点在椭圆内matlab代码,matlab将椭圆拟合到给定的点集

Fitting an ellipse to a given set of points

The most accurate and robust fit minimizes geometric (orthogonal) distances from the observed points to the fitting curve. Our goal is to minimize the sum of squares of orthogonal distances. The Levenberg-Marquardt algorithm requires the computation of the distances and their derivatives with respect to the ellipse parameters. So this method is generated by using implicit differentiation for computing

Jacobian matrix.

Usage: [ParG,RSS,iters] = fit_ellipseLMG(XY,ParGini,LambdaIni)

Child functions:

Residuals_ellipse(from previous submission) , JmatrixLMG (included in the main function)

Input:

XY:given pointsi=1 to n

ParGini = [Center(1:2), Axes(1:2),Angle]'

LambdaIni: the initial value of the control parameter Lambda

Output:

ParG: parameter vector of the ellipse found

RSS: the Residual Sum of Squares (the sum of squares of the distances)

iters:# of iterations