用球座標計算
注意到被積函數與區域都具有徑向對稱性。 被積函數 與邊界 都只依賴到原點的距離。這就是該用球座標的明確訊號:在直角座標下這個積分需要層層套疊的根號上下限,換成球座標後上下限全是常數。
設定球座標代換,別漏掉 Jacobian。 取
被積函數縮成 ,而體積元素帶著 Jacobian 因子
不是可有可無的形式細節 —— 它是座標變換在局部的體積伸縮因子,漏掉它是這個積分最常見的出錯方式。
寫出上下限。 單位球可描述為
注意 (極角,自 軸往下量)只跑到 ,而 (方位角)跑完整的 。把這兩個範圍對調會把球重複算一次。
組出積分並拆開。 被積函數乘上 Jacobian 得 ,所以
因為被積函數可分解成( 的函數)( 的函數)(對 為常數),而且上下限都是常數,這個三重積分就拆成三個單變數積分的乘積。
分別算出三個因子。
相乘:
檢查答案的量級。 球的體積是 ,被積函數 從球心的 變到球面的 ,所以積分必須嚴格落在 與 之間。 在球上的平均值是 ,合理地大於 ,因為球的體積大部分集中在靠外的殼層。在外接正方體上做蒙地卡羅估計得到 ,印證了 。
Need to solve a different problem like this? Open the solver →