B样条曲线

typedef struct {
        int x;
        int y;
        bool w;
    }gVector;
gVector vector[100];
double x0,y0,x2,y2;
void CQuadBzrDlg::OnGene2Bzr()
{
    // TODO: Add your control notification handler code here
   
    double t=0;
    double step=0.01;
    int k=0;//标识点数
    int x1=60;
    int y1=40;

    x0=30;
    y0=200;
   
    x2=200;
    y2=100;

    CDC*pDC=GetDC();   
    CPen pen,*poldPen;   
    pen.CreatePen(PS_SOLID,1,RGB(0,0,0));
    poldPen=pDC->SelectObject(&pen);
   
    vector[k].x =((1-t)*(1-t)*x0)+(2*t*(1-t)*x1)+(t*t*x2);
    vector[k].y =((1-t)*(1-t)*y0)+(2*t*(1-t)*y1)+(t*t*y2);
    pDC->MoveTo(vector[k].x ,vector[k].y );
    do{
       
        t=t+step;
        k++;
        vector[k].x =((1-t)*(1-t)*x0)+(2*t*(1-t)*x1)+(t*t*x2);
        vector[k].y =((1-t)*(1-t)*y0)+(2*t*(1-t)*y1)+(t*t*y2);
        pDC->LineTo(vector[k].x ,vector[k].y );
        pDC->MoveTo(vector[k].x ,vector[k].y );
    }while(t<=1);
    pDC->SelectObject(poldPen);
    ReleaseDC(pDC);
}

void CQuadBzrDlg::On2fitted()
{
    // TODO: Add your control notification handler code here
    //先计算P1.x
   
    double t=0.0;
    double a=0.0;
    double b=0.0;
    double c=0.0;
    double d=0.0;

    double x1,y1;
   
    int i;
   
    for(i=0;i<=100;i++)
    {
        a+=vector[i].x *(2*t*(1-t));
        t=t+0.01;
    }
    t=0.0;
    for(i=0;i<=100;i++)
    {       
        b+=x0*((1-t)*(1-t))*2*t*(1-t);
        t=t+0.01;
    }
    t=0.0;
    for(i=0;i<=100;i++)
    {       
        c+=x2*(t*t)*(2*t*(1-t));
        t=t+0.01;
    }
    t=0.0;
    for(i=0;i<=100;i++)
    {
        d+=2*t*(1-t)*2*t*(1-t);
        t=t+0.01;
    }
    x1=(a-b-c)/d;
    if((x1-(double)((int)x1))*10.0>5.0)
        x1=x1+1;

    //计算y1
    t=0.0;
    a=b=c=0.0;
    for(i=0;i<=100;i++)
    {
        a+=vector[i].y *(2*t*(1-t));
        t=t+0.01;
    }
    t=0.0;
    for(i=0;i<=100;i++)
    {       
        b+=y0*((1-t)*(1-t))*2*t*(1-t);
        t=t+0.01;
    }
    t=0.0;
    for(i=0;i<=100;i++)
    {       
        c+=y2*(t*t)*(2*t*(1-t));
        t=t+0.01;
    }
    t=0.0;
   
    y1=(a-b-c)/d;
    if((y1-(double)((int)y1))*10.0>5.0)
        y1=y1+1;


    //画出拟合曲线
    CDC*pDC=GetDC();   
    CPen pen,*poldPen;   
    pen.CreatePen(PS_SOLID,1,RGB(255,0,0));
    poldPen=pDC->SelectObject(&pen);
    //pDC->SelectObject(poldPen);
   // vector[k].x =((1-t)*(1-t)*x0)+(2*t*(1-t)*x1)+(t*t*x2);
//    vector[k].y =((1-t)*(1-t)*y0)+(2*t*(1-t)*y1)+(t*t*y2);
    pDC->MoveTo(((1-t)*(1-t)*x0)+(2*t*(1-t)*x1)+(t*t*x2) ,((1-t)*(1-t)*y0)+(2*t*(1-t)*y1)+(t*t*y2));
    do{
       
        t=t+0.01;
    //    k++;
    //    vector[k].x =((1-t)*(1-t)*x0)+(2*t*(1-t)*x1)+(t*t*x2);
    //  vector[k].y =((1-t)*(1-t)*y0)+(2*t*(1-t)*y1)+(t*t*y2);
        pDC->LineTo(((1-t)*(1-t)*x0)+(2*t*(1-t)*x1)+(t*t*x2) ,((1-t)*(1-t)*y0)+(2*t*(1-t)*y1)+(t*t*y2) );
    }while(t<=1);
    ReleaseDC(pDC);
}