www.pudn.com > OPF.zip > Imbalance.m, change:2013-03-31,size:2127b


%% *************************************************************** 
%             作者:陆镛 
%             原创日期:2012年11月29号 
%             修改日期:2012年12月2号 
%             函数说明:形成各变量的不平衡量。 
%% *************************************************************** 
function [Ly,Lz,Lw,Lxx,Llu,Luu,Lx] = Imbalance(miu,NodeNum,GenNum,GenI,GenA,GenB,... 
    GenPmin,GenPmax,PvNum,PvI,PvQmin,PvQmax,NodePl,NodeQl,NodePg,NodeQg,... 
    EqualityJmatrix,InequalityJmatrix,Ycos,Ysin,Pg,Qr,AmplitudeU,... 
    l,u,w,z,y,CapacityRef) 
%% 变量说明:****************************************************** 
% 输入参数说明:miu:对偶间隙  NodeNum:节点数   GenNum:发电机节点数  GenI:发电机节点号   
% GenA:耗量特性二次项系数  GenB:耗量特性一次项系数  GenPmin:有功出力下限  GenPmax:有功出力上限 
% PvNum:Pv节点数  PvI:PV节点i  PvQmin:无功出力下限  PvQmax:无功出力上限  
% NodePg:有功出力  NodeQg:无功出力  NodePl:有功负荷  NodeQl:无功负荷 
% EqualityJmatrix:等式约束矩阵  InequalityJmatrix:不等式约束矩阵  CapacityRef:基准容量 
% Ycos:导纳矩阵幅值点乘AngleUij减去导纳矩阵相角的余弦值  Ysin:导纳矩阵幅值点乘AngleUij减去导纳矩阵相角的正弦值  
% Pg:有功出力  Qr:无功出力  AmplitudeU:电压幅值  l、u:松弛变量  w、z、y:拉格朗日乘子 
% 输出参数说明:Ly、Lz、Lw、Lxx、Llu、Luu 
%% *************************************************************** 
DiagU=sparse(1:NodeNum,1:NodeNum,AmplitudeU,NodeNum,NodeNum); 
NodePg(GenI)=0; 
NodeQg(PvI)=0; 
ImbalanceP=DiagU*Ycos*AmplitudeU-sparse(GenI,1,Pg,NodeNum,1)-... 
    (NodePg-NodePl)/CapacityRef;  
ImbalanceQ=DiagU*Ysin*AmplitudeU-sparse(PvI,1,Qr,NodeNum,1)-... 
    (NodeQg-NodeQl)/CapacityRef;  
%% 
Lx=sparse(1:GenNum,ones(GenNum,1),2*GenA.*Pg+GenB,GenNum+PvNum+2*NodeNum... 
    ,1)-EqualityJmatrix*y-InequalityJmatrix*(z+w); 
Ly=[ImbalanceP;ImbalanceQ]; 
Lz=[Pg;Qr;AmplitudeU]-l-[GenPmin/CapacityRef;PvQmin/CapacityRef;0.9*... 
    ones(NodeNum,1)]; 
Lw=[Pg;Qr;AmplitudeU]+u-[GenPmax/CapacityRef;PvQmax/CapacityRef;1.1*... 
    ones(NodeNum,1)]; 
Llu=l.*z-miu; 
Luu=u.*w+miu; 
%% 
Lxx=Lx+InequalityJmatrix*((1./l).*(Llu+z.*Lz)+(1./u).*(Luu-w.*Lw)); 
end