用JavaScript玩转游戏物理(一)运动学模拟与粒子系统

系列简介

也许，三百年前的艾萨克·牛顿爵士(Sir Issac Newton, 1643-1727)并没幻想过，物理学广泛地应用在今天许多游戏、动画中。为什么在这些应用中要使用物理学？笔者认为，自我们出生以来，一直感受着物理世界的规律，意识到物体在这世界是如何"正常移动"，例如射球时球为抛物线(自旋的球可能会做成弧线球) 、石子系在一根线的末端会以固定频率摆动等等。要让游戏或动画中的物体有真实感，其移动方式就要符合我们对"正常移动"的预期。

今天的游戏动画应用了多种物理模拟技术，例如运动学模拟(kinematics simulation)、刚体动力学模拟(rigid body dynamics simulation)、绳子/布料模拟(string/cloth simulation)、柔体动力学模拟(soft body dynamics simulation)、流体动力学模拟(fluid dynamics simulation)等等。另外碰撞侦测(collision detection)是许多模拟系统里所需的。

本系列希望能介绍一些这方面最基础的知识，继续使用JavaScript做例子，以即时互动方式体验。

本文简介

作为系列第一篇，本文介绍最简单的运动学模拟，只有两条非常简单的公式。运动学模拟可以用来模拟很多物体运动(例如马里奥的跳跃、炮弹等)，本文将会配合粒子系统做出一些视觉特效(粒子系统其实也可以用来做游戏的玩法，而不单是视觉特效)。

运动学模拟

运动学(kinematics)研究物体的移动，和动力学(dynamics)不同之处，在于运动学不考虑物体的质量(mass)/转动惯量(moment of inertia)，以及不考虑加之于物体的力(force )和力矩(torque)。

我们先回忆牛顿第一运动定律:

当物体不受外力作用，或所受合力为零时，原先静止者恒静止，原先运动者恒沿着直线作等速度运动。该定律又称为「惯性定律」。

此定律指出，每个物体除了其位置(position)外，还有一个线性速度(linear velocity)的状态。然而，只模拟不受力影响的物体并不有趣。撇开力的概念，我们可以用线性加速度(linear acceleration)去影响物体的运动。例如，要计算一个自由落体在任意时间t的y轴座标，可以使用以下的分析解(analytical solution):

当中，和分别是t=0时的y轴起始座标和速度，而g则是重力加速度(gravitational acceleration)。

这分析解虽然简单，但是有一些缺点，例如g是常数，在模拟过程中不能改变；另外，当物体遇到障碍物，产生碰撞时，这公式也很难处理这种不连续性(discontinuity) 。

在计算机模拟中，通常需要计算连续的物体状态。用游戏的用语，就是计算第一帧的状态、第二帧的状态等等。设物体在任意时间t的状态：位置矢量为、速度矢量为、加速度矢量为。我们希望从时间的状态，计算下一个模拟时间的状态。最简单的方法，是采用欧拉方法(Euler method)作数值积分(numerical integration):

欧拉方法非常简单，但有准确度和稳定性问题，本文会先忽略这些问题。本文的例子采用二维空间，我们先实现一个JavaScript二维矢量类:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17

// Vector2.js Vector2 = function (x, y) { this .x = x; this .y = y; };   Vector2.prototype = {      copy : function () { return  new  Vector2( this .x, this .y); },      length : function () { return  Math.sqrt( this .x * this .x + this .y * this .y); },      sqrLength : function () { return  this .x * this .x + this .y * this .y; },      normalize : function () { var  inv = 1/ this .length(); return  new  Vector2( this .x * inv, this .y * inv); },      negate : function () { return  new  Vector2(- this .x, - this .y); },      add : function (v) { return  new  Vector2( this .x + v.x, this .y + v.y); },