options.auxdata = struct( ...
'om', om, ...
'Ybus', Ybus, ...
'Yf', Yf(il,:), ...
'Yt', Yt(il,:), ...
'mpopt', mpopt, ...
'il', il, ...
'A', A, ...
'nA', nA, ...
'neqnln', neqnlnnumber, ...
'niqnln', niqnlnnumber, ...
'Js', Js, ...
'Hs', Hs );
%% define variable and constraint bounds
options.lb = xmin;
options.ub = xmax;
if gTOTAL_Scenarios==1
options.cl = [zeros(2*nb*gTOTAL_Scenarios, 1); -Inf*ones(2*nl2*gTOTAL_Scenarios, 1); l];
options.cu = [zeros(2*nb*gTOTAL_Scenarios, 1); zeros(2*nl2*gTOTAL_Scenarios, 1); u];
else
options.cl = [zeros(2*nb*gTOTAL_Scenarios, 1); -Inf*ones(2*nl2*gTOTAL_Scenarios+2, 1); l];
options.cu = [zeros(2*nb*gTOTAL_Scenarios, 1); zeros(2*nl2*gTOTAL_Scenarios+2, 1); u];
end
%% assign function handles
%需要提供目标函数,目标函数梯度,约束,约束雅克比矩阵,约束海森矩阵