資源簡介
對正定系數矩陣的特殊求解方程方法。最終結果為矩陣L,并由它以及它的轉置矩陣可以求解方程
代碼片段和文件信息
function?[?lyx?]?=?untitled5(?ab?)
y=[];x=[];
n=size(a1);
for?j=1:n
????c=0;
????for?k=1:j-1
????????c=c+l(jk)^2;
????end
????l(jj)=sqrt(a(jj)-c);
????for?i=j+1:n
????????if?j==1
????????????l(ij)=a(ij)/l(jj);
????????else
????????????c=0;
????????????for?k=1:j-1
????????????????c=c+l(ik)*l(jk);
????????????end
????????????l(ij)=(a(ij)-c)/l(jj);
????????end
????
評論
共有 條評論