单自由度系统自由振动根据物理模型,列得系统得运动方程:或 对系统参数取: m=1 , c=7 , k=9假定有初始条件: , , 列得系统MATLAB程序:Syms k m Q Wn Wd y X E t x0 v0 ;t=0:0、01:10;x0=0;%初始位移v0=4;%初始速度E=7/6;%阻尼k=9;%弹性系数m=1;%质量Wn=sqrt(k/m);%固有频率Wd=sqrt(E*Wn)2-(Wn)2);%阻尼自然频率X=sqrt(x02+(v0+E*Wn*x0)2/Wd2);%初始振幅Q=atan(v0+E*Wn*x0)、/(x0*Wd);%初始相角y=X*exp(-E*Wn*t)、*cos(Wd*t-Q);plot(t,y)画得Simulink程序框图:运行后分别得到仿真图: