www.pudn.com > OPF.zip > Initialization.m, change:2013-01-24,size:2179b


%% *************************************************************** 
%             作者:陆镛 
%             原创日期:2012年11月24号 
%             修改日期:2012年12月1号 
%             函数说明:松弛变量、拉格朗日乘子、各变量初始化。 
%% *************************************************************** 
function [l,u,z,w,y,Pg,Qr,AmplitudeU,AngleU] = Initialization(NodeNum,PvNum... 
    ,PvQmin,PvQmax,GenNum,GenPmin,GenPmax,CapacityRef) 
%% 变量说明:****************************************************** 
% 输入参数说明: NodeNum:节点数   PvNum:PV节点数  PvQmin:无功出力下限  PvQmax:无功出力上限  
% GenNum:发电机节点数  GenPmin:有功出力下限  GenPmax:有功出力上限  CapacityRef:基准容量 
% 输出参数说明: l、u:松弛变量  z、w、y:拉格朗日乘子  Pg:有功出力  Qr:无功出力 
% AngleU:电压相角   AmplitudeU:电压幅值 
%% *************************************************************** 
l=ones(NodeNum+GenNum+PvNum,1);                                             %设置松弛变量l初值1。 
u=ones(NodeNum+GenNum+PvNum,1);                                             %设置松弛变量u初值1。 
%% *************************************************************** 
z=ones(NodeNum+GenNum+PvNum,1);                                             %设置拉格朗日乘子z初值1。 
w=-ones(NodeNum+GenNum+PvNum,1);                                            %设置拉格朗日乘子w初值-1。 
% y=[1.e-10*ones(NodeNum,1);-1.e-10*ones(NodeNum,1)];                         %设置拉格朗日乘子y初值。 
y=sparse(2*NodeNum,1); 
%% *************************************************************** 
Pg=(GenPmin+GenPmax)/2/CapacityRef;                                         %设置有功出力初值,取上限和下限的平均值。 
Qr=(PvQmin+PvQmax)/2/CapacityRef;                                           %设置无功出力初值,取上限和下限的平均值。 
%% ***************************************************************                                                                                                              
AmplitudeU=ones(NodeNum,1);                                                 %设置幅值初值1。        
AngleU=sparse(NodeNum,1);                                                   %设置相角初值0。 
end