What does this error mean? (lsqcurvefit)

2 次查看（过去 30 天）
Vipultomar
Vipultomar 2017-1-15
编辑： Walter Roberson 2018-1-18
I have no idea what line291 is saying here. Seems pretty much advanced stuff. It came up when I was trying running a bunch of lsqcurve fit using Trust region algorithm (due to lb and ub, otherwise would have preferred Levenberg algorithm). I am iterating my input parameters myself too via for loop, as it seemed that iniital parameters were affecting final outcome too much just to see where it goes. and give best solution. But after 5 mins this came up.
Error using trdog>quad1d (line 291)
Square root error in trdog/quad1d.
Error in trdog (line 180)
[nss,tau] = quad1d(nss,ssssave,delta);
Error in snls (line 316)
[sx,snod,qp,posdef,pcgit,Z] = trdog(x,g,A,D,delta,dv,...
Error in lsqncommon (line 150)
[xC,FVAL,LAMBDA,JACOB,EXITFLAG,OUTPUT,msgData]=...
Error in lsqcurvefit (line 253)
[xCurrent,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ...
Error in lambda7 (line 211)
[phiLSQ,resnorm,fval,output]= lsqcurvefit( @(phi,x) lambdafn7(phi,x),initial_param,x,y,lb,ub,options);
  1 个评论
Star Strider
Star Strider 2017-1-15
‘What does this error mean?’
We have no idea. We don’t have your code to refer to.

请先登录，再进行评论。

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

采纳的回答

dpb
dpb 2017-1-15
I'd venture it got into a negative parameter region and hence got an imaginary value from a square root that was trapped by the internal routine quad1d. Well, let's just look and see...
Indeed,
function[nx,tau] = quad1d(x,ss,delta)
%QUAD1D 1D quadratic zero finder for trust region step
% [nx,tau] = quad1d(x,ss,delta) tau is min(1,step-to-zero)
% of a 1-D quadratic ay^2 + b*y + c.
...
% Algorithm:
% numer = -(b + sign(b)*sqrt(b^2-4*a*c));
% root1 = numer/(2*a);
% root2 = c/(a*root1); % because root2*root1 = (c/a);
...
numer = -(b + sign(b)*sqrt(b^2-4*a*c));
r1 = numer/(2*a);
r2 = c/(a*r1);
tau = max(r1,r2);
tau = min(1,tau);
if tau <= 0,
error(message('optimlib:trdog:SqrRt'));
end
nx = tau*x;
So indeed that's the cause. As Star says, not much we can say other than that w/o actual code and data to go with it but looks like your modifications didn't help in this case, anyway. What if you remove those?
  5 个评论
显示 3更早的评论
dpb
dpb 2017-1-15
"Btw which part did you want me to remove?"
The part that is related to ..." I am iterating my input parameters myself too". Wondered if that were buggering up things.
Vipultomar
Vipultomar 2017-1-15
Yes. That was the actual problem. As I am playing around with 10 independent parameters, it seemed very tedious and unsystematic to change them manually for initial guess. Which is why I am changing those parameters in systematic way near my expected solution. As there were multiple (nested) for loops, things got bit messy in the end, resulting in unexpected values of intial parameters. Now I have sorted them out and it is running fine. Thanks for help.

请先登录，再进行评论。

更多回答（1 个）

satya kothapalli
satya kothapalli 2018-1-18
编辑：Walter Roberson 2018-1-18
HI Vipul,
I also facing the same problem with my model code, as you have described in the first comment. It is interesting that you have sorted out multiple loops and identified the unexpected initial parameters in the end. Can you describe exactly what changes you have made to run smoothly. It would help me a lot.
Thank you, Regards, Satya.

类别

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

标签

Community Treasure Hunt

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

Start Hunting!

Translated by