Calculate residuals as a single number in lsqnonlin function

3 次查看(过去 30 天)
I solve a nonlinear curve fitting problem using the lsqnonlin Matlab function.
I am interested to check the residual of my fitting each time so I am returning the residual each time I call the lsqnonlin function. I wanted to calculate a single number of the residual to see if the fits gets better after doing some changes in my data. My data are 3D images. Each time I calculated the sum of residuals using the sum() function. I ended up though with very large numbers, e.g. 7.3507e+04. Is there a better way to represent the residual in a single number than calculating the sum?
  1 个评论
Matt J
Matt J 2018-10-24
It is not clear why you think you need a "better" way. If you want the numbers to be smaller just scale down your objective function, e.g.,
objfun_new =@(x) objfun_old(x)/1e4;

请先登录,再进行评论。

回答(2 个)

Rik
Rik 2018-10-24
Yes, the sum of squares or the RMS:
SSq=sum(res.^2);
RMS=sqrt(mean(res.^2));
The downside of using the plain sum is that a large negative and a large positive will cancel out, which is clearly not what you want. Whether or not a number is large should not be a factor in your decision, just that you want to minimize it.
  1 个评论
Rik
Rik 2018-10-30
Did this suggestion solve your problem? If so, please consider marking it as accepted answer. It will make it easier for other people with the same question to find an answer. If this didn't solve your question, please comment with what problems you are still having.

请先登录,再进行评论。


Torsten
Torsten 2018-10-24
If you use the usual naming conventions of lsqnonlin, the residual is "sum(F.^2)"
Best wishes
Torsten.

类别

Help CenterFile Exchange 中查找有关 Solver Outputs and Iterative Display 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by