dpq.m

电力系统计算节点数不多的潮流程序,数度比较快

function delPQ=dPQ(V,ang,Y,bus,nb,m)

for i=1:nb,
    for j=1:nb,
        G(i,j)=real(Y(i,j));
        B(i,j)=imag(Y(i,j));
        angl(i,j)=ang(i,1)-ang(j,1);
    end
end

for i=1:nb-1,
    x=0;
    for j=1:nb,
        x=x+V(i,1)*V(j,1)*(G(i,j)*cos(angl(i,j))+B(i,j)*sin(angl(i,j)));
    end
    delP(i,1)=bus(i,4)-x;
end

for i=1:m,
    y=0;
  
    for j=1:nb,
        y=y+V(i,1)*V(j,1)*(G(i,j)*sin(angl(i,j))-B(i,j)*cos(angl(i,j)));
    end
    delQ(i,1)=bus(i,5)-y;
end
        
 
delPQ=[delP;delQ];