卡在這題?用 AI 解開它
Calculus · real student question

計算在整個 R³ 上 exp(−α(x² + y² + z²)) 乘上下列和的三重積分:a cos(xξ + yη + zζ) 加 b sin(xξ + yη + zζ) 加 c(x² + y² + z²)。

Question

設 α>0\alpha > 0 且 (ξ,η,ζ)∈R3(\xi, \eta, \zeta) \in \mathbb{R}^3,計算

I=∭R3e−α(x2+y2+z2)(acos⁡(xξ+yη+zζ)+bsin⁡(xξ+yη+zζ)+c (x2+y2+z2))dx dy dzI = \iiint_{\mathbb{R}^3} e^{-\alpha(x^2+y^2+z^2)}\Bigl(a\cos(x\xi + y\eta + z\zeta) + b\sin(x\xi + y\eta + z\zeta) + c\,(x^2+y^2+z^2)\Bigr) dx\,dy\,dz

Step-by-step solution

  1. 改用向量記號,並拆成三個積分。 記 x=(x,y,z)\mathbf{x} = (x,y,z)、k=(ξ,η,ζ)\mathbf{k} = (\xi,\eta,\zeta)、r=∥x∥r = \|\mathbf{x}\|、k=∥k∥k = \|\mathbf{k}\|。則 xξ+yη+zζ=x⋅kx\xi + y\eta + z\zeta = \mathbf{x}\cdot\mathbf{k},而

    I=aI1+bI2+cI3I = aI_1 + bI_2 + cI_3

    其中 I1=∭e−αr2cos⁡(x⋅k) d3xI_1 = \iiint e^{-\alpha r^2}\cos(\mathbf{x}\cdot\mathbf{k})\,d^3x,I2I_2 是把餘弦換成正弦,I3=∭e−αr2r2 d3xI_3 = \iiint e^{-\alpha r^2} r^2\,d^3x。條件 α>0\alpha > 0 正是讓這三個積分都收斂的關鍵。

  2. 用奇偶性幹掉正弦項。 在 x↦−x\mathbf{x} \mapsto -\mathbf{x} 之下,高斯因子 e−αr2e^{-\alpha r^2} 不變,而 sin⁡(x⋅k)\sin(\mathbf{x}\cdot\mathbf{k}) 變號,所以被積函數在對原點對稱的區域上是奇的:

    I2=0I_2 = 0

    因此係數 bb 根本不會出現在答案裡。

  3. 把餘弦項當成傅立葉轉換來算。 標準的三維高斯轉換是

    ∭R3e−αr2eix⋅k d3x=(πα)3/2e−k2/(4α)\iiint_{\mathbb{R}^3} e^{-\alpha r^2} e^{i\mathbf{x}\cdot\mathbf{k}}\,d^3x = \left(\frac{\pi}{\alpha}\right)^{3/2} e^{-k^2/(4\alpha)}

    (配方之後,每一個直角座標方向各貢獻 π/α e−ki2/4α\sqrt{\pi/\alpha}\,e^{-k_i^2/4\alpha})。取實部——虛部恰好就是 I2=0I_2 = 0 ——得到

    I1=(πα)3/2e−k2/(4α)I_1 = \left(\frac{\pi}{\alpha}\right)^{3/2} e^{-k^2/(4\alpha)}

  4. 用球座標算 r² 那一項。 被積函數只依賴 rr,所以 d3x=4πr2 drd^3x = 4\pi r^2\,dr:

    I3=4π∫0∞r4e−αr2 drI_3 = 4\pi\int_0^{\infty} r^4 e^{-\alpha r^2}\,dr

    代換 u=αr2u = \alpha r^2 把它化成 Gamma 積分:

    ∫0∞r4e−αr2dr=Γ(5/2)2α5/2=3π8α5/2\int_0^\infty r^4 e^{-\alpha r^2}dr = \frac{\Gamma(5/2)}{2\alpha^{5/2}} = \frac{3\sqrt{\pi}}{8\alpha^{5/2}}

    這裡用到 Γ(5/2)=34π\Gamma(5/2) = \tfrac{3}{4}\sqrt{\pi}。因此 I3=3π3/22α5/2I_3 = \dfrac{3\pi^{3/2}}{2\alpha^{5/2}}。

  5. 組出結果。

    I=a(πα)3/2e−(ξ2+η2+ζ2)/(4α)+3c π3/22α5/2I = a\left(\frac{\pi}{\alpha}\right)^{3/2} e^{-(\xi^2+\eta^2+\zeta^2)/(4\alpha)} + \frac{3c\,\pi^{3/2}}{2\alpha^{5/2}}

    注意這個結構:含 aa 的部分隨 kk 衰減(在頻率空間裡也是個高斯),而含 cc 的部分則是與 k\mathbf{k} 無關的純常數。

  6. 數值驗證。 取 α=0.7\alpha = 0.7、k=(0.3,−0.5,1.1)\mathbf{k} = (0.3, -0.5, 1.1)、a=1.3a = 1.3、b=2.1b = 2.1、c=0.9c = 0.9,封閉形式給出 25.442094725.4420947。在 [−6,6]3[-6,6]^3 上用 9090 個節點的高斯–勒壤得網格(高斯因子超出這個範圍後已可忽略)同樣得到 25.442094725.4420947 ——吻合到八位數字,也印證了 b=2.1b = 2.1 那一項毫無貢獻。

Answer

a(πα)3/2e−(ξ2+η2+ζ2)/(4α)+3c π3/22α5/2a\left(\frac{\pi}{\alpha}\right)^{3/2} e^{-(\xi^2+\eta^2+\zeta^2)/(4\alpha)} + \frac{3c\,\pi^{3/2}}{2\alpha^{5/2}}

Need to solve a different problem like this? Open the solver →