北京交通大学matlab上机实验3.doc

第三次上机试验报告 第一题：在同一坐标系下绘制t2，- t2，t2sint在内的曲线图；选择合适的“曲线线型”，“曲线颜色”、“标记符号”选项，并在图形上加注坐标轴名和图名。 t=0:0.1:2*pi;y1=t.^2;y2=-y1;y3=y1.*sin(t); plot(t,y1,--r,t,y2,:c,t,y3,gd) xlabel(t),ylabel(y),title(曲线图) 第二题：选择合适的的范围，将同一图形窗口分割成2行2列绘制下列4幅极坐标图 解： t=0:0.1:8*pi; y1=1.0013.*t.^2; y2=cos(3.5*t); y3=sin(t)./t; y4=1-(cos(7.*t)).^3; subplot(2,2,1),polar(t,y1); subplot(2,2,2),polar(t,y2); subplot(2,2,3),polar(t,y3); subplot(2,2,4),polar(t,y4); 第三题：用鼠标左键在图形窗口上取5个点，在每个点的位置处写出一个字符串来显式鼠标点的横坐标值，然后将这些点连成折线。 axis([0,5,0,5]); hold on; box on; x=[];y=[]; while(1) [x1,y1,button]=ginput(1); if(button~=1)break; end plot(x1,y1,+); x=[x,x1];y=[y,y1]; s=num2str(x1); text(x1,y1,s); end line(x,y);hold off： 第四题：求积分，其中 f1=inline(log(x.^2),x); f2=inline(log(16)./(2+sin((x+1).*pi)),x); [s1,kk]=quad(f1,0,2); [s2,kk]=quad(f2,2,4); s=s1+s2 s =1.9741 第五题：五、 求解方程x5+6x4-3x2=10的5个根，并将其位置用五角星符号标记在复平面上，要求横纵坐标轴刻度等长，注明虚轴和实轴，在title位置上写出方程。 解： a=[1,6,0,-3,0,-10]; b=roots(a); b1=real(b); bi=imag(b); axis([-4,4,-4,4]) hold on plot(b1,bi,p) 第六题：某班同学成绩已经存放在矩阵A中，每行为某一位同学的数据，第1列为学号，第2列至第4列为其三门课程的成绩，试编程按照3门课平均成绩由小到大的顺序重新排列成绩表，并放在矩阵B中。 A=95 73 8884 77 8066 80 7292 93 59] for i=1:4 c(i,2)=i; c(i,1)=(A(i,2)+A(i,3)+A(i,4))/3; end for i=1:4 for j=1:(4-i) if c(j+1,1)c(j,1) temp=c(j,:); c(j,:)=c(j+1,:); c(j+1,:)=temp; end end end for i=1:4 B(i,:)=A(c(i,2),:); end B = 66 80 72 84 77 80 92 93 59 95 73 88 第七题：解如下微分方程，并绘制出图像y(x) 解： function xdot= f1(t,x) xdot=zeros(2,1); xdot(1)=x(2); xdot(2)=2*(x(2).^2); end t0=0 t0 = 0 tf=10 tf = 10 x0=[0 -1] x0 = 0 -1 [T,X]=ode45(f1,[t0 tf],x0); plot(T,X)