简单好用的最小二乘椭圆拟合算法---MATLAB实现

一般椭圆公式可以写为：

MATLAB代码为（这只是一个功能函数，需要添加输入才能运行）：

%*******************************

%    作者:James_Ray_Murpy

% 函数说明：V为输入变量，是一个Nx2的矩阵，第一列和第二列分别代表x和y的值。

%                  alpha表示一般椭圆公式的系数，b表示公式中的belta。

function [alpha,b] = LSF(V)
N = size(V,1);
b = max(abs(V(:)))/255;
x = V(:,1);
y = V(:,2);
Chi = [x.^2,2*x.*y,y.^2,2*b*x,2*b*y,b^2*ones(N,1)];

 opt.issym = true;
 X = Chi'*Chi/N;
 [alpha,~] = eigs(X,1,'sm',opt);
end