改进的cholesky matlab,矩阵的Cholesky分解的Matlab实现

矩阵的Cholesky分解的Matlab实现

下次再补上改进的%Cholesky分解法,

%Cholesky分解法

function [X]=m_chol(A,b)

[N,N]=size(A);

X=zeros(N,1);

Y=zeros(N,1);

for i=1:N

A(i,i)=sqrt(A(i,i)-A(i,1:i-1)*A(i,1:i-1)');

if A(i,i)==0

'A is singular. no unique solution'

break

end

for j=i+1:N

A(j,i)=(A(j,i)-A(j,1:i-1)*A(i,1:i-1)')/A(i,i);

end

end

A

b

%前代法

for j=1:N

Y(j)=(b(j)-A(j,1:j-1)*Y(1:j-1))/A(j,j);

end

Y

%

A=A'

for k=N:-1:1

X(k)=(Y(k)-A(k,k+1:N)*X(k+1:N))/A(k,k);

end

Cholesky分解

如果矩阵X是对称正定的,则Cholesky分解将矩阵X分解成一个下三角矩阵和上三角矩阵的乘积。设上三角矩阵为R,则下三角矩阵为其转置,即X=R'R。MATLAB函数chol(X)用于对矩阵X进行Cholesky分解,其调用格式为:

R=chol(X):产生一个上三角阵R,使R'R=X。若X为非对称正定,则输出一个出错信息。

  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值