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

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

《机电系统动力学仿真matlab课后答案刘白雁.docx》由会员分享，可在线阅读，更多相关《机电系统动力学仿真matlab课后答案刘白雁.docx（29页珍藏版）》请在冰点文库上搜索。

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

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

机工1004班201003130150宋康

习题一

1.3

Abcwu_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）,boxon

[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;

clearall;

closeall;

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