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

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

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

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

x

我正在做一个课题，列出非线性方程组如下：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 发表

 登录/注册后可看大图

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

不一定，刚性仅仅是总原因
如果是刚性造成的可以用带算法阻尼的算法进行计算

账号		自动登录	找回密码
密码			我要加入

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

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

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