MATLAB中如何为lsqcurvefit设置参数上限总和≤100的约束？

核心思路：利用原生线性不等式约束 #

lsqcurvefit本身支持线性不等式约束（格式为A*x ≤ b），刚好能直接实现a+b+c ≤100的要求，配合非负下限约束，这是最简洁的方案。

具体实现步骤 #

1. 定义约束条件

   • 下限约束：lb = [0, 0, 0]，直接限制每个参数非负
   • 总和约束：构造矩阵A = [1 1 1]、向量b = 100，对应数学关系a + b + c ≤ 100
   • 上限约束：ub = [inf, inf, inf]，单个参数无需额外上限，由总和约束间接限制

2. 调用lsqcurvefit示例
   假设你的拟合函数为myfun，初始参数猜测值是x0 = [30, 30, 30]，输入数据为xdata和ydata，代码如下：

% 定义约束
lb = [0, 0, 0];
ub = [inf, inf, inf];
A = [1 1 1];
b = 100;

% 调用拟合函数，可根据需求添加优化选项
x_fit = lsqcurvefit(@myfun, x0, xdata, ydata, lb, ub, optimoptions('lsqcurvefit', 'Display', 'iter'));

MATLAB优化工具箱会自动处理线性不等式约束，确保最终拟合出的参数满足a+b+c ≤100且每个参数≥0。

关于之前尝试的问题说明 #

• 在ode45中加if语句：这是在微分方程求解阶段限制变量，而非对拟合参数施加约束，无法保证最终参数满足总和要求。
• 设单个参数上限为100：这种设置下参数总和可能远超100（比如三个参数都取100，总和300），不符合你的需求。

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