北京交通大学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');
仿真结果：。