星系模拟器开发日志(一) 如何科学地用C++画图

代码下载地址：

http://pan.baidu.com/s/1eQjiETc

2015.8.11更新：

最近突然有一个想法：写一个程序，用来模拟太阳系的行星运动，甚至是任意星球的运动。感觉这个想法非常excited，所以就准备开始写。程序的名字就叫“星系模拟器”吧，或者也可以称作“拉普拉斯的长者？”，英文名Solar Simulator

为了避免写完后过一个月看不懂代码的悲剧重演，我准备把整个开发过程都记在这里。

工具：

①一只NOIP选手

②一台电脑

③稍有的物理学常识

2015.8.13更新：

首先说一下思路。

①每个物体都是一个质点，只考虑万有引力

②先假设所有物体都在二维平面上运动，三维的以后再说

③每个物体都用一个小点表示

④怎么模拟运动呢？最早想的是能不能用方程，然后我突然想到“三体问题没有解析解”这回事……果断大模拟

先解决“在屏幕上显示图形”这回事。毕竟作为一只NOIP选手，我以前写过的都是“黑白屏幕的傻X程序”。

先找了半天，找了两个备选方案：一是按着《C++图形与游戏编程基础》用Dark GDK，但它只兼容vs2008，而且功能非常受限。二是用DirectX，但找了找教程发现相当于学习windows编程，上来就给一堆乱码般的函数和类名跪了，写个显示亮点的小程序我至于吗我……

后来经光神指点，找到一个好东西：Easy Graphics Engine（EGE），主页http://tcgraphics.sourceforge.net/，在VS2015下的安装教程见http://www.jianshu.com/p/b12163e5a0b7

环境已备好（VS2015 Community+EGE），可以开始了。

EGE的函数名称啥的都很直白，我也不介绍了，主页里写的非常清楚。先写个“画东西”的小程序压压惊吧：

#include "stdafx.h"
#include <windows.h>
#include <iostream>
#include <graphics.h>
#include <time.h>
#include <cstdio>
#include <cmath>
using namespace std;int main() {//初始化一个200*200的图像区域initgraph(200, 200);//设置绘图颜色setcolor(WHITE);//画空心圆circle(100, 100, 3);//设置填充颜色setfillcolor(WHITE);//填充圆floodfill(100, 100,WHITE);//等待按任意键getch();//关闭图像closegraph();return 0;
}

效果图：

我们将用这样的“亮点”表示物体，当然屏幕肯定比200*200大。

EGE中，坐标(x,y)代表距离左边x，距离上边y，不要搞反了

EGE还可以画像素。画一个抛物线作为练习：

#include "stdafx.h"
#include <windows.h>
#include <iostream>
#include <graphics.h>
#include <time.h>
#include <cstdio>
#include <cmath>
using namespace std;#pragma warning(disable:4996)//无视掉freopen的警告const int gsize = 600;int main() {//打开图像initgraph(gsize, gsize);//绘图for (float x = 0, y = 0 ; y < gsize ; x+=0.001, y = x*x/500) {putpixel(int(x+0.5), int(y+0.5), WHITE);//画像素}//关闭图像getch();closegraph();return 0;
}

方程y=x^2/500

效果图：

接下来该画动图了。按照EGE官网上的教程，写了一个画动图的程序：

#include "stdafx.h"
#include <windows.h>
#include <iostream>
#include <graphics.h>
#include <time.h>
#include <cstdio>
#include <cmath>
using namespace std;#pragma warning(disable:4996)//无视掉freopen的警告void mainloop()//主循环
{//将线条和填充颜色设为白色setcolor(WHITE);setfillcolor(WHITE);float x = 0, y = 0;for (; is_run(); delay_fps(60))//每秒60帧{//cleardevice();//清屏y = x*x / 500.0;fillellipsef(x, y, 3, 3);//画以(x,y)为中心,长短轴均为3的椭圆x += 1.0;//更新x}
}int main(void)
{//设置初始化图形,差不多就是默认setinitmode(INIT_DEFAULT | INIT_NOFORCEEXIT);//初始化窗口initgraph(600, 600);//初始化随机种子randomize();//设置更新窗口模式,为手动模式setrendermode(RENDER_MANUAL);//主循环mainloop();//关闭窗口closegraph();return 0;
}

这是动态绘制上述抛物线的程序。

注意，每一帧绘制完成后我没有清屏，这意味着它将绘制出整条轨迹。

画到一半的效果图：