computing.m

来自「前推潮流算法」· M 代码 · 共 29 行

% end to root,U to P and Q
n=11;
for t=n:-1:1
    if t==1 k=0,j1=2,j2=9,qiantui1;
    elseif t==2 k=1,j1=7,j2=3,qiantui1;
    elseif t==3 k=2,j1=4,j2=0,qiantui3;
    elseif t==4 k=3,j1=5,j2=6,qiantui1;
    elseif t==5 k=4,j1=0,j2=0,qiantui2;
    elseif t==6 k=4,j1=0,j2=0,qiantui2;
    elseif t==7 k=2,j1=8,j2=0,qiantui3;
    elseif t==8 k=7,j1=0,j2=0,qiantui2;
    elseif t==9 k=1,j1=10,j2=0,qiantui3;
    elseif t==10 k=9,j1=11,j2=0,qiantui3;
    elseif t==11k=10,j1=0,j2=0,qiantui2; 
    else break;
    end
end
% root to end, P and Q to U
m=11;
U1=U;
for b=1:m
    if b==1 
    Iab(b)=(Pkt(b)-i*Qkt(b))/U0,
    U(b)=U0-Iab(b)*(r(b)+i*x(b))/1000;
    else
    Iab(b)=(Pkt(b)-i*Qkt(b))/conj(U(b-1)),
    U(b)=U(b-1)-Iab(b)*(r(b)+i*x(b))/1000;
    end
end