matlab实现最小二乘法拟合椭圆方程

用matlab实现的最小二乘法拟合椭圆方程，实验结果已验证，具体基础知识和推导过程可见http://blog.sina.com.cn/s/blog_471e6c930102x96q.html

function ellipse_paras = ellipsefit（XYn）
% 椭圆方程：x^2+Axy+By^2+Cx+Dy+E=0
% 采用最小二乘法进行拟合椭圆
% Input: X --- a vector of x measurements
%        Y --- a vector of y measurements
%        n --- the number of measurements

% M*[A B C D E]‘ = N

%初始化椭圆方程结果
ellipse_paras.A=0;
ellipse_paras.B=0;
ellipse_paras.C=0;
ellipse_paras.D=0;
ellipse_paras.E=0;


x=X;
y=Y;
xy=x.*y;
x2=x.*x;
y2=y.*y;
x3=x2.*x;
y3=y2.*y;
x2y=x2.*y;
xy2=x.*y2;
x4=x2.*x2;
y4=y2.*y2;
x3y=x3.*y;