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

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

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

-5

5

10

-10

-5

5

10

-0.5

0.5

1

x

y

z

1.7

（1）

（2）

1.11

-5

5

-50

5

x

y

z

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

-0.5

0.5

1

-0.5

0.5

1

02

4

6

8

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