想用 MATLAB 解决优化问题最快的方法是直接使用官方优化算法工具箱Optimization Toolbox。它解决的不是某一个具体问题而是一整类问题线性规划、带线性或非线性约束的优化、无约束优化、二次规划、最小二乘数据拟合、非线性方程组求解全部封装成统一函数接口。你不需要从零写迭代算法不需要手搓梯度下降甚至不需要自己判断该用哪种数值方法只需要把目标函数、约束、初始点按 MATLAB 的规则传进去求解器就会自动处理算法选择、迭代终止、数值稳定性这些细节。从工程使用角度来看这个工具箱最值得关注的几个点第一覆盖面广从一元函数最小值到多变量强非线性约束问题都有对应求解器第二接口风格统一所有求解器都遵循“目标函数 约束 选项”的调用形式学会一个就能快速迁移到其他问题第三扩展性强可以结合遗传算法、模拟退火、模式搜索等全局优化方法处理多峰、非光滑、离散变量问题也能通过 for/parfor 批量跑大量初值和多组数据。下面我会从环境检查、工具箱启用、常用求解器、批量任务、性能调优到排错给出一套可以照跑的完整流程。上一期提到的遗传算法工具箱和这里的优化工具箱是同一条技术体系里的不同层次建议结合起来看。1. 核心能力速览能力项说明项目类型MathWorks 官方 MATLAB 优化工具箱主要功能线性规划、二次规划、无约束/约束非线性优化、最小二乘拟合、非线性方程组求解、多目标优化常用求解器fminunc、fmincon、linprog、intlinprog、quadprog、fminbnd、fsolve、lsqcurvefit、lsqnonlin全局优化求解器ga、patternsearch、simulannealbnd、multistart、globalsearch需配合全局优化能力使用调用方式命令行函数调用、.m 脚本调用、实时编辑器 Optimize 任务、传统 App 界面是否支持批量支持可结合 for/parfor 批量求解多个初值、多组数据是否支持并行支持可配合 Parallel Computing Toolbox 使用硬件要求普通办公电脑即可运行大规模问题建议配置充足内存支持平台Windows、Linux、macOS 上的 MATLAB许可要求商业授权必须使用正版许可这里需要单独说明一下版本差异。在新版本 MATLAB 中全局优化工具箱和优化工具箱已经合并ga、patternsearch 这些求解器直接通过优化工具箱提供。如果你使用的是老版本可能还需要单独安装 Global Optimization Toolbox 才能用遗传算法相关函数。所以后面代码里用到 ga 时我会明确标注这一点。2. 适用场景与使用边界这个工具箱适合解决以下类型的问题参数辨识与数据拟合比如根据实验数据反推模型参数对应 lsqcurvefit、lsqnonlin。工程设计优化比如结构参数、尺寸、控制参数寻优对应 fmincon、fminunc。资源调度与路径规划中的线性规划、整数规划对应 linprog、intlinprog。求解复杂的非线性方程组比如多元隐式方程、平衡点计算对应 fsolve。需要跳出局部最优的全局优化任务配合 ga、multistart、patternsearch 使用。但也要说清楚边界。这个工具箱不是万能建模软件它解决的问题类型有明确数学形式如果问题不是凸优化fminunc、fmincon 这类局部求解器只能保证找到局部最优不保证全局最优。如果变量是离散整数普通 fmincon 不能直接处理要用 intlinprog 或 ga。如果目标函数无法写成可计算的 MATLAB 函数比如完全依赖外部软件仿真那每次求解都需要调用外部程序效率会明显下降。如果问题本身没有可行解求解器会返回不可行状态这时候需要回头检查约束是否自相矛盾。最后是合规边界。这个工具箱是商业软件需要从 MathWorks 官方渠道获得授权许可。不要使用任何非官方激活方式容易出现函数缺失、许可证报错也存在安全隐患。如果你用 MATLAB 处理真实业务数据、人脸数据、内部图纸或受版权保护的素材要确认自己有权使用这些数据。特别是把优化结果发布到论文、报告或商用产品时务必核验数据的来源和授权链条。优化代码本身是一个工具工具怎么用、用在哪里由使用者负责。3. 环境准备与前置条件开始写代码之前先确认 MATLAB 环境是否已经具备优化工具箱。打开 MATLAB 命令窗口执行下面这几条命令。% 查看优化工具箱版本装了会显示版本号和描述信息 ver(optim)% 检测许可证是否可用返回 1 表示可用返回 0 表示不可用 license(test, Optimization_Toolbox)% 检查关键求解器是否在搜索路径上 which fmincon which linprog如果ver(optim)显示出版本信息说明工具箱已安装。如果license(test, Optimization_Toolbox)返回 0说明安装过但许可证没有覆盖当前用户需要联系 MATLAB 授权管理员开通。如果which fmincon返回空说明工具箱不在 MATLAB 搜索路径中。老版本环境还需要额外检查全局优化工具箱% 老版本中遗传算法等函数由全局优化工具箱提供 ver(globaloptim) license(test, Global_Optimization_Toolbox)对于一般的学习型和科研型优化问题普通办公电脑即可满足要求。复杂工程优化问题建议内存 16GB 以上并在求解时关注迭代日志和内存占用。MATLAB 本身支持 Windows、Linux、macOS 三大平台优化工具箱在这三个平台上的行为完全一致。还要处理好工作目录。建议所有优化代码放在独立文件夹里用cd切换过去或者把这个文件夹加入搜索路径。cd(D:\matlab_opt_project) addpath(D:\matlab_opt_project)4. 安装部署与工具箱启用如果你的 MATLAB 环境里没有优化工具箱要通过官方渠道补装。常见方式有两种。第一种是学校的校园授权或企业授权通常管理员会统一开通工具箱你在 MATLAB 里执行ver(optim)能直接看到。第二种是个人授权可以在 MATLAB 主页选项卡中打开“附加功能”管理器找到“管理附加功能”从 MathWorks 官网获取新增产品。也可以重新运行 MATLAB 安装程序在产品选择界面勾选 Optimization Toolbox。安装完以后需要重启 MATLAB 才能生效。重启后不要急着写复杂代码先跑一个最小示例验证工具箱可用性。% 一元函数最小值求解验证工具箱核心功能 fun (x) (x - 3).^2 2; [x, fval] fminbnd(fun, -10, 10); fprintf(最优解 x %.6f目标函数值 %.6f\n, x, fval);如果命令窗口输出类似“最优解 x 3.000000目标函数值 2.000000”的内容说明工具箱已经正常工作。这里必须强调整个过程只使用正版许可。有些破解版 MATLAB 表面能运行基础功能但到调用 fmincon、linprog 这类工具箱函数时会出现“Undefined function”或许可证弹窗很难排查。从官方渠道获取授权是最稳妥的路线。5. 功能测试与效果验证下面按照不同问题类型逐一跑通常用求解器。所有代码都可以直接复制到 MATLAB 中运行每个示例都会给出判断成功与否的标准。5.1 无约束优化fminuncfminunc 用于求解无约束多变量函数的最小值典型场景是参数估计、模型校准。这里用经典的 Rosenbrock 函数测试函数表达式为 f(x) 100×(x2 - x1²)² (1 - x1)²全局最小值在 (1, 1) 处目标函数值为 0。% 目标函数Rosenbrock 函数 fun (x) 100 * (x(2) - x(1)^2)^2 (1 - x(1))^2; % 初始点 x0 [-1.2; 1]; % 调用求解器 [x, fval, exitflag, output] fminunc(fun, x0); fprintf(x [%.4f, %.4f]\n, x(1), x(2)); fprintf(目标函数值 %.6f\n, fval); fprintf(退出标志 %d\n, exitflag); fprintf(迭代次数 %d\n, output.iterations);运行成功后x 应该接近 [1, 1]目标函数值接近 0退出标志为 1。如果 exitflag 为 0说明达到最大迭代次数或函数评估次数上限需要增大 MaxFunctionEvaluations 或 MaxIterations。fminunc 适用的前提是目标函数连续可导。如果目标函数不可导、不连续或者存在大量局部极值点fminunc 不一定能稳定求解这时要改用 ga 或 patternsearch。5.2 线性规划与整数规划linprog、intlinprog线性规划是实际工程中最常见的规划类型。linprog 求解形式为 min f×x满足线性不等式 A×x ≤ b、线性等式 Aeq×x beq、边界 lb ≤ x ≤ ub。下面是一个两变量线性规划示例。% 目标函数系数即最小化 -x1 - 2x2 f [-1; -2]; % 线性不等式约束 A*x b A [1, 1; -1, 1; 0, 1]; b [6; -2; 3]; % 无等式约束 Aeq []; beq []; % 变量下界和上界 lb [0; 0]; ub []; % 求解 [x, fval, exitflag] linprog(f, A, b, Aeq, beq, lb, ub); fprintf(x [%.4f, %.4f]\n, x(1), x(2)); fprintf(目标函数值 %.4f\n, fval); fprintf(退出标志 %d\n, exitflag);该问题的理论最优解是 x [4, 2]目标函数值 -8。如果输出的 exitflag 为 1说明求解成功。linprog 返回结果后要主动检查约束是否满足% 验证不等式约束 disp(A * x); disp(b);如果 A×x 明显超出 b说明问题建模有误或者求解器返回了不可行解。如果决策变量需要取整数比如调度问题里的台数、人数、设备数量用 intlinprog 替代 linprog。调用时需要多传一个整数变量索引向量。% 假设 x1 为整数x2 为连续变量 intcon 1; [x_int, fval_int] intlinprog(f, intcon, A, b, Aeq, beq, lb, ub);5.3 带约束的非线性优化fminconfmincon 是优化工具箱里使用频率最高的求解器。它可以同时处理线性约束、边界约束和非线性约束。调用格式为[x, fval, exitflag] fmincon(fun, x0, A, b, Aeq, beq, lb, ub, nonlcon, options);其中 nonlcon 是返回非线性约束的函数句柄。非线性约束函数必须返回两个输出c 是不等式约束需要满足 c(x) ≤ 0ceq 是等式约束需要满足 ceq(x) 0。这里用一个带非线性约束的测试问题最小化 (x1-3)² (x22)²要求 x1²x2² ≤ 4且 x1x2 1。先把非线性约束写成单独函数文件或者直接用局部函数。为了简单在脚本中使用函数文件方式% 主脚本 fun (x) (x(1)-3)^2 (x(2)2)^2; x0 [0; 0]; A []; b []; Aeq []; beq []; lb [-5; -5]; ub [5; 5]; nonlcon circle_constraint; [x, fval, exitflag] fmincon(fun, x0, A, b, Aeq, beq, lb, ub, nonlcon); fprintf(x [%.4f, %.4f]\n, x(1), x(2)); fprintf(目标函数值 %.4f\n, fval); fprintf(退出标志 %d\n, exitflag);function [c, ceq] circle_constraint(x) % 不等式约束x1^2 x2^2 4 c x(1)^2 x(2)^2 - 4; % 等式约束x1 x2 1 ceq x(1) x(2) - 1; end判断成功的标准是exitflag 为 1且验证约束满足程度。[c, ceq] circle_constraint(x); fprintf(不等式约束残差 %.6f\n, max(c)); fprintf(等式约束残差 %.6f\n, norm(ceq));如果 max(c) 大于 1e-6 或 norm(ceq) 过大说明约束没有被满足。常见原因包括收敛容忍度设置过松、初始点离可行域太远、或者问题本身无解。5.4 非线性方程组求解fsolvefsolve 用于求解形式为 F(x) 0 的非线性方程组。典型场景是化学平衡计算、机构运动学分析、隐式方程求反函数。这里求解一个简单方程组x1² x2² 1x1 - x2 0也就是 F(x) [x1² x2² - 1; x1 - x2]。% 定义方程组输入 x 为列向量 fun (x) [x(1)^2 x(2)^2 - 1; x(1) - x(2)]; % 初始猜测 x0 [0.5; 2]; % 求解 [x, fval, exitflag, output] fsolve(fun, x0); fprintf(x [%.6f, %.6f]\n, x(1), x(2)); fprintf(函数残差范数 %.6e\n, norm(fval)); fprintf(退出标志 %d\n, exitflag);该方程组的真实解是 x [√2/2, √2/2] 和 x [-√2/2, -√2/2]。从不同初始点出发可能收敛到不同解。fsolve 强调初始点的重要性因为非线性方程组通常有多个解。判断成功的标准是norm(fval)足够小比如小于 1e-6且 exitflag 为 1。如果 exitflag 为 0 或负数优先调整 x0然后调整 TolFun 和 MaxFunctionEvaluations。5.5 曲线拟合与参数辨识lsqcurvefitlsqcurvefit 是工程里使用频率非常高的函数。它解决的是最小二乘拟合问题给定一组自变量 t 和观测值 y找到模型参数 p使模型输出 model(p, t) 与 y 的误差平方和最小。以下代码用指数衰减模型拟合带噪声数据。% 生成模拟数据 rng(0); % 固定随机种子保证结果可复现 t linspace(0, 5, 20); y 2.5 * exp(-0.8 * t) 0.3 * randn(20, 1); % 定义模型函数y p(1) * exp(p(2) * t) model (p, t) p(1) * exp(p(2) * t); % 初始参数猜测 p0 [1; -0.1]; % 拟合 [p, resnorm, residual, exitflag] lsqcurvefit(model, p0, t, y); fprintf(p [%.4f, %.4f]\n, p(1), p(2)); fprintf(残差平方和 %.6f\n, resnorm); fprintf(退出标志 %d\n, exitflag); % 预测并对比 y_fit model(p, t); plot(t, y, o, t, y_fit, -); legend(观测数据, 拟合曲线);真实参数为 p [2.5, -0.8]。因为有随机噪声拟合结果会有微小偏差只要 p(1) 接近 2.5、p(2) 接近 -0.8说明拟合成功。resnorm 越小说明拟合效果越好。5.6 全局优化与遗传算法ga当目标函数是多峰的、非光滑的或者有大量局部极值点时局部求解器 fmincon 很容易陷入局部最优。这时可以用遗传算法 ga 做全局搜索。这个求解器在老版本中由全局优化工具箱提供新版本中已经合并进优化工具箱。仍然用 5.3 节的约束问题测试 gafun (x) (x(1)-3)^2 (x(2)2)^2; lb [-5; -5]; ub [5; 5]; A []; b []; Aeq []; beq []; % 非线性约束函数不能直接用句柄ga 需要自定义函数 nonlcon circle_constraint; % 遗传算法优化 [x_ga, fval_ga, exitflag_ga] ga(fun, 2, A, b, Aeq, beq, lb, ub, nonlcon); fprintf(ga 结果x [%.4f, %.4f]\n, x_ga(1), x_ga(2)); fprintf(目标函数值 %.4f\n, fval_ga); fprintf(退出标志 %d\n, exitflag_ga);注意ga 是随机算法每次运行结果可能不同。如果想要结果稳定可以设置随机数种子或增加种群规模和代数。rng(0); options_ga optimoptions(ga, PopulationSize, 200, MaxGenerations, 500); [x_ga, fval_ga] ga(fun, 2, A, b, Aeq, beq, lb, ub, nonlcon, options_ga);判断 ga 是否成功的标准不是单纯看 exitflag还要看目标函数值是否比 fmincon 的结果更优或 
长河君 · 2026/8/31 1:32:32 · 阅读 1.2k MATLAB优化工具箱实战指南：线性规划、非线性拟合与全局优化
想用 MATLAB 解决优化问题最快的方法是直接使用官方优化算法工具箱Optimization Toolbox。它解决的不是某一个具体问题而是一整类问题线性规划、带线性或非线性约束的优化、无约束优化、二次规划、最小二乘数据拟合、非线性方程组求解全部封装成统一函数接口。你不需要从零写迭代算法不需要手搓梯度下降甚至不需要自己判断该用哪种数值方法只需要把目标函数、约束、初始点按 MATLAB 的规则传进去求解器就会自动处理算法选择、迭代终止、数值稳定性这些细节。从工程使用角度来看这个工具箱最值得关注的几个点第一覆盖面广从一元函数最小值到多变量强非线性约束问题都有对应求解器第二接口风格统一所有求解器都遵循“目标函数 约束 选项”的调用形式学会一个就能快速迁移到其他问题第三扩展性强可以结合遗传算法、模拟退火、模式搜索等全局优化方法处理多峰、非光滑、离散变量问题也能通过 for/parfor 批量跑大量初值和多组数据。下面我会从环境检查、工具箱启用、常用求解器、批量任务、性能调优到排错给出一套可以照跑的完整流程。上一期提到的遗传算法工具箱和这里的优化工具箱是同一条技术体系里的不同层次建议结合起来看。1. 核心能力速览能力项说明项目类型MathWorks 官方 MATLAB 优化工具箱主要功能线性规划、二次规划、无约束/约束非线性优化、最小二乘拟合、非线性方程组求解、多目标优化常用求解器fminunc、fmincon、linprog、intlinprog、quadprog、fminbnd、fsolve、lsqcurvefit、lsqnonlin全局优化求解器ga、patternsearch、simulannealbnd、multistart、globalsearch需配合全局优化能力使用调用方式命令行函数调用、.m 脚本调用、实时编辑器 Optimize 任务、传统 App 界面是否支持批量支持可结合 for/parfor 批量求解多个初值、多组数据是否支持并行支持可配合 Parallel Computing Toolbox 使用硬件要求普通办公电脑即可运行大规模问题建议配置充足内存支持平台Windows、Linux、macOS 上的 MATLAB许可要求商业授权必须使用正版许可这里需要单独说明一下版本差异。在新版本 MATLAB 中全局优化工具箱和优化工具箱已经合并ga、patternsearch 这些求解器直接通过优化工具箱提供。如果你使用的是老版本可能还需要单独安装 Global Optimization Toolbox 才能用遗传算法相关函数。所以后面代码里用到 ga 时我会明确标注这一点。2. 适用场景与使用边界这个工具箱适合解决以下类型的问题参数辨识与数据拟合比如根据实验数据反推模型参数对应 lsqcurvefit、lsqnonlin。工程设计优化比如结构参数、尺寸、控制参数寻优对应 fmincon、fminunc。资源调度与路径规划中的线性规划、整数规划对应 linprog、intlinprog。求解复杂的非线性方程组比如多元隐式方程、平衡点计算对应 fsolve。需要跳出局部最优的全局优化任务配合 ga、multistart、patternsearch 使用。但也要说清楚边界。这个工具箱不是万能建模软件它解决的问题类型有明确数学形式如果问题不是凸优化fminunc、fmincon 这类局部求解器只能保证找到局部最优不保证全局最优。如果变量是离散整数普通 fmincon 不能直接处理要用 intlinprog 或 ga。如果目标函数无法写成可计算的 MATLAB 函数比如完全依赖外部软件仿真那每次求解都需要调用外部程序效率会明显下降。如果问题本身没有可行解求解器会返回不可行状态这时候需要回头检查约束是否自相矛盾。最后是合规边界。这个工具箱是商业软件需要从 MathWorks 官方渠道获得授权许可。不要使用任何非官方激活方式容易出现函数缺失、许可证报错也存在安全隐患。如果你用 MATLAB 处理真实业务数据、人脸数据、内部图纸或受版权保护的素材要确认自己有权使用这些数据。特别是把优化结果发布到论文、报告或商用产品时务必核验数据的来源和授权链条。优化代码本身是一个工具工具怎么用、用在哪里由使用者负责。3. 环境准备与前置条件开始写代码之前先确认 MATLAB 环境是否已经具备优化工具箱。打开 MATLAB 命令窗口执行下面这几条命令。% 查看优化工具箱版本装了会显示版本号和描述信息 ver(optim)% 检测许可证是否可用返回 1 表示可用返回 0 表示不可用 license(test, Optimization_Toolbox)% 检查关键求解器是否在搜索路径上 which fmincon which linprog如果ver(optim)显示出版本信息说明工具箱已安装。如果license(test, Optimization_Toolbox)返回 0说明安装过但许可证没有覆盖当前用户需要联系 MATLAB 授权管理员开通。如果which fmincon返回空说明工具箱不在 MATLAB 搜索路径中。老版本环境还需要额外检查全局优化工具箱% 老版本中遗传算法等函数由全局优化工具箱提供 ver(globaloptim) license(test, Global_Optimization_Toolbox)对于一般的学习型和科研型优化问题普通办公电脑即可满足要求。复杂工程优化问题建议内存 16GB 以上并在求解时关注迭代日志和内存占用。MATLAB 本身支持 Windows、Linux、macOS 三大平台优化工具箱在这三个平台上的行为完全一致。还要处理好工作目录。建议所有优化代码放在独立文件夹里用cd切换过去或者把这个文件夹加入搜索路径。cd(D:\matlab_opt_project) addpath(D:\matlab_opt_project)4. 安装部署与工具箱启用如果你的 MATLAB 环境里没有优化工具箱要通过官方渠道补装。常见方式有两种。第一种是学校的校园授权或企业授权通常管理员会统一开通工具箱你在 MATLAB 里执行ver(optim)能直接看到。第二种是个人授权可以在 MATLAB 主页选项卡中打开“附加功能”管理器找到“管理附加功能”从 MathWorks 官网获取新增产品。也可以重新运行 MATLAB 安装程序在产品选择界面勾选 Optimization Toolbox。安装完以后需要重启 MATLAB 才能生效。重启后不要急着写复杂代码先跑一个最小示例验证工具箱可用性。% 一元函数最小值求解验证工具箱核心功能 fun (x) (x - 3).^2 2; [x, fval] fminbnd(fun, -10, 10); fprintf(最优解 x %.6f目标函数值 %.6f\n, x, fval);如果命令窗口输出类似“最优解 x 3.000000目标函数值 2.000000”的内容说明工具箱已经正常工作。这里必须强调整个过程只使用正版许可。有些破解版 MATLAB 表面能运行基础功能但到调用 fmincon、linprog 这类工具箱函数时会出现“Undefined function”或许可证弹窗很难排查。从官方渠道获取授权是最稳妥的路线。5. 功能测试与效果验证下面按照不同问题类型逐一跑通常用求解器。所有代码都可以直接复制到 MATLAB 中运行每个示例都会给出判断成功与否的标准。5.1 无约束优化fminuncfminunc 用于求解无约束多变量函数的最小值典型场景是参数估计、模型校准。这里用经典的 Rosenbrock 函数测试函数表达式为 f(x) 100×(x2 - x1²)² (1 - x1)²全局最小值在 (1, 1) 处目标函数值为 0。% 目标函数Rosenbrock 函数 fun (x) 100 * (x(2) - x(1)^2)^2 (1 - x(1))^2; % 初始点 x0 [-1.2; 1]; % 调用求解器 [x, fval, exitflag, output] fminunc(fun, x0); fprintf(x [%.4f, %.4f]\n, x(1), x(2)); fprintf(目标函数值 %.6f\n, fval); fprintf(退出标志 %d\n, exitflag); fprintf(迭代次数 %d\n, output.iterations);运行成功后x 应该接近 [1, 1]目标函数值接近 0退出标志为 1。如果 exitflag 为 0说明达到最大迭代次数或函数评估次数上限需要增大 MaxFunctionEvaluations 或 MaxIterations。fminunc 适用的前提是目标函数连续可导。如果目标函数不可导、不连续或者存在大量局部极值点fminunc 不一定能稳定求解这时要改用 ga 或 patternsearch。5.2 线性规划与整数规划linprog、intlinprog线性规划是实际工程中最常见的规划类型。linprog 求解形式为 min f×x满足线性不等式 A×x ≤ b、线性等式 Aeq×x beq、边界 lb ≤ x ≤ ub。下面是一个两变量线性规划示例。% 目标函数系数即最小化 -x1 - 2x2 f [-1; -2]; % 线性不等式约束 A*x b A [1, 1; -1, 1; 0, 1]; b [6; -2; 3]; % 无等式约束 Aeq []; beq []; % 变量下界和上界 lb [0; 0]; ub []; % 求解 [x, fval, exitflag] linprog(f, A, b, Aeq, beq, lb, ub); fprintf(x [%.4f, %.4f]\n, x(1), x(2)); fprintf(目标函数值 %.4f\n, fval); fprintf(退出标志 %d\n, exitflag);该问题的理论最优解是 x [4, 2]目标函数值 -8。如果输出的 exitflag 为 1说明求解成功。linprog 返回结果后要主动检查约束是否满足% 验证不等式约束 disp(A * x); disp(b);如果 A×x 明显超出 b说明问题建模有误或者求解器返回了不可行解。如果决策变量需要取整数比如调度问题里的台数、人数、设备数量用 intlinprog 替代 linprog。调用时需要多传一个整数变量索引向量。% 假设 x1 为整数x2 为连续变量 intcon 1; [x_int, fval_int] intlinprog(f, intcon, A, b, Aeq, beq, lb, ub);5.3 带约束的非线性优化fminconfmincon 是优化工具箱里使用频率最高的求解器。它可以同时处理线性约束、边界约束和非线性约束。调用格式为[x, fval, exitflag] fmincon(fun, x0, A, b, Aeq, beq, lb, ub, nonlcon, options);其中 nonlcon 是返回非线性约束的函数句柄。非线性约束函数必须返回两个输出c 是不等式约束需要满足 c(x) ≤ 0ceq 是等式约束需要满足 ceq(x) 0。这里用一个带非线性约束的测试问题最小化 (x1-3)² (x22)²要求 x1²x2² ≤ 4且 x1x2 1。先把非线性约束写成单独函数文件或者直接用局部函数。为了简单在脚本中使用函数文件方式% 主脚本 fun (x) (x(1)-3)^2 (x(2)2)^2; x0 [0; 0]; A []; b []; Aeq []; beq []; lb [-5; -5]; ub [5; 5]; nonlcon circle_constraint; [x, fval, exitflag] fmincon(fun, x0, A, b, Aeq, beq, lb, ub, nonlcon); fprintf(x [%.4f, %.4f]\n, x(1), x(2)); fprintf(目标函数值 %.4f\n, fval); fprintf(退出标志 %d\n, exitflag);function [c, ceq] circle_constraint(x) % 不等式约束x1^2 x2^2 4 c x(1)^2 x(2)^2 - 4; % 等式约束x1 x2 1 ceq x(1) x(2) - 1; end判断成功的标准是exitflag 为 1且验证约束满足程度。[c, ceq] circle_constraint(x); fprintf(不等式约束残差 %.6f\n, max(c)); fprintf(等式约束残差 %.6f\n, norm(ceq));如果 max(c) 大于 1e-6 或 norm(ceq) 过大说明约束没有被满足。常见原因包括收敛容忍度设置过松、初始点离可行域太远、或者问题本身无解。5.4 非线性方程组求解fsolvefsolve 用于求解形式为 F(x) 0 的非线性方程组。典型场景是化学平衡计算、机构运动学分析、隐式方程求反函数。这里求解一个简单方程组x1² x2² 1x1 - x2 0也就是 F(x) [x1² x2² - 1; x1 - x2]。% 定义方程组输入 x 为列向量 fun (x) [x(1)^2 x(2)^2 - 1; x(1) - x(2)]; % 初始猜测 x0 [0.5; 2]; % 求解 [x, fval, exitflag, output] fsolve(fun, x0); fprintf(x [%.6f, %.6f]\n, x(1), x(2)); fprintf(函数残差范数 %.6e\n, norm(fval)); fprintf(退出标志 %d\n, exitflag);该方程组的真实解是 x [√2/2, √2/2] 和 x [-√2/2, -√2/2]。从不同初始点出发可能收敛到不同解。fsolve 强调初始点的重要性因为非线性方程组通常有多个解。判断成功的标准是norm(fval)足够小比如小于 1e-6且 exitflag 为 1。如果 exitflag 为 0 或负数优先调整 x0然后调整 TolFun 和 MaxFunctionEvaluations。5.5 曲线拟合与参数辨识lsqcurvefitlsqcurvefit 是工程里使用频率非常高的函数。它解决的是最小二乘拟合问题给定一组自变量 t 和观测值 y找到模型参数 p使模型输出 model(p, t) 与 y 的误差平方和最小。以下代码用指数衰减模型拟合带噪声数据。% 生成模拟数据 rng(0); % 固定随机种子保证结果可复现 t linspace(0, 5, 20); y 2.5 * exp(-0.8 * t) 0.3 * randn(20, 1); % 定义模型函数y p(1) * exp(p(2) * t) model (p, t) p(1) * exp(p(2) * t); % 初始参数猜测 p0 [1; -0.1]; % 拟合 [p, resnorm, residual, exitflag] lsqcurvefit(model, p0, t, y); fprintf(p [%.4f, %.4f]\n, p(1), p(2)); fprintf(残差平方和 %.6f\n, resnorm); fprintf(退出标志 %d\n, exitflag); % 预测并对比 y_fit model(p, t); plot(t, y, o, t, y_fit, -); legend(观测数据, 拟合曲线);真实参数为 p [2.5, -0.8]。因为有随机噪声拟合结果会有微小偏差只要 p(1) 接近 2.5、p(2) 接近 -0.8说明拟合成功。resnorm 越小说明拟合效果越好。5.6 全局优化与遗传算法ga当目标函数是多峰的、非光滑的或者有大量局部极值点时局部求解器 fmincon 很容易陷入局部最优。这时可以用遗传算法 ga 做全局搜索。这个求解器在老版本中由全局优化工具箱提供新版本中已经合并进优化工具箱。仍然用 5.3 节的约束问题测试 gafun (x) (x(1)-3)^2 (x(2)2)^2; lb [-5; -5]; ub [5; 5]; A []; b []; Aeq []; beq []; % 非线性约束函数不能直接用句柄ga 需要自定义函数 nonlcon circle_constraint; % 遗传算法优化 [x_ga, fval_ga, exitflag_ga] ga(fun, 2, A, b, Aeq, beq, lb, ub, nonlcon); fprintf(ga 结果x [%.4f, %.4f]\n, x_ga(1), x_ga(2)); fprintf(目标函数值 %.4f\n, fval_ga); fprintf(退出标志 %d\n, exitflag_ga);注意ga 是随机算法每次运行结果可能不同。如果想要结果稳定可以设置随机数种子或增加种群规模和代数。rng(0); options_ga optimoptions(ga, PopulationSize, 200, MaxGenerations, 500); [x_ga, fval_ga] ga(fun, 2, A, b, Aeq, beq, lb, ub, nonlcon, options_ga);判断 ga 是否成功的标准不是单纯看 exitflag还要看目标函数值是否比 fmincon 的结果更优或