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

61 views (last 30 days)
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.

Accepted Answer

Torsten
Torsten on 11 Nov 2024 at 13:58
Edited: Torsten on 11 Nov 2024 at 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 Comments
Torsten
Torsten on 11 Nov 2024 at 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 on 12 Nov 2024 at 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.

Sign in to comment.

More Answers (0)

Products


Release

R2023b

Community Treasure Hunt

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

Start Hunting!