factorization_cholesky.m

matlab基础程序 便于初学者学习

% factorization_cholesky.m
% Cholesky分解示例
% 产生n阶pascal矩阵,pascal矩阵为pascal三角阵发展而来，为对称正定阵
disp('正定对称阵的分解：')
n=6;
X=pascal(n)
R=chol(X)
disp('非正定对称阵的分解：')
X(n,n)=X(n,n)-1;
X
[R, p]=chol(X)
Xn_1=R'*R
disp('利用cholesky分解结果计算线性方程组')
A=pascal(6);
b=[1 3 5 7 9 4]';
x=A\b
R=chol(A);
xc=R\(R'\b)