声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2041|回复: 3

[编程技巧] 高维时变非线性微分方程求解

[复制链接]
发表于 2009-2-23 19:47 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
请大家帮忙看一下,这是有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

donglixue.m

2.66 KB, 下载次数: 24

主程序

fangcheng.m

75.11 KB, 下载次数: 20

调用程序

回复
分享到:

使用道具 举报

发表于 2009-2-23 21:16 | 显示全部楼层

回复 楼主 greatchina 的帖子

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

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

评分

1

查看全部评分

发表于 2014-12-18 16:53 | 显示全部楼层
发表于 2015-10-14 10:02 | 显示全部楼层
本帖最后由 oneonly 于 2015-10-14 10:04 编辑

楼主您好,这个问题是怎么解决的,我的是变系数的112个二阶微分方程组,和您这个有几分相似,解的时候遇到好多问题。。。能指点一二吗,谢谢谢谢http://forum.vibunion.com/thread-137719-1-1.html
%B$T}XGJVFY92(_0C7VX7`B.png
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-3 10:07 , Processed in 0.076203 second(s), 26 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表