Cholesky分解（Matlab代码实现）

 👨‍🎓个人主页：研学社的博客 

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

​

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

💥1 概述

Cholesky分解法又叫平方根法，是求解对称正定线性方程组最常用的方法之一。对于一般矩阵，为了消除LU分解的局限性和误差的过分积累，采用了选主元的方法，但对于对称正定矩阵而言，选主元是不必要的。

更详细知识点，这个博主讲解很清晰，超棒：

矩阵Cholesky分解，及其在求解线性方程组、矩阵逆的应用

📚2 运行结果

 

部分代码：

for i=1:N
    W3(:,i)=randn(2,1); 
end
for i=1:N
    WW3(:,i)=wblrnd(9.0,2.15,2,1); 
end

Z3=R3*W3;
ccc03=corrcoef(W3');
ccc03
ccc13=corrcoef(Z3');
ccc13

aa3=zeros(2,N);
Ls3=zeros(2,N);
for p=1:2
hig=max(Z3(p,:));
k=1;
for i=1:N
[b c]=min(Z3(p,:));
Ls3(p,c)=k;
k=k+1;
Z3(p,c)=hig+1;
end
end
LLs3=Ls3;

x=zeros(2,N);
for i=1:N
    tt=(i-0.5)/N;
    x3(1,i)=wblinv(tt,9.0,2.15);
    x3(2,i)=wblinv(tt,9.0,2.15);   
end

for i=1:2
y3=x3(i,:);
bb3=LLs3(i,:);
hig1=max(y3);
hig2=max(bb3);

for p=1:N
[C1 D1]=min(y3);
[C2 D2]=min(LLs3(i,:));
bb3(1,D2)=C1;
y3(1,D1)=hig1+1;
LLs3(i,D2)=hig2+1;
end
yy3(i,:)=bb3;
end

🎉3 参考文献

部分理论来源于网络，如有侵权请联系删除。

[1]李辉,郝如江.基于不完全Cholesky分解相关熵双谱的轴承故障诊断[J].振动与冲击,2022,41(11):123-132

[2]笪涵,胡圣波.基于Cholesky矩阵分解的贝叶斯压缩感知信号处理[J].贵州师范大学学报：自然科学版,2021,39(1):72-76

🌈4 Matlab代码实现