第八天 矩阵分界 稀疏矩阵 奥力给！！！

 1.cholesky分解又称平方根法

用chol（）函数实现

R=chol（X）：X为正定矩阵，R是上三角矩阵 X=RR‘。X不是正定的，则返回错误信息

[R,p]=chol(X) 返回两个参数，不返回错误信息。当X为正定矩阵时，返回上三角矩阵R满足X=RR’，且p=0;当是非正定矩阵时，返回值p是正整数，R是上三角矩阵，阶数为p-1，X（1：p-1，1：p-1）=RR‘

2.一般方阵高斯消去法分解又称LU分解

将任意一个方阵A分解为一个下三角矩阵L和一个上三角矩阵U的乘积，A=LU。

用lu（）函数来实现

[L,U]=lu(x)   x为方阵满足上述关系

[L,U,P]=lu(x),P为置换矩阵PX=LU

Y=lu（x）上三角和下三角矩阵结合，对角元素为上三角矩阵Y=L+U-I。I为单位矩阵。置换矩阵p丢失

上面两种方法都大大降低了计算Ax=b的量

LUx=b   RR’x=b     x=.....除一下就行了

luinc（）做不完全LU分解