[GL]行星运行2——功能实现

1.项目描述

　　实现太阳系的模拟，包括天体的公转和自转。

2. 实现思路

　　

3. 功能实现

　　打开VS2010，新建C++ Win32控制台空项目；新建Plant.cpp文件；添加文件夹Data，将位图文件拷贝到文件夹中并添加到Data中。

　　工程属性->配置（设为‘所有配置’）->配置属性->C/C++->预处理器定义(设为GLUT_BUILDING_LIB)

　　链接器->输入->附加依赖项中添加glu32.lib;glut32.lib;glut.lib;GLaux.lib;opengl32.lib

　　添加Plant.cpp代码：

#include <windows.h>
#include <stdio.h>
#include <math.h>
#include <GL/gl.h>
#include <GL/glu.h>
#include <GL/glut.h>
#include "gl/glaux.h"
#define PI                3.141592654
#define DEG2RAD(angle)    angle*(PI/180)

static GLdouble viewPoint[]= {1.6,8.5,-0.3};
static double Eyear = 0,Vyear=0,Meryear=0,Maryear=0, Jyear=0,Syear=0,Uyear=0,Nyear=0,
              month=0,
              Eday = 0,Merday=0,Vday=360,Marday=0,Jday=0,Sday=0,Uday=0,Nday=0,
              angle,r;
static double MerDSpeed=0.72,VDSpeed=0.15,EDSpeed=36.0,MarDSpeed=35.0,JDSpeed=90.0,SDSpeed=81.0,UDSpeed=49.6,NDSpeed=53.6,
           monSpeed=12,
           MerYSpeed=4,VYSpeed=1.6,EYSpeed=1,MarYSpeed=0.5,JYSpeed=0.09,SYSpeed=0.03,UYSpeed=0.01,NYSpeed=0.006;
GLuint texture[10];


void AutoRun(void)
{   /*自转控制*/
    Merday += MerDSpeed;
    if( Merday > 360.0 ) Merday -= 360.0;
    Vday-=VDSpeed;
    if( Vday <= 0 ) Vday += 360.0;    // 金星是逆向自转
    Eday += EDSpeed;
    if( Eday > 360.0 ) Eday -= 360.0;
    Marday+=MarDSpeed;
    if( Marday > 360.0 ) Marday -= 360.0;
    Jday+=JDSpeed;
    if( Jday > 360.0 ) Jday -= 360.0;
    Sday+=SDSpeed;
    if( Sday > 360.0 ) Sday -= 360.0;
    Uday+=UDSpeed;
    if( Uday > 360.0 ) Uday -= 360.0;
    Nday+=NDSpeed;
    if( Nday > 360.0 ) Nday -= 360.0;
    /*公转控制*/
    Meryear+=MerYSpeed;
    if( Meryear > 360.0 ) Meryear -= 360.0;
    Vyear += VYSpeed;
    if( Vyear > 360.0 ) Vyear -= 360.0;
    Eyear += EYSpeed;
    if( Eyear > 360.0 ) Eyear -= 360.0;
    Maryear += MarYSpeed;
    if( Maryear > 360.0 ) Maryear -= 360.0;
    Jyear += JYSpeed;
    if( Jyear > 360.0 ) Jyear -= 360.0;
    Syear += SYSpeed;
    if( Syear > 360.0 ) Syear -= 360.0;
    Uyear += UYSpeed;
    if( Uyear > 360.0 ) Uyear -= 360.0;
    Nyear += NYSpeed;
    if( Nyear > 360.0 ) Nyear -= 360.0;

    month += monSpeed;
    if( month > 360.0 ) month -= 360.0;
    glutPostRedisplay();
}

/*绘制球体*/
void gltDrawSphere(GLfloat fRadius, GLint iSlices, GLint iStacks)
    {
    GLfloat drho = (GLfloat)(3.141592653589) / (GLfloat) iStacks;
    GLfloat dtheta = 2.0f * (GLfloat)(3.141592653589) / (GLfloat) iSlices;
    GLfloat ds = 1.0f / (GLfloat) iSlices;
    GLfloat dt = 1.0f / (GLfloat) iStacks;
    GLfloat t = 1.0f;    
    GLfloat s = 0.0f;
    GLint i, j;     

    for (i = 0; i < iStacks; i++) 
        {
        GLfloat rho = (GLfloat)i * drho;
        GLfloat srho = (GLfloat)(sin(rho));
        GLfloat crho = (GLfloat)(cos(rho));
        GLfloat srhodrho = (GLfloat)(sin(rho + drho));
        GLfloat crhodrho = (GLfloat)(cos(rho + drho));

        glBegin(GL_TRIANGLE_STRIP);
        s = 0.0f;
        for ( j = 0; j <= iSlices; j++) 
            {
            GLfloat theta = (j == iSlices) ? 0.0f : j * dtheta;
            GLfloat stheta = (GLfloat)(-sin(theta));
            GLfloat ctheta = (GLfloat)(cos(theta));

            GLfloat x = stheta * srho;
            GLfloat y = ctheta * srho;
            GLfloat z = crho;
            
            glTexCoord2f(s, t);
            glNormal3f(x, y, z);
            glVertex3f(x * fRadius, y * fRadius, z * fRadius);

            x = stheta * srhodrho;
            y = ctheta * srhodrho;
            z = crhodrho;
            glTexCoord2f(s, t - dt);
            s += ds;
            glNormal3f(x, y, z);
            glVertex3f(x * fRadius, y * fRadius, z * fRadius);
            }
        glEnd();

        t -= dt;
        }
    }

/*从文件获取BMP图片*/
AUX_RGBImageRec *LoadBMP(char *Filename)            // 载入位图图象
{
    FILE *File=NULL;                                    // 文件句柄
    if (!Filename)                                        // 确保文件名已提供
    {
        return NULL;                                    //如果没提供，返回 NULL

    }
    File=fopen(Filename,"r");                            //尝试打开文件
    if (File)                                            // 文件存在与否
    {
        fclose(File);                                    // 关闭句柄
        return auxDIBImageLoad(Filename);                //载入位图并返回指针
    }
    return NULL;                                        // 如果载入失败，返回 NULL
}
int LoadGLTextures()                                    // 载入位图并转换成纹理
{
    int Status=FALSE;                                    // Status 状态指示器
    int i;
    AUX_RGBImageRec *TextureImage[10];                    // 创建纹理的存储空间
    char bmp[10][30]={
        "Data/Sun.bmp",
        "Data/Mercury.bmp",
        "Data/Venus.bmp",
        "Data/earth.bmp",
        "Data/moon.bmp",
        "Data/Mars.bmp",
        "Data/Jutiper.bmp",
        "Data/Saturn.bmp",
        "Data/Uranus.bmp",
        "Data/Neptune.bmp",

    };
    TextureImage[0]= LoadBMP(bmp[0]);
    TextureImage[1]=LoadBMP(bmp[1]);
    TextureImage[2]=LoadBMP(bmp[2]);
    TextureImage[3]=LoadBMP(bmp[3]);
    TextureImage[4]=LoadBMP(bmp[4]);
    TextureImage[5]=LoadBMP(bmp[5]);
    TextureImage[6]=LoadBMP(bmp[6]);
    TextureImage[7]=LoadBMP(bmp[7]);
    TextureImage[8]=LoadBMP(bmp[8]);
    TextureImage[9]=LoadBMP(bmp[9]);
    for(i=0;i<10;i++)
    {
    memset(TextureImage[i],0,sizeof(void *)*1);               // 将指针设置为Null
    /* 载入BMP图片，并查错*/
    if (TextureImage[i]=LoadBMP(bmp[i]))
    {
        Status=TRUE;                                    // 将状态设置为TRUE
        glGenTextures(1, &texture[i]);                    // 创建纹理
        glBindTexture(GL_TEXTURE_2D, texture[i]);       // 根据不同Bitmap图片的数据产生不同的纹理
        /*定义纹理*/
        glTexImage2D(GL_TEXTURE_2D, 0, 3, TextureImage[i]->sizeX, TextureImage[i]->sizeY, 0, GL_RGB, GL_UNSIGNED_BYTE, TextureImage[i]->data);
        /* 控制滤波*/
        glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_WRAP_S, GL_CLAMP);
        glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_WRAP_T, GL_CLAMP);
        glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);
        glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);
    }
    if (TextureImage[i])                                    // 判断纹理是否存在
    {
        if (TextureImage[i]->data)                            // 判断纹理图像是否存在
        {
            free(TextureImage[i]->data);                    // 释放纹理图像占用的内存
        }
        free(TextureImage[i]);                                // 释放图像结构
    }
    }
    
    return Status;                                        // 返回状态
}

void display(void)
{   /*材质参数*/
    GLfloat sun_mat_ambient[4]={1.0,1.0,1.0,0.0};
    GLfloat sun_mat_diffuse[4]={1.0,1.0,0.5,1.0};
    GLfloat sun_mat_specular[4]={1.0,1.0,1.0,1.0};
    GLfloat sun_mat_shininess[]={10.0};
    GLfloat sun_mat_emission[4]={0.1,0.1,0.1,1.0};
    
    GLfloat mat_ambient[4]={0.2,0.2,0.2,1.0};
    GLfloat mat_diffuse[4]={1.0,1.0,1.0,1.0};
    GLfloat mat_specular[4]={0.5,0.5,0.5,1.0};
    GLfloat mat_shininess[]={5.0};
    
    GLfloat light_position[] = { 0.0, 0.0, 0.0, 1.0 };          // 光源位置设置

    glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);         // 清空屏幕和深度缓存
    glLoadIdentity();                                           // 重置当前的模型观察矩阵

    gluLookAt(0,viewPoint[0],viewPoint[1],0,0,0,viewPoint[2],1,0);       // 观察点设置

     /* draw sun */
    glPushMatrix();
    glLightfv(GL_LIGHT0, GL_POSITION, light_position);            // 创建光源
    glMaterialfv(GL_FRONT_AND_BACK,GL_AMBIENT,sun_mat_ambient);   // 材质设定
    glMaterialfv(GL_FRONT_AND_BACK,GL_DIFFUSE,sun_mat_diffuse);
    glMaterialfv(GL_FRONT_AND_BACK,GL_SPECULAR,sun_mat_specular);
    glMaterialfv(GL_FRONT_AND_BACK,GL_SHININESS,sun_mat_shininess);
    glMaterialfv(GL_FRONT_AND_BACK,GL_EMISSION,sun_mat_emission);
    glBindTexture(GL_TEXTURE_2D, texture[0]);                     // 绑定纹理
        glRotatef((GLfloat) Eyear, 0.0, 1.0, 0.0);
        glRotatef(90,1,0.0,0.0);
        gltDrawSphere(1.0,200,200);             //绘制球体
    glPopMatrix();

/* draw Mercury */
    glPushMatrix();
    glMaterialfv(GL_FRONT_AND_BACK,GL_AMBIENT,mat_ambient);   // 材质设定
    glMaterialfv(GL_FRONT_AND_BACK,GL_DIFFUSE,mat_diffuse);
    glMaterialfv(GL_FRONT_AND_BACK,GL_SPECULAR,mat_specular);
    glMaterialfv(GL_FRONT_AND_BACK,GL_SHININESS,mat_shininess);
    glBindTexture(GL_TEXTURE_2D, texture[1]);                 // 绑定纹理
        glRotatef((GLfloat) Meryear, 0.0, 1.0, 0.0);   //公转
        glTranslatef (0.0, 0.0, 1.8);           //与太阳的距离
        glRotatef((GLfloat) Merday, 0.0, 1.0, 0.0);    //自转
        glRotatef(90,1.0,0.0,0.0);
        gltDrawSphere(0.2,200,200);        
   glPopMatrix();

/*画水星轨道*/
   glDisable(GL_TEXTURE_2D);
   glBegin(GL_LINE_LOOP);
   glColor3f(1.0,1.0,1.0);
    for(angle=0;angle<=360;angle++)
            glVertex3f(1.8*sin(DEG2RAD(angle)),0,1.8*cos(DEG2RAD(angle)));        
   glEnd();
   glEnable(GL_TEXTURE_2D);

/* draw Venus */
    glPushMatrix();
    glMaterialfv(GL_FRONT_AND_BACK,GL_AMBIENT,mat_ambient);   // 材质设定
    glMaterialfv(GL_FRONT_AND_BACK,GL_DIFFUSE,mat_diffuse);
    glMaterialfv(GL_FRONT_AND_BACK,GL_SPECULAR,mat_specular);
    glMaterialfv(GL_FRONT_AND_BACK,GL_SHININESS,mat_shininess);
    glBindTexture(GL_TEXTURE_2D, texture[2]);                 // 绑定纹理
        glRotatef((GLfloat) Vyear, 0.0, 1.0, 0.0);
        glTranslatef (0.0, 0.0, 2.5);           //与太阳的距离
        glRotatef((GLfloat) Vday, 0.0, 1.0, 0.0);
        glRotatef(90,1.0,0.0,0.0);
        gltDrawSphere(0.22,200,200);        
    glPopMatrix();

/*画金星轨道*/
    glDisable(GL_TEXTURE_2D);
    glBegin(GL_LINE_LOOP);
    glColor3f(1.0,1.0,1.0);
    for(angle=0;angle<=360;angle++)
            glVertex3f(2.5*sin(DEG2RAD(angle)),0,2.5*cos(DEG2RAD(angle)));        
    glEnd();
    glEnable(GL_TEXTURE_2D);

/* draw Earth */
    glPushMatrix();
    glMaterialfv(GL_FRONT_AND_BACK,GL_AMBIENT,mat_ambient);   // 材质设定
    glMaterialfv(GL_FRONT_AND_BACK,GL_DIFFUSE,mat_diffuse);
    glMaterialfv(GL_FRONT_AND_BACK,GL_SPECULAR,mat_specular);
    glMaterialfv(GL_FRONT_AND_BACK,GL_SHININESS,mat_shininess);
    glBindTexture(GL_TEXTURE_2D, texture[3]);                 // 绑定纹理
        glRotatef((GLfloat) Eyear, 0.0, 1.0, 0.0);
        glTranslatef (0.0, 0.0, 3.2);           //与太阳的距离
        glRotatef((GLfloat) Eday, 0.0, 1.0, 0.0);
        glRotatef(90,1.0,0.0,0.0);
        gltDrawSphere(0.2,200,200);        
    glPopMatrix();

/*画地球轨道*/
    glDisable(GL_TEXTURE_2D);
    glBegin(GL_LINE_LOOP);
    glColor3f(1.0,1.0,1.0);
    for(angle=0;angle<=360;angle++)
            glVertex3f(3.2*sin(DEG2RAD(angle)),0,3.2*cos(DEG2RAD(angle)));        
    glEnd();
    glEnable(GL_TEXTURE_2D);

/* draw moon */
    glPushMatrix();
    glMaterialfv(GL_FRONT_AND_BACK,GL_AMBIENT,mat_ambient);   // 材质设定
    glMaterialfv(GL_FRONT_AND_BACK,GL_DIFFUSE,mat_diffuse);
    glMaterialfv(GL_FRONT_AND_BACK,GL_SPECULAR,mat_specular);
    glMaterialfv(GL_FRONT_AND_BACK,GL_SHININESS,mat_shininess);
    glBindTexture(GL_TEXTURE_2D, texture[4]);                 // 绑定纹理
        glRotatef((GLfloat) Eyear, 0.0, 1.0, 0.0);
        glTranslatef (0.0, 0.0, 3.2);           //与太阳的距离
        glRotatef((GLfloat) month, 0.0, 1.0, 0.0);
        glTranslatef (0.0, 0.0, 0.3);          //与地球的距离
        glRotatef((GLfloat) Eday, 0.0, 1.0, 0.0);
        glRotatef(90,1.0,0.0,0.0);
        glColor3f (1.0, 1.0, 1.0);
        gltDrawSphere(0.06,100,100);           
   glPopMatrix();

/* draw Mars */
       glPushMatrix();
    glMaterialfv(GL_FRONT_AND_BACK,GL_AMBIENT,mat_ambient);   // 材质设定
    glMaterialfv(GL_FRONT_AND_BACK,GL_DIFFUSE,mat_diffuse);
    glMaterialfv(GL_FRONT_AND_BACK,GL_SPECULAR,mat_specular);
    glMaterialfv(GL_FRONT_AND_BACK,GL_SHININESS,mat_shininess);
    glBindTexture(GL_TEXTURE_2D, texture[5]);                 // 绑定纹理
        glRotatef((GLfloat)  Maryear, 0.0, 1.0, 0.0);
        glTranslatef (0.0, 0.0, 3.9);           //与太阳的距离
        glRotatef((GLfloat) Marday, 0.0, 1.0, 0.0);
        glRotatef(90,1.0,0.0,0.0);
        gltDrawSphere(0.15,200,200);        
    glPopMatrix();

/*画火星轨道*/
    glDisable(GL_TEXTURE_2D);
    glBegin(GL_LINE_LOOP);
    glColor3f(1.0,1.0,1.0);
    for(angle=0;angle<=360;angle++)
            glVertex3f(3.9*sin(DEG2RAD(angle)),0,3.9*cos(DEG2RAD(angle)));        
    glEnd();
    glEnable(GL_TEXTURE_2D);

/* draw Jupiter*/
    glPushMatrix();
    glMaterialfv(GL_FRONT_AND_BACK,GL_AMBIENT,mat_ambient);   // 材质设定
    glMaterialfv(GL_FRONT_AND_BACK,GL_DIFFUSE,mat_diffuse);
    glMaterialfv(GL_FRONT_AND_BACK,GL_SPECULAR,mat_specular);
    glMaterialfv(GL_FRONT_AND_BACK,GL_SHININESS,mat_shininess);
    glBindTexture(GL_TEXTURE_2D, texture[6]);                 // 绑定纹理
        glRotatef((GLfloat)  Jyear, 0.0, 1.0, 0.0);
        glTranslatef (0.0, 0.0, 4.5);           //与太阳的距离
        glRotatef((GLfloat) Jday, 0.0, 1.0, 0.0);
        glRotatef(90,1.0,0.0,0.0);
        gltDrawSphere(0.3,200,200);        
    glPopMatrix();

/*画木星轨道*/
    glDisable(GL_TEXTURE_2D);
    glBegin(GL_LINE_LOOP);
    glColor3f(1.0,1.0,1.0);
    for(angle=0;angle<=360;angle++)
            glVertex3f(4.5*sin(DEG2RAD(angle)),0,4.5*cos(DEG2RAD(angle)));        
    glEnd();
    glEnable(GL_TEXTURE_2D);

/* draw Saturn*/
    glPushMatrix();
    glMaterialfv(GL_FRONT_AND_BACK,GL_AMBIENT,mat_ambient);   // 材质设定
    glMaterialfv(GL_FRONT_AND_BACK,GL_DIFFUSE,mat_diffuse);
    glMaterialfv(GL_FRONT_AND_BACK,GL_SPECULAR,mat_specular);
    glMaterialfv(GL_FRONT_AND_BACK,GL_SHININESS,mat_shininess);
    glBindTexture(GL_TEXTURE_2D, texture[7]);                 // 绑定纹理
        glRotatef((GLfloat)  Syear, 0.0, 1.0, 0.0);
        glTranslatef (0.0, 0.0, 5.5);           //与太阳的距离
        glRotatef((GLfloat) Sday, 0.0, 1.0, 0.0);
        glRotatef(90,1.0,0.0,0.0);
        gltDrawSphere(0.26,200,200);        
    glPopMatrix();

/*画土星的光环*/
    glPushMatrix();
    glRotatef((GLfloat)  Syear, 0.0, 1.0, 0.0);
    glTranslatef (0.0, 0.0, 5.5);
    glBegin(GL_LINE_LOOP);
    for(r=0.32;r<=0.4;r+=0.006)
    for(angle=0;angle<=360;angle++)
            glVertex3f(r*sin(DEG2RAD(angle)),0,r*cos(DEG2RAD(angle)));        
    glEnd();
    glPopMatrix();

/*画土星轨道*/
    glDisable(GL_TEXTURE_2D);
    glBegin(GL_LINE_LOOP);
    glColor3f(1.0,1.0,1.0);
    for(angle=0;angle<=360;angle++)
            glVertex3f(5.5*sin(DEG2RAD(angle)),0,5.5*cos(DEG2RAD(angle)));        
    glEnd();
    glEnable(GL_TEXTURE_2D);

/* draw Uranus*/
    glPushMatrix();
    glMaterialfv(GL_FRONT_AND_BACK,GL_AMBIENT,mat_ambient);   // 材质设定
    glMaterialfv(GL_FRONT_AND_BACK,GL_DIFFUSE,mat_diffuse);
    glMaterialfv(GL_FRONT_AND_BACK,GL_SPECULAR,mat_specular);
    glMaterialfv(GL_FRONT_AND_BACK,GL_SHININESS,mat_shininess);
    glBindTexture(GL_TEXTURE_2D, texture[8]);                 // 绑定纹理
        glRotatef((GLfloat)  Uyear, 0.0, 1.0, 0.0);
        glTranslatef (0.0, 0.0, 6.3);             //与太阳的距离
        glRotatef((GLfloat) Uday, 0.0, 1.0, 0.0);
        glRotatef(90,1.0,0.0,0.0);
        gltDrawSphere(0.24,200,200);        
    glPopMatrix();

/*画天王星轨道*/
    glDisable(GL_TEXTURE_2D);
    glBegin(GL_LINE_LOOP);
    glColor3f(1.0,1.0,1.0);
    for(angle=0;angle<=360;angle++)
            glVertex3f(6.3*sin(DEG2RAD(angle)),0,6.3*cos(DEG2RAD(angle)));        
    glEnd();
    glEnable(GL_TEXTURE_2D);

/* draw Neptune*/
    glPushMatrix();
    glMaterialfv(GL_FRONT_AND_BACK,GL_AMBIENT,mat_ambient);   // 材质设定
    glMaterialfv(GL_FRONT_AND_BACK,GL_DIFFUSE,mat_diffuse);
    glMaterialfv(GL_FRONT_AND_BACK,GL_SPECULAR,mat_specular);
    glMaterialfv(GL_FRONT_AND_BACK,GL_SHININESS,mat_shininess);
    glBindTexture(GL_TEXTURE_2D, texture[9]);                 // 绑定纹理
        glRotatef((GLfloat)  Nyear, 0.0, 1.0, 0.0);
        glTranslatef (0.0, 0.0, 7);           //与太阳的距离
        glRotatef((GLfloat) Nday, 0.0, 1.0, 0.0);
        glRotatef(90,1.0,0.0,0.0);
        gltDrawSphere(0.23,200,200);        
    glPopMatrix();

/*画海王星轨道*/
    glDisable(GL_TEXTURE_2D);
    glBegin(GL_LINE_LOOP);
    glColor3f(1.0,1.0,1.0);
    for(angle=0;angle<=360;angle++)
            glVertex3f(7*sin(DEG2RAD(angle)),0,7*cos(DEG2RAD(angle)));        
    glEnd();
    glEnable(GL_TEXTURE_2D);

    glutSwapBuffers();    
    glFlush();
}

void myinit(void) 
{
    
    GLfloat light_ambient[]={1.0, 1.0, 1.0, 1.0};     //环境光
    GLfloat light_diffuse[]={1.0, 1.0, 1.0, 1.0};     //漫反射光
    GLfloat light_specular[]={1.0, 1.0, 1.0, 1.0};    //镜面反射光

    //glClearColor(0.5f, 0.5f, 0.5f, 1.0f);              //设置背景颜色
    glLightfv(GL_LIGHT0, GL_AMBIENT, light_ambient);
    glLightfv(GL_LIGHT0, GL_DIFFUSE, light_diffuse);
    glLightfv(GL_LIGHT0, GL_SPECULAR, light_specular);

    if (LoadGLTextures())                            // 调用纹理载入子例程 
    glEnable(GL_LIGHTING);                          // 开启光照效果
    glEnable(GL_LIGHT0);
    glEnable(GL_TEXTURE_2D);                        //启动纹理映射 
    glShadeModel(GL_SMOOTH);                        // 启用阴影平滑
    
    glTexEnvf(GL_TEXTURE_ENV, GL_TEXTURE_ENV_MODE, GL_MODULATE);   // 说明映射方式
    glDepthFunc(GL_LESS);
    glClearDepth(1.0f);                                // 深度缓存创建
    glEnable(GL_DEPTH_TEST);                        // 开启深度测试
    glDepthFunc(GL_LEQUAL);        

}


void myReshape(GLsizei w, GLsizei h)
{
    glViewport(0, 0, w, h);
    glMatrixMode(GL_PROJECTION);
    glLoadIdentity();
    gluPerspective(800.0, (GLfloat) w/(GLfloat) h, 1.0, 20.0);
    glMatrixMode(GL_MODELVIEW);
    glLoadIdentity();
    glTranslatef (0.0, 0.0, -5.0);
}

void View_SpeedControl(unsigned char key ,int x,int y)
{
    switch(key)
    {

    case'z':
        viewPoint[0]+= 0.1;
        break;

    case'x':
        viewPoint[0]-= 0.1;
        break;

    case'c':
        viewPoint[1]+=0.1;
        break;

    case'v':
        viewPoint[1]-=0.1;
        break;

    case'b':
        viewPoint[2]+=0.1;
        break;

    case'n':
        viewPoint[2]-=0.1;
        break;

    case'=':
        {  MerDSpeed*=1.5,VDSpeed*=1.5,EDSpeed*=1.5,MarDSpeed*=1.5,JDSpeed*=1.5,SDSpeed*=1.5,UDSpeed*=1.5,NDSpeed*=1.5,
           monSpeed*=2,
           MerYSpeed*=2,VYSpeed*=2,EYSpeed*=2,MarYSpeed*=2,JYSpeed*=2,SYSpeed*=2,UYSpeed*=2,NYSpeed*=2;};
        break;
    case'-':
        {  MerDSpeed/=1.5,VDSpeed/=1.5,EDSpeed/=1.5,MarDSpeed/=1.5,JDSpeed/=1.5,SDSpeed/=1.5,UDSpeed/=1.5,NDSpeed/=1.5,
           monSpeed/=1.5,
           MerYSpeed/=2,VYSpeed/=2,EYSpeed/=2,MarYSpeed/=2,JYSpeed/=2,SYSpeed/=2,UYSpeed/=2,NYSpeed/=2;};
        break;

    default:
        break;

    }
    display();
}
int main(int argc, char** argv)
{
    glutInit(&argc,argv);
    glutInitDisplayMode(GLUT_DOUBLE|GLUT_RGB | GLUT_DEPTH);
    glutInitWindowSize(800,600);
    glutCreateWindow("solar system");
    myinit();
    glutDisplayFunc(display);
    glutReshapeFunc(myReshape);
    glutKeyboardFunc(View_SpeedControl);
    glutIdleFunc(AutoRun);
    glutMainLoop();

}

运行结果：