w=2*pi; x0=0;v0=0; h=H/m;
b=h/(wn^2-w^2);
A=sqrt(x0^2+(v0-b*w)^2/wn^2); a=atan(wn*x0/(v0-b*w)); t=0:0.005:5;
x=A*sin(wn*t+a)+b*sin(2*pi*t); Plot(t,x)；
所得曲线如下：
第 20 页 共 29 页
3.1.3 simulink仿真：
3.2有阻尼
3.2.1Soloidworks运动仿真
第 21 页 共 29 页
3.2.2 simulink仿真：
3.3.3 MATLAB编程
运动方程如下：
由初始条件可求得A和θ。 程序如下：
m=1.6632; k=15; H=1;
wn=sqrt(k/m); w=2*pi; h=H/m; c=7;
n=c/(2*m);
wd=sqrt(wn^2-n^2);
b=h/(sqrt((wn^2-w^2)^2+4*n^2*w^2)); o=atan(2*n*w/(wn^2-w^2)); y=sin(o); y1=cos(o);
A=-sqrt((n*b*y-b*w*y1)^2/wd^2+b^2*y^2); a=atan(wd*y/(n*y-w*y1));
第 22 页 共 29 页
t=0:0.05:5;
x=-(A*(exp(-n*t)).*sin(wd*t+a)+b*sin(w*t-o)); Plot(t,x);
所得曲线如下：
三、两自由度振动
给定两自由度系统参数，对两自由度自由振动系统进行计算，绘
制响应曲线，并进行仿真；对其中末端物体上施加一简谐力，绘制受迫振动曲线，并进行仿真。
要求：编程计算软件不限，只要绘出曲线即可。仿真使用solidwoks
第 23 页 共 29 页
和matlab/simulink两种软件同时进行,计算结果和仿真结果一致；
1、自由振动
1.1 Soloidworks运动仿真
1.2 simulink仿真
出图如下：
第 24 页 共 29 页