%**************************************************************************
%（第四題）求矩陣的LDLT分解及cholesky分解--3111312053龔小林?碩1080
%**************************************************************************
clear
clc
format?short
%待分解計算矩陣A
for?i=1:10
????for?j=1:10
????????if?j==i
????????????A（ij）=i;
????????else
????????????A（ij）=min（ij）-2;
????????end
????end
end
disp（‘待分解矩陣A為‘）
A
AG=A;
%矩陣LDLT分解
for?k=1:10
????for?p=1:k-1
????????r（p）=A（pp）*A（kp）;
????end
????sum=0;
????for?p=1:k-1
????????sum=sum+A（kp）*r（p）;??
????end
????????A（kk）=A（kk）-sum;
????if?A（kk）==0
????????disp（‘錯誤‘）
????????break;
????end
????D（k