这是一个快速和非迭代椭圆拟合算法 . 用法： A = EllipseDirectFit(XY) 输入: XY(n,2) 数组是n个点的坐标 x(i)=XY(i,1), y(i)=XY(i,2) 输出: A = [a b c d e f]' 时椭圆拟合的系数向量其方程方程为：: ax^2 + bxy + cy^2 + dx + ey + f = 0,其中A是被归一化的 ||A||=1 可以转换输出的几何参数，比如（半轴，中心等），具体的理论公式在http://www.mathworks.com/
[2010 PAMI] Object Detection with Discriminatively Trained Part Based Models Pedro F. Felzenszwalb, Ross B.Girshick, David McAllester and Deva Ramanan 翻译见：http://blog.csdn.net/masibuaa/article/details/17924671
本代码是由Loy G, Zelinsky A. 2003年发表在pami上面的方法FRST的实现代码。 Loy G, Zelinsky A. Fast radial symmetry for detecting points of interest[J]. Pattern Analysis and Machine Intelligence, IEEE Transactions on, 2003, 25(8): 959-973.
