单自由度系统自由振动matlab程序

单自由度系统自由振动根据物理模型，列得系统的运动方程：0mxcxkx或220nnxxx2/ncm2/nkm2222dnn对系统参数取：m=1,c=7，k=9假定有初始条件：4x，0x7/6，3n列的系统MATLAB程序：SymskmQWnWdyXEtx0v0;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(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程序框图：运行后分别得到仿真图: