新人求教多自由度振动问题（Runge-Kutta法）

Jonahting

我最近建了一个三个自由的弹簧阻尼模型，列出了相关的三个二阶微分方程，尝试用Runge-Kutta法求解，但结果却不对：
这是该模型所受到的外部激励，，是个分段函数。
0.*(t>=0&t<10)+192000*((2*0.06*(40/3.6).*(t-10))/0.4).*(t>=10&t<10.018)+192000*((2*0.06*(0.4-(40/3.6).*(t-10)))/0.4).*(t>=10.018&t<10.036)+0.*(t>=10.036&t<=30)




这是主程序，激励，填在下面代码的引号内。
t=0:0.01:30;
[url=mailto:odefun=@(t,z)[z(2);(‘                 ’214000*z(1)+22000*z(3)-1600*z(2)+1600*z(4))/45.4;z(4);(22000*z(1)-29500*z(3)+7500*z(5)+1600*z(2)-1900*z(4)+300*z(6))/318;z(6);(7500*z(3)-7500*z(5)+300*z(4)-300*z(6))/70]odefun=@(t,z)[z(2);(‘                 ’214000*z(1)+22000*z(3)-1600*z(2)+1600*z(4))/45.4;z(4);(22000*z(1)-29500*z(3)+7500*z(5)+1600*z(2)-1900*z(4)+300*z(6))/318;z(6);(7500*z(3)-7500*z(5)+300*z(4)-300*z(6))/70[/url]];
[t,z]=ode45(odefun,t,[0;0;0;0;0;0;]);
plot(t,z(:,1),'r-')
出图结果是一条直线，不合理。。。谢谢大家。烦了好多天了