opencv最小二成拟合椭圆算法

//传入样本点，返回椭圆的5个参数

vector<double> getEllipsepar(vector<CvPoint> vec_point)
{
    vector<double> vec_result;
    double  x3y1 = 0,x1y3= 0,x2y2= 0,yyy4= 0, xxx3= 0,xxx2= 0,x2y1= 0,yyy3= 0,x1y2= 0 ,yyy2= 0,x1y1= 0,xxx1= 0,yyy1= 0;

 

    int N = vec_point.size();
    cout << N << endl;
    for (int m_i = 0;m_i < N ;++m_i )
    {
        double xi = vec_point[m_i].x ;
        double yi = vec_point[m_i].y;
        x3y1   += xi*xi*xi*yi ;
        x1y3   += xi*yi*yi*yi;
        x2y2   += xi*xi*yi*yi; ;
        yyy4   +=yi*yi*yi*yi;
        xxx3   += xi*xi*xi ;
        xxx2   += xi*xi ;
        x2y1 += xi*xi*yi;

        x1y2 += xi*yi*yi;
        yyy2 += yi*yi;