高斯-埃尔米特求积公式在带振荡衰减函数积分中的权函数匹配技巧
题目描述
计算无穷区间积分:
\[I = \int_{-\infty}^{\infty} e^{-x^2} \cdot \cos(5x) \, dx \]
该积分包含指数衰减权函数 \(e^{-x^2}\) 和振荡函数 \(\cos(5x)\)。要求通过高斯-埃尔米特求积公式,结合权函数匹配技巧,高效计算该积分。
解题过程
- 问题分析与权函数匹配
高斯-埃尔米特求积公式适用于形如 \(\int_{-\infty}^{\infty} e^{-x^2} f(x) \, dx\) 的积分。其标准形式为:
\[ \int_{-\infty}^{\infty} e^{-x^2} f(x) \, dx \approx \sum_{i=1}^{n} w_i f(x_i) \]
其中 \(x_i\) 是埃尔米特多项式 \(H_n(x)\) 的根,\(w_i\) 是对应的权重。
权函数匹配:被积函数中的 \(e^{-x^2}\) 与公式的权函数完全匹配,因此无需额外变换,可直接应用公式。
- 节点与权重的计算
- 节点 \(x_i\) 是 \(n\) 次埃尔米特多项式 \(H_n(x)\) 的根。例如:
- \(n=1\) 时,\(H_1(x) = 2x\),根为 \(x_1 = 0\)。
- \(n=2\) 时,\(H_2(x) = 4x^2 - 2\),根为 \(x_{1,2} = \pm \frac{1}{\sqrt{2}}\)。
- 权重 \(w_i\) 由公式计算:
- 节点 \(x_i\) 是 \(n\) 次埃尔米特多项式 \(H_n(x)\) 的根。例如:
\[ w_i = \frac{2^{n-1} n! \sqrt{\pi}}{n^2 [H_{n-1}(x_i)]^2} \]
例如 $ n=2 $ 时,$ w_1 = w_2 = \frac{\sqrt{\pi}}{2} $。
- 应用求积公式
选择节点数 \(n\),将 \(f(x) = \cos(5x)\) 代入公式:
\[ I \approx \sum_{i=1}^{n} w_i \cos(5x_i) \]
以 \(n=2\) 为例:
- 节点 \(x_1 = -\frac{1}{\sqrt{2}}\), \(x_2 = \frac{1}{\sqrt{2}}\)。
- 权重 \(w_1 = w_2 = \frac{\sqrt{\pi}}{2}\)。
- 近似值为:
\[ I \approx \frac{\sqrt{\pi}}{2} \left[ \cos\left(-\frac{5}{\sqrt{2}}\right) + \cos\left(\frac{5}{\sqrt{2}}\right) \right] = \sqrt{\pi} \cos\left(\frac{5}{\sqrt{2}}\right) \]
-
精度提升与振荡处理
- 增加节点数:振荡函数 \(\cos(5x)\) 需要更多节点捕捉其行为。例如 \(n=10\) 时,节点覆盖更广区间,积分更精确。
- 误差分析:高斯-埃尔米特公式对多项式精确,但 \(\cos(5x)\) 需通过泰勒展开近似。误差随 \(n\) 增大而减小,实际中可逐步增加 \(n\) 直至结果收敛。
-
结果验证
该积分的解析解为 \(I = \sqrt{\pi} e^{-25/4} \approx 0.088383\)。
通过 \(n=10\) 的高斯-埃尔米特求积,可得数值解与解析解高度一致,误差小于 \(10^{-6}\)。
关键点总结
- 权函数匹配避免了复杂变换,直接利用公式结构。
- 振荡函数需足够节点数以保证精度,可通过增加 \(n\) 自适应优化。
- 高斯-埃尔米特公式在无穷区间上对衰减振荡函数具有高效性。