机电系统动力学仿真实验报告 机工1004 班 201003130150宋康 习题一 1.3 Abc w u_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.51xyz 1.7 （1） （2） 1 .1 1 -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.8 3.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) -5 005 0-5 005 01 0 000 .0 20 .0 40 .0 60 .0 8 -5 005 0-5 005 01 0 000 .10 .20 .30 .40 .5 05 01 0 01 5 005 01 0 01 5 000 .10 .20 .30 .40 .5 3.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
