几种矩阵分解算法的matlab实现

%矩阵分解
%Cholesky分解
%产生n阶pascal矩阵，pascal矩阵由pascal三角阵发展而来，为对称正定阵
disp（‘正定对称阵的分解‘）
n=6;
X=pascal（n）
R=chol（X）
disp（‘非正定对称阵的分解‘）
X（nn）=X（nn）-1;
X
[Rp]=chol（X）
Xn_1=R‘*R
disp（‘利用pascal分解结果计算线性方程组‘）
A=pascal（6）;
b=[1 3 5 7 9 4]‘;
x=A\b
R=chol（A）;
xc=R\（R‘\b）



%不完全cholesky分解
S=sparse（[1  0   3  0;
          0  25  0  30;
          3  0   9  0 ;
          0  30  0  661]）;
      disp（‘完全Cholesky分解‘）
      [Rp]=cholinc（S‘0‘）
      Rfull=full（R）
      Srfull=Rfull‘*Rfull
      disp（‘不完全Cholesky分解‘）
      Rinf=cholinc（S‘inf‘）
      Rinf_full=full（Rinf）
      Srinf_full=Rinf_full*Rinf_full
      
      
      
      
      
      %Lu矩阵分解
      A=[1 2 3
          4 5 6
          7 8 0];
      disp（‘Lu分解法1‘）;
      [L1U]=lu（A）
      X=L1*U;
      disp