node.m

来自「电力系统计算节点数不多的潮流程序,数度比较快」· M 代码 · 共 55 行

function [S,SS,delS]=node(V,line,ang,nb)
S=zeros(nb);
SS=zeros(nb);
delS=zeros(nb);
[nl,ml]=size(line);
for k=1:nl,
    p=line(k,1);
    q=line(k,2);
    Zt=line(k,3)+j*line(k,4);
    Yt=1/Zt;
    Ym=line(k,5)+j*line(k,6);
    K=line(k,7);

    
        if (K==0)&(q~=0)
               
                 I(p,q)=V(p,1)*(cos(ang(p,1))+j*sin(ang(p,1)))*(Yt+Ym)-V(q,1)*(cos(ang(q,1))+j*sin(ang(q,1)))*Yt;
                 I(q,p)=V(q,1)*(cos(ang(q,1))+j*sin(ang(q,1)))*(Yt+Ym)-V(p,1)*(cos(ang(p,1))+j*sin(ang(p,1)))*Yt;
                 S(p,q)=V(p,1)*(cos(ang(p,1))+j*sin(ang(p,1)))*conj(I(p,q));
                 SS(q,p)=V(q,1)*(cos(ang(q,1))+j*sin(ang(q,1)))*conj(I(q,p));
                 delS(p,q)=S(p,q)+SS(q,p);
             
         end 
             if K>0
                 
                  I(p,q)=V(p,1)*(cos(ang(p,1))+j*sin(ang(p,1)))*(Yt/K+Ym+Yt*(K-1)/K)-V(q,1)*(cos(ang(q,1))+j*sin(ang(q,1)))*Yt/K;
                  I(q,p)=V(q,1)*(cos(ang(q,1))+j*sin(ang(q,1)))*(Yt/K+Yt*(1-K)/K/K)-V(p,1)*(cos(ang(p,1))+j*sin(ang(p,1)))*Yt/K;
                   S(p,q)=V(p,1)*(cos(ang(p,1))+j*sin(ang(p,1)))*conj(I(p,q));
                  SS(q,p)=V(q,1)*(cos(ang(q,1))+j*sin(ang(q,1)))*conj(I(q,p));
                  delS(p,q)=S(p,q)+SS(q,p);
              
              end
        
            if K<0
                K=-K;
                  I(p,q)=V(p,1)*(cos(ang(p,1))+j*sin(ang(p,1)))*(Yt*K+Ym+Yt*(1-K))-V(q,1)*(cos(ang(q,1))+j*sin(ang(q,1)))*Yt*K;
                  I(q,p)=V(q,1)*(cos(ang(q,1))+j*sin(ang(q,1)))*(Yt*K+Yt*(K*K-K))-V(p,1)*(cos(ang(p,1))+j*sin(ang(p,1)))*Yt*K;
                  S(p,q)=V(p,1)*(cos(ang(p,1))+j*sin(ang(p,1)))*conj(I(p,q));
                  SS(q,p)=V(q,1)*(cos(ang(q,1))+j*sin(ang(q,1)))*conj(I(q,p));
                  delS(p,q)=S(p,q)+SS(q,p);
              
              end
        
            if (K==0)&(q==0)
                 q=p;
                   S(p,q)=V(p,1)*V(p,1)*conj(Ym);
                   SS(q,p)=-S(p,q);
                   delS(p,q)=S(p,q)+SS(q,p);
               end
           end