请教各位求解12变量非线性方程组!
我正在做一个课题,列出非线性方程组如下:MX''+CX'+KX=F,其中M、C、K分别为12*12的矩阵,F为12*1的列向量。我用的是MATLAB调用ode45,可是结果只有前一个周期的图像,后面就发散了,不知为什么呢?我看别人都用五阶变步长自适应龙格库塔进行求解,结果都很理想,我一直很纳闷,问题一直得不到解决,不知是参数问题还是算法有问题,请教各位大虾帮忙!!非常感谢!!![ 本帖最后由 eight 于 2007-12-5 11:16 编辑 ] 那你也用五阶变步长自适应龙格库塔解啊 有可能是参数问题,用ode45最容易出现的问题就是发散,刚性方程
是参数问题还是算法问题?
是我参数有问题导致发散,还是ode45函数导致了发散??参数都是参考别人数据得来的啊!是不是换种工具或算法?能给点儿建议吗?谢谢!!! 请将你的代码及问题公式贴一下,否则只能是空谈.
这是具体参数和程序!请各位帮忙看看问题出在哪儿???
改变积分区间时各图形
请各位分析下原因!
非常感谢,这个难题我已经很长时间了,都得不到解决![ 本帖最后由 eight 于 2007-12-7 19:44 编辑 ] 可以用ode15s试试看,这个对刚性方程比较有效。
回复 #4 fineshang 的帖子
可能是算法的原因,但是之前你最好看一下参数 我曾经也遇到类似的问题,仔细看看吧,可能程序错了回复 #7 fineshang 的帖子
如果你的程序一直没变,而改变了积分区间。那么你画得图都是正确的,随着积分区间的增加,结果也在递增。
特别是40以后,仔细看一下几幅图y轴的范围你就可以明白了:-3--3;-5--20;-1--7*10^5
因此单位不同时,图形是没有可比性的 请将公式贴一下。
------------------------------------------------------------------
看了一下你的代码,目前只能作如下判断:
%%%
变化tspan, tfinal=35可明显看到长时间趋于发散。
变化初值, 若Y0=rand(24,1), tfinal=10就可明显看到趋于发散,而且行为也与Y0=0很不一样!
所以如果没有输入错误的话,很可能这个系统不稳定。自设定时间步长也无改善.
------------------------
另外,ode45已经采用了变步长求解,所以原则上ode45不能求解的问题,你自己编程用5阶RK求解的话,情况应该不会有什么改善,不过你还是可以试试。
谢谢各位的解答与讨论!
请遇到过这类问题的朋友给出解决办法?谢谢~其实我已经仔细检查过系统参数了
好像没有问题
会继续关注这个问题
我做的是行星齿轮的动力学分析
请做过这方面的朋友多多讨论交流!!! 同样遇到这样的问题
页:
[1]