function?[varargout]=ellipsefit（xy）
%ELLIPSEFIT?Stable?Direct?Least?Squares?Ellipse?Fit?to?Data.
%?[XcYcABPhiP]=ELLIPSEFIT（XY）?finds?the?least?squares?ellipse?that
%?best?fits?the?data?in?X?and?Y.?X?and?Y?must?have?at?least?5?data?points.
%?Xc?and?Yc?are?the?x-?and?y-axis?center?of?the?ellipse?respectively.
%?A?and?B?are?the?major?and?minor?axis?of?the?ellipse?respectively.
%?Phi?is?the?radian?angle?of?the?major?axis?with?respect?to?the?x-axis.
%?P?is?a?vector?containing?the?general?conic?parameters?of?the?ellipse.
%?The?conic?representation?of?the?ellipse?is?given?by:
%
%?P（1）*x^2?+?P（2）*x*y?+?P（3）*y^2?+?P（4）*x?+?P（5）*y?+?P（6）?=?0
%
%?S=ELLIPSEFIT（XY）?returns?the?output?data?in?a?structure?with?field?names
%?equal?to?the?variable?names?given?above?e.g.?S.Xc?S.Yc?S.A?S.B
%?S.Phi?and?S.P
%
%?Reference:?R.?Halif?and?J.?Flusser?“Numerically?Stable?Direct?Least
%?Squares?FItting?of?Ellipses“?Department?of?Software?Engineering?Charles
%?University?Czech?Republic?2000.

%?Conversion?from?conic?to?conventional?ellipse?equation?inspired?by
%?fit_ellipse.m?on?MATLAB?Central

%?D.C.?Hanselman?University?of?Maine?Orono?ME?04469
%?Mastering?MATLAB?7
%?2005-02-28
%?Rotation?angle?fixed?2005-08-09

%--------------------------------------------------------------------------
x=x（:）;?%?convert?data?to?column?vectors
y=y（:）;
if?numel（x）~=numel（y）?||?numel（x）<5
???error（‘X?and?Y?Must?be?the?Same?Length?and?Contain?at?Least?5?Values.‘）