OpenGL画蜗型线、心形线、三叶曲线、四叶曲线、螺旋线

这里讲的是一些有趣而美丽的曲线，蜗型线、心形线、三叶曲线、四叶曲线、螺旋线等，都是由圆公式的一些特殊的变化带来的，我们使用了参数极坐标方程来计算曲线路径的点，这些点用做显示弧的逼近折现中直线段的端点。这些曲线主要都是通过圆半径r的变化来生成。通过元半径的不同，可生成这些图形以及其他各种各样的图形。
下面介绍这些曲线的的原始公式：
首先是圆，众所周知，圆是具有固定圆心（x0，y0）和固定半径r的图形。
x=x0+r*cos(theta)
y=y0+r*sin(theta）(theta是圆心圆上的点与x轴所成的角度)
我们要看的图形都来自圆半径的变化来实现的，x和y的坐标都遵循以上公式。一下公式都使用了参数极坐标方程来计算曲线路径。附有图。
1、蜗形线：
r=a*cos（theta）+b（a、b是常数，theta是角度）

2、心形线：
r=a*（1+cos（theta））（a是常数，theta是角度）

3、三叶曲线：
r=a*cos（3*theta）（a是常数，theta是角度）

4、四叶曲线：
r=a*cos(2*theta)（a是常数，theta是角度）

5、螺旋线：
r=a*theta（a是常数，theta是角度）

附上代码：

#include "stdafx.h"#include<GL/glut.h>#include<stdlib.h>#include<math.h>#include<iostream>using namespace std;//screenPt:定义 屏幕上的一个点。struct screenPt{    GLint x;    GLint y;};//curveName:要画的曲线的名称，1-蜗型线，2-心形线，3-三叶曲线，4-四叶曲线，5-螺旋线。typedef enum{ limacon = 1, cardioid, threeLeaf, fourLeaf, spiral } curveName;GLsizei winWidth = 600, winHeight = 500;//init//init:初始化屏幕。void init(void){    glClearColor(1.0, 1.0, 1.0, 1.0);//设置显示窗口的背景维白色。    glMatrixMode(GL_PROJECTION);    gluOrtho2D(0.0, 200.0, 0.0, 150.0);}//lineSegment:画出pt1和pt2之间的直线。void lineSegment(screenPt pt1, screenPt pt2){    glBegin(GL_LINES);        glVertex2i(pt1.x, pt1.y);        glVertex2i(pt2.x, pt2.y);    glEnd();}//drawCurve:根据curveNum的值画出对应曲线。void drawCurve(GLint curveNum){    const GLdouble twoPi = 6.283185;    const GLint a = 175, b = 60;    GLfloat r, theta, dtheta = 1.0 / float(a);    GLint x0 = 200, y0 = 250;    screenPt curvePt[2];    glColor3f(0.0,0.0,0.0);    //curvePt[0]是逼近折线的起始位置。    //默认值是：（x0，y0）。    curvePt[0].x = x0;    curvePt[0].y = y0;    switch(curveNum)    {//计算开始逼近折线的起始位置，即theta为0时的曲线上点的x值。        case limacon: curvePt[0].x += a + b; break;        case cardioid: curvePt[0].x += a + a; break;        case threeLeaf: curvePt[0].x += a; break;        case fourLeaf:curvePt[0].x += a; break;        case spiral: break;        default: break;    }    theta = dtheta;    while(theta < twoPi)    {        switch(curveNum)        {            case limacon: r = a* cos(theta) + b; break;            case cardioid: r = a * (1 + cos(theta)); break;            case threeLeaf: r = a * cos(3 * theta) ; break;            case fourLeaf:r = a * cos(2 * theta); break;            case spiral: r = (a / 20.0) * theta; break;            default: break;        }        curvePt[1].x = x0 + r * cos(theta);        curvePt[1].y = y0 + r * sin(theta);        lineSegment(curvePt[0], curvePt[1]);        curvePt[0].x = curvePt[1].x;        curvePt[0].y = curvePt[1].y;        theta += dtheta;    }}//displayFcn:显示选项,根据提示选择要显示的曲线。void displayFcn(void){    GLint curveNum;    glClear(GL_COLOR_BUFFER_BIT);    cout << "\nEnter the integer value corresponding to\n";    cout << "one of the following curve names.\n";    cout << "Press any other key to exit.\n";    cout << "\n1-limacon, 2-cardioid, 3-threeLeaf, 4-fourLeaf, 5-spiral:";    cin >> curveNum;    if(curveNum == 1 || curveNum == 2 || curveNum == 3 || curveNum == 4 || curveNum == 5)    {        drawCurve(curveNum);    }    else         exit(0);    glFlush();}//重定义窗口参数void winReshapeFcn(GLint newWidth, GLint newHeight){    glMatrixMode(GL_PROJECTION);    glLoadIdentity();    gluOrtho2D(0.0, (GLdouble) newWidth, 0.0, (GLdouble) newHeight);    glClear(GL_COLOR_BUFFER_BIT);}void main(int argc, char** argv){    glutInit(&argc, argv);//初始化GLUT。    glutInitDisplayMode(GLUT_SINGLE | GLUT_RGB);//设置显示模式，使用单缓存和RGB颜色模型。    glutInitWindowPosition(100, 100);//设置窗口位置。    glutInitWindowSize(winWidth, winHeight);//设置窗口大小。    glutCreateWindow("Draw Curves");//创造窗口。    init();//执行初始化程序。    glutDisplayFunc(displayFcn);//将图形信息送往窗口显示。    glutReshapeFunc(winReshapeFcn);//重定义窗口参数。    glutMainLoop();//循环执行。}