矩阵分解LDLfunction [l,d]=ldlt(a)
%by swell2005
n=length(a(:,1));
l(1,1)=1;d(1,1)=a(1,1);
l(2,1)=a(2,1)/d(1,1);
l(2,2)=1;d(2,2)=a(2,2)-l(2,1)*d(1,1)*l(2,1);
for i=3:n,
l(i,1)=a(i,1)/d(1,1);
l(i,i)=1;
for j=2:i-1,
l(i,j)=(a(i,j)-l(i,1:j-1)*d(1:j-...
function [l,d]=ldlt(a)
%by swell2005
n=length(a(:,1));
l(1,1)=1;d(1,1)=a(1,1);
l(2,1)=a(2,1)/d(1,1);
l(2,2)=1;d(2,2)=a(2,2)-l(2,1)*d(1,1)*l(2,1);
for i=3:n,
l(i,1)=a(i,1)/d(1,1);
l(i,i)=1;
for j=2:i-1,
l(i,j)=(a(i,j)-l(i,1:j-1)*d(1:j-1,1:j-1)*l(j,1:j-1)')/d(j,j);
d(i,i)=a(i,i)-l(i,1:i-1)*d(1:i-1,1:i-1)*l(i,1:i-1)';
end
end
本文档为【矩阵分解LDL】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑,
图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。