机电系统动力学仿真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')

z

y

-10

x

10

1.7

（1）

（2）

1.11

z

y

-10

x

10

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

-50

50

-50

50

150

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

-1

1

习题四 4.1

习题五

5.1

5.4

5.5

5.6

5.7

习题六 6.2

6.3

6.4