MATLAB 中的约束最小化

Question

我想解决一个受约束的最小化问题，我正在寻求一些有关如何构建代码的帮助。

我知道fmincon 是我应该通过使用@mycon 参数来使用的，但我正在努力使其适应我的情况。任何建议将不胜感激。

这些是我的文件（a 和 b 是预定义的参数）：

1. f1.m

   function [y1, y2, y3]=f1(x1, x2, a)
   ...
   end
   

2. f2.m

   function w1=f2(x1, x2, y2, y3, b)
   ...
   end
   

我要编码的问题：

分钟 y1 w.r.t x1, x2
这样y1<=w1

Answer 1

您可以按如下方式使用fmincon：

x = fmincon(@(x) f1(x(1), x(2), a), [x1_start x2_start], [], [], [], [], [], [], @(x) mycon(x(1), x(2), a, b));
x1 = x(1)
x2 = x(2)

mycon 定义为：

% C <= 0 and Ceq == 0
function [C, Ceq] = mycon(x1, x2, a, b)
    [y1, y2, y3] = f1(x1, x2, a);
    C = y1 - f2(x1, x2, y2, y3, b);
    Ceq = [];
end

和x1_start 和x2_start 在Matlab 使用的迭代求解器中x1 和x2 的起始值。

请查看 matlab 文档和示例以获取更多信息。