Non-linear least squares estimation with linear constraints in Simulink?

58 次查看（过去 30 天）
Benjamin Pepper
Benjamin Pepper 2024-11-11，12:40
评论： Benjamin Pepper 2024-11-12，15:25
I'm trying to implement a non-linear least squares solver in Simulink, with linear constraints. However, it seems that only the 'levenberg-marquardt' algorithm can be used in Simulink, which does not have the facility to add linear constraints.
I get the following error message when I try to use the 'interior-point' algorithm:
OPTIMOPTIONS property 'Algorithm' must be one of the following values for lsqcurvefit:
'levenberg-marquardt'
My code is:
options = optimoptions('lsqcurvefit', 'Algorithm', 'interior-point');
x = lsqcurvefit(fun, x0, Q, P, lb, ub, options);
This problem does not occur when I use the lsqcurvefit function in a regular MATLAB script.
Do you have any advice on how I could implement a non-linear least squares solver with linear constraints in Simulink?
Many thanks.

请先登录，再回答此问题。

采纳的回答

Torsten
Torsten 2024-11-11，13:58
编辑：Torsten 2024-11-11，14:16
I vaguely remember that using "lsqcurvefit" within Simulink requires to use the Levenberg-Marquardt algorithm. I can't answer why this is the case (if you are interested, you should ask Technical Support), but you simply have to accept this. But you only seem to have lower and upper bounds on the variables
x = lsqcurvefit(fun, x0, Q, P, lb, ub, options);
- so why don't you try Levenberg-Marquard ?
If this algorithm within "lsqcurvefit" doesn't succeed, what about using "fmincon" instead ?
  3 个评论
显示 1更早的评论
Torsten
Torsten 2024-11-11，19:07
If you use
sum((yexp-ysim).^2)
as objective function for "fmincon", you should get the same solution as if you specify
yexp-ysim
in "lsqnonlin" or
ysim
in "lsqcurvefit".
I don't know if special algorithms are excluded if you use "fmincon" within Simulink. Just test it.
Benjamin Pepper
Benjamin Pepper 2024-11-12，15:25
Have had some success with 'fmincon'. Seems it can be used in Simulink with 'sqp' or 'sqp-legacy' algorithms. The estimates can be good, but not as good or robust as 'lsqcurvefit'. Might try to find some C code for a non-linear least squares estimator with linear constraints, and implement that directly in Simulink.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Solver Outputs and Iterative Display 的更多信息

标签

产品


版本

R2023b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by