powerlosses.m

这个程序是采用模块化的方式对潮流进行求解

function PowerLosses()
global bus line n m nodenum Y V delta myf;

[nl,ml]=size(line);
S_ij=zeros(nl,1);
S_ji=zeros(nl,1);
dS_ij=zeros(nl,1);

for t=1:nl
    Zt=line(t,3)+j*line(t,4);
    Yt=1/Zt;
    Ym=line(t,5)+j*line(t,6);
    I=line(t,1);
    J=line(t,2);
    K=line(t,7);
    
    if (K==0)&(J~=0)                  % 线路为普通线路 且不是接地端
        S_ij(t)=V(I)^2*conj(Yt+Ym)-V(I)*exp(j*delta(I))*conj(V(J)*exp(j*delta(J)))*conj(Yt);
        S_ji(t)=V(J)^2*conj(Yt+Ym)-V(J)*exp(j*delta(J))*conj(V(I)*exp(j*delta(I)))*conj(Yt);
    end
    
     if (K==0)&(J==0)                 % 线路为接地支路
         S_ij(t)=V(I)^2*conj(Ym);
     end
     
     if K>0                           % 变压器线路: Zt和Ym为折算到i侧的值,K在j侧
        S_ij(t) = V(I)^2 * conj(Yt+Ym) - V(I)*exp(j*delta(I)) * conj(V(J)*exp(j*delta(J)))*conj(Yt/K);
        S_ji(t) = V(J)^2 * conj(Yt/K/K) -V(J)*exp(j*delta(J)) * conj(V(I)*exp(j*delta(I)))*conj(Yt/K);
    end
    
     if K<0                           % 变压器线路: Zt和Ym为折算到K侧的值,K在i侧
        S_ij(t) = V(I)^2 * conj(Yt+Ym) + V(I)*exp(j*delta(I)) * conj(V(J)*exp(j*delta(J)))*conj(Yt*K);
        S_ji(t) = V(J)^2 * conj(Yt*K*K) + V(J)*exp(j*delta(J)) * conj(V(I)*exp(j*delta(I)))*conj(Yt*K);% 到底是加还是减？？？
    end
end
         
dS_ij = S_ij + S_ji;
  
fprintf(myf,'\n');
fprintf(myf,'线路计算结果：\n');
fprintf(myf,'节点I    节点J                  线路功率S(I,J)                   线路功率S(J,I)                  线路损耗dS(I,J)\n');
 
for t = 1:nl
    I=line(t,1);                      
    J=line(t,2);
    fprintf(myf,' %4d     %4d        %12.6f+j%8.6f        %12.6f+j%12.6f         %12.6f+j%12.6f\n',renode(I),renode(J),real(S_ij(t)),imag(S_ij(t)),real(S_ji(t)),imag(S_ji(t)),real(dS_ij(t)),imag(dS_ij(t)));
end