About implementation of lsqnonlin

3 次查看（过去 30 天）
SreeHarish Muppirisetty
评论： Matt J 2020-5-12
Hi
I have time-course data variable 'prodconcnumbers', size 18000X22 meaning 18000 data points for each of the 22 molecules in my model.
I have 30 parameters to be fit.
I invoke fitting using lsqnonlin as follows:
karbitraryvalues (My initial value parameters)=row vector of some 30 starting values.
lb and ub were again assigned bounds.
lsqnonlinoptions=optimset('Algorithm','Levenberg-Marquardt');
[kestimates,errorperspecies]= lsqnonlin(@lsqnonlin_fun, karbitraryvalues,lb,ub,lsqnonlinoptions);
So, first I built my objfunction as follows:
for eachspecies from 1 to 22
errorpermolecule(eachspecies)=prodconcnumbers(:,eachspecies)-estimatedconcnumbers(:,eachspecies);
end
So, in conclusion, my obj.function I want to minimize errorperspecies is of length 22 where as I have 30 parameters to b estimated using this lsqnonlin.
I am getting the following error:
??? Error using ==> lsqncommon at 102
Levenberg-Marquardt and Gauss-Newton algorithms do not handle bound constraints and trust-region-reflective algorithm requires at least as
many equations as variables; aborting.
What mistake am I doing? Appreciate your time and help.
  2 个评论
Matt J
Matt J 2014-5-11
lsqcurvefit might be simpler here than lsqnonlin. Your data comparison is simple subtraction.
Matt J
Matt J 2020-5-12
Abdelwahab Afifi's comment moved here:
How you define your objective function interms of the variables' vector of length 30 ?

请先登录，再进行评论。

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

采纳的回答

Matt J
Matt J 2014-5-11
Since you have bound constraints, there was no point choosing the Levenberg-Marquardt algorithm, since it will not handle those. When it saw that you were applying bounds, lsqnonlin defaulted back to the trust-region-reflective algorithm.
The other problem was that you are trying to fit 30 parameters with only 22 data points. The trust-region-reflective algorithm cannot process under-determined problems, but even if it could, you would get poor results because your problem is undetermined.
Find more data points with which to do the fit.
  4 个评论
显示 2更早的评论
Matt J
Matt J 2014-5-12
The code you've shown suggests that your objective function is only generating one error term per molecule for a total of 22 error measurements. If you really have 18000*22 measurements, couldn't you be returning the error for each of those?
SreeHarish Muppirisetty
SreeHarish Muppirisetty 2014-5-12
Thank you, that is what I already started doing with implementation of 'lsqcurve' which you suggested, i.e., calculating error at each of 18000 time points for all the 22 species. So, the order of my error (obj.function value) will be also 18000X22 like that of data. Will see how it works. Will close this question later.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Nonlinear ARX Models 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by