%矩陣分解
%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