高维时变非线性微分方程求解
请大家帮忙看一下，这是有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)=.
Error in ==> donglixue at 52
y0=y(length(y),:);%%%length(y),y的最后一行，所有列付给y0
回复 楼主 greatchina 的帖子
楼主的这个问题有点复杂，运动的状态随时间的变化而变化
有人有用过stateflow加simulink解决方程问题的吗
有的话，就是你可以帮楼主啊 {:{10}:} 本帖最后由 oneonly 于 2015-10-14 10:04 编辑
楼主您好，这个问题是怎么解决的，我的是变系数的112个二阶微分方程组，和您这个有几分相似，解的时候遇到好多问题。。。能指点一二吗，谢谢谢谢http://forum.vibunion.com/thread-137719-1-1.html
页:
[1]