loss.m

利用该MATLAB程序可实现对一般电力系统的节点、支路及其损耗的潮流计算。

function Pf=loss(bus,line)
%计算线路损耗、线路潮流
[nl,ml]=size(line);
Pf=zeros(nl,5);
for k=1:nl
    I=line(k,1);
    J=line(k,2);
    Zt=line(k,3)+i*line(k,4);
     if Zt==0
        Yt=inf;
     else
        Yt=1/Zt;
     end
    Ym=line(k,5)+i*line(k,6);
    K=line(k,7);
    
    if (K==0)&(J~=0)   %普通线路潮流       
       S(I,J)=bus(I,2)^2*(conj(Yt)+conj(Ym))-bus(I,2)*(cos(bus(I,3))+i*sin(bus(I,3)))*bus(J,2)*(cos(bus(J,3))-i*sin(bus(J,3)))*conj(Yt);
       S(J,I)=bus(J,2)^2*(conj(Yt)+conj(Ym))-bus(J,2)*(cos(bus(J,3))+i*sin(bus(J,3)))*bus(I,2)*(cos(bus(I,3))-i*sin(bus(I,3)))*conj(Yt);
       delS(I,J)=S(I,J)+S(J,I);       
    end
    
    if(K==0)&(J==0)  %对地支路潮流
        J=5;
        S(I,5)=bus(I,2)^2*conj(Ym);
    end
    
    if K>0   %变压器支路k>0时的潮流
       S(I,J)=bus(I,2)^2*(conj(Ym+Yt*(1-1/K))+conj(Yt/K))-bus(I,2)*(cos(bus(I,3))+i*sin(bus(I,3)))*bus(J,2)*(cos(bus(J,3))-i*sin(bus(J,3)))*conj(Yt/K);
       S(J,I)=bus(J,2)^2*(conj(Yt))/K^2-bus(J,2)*(cos(bus(J,3))+i*sin(bus(J,3)))*bus(I,2)*(cos(bus(I,3))-i*sin(bus(I,3)))*conj(Yt/K);
       delS(I,J)=S(I,J)+S(J,I);    
    end
    
    if K<0    %变压器支路k<0时的潮流
       S(I,J)=bus(I,2)^2*(conj(Ym+Yt))+bus(I,2)*(cos(bus(I,3))+i*sin(bus(I,3)))*bus(J,2)*(cos(bus(J,3))-i*sin(bus(J,3)))*conj(Yt*K);
       S(J,I)=bus(J,2)^2*(conj(Yt))*K^2+bus(J,2)*(cos(bus(J,3))+i*sin(bus(J,3)))*bus(I,2)*(cos(bus(I,3))-i*sin(bus(I,3)))*conj(Yt*K);
       delS(I,J)=S(I,J)+S(J,I);       
    end
      if  J==5&Zt==0
          Sp=[line(k,1) line(k,2) S(I,5) 0  S(I,5)];      
      else
       Sp=[line(k,1) line(k,2) S(I,J) S(J,I) delS(I,J)];
      end
      Pf(k,:)=Sp;
end