Matlab结合fmincon与ODE求解x3最大值的优化方法咨询

核心逻辑说明 #

fmincon默认求解最小值，我们只需要将目标函数设置为-x3，对该值求最小值即可等价于求x3的最大值，同时配置好参数边界约束和x3的区间约束即可完成优化。

具体实现步骤 #

1. 封装优化目标函数 #

将已有ODE求解逻辑封装为可被优化器调用的函数，输入为待优化的4个参数组成的向量，输出为目标函数值、非线性约束值。
示例代码如下，你可以根据自己的Scrpt1函数和求解逻辑修改对应部分：

function [f, c, ceq] = opt_obj(params)
    % 参数映射：params(1)=kst, params(2)=x1, params(3)=x5, params(4)=xo
    kst = params(1);
    x1 = params(2);
    x5 = params(3);
    xo = params(4);
    
    % 调用你已有的ode15s求解逻辑，替换成你自己的tspan、初始值y0等参数
    [t,y] = ode15s(@Scrpt1, tspan, y0, [], kst, x1, x5, xo);
    
    % 提取x3值，这里假设x3是求解结果y的第3列，根据你自己的变量顺序调整
    x3 = y(:,3);
    
    % 配置x3的区间约束：要求x3处于[-8e-4, 2e-4]，非线性不等式约束要求c<=0
    c(1) = -8e-4 - x3; % 等价于x3 >= -8e-4
    c(2) = x3 - 2e-4; % 等价于x3 <= 2e-4
    ceq = []; % 无等式约束时留空
    
    % 目标函数：如果要取最终时刻x3的最大值，就取x3(end)，要取整个时域x3最大值就替换为max(x3)
    f = -x3(end);
end

2. 配置优化参数 #

% 4个参数的上下界，顺序对应kst、x1、x5、xo
lb = [0.1, 1e-9, 4, 4];
ub = [2, 1e-6, 15, 10];

% 初始值可以选参数区间的中间值，也可以自己设置合理值
x0 = [1, 5e-7, 9.5, 7];

% 优化选项配置，选择sqp算法适配带ODE的非线性优化，打开迭代显示方便查看过程
options = optimoptions('fmincon','Display','iter','Algorithm','sqp');

3. 调用fmincon求解 #

[x_opt, fval, exitflag] = fmincon(@opt_obj, x0, [], [], [], [], lb, ub, [], options);

% 解析最优结果
kst_opt = x_opt(1);
x1_opt = x_opt(2);
x5_opt = x_opt(3);
xo_opt = x_opt(4);
x3_max = -fval; % 转换得到x3的最大值

替代方案 #

如果模型存在较多局部最优导致fmincon收敛到局部解，可以改用Matlab全局优化工具箱的patternsearch（模式搜索）或者ga（遗传算法），调用逻辑和上述配置一致，仅替换求解器函数即可。

注意事项 #

• 确保你的Scrpt1函数支持传入kst、x1、x5、xo四个参数，不要硬写参数值
• 如果只要求最终时刻的x3满足区间约束，把约束里的x3替换为x3(end)即可，减少约束计算量
• 出现不收敛情况时可以调整初始值，或者适当修改优化选项的终止容差参数

内容的提问来源于stack exchange，提问作者user16739361