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

机电系统动力学仿真实验报告

机工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.7
（1）

（2）

1.11

习题二
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')

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)

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-')

习题四
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