如何在三个点间插值平滑粒子发射器的运动路径？

问题描述 #

我有一个移动的粒子发射器，每帧发射粒子，已通过**子帧插值（sub-frame interpolation）**实现粒子帧间过渡。发射器逐帧移动时，粒子发射路径呈线性，速度设为0时会得到预期的锯齿状路径。

我想通过当前路径点及前两个路径点这三个点来平滑路径，尝试了一段声称是**二次B样条插值（Quadratic BSpline Interpolation）**的代码，但输出完全混乱。我的世界空间范围为-1.0到1.0，原点在中心。需要解决：

• 是不是选了错误的方法？
• 若算法正确，可能存在哪些错误？
• 有没有更简单的三点平滑路径方法？

附尝试的代码：

//Quadratic BSpline Interpolation
// t = position between 0.0-1.0
double t2 = t * t;

// Calculate coefficients
double a0 = 0.5 * emitposlastx - emitposx + 0.5 * emitposlastx2;
double a1 = emitposx - emitposlastx;
double a2 = 0.5 * (emitposlastx2 + emitposlastx);

// Interpolate x coordinate
posx = a0 * t2 + a1 * t + a2;

a0 = 0.5 * emitposlasty - emitposy + 0.5 * emitposlasty2;
a1 = emitposy - emitposlasty;
a2 = 0.5 * (emitposlasty2 + emitposlasty);

// Interpolate y coordinate
posy = a0 * t2 + a1 * t + a2;

a0 = 0.5 * emitposlastz - emitposz + 0.5 * emitposlastz2;
a1 = emitposz - emitposlastz;
a2 = 0.5 * (emitposlastz2 + emitposlastz);

// Interpolate z coordinate
posz = a0 * t2 + a1 * t + a2;

问题分析 #

1. 当前二次B样条代码的核心错误 #

你的代码里的二次B样条公式完全错误，尤其是a1的计算。标准二次B样条的三个控制点（P0=前前帧位置、P1=前帧位置、P2=当前帧位置）对应的插值公式为：

P(t) = 0.5*[(1-t)²*P0 + (-2t²+2t+1)*P1 + t²*P2]

拆解为多项式形式后，正确的系数计算应该是：

• a0 = 0.5*P0 - P1 + 0.5*P2
• a1 = -P0 + P1
• a2 = 0.5*P0 + 0.5*P1

对比你的代码：

• 你把a1写成了emitposx - emitposlastx（也就是P2-P1），但正确的a1应该是emitposlastx - emitposlastx2（P1-P0），这直接导致了曲线方向和形态完全混乱。
• 另外要确认变量对应关系：emitposlastx2必须是P0，emitposlastx是P1，emitposx是P2，顺序不能搞反。

2. 其他可能的错误点 #

• t的取值逻辑错误：如果t是当前帧内的子帧进度（0到1），二次B样条的这段曲线是从(P0+P1)/2到(P1+P2)/2，而非直接连接P1到P2，这可能和你预期的帧间过渡不符。
• 变量赋值错误：比如把当前帧位置误设为前前帧位置，导致曲线方向完全偏离。
• t超出0-1范围：如果t的取值不小心超过1或者小于0，二次B样条会出现失控的延伸，超出你的世界空间范围。

替代解决方案 #

1. Catmull-Rom样条（推荐） #

Catmull-Rom样条会直接经过所有控制点，完美适配帧间路径平滑的需求，计算也简单。如果只有三个点（P0、P1、P2），可以用简化版的二次Catmull-Rom，或者补一个虚拟点（比如把P0复制为P1的前导点，P3复制为P2的后续点）。

简化版代码（基于三个点，t∈[0,1]对应从P1到P2的过渡）：

// 计算切线，用前后点的差值
double tangent1x = 0.5 * (emitposx - emitposlastx2); // P2-P0的一半
double tangent2x = tangent1x; // 没有下一帧的话，复用前切线

double t2 = t * t;
double t3 = t2 * t;

// Hermite插值实现Catmull-Rom
posx = (2*t3 - 3*t2 + 1) * emitposlastx + 
       (t3 - 2*t2 + t) * tangent1x + 
       (-2*t3 + 3*t2) * emitposx + 
       (t3 - t2) * tangent2x;

// y、z轴同理

2. 二次贝塞尔曲线 #

二次贝塞尔曲线只需要一个控制点，你可以用(emitposlastx2 + emitposlastx)/2（P0和P1的中点）作为控制点，连接P1到P2，实现平滑过渡：

double controlx = 0.5 * (emitposlastx2 + emitposlastx);
posx = (1-t)*(1-t)*emitposlastx + 2*(1-t)*t*controlx + t*t*emitposx;

// y、z轴同理

3. 简单加权平滑（快速方案） #

如果不需要严格的曲线，只是快速弱化锯齿，可以用动态加权平均：

• 当t=0（帧开始）：位置=P1
• 当t=0.5（帧中间）：位置=0.2P0 + 0.6P1 + 0.2*P2
• 当t=1（帧结束）：位置=P2

这种方法实现成本极低，但平滑效果不如样条曲线自然。

内容的提问来源于stack exchange，提问作者Rich95