声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2359|回复: 14

[综合讨论] 请教各位求解12变量非线性方程组!

[复制链接]
发表于 2007-12-5 10:56 | 显示全部楼层 |阅读模式

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

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

x
我正在做一个课题,列出非线性方程组如下:MX''+CX'+KX=F,其中M、C、K分别为12*12的矩阵,F为12*1的列向量。我用的是MATLAB调用ode45,可是结果只有前一个周期的图像,后面就发散了,不知为什么呢?我看别人都用五阶变步长自适应龙格库塔进行求解,结果都很理想,我一直很纳闷,问题一直得不到解决,不知是参数问题还是算法有问题,请教各位大虾帮忙!!非常感谢!!!

[ 本帖最后由 eight 于 2007-12-5 11:16 编辑 ]
回复
分享到:

使用道具 举报

发表于 2007-12-5 13:59 | 显示全部楼层
那你也用五阶变步长自适应龙格库塔解啊
发表于 2007-12-5 14:26 | 显示全部楼层
有可能是参数问题,用ode45最容易出现的问题就是发散,刚性方程
 楼主| 发表于 2007-12-6 16:54 | 显示全部楼层

是参数问题还是算法问题?

是我参数有问题导致发散,还是ode45函数导致了发散??参数都是参考别人数据得来的啊!是不是换种工具或算法?能给点儿建议吗?
谢谢!!!
发表于 2007-12-6 17:24 | 显示全部楼层
请将你的代码及问题公式贴一下,否则只能是空谈.
 楼主| 发表于 2007-12-7 18:17 | 显示全部楼层

这是具体参数和程序!请各位帮忙看看问题出在哪儿???

CarryOut.m

1.75 KB, 下载次数: 17

运行结果程序

DYD.m

3.47 KB, 下载次数: 19

时变参数矩阵

 楼主| 发表于 2007-12-7 18:28 | 显示全部楼层

改变积分区间时各图形

积分区间t=0~20

积分区间t=0~20

积分区间t=0~30

积分区间t=0~30

积分区间t=0~50

积分区间t=0~50
 楼主| 发表于 2007-12-7 18:30 | 显示全部楼层

请各位分析下原因!

非常感谢,这个难题我已经很长时间了,都得不到解决!

[ 本帖最后由 eight 于 2007-12-7 19:44 编辑 ]
发表于 2007-12-7 19:13 | 显示全部楼层
可以用ode15s试试看,这个对刚性方程比较有效。
发表于 2007-12-7 22:05 | 显示全部楼层

回复 #4 fineshang 的帖子

可能是算法的原因,但是之前你最好看一下参数
发表于 2007-12-7 22:23 | 显示全部楼层
我曾经也遇到类似的问题,仔细看看吧,可能程序错了
发表于 2007-12-7 22:37 | 显示全部楼层

回复 #7 fineshang 的帖子

如果你的程序一直没变,而改变了积分区间。

那么你画得图都是正确的,随着积分区间的增加,结果也在递增。

特别是40以后,仔细看一下几幅图y轴的范围你就可以明白了:-3--3;-5--20;-1--7*10^5

因此单位不同时,图形是没有可比性的
发表于 2007-12-8 15:09 | 显示全部楼层
请将公式贴一下。
------------------------------------------------------------------
看了一下你的代码,目前只能作如下判断:
%%%
变化tspan, tfinal=35可明显看到长时间趋于发散。
变化初值, 若Y0=rand(24,1), tfinal=10就可明显看到趋于发散,而且行为也与Y0=0很不一样!
所以如果没有输入错误的话,很可能这个系统不稳定。自设定时间步长也无改善.
------------------------
另外,ode45已经采用了变步长求解,所以原则上ode45不能求解的问题,你自己编程用5阶RK求解的话,情况应该不会有什么改善,不过你还是可以试试。

评分

1

查看全部评分

 楼主| 发表于 2007-12-9 10:14 | 显示全部楼层

谢谢各位的解答与讨论!

请遇到过这类问题的朋友给出解决办法?谢谢~
其实我已经仔细检查过系统参数了
好像没有问题
会继续关注这个问题
我做的是行星齿轮的动力学分析
请做过这方面的朋友多多讨论交流!!!
发表于 2012-11-29 11:36 | 显示全部楼层
同样遇到这样的问题
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-29 02:37 , Processed in 0.090301 second(s), 23 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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