北京交通大学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幅极坐标图

23(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 中。

99234012

957388992340338477809923400966807299234067

9293

59A ?????

?=???

???

仿真程序:

%实验三，题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');

仿真结果：