原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
. .
机电系统动力学仿真实验报告
机工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
文档评论（0）