算法机制
如上图所示，实现的是地球（绿色）围绕太阳（红色）公转，月亮（黄色）围绕地球转，同时地球也在自转的动画。具体原理如下：
首先初始化一块画布，开启深度测试，虽然当下视角看不到遮盖效果，但是在改变观测视角的情况下，还是有必要开启深度测试，避免天体之间的前后关系混乱。
利用到了openGL里的glMatrixMode(GL_PROJECTION)函数将矩阵设置成投影模式，满足绘图所需执行的矩阵变换，glLoadIdentity()函数将当前矩阵设为单位矩阵，因为在变换前需要将矩阵设置为单位矩阵，也用到了gluPerspective(GLdouble fovy,GLdouble aspect,GLdouble zNear,GLdouble zFar)透视投影的函数来设置 眼睛上下睁开的幅度 、实际窗口的横纵比、 近截面和远截面。以及glMatrixMode(GL_MODELVIEW)将矩阵设置成视图模式，满足绘图所需执行的矩阵变换。gluLookAt()来设置观测点坐标，观测物体坐标，以及观测方向。
在绘图函数draw（）里，我定义一个全局变量day，来控制地球绕太阳旋转的角度，以及月亮绕地球旋转的角度，以及控制地球自转。太阳，地球和月亮是利用OpenGL里自带的画球函数glutWireSphere来绘制的。通过glRotatef()函数来设置天体绕中心天体旋转的角度，以及绕天体运动的轴，我选择的是都绕Z轴正方向旋转。glTranslatef()来使绘制的天体按设置的方向和距离平移。
在main函数里，首先令变量day=100，然后每次调用draw函数时，变量day都加一，day大于等于360时，重新设置为零，通过不断的绘制图形以此来实现天体动态运动的效果。然而即使开启双缓冲区，即glutSwapBuffers()函数，我发现天体运动的速度还是过快，最后只好使用time.sleep()函数，来强制限制图形绘制显示的速度，最后实验效果较为理想。
代码
from OpenGL.GL import *
from OpenGL.GLUT import *
from OpenGL.raw.GLU import *
import time
def init():
glEnable(GL_DEPTH_TEST)
glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT)
glViewport(0, 0, 200, 200) # 定义视区大小
glMatrixMode(GL_PROJECTION)
glLoadIdentity()
gluPerspective(60, 100 / 100, 1, 20)
glMatrixMode(GL_MODELVIEW)
glLoadIdentity()
gluLookAt(0.0, 0.0, 5.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0)
def draw():
global day
day = day + 1
if day >= 360:
day = 0
glEnable(GL_DEPTH_TEST)
glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT)
# 太阳
glPushMatrix()
glColor3f(1.0, 0.0, 0.0)
glutWireSphere(1.0, 20, 16)
# 地球
glRotatef(day, 0.0, 0.0, 1.0)
glTranslatef(2.0, 0.0, 0.0)
glRotatef(day / 1 * 360, 0.0, 0.0, 1.0)
glPushMatrix()
glMatrixMode(GL_MODELVIEW)
glRotatef(60.0, 0.0, 1.0, 1.0)
glColor3f(0.0, 1.0, 0.0)
glutWireSphere(0.2, 10, 8)
glPopMatrix()
# 月亮
glRotatef(day / 30 * 360, 0.0, 0.0, 1.0)
glTranslatef(0.4, 0.0, 0.0)
glColor3f(1.0, 1.0, 0.0)
glutWireSphere(0.05, 10, 8)
glPopMatrix()
glutSwapBuffers()
glFlush
time.sleep(0.05)
if __name__ == "__main__":
day = 100
glutInit() # 1. 初始化glut库
glutInitDisplayMode(GLUT_SINGLE | GLUT_RGB | GLUT_DOUBLE | GLUT_DOUBLE)
glutInitWindowPosition(100, 100)
glutInitWindowSize(400, 400)
glutCreateWindow("run")
init()
glutDisplayFunc(draw) # 3. 注册回调函数draw()
glutIdleFunc(draw)
glutMainLoop()
