机电系统动力学仿真matlab课后答案刘白雁课件.docVIP

机电系统动力学仿真实验报告 机工 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) 1 0.5 z 0 -0.5 10 5 10 5 0 0 -5 -5 -10 -10 y x 1.7 （ 1） （ 2） 1.11 1 0.5 z 0 -0.5 10 5 10 5 0 0 -5 -5 -10 -10 y x 习题二 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) 1.6 =0.2 =0.4 1.4 =0.6 =0.8 1.2 1 0.8 0.6 0.4 0.2 0 0 2 4 6 8 10 12 14 16 18 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) 0.08 0.06 0.04 0.02 0 100 50 50 0 0 -50 -50 0.5 0.4 0.3 0.2 0.1 0 100 50 50 0 0 -50 -50 0.5 0.4 0.3 0.2 0.1 0 150 150 100 100 50 50 0 0 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-) 8 6 4 2 0 1 0.5 1 0.5 0 0 -0.5 -0.5 -1 -1 习题四 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