Skip to content

混合大型运算的计算方式

摘要

在高等数学和数学分析中,我们经常面对同时涉及求和(Summation)求积(Product)积分(Integration)级数(Series) 的混合运算。这类运算的难点在于:直接硬算往往不可行,必须运用次序交换、级数展开、渐进估计等高级技巧来化简。本文系统梳理混合大型运算的核心策略与经典方法,每节配套详细例题与完整推导,帮助读者建立从识别问题类型到选择合适工具的完整思路链。


一、混合运算的基本策略

1.1 混合运算的分类

在大型运算中,常见的混合类型主要有三种:

类型形式典型场景
求和 + 求积$\sum \prod$或$\prod \sum$多项式展开、韦达定理应用
求和 + 积分$\sum \int$或$\int \sum$级数和、含参积分
求积 + 积分$\prod \int$或$\int \prod$无穷乘积、概率论中的矩

1.2 四大基本策略

处理混合运算有四条核心思路,层层递进:

策略一:分离变量 (Separation of Variables) 将混合运算拆分为独立的部分分别处理。例如 $\sum_{k=1}^n \int_a^b f_k(x) g(x)\thinspace{}dx$中,若$g(x)$与$k$无关,则$g(x)$ 可提到求和号外。

策略二:交换次序 (Interchange of Order) 在满足一定条件下(一致收敛、绝对可积等),交换 $\sum$与$\int$、$\sum$与$\sum$、$\int$与$\int$ 的次序,使运算变得可行。

策略三:递推降阶 (Recursive Reduction) 利用递推关系将高维/高阶问题降解为低维/低阶问题,逐层求解。

策略四:渐进估计 (Asymptotic Estimation) 当精确解不可求时,用渐进方法(大O、拉普拉斯方法、鞍点法)获得主导项。

1.3 运算优先级与符号规范

运算优先级

混合运算中,先内后外是基本原则:

  1. 最内层的积分/求和/求积运算符先计算
  2. 若运算符之间可以交换次序,优先选择使计算简单的顺序
  3. 先化简再计算:能用恒等式化简的部分先化简

符号规范:

  • 求和:$\displaystyle\sum_{k=m}^{n} a_k$,无穷和记 $\displaystyle\sum_{k=1}^{\infty} a_k$
  • 求积:$\displaystyle\prod_{k=1}^{n} a_k$
  • 积分:$\displaystyle\int_a^b f(x)\thinspace{}dx$
  • 多重运算:用括号明确优先级,如 $\displaystyle\sum_{k=1}^{n}\left(\int_0^1 f_k(x)\thinspace{}dx\right)$

二、求和与积分次序交换

2.1 有限和与积分的交换

对于有限和,积分与求和总是可以直接交换:

$$ \boxed{\int_a^b \sum_{k=1}^{n} f_k(x)\thinspace{}dx = \sum_{k=1}^{n} \int_a^b f_k(x)\thinspace{}dx} $$

这是因为有限和的积分等于积分的有限和,无需任何附加条件。

例题 2.1:计算 $\displaystyle\int_0^1 \sum_{k=1}^{n} k x^{k-1}\thinspace{}dx$

方法一:先积分后求和

$$ \int_0^1 \sum_{k=1}^{n} k x^{k-1}\thinspace{}dx = \int_0^1 \left(1 + 2x + 3x^2 + \cdots + n x^{n-1}\right) dx $$

注意到 $\sum_{k=1}^{n} k x^{k-1} = \frac{d}{dx}\left(\sum_{k=0}^{n} x^k\right) = \frac{d}{dx}\left(\frac{1-x^{n+1}}{1-x}\right)$。

计算导数: $$ \frac{d}{dx}\left(\frac{1-x^{n+1}}{1-x}\right) = \frac{-(n+1)x^n(1-x) + (1-x^{n+1})}{(1-x)^2} $$

直接积分比较复杂,换方法二。

方法二:先求和后积分(交换次序)

$$ \begin{aligned} \int_0^1 \sum_{k=1}^{n} k x^{k-1}\thinspace{}dx &= \sum_{k=1}^{n} k \int_0^1 x^{k-1}\thinspace{}dx \newline &= \sum_{k=1}^{n} k \cdot \left[\frac{x^k}{k}\right]_0^1 \newline &= \sum_{k=1}^{n} k \cdot \frac{1}{k} = \sum_{k=1}^{n} 1 = n \end{aligned} $$

两种方法结果一致,但交换次序后计算极其简单。这说明了选择正确的运算顺序的重要性。

2.2 无穷级数与积分的交换(逐项积分定理)

对于无穷级数,交换积分与求和需要一致收敛(Uniform Convergence)的条件。

逐项积分定理 (Term-by-Term Integration Theorem)

若 $\lbrace f_n(x)\rbrace $在$[a,b]$上连续,且$\sum f_n(x)$在$[a,b]$ 上一致收敛,则: $$ \int_a^b \sum_{n=1}^{\infty} f_n(x)\thinspace{}dx = \sum_{n=1}^{\infty} \int_a^b f_n(x)\thinspace{}dx $$

幂级数的特殊性质:幂级数 $\sum a_n x^n$ 在其收敛区间内内闭一致收敛,因此可以在收敛区间内逐项积分和逐项求导。

例题 2.2:利用 $\displaystyle\sum_{n=0}^{\infty} x^n = \frac{1}{1-x}$逐项积分求$\ln(1+x)$ 的级数展开

推导

已知几何级数:$\displaystyle\sum_{n=0}^{\infty} t^n = \frac{1}{1-t}$,$|t| < 1$。

将 $t$替换为$-t$: $$ \sum_{n=0}^{\infty} (-1)^n t^n = \frac{1}{1+t}, \quad |t| < 1 $$

两边从 $0$到$x$ 积分($|x| < 1$):

$$ \begin{aligned} \int_0^x \frac{1}{1+t}\thinspace{}dt &= \int_0^x \sum_{n=0}^{\infty} (-1)^n t^n\thinspace{}dt \newline &= \sum_{n=0}^{\infty} (-1)^n \int_0^x t^n\thinspace{}dt \quad \text{(逐项积分)} \newline &= \sum_{n=0}^{\infty} (-1)^n \frac{x^{n+1}}{n+1} \newline &= x - \frac{x^2}{2} + \frac{x^3}{3} - \frac{x^4}{4} + \cdots \end{aligned} $$

而左边 $\int_0^x \frac{1}{1+t}\thinspace{}dt = \ln(1+x)$,故:

$$ \boxed{\ln(1+x) = \sum_{n=1}^{\infty} \frac{(-1)^{n-1} x^n}{n}, \quad -1 < x \leq 1} $$

技巧总结

先找一个已知的简单级数(如几何级数),通过代换变量逐项积分或求导 → 得到新函数的级数展开。这是处理级数求和问题最常用的套路之一。

2.3 积分号下求和(Fubini 定理的离散版本)

对于非负函数,求和与积分的交换更加自由。

Tonelli 定理(离散版本)

若 $f_n(x) \geq 0$对所有$n$和$x$ 成立,则求和与积分可无条件交换: $$ \sum_{n=1}^{\infty} \int f_n(x)\thinspace{}dx = \int \sum_{n=1}^{\infty} f_n(x)\thinspace{}dx $$ 两边同时有限或同时无穷。

例题 2.3:计算 $\displaystyle\sum_{n=1}^{\infty} \int_0^1 x^n \ln x\thinspace{}dx$

:先计算内层积分。

$$ \int_0^1 x^n \ln x\thinspace{}dx $$

使用分部积分法(Integration by Parts):令 $u = \ln x$,$dv = x^n dx$,则 $du = \frac{1}{x}dx$,$v = \frac{x^{n+1}}{n+1}$。

$$ \begin{aligned} \int_0^1 x^n \ln x\thinspace{}dx &= \left[\frac{x^{n+1}}{n+1} \ln x\right]_0^1 - \int_0^1 \frac{x^{n+1}}{n+1} \cdot \frac{1}{x}\thinspace{}dx \newline &= 0 - \int_0^1 \frac{x^n}{n+1}\thinspace{}dx \newline &= -\frac{1}{(n+1)^2} \end{aligned} $$

注意:$\lim_{x \to 0^+} x^{n+1} \ln x = 0$(因为幂函数压倒对数函数)。

现在对 $n$ 求和:

$$ \sum_{n=1}^{\infty} \int_0^1 x^n \ln x\thinspace{}dx = -\sum_{n=1}^{\infty} \frac{1}{(n+1)^2} = -\left(\frac{1}{2^2} + \frac{1}{3^2} + \cdots\right) $$

已知 $\sum_{n=1}^{\infty} \frac{1}{n^2} = \frac{\pi^2}{6}$(巴塞尔问题),故:

$$ \sum_{n=1}^{\infty} \frac{1}{(n+1)^2} = \frac{\pi^2}{6} - 1 $$

因此原式 $= -\left(\frac{\pi^2}{6} - 1\right) = 1 - \frac{\pi^2}{6}$。

验证(交换次序)

也可先求和再积分:$\int_0^1 \ln x \sum_{n=1}^{\infty} x^n\thinspace{}dx = \int_0^1 \ln x \cdot \frac{x}{1-x}\thinspace{}dx$,结果相同。


三、求和与求积的相互转换

3.1 和转积:对数转换法

求和与求积之间的桥梁是对数函数:

$$ \boxed{\prod_{k=1}^{n} a_k = \exp\left(\sum_{k=1}^{n} \ln a_k\right)} $$

这一转换将乘积问题转化为求和问题,从而可以运用丰富的求和技术。

例题 3.1:计算 $\displaystyle\lim_{n \to \infty} \prod_{k=1}^{n} \left(1 + \frac{k}{n^2}\right)$

:取对数转和:

$$ \ln P_n = \ln \prod_{k=1}^{n} \left(1 + \frac{k}{n^2}\right) = \sum_{k=1}^{n} \ln\left(1 + \frac{k}{n^2}\right) $$

当 $n$ 很大时,$\frac{k}{n^2}$很小,使用近似$\ln(1+t) \approx t - \frac{t^2}{2}$(当 $t \to 0$):

$$ \begin{aligned} \ln P_n &= \sum_{k=1}^{n} \left[\frac{k}{n^2} - \frac{k^2}{2n^4} + O\left(\frac{k^3}{n^6}\right)\right] \newline &= \frac{1}{n^2} \sum_{k=1}^{n} k - \frac{1}{2n^4} \sum_{k=1}^{n} k^2 + \cdots \newline &= \frac{1}{n^2} \cdot \frac{n(n+1)}{2} - \frac{1}{2n^4} \cdot \frac{n(n+1)(2n+1)}{6} + \cdots \end{aligned} $$

简化主导项:

$$ \begin{aligned} \ln P_n &= \frac{n+1}{2n} - \frac{(n+1)(2n+1)}{12n^3} + \cdots \newline &= \frac{1}{2} + \frac{1}{2n} + O\left(\frac{1}{n^2}\right) \end{aligned} $$

因此:

$$ \lim_{n \to \infty} \ln P_n = \frac{1}{2} $$

$$ \boxed{\lim_{n \to \infty} \prod_{k=1}^{n} \left(1 + \frac{k}{n^2}\right) = e^{1/2} = \sqrt{e}} $$

3.2 积转和:部分分式与韦达定理

多项式 $\prod_{k=1}^n (x - a_k)$ 展开为幂和形式时,系数与根的对称多项式有关。

韦达定理 (Vieta's Formulas)

若 $P(x) = \prod_{k=1}^n (x - a_k) = x^n - e_1 x^{n-1} + e_2 x^{n-2} - \cdots + (-1)^n e_n$,则 $$ e_1 = \sum a_k,\quad e_2 = \sum_{i<j} a_i a_j,\quad e_3 = \sum_{i<j<k} a_i a_j a_k,\ \cdots $$ 即 $e_m$是根的所有$m$ 元乘积之和(初等对称多项式)。

例题 3.2:将 $\displaystyle\prod_{k=1}^{n} (x - k)$展开为$\sum$ 形式

直接展开: $$ \prod_{k=1}^{n} (x - k) = x^n - S_1 x^{n-1} + S_2 x^{n-2} - \cdots + (-1)^n S_n $$

其中 $S_1 = \sum_{k=1}^{n} k = \frac{n(n+1)}{2}$,

$S_2 = \sum_{1 \leq i < j \leq n} ij$。利用 $(\sum k)^2 = \sum k^2 + 2\sum_{i<j} ij$,得:

$$ S_2 = \frac{1}{2}\left[\left(\frac{n(n+1)}{2}\right)^2 - \frac{n(n+1)(2n+1)}{6}\right] $$

更高阶系数可使用牛顿恒等式(Newton's Identities)由幂和 $p_m = \sum k^m$ 递推得到。

3.3 混合运算中的同步化简

例题 3.3:计算 $\displaystyle\prod_{k=1}^{n} \frac{2k-1}{2k}$当$n \to \infty$ 时的渐进行为

:先转换为阶乘形式,再用斯特林公式:

$$ \begin{aligned} P_n &= \prod_{k=1}^{n} \frac{2k-1}{2k} = \frac{1 \cdot 3 \cdot 5 \cdots (2n-1)}{2 \cdot 4 \cdot 6 \cdots (2n)} \newline &= \frac{(2n)!}{(2 \cdot 4 \cdots 2n)^2} \cdot \frac{2 \cdot 4 \cdots 2n}{2 \cdot 4 \cdots 2n \cdot (\text{修正})} \end{aligned} $$

更直接的方法是利用 Wallis 乘积的已知结论:

$$ \prod_{k=1}^{n} \frac{2k-1}{2k} \sim \frac{1}{\sqrt{\pi n}} \quad (n \to \infty) $$

推导可用取对数转和 + 积分近似,最终得到上述渐进表达式。


四、含参积分与求和

4.1 积分号下求导法 (Differentiation Under the Integral Sign)

Leibniz 积分法则

$$ \frac{d}{d\alpha} \int_{a(\alpha)}^{b(\alpha)} f(x,\alpha)\thinspace{}dx = \int_{a(\alpha)}^{b(\alpha)} \frac{\partial f}{\partial \alpha}\thinspace{}dx + f(b(\alpha),\alpha) b'(\alpha) - f(a(\alpha),\alpha) a'(\alpha) $$

对于定限积分($a,b$ 为常数),简化为: $$ \frac{d}{d\alpha} \int_a^b f(x,\alpha)\thinspace{}dx = \int_a^b \frac{\partial f}{\partial \alpha}\thinspace{}dx $$

例题 4.1:利用含参积分求 $\displaystyle\int_0^{\infty} \frac{\sin x}{x}\thinspace{}dx$

:引入参数 $\alpha$,定义:

$$ I(\alpha) = \int_0^{\infty} e^{-\alpha x} \frac{\sin x}{x}\thinspace{}dx, \quad \alpha \geq 0 $$

目标:求 $I(0) = \int_0^{\infty} \frac{\sin x}{x}\thinspace{}dx$。

对 $\alpha$ 求导(积分号下求导):

$$ I'(\alpha) = \int_0^{\infty} \frac{\partial}{\partial \alpha}\left(e^{-\alpha x} \frac{\sin x}{x}\right) dx = -\int_0^{\infty} e^{-\alpha x} \sin x\thinspace{}dx $$

计算此积分(两次分部积分或利用 $\sin x = \operatorname{Im}(e^{ix})$):

$$ \int_0^{\infty} e^{-\alpha x} \sin x\thinspace{}dx = \frac{1}{\alpha^2 + 1} $$

所以 $I'(\alpha) = -\frac{1}{\alpha^2 + 1}$。

积分回去:

$$ I(\alpha) = -\arctan \alpha + C $$

当 $\alpha \to \infty$ 时,$e^{-\alpha x} \to 0$($x > 0$),故 $I(\infty) = 0$,得 $C = \frac{\pi}{2}$。

因此 $I(0) = \frac{\pi}{2}$:

$$ \boxed{\int_0^{\infty} \frac{\sin x}{x}\thinspace{}dx = \frac{\pi}{2}} $$

4.2 积分表示求和

很多经典的求和可以用积分来表示,从而借助积分的工具来求值。

调和数(Harmonic Numbers)的积分表示

$$ \boxed{H_n = \sum_{k=1}^{n} \frac{1}{k} = \int_0^1 \frac{1 - x^n}{1 - x}\thinspace{}dx} $$

推导

$$ \frac{1 - x^n}{1 - x} = 1 + x + x^2 + \cdots + x^{n-1} = \sum_{k=0}^{n-1} x^k $$

积分: $$ \int_0^1 \sum_{k=0}^{n-1} x^k\thinspace{}dx = \sum_{k=0}^{n-1} \int_0^1 x^k\thinspace{}dx = \sum_{k=0}^{n-1} \frac{1}{k+1} = \sum_{k=1}^{n} \frac{1}{k} = H_n $$

二项式系数和的积分表示(利用柯西积分公式的离散版本):

$$ \sum_{k=0}^{n} \binom{n}{k} = \frac{1}{2\pi} \int_0^{2\pi} (1 + e^{i\theta})^n e^{-i\theta \cdot 0}\thinspace{}d\theta $$

例题 4.2:用积分表示证明 $\displaystyle\sum_{k=0}^{n} \binom{n}{k} = 2^n$

证明:利用二项式定理:

$$ (1 + x)^n = \sum_{k=0}^{n} \binom{n}{k} x^k $$

令 $x = 1$,直接得到 $\sum_{k=0}^{n} \binom{n}{k} = 2^n$。

说明

虽然此题直接代值即可,但在更复杂的场合(如带权重的组合和),积分表示法(利用正交性 $\int_0^{2\pi} e^{im\theta} e^{-in\theta} d\theta = 2\pi \delta_{mn}$)是极其有力的工具。


五、级数展开法

5.1 泰勒级数展开求和 (Taylor Series)

常见函数的泰勒展开是将求和式转化为闭形式(Closed Form)的利器。

函数泰勒级数收敛半径
$e^x$$\displaystyle\sum_{n=0}^{\infty} \frac{x^n}{n!}$$\infty$
$\sin x$$\displaystyle\sum_{n=0}^{\infty} \frac{(-1)^n x^{2n+1}}{(2n+1)!}$$\infty$
$\cos x$$\displaystyle\sum_{n=0}^{\infty} \frac{(-1)^n x^{2n}}{(2n)!}$$\infty$
$\ln(1+x)$$\displaystyle\sum_{n=1}^{\infty} \frac{(-1)^{n-1} x^n}{n}$$(-1, 1]$
$\arctan x$$\displaystyle\sum_{n=0}^{\infty} \frac{(-1)^n x^{2n+1}}{2n+1}$$[-1, 1]$
$\frac{1}{1-x}$$\displaystyle\sum_{n=0}^{\infty} x^n$$(-1, 1)$

例题 5.1:求 $\displaystyle\sum_{n=1}^{\infty} \frac{1}{n \cdot 2^n}$

:观察级数形式 $\sum \frac{x^n}{n}$,取 $x = \frac{1}{2}$。

已知 $\displaystyle-\ln(1-x) = \sum_{n=1}^{\infty} \frac{x^n}{n}$,$|x| < 1$。

令 $x = \frac{1}{2}$:

$$ \sum_{n=1}^{\infty} \frac{1}{n \cdot 2^n} = -\ln\left(1 - \frac{1}{2}\right) = -\ln\frac{1}{2} = \ln 2 $$

例题 5.2:求 $\displaystyle\sum_{n=0}^{\infty} \frac{(-1)^n}{(2n+1) \cdot 3^n}$

:已知 $\displaystyle\arctan x = \sum_{n=0}^{\infty} \frac{(-1)^n x^{2n+1}}{2n+1}$。

令 $x = \frac{1}{\sqrt{3}}$,则 $x^{2n+1} = \frac{1}{3^n \sqrt{3}}$:

$$ \sum_{n=0}^{\infty} \frac{(-1)^n}{(2n+1) \cdot 3^n \sqrt{3}} = \arctan\frac{1}{\sqrt{3}} = \frac{\pi}{6} $$

因此:

$$ \boxed{\sum_{n=0}^{\infty} \frac{(-1)^n}{(2n+1) \cdot 3^n} = \frac{\pi \sqrt{3}}{6}} $$

5.2 傅里叶级数展开求和 (Fourier Series)

傅里叶级数是求特殊常数项级数和的经典方法。

例题 5.3:利用傅里叶级数证明 $\displaystyle\sum_{n=1}^{\infty} \frac{1}{n^2} = \frac{\pi^2}{6}$(巴塞尔问题)

推导:考虑函数 $f(x) = x^2$在$[-\pi, \pi]$ 上的傅里叶级数。

$f(x) = x^2$ 是偶函数,傅里叶余弦展开:

$$ f(x) = \frac{a_0}{2} + \sum_{n=1}^{\infty} a_n \cos(nx) $$

计算系数:

$$ a_0 = \frac{2}{\pi} \int_0^{\pi} x^2\thinspace{}dx = \frac{2}{\pi} \cdot \frac{\pi^3}{3} = \frac{2\pi^2}{3} $$

$$ a_n = \frac{2}{\pi} \int_0^{\pi} x^2 \cos(nx)\thinspace{}dx $$

分部积分两次:

$$ \begin{aligned} \int_0^{\pi} x^2 \cos(nx)\thinspace{}dx &= \left[\frac{x^2 \sin(nx)}{n}\right]_0^{\pi} - \frac{2}{n} \int_0^{\pi} x \sin(nx)\thinspace{}dx \newline &= 0 - \frac{2}{n} \left(\left[-\frac{x \cos(nx)}{n}\right]_0^{\pi} + \frac{1}{n} \int_0^{\pi} \cos(nx)\thinspace{}dx\right) \newline &= \frac{2}{n^2} \left[ x \cos(nx) \right]_0^{\pi} - 0 \newline &= \frac{2}{n^2} \cdot \pi \cos(n\pi) = \frac{2\pi (-1)^n}{n^2} \end{aligned} $$

因此 $a_n = \frac{2}{\pi} \cdot \frac{2\pi (-1)^n}{n^2} = \frac{4(-1)^n}{n^2}$。

傅里叶级数为:

$$ x^2 = \frac{\pi^2}{3} + 4 \sum_{n=1}^{\infty} \frac{(-1)^n}{n^2} \cos(nx), \quad x \in [-\pi, \pi] $$

在 $x = \pi$处(傅里叶级数收敛于$\frac{f(\pi^+) + f(\pi^-)}{2} = \pi^2$):

$$ \pi^2 = \frac{\pi^2}{3} + 4 \sum_{n=1}^{\infty} \frac{(-1)^n}{n^2} \cos(n\pi) $$

由于 $\cos(n\pi) = (-1)^n$,得:

$$ \pi^2 = \frac{\pi^2}{3} + 4 \sum_{n=1}^{\infty} \frac{1}{n^2} $$

解出:

$$ \boxed{\sum_{n=1}^{\infty} \frac{1}{n^2} = \frac{\pi^2}{6}} $$

例题 5.4:证明 $\displaystyle\sum_{n=1}^{\infty} \frac{(-1)^{n+1}}{n^2} = \frac{\pi^2}{12}$

在 $x = 0$ 处代入傅里叶级数:

$$ 0 = \frac{\pi^2}{3} + 4 \sum_{n=1}^{\infty} \frac{(-1)^n}{n^2} $$

故 $\displaystyle\sum_{n=1}^{\infty} \frac{(-1)^n}{n^2} = -\frac{\pi^2}{12}$,即 $\displaystyle\sum_{n=1}^{\infty} \frac{(-1)^{n+1}}{n^2} = \frac{\pi^2}{12}$。

5.3 幂级数的逐项积分与逐项求导

例题 5.5:求 $\displaystyle\sum_{n=1}^{\infty} n x^n$ 的闭形式(Closed Form)

:从几何级数出发:

$$ \sum_{n=0}^{\infty} x^n = \frac{1}{1-x}, \quad |x| < 1 $$

两边对 $x$ 求导:

$$ \sum_{n=1}^{\infty} n x^{n-1} = \frac{1}{(1-x)^2} $$

乘以 $x$:

$$ \boxed{\sum_{n=1}^{\infty} n x^n = \frac{x}{(1-x)^2}, \quad |x| < 1} $$

例题 5.6:求 $\displaystyle\sum_{n=1}^{\infty} \frac{x^n}{n}$ 的闭形式

:从 $\frac{1}{1-x} = \sum_{n=0}^{\infty} x^n$出发,两边从$0$到$x$ 积分:

$$ \int_0^x \frac{1}{1-t}\thinspace{}dt = \sum_{n=0}^{\infty} \frac{x^{n+1}}{n+1} = \sum_{n=1}^{\infty} \frac{x^n}{n} $$

左边 = $-\ln(1-x)$,故:

$$ \boxed{\sum_{n=1}^{\infty} \frac{x^n}{n} = -\ln(1-x), \quad -1 \leq x < 1} $$


六、渐进分析与估计

当精确值无法求得时,渐进分析(Asymptotic Analysis)提供主导项近似。

6.1 渐进记号

记号定义含义
$f(n) = O(g(n))$$\exists C, N: \Vert{}f(n)\Vert \leq C\Vert{}g(n)\Vert, \forall n \geq N$$f$增长不超过$g$
$f(n) = o(g(n))$$\lim_{n\to\infty} f(n)/g(n) = 0$$f$远小于$g$
$f(n) \sim g(n)$$\lim_{n\to\infty} f(n)/g(n) = 1$$f$与$g$ 渐进等价

6.2 欧拉-麦克劳林求和公式 (Euler-Maclaurin Formula)

这是连接求和与积分的核心公式:

欧拉-麦克劳林求和公式

$$ \sum_{k=1}^{n} f(k) = \int_1^n f(x)\thinspace{}dx + \frac{f(1)+f(n)}{2} + \sum_{j=1}^{m} \frac{B_{2j}}{(2j)!}\left[f^{(2j-1)}(n) - f^{(2j-1)}(1)\right] + R_m $$ 其中 $B_{2j}$ 是伯努利数(Bernoulli Numbers),$R_m$ 是余项。

伯努利数前几项:$B_2 = \frac{1}{6}$,$B_4 = -\frac{1}{30}$,$B_6 = \frac{1}{42}$,$B_8 = -\frac{1}{30}$,奇数项(除 $B_1$ 外)为零。

例题 6.1:用欧拉-麦克劳林公式估计 $\displaystyle\sum_{k=1}^{n} \sqrt{k}$ 的渐进表达式

:$f(x) = \sqrt{x} = x^{1/2}$,取 $m = 1$。

计算各项:

  1. 积分项:$\displaystyle\int_1^n x^{1/2}\thinspace{}dx = \frac{2}{3}\left(n^{3/2} - 1\right) = \frac{2}{3} n^{3/2} - \frac{2}{3}$

  2. 边界项:$\frac{f(1)+f(n)}{2} = \frac{1+\sqrt{n}}{2} = \frac{1}{2} \sqrt{n} + \frac{1}{2}$

  3. 一阶导数修正($B_2 = \frac{1}{6}$):

    $f'(x) = \frac{1}{2\sqrt{x}}$,所以 $$ \frac{B_2}{2!}\left[f'(n) - f'(1)\right] = \frac{1/6}{2}\left(\frac{1}{2\sqrt{n}} - \frac{1}{2}\right) = \frac{1}{24\sqrt{n}} - \frac{1}{24} $$

综合得到:

$$ \sum_{k=1}^{n} \sqrt{k} = \frac{2}{3} n^{3/2} + \frac{1}{2}\sqrt{n} + \left(-\frac{2}{3} + \frac{1}{2} - \frac{1}{24}\right) + \frac{1}{24\sqrt{n}} + \cdots $$

即:

$$ \boxed{\sum_{k=1}^{n} \sqrt{k} = \frac{2}{3} n^{3/2} + \frac{1}{2}\sqrt{n} - \frac{5}{24} + \frac{1}{24\sqrt{n}} + O\left(\frac{1}{n^{3/2}}\right)} $$

6.3 拉普拉斯方法 (Laplace's Method)

用于估计形如 $\int e^{n f(x)} dx$ 的大参数积分。

拉普拉斯方法

若 $f(x)$在$[a,b]$内有唯一最大值点$x_0$($f'(x_0)=0, f''(x_0)<0$),则当 $n \to \infty$: $$ \int_a^b e^{n f(x)}\thinspace{}dx \sim e^{n f(x_0)} \sqrt{\frac{2\pi}{-n f''(x_0)}} $$

例题 6.2:推导斯特林公式 (Stirling's Formula) 的首项

推导:$n! = \Gamma(n+1) = \int_0^{\infty} x^n e^{-x}\thinspace{}dx$。

令 $x = n t$,则:

$$ n! = n^{n+1} \int_0^{\infty} t^n e^{-nt}\thinspace{}dt = n^{n+1} \int_0^{\infty} e^{n(\ln t - t)}\thinspace{}dt $$

令 $f(t) = \ln t - t$,则 $f'(t) = \frac{1}{t} - 1 = 0$给出$t_0 = 1$,且 $f''(1) = -1$,$f(1) = -1$。

应用拉普拉斯方法:

$$ n! \sim n^{n+1} \cdot e^{-n} \cdot \sqrt{\frac{2\pi}{n \cdot 1}} = \sqrt{2\pi n}\left(\frac{n}{e}\right)^n $$

即斯特林公式:

$$ \boxed{n! \sim \sqrt{2\pi n}\left(\frac{n}{e}\right)^n \quad (n \to \infty)} $$

6.4 鞍点法与最速下降法简介

鞍点法(Saddle Point Method)是拉普拉斯方法在复平面上的推广。对于形如:

$$ \int_C e^{n f(z)}\thinspace{}dz $$

的围道积分,选择通过鞍点($f'(z_0) = 0$)的最速下降路径(Steepest Descent Path),将积分近似为鞍点处的高斯型积分。该方法在组合数学(生成函数系数提取)和数学物理中有广泛应用。


七、大型运算的数值计算策略

7.1 递推关系 (Recurrence Relations)

利用递推避免重复计算大指数幂和。

例题 7.1:用递推计算 $\displaystyle\sum_{k=1}^{n} k^m$(伯努利递推)

递推公式(参见 幂和公式及推导):

$$ S_m(n) = \frac{1}{m+1} \left[(n+1)^{m+1} - 1 - \sum_{j=0}^{m-1} \binom{m+1}{j} S_j(n)\right] $$

从 $S_0(n) = n$开始,可逐步递推计算任意$m$ 的幂和。

算法优势

递推只需 $O(m^2)$时间即可求得$S_0, S_1, \dots, S_m$的值,而直接求和需要$O(mn)$。

7.2 分治策略与快速算法

快速傅里叶变换(FFT)加速卷积求和

对于卷积形式的求和 $c_k = \sum_{i+j=k} a_i b_j$(离散卷积),直接计算需要 $O(n^2)$,而使用 FFT 可降至 $O(n \log n)$。

二分法在数值求和中的应用

对于单调级数,可通过积分比较进行二分截断。例如估计 $\sum_{k=n+1}^{\infty} \frac{1}{k^2}$:

$$ \int_{n+1}^{\infty} \frac{dx}{x^2} < \sum_{k=n+1}^{\infty} \frac{1}{k^2} < \int_n^{\infty} \frac{dx}{x^2} $$

即 $\frac{1}{n+1} < R_n < \frac{1}{n}$,给出精确的误差界。

7.3 误差控制

截断误差

用有限和 $\sum_{k=1}^{N} a_k$逼近无穷和$\sum_{k=1}^{\infty} a_k$时,余项$R_N = \sum_{k=N+1}^{\infty} a_k$ 需要估计。

常见估计方法

  1. 积分判别法:若 $a_k = f(k)$且$f(x)$单调递减,则$\int_{N+1}^{\infty} f(x)dx \leq R_N \leq \int_N^{\infty} f(x)dx$
  2. 交错级数:若为交错递减,则 $\Vert{}R_N\Vert \leq \Vert{}a_{N+1}\Vert$
  3. 比值法:若 $\left\Vert\frac{a_{k+1}}{a_k}\right\Vert \leq r < 1$,则几何级数给出余项上界

八、综合例题

综合题 1:极限与积分近似

求 $\displaystyle\lim_{n \to \infty} \frac{1}{n} \sum_{k=1}^{n} \sqrt{\frac{n}{k}}$

:化简求和式:

$$ \frac{1}{n} \sum_{k=1}^{n} \sqrt{\frac{n}{k}} = \frac{1}{n} \sum_{k=1}^{n} \frac{\sqrt{n}}{\sqrt{k}} = \frac{1}{\sqrt{n}} \sum_{k=1}^{n} \frac{1}{\sqrt{k}} $$

改写为黎曼和的形式:

$$ \frac{1}{\sqrt{n}} \sum_{k=1}^{n} \frac{1}{\sqrt{k}} = \frac{1}{n} \sum_{k=1}^{n} \frac{1}{\sqrt{k/n}} $$

识别为黎曼积分:

$$ \lim_{n \to \infty} \frac{1}{n} \sum_{k=1}^{n} \frac{1}{\sqrt{k/n}} = \int_0^1 \frac{1}{\sqrt{x}}\thinspace{}dx = [2\sqrt{x}]_0^1 = 2 $$

$$ \boxed{\lim_{n \to \infty} \frac{1}{n} \sum_{k=1}^{n} \sqrt{\frac{n}{k}} = 2} $$

综合题 2:交换求和与积分 + 级数展开

计算 $\displaystyle I = \int_0^1 \left(\sum_{k=1}^{\infty} \frac{x^k}{k}\right) dx$

:内层是已知级数 $\sum_{k=1}^{\infty} \frac{x^k}{k} = -\ln(1-x)$。

但也可以交换次序验证:

$$ \begin{aligned} I &= \int_0^1 \sum_{k=1}^{\infty} \frac{x^k}{k}\thinspace{}dx = \sum_{k=1}^{\infty} \frac{1}{k} \int_0^1 x^k\thinspace{}dx \newline &= \sum_{k=1}^{\infty} \frac{1}{k(k+1)} = \sum_{k=1}^{\infty} \left(\frac{1}{k} - \frac{1}{k+1}\right) \end{aligned} $$

这是裂项相消(Telescoping):

$$ I = \left(1 - \frac{1}{2}\right) + \left(\frac{1}{2} - \frac{1}{3}\right) + \left(\frac{1}{3} - \frac{1}{4}\right) + \cdots = 1 $$

用 $-\ln(1-x)$ 直接积分验证:

$$ I = \int_0^1 -\ln(1-x)\thinspace{}dx = \int_0^1 -\ln t\thinspace{}dt \quad (t = 1-x) $$

分部积分:$\int_0^1 -\ln t\thinspace{}dt = [-t\ln t + t]_0^1 = 1$。

$$ \boxed{I = 1} $$

综合题 3:乘积极限 → 取对数 → 求和 → 积分近似

求 $\displaystyle\lim_{n \to \infty} \prod_{k=1}^{n} \left(1 + \frac{k}{n^2}\right)$

此题在 3.1 节已详细推导,结果为 $\sqrt{e}$。

核心思路:乘积 → 取对数变求和 → 泰勒展开 → 黎曼和 → 积分 → 指数还原。这是处理乘积极限的标准范式。

综合题 4:证明 $\displaystyle\sum_{n=1}^{\infty} \frac{1}{n^2} = \frac{\pi^2}{6}$(巴塞尔问题)

此题在 5.2 节已有完整傅里叶级数推导。另有一种使用帕塞瓦尔恒等式(Parseval's Identity)的方法:

对 $f(x) = x$在$[-\pi, \pi]$上的傅里叶正弦展开$f(x) = 2\sum_{n=1}^{\infty} \frac{(-1)^{n+1}}{n} \sin(nx)$,应用帕塞瓦尔恒等式:

$$ \frac{1}{\pi} \int_{-\pi}^{\pi} x^2\thinspace{}dx = \sum_{n=1}^{\infty} \frac{4}{n^2} $$

左边 $= \frac{1}{\pi} \cdot \frac{2\pi^3}{3} = \frac{2\pi^2}{3}$,故 $\sum \frac{1}{n^2} = \frac{\pi^2}{6}$。

综合题 5:估计 $\displaystyle\sum_{k=1}^{n} \sqrt{k}$ 的渐进表达式

此题在 6.2 节用欧拉-麦克劳林公式已完整推导,结果为:

$$ \sum_{k=1}^{n} \sqrt{k} = \frac{2}{3} n^{3/2} + \frac{1}{2}\sqrt{n} + \zeta\left(-\frac{1}{2}\right) + o(1) $$

其中常数项与黎曼 $\zeta$ 函数在负半整数点的值有关。


知识图谱链接

相关笔记关联说明
求和式的计算方式求和的基础方法与技巧
求积式的计算方式求积的基础方法与 $\prod$转$\sum$
数列与微分方程核心方法:原理、推导及线性代数本质递推与特征根法的统一视角
微积分积分号下求导、分部积分
定积分与积分应用黎曼和、含参积分
幂和公式及推导伯努利数、欧拉-麦克劳林公式
不等式证明方法综述渐进估计中的放缩技巧
微分方程积分因子法与含参积分的联系
伯努利不等式乘积求和中的不等式放缩

方法总结

混合大型运算的核心是次序交换表示转化。面对一个混合运算问题,可按照以下思路链分析:

  1. 能否交换 $\sum$与$\int$ 的次序?(有限和直接交换 / 无穷和检查一致收敛)
  2. 能否用级数展开将求和或积分转化为已知函数?(泰勒级数 / 傅里叶级数)
  3. 乘积能否取对数化为求和?($\prod \to \sum \ln$)
  4. 能否引入参数将问题嵌入更一般的含参积分?(积分号下求导)
  5. 若精确解不可求,能否得到渐进表达式?(欧拉-麦克劳林 / 拉普拉斯方法)

掌握这五条思路链,即可应对绝大多数混合大型运算问题。

基于 Obsidian 整理 · 由 VitePress 构建