北京交通大学MATLAB实验3

北京交通大学MATLAB 实验3

实验内容

一、 在同一坐标系下绘制t 2，- t 2，t 2sin t 在[0,2]t π∈内的曲线图；选择合适的“曲线线型”，“曲线颜色”、“标记符号”选项，并在图形上加注坐标轴名和图名。 仿真程序：

%实验三，题1

close all;clear all;clc;

t=linspace(0,2*pi,1024);

y1=t.^2;

y2=-t.^2;

y3=(t.^2).*sin(t);

figure(2);plot(t,y1,':y',t,y2,'--g',t,y3,'-r');

legend('y=t^2','y=-t^2','y=t^2*sin(t)');

title('三个函数图形');xlabel('t');ylabel('f(t)');

仿真结果：

二、 选择合适的θ的范围，讲同一图形窗口分割成2行2列绘制下列4幅极坐标图 2

3(1) 1.0013(2)cos(3.5)sin (3)(4)1cos (7)

ρθρρρθθθ

θ====-

仿真程序：

%实验三，题2

theta=linspace(0,16*pi,2048);

subplot(2,2,1);polar(theta,1.0013*theta.^2);

subplot(2,2,2);polar(theta,cos(3.5*theta));

subplot(2,2,3);polar(theta,sin(theta)./theta);

subplot(2,2,4);polar(theta,1-(cos(7*theta)).^3);

仿真结果：

三、 用鼠标左键在图形窗口上取5个点，在每个点的位置处写出一个字符串来显式鼠标

点的横坐标值，然后将这些点连成折线。

仿真程序：

%实验三，题3

close all;clc;clear all;

axis([0,10,0,10]);

hold on;

box on;

x=[];y=[];

while(1)

[x1,y1,button]=ginput(1);

if(button~=1)

break;

end

plot(x1,y1,'ro');

x=[x,x1];

y=[y,y1];

text((x1+0.1),y1,mat2str(x1));

end

line(x,y);

hold off;

仿真结果：

四、求积分

4

()

a f x dx

=⎰，其中

2

ln,

()ln16

,2

2sin(1)

2

x if x

f x

if x

xπ

⎧

⎪

=⎨

>

⎪++

⎩

≤

1.9741

仿真程序：

%实验三，题4

close all;clear all;clc;

syms y1 y2 x;

y1=int(log(x^2),x,0,2);

y2=int(log(16)/(2+sin((x+1)*pi)),x,2,4);

y=y1+y2;

disp(double(y));

仿真结果：

五、求解方程x5+6x4-3x2=10的5个根，并将其位置用五角星符号标记在复平面上，要

求横纵坐标轴刻度等长，注明虚轴和实轴，在title位置上写出方程。

仿真程序：

%实验三，题5

clc;close all;clear all;

factor=[1 6 0 -3 0 -10];

ro=roots(factor);

disp(ro);

scatter(real(ro),imag(ro),'p');

axis equal;

title('x^5+6x^4-3x^2=10');

xlabel('real axis');ylabel('imagine axis(j)');

仿真结果：

六、 某班同学成绩已经存放在矩阵A 中，每行为某一位同学的数据，第1列为学号，

第2列至第4列为其三门课程的成绩，试编程按照3门课平均成绩由小到大的顺序重新排列成绩表，并放在矩阵B 中。

99234012957388992340338477809923400966807299234067929359A ⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦

仿真程序:

%实验三，题6

clc;clear all;close all;

A=[99234012,95,73,88;

99234033,84,77,80;

99234009,66,80,72;

99234067,92,93,59];

[y,index]=sort(mean(A(:,2:end),2));

B=A(index,:);

disp(B);

仿真结果：

七、 解如下微分方程，并绘制出图像y (x ) 2002()0,|0[0,1,|]1,x x y y y y x ==-===∈- 仿真程序：

function xdot=f_03(t,x)

xdot=[x(2);2*(x(2))^2];

t0=0;tf=1;

x0=[0,-1];

[t,x]=ode45('f_03',[t0,tf],x0);

plot(t,x(:,1));

title('y=f(x)');

xlabel('x');ylabel('y');

仿真结果：