数学之美--贝塞尔曲线

自从上个世纪60年代，雷诺汽车公司第一次把由手工设计车体(粘土)的任务转到由计算机来完成，二维的贝塞尔曲线就成了计算机图形学中最有用的曲线之一(继直线和椭圆之后)。在PostScript中，所有曲线都用贝塞尔曲线表示——椭圆线也用贝塞尔曲线来逼近，贝塞尔曲线也用于定义PostScript字体的字符轮廓。今天的我们要感谢Pierm Bezier，是他通过一些数学的计算和推导，最后找到了这套近乎完美的曲线公式，它的作用毫不逊于数学当中另一个奇迹——黄金分割点。不仅仅如此，其实自然界中，看似纷繁复杂的万事万物没有哪一个不蕴含着它自身的规律，或许它们就是一个一个完美的奇迹，等着用数学去发现。

     在说贝塞尔曲线之前，我首先要说一说“样条”这个概念。“样条”这个词以前指的是一片木头、橡皮或者金属，用来在纸上画曲线。比如说，如果你有一些分开的图上的点，想要在它们之间画一条曲线(内插或者外插)，首先将这些点描在绘图纸上，然后，将样条定在这些点上，并用铅笔沿着样条绕着这些点弯曲的方向画曲线。现在，“样条”指的是数学公式。用贝塞尔样条(即公式)画线，我们就能得到贝塞尔曲线。

     贝塞尔曲线之所以被广泛的应用于计算机的辅助设计，是因为贝塞尔样条(公式)具有以下几个重要的特点：

            第一，贝塞尔样条常常比较具有美感，这是一个公认的主观评价；

            第二，经过不断的调整，贝塞尔样条可以逼近任意的形状；

            第三，贝塞尔样条总是由其两个端点开始和结束的，因此它非常好控制；

            第四，贝塞尔样条没有奇点，这在计算机的设计中是非常重要的。事实上，贝塞尔曲线总是受限于由端点和控制点连接而成的四边形(称作“凸包”)。

            第五，贝塞尔曲线总是与第一个控制点到起点的直线相切，并保持同一方向；同时，也与第二个控制点到终点的直线相切，并保持同一方向。

     一条二维的贝塞尔样条是由4个点定义的，两个端点和两个控制点。假设起点是(x₀,y₀)，终点是(x₃,y₃)，两个控制点分别是(x₁,y₁)和(x₂,y₂)，下面是贝塞尔样条的参数方程：

                       x(t) = (1-t)³x₀ + 3t(1-t)²x₁+ 3t²(1-t)x₂ + t³x₃

                       y(t) = (1-t)³y₀ + 3t(1-t)²y₁+ 3t²(1-t)y₂ + t³y₃

            其中，t 的值从0到1变化，这样就可以画出曲线。

     在实际的编程中，即使不知道上面的公式也可以使用贝塞尔样条。在Window API中，有现成的函数供调用，利用该函数可以画一条或多条连接的贝塞尔样条：PolyBezier(hdc,apt,iCount)。这个函数中，hdc是设备描述表句柄；apt是POINT结构的数组，前四个点依次给出贝塞尔曲线的起点、第一个控制点、第二个控制点和终点，此后的每一条贝塞尔曲线只需给出三个点，因为后一条贝塞尔曲线的起点就是前一条贝塞尔曲线的终点，如此类推；iCount是数组中点的个数，它的值等于1加上所绘制的曲线条数的三倍。

     另外要注意，在画一系列的贝塞尔样条时，只有当第一条贝塞尔曲线的第二个控制点、第一条贝塞尔曲线的终点(也是第二条贝塞尔曲线的起点)和第二条贝塞尔曲线的第一个控制点线性相关时，也就是这三个点在同一条直线上时，曲线在连接点处才是光滑的。

     下面是一个交互式的程序，利用上面的这个函数画了一条二维的贝塞尔曲线。程序中，这条贝塞尔曲线的两个控制点，可以通过按住鼠标左键或右键拖动鼠标分别进行改动。程序的源码如下：

///  
// Bezier.c文件  
 
#include <windows.h>  
  
LRESULT CALLBACK WndProc(HWND,UINT,WPARAM,LPARAM);  
void DrawBezier(HDC hdc,POINT apt[]);  
  
int WINAPI WinMain(HINSTANCE hInstance,HINSTANCE hPrevInstance,  
       LPSTR lpCmdLine,int nShowCmd)  
{  
 static TCHAR szAppName[] = TEXT("Bezier");  
 HWND hwnd;  
 MSG  msg;  
 WNDCLASSEX wndclass;  
  
 wndclass.cbSize = sizeof(WNDCLASSEX);  
 wndclass.style = CS_HREDRAW | CS_VREDRAW;  
 wndclass.cbClsExtra = 0;  
 wndclass.cbWndExtra = 0;  
 wndclass.lpszClassName = szAppName;  
 wndclass.lpfnWndProc = WndProc;  
 wndclass.hInstance = hInstance;  
 wndclass.hIcon = LoadIcon(NULL,IDI_APPLICATION);  
 wndclass.hCursor = LoadCursor(NULL,IDC_ARROW);  
 wndclass.hbrBackground = (HBRUSH)GetStockObject(WHITE_BRUSH);  
 wndclass.lpszMenuName = NULL;  
 wndclass.hIconSm = NULL;  
  
 RegisterClassEx(&wndclass);  
  
 hwnd = CreateWindowEx(WS_EX_CLIENTEDGE,szAppName,TEXT("Bezier"),WS_OVERLAPPEDWINDOW,  
  CW_USEDEFAULT,CW_USEDEFAULT,CW_USEDEFAULT,CW_USEDEFAULT,NULL,NULL,  
  hInstance,NULL);  
  
 ShowWindow(hwnd,nShowCmd);  
 UpdateWindow(hwnd);  
  
 while(GetMessage(&msg,NULL,0,0))  
 {  
  TranslateMessage(&msg);  
  DispatchMessage(&msg);  
 }  
 return msg.wParam;  
}  
  
LRESULT CALLBACK WndProc(HWND hwnd,UINT message,WPARAM wParam,LPARAM lParam)  
{  
 static POINT apt[4];  
 HDC hdc;  
 int cxClient,cyClient;  
 PAINTSTRUCT ps;  
  
 switch(message)  
 {  
 case WM_SIZE:  
  cxClient = LOWORD(lParam); //低16位为客户区宽度  
  cyClient = HIWORD(lParam); //高16位为客户区高度  
   
  apt[0].x = cxClient / 4; //第一个端点，也称为起始点  
  apt[0].y = cyClient / 2;  
    
  apt[1].x = cxClient / 2; //第一个控制点  
  apt[1].y = cyClient / 4;  
  
  apt[2].x = cxClient / 2; //第二个控制点  
  apt[2].y = 3*cyClient / 4;  
  
  apt[3].x = 3*cxClient / 4; //第二个端点，也称为终点  
  apt[3].y = cyClient / 2;  
  
  return 0;  
  
 case WM_LBUTTONDOWN:  
 case WM_RBUTTONDOWN:  
 case WM_MOUSEMOVE:  
  if(wParam & MK_LBUTTON || wParam & MK_RBUTTON)  
  {  
   hdc = GetDC(hwnd);  
   SelectObject(hdc,GetStockObject(WHITE_PEN));  
   DrawBezier(hdc,apt);   //使用白色的笔插除以前的线  
  
   if(wParam & MK_LBUTTON)  
   {  
    apt[1].x = LOWORD(lParam); //低16位为鼠标点的x值  
    apt[1].y = HIWORD(lParam); //高16位为鼠标点的y值  
   }  
  
   if(wParam & MK_RBUTTON)  
   {  
    apt[2].x = LOWORD(lParam);   
    apt[2].y = HIWORD(lParam);   
   }  
  
   SelectObject(hdc,GetStockObject(BLACK_PEN));  
   DrawBezier(hdc,apt);   //使用黑色的笔重新画  
   ReleaseDC(hwnd,hdc);  
  }  
  return 0;  
  
 case WM_PAINT:  
  InvalidateRect(hwnd,NULL,TRUE);  
  hdc = BeginPaint(hwnd,&ps);  
  DrawBezier(hdc,apt);  
  EndPaint(hwnd,&ps);  
  return 0;  
  
 case WM_DESTROY:  
  PostQuitMessage(0);  
  return 0;  
 }  
 return DefWindowProc(hwnd,message,wParam,lParam);  
}  
  
void DrawBezier(HDC hdc,POINT apt[])  
{  
 PolyBezier(hdc,apt,4); //Windows提供的画一条贝塞尔曲线的函数  
   
 MoveToEx(hdc,apt[0].x,apt[0].y,NULL); //画起始点和第一个控制点的连线  
 LineTo(hdc,apt[1].x,apt[1].y);  
  
 MoveToEx(hdc,apt[2].x,apt[2].y,NULL); //画第二个控制点和终点的连线  
 LineTo(hdc,apt[3].x,apt[3].y);  
}