机电系统动力学仿真实验报告机工 1004 班 201003130150宋康习题一1.3 Abc wu_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.2 2.4 2.5 习题三3.1 clf,clc,clear t=(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.2 clf,clear x=-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 on 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.505010015005010015000.10.20.30.40.53.3 clc; clear all; close all; 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.4 5.5 5.6 5.7 习题六6.2 6.3 6.4 习题七7.1 （1）临界稳定（2）不稳定（3）不稳定（4）不稳定7.2 K=36 7.3 7.4