高分求一道MATLAB 题 答案，题不难，MATLAB高手入内，在线等！！
有一道MATLAB题，没标准答案，求掌握MATLAB的童鞋帮忙给个答案，本人在线等待。题内容在附图里面：
某弹簧-质量-阻尼器系统如图所示。
振体质量为m，弹簧的线性系数为k，阻尼器的阻尼系数为c，用x(t)表示系统的位移，系统受一外力F(t)=X0kf(t)作用， 为固有频率，τ = ωnt 。
试求：
（1）写出系统的运动方程；
（2）将系统的运动方程转化为状态空间形式；
（3） 考虑当系统为线性系统时，在同一座标下，求当0≤τ≤40范围内，阻尼比ζ分别为0.1、0.5、1.0、5.0时的自由振荡响应（即x(t)与τ的关系）。
要求：写出解题过程、相应的MatLab程序脚本，并附相应的仿真曲线图。
谢谢了，本人积分4000+，只要帮我解决了这道题，我积分你要多少全给你。
如果有做好了的，发到我的QQ邮箱吧，[email protected]
我一直在线的，Hi我，或者加我QQ 584195364都可以的，谢谢了，悬赏只能100分，收到答案我追加分。
希望能把解题思路，MATLAB程序、相应的图线（截图）一起发给我。不甚感激！
figure('pos',[120,130,300,250],'color','w')%设置图形句柄
axes('pos',[0.1 0.15 0.8 0.72]);%设置坐标轴句柄
syms s%设置符号s
for zeta=[0.1,0.5,1.0,5.0];%zeta为阻尼比，设置循环计算，并绘制图形
wn=0.15;%设置固有频率为0.15
ezplot(ilaplace(wn^2/s/(s^2+2*zeta*wn*s+wn^2)),[0 100]);hold on;%wn^2/(s^2+2*zeta*wn*s+wn^2为系统传递函数，wn^2/s/(s^2+2*zeta*wn*s+wn^2为阶跃响应下的传递函数，仿真时间为100s
if zeta==.1,text(18,1.85,'0.1');end;%标出0.1的曲线
if zeta==.5,text(20,1.25,'0.5');end;%标出0.5的曲线
if zeta==1,text(25.5,0.8,'1.0');end;%标出1的曲线
if zeta==5,text(30,0.2,'5.0');end;%标出5的曲线
end;
grid;%打开grid
axis([0 100 0 2]);title('阻尼比:0.1,0.5,1.0,5.0')%设置坐标轴长度及标题
微分方程及状态空间用见所插入的图片。因为只能插入一张图片，所以结果曲线无法再次插入了，但是只要运行以上程序，就会得到结果曲线。
