中心点画圆算法

参考网址：
http://wenku.baidu.com/link?url=xlbhkU71OkL1PTcf-q01q6FGT1S10S9S91UbkF-Ddf1XvIL4lc6gbAlMIlTvD37iaIcxBlIIwtQpsnBTCw4G0cs14PMSLFNJpSesvgGpVBW

圆心在原点（0,0），半径为R。首先画八分之一的圆，如下图所示：

只要画出右图中的八分之一的圆，利用圆的对称性，画出剩余7/8的的圆。
圆的方程：
改写成如下的形式：

任何点(x,y)的相对位置可由圆函数的符号来检测：
1、F(x,y) = 0，则表示点(x,y)在圆上；
2、F(x,y) > 0，则表示点(x,y)在圆的外侧；
3、F(x,y) < 0，则表示点(x,y)在圆的内侧；

如上图所示，两个弧A和B，假设当前的点为Pi(Xi,Yi)，那么下面确定是使用点P1还是点P4。其中M点为P1和P4的中点。
M点的坐标为（Xi+1,Yi-1/2）。
分情况讨论：
1、F(M)=0，则M点在圆上，此时可以取P1点或者P4点。
2、F(M)>0，则此时如弧线A所示，此时M点在弧线的外侧，那么比较P1和P4，当然是P4距离理想点P3更近，所以此时选择P4点。
3、F(M)<0，则此时如弧线B所示，此时M点在弧线的内侧，那么比较P1和P4，当然是
P1点距离理想点P2更近，所以此时选择P1点。

M点(P1和P4的中点)的决策函数为：

这里分类讨论di+1的决策函数。
情况1：当di小于0，此时选择P1点，那么当选择P1点作为最终绘制的点时，那么把他看做这里的Pi，那么下一个中点为M’：

为什么不往P1的右上方找中心点呢？因为遵从弧线的方向，中心点必然为M’。
此时：M’的坐标为：
代入我们的式子F（x,y），得到：

此时我们知道了di以及di+1，那么用di+1减去di，得到：

情况2：当di大于等于0时，此时选择P4点，那么当选择P4点作为最终绘制的点时，那么把他看作这里的Pi，那么下一个中点为M”：

此时：M”的坐标为：
代入我们的式子F(x,y)，得到：

此时我们知道了di以及di+1，那么用di+1减去di，得到：

我们利用递推计算这八分之一圆弧上的每个点，每次迭代需要两步代理：
（1）用前一次迭代算出的决策变量的符号来决定本次选择的点；
（2）对本次选择的点，重新递推计算出新的决策变量的值；

接下的问题，就是确定d0，如下图所示。对于初始的点（0，R），顺时针生成八分之一圆，下一个中心点M的坐标是：

画八分之一的圆：

圆的对称性，得到8个对称点。

代码如下：

#include "stdafx.h"
#include <windows.h>
#include <GL/glu.h>
#include <GL/gl.h>
#include <GL/glut.h>
#include<stdlib.h>  
#include<math.h>

// 圆的对称性，算出8个对称点
void DrawPoints(double x0, double y0, double x, double y)
{
    glVertex3f(x0 + x, y0 + y, 0);
    glVertex3f(x0 - x, y0 + y, 0);
    glVertex3f(x0 + x, y0 - y, 0);
    glVertex3f(x0 - x, y0 - y, 0);

    glVertex3f(x0 + y, y0 + x, 0);
    glVertex3f(x0 - y, y0 + x, 0);
    glVertex3f(x0 + y, y0 - x, 0);
    glVertex3f(x0 - y, y0 - x, 0);
}

//八分之一圆
void DrawCircle2()
{
    glBegin(GL_POINTS);
    double r = 400;
    double x = 0;
    double y = r;
    double d = 5 / 4 - r;
    while (x <= y)
    {
        if (d >= 0)
        {
            y = y - 1;
            d = d + 2 * (x - y) + 5;
        }
        else
        {
            d = d + 2 * x + 3;
        }
        x = x + 1;
        DrawPoints(0, 0, x, y);
    }
    glEnd();
    glFlush();
}

//主循环方法
void display(void)
{
    glClear(GL_COLOR_BUFFER_BIT);
    glColor3f(1.0f, 0.0f, 1.0f);
    gluOrtho2D(-500,500,-500,500); //屏幕的坐标左下角为(-500,500)，右上角为(500,500)。
    glPointSize(3);
    DrawCircle2();
}

int main(int argc, char* argv[])
{
    glutInit(&argc, argv);
    glutInitDisplayMode(GLUT_SINGLE | GLUT_RGB);
    glutInitWindowSize(400, 400);
    glutInitWindowPosition(200, 200);
    glutCreateWindow("Circle");
    glutDisplayFunc(display);
    glutMainLoop();
    return(0);
}

结果：