一般椭圆公式可以写为：
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