Jonahting 发表于 2015-4-27 21:37

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

我最近建了一个三个自由的弹簧阻尼模型,列出了相关的三个二阶微分方程,尝试用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;
odefun=@(t,z)];
=ode45(odefun,t,);
plot(t,z(:,1),'r-')
出图结果是一条直线,不合理。。。谢谢大家。烦了好多天了
页: [1]
查看完整版本: 新人求教多自由度振动问题(Runge-Kutta法)