dpq.m

来自「这个程序是采用模块化的方式对潮流进行求解」· M 代码 · 共 27 行

function [delP,delQ] =dPQ(Y,bus)
global   G B V delta P Q m  n
delP = zeros(n-1,1);
delQ = zeros(m,1);

    PPP=0;
    for i=1:n-1
        PPP=0;
        for j=1:n
        PPP=PPP+V(i)*V(j)*(G(i,j)*cos(delta(i)-delta(j))+B(i,j)*sin(delta(i)-delta(j)));   
        end
        delP(i)=P(i)-PPP;      % 计算有功不平衡量 
    end    
    

    
   
    for i=1:m
        PPP=0;
        for j=1:n
         PPP= PPP+V(i)*V(j)*(G(i,j)*sin(delta(i)-delta(j))-B(i,j)*cos(delta(i)-delta(j)));
        end
        delQ(i)=Q(i)-PPP;       %  计算无功不平衡量
    end