看完不会?让 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 就把它变成一个伽马积分:

    ∫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 →