dpq.m

来自「计算机潮流计算 应用的是牛拉法」· M 代码 · 共 20 行

function x=dPQ(Y,bus,n,m)
dP=zeros((n-1),1);
dQ=zeros(m,1);
for i=1:(n-1)
    sgmP=0;
    for j=1:n
   sgmP=sgmP+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)));
    end
   dP(i,1)=bus(i,4)-bus(i,2)*sgmP;
end
for i=1:m
    sgmQ=0;
    for j=1:n
   sgmQ=sgmQ+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)));
    end
   dQ(i,1)=bus(i,5)-bus(i,2)*sgmQ;
end
x=[dP;dQ];