您现在的位置：程式碼之家 >> 編程語言 > >> C語言 >> C++ >> C++ 入門

C++ 實時水波模擬：從物理特性到高效算法實現

編輯：C++ 入門

引言

水波模擬是計算機圖形學中的經典問題。要在實時應用中呈現逼真的水面效果，必須在物理真實感與計算效率之間取得平衡。本文將詳細介紹如何利用 C++ 實現一個高效的實時水波模擬系統，涵蓋波的物理特性、數學建模、數據結構設計以及渲染技術。

波的物理特性分析

在開始編程之前，我們需要理解水波的幾個核心物理特性：

擴散性：水波上的任意一點都會以自己為中心向四周擴散，環狀波紋的形成是因為內部擴散的對稱相互抵消。
衰減性：水具有阻尼特性，否則波濤將永不停息。
折射效應：水面傾斜導致光線折射，我們看到的水底景物位置會發生偏移，這正是我們感知波浪形狀的原因。
反射與繞射：波遇障礙物反射，通過縫隙或繞過障礙物時發生繞射。

核心算法：快速近似波動方程

傳統正弦函數或精確物理公式計算量過大，無法滿足實時渲染（≥20 FPS）的需求。我們採用一種基於差分的快速近似算法，僅使用加減法和位移運算。

雙緩衝區設計

建立兩個與水面圖像等大的數組 buf1 和 buf2，分別存儲前一時刻和後一時刻每個像素點的波幅（能量）數據。初始狀態下水面平靜，所有值為 0。

波幅遞推公式推導

假設某點下一時刻的振幅受其上下左右四鄰域點及自身當前振幅影響，且四鄰域影響機會均等：

X0' = a(X1+X2+X3+X4) + bX0

在無阻尼理想條件下，系統總能量守恆，即所有點振幅之和不變，可得 4a + b = 1。取最簡解 a = 1/2, b = -1，因 1/2 可用右移運算 >>1 實現，避免乘除。最終公式為：

X0' = (X1+X2+X3+X4)/2 - X0

引入實際阻尼，每次計算後波幅按比例衰減，經測試 1/32 (即 >>5) 為宜。

核心計算代碼

void RippleSpread() {
    for (int i = BACKWIDTH; i < BACKWIDTH*BACKHEIGHT - BACKWIDTH; i++) {
        // 波能擴散
        buf2[i] = ((buf1[i-1] + buf1[i+1] + buf1[i-BACKWIDTH] + buf1[i+BACKWIDTH]) >> 1) - buf2[i];
        // 波能衰減
        buf2[i] -= buf2[i] >> 5;
    }
    // 交換緩衝區
    short *ptmp = buf1;
    buf1 = buf2;
    buf2 = ptmp;
}

渲染：基於折射的像素偏移

利用水面折射原理，水面傾斜導致視線偏移。我們用相鄰點波幅之差近似偏移量：

xoff = buf1[k-1] - buf1[k+1];  // 左右波幅差
yoff = buf1[k-BACKWIDTH] - buf1[k+BACKWIDTH];  // 上下波幅差

根據偏移量從原始圖像采樣像素到渲染目標，實現水面折射視覺效果。渲染循環中需進行邊界檢查，防止越界訪問。

波源激發：投石入水

在指定位置製造負脈衝波源，模擬投石入水。以石頭入水點為圓心，在半徑範圍內將波幅設置為負值（-n, n∈[32,128]），圓形區域近似處理。代碼如下：

void DropStone(int x, int y, int stonesize, int stoneweight) {
    // 邊界檢查
    for (int posx = x-stonesize; posx < x+stonesize; posx++)
        for (int posy = y-stonesize; posy < y+stonesize; posy++)
            if ((posx-x)*(posx-x) + (posy-y)*(posy-y) < stonesize*stonesize)
                buf1[BACKWIDTH*posy + posx] = -stoneweight;
}