lixixi0631032的专栏

原创 方形矩阵的LU分解法

将系数矩阵A转变成等价两个矩阵L和U的乘积 ，其中L和U分别是单位下三角矩阵和上三角矩阵。当A的所有顺序主子式都不为0时，矩阵A可以分解为A=LU。其中L是下三角矩阵，U是上三角矩阵。 根据上图可以得到一般的思路，最下面这个方阵中的元素可表示为A[i,j]=ΣL[i,k]U[k,j]，其中0=<k<=min(i,j).我们应该从最左和最上开始计算，依次向内层推进。 """Lower-Upper (LU) Decomposition.""" # lower–upper (LU) decompo