fineshang 发表于 2007-12-5 10:55

请教各位求解12变量非线性方程组！

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

messenger 发表于 2007-12-5 13:18

这东西发散很正常吧

你看到的都是人家做的结果比较理想的，不理想的也不会给别人看

建议你换别的方法试试

风花雪月 发表于 2007-12-5 16:37

收敛性和你的模型本身有很大的关系

fineshang 发表于 2007-12-7 13:15

一个方程组不收敛的原因是什么呢？

是刚性导致的么？怎么改善啊？请指教！谢谢！！！

风花雪月 发表于 2007-12-11 15:18

原帖由 fineshang 于 2007-12-7 13:15 发表 http://www.chinavib.com/forum/images/common/back.gif
是刚性导致的么？怎么改善啊？请指教！谢谢！！！

不一定，刚性仅仅是总原因
如果是刚性造成的可以用带算法阻尼的算法进行计算
页: [1]
查看完整版本: 请教各位求解12变量非线性方程组！