dpq.m

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

function dPQ=dPQ(Y,bus)
%delP--有功偏差量
%delQ--无功偏差量
%Y--节点导纳矩阵
%bus--节点参数（P,Q,U及相角）矩阵
global n;
global m;
delP=zeros(1,n-1);
delQ=zeros(1,m);
for i=1:(n-1)  %形成delP矩阵(PQ、PV节点)
    P1=0;T=0;
    for j=1:n
        T=bus(i,2)*bus(j,2)*(real(Y(i,j))*cos(bus(i,3)-bus(j,3))+imag(Y(i,j))*sin(bus(i,3)-bus(j,3)));
        P1=P1+T;
    end
    delP(1,i)=bus(i,4)-P1;
end
for i=1:m %形成delQ矩阵(PQ节点)
    Q1=0;T=0;
    for j=1:n
        T=bus(i,2)*bus(j,2)*(real(Y(i,j))*sin(bus(i,3)-bus(j,3))-imag(Y(i,j))*cos(bus(i,3)-bus(j,3)));
        Q1=Q1+T;
    end
    delQ(1,i)=bus(i,5)-Q1;
end
dPQ=[delP,delQ];