高斯列主元素消去法求解矩阵方程
高斯列主元素消去法求解矩阵方程,供参考 附件传错了吧 我也做过高斯列主元素消去法的程序,可是当行列式接近于零的时候,即使阶数不是很大(仅12阶),得到的结果误差也很大(和matlab)相比。行列式如果是10^-200量级的话可能差个10^15倍左右。看来对病态矩阵不太适用。 全部选主元! 全选主元得记录列交换的顺序,太麻烦我就没做。非病态矩阵列主元算的还算准。病态矩阵我估计全选主元也不一定算得准吧。不知道matlab的内部算法是怎样的原理。 搞研究应该不怕麻烦!而且不是很麻烦做个解顺序调换就可以。
全选肯定是准的 除非是奇异阵。 嗯,回来做一个和matlab比较一下。 选主元其实也不算麻烦,只是在循环里交换一下顺序就好了,matlab的精度也没那么好,病态方程一般用奇异值分解做吧。。。很麻烦的说。。。
页:
[1]