make_in_p.m

来自「电力系统考虑网损的经济调度程序」· M 代码 · 共 16 行

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