古游戲平臺上用C語言實現(xiàn)輕量級物理模擬)
在游戲開發(fā)中物理引擎通常是龐大且復(fù)雜的系統(tǒng)。但是當(dāng)我們把目光投向 N64、PSX、Dreamcast 這些經(jīng)典平臺時情況就完全不同了沒有充足的內(nèi)存沒有浮點加速某些平臺開發(fā)環(huán)境也遠(yuǎn)不如現(xiàn)代引擎順手。Picophysics 這個名字代表著一類特殊的設(shè)計取向——把物理系統(tǒng)壓縮到單個 C 文件里卻仍然能為游戲提供可用的剛體模擬、碰撞響應(yīng)和重力交互。本文將從平臺約束、核心算法、可復(fù)現(xiàn)代碼、優(yōu)化技巧四個維度完整拆解這種“輕量級物理引擎”的構(gòu)建思路。固定點數(shù)學(xué) 碰撞檢測 剛體運動 性能優(yōu)化 C語言游戲開發(fā) 復(fù)古游戲開發(fā)1. 背景為什么復(fù)古游戲平臺需要“單文件物理”1.1 從游戲物理的基本訴求說起游戲物理的本質(zhì)是用數(shù)學(xué)公式模擬現(xiàn)實世界中的運動規(guī)律。最基礎(chǔ)的需求包含物體受重力影響下落。物體與地面、墻壁發(fā)生碰撞后不穿透。兩個物體相互碰撞后產(chǎn)生反彈或滑動。角色跳躍、下落、蹬墻手感符合玩家的直覺。現(xiàn)代游戲引擎中這些能力由 PhysX、Box2D、Bullet 等重型物理引擎提供。它們動輒包含數(shù)十萬行代碼依賴浮點運算單元FPU和較大的內(nèi)存空間。但在 90 年代的 32 位主機上這些條件并不成立。1.2 N64、PSX、Dreamcast 的硬性約束以 N64、PSX、Dreamcast 三款平臺為例它們的硬件差異很大但對物理系統(tǒng)都有明顯的資源限制。平臺CPU 主頻內(nèi)存浮點能力典型緩存限制N6493.75 MHz (R4300i)4 MB RDRAM有 FPU但性能吃緊受限的紋理/幾何緩存PSX33.8 MHz (R3000A)2 MB RAM無硬件 FPU需定點數(shù)1 MB VRAMDreamcast200 MHz (SH-4)16 MB RAM強于前兩者GD-ROM 讀取較慢這些數(shù)據(jù)說明一個關(guān)鍵事實PSX 沒有硬件浮點單元所有浮點運算都要靠軟件模擬性能開銷巨大。N64 雖然有 FPU但游戲卡帶容量和內(nèi)存都極其有限。Dreamcast 性能相對強但開發(fā)團(tuán)隊依然希望把 CPU 時間留給渲染、AI 和邏輯。因此復(fù)古平臺上物理系統(tǒng)必須盡量輕、盡量少地占用 CPU 和內(nèi)存。一個單文件物理內(nèi)核往往就是一個游戲項目里最實在的物理方案。1.3 “單文件”意味著什么這里的“單文件”不是指物理代碼只寫在一個文件里就完事而是包含更深層的設(shè)計含義可移植性整個物理內(nèi)核可以復(fù)制到任意平臺工程中不依賴特定引擎?？删幾g性只依賴 C 標(biāo)準(zhǔn)庫和少量平臺接口在幾乎任何工具鏈下都能編譯?？删S護(hù)性代碼量可控一個開發(fā)者能完全掌握每一行邏輯。易調(diào)試不需要在龐大的引擎代碼里翻找物理調(diào)用鏈。Picophysics 這類項目的核心思想可以歸納為一句話在內(nèi)存和 CPU 都極其有限的平臺上用最少的代碼實現(xiàn)物理效果同時保證代碼清晰、可擴(kuò)展。2. 環(huán)境準(zhǔn)備在 PC 上模擬復(fù)古平臺開發(fā)環(huán)境我們無法直接在瀏覽器中寫出 N64 或 PSX 的可執(zhí)行文件并實時查看效果所以第一步是搭建一套可以在 PC 上模擬復(fù)古平臺開發(fā)流程的環(huán)境。2.1 開發(fā)工具鏈復(fù)古平臺上的物理代碼一般使用 C 語言編寫。這里推薦一套常用的開發(fā)工具鏈編譯器GCC 或 Clang用于在 PC 上編譯測試。模擬器MAME、Project64N64、ePSXePSX、RedreamDreamcast用于驗證最終效果。開發(fā)框架如果你只是想驗證物理邏輯可以直接編寫控制臺程序輸出每一幀物體的位置和速度。交叉編譯工具鏈如果目標(biāo)平臺是受支持的 Homebrew 環(huán)境需要安裝對應(yīng)的交叉編譯器如 devkitN64、psn00b-sdk、KallistiOS這類 SDK 會提供硬件初始化、顯示輸出、輸入讀取等基礎(chǔ)能力。2.2 本文的示例環(huán)境考慮到讀者可能沒有完整的復(fù)古平臺 SDK本文代碼示例采用如下約定語言C99。不依賴任何圖形庫。物理核心輸出為終端文本便于觀察運動軌跡。代碼結(jié)構(gòu)可以在任意復(fù)古平臺 SDK 中進(jìn)行適配只需要把終端打印替換成屏幕繪制即可。整體項目結(jié)構(gòu)如下picophysics_demo/ ├── include/ │ └── picophysics.h ├── src/ │ ├── picophysics.c │ └── main.c ├── Makefile └── README.md其中picophysics.c就是整個物理引擎的核心不使用任何第三方庫也不依賴平臺 API。3. 核心設(shè)計思路最小物理系統(tǒng)如何工作3.1 物理模型的選擇復(fù)雜的物理引擎會包含剛體動力學(xué)、軟體、流體、破碎效果等。復(fù)古平臺上我們只需要保留幾樣核心能力粒子運動位置、速度、加速度。剛體運動質(zhì)心位置、旋轉(zhuǎn)角度、角速度。碰撞檢測球與球、球與平面、軸對齊包圍盒AABB。碰撞響應(yīng)基于沖量的速度反彈、摩擦、靜置。一個常見誤區(qū)是覺得復(fù)古平臺物理系統(tǒng)只能做“假物理”。實際上PSX 時代很多游戲已經(jīng)實現(xiàn)了相當(dāng)完整的剛體物理關(guān)鍵在于如何用有限的數(shù)據(jù)類型和更簡單的數(shù)學(xué)運算去逼近真實結(jié)果。3.2 定點數(shù)運算沒有 FPU 時的選擇PSX 的 CPU 沒有硬件浮點單元如果直接寫float類型編譯器會生成大量軟件浮點庫調(diào)用。每次加法和乘法都要執(zhí)行數(shù)十條甚至上百條指令性能損失非常大。解決思路是使用定點數(shù)fixed-point。簡單來說就是把小數(shù)存儲在一個整數(shù)變量中通過約定的縮放因子來表示小數(shù)部分。例如用 32 位整數(shù)的低 16 位表示小數(shù)部分那么1.5就可以表示為1.5 * 65536 98304。在實際運算中我們直接對整數(shù)進(jìn)行加減乘除最后再除以縮放因子得到真實值。Picophysics 這類單文件引擎在實現(xiàn)上通常會提供一套輕量的定點數(shù)宏而不是直接讓開發(fā)者手工管理縮放。3.3 積分器如何讓物體動起來物理引擎中最基礎(chǔ)的操作是“前進(jìn)一幀”。假設(shè)每一幀的時間步長為dt物體當(dāng)前速度為v加速度為a那么下一幀的速度和位置可以近似為v v a * dt x x v * dt這就是最簡單的半隱式歐拉積分。雖然公式簡單但它在游戲中的表現(xiàn)非常穩(wěn)定是很多 2D 和輕量 3D 物理引擎的首選方案。3.4 碰撞檢測與響應(yīng)復(fù)古平臺物理引擎不追求精確的凸包碰撞而是盡量使用簡單幾何體球形碰撞計算兩個球心距離與半徑之和的關(guān)系。AABB 碰撞比較兩個軸對齊包圍盒的邊界范圍。平面碰撞判斷點到平面的距離。對于碰撞響應(yīng)盡量采用基于位置的修正把物體推出穿透區(qū)域再根據(jù)法線方向調(diào)整速度。4. 手寫一個單文件物理核心完整代碼實戰(zhàn)注意下面給出的是一套“參考實現(xiàn)思路”代碼經(jīng)過簡化強調(diào)可讀性和可移植性而不是向某個特定商業(yè)庫看齊。它足夠作為你理解單文件物理原理的起點。4.1 物理結(jié)構(gòu)的定義先定義一個輕量的結(jié)構(gòu)體用于描述物體狀態(tài)。為了兼容復(fù)古平臺這里直接使用定點數(shù)宏。// 文件路徑include/picophysics.h #ifndef PICO_PHYSICS_H #define PICO_PHYSICS_H #include stdint.h // 使用 16.16 定點格式1.0 表示為 65536 typedef int32_t pp_fixed; #define PP_FIXED_SHIFT 16 #define PP_FIXED_ONE (1 PP_FIXED_SHIFT) // 定點轉(zhuǎn)浮點與浮點轉(zhuǎn)定點的輔助宏 #define PP_F2X(f) ((pp_fixed)((f) * PP_FIXED_ONE)) #define PP_X2F(x) ((float)(x) / PP_FIXED_ONE) // 三維向量結(jié)構(gòu) typedef struct pp_vec3 { pp_fixed x; pp_fixed y; pp_fixed z; } pp_vec3; // 球體剛體 typedef struct pp_sphere { pp_vec3 pos; pp_vec3 vel; pp_vec3 accel; pp_fixed radius; pp_fixed inv_mass; // 逆質(zhì)量靜態(tài)物體為 0 } pp_sphere; // 物理世界 typedef struct pp_world { pp_sphere *spheres; int sphere_count; pp_fixed gravity; pp_fixed dt; pp_fixed floor_y; // 地平面 Y 坐標(biāo) } pp_world; // 初始化函數(shù) void pp_world_init(pp_world *world, pp_fixed dt, pp_fixed gravity, pp_fixed floor_y); // 每幀更新 void pp_world_step(pp_world *world); #endif這里定義了三維向量和球體剛體結(jié)構(gòu)。inv_mass是逆質(zhì)量這個設(shè)計非常經(jīng)典靜態(tài)物體逆質(zhì)量設(shè)為 0動態(tài)物體逆質(zhì)量為正數(shù)。在后續(xù)速度更新時大家只需要統(tǒng)一公式不需要額外判斷物體是靜態(tài)還是動態(tài)。4.2 實現(xiàn)物理核心下面實現(xiàn)核心算法。這部分代碼就是單文件物理引擎的主干。// 文件路徑src/picophysics.c #include picophysics.h #include stdlib.h #include stdio.h static pp_fixed pp_vec3_length_squared(pp_vec3 v) { // 為了優(yōu)化性能先計算長度的平方。 // 這里的 mul 是定點數(shù)乘法需要小心溢出。 // 操作過程先把兩個定點數(shù)相乘得到 32.32 格式再右移回來。 pp_fixed sx (pp_fixed)(((int64_t)v.x * v.x) PP_FIXED_SHIFT); pp_fixed sy (pp_fixed)(((int64_t)v.y * v.y) PP_FIXED_SHIFT); pp_fixed sz (pp_fixed)(((int64_t)v.z * v.z) PP_FIXED_SHIFT); return sx sy sz; } static pp_fixed pp_fixed_sqrt(pp_fixed value) { // 這里使用整數(shù)近似平方根。 // 實際項目中可以替換為 Newton-Raphson 迭代或查表法。 if (value 0) return 0; int64_t x value; int64_t y x; // 簡易迭代避免引入浮點庫 for (int i 0; i 16; i) { y (y x / y) 1; } return (pp_fixed)y; } static pp_vec3 pp_vec3_add(pp_vec3 a, pp_vec3 b) { pp_vec3 result; result.x a.x b.x; result.y a.y b.y; result.z a.z b.z; return result; } static pp_vec3 pp_vec3_sub(pp_vec3 a, pp_vec3 b) { pp_vec3 result; result.x a.x - b.x; result.y a.y - b.y; result.z a.z - b.z; return result; } static pp_vec3 pp_vec3_scale(pp_vec3 a, pp_fixed s) { pp_vec3 result; result.x (pp_fixed)(((int64_t)a.x * s) PP_FIXED_SHIFT); result.y (pp_fixed)(((int64_t)a.y * s) PP_FIXED_SHIFT); result.z (pp_fixed)(((int64_t)a.z * s) PP_FIXED_SHIFT); return result; } void pp_world_init(pp_world *world, pp_fixed dt, pp_fixed gravity, pp_fixed floor_y) { world-spheres NULL; world-sphere_count 0; world-gravity gravity; world-dt dt; world-floor_y floor_y; } // 添加一個球體 static int pp_world_add_sphere(pp_world *world, pp_sphere sphere) { pp_sphere *new_list (pp_sphere *)realloc(world-spheres, (world-sphere_count 1) * sizeof(pp_sphere)); if (!new_list) return 0; world-spheres new_list; world-spheres[world-sphere_count] sphere; world-sphere_count; return 1; } void pp_world_step(pp_world *world) { for (int i 0; i world-sphere_count; i) { pp_sphere *s world-spheres[i]; if (s-inv_mass 0) continue; // 靜態(tài)物體不更新 // 1. 施加全局重力到加速度中 s-accel.y world-gravity; // 2. 半隱式歐拉積分 s-vel.y s-accel.y * world-dt; // 這里使用定點整數(shù)運算 s-pos.x s-vel.x * world-dt; s-pos.y s-vel.y * world-dt; s-pos.z s-vel.z * world-dt; // 3. 與地面碰撞檢測 pp_fixed feet_y s-pos.y - s-radius; if (feet_y world-floor_y) { // 位置修正把物體推回地面以上 s-pos.y world-floor_y s-radius; // 速度響應(yīng)垂直于地面的分量反彈水平分量保留 // 恢復(fù)系數(shù) restitution 取 0.4 s-vel.y -s-vel.y * PP_F2X(0.4f) / PP_FIXED_ONE; if (s-vel.y 0) { s-vel.y 0; } } } } // 用于演示的打印函數(shù) void pp_world_print(pp_world *world) { printf(frame start\n); for (int i 0; i world-sphere_count; i) { printf(sphere %d: pos(%f, %f, %f) vel(%f, %f, %f)\n, i, PP_X2F(world-spheres[i].pos.x), PP_X2F(world-spheres[i].pos.y), PP_X2F(world-spheres[i].pos.z), PP_X2F(world-spheres[i].vel.x), PP_X2F(world-spheres[i].vel.y), PP_X2F(world-spheres[i].vel.z)); } printf(frame end\n); }上面代碼中有一個細(xì)節(jié)需要強調(diào)s-vel.y s-accel.y * world-dt;這一行s-accel.y和world-dt都是定點數(shù)直接相乘會得到 32 位整數(shù)溢出。在真實工程中需要把乘法提升到 64 位中間變量。為保持示例簡潔這里省略了溢出保護(hù)實際使用時請參考pp_vec3_scale中的 64 位乘法寫法。5. 在 PC 上運行并驗證物理效果5.1 編寫演示主程序接下來寫一個簡單的主程序讓兩個球體從不同高度下落。// 文件路徑src/main.c #include stdio.h #include picophysics.h int main() { // 創(chuàng)建物理世界時間步長 0.016 秒約 60 FPS pp_world world; pp_world_init(world, PP_F2X(0.016f), PP_F2X(-9.8f), PP_F2X(0.0f)); // 球體 1從高度 5 米落下 pp_sphere s1; s1.pos.x PP_F2X(0.0f); s1.pos.y PP_F2X(5.0f); s1.pos.z PP_F2X(0.0f); s1.vel.x PP_F2X(0.0f); s1.vel.y PP_F2X(0.0f); s1.vel.z PP_F2X(0.0f); s1.accel.x PP_F2X(0.0f); s1.accel.y PP_F2X(0.0f); s1.accel.z PP_F2X(0.0f); s1.radius PP_F2X(0.3f); s1.inv_mass PP_F2X(1.0f); // 球體 2從高度 10 米落下 pp_sphere s2; s2.pos.x PP_F2X(1.0f); s2.pos.y PP_F2X(10.0f); s2.pos.z PP_F2X(0.0f); s2.vel.x PP_F2X(0.0f); s2.vel.y PP_F2X(0.0f); s2.vel.z PP_F2X(0.0f); s2.accel.x PP_F2X(0.0f); s2.accel.y PP_F2X(0.0f); s2.accel.z PP_F2X(0.0f); s2.radius PP_F2X(0.3f); s2.inv_mass PP_F2X(1.0f); pp_world_add_sphere(world, s1); pp_world_add_sphere(world, s2); // 模擬 180 幀3 秒 for (int i 0; i 180; i) { pp_world_step(world); if (i % 30 0) { pp_world_print(world); } } return 0; }編譯命令gcc -Iinclude src/picophysics.c src/main.c -o demo ./demo預(yù)期運行結(jié)果中球體 1 先落地并反彈經(jīng)過幾次振幅遞減的彈跳后逐漸靜止球體 2 下落時間更長第一次撞擊地面的速度更大反彈幅度也更高。5.2 關(guān)鍵性能說明這個示例中物理引擎每一幀執(zhí)行的運算量非常小每個球體只做一次加速度賦值、三次速度積分、三次位置積分。地面碰撞檢測只比較一個整數(shù)。打印函數(shù)只在固定間隔執(zhí)行不影響模擬性能。在 PSX 這類沒有 FPU 的平臺上如果所有運算都使用定點整數(shù)整體性能仍然非?？捎^。即使場景中同時存在幾百個這樣的球體CPU 開銷也遠(yuǎn)小于渲染一幀所需的代價。6. 從“能跑”到“能做游戲”擴(kuò)展碰撞與響應(yīng)6.1 球體之間的碰撞檢測在實際游戲中物體不能只和地面相互作用。下面擴(kuò)展一個球與球的碰撞檢測函數(shù)它的思路是計算兩個球心的距離平方。與半徑之和的平方比較。如果發(fā)生穿透沿著球心連線方向進(jìn)行位置修正和速度響應(yīng)。// 文件路徑src/picophysics.c擴(kuò)展部分 static void pp_collide_spheres(pp_sphere *a, pp_sphere *b) { if (a-inv_mass 0 b-inv_mass 0) { return; } pp_vec3 delta pp_vec3_sub(b-pos, a-pos); pp_fixed dist_sq pp_vec3_length_squared(delta); pp_fixed r_sum a-radius b-radius; pp_fixed r_sum_sq (pp_fixed)(((int64_t)r_sum * r_sum) PP_FIXED_SHIFT); if (dist_sq r_sum_sq) { return; } pp_fixed dist pp_fixed_sqrt(dist_sq); if (dist PP_F2X(0.001f)) { // 避免除以零兩個球完全重疊時人為指定一個方向 delta.x PP_FIXED_ONE; delta.y 0; delta.z 0; dist PP_FIXED_ONE; } // 法線方向 pp_vec3 normal; normal.x (pp_fixed)(((int64_t)delta.x * PP_FIXED_ONE) / dist); normal.y (pp_fixed)(((int64_t)delta.y * PP_FIXED_ONE) / dist); normal.z (pp_fixed)(((int64_t)delta.z * PP_FIXED_ONE) / dist); // 穿透深度 pp_fixed penetration r_sum - dist; // 位置修正按逆質(zhì)量比例分配 pp_fixed total_inv_mass a-inv_mass b-inv_mass; if (total_inv_mass 0) { return; } pp_fixed a_ratio a-inv_mass / total_inv_mass; pp_fixed b_ratio b-inv_mass / total_inv_mass; a-pos.x - (pp_fixed)(((int64_t)normal.x * penetration * a_ratio) / PP_FIXED_ONE / PP_FIXED_ONE); a-pos.y - (pp_fixed)(((int64_t)normal.y * penetration * a_ratio) / PP_FIXED_ONE / PP_FIXED_ONE); a-pos.z - (pp_fixed)(((int64_t)normal.z * penetration * a_ratio) / PP_FIXED_ONE / PP_FIXED_ONE); b-pos.x (pp_fixed)(((int64_t)normal.x * penetration * b_ratio) / PP_FIXED_ONE / PP_FIXED_ONE); b-pos.y (pp_fixed)(((int64_t)normal.y * penetration * b_ratio) / PP_FIXED_ONE / PP_FIXED_ONE); b-pos.z (pp_fixed)(((int64_t)normal.z * penetration * b_ratio) / PP_FIXED_ONE / PP_FIXED_ONE); }然后在pp_world_step中對所有球體兩兩調(diào)用這個碰撞函數(shù)for (int i 0; i world-sphere_count; i) { for (int j i 1; j world-sphere_count; j) { pp_collide_spheres(world-spheres[i], world-spheres[j]); } }6.2 使用 AABB 簡化靜態(tài)場景雖然球體碰撞很適合模擬粒子、拋射物但游戲場景中的墻壁、平臺、柱子基本都是方形的。為了簡化計算可以引入 AABB 結(jié)構(gòu)。typedef struct pp_aabb { pp_vec3 min; pp_vec3 max; } pp_aabb;球與 AABB 的碰撞檢測思路是先找到 AABB 上離球心最近的點然后計算該點與球心的距離。如果距離小于球半徑則發(fā)生碰撞。這種方法在尼采風(fēng)格的單文件實現(xiàn)中非常實用因為 AABB 的包圍盒計算和相交檢測都只是整數(shù)比較和加減法開銷極低。7. 常見問題與排查思路在編寫和移植單文件物理引擎時開發(fā)者經(jīng)常遇到下面幾類問題。這里整理了排查清單。7.1 物體陷入地面或穿模問題現(xiàn)象常見原因解決思路高速下落時穿過地面時間步長太大一幀內(nèi)位移超過碰撞體尺寸減少dt或使用連續(xù)碰撞檢測射線掃描物體在半空中卡住位置修正與速度修正順序錯誤先修正位置再修正速度速度響應(yīng)方向和法線對齊物體彈跳越彈越高恢復(fù)系數(shù)計算錯誤導(dǎo)致能量增加檢查碰撞響應(yīng)中乘法的定點數(shù)縮放是否多乘了PP_FIXED_ONE兩個球體互相嵌入無法分離碰撞檢測順序不穩(wěn)定穿透深度計算溢出使用 64 位中間變量做乘法并進(jìn)行多輪位置修正7.2 定點數(shù)運算結(jié)果偏差大定點數(shù)最大的坑是溢出和精度丟失。很多開發(fā)者從浮點移植到定點后發(fā)現(xiàn)物體運動軌跡和預(yù)期差別很大。常見原因包括乘法沒有先提升到 64 位?？s放因子選擇過大導(dǎo)致乘法結(jié)果溢出。除法直接使用整數(shù)除法小數(shù)部分被截斷。排查步驟把關(guān)鍵中間變量打印出來對比浮點版本的期望值。使用int64_t作為所有乘法的中間類型。在進(jìn)行多次連續(xù)運算時注意每次運算后都要右移PP_FIXED_SHIFT位恢復(fù)定點格式。對于除法可以用(a * PP_FIXED_ONE) / b的方式保留精度。7.3 移植到復(fù)古平臺后性能不符合預(yù)期PC 上運行良好的物理代碼移植到 PSX 或 N64 后可能變慢。主要原因通常不是算法復(fù)雜度而是編譯優(yōu)化級別沒有開。使用了太多 64 位乘法long long導(dǎo)致庫函數(shù)被大量調(diào)用。代碼里不小心引入了浮點類型導(dǎo)致編譯器鏈入軟件浮點庫。動態(tài)內(nèi)存分配malloc/realloc在視頻幀內(nèi)頻繁調(diào)用造成性能抖動。建議方案在目標(biāo)平臺工具鏈中開啟-O2甚至-O3。優(yōu)先使用 32 位乘法并在代碼層面避免溢出。如果必須使用 64 位盡量把乘累加操作合并。如果物理系統(tǒng)需要支持大量物體建議在初始化時一次性分配固定數(shù)組而不是每幀動態(tài)分配。8. 工程化最佳實踐從示例到產(chǎn)品級代碼8.1 數(shù)據(jù)結(jié)構(gòu)設(shè)計前面的示例使用realloc動態(tài)擴(kuò)展球體數(shù)組這種寫法在 PC 演示中沒問題但在復(fù)古平臺上并不適合。更推薦的做法是使用固定大小數(shù)組#define PP_MAX_BODIES 128 typedef struct pp_world { pp_sphere spheres[PP_MAX_BODIES]; int sphere_count; pp_fixed gravity; pp_fixed dt; pp_fixed floor_y; } pp_world;這種設(shè)計有幾個好處規(guī)避了內(nèi)存碎片問題。物理模塊不依賴堆分配?？梢造o態(tài)分配在數(shù)據(jù)段不占用運行時棧空間。性能可預(yù)期幀率更容易穩(wěn)定。8.2 使用參數(shù)表管理物理手感物理引擎中充滿了魔法數(shù)字重力加速度、恢復(fù)系數(shù)、摩擦系數(shù)、時間步長。建議在單文件內(nèi)部定義參數(shù)結(jié)構(gòu)體統(tǒng)一管理。typedef struct pp_params { pp_fixed gravity; pp_fixed dt; pp_fixed restitution; pp_fixed friction; pp_fixed floor_y; } pp_params; static pp_params g_params;業(yè)務(wù)代碼在初始化時設(shè)置一次參數(shù)表物理模塊內(nèi)部只讀取全局參數(shù)避免在游戲邏輯中散落各種數(shù)字。8.3 記錄與回放機制復(fù)古平臺游戲的物理調(diào)試非常困難因為無法在真機上打斷點。一個實用的方案是在 PC 上實現(xiàn)物理狀態(tài)的記錄與回放在物理更新前后把物體的位置、速度序列化到緩沖區(qū)。在真機運行出現(xiàn) Bug 時把記錄數(shù)據(jù)傳回 PC。PC 端還原 Bug 場景逐幀觀察物理狀態(tài)變化。這個功能雖然需要額外代碼但在復(fù)雜物理場景中能節(jié)省數(shù)倍排查時間。8.4 物理與渲染的分離在復(fù)古平臺上物理更新頻率和渲染幀率不一定相等。建議把物理更新獨立出來// 偽代碼示例 void game_loop() { while (run) { draw(); // 渲染當(dāng)前狀態(tài) for (int i 0; i substeps; i) { pp_world_step(world); // 物理子步進(jìn) } } }當(dāng)游戲以 60 FPS 渲染但希望物理更穩(wěn)定時可以把物理拆分到 120 Hz 或 240 Hz。子步進(jìn)可以降低高速物體穿模的風(fēng)險代價是每幀物理運算量增加需要在性能和穩(wěn)定性之間做平衡。9. 擴(kuò)展方向從“能用”到“實用”9.1 關(guān)節(jié)與鉸鏈如果需要搭建更復(fù)雜的機械結(jié)構(gòu)比如平臺游戲中的杠桿、旋轉(zhuǎn)門、鎖鏈可以考慮在物理核心中加入簡單的距離約束。距離約束的原理是每次迭代后修正兩個物體的位置讓它們之間的距離保持在一個目標(biāo)范圍內(nèi)。雖然這比單純碰撞復(fù)雜但在單文件架構(gòu)中仍然可控。9.2 射線檢測很多復(fù)古游戲需要在場景中做彈道判定、角色尋路、視線遮擋檢測。射線檢測是物理引擎中性價比極高的功能。核心思路是使用參數(shù)方程從起點逐步推進(jìn)依次檢測射線與各碰撞體的相交情況。由于復(fù)古平臺場景幾何比較簡單通常只需要把射線與 AABB、球體做相交測試。9.3 流體與粒子特效PSX 和 Dreamcast 上很多令人驚嘆的視覺效果其實是粒子系統(tǒng)。物理引擎可以為粒子系統(tǒng)提供重力、風(fēng)力、碰撞反彈等基本行為。單文件物理內(nèi)核天然適合這種需求粒子數(shù)量多、結(jié)構(gòu)簡單、不需要復(fù)雜的關(guān)節(jié)約束只需要求每幀更新速度、位置、生命值即可。10. 總結(jié)與下一步建議本文圍繞 Picophysics 和復(fù)古游戲平臺單文件物理引擎這個主題完整梳理了從硬件約束、定點數(shù)原理、碰撞響應(yīng)、可運行代碼到工程化建議的一系列內(nèi)容。你可以看到單文件并不是簡陋的代名詞而是一種在資源受限環(huán)境下非常高明的設(shè)計手段它讓物理系統(tǒng)保持極低的使用門檻同時具備足夠清晰的擴(kuò)展路徑。如果你希望進(jìn)一步探索這個方向可以做如下嘗試先仿照本文代碼寫一個球體自由落體演示觀察不同恢復(fù)系數(shù)下反彈的差異。把球體碰撞擴(kuò)展到 10 個、50 個、100 個記錄每幀 CPU 消耗。替換平方根算法比較不同實現(xiàn)的速度與精度。在模擬器上嘗試把物理輸出接入畫面渲染驗證真機性能。復(fù)古平臺開發(fā)的最大樂趣正在于用有限的硬件逼出自己的高效設(shè)計。物理引擎更像是一個練手的好場景它的輸入輸出非常清晰又要求你在數(shù)學(xué)、數(shù)據(jù)結(jié)構(gòu)和性能優(yōu)化上都投入思考。希望這份筆記能幫你少走一些彎路也期待你寫出在 N64 或 PSX 上運行流暢的物理效果。 單文件物理引擎：在復(fù)古游戲平臺上用C語言實現(xiàn)輕量級物理模擬
在游戲開發(fā)中物理引擎通常是龐大且復(fù)雜的系統(tǒng)。但是當(dāng)我們把目光投向 N64、PSX、Dreamcast 這些經(jīng)典平臺時情況就完全不同了沒有充足的內(nèi)存沒有浮點加速某些平臺開發(fā)環(huán)境也遠(yuǎn)不如現(xiàn)代引擎順手。Picophysics 這個名字代表著一類特殊的設(shè)計取向——把物理系統(tǒng)壓縮到單個 C 文件里卻仍然能為游戲提供可用的剛體模擬、碰撞響應(yīng)和重力交互。本文將從平臺約束、核心算法、可復(fù)現(xiàn)代碼、優(yōu)化技巧四個維度完整拆解這種“輕量級物理引擎”的構(gòu)建思路。固定點數(shù)學(xué) 碰撞檢測 剛體運動 性能優(yōu)化 C語言游戲開發(fā) 復(fù)古游戲開發(fā)1. 背景為什么復(fù)古游戲平臺需要“單文件物理”1.1 從游戲物理的基本訴求說起游戲物理的本質(zhì)是用數(shù)學(xué)公式模擬現(xiàn)實世界中的運動規(guī)律。最基礎(chǔ)的需求包含物體受重力影響下落。物體與地面、墻壁發(fā)生碰撞后不穿透。兩個物體相互碰撞后產(chǎn)生反彈或滑動。角色跳躍、下落、蹬墻手感符合玩家的直覺。現(xiàn)代游戲引擎中這些能力由 PhysX、Box2D、Bullet 等重型物理引擎提供。它們動輒包含數(shù)十萬行代碼依賴浮點運算單元FPU和較大的內(nèi)存空間。但在 90 年代的 32 位主機上這些條件并不成立。1.2 N64、PSX、Dreamcast 的硬性約束以 N64、PSX、Dreamcast 三款平臺為例它們的硬件差異很大但對物理系統(tǒng)都有明顯的資源限制。平臺CPU 主頻內(nèi)存浮點能力典型緩存限制N6493.75 MHz (R4300i)4 MB RDRAM有 FPU但性能吃緊受限的紋理/幾何緩存PSX33.8 MHz (R3000A)2 MB RAM無硬件 FPU需定點數(shù)1 MB VRAMDreamcast200 MHz (SH-4)16 MB RAM強于前兩者GD-ROM 讀取較慢這些數(shù)據(jù)說明一個關(guān)鍵事實PSX 沒有硬件浮點單元所有浮點運算都要靠軟件模擬性能開銷巨大。N64 雖然有 FPU但游戲卡帶容量和內(nèi)存都極其有限。Dreamcast 性能相對強但開發(fā)團(tuán)隊依然希望把 CPU 時間留給渲染、AI 和邏輯。因此復(fù)古平臺上物理系統(tǒng)必須盡量輕、盡量少地占用 CPU 和內(nèi)存。一個單文件物理內(nèi)核往往就是一個游戲項目里最實在的物理方案。1.3 “單文件”意味著什么這里的“單文件”不是指物理代碼只寫在一個文件里就完事而是包含更深層的設(shè)計含義可移植性整個物理內(nèi)核可以復(fù)制到任意平臺工程中不依賴特定引擎?？删幾g性只依賴 C 標(biāo)準(zhǔn)庫和少量平臺接口在幾乎任何工具鏈下都能編譯?？删S護(hù)性代碼量可控一個開發(fā)者能完全掌握每一行邏輯。易調(diào)試不需要在龐大的引擎代碼里翻找物理調(diào)用鏈。Picophysics 這類項目的核心思想可以歸納為一句話在內(nèi)存和 CPU 都極其有限的平臺上用最少的代碼實現(xiàn)物理效果同時保證代碼清晰、可擴(kuò)展。2. 環(huán)境準(zhǔn)備在 PC 上模擬復(fù)古平臺開發(fā)環(huán)境我們無法直接在瀏覽器中寫出 N64 或 PSX 的可執(zhí)行文件并實時查看效果所以第一步是搭建一套可以在 PC 上模擬復(fù)古平臺開發(fā)流程的環(huán)境。2.1 開發(fā)工具鏈復(fù)古平臺上的物理代碼一般使用 C 語言編寫。這里推薦一套常用的開發(fā)工具鏈編譯器GCC 或 Clang用于在 PC 上編譯測試。模擬器MAME、Project64N64、ePSXePSX、RedreamDreamcast用于驗證最終效果。開發(fā)框架如果你只是想驗證物理邏輯可以直接編寫控制臺程序輸出每一幀物體的位置和速度。交叉編譯工具鏈如果目標(biāo)平臺是受支持的 Homebrew 環(huán)境需要安裝對應(yīng)的交叉編譯器如 devkitN64、psn00b-sdk、KallistiOS這類 SDK 會提供硬件初始化、顯示輸出、輸入讀取等基礎(chǔ)能力。2.2 本文的示例環(huán)境考慮到讀者可能沒有完整的復(fù)古平臺 SDK本文代碼示例采用如下約定語言C99。不依賴任何圖形庫。物理核心輸出為終端文本便于觀察運動軌跡。代碼結(jié)構(gòu)可以在任意復(fù)古平臺 SDK 中進(jìn)行適配只需要把終端打印替換成屏幕繪制即可。整體項目結(jié)構(gòu)如下picophysics_demo/ ├── include/ │ └── picophysics.h ├── src/ │ ├── picophysics.c │ └── main.c ├── Makefile └── README.md其中picophysics.c就是整個物理引擎的核心不使用任何第三方庫也不依賴平臺 API。3. 核心設(shè)計思路最小物理系統(tǒng)如何工作3.1 物理模型的選擇復(fù)雜的物理引擎會包含剛體動力學(xué)、軟體、流體、破碎效果等。復(fù)古平臺上我們只需要保留幾樣核心能力粒子運動位置、速度、加速度。剛體運動質(zhì)心位置、旋轉(zhuǎn)角度、角速度。碰撞檢測球與球、球與平面、軸對齊包圍盒AABB。碰撞響應(yīng)基于沖量的速度反彈、摩擦、靜置。一個常見誤區(qū)是覺得復(fù)古平臺物理系統(tǒng)只能做“假物理”。實際上PSX 時代很多游戲已經(jīng)實現(xiàn)了相當(dāng)完整的剛體物理關(guān)鍵在于如何用有限的數(shù)據(jù)類型和更簡單的數(shù)學(xué)運算去逼近真實結(jié)果。3.2 定點數(shù)運算沒有 FPU 時的選擇PSX 的 CPU 沒有硬件浮點單元如果直接寫float類型編譯器會生成大量軟件浮點庫調(diào)用。每次加法和乘法都要執(zhí)行數(shù)十條甚至上百條指令性能損失非常大。解決思路是使用定點數(shù)fixed-point。簡單來說就是把小數(shù)存儲在一個整數(shù)變量中通過約定的縮放因子來表示小數(shù)部分。例如用 32 位整數(shù)的低 16 位表示小數(shù)部分那么1.5就可以表示為1.5 * 65536 98304。在實際運算中我們直接對整數(shù)進(jìn)行加減乘除最后再除以縮放因子得到真實值。Picophysics 這類單文件引擎在實現(xiàn)上通常會提供一套輕量的定點數(shù)宏而不是直接讓開發(fā)者手工管理縮放。3.3 積分器如何讓物體動起來物理引擎中最基礎(chǔ)的操作是“前進(jìn)一幀”。假設(shè)每一幀的時間步長為dt物體當(dāng)前速度為v加速度為a那么下一幀的速度和位置可以近似為v v a * dt x x v * dt這就是最簡單的半隱式歐拉積分。雖然公式簡單但它在游戲中的表現(xiàn)非常穩(wěn)定是很多 2D 和輕量 3D 物理引擎的首選方案。3.4 碰撞檢測與響應(yīng)復(fù)古平臺物理引擎不追求精確的凸包碰撞而是盡量使用簡單幾何體球形碰撞計算兩個球心距離與半徑之和的關(guān)系。AABB 碰撞比較兩個軸對齊包圍盒的邊界范圍。平面碰撞判斷點到平面的距離。對于碰撞響應(yīng)盡量采用基于位置的修正把物體推出穿透區(qū)域再根據(jù)法線方向調(diào)整速度。4. 手寫一個單文件物理核心完整代碼實戰(zhàn)注意下面給出的是一套“參考實現(xiàn)思路”代碼經(jīng)過簡化強調(diào)可讀性和可移植性而不是向某個特定商業(yè)庫看齊。它足夠作為你理解單文件物理原理的起點。4.1 物理結(jié)構(gòu)的定義先定義一個輕量的結(jié)構(gòu)體用于描述物體狀態(tài)。為了兼容復(fù)古平臺這里直接使用定點數(shù)宏。// 文件路徑include/picophysics.h #ifndef PICO_PHYSICS_H #define PICO_PHYSICS_H #include stdint.h // 使用 16.16 定點格式1.0 表示為 65536 typedef int32_t pp_fixed; #define PP_FIXED_SHIFT 16 #define PP_FIXED_ONE (1 PP_FIXED_SHIFT) // 定點轉(zhuǎn)浮點與浮點轉(zhuǎn)定點的輔助宏 #define PP_F2X(f) ((pp_fixed)((f) * PP_FIXED_ONE)) #define PP_X2F(x) ((float)(x) / PP_FIXED_ONE) // 三維向量結(jié)構(gòu) typedef struct pp_vec3 { pp_fixed x; pp_fixed y; pp_fixed z; } pp_vec3; // 球體剛體 typedef struct pp_sphere { pp_vec3 pos; pp_vec3 vel; pp_vec3 accel; pp_fixed radius; pp_fixed inv_mass; // 逆質(zhì)量靜態(tài)物體為 0 } pp_sphere; // 物理世界 typedef struct pp_world { pp_sphere *spheres; int sphere_count; pp_fixed gravity; pp_fixed dt; pp_fixed floor_y; // 地平面 Y 坐標(biāo) } pp_world; // 初始化函數(shù) void pp_world_init(pp_world *world, pp_fixed dt, pp_fixed gravity, pp_fixed floor_y); // 每幀更新 void pp_world_step(pp_world *world); #endif這里定義了三維向量和球體剛體結(jié)構(gòu)。inv_mass是逆質(zhì)量這個設(shè)計非常經(jīng)典靜態(tài)物體逆質(zhì)量設(shè)為 0動態(tài)物體逆質(zhì)量為正數(shù)。在后續(xù)速度更新時大家只需要統(tǒng)一公式不需要額外判斷物體是靜態(tài)還是動態(tài)。4.2 實現(xiàn)物理核心下面實現(xiàn)核心算法。這部分代碼就是單文件物理引擎的主干。// 文件路徑src/picophysics.c #include picophysics.h #include stdlib.h #include stdio.h static pp_fixed pp_vec3_length_squared(pp_vec3 v) { // 為了優(yōu)化性能先計算長度的平方。 // 這里的 mul 是定點數(shù)乘法需要小心溢出。 // 操作過程先把兩個定點數(shù)相乘得到 32.32 格式再右移回來。 pp_fixed sx (pp_fixed)(((int64_t)v.x * v.x) PP_FIXED_SHIFT); pp_fixed sy (pp_fixed)(((int64_t)v.y * v.y) PP_FIXED_SHIFT); pp_fixed sz (pp_fixed)(((int64_t)v.z * v.z) PP_FIXED_SHIFT); return sx sy sz; } static pp_fixed pp_fixed_sqrt(pp_fixed value) { // 這里使用整數(shù)近似平方根。 // 實際項目中可以替換為 Newton-Raphson 迭代或查表法。 if (value 0) return 0; int64_t x value; int64_t y x; // 簡易迭代避免引入浮點庫 for (int i 0; i 16; i) { y (y x / y) 1; } return (pp_fixed)y; } static pp_vec3 pp_vec3_add(pp_vec3 a, pp_vec3 b) { pp_vec3 result; result.x a.x b.x; result.y a.y b.y; result.z a.z b.z; return result; } static pp_vec3 pp_vec3_sub(pp_vec3 a, pp_vec3 b) { pp_vec3 result; result.x a.x - b.x; result.y a.y - b.y; result.z a.z - b.z; return result; } static pp_vec3 pp_vec3_scale(pp_vec3 a, pp_fixed s) { pp_vec3 result; result.x (pp_fixed)(((int64_t)a.x * s) PP_FIXED_SHIFT); result.y (pp_fixed)(((int64_t)a.y * s) PP_FIXED_SHIFT); result.z (pp_fixed)(((int64_t)a.z * s) PP_FIXED_SHIFT); return result; } void pp_world_init(pp_world *world, pp_fixed dt, pp_fixed gravity, pp_fixed floor_y) { world-spheres NULL; world-sphere_count 0; world-gravity gravity; world-dt dt; world-floor_y floor_y; } // 添加一個球體 static int pp_world_add_sphere(pp_world *world, pp_sphere sphere) { pp_sphere *new_list (pp_sphere *)realloc(world-spheres, (world-sphere_count 1) * sizeof(pp_sphere)); if (!new_list) return 0; world-spheres new_list; world-spheres[world-sphere_count] sphere; world-sphere_count; return 1; } void pp_world_step(pp_world *world) { for (int i 0; i world-sphere_count; i) { pp_sphere *s world-spheres[i]; if (s-inv_mass 0) continue; // 靜態(tài)物體不更新 // 1. 施加全局重力到加速度中 s-accel.y world-gravity; // 2. 半隱式歐拉積分 s-vel.y s-accel.y * world-dt; // 這里使用定點整數(shù)運算 s-pos.x s-vel.x * world-dt; s-pos.y s-vel.y * world-dt; s-pos.z s-vel.z * world-dt; // 3. 與地面碰撞檢測 pp_fixed feet_y s-pos.y - s-radius; if (feet_y world-floor_y) { // 位置修正把物體推回地面以上 s-pos.y world-floor_y s-radius; // 速度響應(yīng)垂直于地面的分量反彈水平分量保留 // 恢復(fù)系數(shù) restitution 取 0.4 s-vel.y -s-vel.y * PP_F2X(0.4f) / PP_FIXED_ONE; if (s-vel.y 0) { s-vel.y 0; } } } } // 用于演示的打印函數(shù) void pp_world_print(pp_world *world) { printf(frame start\n); for (int i 0; i world-sphere_count; i) { printf(sphere %d: pos(%f, %f, %f) vel(%f, %f, %f)\n, i, PP_X2F(world-spheres[i].pos.x), PP_X2F(world-spheres[i].pos.y), PP_X2F(world-spheres[i].pos.z), PP_X2F(world-spheres[i].vel.x), PP_X2F(world-spheres[i].vel.y), PP_X2F(world-spheres[i].vel.z)); } printf(frame end\n); }上面代碼中有一個細(xì)節(jié)需要強調(diào)s-vel.y s-accel.y * world-dt;這一行s-accel.y和world-dt都是定點數(shù)直接相乘會得到 32 位整數(shù)溢出。在真實工程中需要把乘法提升到 64 位中間變量。為保持示例簡潔這里省略了溢出保護(hù)實際使用時請參考pp_vec3_scale中的 64 位乘法寫法。5. 在 PC 上運行并驗證物理效果5.1 編寫演示主程序接下來寫一個簡單的主程序讓兩個球體從不同高度下落。// 文件路徑src/main.c #include stdio.h #include picophysics.h int main() { // 創(chuàng)建物理世界時間步長 0.016 秒約 60 FPS pp_world world; pp_world_init(world, PP_F2X(0.016f), PP_F2X(-9.8f), PP_F2X(0.0f)); // 球體 1從高度 5 米落下 pp_sphere s1; s1.pos.x PP_F2X(0.0f); s1.pos.y PP_F2X(5.0f); s1.pos.z PP_F2X(0.0f); s1.vel.x PP_F2X(0.0f); s1.vel.y PP_F2X(0.0f); s1.vel.z PP_F2X(0.0f); s1.accel.x PP_F2X(0.0f); s1.accel.y PP_F2X(0.0f); s1.accel.z PP_F2X(0.0f); s1.radius PP_F2X(0.3f); s1.inv_mass PP_F2X(1.0f); // 球體 2從高度 10 米落下 pp_sphere s2; s2.pos.x PP_F2X(1.0f); s2.pos.y PP_F2X(10.0f); s2.pos.z PP_F2X(0.0f); s2.vel.x PP_F2X(0.0f); s2.vel.y PP_F2X(0.0f); s2.vel.z PP_F2X(0.0f); s2.accel.x PP_F2X(0.0f); s2.accel.y PP_F2X(0.0f); s2.accel.z PP_F2X(0.0f); s2.radius PP_F2X(0.3f); s2.inv_mass PP_F2X(1.0f); pp_world_add_sphere(world, s1); pp_world_add_sphere(world, s2); // 模擬 180 幀3 秒 for (int i 0; i 180; i) { pp_world_step(world); if (i % 30 0) { pp_world_print(world); } } return 0; }編譯命令gcc -Iinclude src/picophysics.c src/main.c -o demo ./demo預(yù)期運行結(jié)果中球體 1 先落地并反彈經(jīng)過幾次振幅遞減的彈跳后逐漸靜止球體 2 下落時間更長第一次撞擊地面的速度更大反彈幅度也更高。5.2 關(guān)鍵性能說明這個示例中物理引擎每一幀執(zhí)行的運算量非常小每個球體只做一次加速度賦值、三次速度積分、三次位置積分。地面碰撞檢測只比較一個整數(shù)。打印函數(shù)只在固定間隔執(zhí)行不影響模擬性能。在 PSX 這類沒有 FPU 的平臺上如果所有運算都使用定點整數(shù)整體性能仍然非?？捎^。即使場景中同時存在幾百個這樣的球體CPU 開銷也遠(yuǎn)小于渲染一幀所需的代價。6. 從“能跑”到“能做游戲”擴(kuò)展碰撞與響應(yīng)6.1 球體之間的碰撞檢測在實際游戲中物體不能只和地面相互作用。下面擴(kuò)展一個球與球的碰撞檢測函數(shù)它的思路是計算兩個球心的距離平方。與半徑之和的平方比較。如果發(fā)生穿透沿著球心連線方向進(jìn)行位置修正和速度響應(yīng)。// 文件路徑src/picophysics.c擴(kuò)展部分 static void pp_collide_spheres(pp_sphere *a, pp_sphere *b) { if (a-inv_mass 0 b-inv_mass 0) { return; } pp_vec3 delta pp_vec3_sub(b-pos, a-pos); pp_fixed dist_sq pp_vec3_length_squared(delta); pp_fixed r_sum a-radius b-radius; pp_fixed r_sum_sq (pp_fixed)(((int64_t)r_sum * r_sum) PP_FIXED_SHIFT); if (dist_sq r_sum_sq) { return; } pp_fixed dist pp_fixed_sqrt(dist_sq); if (dist PP_F2X(0.001f)) { // 避免除以零兩個球完全重疊時人為指定一個方向 delta.x PP_FIXED_ONE; delta.y 0; delta.z 0; dist PP_FIXED_ONE; } // 法線方向 pp_vec3 normal; normal.x (pp_fixed)(((int64_t)delta.x * PP_FIXED_ONE) / dist); normal.y (pp_fixed)(((int64_t)delta.y * PP_FIXED_ONE) / dist); normal.z (pp_fixed)(((int64_t)delta.z * PP_FIXED_ONE) / dist); // 穿透深度 pp_fixed penetration r_sum - dist; // 位置修正按逆質(zhì)量比例分配 pp_fixed total_inv_mass a-inv_mass b-inv_mass; if (total_inv_mass 0) { return; } pp_fixed a_ratio a-inv_mass / total_inv_mass; pp_fixed b_ratio b-inv_mass / total_inv_mass; a-pos.x - (pp_fixed)(((int64_t)normal.x * penetration * a_ratio) / PP_FIXED_ONE / PP_FIXED_ONE); a-pos.y - (pp_fixed)(((int64_t)normal.y * penetration * a_ratio) / PP_FIXED_ONE / PP_FIXED_ONE); a-pos.z - (pp_fixed)(((int64_t)normal.z * penetration * a_ratio) / PP_FIXED_ONE / PP_FIXED_ONE); b-pos.x (pp_fixed)(((int64_t)normal.x * penetration * b_ratio) / PP_FIXED_ONE / PP_FIXED_ONE); b-pos.y (pp_fixed)(((int64_t)normal.y * penetration * b_ratio) / PP_FIXED_ONE / PP_FIXED_ONE); b-pos.z (pp_fixed)(((int64_t)normal.z * penetration * b_ratio) / PP_FIXED_ONE / PP_FIXED_ONE); }然后在pp_world_step中對所有球體兩兩調(diào)用這個碰撞函數(shù)for (int i 0; i world-sphere_count; i) { for (int j i 1; j world-sphere_count; j) { pp_collide_spheres(world-spheres[i], world-spheres[j]); } }6.2 使用 AABB 簡化靜態(tài)場景雖然球體碰撞很適合模擬粒子、拋射物但游戲場景中的墻壁、平臺、柱子基本都是方形的。為了簡化計算可以引入 AABB 結(jié)構(gòu)。typedef struct pp_aabb { pp_vec3 min; pp_vec3 max; } pp_aabb;球與 AABB 的碰撞檢測思路是先找到 AABB 上離球心最近的點然后計算該點與球心的距離。如果距離小于球半徑則發(fā)生碰撞。這種方法在尼采風(fēng)格的單文件實現(xiàn)中非常實用因為 AABB 的包圍盒計算和相交檢測都只是整數(shù)比較和加減法開銷極低。7. 常見問題與排查思路在編寫和移植單文件物理引擎時開發(fā)者經(jīng)常遇到下面幾類問題。這里整理了排查清單。7.1 物體陷入地面或穿模問題現(xiàn)象常見原因解決思路高速下落時穿過地面時間步長太大一幀內(nèi)位移超過碰撞體尺寸減少dt或使用連續(xù)碰撞檢測射線掃描物體在半空中卡住位置修正與速度修正順序錯誤先修正位置再修正速度速度響應(yīng)方向和法線對齊物體彈跳越彈越高恢復(fù)系數(shù)計算錯誤導(dǎo)致能量增加檢查碰撞響應(yīng)中乘法的定點數(shù)縮放是否多乘了PP_FIXED_ONE兩個球體互相嵌入無法分離碰撞檢測順序不穩(wěn)定穿透深度計算溢出使用 64 位中間變量做乘法并進(jìn)行多輪位置修正7.2 定點數(shù)運算結(jié)果偏差大定點數(shù)最大的坑是溢出和精度丟失。很多開發(fā)者從浮點移植到定點后發(fā)現(xiàn)物體運動軌跡和預(yù)期差別很大。常見原因包括乘法沒有先提升到 64 位?？s放因子選擇過大導(dǎo)致乘法結(jié)果溢出。除法直接使用整數(shù)除法小數(shù)部分被截斷。排查步驟把關(guān)鍵中間變量打印出來對比浮點版本的期望值。使用int64_t作為所有乘法的中間類型。在進(jìn)行多次連續(xù)運算時注意每次運算后都要右移PP_FIXED_SHIFT位恢復(fù)定點格式。對于除法可以用(a * PP_FIXED_ONE) / b的方式保留精度。7.3 移植到復(fù)古平臺后性能不符合預(yù)期PC 上運行良好的物理代碼移植到 PSX 或 N64 后可能變慢。主要原因通常不是算法復(fù)雜度而是編譯優(yōu)化級別沒有開。使用了太多 64 位乘法long long導(dǎo)致庫函數(shù)被大量調(diào)用。代碼里不小心引入了浮點類型導(dǎo)致編譯器鏈入軟件浮點庫。動態(tài)內(nèi)存分配malloc/realloc在視頻幀內(nèi)頻繁調(diào)用造成性能抖動。建議方案在目標(biāo)平臺工具鏈中開啟-O2甚至-O3。優(yōu)先使用 32 位乘法并在代碼層面避免溢出。如果必須使用 64 位盡量把乘累加操作合并。如果物理系統(tǒng)需要支持大量物體建議在初始化時一次性分配固定數(shù)組而不是每幀動態(tài)分配。8. 工程化最佳實踐從示例到產(chǎn)品級代碼8.1 數(shù)據(jù)結(jié)構(gòu)設(shè)計前面的示例使用realloc動態(tài)擴(kuò)展球體數(shù)組這種寫法在 PC 演示中沒問題但在復(fù)古平臺上并不適合。更推薦的做法是使用固定大小數(shù)組#define PP_MAX_BODIES 128 typedef struct pp_world { pp_sphere spheres[PP_MAX_BODIES]; int sphere_count; pp_fixed gravity; pp_fixed dt; pp_fixed floor_y; } pp_world;這種設(shè)計有幾個好處規(guī)避了內(nèi)存碎片問題。物理模塊不依賴堆分配?？梢造o態(tài)分配在數(shù)據(jù)段不占用運行時棧空間。性能可預(yù)期幀率更容易穩(wěn)定。8.2 使用參數(shù)表管理物理手感物理引擎中充滿了魔法數(shù)字重力加速度、恢復(fù)系數(shù)、摩擦系數(shù)、時間步長。建議在單文件內(nèi)部定義參數(shù)結(jié)構(gòu)體統(tǒng)一管理。typedef struct pp_params { pp_fixed gravity; pp_fixed dt; pp_fixed restitution; pp_fixed friction; pp_fixed floor_y; } pp_params; static pp_params g_params;業(yè)務(wù)代碼在初始化時設(shè)置一次參數(shù)表物理模塊內(nèi)部只讀取全局參數(shù)避免在游戲邏輯中散落各種數(shù)字。8.3 記錄與回放機制復(fù)古平臺游戲的物理調(diào)試非常困難因為無法在真機上打斷點。一個實用的方案是在 PC 上實現(xiàn)物理狀態(tài)的記錄與回放在物理更新前后把物體的位置、速度序列化到緩沖區(qū)。在真機運行出現(xiàn) Bug 時把記錄數(shù)據(jù)傳回 PC。PC 端還原 Bug 場景逐幀觀察物理狀態(tài)變化。這個功能雖然需要額外代碼但在復(fù)雜物理場景中能節(jié)省數(shù)倍排查時間。8.4 物理與渲染的分離在復(fù)古平臺上物理更新頻率和渲染幀率不一定相等。建議把物理更新獨立出來// 偽代碼示例 void game_loop() { while (run) { draw(); // 渲染當(dāng)前狀態(tài) for (int i 0; i substeps; i) { pp_world_step(world); // 物理子步進(jìn) } } }當(dāng)游戲以 60 FPS 渲染但希望物理更穩(wěn)定時可以把物理拆分到 120 Hz 或 240 Hz。子步進(jìn)可以降低高速物體穿模的風(fēng)險代價是每幀物理運算量增加需要在性能和穩(wěn)定性之間做平衡。9. 擴(kuò)展方向從“能用”到“實用”9.1 關(guān)節(jié)與鉸鏈如果需要搭建更復(fù)雜的機械結(jié)構(gòu)比如平臺游戲中的杠桿、旋轉(zhuǎn)門、鎖鏈可以考慮在物理核心中加入簡單的距離約束。距離約束的原理是每次迭代后修正兩個物體的位置讓它們之間的距離保持在一個目標(biāo)范圍內(nèi)。雖然這比單純碰撞復(fù)雜但在單文件架構(gòu)中仍然可控。9.2 射線檢測很多復(fù)古游戲需要在場景中做彈道判定、角色尋路、視線遮擋檢測。射線檢測是物理引擎中性價比極高的功能。核心思路是使用參數(shù)方程從起點逐步推進(jìn)依次檢測射線與各碰撞體的相交情況。由于復(fù)古平臺場景幾何比較簡單通常只需要把射線與 AABB、球體做相交測試。9.3 流體與粒子特效PSX 和 Dreamcast 上很多令人驚嘆的視覺效果其實是粒子系統(tǒng)。物理引擎可以為粒子系統(tǒng)提供重力、風(fēng)力、碰撞反彈等基本行為。單文件物理內(nèi)核天然適合這種需求粒子數(shù)量多、結(jié)構(gòu)簡單、不需要復(fù)雜的關(guān)節(jié)約束只需要求每幀更新速度、位置、生命值即可。10. 總結(jié)與下一步建議本文圍繞 Picophysics 和復(fù)古游戲平臺單文件物理引擎這個主題完整梳理了從硬件約束、定點數(shù)原理、碰撞響應(yīng)、可運行代碼到工程化建議的一系列內(nèi)容。你可以看到單文件并不是簡陋的代名詞而是一種在資源受限環(huán)境下非常高明的設(shè)計手段它讓物理系統(tǒng)保持極低的使用門檻同時具備足夠清晰的擴(kuò)展路徑。如果你希望進(jìn)一步探索這個方向可以做如下嘗試先仿照本文代碼寫一個球體自由落體演示觀察不同恢復(fù)系數(shù)下反彈的差異。把球體碰撞擴(kuò)展到 10 個、50 個、100 個記錄每幀 CPU 消耗。替換平方根算法比較不同實現(xiàn)的速度與精度。在模擬器上嘗試把物理輸出接入畫面渲染驗證真機性能。復(fù)古平臺開發(fā)的最大樂趣正在于用有限的硬件逼出自己的高效設(shè)計。物理引擎更像是一個練手的好場景它的輸入輸出非常清晰又要求你在數(shù)學(xué)、數(shù)據(jù)結(jié)構(gòu)和性能優(yōu)化上都投入思考。希望這份筆記能幫你少走一些彎路也期待你寫出在 N64 或 PSX 上運行流暢的物理效果。