make_in_q.m

电力系统考虑网损的经济调度程序

function [dQ]=make_in_Q(A,G,B)
[k1,k2]=size(A);
[p1,p2]=size(G);
dQ=[];
for m=1:k1
    temp_q=0;
    for n=1:k1
        temp_q=temp_q+A(m,5)*A(n,5)*(G(m,n)*sin(A(m,6)-A(n,6))-B(m,n)*cos(A(m,6)-A(n,6)));
    end
    dQ=[dQ;temp_q];
end
for p=1:k1
  dQ(p,1)=A(p,4)-dQ(p,1);
end
dQ;
return;