function dPQ=dPQ(Y,bus)
%delP--有功偏差量
%delQ--無功偏差量
%Y--節(jié)點(diǎn)導(dǎo)納矩陣
%bus--節(jié)點(diǎn)參數(shù)（P,Q,U及相角）矩陣
global n;
global m;
delP=zeros(1,n-1);
delQ=zeros(1,m);
for i=1:(n-1)  %形成delP矩陣(PQ、PV節(jié)點(diǎn))
    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節(jié)點(diǎn))
    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];