fitting a circle with fitnlm

Question

Hoschang Noori 2020-7-4

2
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/559322-fitting-a-circle-with-fitnlm

编辑： Star Strider 2020-7-4

采纳的回答： Star Strider

coordinates.txt

在 MATLAB Online 中打开

Hello, I tried to fit following data to a circle with help of fitnlm

i used following circle function to determine the center of the circle and the radius:

X = readmatrix('coordinates.txt') %you can find the coordinats.txt in the attachments
circlefun = @(X, b) (X(:,1).^2 + X(:,2).^2 + b(1)*X(:,1) + b(2)*X(:,2) + b(3));
y = zeros(length(X(:,1)),1);
beta0 = [0 0 400];
mdl = fitnlm(X,y,circlefun, beta0)

it gives me following output:

Error using nlinfit (line 219)
MODELFUN must be a function that returns a vector of fitted values the same size as Y (1000-by-1).  The model function you provided returned a result that was 1-by-1.
One common reason for a size mismatch is using matrix operators (*, /, ^) in your function instead of the corresponding elementwise operators (.*, ./, .^).
Error in NonLinearModel/fitter (line 1127)
                    nlinfit(X,y,F,b0,opts,wtargs{:},errormodelargs{:});
Error in classreg.regr.FitObject/doFit (line 94)
            model = fitter(model);
Error in NonLinearModel.fit (line 1434)
            model = doFit(model);
Error in fitnlm (line 99)
model = NonLinearModel.fit(X,varargin{:});

What I already tried: I used Newton Method for multivariables and the results were quite good in my opinion, but i would like to understand how i can use fitnlm.

Thanks in advance!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Star Strider 2020-7-4

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/559322-fitting-a-circle-with-fitnlm#answer_460883

在 MATLAB Online 中打开

The arguments to ‘circlefun’ are reversed. The parameter vector must always be the first argument, and the independent variable vector (or matrix) must always be the second.

This works:

circlefun = @(b,X) (X(:,1).^2 + X(:,2).^2 + b(1)*X(:,1) + b(2)*X(:,2) + b(3));

and produces:

mdl = 
Nonlinear regression model:
    y ~ (x1^2 + x2^2 + b1*x1 + b2*x2 + b3)
Estimated Coefficients:
           Estimate        SE        tStat       pValue  
          ___________    _______    _______    __________
    b1        -156.27    0.79825    -195.76             0
    b2         51.732     5.9203     8.7381    9.8825e-18
    b3    -3.4442e+05     3394.1    -101.48             0
Number of observations: 1000, Error degrees of freedom: 997
Root Mean Squared Error: 1.08e+03
R-Squared: -Inf,  Adjusted R-Squared -Inf
F-statistic vs. constant model: 0, p-value = 1

.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Image Analyst 2020-7-4

How do you plot the fitted circle once you have the model?

Star Strider 2020-7-4

编辑：Star Strider 2020-7-4

在 MATLAB Online 中打开

Since it is an implicit equation, the intent appears to be to estimate the parameters. Everything else would be derived from them.

EDIT — (4 Jul 2020 at 17:55)

The plot is relatively straightforward:

B = mdl.Coefficients.Estimate;
Xm = -B(1)/2;
Ym = -B(2)/2;
R = sqrt((Xm^2 + Ym^2) - B(3));
A = atan2(X(:,2)-Ym, X(:,1)-Xm);
Ycir = R*sin(A) + Ym;
Xcir = R*cos(A) + Xm;
figure
plot(X(:,1), X(:,2), 'p')
hold on
plot(Xcir, Ycir, '-r', 'LineWidth',2)
hold off
grid
axis('equal')

producing:

.

请先登录，再进行评论。

Answer 2

Image Analyst 2020-7-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/559322-fitting-a-circle-with-fitnlm#answer_460925

You might also be interested in the FAQ: https://matlab.fandom.com/wiki/FAQ#How_can_I_fit_a_circle_to_a_set_of_XY_data.3F

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Hoschang Noori 2020-7-4

编辑：Hoschang Noori 2020-7-4

Thanks, I know the FAQ, i saw some posts from you about this elsewhere. :D The solution presented there works quite well for circles with few outliers, but it seems not to be robust enough for my purposes since i have some noise + some more circles in my data.

请先登录，再进行评论。