《机电系统动力学仿真matlab课后问题详解 刘白雁》由会员分享,可在线阅读,更多相关《机电系统动力学仿真matlab课后问题详解 刘白雁(28页珍藏版)》请在金锄头文库上搜索。
1、机电系统动力学仿验报告机工1004班 201003130150宋康习题一1.3Abc wu_2004 两个是合法的。1.4123clear;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)1.7121.11习题二2.12.22.42.5习题三3.1clf,clc,cleart=(0:18/100:18);xi=0.2:0.2:0.8;beita=sqrt(1-xi.2);si
2、ta=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)3.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),box onX,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)3.3clc;clear all;close all;t=0:pi/50:2*pi;x=sin(t);y=cos(t);z=t;plot3(x,y,z,g-)习题四4.1习题五5.15.45.55.65.7习题六6.26.36.4习题七7.11临界稳定2不稳定3不稳定4不稳定7.2K=367.37.4