ipopt 使用总结

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