fmincon函数的代码例子 最小二乘拟合

以下是利用fmincon函数进行最小二乘拟合的代码示例：

假设我们要拟合的函数为 y= a*sin(b*x)+c*cos(d*x)，其中a，b，c，d为拟合参数，x和y为已知数据。 首先，我们需要定义一个目标函数，即待优化的函数。对于最小二乘拟合，我们的目标函数应该是残差平方和，即：

fun = @(x) (x(1)*sin(x(2)*xdata) + x(3)*cos(x(4)*xdata) - ydata).^2;

其中，xdata和ydata为已知数据，x为待拟合的参数。

然后，我们需要指定参数的初始值和参数的上下限。这可以通过定义一个包含参数初始值和参数上下限的结构体来实现，如下所示：

x0 = [1, 1, 1, 1]; % 初始值
lb = [0, 0, 0, 0]; % 下限
ub = [Inf, Inf, Inf, Inf]; % 上限
options = optimoptions('fmincon','Algorithm','interior-point',... % 设置优化选项
'MaxFunctionEvaluations',5000,...
'MaxIterations',5000);

params = fmincon(fun,x0,[],[],[],[],lb,ub,[],options); % 进行拟合

最后，我们可以得到最优的参数估计params，从而得到拟合函数。