单自由度系统自由振动根据物理模型，列得系统的运动方程：或 对系统参数取： m=1 , c=7 ， k=9假定有初始条件： ， ， 列的系统 MATLAB 程序：Syms k m Q Wn Wd y X E t x0 v0 ;t=0::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(x0^2+(v0+E*Wn*x0)^2/Wd^2);%初始振幅Q=atan((v0+E*Wn*x0)./(x0*Wd));%初始相角y=X*exp(-E*Wn*t).*cos(Wd*t-Q);plot(t,y)画得 Simulink 程序框图：运行后分别得到仿真图:
