OpenGL实现的烟花粒子系统

OpenGL实现的烟花粒子系统

最近在学习OpenGL，看到一个网友写的烟花系统非常漂亮，自己也实现一下，作为学习的练习。代码分别在VC和Android上实现。

基本原理：

物理上，粒子首先向上做匀减速直线运动，s = v * t，v = a * t；x、z方向速度是0。到最高点之后，360度均匀炸开成24束小粒子，x、z方向做匀速直线运动，vx = r * sin(radian)，vz = r * cos(radian)；y方向做匀加速直线运动。拖尾的实现原理是每个粒子后面跟一串粒子，就像一只鸡妈妈后面领了大小一队鸡娃娃，粒子娃娃的位置依次是粒子妈妈最近N次运动经过的位置（N表示后面跟的粒子娃娃个数），粒子娃娃的大小分别为：size_baby = size_mother * (1 - (float)n/N)；类似地，粒子娃娃的透明效果：alpha_baby = alpha_mother * (1 - (float)n/N)，很简单，但是实现的尾迹效果不错呢，如果追求更好的效果，可以试试一些递减的非线性函数，比如1 - sin(3.141593f/2*n/N)。

数据结构：

1. const int MAX_FIRE = 5;       // 最多5个烟花   
2. const int MAX_PARTICLES = 24; // 每个烟花升空散开的小烟花个数   
3. const int MAX_TAIL = 30;      // 烟花尾迹数   
4.   
5. typedef struct {  
6.    float r, g, b;      /* color */  
7.    float x, y, z;      /* position  */  
8.    float xs, ys, zs;   /* speed  */  
9.    float xg, yg, zg;   /* gravity  */  
10.    boolean up;         /* up or down */  
11. } Particle;  
12.   
13. typedef struct {  
14.    Particle particle[MAX_PARTICLES][MAX_TAIL]; // 烟花系统数组   
15.    float life, fade, rad; // 生命，衰减速度，x-z平面上的运动速度   
16. } Fire;  
17.   
18. Fire fire[MAX_FIRE];  

初始化：

1. for(int loop = 0; loop < MAX_FIRE; loop++) {  
2.    resetFire(loop);  
3. }//for loop end   
4.      
5. void resetFire(int loop) {  
6.    // init position   
7.    float xtemp = rand()%30 - 15.f;  
8.    float ytemp = -1*rand()%5 - 15.f;//8.f;   
9.    float ztemp = -1*rand()%5 - 15.f;//100.f;   
10.    float speed = rand()%5 + 15.f;  
11.   
12.    fire[loop].life = 1.5f;//1.0f;   
13.    fire[loop].fade = (float) ((rand()%100)/20000 + 0.002);  
14.    fire[loop].rad  = rand()%3 + 4.0f;  
15.   
16.    for (int loop1 = 0; loop1 < MAX_PARTICLES; loop1++) {  
17.       Particle* pat = &fire[loop].particle[loop1][0];  
18.       //初始颜色   
19.       pat->r = 1.0f; pat->g = 1.0f; pat->b = 1.0f;  
20.       //初始位置   
21.       pat->x = xtemp; pat->y = ytemp; pat->z = ztemp;  
22.       //初始速度   
23.       pat->xs = 0.0f; pat->ys = speed; pat->zs = 0.0f;  
24.       //初始加速度   
25.       pat->xg = 0.0f; pat->yg = -5.f; pat->zg = 0.0f;  
26.       pat->up = true;  
27.   
28.       //尾部初始化   
29.       for(int loop2 = 1; loop2 < MAX_TAIL; loop2++) {  
30.          pat = &fire[loop].particle[loop1][loop2];  
31.          pat->x = fire[loop].particle[loop1][0].x;  
32.          pat->y = fire[loop].particle[loop1][0].y;  
33.          pat->z = fire[loop].particle[loop1][0].z;  
34.       } //for loop2 end   
35.    }//for loop1 end   
36. }