資源簡(jiǎn)介
高等電力網(wǎng)絡(luò)節(jié)點(diǎn)導(dǎo)納矩陣LDU分解MATLAB實(shí)現(xiàn)代碼,該代碼簡(jiǎn)潔清晰,可以適用于任何形式的節(jié)點(diǎn)導(dǎo)納矩陣和節(jié)點(diǎn)阻抗矩陣的LDU分解,為下一步的前代回代求解做好數(shù)學(xué)準(zhǔn)備。
代碼片段和文件信息
function?[LDU]?=?calldu(?a?)
%對(duì)稱矩陣LDU因子分解
b=size(a);
m=b(1);
n=b(2);
L=eye(n);
D=eye(n);
U=eye(n);
for?j=1:n;
????for?i=1:j-1
????????U(ij)=a(ij);
????????for?k=1:i-1
???????????U(ij)=U(ij)-D(kk)*L(ik)*U(kj);
????????end
????????U(ij)=U(ij)/D(ii);
?
評(píng)論
共有 條評(píng)論