题号：72编程计算潮流（电力系统分析）节点数：3 支路数：3 计算精度：0.00010支路1：0.0

网友回答

n=5;nl=5;isb=1;
pr=0.00001; 
B1=[1 2 0.03i 0 1.05 0;2 3 0.08+0.3i 0.5i 1 0;2 4 0.1+0.35i 0 1 0;3 4 0.04+0.25i 0.5i 1 0;3 5 0.015i 0 1.05 1];
B2=[0 0 1.05 1.05 0 1;0 3.7+1.3i 1.05 0 0 2;0 2+1i 1.05 0 0 2;0 1.6+0.8i 1.05 0 0 2;5 0 1.05 1.05 0 3];
X=[1 0;2 0;3 0;4 0;5 0];
na=3;Y=zeros(n);YI=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);O=zeros(1,n);
for i=1:n
    if X(i,2)~=0;
        p=X(i,1);
        Y(p,p)=1./X(i,2);
    end
endfor i=1:n
l    if B1(i,6)==0
        p=B1(i,1);q=B1(i,2);
    else p=B1(i,2);q=B1(i,1);
    end
    Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));
    YI(p,q)=YI(p,q)-1./B1(i,3);
    Y(q,p)=Y(p,q);
    YI(q,p)=YI(p,q);
    Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;
    YI(q,q)=YI(q,q)+1./B1(i,3);
    Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;
    YI(p,p)=YI(p,p)+1./B1(i,3);
endG=real(Y);B=imag(YI);BI=imag(Y);for i=1:n
    S(i)=B2(i,1)-B2(i,2);    BI(i,i)=BI(i,i)+B2(i,5);endP=real(S);Q=imag(S);for i=1:n
    e(i)=real(B2(i,3));    f(i)=imag(B2(i,3));    V(i)=B2(i,4);endfor i=1:n
    if B2(i,6)==2        V(i)=sqrt(e(i)^2+f(i)^2);        O(i)=atan(f(i)./e(i));    end
endfor i=2:n    if i==n&nbs