高等电力网络节点导纳矩阵LDU分解MATLAB实现代码

function [LDU] = calldu（ a ）
%对称矩阵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）;