get_pq.m

来自「电力系统计算节点数不多的潮流程序,数度比较快」· M 代码 · 共 24 行

function [P,Q]=get_PQ(V,ang,Y,nb,m,bus)
Q(:,1)=bus(:,5);
P(:,1)=bus(:,4);
for p=1:nb,
    for q=1:nb,
        G(p,q)=real(Y(p,q));
        B(p,q)=imag(Y(p,q));
        angl(p,q)=ang(p,1)-ang(q,1);
    end
end
for p=m+1:nb-1,
    Q(p,1)=0;
    for q=1:nb,
        Q(p,1)=Q(p,1)+V(p,1)*V(q,1)*(G(p,q)*sin(angl(p,q))-B(p,q)*cos(angl(p,q)));
    end
end
P(nb,1)=0;
Q(nb,1)=0;
for q=1:nb,
    P(nb,1)=P(nb,1)+V(nb,1)*V(q,1)*(G(nb,q)*cos(angl(nb,q))+B(nb,q)*sin(angl(nb,q)));
    Q(nb,1)=Q(nb,1)+V(nb,1)*V(q,1)*(G(nb,q)*sin(angl(nb,q))-B(nb,q)*cos(angl(nb,q)));
end