避免大M取值过大引起的数值问题

在数学建模当中,常常会见到大M法,它之所以叫大M法,是因为它涉及到一个(绝对值)较大的系数M,这个大M的值应大于约束中的连续变量或者约束表达式可能取到的任何合理值,M值取过大往往会造成优化问题求解的不稳定性。举以下例子:

x ≤ 1 0 6 y x ≥ 0 y ∈ { 0 , 1 } x\leq 10^6y\\x\geq 0\\ y\in \{0,1\} x≤106yx≥0y∈{0,1}

大M约束通常用在将二元变量的信息传播给连续变量,如上述例子,只有当 y = 1 y=1 y=1 时, x x x 才能取到非0的整数值,假如求解器的整数容忍误差为 1 0 − 5 10^{-5} 10−5(详见相关文章),则 y = 0.000009999 y=0.000009999 y=0.000009999 满足整数条件,被视为等于0,但此时代入上述约束后,得到 x ≤ 9.999 x\leq 9.999 x≤9.999, x x x 最大可以取到 9 的正整数值,并不符合一开始说的只有 y = 1 y=1 y=1 时, x x x 才能取到非 0 整数值的关系。尽管我们能够调整求解器的整数容忍误差,但当大M的取值更大时,就不能避免这种情况。

因此,另一个常用的方式是通过额外的信息来收紧大M的值(详见相关文章),对于整数容忍误差为 1 0 − 5 10^{-5} 10−5 的求解器而言,大M的值比 1 0 − 5 10^{-5} 10−5 小就能实现目标,例如上述约束,在如下形式:

x ≤ 1 0 3 y x ≥ 0 y ∈ { 0 , 1 } x\leq 10^3y\\x\geq 0\\ y\in \{0,1\} x≤103yx≥0y∈{0,1}

即使 y = 0.000009999 y=0.000009999 y=0.000009999,约束也仅允许 x ≤ 0.009999 x\leq 0.009999 x≤0.009999,满足约束的初始意图。

特别的,如果我们自己用大M法实在是难以收紧大M的取值,不得不用较大的值进行约束,则可以利用一些求解器自带的SOS约束函数(当 y = 0 ⇒ x = 0 y=0\Rightarrow x=0 y=0⇒x=0),这些SOS约束函数在求解器内部也是通过大M法进行转化,且会以一定的额外求解时间为代价。

相关推荐
大江东去浪淘尽千古风流人物17 小时前
【MASt3R-SLAM】从DUSt3R到MASt3R-SLAM:基于3D重建先验的实时稠密SLAM系统演进与深度解析
计算机视觉·数学建模·参数化模型·smpl·手部重建·mano·3d手部模型
San813_LDD18 小时前
[数据结构]图:邻接矩阵
数学建模
大江东去浪淘尽千古风流人物2 天前
【MANO】参数化三维手部模型:从1000次扫描到通用手部重建的数学原理与工程实践
计算机视觉·数学建模·参数化模型·smpl·手部重建·mano·3d手部模型
CC数学建模3 天前
2026年江西省研究生数学建模竞赛1题:空间数据分析中的过拟合识别完整思路、代码、模型、文章,全网首发高质量分享!
python·算法·数学建模
逻辑君3 天前
Foresight研究报告【20260023】
人工智能·深度学习·机器学习·数学建模
茜茜数模3 天前
2026年江西省研究生数学建模竞赛赛题1 全套获奖级解决方案 全套成品word论文+代码+思路+结果图表
数学建模
一只小小的土拨鼠3 天前
【深度解析】2026 年江西省研究生数学建模竞赛赛题 1:空间数据分析中的过拟合识别——地理探测器 q 统计量、OPGD 与过拟合检验框架
数学建模
阿里matlab建模师3 天前
【机场停机位分配】matlab实现基于遗传算法的机场停机位分配优化研究
开发语言·算法·数学建模·matlab·全国大学生数学建模竞赛
2301_764441334 天前
基于Stackelberg博弈的分散式库存模型
python·算法·数学建模
ECT-OS-JiuHuaShan5 天前
什么是对和错?——“有针对性定义域的逻辑值的真伪”:认识论终极追问的公理化裁决
数据库·人工智能·算法·机器学习·数学建模