高斯列主元素消去法求解矩阵方程

zxl1334 · 发表于 2009-3-3 12:11

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

x

高斯列主元素消去法求解矩阵方程，供参考

logxing · 发表于 2009-3-3 13:21

附件传错了吧

logxing · 发表于 2009-3-3 14:26

我也做过高斯列主元素消去法的程序，可是当行列式接近于零的时候，即使阶数不是很大（仅12阶），得到的结果误差也很大（和matlab）相比。行列式如果是10^-200量级的话可能差个10^15倍左右。看来对病态矩阵不太适用。

lq12131010 · 发表于 2009-3-3 17:25

全部选主元！

logxing · 发表于 2009-3-4 08:17

全选主元得记录列交换的顺序，太麻烦我就没做。非病态矩阵列主元算的还算准。病态矩阵我估计全选主元也不一定算得准吧。不知道matlab的内部算法是怎样的原理。

lq12131010 · 发表于 2009-3-4 10:28

搞研究应该不怕麻烦！
而且不是很麻烦做个解顺序调换就可以。
全选肯定是准的除非是奇异阵。

logxing · 发表于 2009-3-4 12:11

嗯，回来做一个和matlab比较一下。

lilili606 · 发表于 2009-3-5 22:20

选主元其实也不算麻烦，只是在循环里交换一下顺序就好了，matlab的精度也没那么好，病态方程一般用奇异值分解做吧。。。很麻烦的说。。。

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

[Fortran] 高斯列主元素消去法求解矩阵方程