[运算] ode解震动方程式 时间拉长有误差

楼主: dooksp (颓废的人)   2017-07-05 11:01:06
各位好
目前在利用ode解简单的力学问题
方程式为 F=-kx
当我的dt设置的范围再0至10(时间很短的情况下)曲线在+-1震荡
但当我dt设置范围改为0至100时(时间拉长了)曲线后来就不在+-1震荡了
想询问是否有方法解决这个问题
感谢
以下是我写
function dydt=odefcn(t,y,m,b,k,w)
dydt=zeros(2,1);
dydt(1)=y(2);
dydt(2)=-(k/m)*y(1);
作者: sin55688 (单手挑藏獒)   2017-07-05 20:31:00
依旧在+-1之间阿
楼主: dooksp (颓废的人)   2017-07-07 10:47:00
又改成 dt=0:0.01:100吗

Links booklink

Contact Us: admin [ a t ] ucptt.com