实用标准文案精彩文档机电系统动力学仿真实验报告机工1004班201003130150宋康习题一1.3Abcwu_2004两个是合法的。1.4（1）（2）（3）clear;x=-8:0.5:8;y=x';X=ones(size(y))*x;Y=y*ones(size(x));实用标准文案精彩文档R=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R;mesh(X,Y,Z);colormap(hot)xlabel('x'),ylabel('y'),zlabel('z')-10-50510-10-50510-0.500.51xyz1.7（1）（2）实用标准文案精彩文档1.11-10-50510-10-50510-0.500.51xyz习题二2.1实用标准文案精彩文档2.22.4实用标准文案精彩文档2.5实用标准文案精彩文档习题三3.1clf,clc,cleart=(0:18/100:18)';xi=0.2:0.2:0.8;beita=sqrt(1-xi.^2);sita=atan(beita./xi);y=1-exp(-t*xi).*sin(t*beita+ones(101,1)*sita)./(ones(101,1)*beita);figure(1)plot(t,y(:,1),'b:',t,y(:,2),'r-.',t,y(:,3),'g-',t,y(:,4),'ko')实用标准文案精彩文档legend('\xi=0.2','\xi=0.4','\xi=0.6','\xi=0.8')02468101214161800.20.40.60.811.21.41.6=0.2=0.4=0.6=0.83.2clf,clearx=-50:1:50;y=-40:1:60;z=1./(sqrt((1-x).^2+y.^2)+sqrt((1+x).^2+y.^2));figure(1)plot3(x,y,z),boxon[X,Y]=meshgrid(x,y);Z=1./(sqrt((1-X).^2+Y.^2)+sqrt((1+X).^2+Y.^2));figure(2)mesh(X,Y,Z)figure(3)surf(Z)实用标准文案精彩文档-50050-5005010000.020.040.060.08-50050-5005010000.10.20.30.40.5实用标准文案精彩文档05010015005010015000.10.20.30.40.53.3clc;clearall;closeall;t=0:pi/50:2*pi;x=sin(t);y=cos(t);z=t;plot3(x,y,z,'g-')实用标准文案精彩文档-1-0.500.51-1-0.500.5102468习题四4.1实用标准文案精彩文档习题五5.1实用标准文案精彩文档5.45.5实用标准文案精彩文档5.6实用标准文案精彩文档5.7习题六6.2实用标准文案精彩文档6.3实用标准文案精彩文档6.4实用标准文案精彩文档实用标准文案精彩文档习题七7.1（1）临界稳定（2）实用标准文案精彩文档不稳定（3）实用标准文案精彩文档不稳定（4）实用标准文案精彩文档不稳定7.2K=36实用标准文案精彩文档实用标准文案精彩文档实用标准文案精彩文档7.3实用标准文案精彩文档实用标准文案精彩文档7.4实用标准文案精彩文档实用标准文案精彩文档
