高斯消元法线性方程组python_线性方程组的直接解法

我们学过高等代数(或者线性代数）知道，求解线性方程组可以用克莱姆法则，但是在实际交给机器计算时，因为克莱姆法则需要计算乘除(n-1)（n+1）!+n次，计算量太复杂。我们需要寻找更好的数值解法。

高斯消元法

它包括两部分：消元和回代。

它的乘除计算量约为

。但是这个方法有一个缺陷就是可能分母会是接近于0的数，导致除出来的数误差大。所以我们需要对它进行改进。出现了列主元消去法。

列主元消去法

它只是在高斯消元法的基础上，进行了找每列绝对值最大的数作为除数的步骤。

因此运算量与高斯消去法相当。

直接三角分解法（LU法）

转化为求解两个简单三角形方程问题。

则求解下面两个方程即可

上述A的分解也叫Doolittle分解。它的乘除计算量大约

,与高斯消元法相当。

好处：若实现了A的D分解，且L、U保存在A的相应位置。则对于系数相同，只是b不同的方程组而言，只需计算

次了。

平方根法

针对的是正定对称矩阵。

对于对称矩阵，若它的所有顺序主子式均不为0，则A可唯一分解为

。

对于正定对称矩阵，说明上述D中的对角线上的元素都是正的。

若A是正定对称矩阵，则存在一个实的非奇异下三角矩阵L,使得

，当限定L的对角元素为正时，这种分解是唯一的。

平方根法约需

乘除次，大约为LU法的一半。

追赶法

追赶法解决的是系数矩阵满足上述条件的方程组

它将A分解为下面的形式

计算量为5n-4次乘除。后面均为高斯消元法的变形。