完整版机电系统动力学仿真matlab课后答案刘白雁Word格式文档下载.docx

完整版机电系统动力学仿真matlab课后答案刘白雁Word格式文档下载.docx

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

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

3）,'

g-'

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:

60;

z=1./（sqrt（（1-x）.^2+y.^2）+sqrt（（1+x）.^2+y.^2））;

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

习题四

4.1

习题五

5.1

5.4

5.5

5.6

5.7

习题六

6.2

6.3

6.4

习题七

7.1

临界稳定

不稳定

（4）

7.2

K=36

7.3

7.4