好的,我们来聊聊如何在 R^n 上计算这个积分:$int_{mathbb{R}^n} |x|^2 e^{|x|^2} dx$。这个积分在很多领域,比如统计物理、量子力学以及概率论中都非常常见。
首先,让我们先分解一下被积函数 $|x|^2 e^{|x|^2}$。
$|x|^2$: 在 R^n 中,如果 $x = (x_1, x_2, dots, x_n)$,那么 $|x|^2 = x_1^2 + x_2^2 + dots + x_n^2$。
$e^{|x|^2}$: 这是一个高斯函数的形式,它在原点处达到最大值,并随着 $|x|$ 的增大而迅速衰减。
这个积分的计算,最直接也最有效的方法是利用 高斯积分 的知识以及 球坐标系 的思想。
核心思路:将多维积分转化为一维积分
R^n 上的积分直接计算会很困难,特别是当被积函数涉及到 $|x|$(距离原点的距离)时。这时候,我们往往会尝试使用高维的球坐标系来简化计算。
第一步:引入一个辅助的、独立的 R^n 高斯积分
考虑一个非常著名的积分:
$$I_n = int_{mathbb{R}^n} e^{|x|^2} dx$$
这个积分是高斯函数在 R^n 上的积分。我们可以通过将它写成 n 个独立的一维高斯积分的乘积来计算:
$$I_n = int_{mathbb{R}^n} e^{(x_1^2 + x_2^2 + dots + x_n^2)} dx_1 dx_2 dots dx_n$$
$$I_n = left(int_{infty}^{infty} e^{x_1^2} dx_1
ight) left(int_{infty}^{infty} e^{x_2^2} dx_2
ight) dots left(int_{infty}^{infty} e^{x_n^2} dx_n
ight)$$
我们知道一维高斯积分 $int_{infty}^{infty} e^{x^2} dx = sqrt{pi}$。
所以,
$$I_n = (sqrt{pi})^n = pi^{n/2}$$
这个结果非常重要,我们将会在后面的计算中用到它。
第二步:对被积函数进行操作,引入 $|x|^2$
我们的目标积分是 $J_n = int_{mathbb{R}^n} |x|^2 e^{|x|^2} dx$。
如何引入 $|x|^2$ 呢?我们可以利用导数性质。考虑对 $I_n$ 做一些“微调”。
一种巧妙的方法是考虑这样一个量:
$$int_{mathbb{R}^n} e^{a|x|^2} dx$$
其中 $a$ 是一个正的常数。通过类似的推导(将 $x_i^2$ 替换为 $ax_i^2$),我们得到:
$$int_{mathbb{R}^n} e^{a|x|^2} dx = left(int_{infty}^{infty} e^{ax^2} dx
ight)^n$$
对于一维积分 $int_{infty}^{infty} e^{ax^2} dx$,我们可以做变量替换 $y = sqrt{a}x$,则 $dy = sqrt{a} dx$,$dx = frac{1}{sqrt{a}} dy$。
$$int_{infty}^{infty} e^{ax^2} dx = int_{infty}^{infty} e^{y^2} frac{1}{sqrt{a}} dy = frac{1}{sqrt{a}} int_{infty}^{infty} e^{y^2} dy = frac{sqrt{pi}}{sqrt{a}}$$
所以,
$$int_{mathbb{R}^n} e^{a|x|^2} dx = left(frac{sqrt{pi}}{sqrt{a}}
ight)^n = frac{pi^{n/2}}{a^{n/2}}$$
现在,我们来看如何利用导数来引入 $|x|^2$。
考虑对关于参数 $a$ 的函数 $f(a) = int_{mathbb{R}^n} e^{a|x|^2} dx$ 进行求导。
$$f'(a) = frac{d}{da} int_{mathbb{R}^n} e^{a|x|^2} dx$$
我们通常可以在积分号下求导(假设收敛性满足条件,这里是满足的):
$$f'(a) = int_{mathbb{R}^n} frac{partial}{partial a} (e^{a|x|^2}) dx$$
$$f'(a) = int_{mathbb{R}^n} (|x|^2) e^{a|x|^2} dx$$
$$f'(a) = int_{mathbb{R}^n} |x|^2 e^{a|x|^2} dx$$
另一方面,我们已经知道 $f(a) = frac{pi^{n/2}}{a^{n/2}} = pi^{n/2} a^{n/2}$。
对 $f(a)$ 关于 $a$ 求导:
$$f'(a) = pi^{n/2} left(frac{n}{2}
ight) a^{n/2 1}$$
$$f'(a) = frac{n}{2} pi^{n/2} a^{(n+2)/2}$$
现在,我们令 $a=1$(因为我们的目标积分是 $e^{|x|^2}$,相当于 $a=1$ 的情况)。
从 $f'(a) = int_{mathbb{R}^n} |x|^2 e^{a|x|^2} dx$ 来看,当 $a=1$ 时:
$$f'(1) = int_{mathbb{R}^n} |x|^2 e^{|x|^2} dx$$
从 $f'(a) = frac{n}{2} pi^{n/2} a^{(n+2)/2}$ 来看,当 $a=1$ 时:
$$f'(1) = frac{n}{2} pi^{n/2} (1)^{(n+2)/2} = frac{n}{2} pi^{n/2}$$
将两者联立,我们就得到了:
$$int_{mathbb{R}^n} |x|^2 e^{|x|^2} dx = frac{n}{2} pi^{n/2}$$
因此,
$$int_{mathbb{R}^n} |x|^2 e^{|x|^2} dx = frac{n}{2} pi^{n/2}$$
这个方法非常简洁高效。
第二种方法:利用球坐标系和期望值
这种方法可能更加直观一些,因为它直接利用了 R^n 的球对称性。
第一步:利用 R^n 的球坐标系
R^n 的球坐标系可以将一个点 $x$ 表示为 $(r, Omega)$,其中 $r = |x|$ 是径向距离,而 $Omega$ 是位于单位球面 $S^{n1}$ 上的角度坐标。
在 R^n 中,体积元 $dx$ 可以写成:
$$dx = r^{n1} dr dOmega$$
其中 $dOmega$ 是单位球面 $S^{n1}$ 上的面积元。
被积函数 $|x|^2 e^{|x|^2}$ 在球坐标系下变成 $r^2 e^{r^2}$。
所以,积分变为:
$$int_{mathbb{R}^n} |x|^2 e^{|x|^2} dx = int_0^infty int_{S^{n1}} r^2 e^{r^2} r^{n1} dr dOmega$$
$$= int_0^infty r^{n+1} e^{r^2} dr int_{S^{n1}} dOmega$$
第二步:计算单位球面的面积
$int_{S^{n1}} dOmega$ 是单位球面 $S^{n1}$ 的总面积。单位球面的面积公式是:
$$A_{n1} = frac{2pi^{n/2}}{Gamma(n/2)}$$
其中 $Gamma(z)$ 是伽马函数。我们知道 $Gamma(1/2) = sqrt{pi}$ 和 $Gamma(1) = 1$。
当 n 为偶数,比如 $n=2k$,则 $n/2 = k$,$Gamma(k) = (k1)!$。所以 $A_{2k1} = frac{2pi^k}{k!} imes (2k1)!$ (这里的记号有点乱,更准确的是球面 $S^{n1}$ 的面积)
更常见的单位球面 $S^{n1}$ 的面积是:
$$A_{n1} = frac{2pi^{n/2}}{Gamma(n/2)}$$
所以积分是:
$$J_n = left(int_0^infty r^{n+1} e^{r^2} dr
ight) left(frac{2pi^{n/2}}{Gamma(n/2)}
ight)$$
第三步:计算径向积分
我们需要计算 $int_0^infty r^{n+1} e^{r^2} dr$。
我们可以做变量替换 $u = r^2$。那么 $du = 2r dr$。
$r = sqrt{u}$,$r^{n+1} = (sqrt{u})^{n+1} = u^{(n+1)/2}$。
$dr = frac{1}{2sqrt{u}} du = frac{1}{2} u^{1/2} du$。
代入积分:
$$int_0^infty r^{n+1} e^{r^2} dr = int_0^infty u^{(n+1)/2} e^{u} left(frac{1}{2} u^{1/2} du
ight)$$
$$= frac{1}{2} int_0^infty u^{(n+1)/2 1/2} e^{u} du$$
$$= frac{1}{2} int_0^infty u^{n/2} e^{u} du$$
这个积分的形式正是伽马函数的定义:$Gamma(z) = int_0^infty t^{z1} e^{t} dt$。
所以,我们的积分是 $frac{1}{2} Gamma(n/2 + 1)$。
第四步:将结果组合起来
现在我们将径向积分和球面面积结合起来:
$$J_n = left(frac{1}{2} Gamma(n/2 + 1)
ight) left(frac{2pi^{n/2}}{Gamma(n/2)}
ight)$$
我们知道伽马函数的一个重要性质是 $Gamma(z+1) = z Gamma(z)$。
所以 $Gamma(n/2 + 1) = frac{n}{2} Gamma(n/2)$。
代入上式:
$$J_n = left(frac{1}{2} cdot frac{n}{2} Gamma(n/2)
ight) left(frac{2pi^{n/2}}{Gamma(n/2)}
ight)$$
$$J_n = frac{n}{4} Gamma(n/2) cdot frac{2pi^{n/2}}{Gamma(n/2)}$$
消去 $Gamma(n/2)$:
$$J_n = frac{n}{4} cdot 2pi^{n/2}$$
$$J_n = frac{n}{2} pi^{n/2}$$
两种方法都得到了相同的结果,$frac{n}{2} pi^{n/2}$。第一种方法利用了对参数求导的技巧,非常巧妙;第二种方法则更加依赖于对 R^n 球坐标系和伽马函数的熟悉程度。
举个例子:n=3 的情况
让我们来验证一下 n=3 的结果。
目标是计算 $int_{mathbb{R}^3} |x|^2 e^{|x|^2} dx$。
根据公式,结果应该是 $frac{3}{2} pi^{3/2}$。
我们来用方法一验证一下:
考虑 $I_3(a) = int_{mathbb{R}^3} e^{a|x|^2} dx = frac{pi^{3/2}}{a^{3/2}}$。
$I_3'(a) = frac{d}{da} left(pi^{3/2} a^{3/2}
ight) = pi^{3/2} left(frac{3}{2}
ight) a^{5/2}$。
同时,$I_3'(a) = int_{mathbb{R}^3} (|x|^2) e^{a|x|^2} dx = int_{mathbb{R}^3} |x|^2 e^{a|x|^2} dx$。
令 $a=1$:
$int_{mathbb{R}^3} |x|^2 e^{|x|^2} dx = pi^{3/2} left(frac{3}{2}
ight) (1)^{5/2} = frac{3}{2} pi^{3/2}$。
所以 $int_{mathbb{R}^3} |x|^2 e^{|x|^2} dx = frac{3}{2} pi^{3/2}$。结果一致。
用方法二验证一下:
$int_{mathbb{R}^3} |x|^2 e^{|x|^2} dx = int_0^infty r^2 e^{r^2} r^{31} dr int_{S^2} dOmega$
$S^2$ 的面积是 $4pi$。
径向积分是 $int_0^infty r^2 e^{r^2} r^2 dr = int_0^infty r^4 e^{r^2} dr$。
令 $u = r^2$, $du = 2r dr$, $dr = frac{1}{2sqrt{u}} du$。
$int_0^infty u^2 e^{u} frac{1}{2} u^{1/2} du = frac{1}{2} int_0^infty u^{3/2} e^{u} du = frac{1}{2} Gamma(5/2)$。
$Gamma(5/2) = frac{3}{2} Gamma(3/2) = frac{3}{2} cdot frac{1}{2} Gamma(1/2) = frac{3}{4} sqrt{pi}$。
所以径向积分是 $frac{1}{2} cdot frac{3}{4} sqrt{pi} = frac{3}{8} sqrt{pi}$。
总积分是 $frac{3}{8} sqrt{pi} cdot 4pi = frac{3}{2} pi^{3/2}$。结果也一致。
总结一下计算步骤:
1. 识别积分形式: 积分包含 $|x|^2$ 和指数项 $e^{|x|^2}$,这提示使用高斯积分和球坐标系。
2. 方法一 (参数求导):
考虑更一般的积分 $int_{mathbb{R}^n} e^{a|x|^2} dx = pi^{n/2} a^{n/2}$。
对关于 $a$ 求导,得到 $int_{mathbb{R}^n} |x|^2 e^{a|x|^2} dx = frac{n}{2} pi^{n/2} a^{(n+2)/2}$。
令 $a=1$,直接得到结果 $int_{mathbb{R}^n} |x|^2 e^{|x|^2} dx = frac{n}{2} pi^{n/2}$。
3. 方法二 (球坐标与伽马函数):
将积分转化为球坐标下的形式:$int_0^infty r^{n+1} e^{r^2} dr int_{S^{n1}} dOmega$。
计算单位球面的面积 $A_{n1} = frac{2pi^{n/2}}{Gamma(n/2)}$。
计算径向积分 $int_0^infty r^{n+1} e^{r^2} dr$,通过变量替换得到 $frac{1}{2} Gamma(n/2 + 1)$。
利用 $Gamma(z+1) = zGamma(z)$ 将结果组合,得到 $frac{n}{2} pi^{n/2}$。
这两种方法都非常强大,而且在处理这类积分时经常会用到。希望这个详细的解释能帮助你理解这个计算过程!