fmincon中目标函数和约束函数的单次计算

Question

在我的程序中，我需要解决优化程序中的边界值问题，该程序也具有非线性约束。使用fmincon 解决问题，我需要BVP 的解决方案来评估目标函数和非线性约束函数。

目前我正在解决目标函数和约束函数中的 BVP。有没有一种更有效的方法，在目标函数中对 BVP 进行一次评估后，我可以将解决方案传递给约束函数，以减少一次 BVP 评估。任何想法

Answer 1

创建一个带有额外输入的函数，用于区分两个调用。这基本上是目标和约束函数的包装，将所有必要的中间输出保存在persistent 变量中：

function varargout = my_obj_con_function(..., option)

    persistent C, Ceq; % ... or any other data you might need

    switch (option)

        case 'obj'

            % objective function 
            [varargout{1}, data] = objfun(...); % 'data' is an example

            % constraint function 
            [C, Ceq] = confcn(data, ...);

        case 'con'

            % You've just computed this -- just return it
            varargout{1} = C;
            varargout{2} = Ceq;
            return 

    end

end

这样使用：

[...] = fmincon(@(x) my_obj_con_function(x, ..., 'obj'), ...
                x0,...
                Lb, Ub, ...
                A, b,...
                Aeq, beq, ...
                @(x) my_obj_con_function(x, ..., 'con'), ...);

这明确假设fmincon 总是调用约束函数在调用目标函数之后，完全相同决策值相同变量。现在我不确定这是否可以保证——您可能希望在不能保证的情况下建立一些保护和机制。