高维时变非线性微分方程求解

greatchina

请大家帮忙看一下，这是有94个二阶非线性时变微分方程组成的方程组，运行结果总是返回NAN，并提示下面的错误，不知道是什么原因，而且各种ODE方法都尝试过，均没有成功。如果将tspan=linspace(t0,tf);改为tspan=linspace(t0,tf,188);（这样改是不符合原理的）就不会出现该错误，但返回值仍然是NAN。
请高手指教，共同学习！

主程序及调用程序见附件

运行错误：

??? Attempted to access y(188,:); index out of bounds because size(y)=[100,188].
Error in ==> donglixue at 52
    y0=y(length(y),:);%%%length(y),y的最后一行，所有列付给y0

xukai871105

楼主的这个问题
有点复杂，运动的状态随时间的变化而变化
有人有用过stateflow加simulink解决方程问题的吗

有的话，就是你可以帮楼主啊

haidi114

oneonly

楼主您好，这个问题是怎么解决的，我的是变系数的112个二阶微分方程组，和您这个有几分相似，解的时候遇到好多问题。。。能指点一二吗，谢谢谢谢http://forum.vibunion.com/thread-137719-1-1.html