---
url: /integral/integral-19/index.md
---
## Part 1 积分应用

### · 曲线的长度

设 $I = \[\alpha, \beta]$ 为区间，映射 $\sigma: I \to \mathbb{R}^2$ 用分量表示为

$$
\sigma(t) = (x(t), y(t)), \quad t \in I.
$$

如果 $x(t), y(t)$ 均为连续函数，则称 $\sigma$ 为 $\mathbb{R}^2$ 上的连续曲线.如果 $x(t), y(t)$ 均可微（连续可微），则称 $\sigma$ 为可微（连续可微）曲线.

设 $\sigma$ 为连续可微曲线，通过分割曲线并用直线段长度之和作逼近，我们可以定义 $\sigma$ 的长度为

$$
L(\sigma) = \int\_\alpha^\beta \left\[(x'(t))^2 + (y'(t))^2\right]^{\frac{1}{2}} dt.
$$

这个公式可以如下推导.首先注意到下面的简单不等式：

$$
|\sqrt{a^2 + b^2} - \sqrt{a^2 + c^2}| \leq |b - c|, \quad \forall a,b,c \in \mathbb{R}.
$$

我们将 $\[\alpha, \beta]$ 分割为 $\alpha = t\_0 < t\_1 < \cdots < t\_n = \beta$，点 $(x(t\_i), y(t\_i))$ 把曲线分成若干段，每一段的长度可以近似地用直线段的长度表示，即

$$
L(\sigma) \approx \sum\_{i=1}^{n} \sqrt{(x(t\_i) - x(t\_{i-1}))^2 + (y(t\_i) - y(t\_{i-1}))^2},
$$

由微分中值定理，存在 $\xi\_i, \eta\_i \in (t\_{i-1}, t\_i)$，使得

$$
x(t\_i) - x(t\_{i-1}) = x'(\xi\_i)(t\_i - t\_{i-1}), \quad y(t\_i) - y(t\_{i-1}) = y'(\eta\_i)(t\_i - t\_{i-1}),
$$

从而有

$$
\sqrt{(x(t\_i) - x(t\_{i-1}))^2 + (y(t\_i) - y(t\_{i-1}))^2} = \sqrt{(x'(\xi\_i))^2 + (y'(\eta\_i))^2} \Delta t\_i.
$$

因为

$$
\left|\sqrt{(x'(\xi\_i))^2 + (y'(\eta\_i))^2} \Delta t\_i - \sqrt{(x'(\xi\_i))^2 + (y'(\xi\_i))^2} \Delta t\_i\right| \leq |y'(\eta\_i) - y'(\xi\_i)| \Delta t\_i,
$$

而

$$
\sum\_{i=1}^{n} |y'(\eta\_i) - y'(\xi\_i)| \Delta t\_i \leq \sum\_{i=1}^{n} \omega\_i(y) \Delta t\_i \to 0, \quad (|\pi| = \max{t\_i - t\_{i-1}} \to 0)
$$

因此有

$$
\begin{aligned}
L(\sigma) &= \lim\_{|\pi| \to 0} \sum\_{i=1}^{n} \sqrt{(x(t\_i) - x(t\_{i-1}))^2 + (y(t\_i) - y(t\_{i-1}))^2} \\\\
&= \lim\_{|\pi| \to 0} \sum\_{i=1}^{n} \sqrt{(x'(\xi\_i))^2 + (y'(\xi\_i))^2} \Delta t\_i \\\\
&= \int\_\alpha^\beta \left\[(x'(t))^2 + (y'(t))^2\right]^{\frac{1}{2}} dt.
\end{aligned}
$$

注：如果 $(x'(t))^2 + (y'(t))^2 \neq 0$，令

$$
s = \phi(t) = \int\_\alpha^t \left\[(x'(u))^2 + (y'(u))^2\right]^{\frac{1}{2}} du, \quad t \in \[\alpha, \beta].
$$

则 $\phi: \[\alpha, \beta] \to \[0, L(\sigma)]$ 是严格单调递增函数，从而可逆，其逆记为 $t = \psi(s)$，$s$ 称为 $\sigma$ 的弧长参数.记 $\tilde{\sigma}(s) = \sigma(\psi(s))$，$s \in \[0, L(\sigma)]$.根据反函数的求导公式易见

$$
|\tilde{\sigma}'(s)| = \sqrt{(\tilde{x}'(s))^2 + (\tilde{y}'(s))^2} = 1.
$$

### · 简单图形面积

(1) 如果 $f > 0$ 为 $\[a, b]$ 上的连续函数，则由 $y = f(x)$, $x = a$, $x = b$ ($a < b$) 与 $y = 0$ 围成的曲边梯形的面积为

$$
S = \int\_a^b f(x) dx.
$$

一般地，当 $f$ 变号时，上式仍有意义，称为代数面积和，而

$$
S = \int\_a^b |f(x)| dx
$$

才是所围面积之和.更一般地，由 $y = f\_2(x)$, $y = f\_1(x)$ 以及 $x = a$, $x = b$ 围成的图形的面积为

$$
S = \int\_a^b |f\_2(x) - f\_1(x)| dx.
$$

(2) 设 $\sigma$ 为平面曲线，由极坐标方程

$$
r = r(\theta), \quad \theta \in \[\alpha, \beta]
$$

给出，其中 $r(\theta)$ 关于 $\theta$ 连续，$\beta - \alpha \leq 2\pi$.则由 $\sigma$, $\theta = \alpha$, $\theta = \beta$ 所围成的图形面积为

$$
S = \lim\_{|\pi| \to 0} \sum\_{i=1}^{m} \frac{1}{2} r^2(\xi\_i) \cdot \Delta \theta\_i = \frac{1}{2} \int\_\alpha^\beta r^2(\theta) d\theta.
$$

这个公式是通过使用扇形的面积和逼近图形面积得到的.

(3) 如果曲线 $\sigma$ 由 $\sigma(t) = (x(t), y(t))$, $t \in \[\alpha, \beta]$ 给出，其中 $y(t) \geq 0$, $x$ 关于 $t$ 单调递增，$x(\[\alpha, \beta]) = \[a, b]$.则 $\sigma$ 与 $x = a$, $x = b$ 以及 $y = 0$ 围成的曲边梯形的面积为

$$
S = \int\_\alpha^\beta y(t)x'(t) dt.
$$

这个公式仍然是通过使用矩形面积之和去逼近曲边梯形得到.一般地，如果只设 $x$ 是单调的，则面积公式为

$$
S = \int\_\alpha^\beta |y(t)x'(t)| dt.
$$

如果 $\sigma$ 除在 $t = \alpha, \beta$ 处以外无自交点，则 $\sigma$ 本身围成的图形的面积为

$$
S = \left| \int\_\alpha^\beta y(t)x'(t) dt \right| = \left| \int\_\alpha^\beta x(t)y'(t) dt \right|,
$$

因为

$$
\begin{aligned}
\int\_\alpha^\beta y(t)x'(t) dt &= y(t)x(t)\Big|*\alpha^\beta - \int*\alpha^\beta y'(t)x(t) dt \\\\
&= - \int\_\alpha^\beta y'(t)x(t) dt,
\end{aligned}
$$

故这个面积公式也可以改写为

$$
S = \frac{1}{2} \left| \int\_\alpha^\beta \[y(t)x'(t) - y'(t)x(t)] dt \right|.
$$

(4) 旋转曲面的面积

设 $\sigma$ 为平面曲线

$$
\sigma(t) = (x(t), y(t)), \quad t \in \[\alpha, \beta], \quad y(t) \geq 0.
$$

$\sigma$ 绕 $x$ 轴旋转所得曲面的面积为

$$
S = \int\_\alpha^\beta 2\pi y(t)\left\[(x'(t))^2 + (y'(t))^2\right]^{\frac{1}{2}} dt.
$$

> 取 $\[\alpha, \beta]$ 的一个分割，在分点 $t\_{i-1}, t\_i$ 之间的曲线段经过旋转后所形成的曲面的面积可以用圆台的面积近似逼近，这一部分圆台的面积为
>
> $$
> \pi(y(t\_{i-1}) + y(t\_i)) \sqrt{(x(t\_i) - x(t\_{i-1}))^2 + (y(t\_i) - y(t\_{i-1}))^2},
> $$
>
> 因此
>
> $$
> S \approx \sum\_{i=1}^{n} \pi(y(t\_{i-1}) + y(t\_i)) \sqrt{(x(t\_i) - x(t\_{i-1}))^2 + (y(t\_i) - y(t\_{i-1}))^2},
> $$
>
> 和曲线弧长公式的推导过程类似，当分割的模趋于零时，我们近似地有
>
> $$
> (y(t\_{i-1}) + y(t\_i)) \approx 2y(\xi\_i), \quad (\xi\_i \in \[t\_{i-1}, t\_i])
> $$
>
> 以及
>
> $$
> \sqrt{(x(t\_i) - x(t\_{i-1}))^2 + (y(t\_i) - y(t\_{i-1}))^2} \approx \sqrt{(x'(\xi\_i))^2 + (y'(\xi\_i))^2} \Delta t\_i,
> $$
>
> 当分割的模趋于零时，近似逼近所引起的这些误差之和趋于零.因此有
>
> $$
> \begin{aligned}
> S &= \lim\_{|\pi| \to 0} \sum\_{i=1}^{n} 2\pi y(\xi\_i)\left\[(x'(\xi\_i))^2 + (y'(\xi\_i))^2\right]^{\frac{1}{2}} \Delta t\_i \\\\
> &= \int\_\alpha^\beta 2\pi y(t)\left\[(x'(t))^2 + (y'(t))^2\right]^{\frac{1}{2}} dt.
> \end{aligned}
> $$

### · 简单立体体积

(1) 平行截面之间的立体体积

设 $\Omega$ 为 $\mathbb{R}^3$ 中一块立体区域，夹在平面 $x = a$ 与 $x = b$ ($a < b$) 之间.记 $S(x)$ 为 $x \in \[a, b]$ 处垂直于 $x$ 轴的平面截 $\Omega$ 的截面面积函数.如果 $S(x)$ 关于 $x$ 连续，则 $\Omega$ 的体积为

$$
V = \int\_a^b S(x) dx.
$$

特别地，如果两块区域 $\Omega\_A$ 和 $\Omega\_B$ 的截面面积函数相等，则其体积相同.这个事实在公元 5 到 6 世纪由祖暅（祖冲之之子）所发现，17 世纪时意大利人 Cavalieri 也发现了这一事实.

(2) 旋转体的体积

设 $f$ 为 $\[a, b]$ 上的连续函数，$\Omega$ 是由平面图形

$$
{(x, y) \mid a \leq x \leq b, ; 0 \leq |y| \leq |f(x)|}
$$

绕 $x$ 轴旋转一周所得旋转体.该旋转体在 $x \in \[a, b]$ 处的截面为圆盘，其面积为

$$
S(x) = \pi f^2(x).
$$

因此 $\Omega$ 的体积为

$$
V = \int\_a^b S(x) dx = \pi \int\_a^b f^2(x) dx.
$$

### · 物理应用

(1) 降落伞的原理

质量为 $m$ 的物体在重力作用下自由下落，下落时所受空气阻力与下落速度成正比，比例常数为 $k$，则由牛顿定律，

$$
mg - kv = m \frac{dv}{dt},
$$

其中，$g$ 为重力加速度，$v$ 为物体的速度，我们选择指向地心的坐标.上面的方程等价于

$$
\frac{d}{dt}(e^{\frac{k}{m}t} v) = g e^{\frac{k}{m}t},
$$

假设初速度为零，则

$$
e^{\frac{k}{m}t} v = g \int\_0^t e^{\frac{k}{m}s} ds = \frac{mg}{k}(e^{\frac{k}{m}t} - 1),
$$

即

$$
v(t) = \frac{mg}{k}(1 - e^{-\frac{k}{m}t}).
$$

特别地，$t \to \infty$ 时 $v(t) \to \dfrac{mg}{k}$，即速度不会增加到无限大.

(2) 第二宇宙速度

从地球表面发射火箭，如果要求火箭无限飞离地球，问：火箭的初速度至少为多大？

根据万有引力定律，在距地心 $x$ 处火箭所受地球引力为

$$
F = GMm x^{-2},
$$

其中，$G$ 为万有引力常数，$M$ 为地球质量，$m$ 为火箭质量.在地球表面，有

$$
GMm R^{-2} = mg,
$$

其中 $R$ 为地球半径.火箭从地面升到距地心 $r$ ($r > R$) 处需要做的功为

$$
\int\_R^r GMm x^{-2} dx = \int\_R^r mg R^2 x^{-2} dx = mg R^2 \left(\frac{1}{R} - \frac{1}{r}\right).
$$

因此，火箭无限飞离地球需要做功

$$
W = \lim\_{r \to \infty} mg R^2 \left(\frac{1}{R} - \frac{1}{r}\right) = mg R.
$$

由能量守恒原理，火箭的初速度至少为 $v\_0$，则

$$
\frac{1}{2} m v\_0^2 = mg R,
$$

因而

$$
v\_0 = \sqrt{2gR} \approx \sqrt{2 \times 9.81 (m/s^2) \times 6.371 \times 10^6 m} \approx 11.2 , (km/s).
$$

(3) 缆绳的工作原理

绳索在日常生活中应用十分广泛，例如在码头上经常用来系住船舶.为什么绳索能拉住大型船舶？下面我们就来作一个力学分析，它揭示了绳索产生巨大拉力的原理.

设一段绳索缠绕在一圆柱体上，绳索一端施以拉力 $f$，绳索与圆柱体之间的摩擦系数为 $k$，如果绳索共绕了 $n$ 圈，在绳索的另一端产生的拉力为 $F$，我们来求 $F$ 的值.

取角度为 $\Delta\theta$ 的一小段绳索，研究其受力状况.设这一段绳索承受圆柱体的正压力为 $\Delta N$，则摩擦力为 $k \Delta N$.这一段绳索两端所受拉力分别为 $F$, $F + \Delta F$，则考虑沿圆柱体外法向和切向这两个方向绳索的受力，得到方程

$$
\begin{cases}
\Delta N = (F + \Delta F) \sin \dfrac{\Delta\theta}{2} + F \sin \dfrac{\Delta\theta}{2}, \\\\
(F + \Delta F) \cos \dfrac{\Delta\theta}{2} = F \cos \dfrac{\Delta\theta}{2} + k \Delta N.
\end{cases}
$$

从方程中消去 $\Delta N$，令 $\Delta\theta \to 0$，得

$$
\frac{dF}{d\theta} = \lim\_{\Delta\theta \to 0} \frac{\Delta F}{\Delta\theta} = kF,
$$

利用积分解得

$$
F(\theta) = f \cdot e^{k\theta}.
$$

当 $\theta = 2n\pi$ 时，$F = f \cdot e^{2kn\pi}$.例如，设摩擦系数 $k = \dfrac{1}{4}$，$n = 6$，$f = 10 kg$，则 $F = 10 e^{3\pi} kg > 100000 kg$.

### · 近似计算

设 $f$ 为 $\[a,b]$ 上的二次连续可微函数，则由微分中值定理或 Taylor 展开我们得到

$$
|f(x) - l(x)| \leq \frac{1}{2} M (x - a)(b - x), \quad \forall x \in \[a,b],
$$

其中，$M = \max\_{x \in \[a,b]} |f''(x)|$，且

$$
l(x) = f(a) + \frac{f(b) - f(a)}{b - a}(x - a), \quad x \in \[a,b].
$$

因此有如下的积分估计

$$
\begin{aligned}
\left| \int\_a^b f(x) dx - \int\_a^b l(x) dx \right|&\leq \frac{1}{2} M \int\_a^b (x - a)(b - x) dx \\\\
&= \frac{1}{12} M (b - a)^3.
\end{aligned}
$$

这也就是 $f$ 在 $\[a,b]$ 上的积分用梯形面积逼近的误差公式.

我们考虑函数 $f = \ln x$ 在 $\[1,n]$ 上的积分.令

$$
\begin{aligned}
A\_n = \int\_1^n \ln x , dx &= x \ln x \Big|\_1^n - \int\_1^n (\ln x)' x dx \\\\
&= n \ln n - n + 1,
\end{aligned}
$$

$$
\begin{aligned}
B\_n &= \frac{1}{2}(\ln 1 + \ln 2) + \frac{1}{2}(\ln 2 + \ln 3) + \cdots + \frac{1}{2}(\ln(n-1) + \ln n)\\\\
&= \ln n! - \frac{1}{2} \ln n,
\end{aligned}
$$

根据上面的误差估计，并注意 $\ln x$ 为凹函数，则有

$$
0 < \int\_k^{k+1} \ln x , dx - \frac{1}{2}(\ln k + \ln(k+1)) < \frac{1}{12 k^2}.
$$

令 $C\_n = A\_n - B\_n$，则 $C\_n$ 是 $n-1$ 次累计误差，它关于 $n$ 是单调递增的.从而

$$
0 < C\_n < \frac{1}{12} \sum\_{k=1}^{n-1} \frac{1}{k^2} < \frac{1}{12} \left(1 + \sum\_{k=1}^\infty \left(\frac{1}{k} - \frac{1}{k+1}\right)\right) = \frac{1}{6},
$$

这说明极限 $\displaystyle\lim\_{n \to \infty} C\_n = C$ 存在，且

$$
\begin{aligned}
0 < C - C\_n &< \frac{1}{12} \sum\_{k=n}^\infty \frac{1}{k^2} \\\\
&< \frac{1}{12} \left\[\frac{1}{n^2} + \frac{1}{n(n+1)} + \frac{1}{(n+1)(n+2)} + \cdots \right] \\\\
&= \frac{1}{12} \left(\frac{1}{n^2} + \frac{1}{n}\right).
\end{aligned}
$$

下面我们来求极限 $C$ 的值.由定义，有

$$
C\_n = A\_n - B\_n = n \ln n - n + 1 - \ln n! + \frac{1}{2} \ln n,
$$

因此

$$
n! = e^{1 - C\_n} n^{n + \frac{1}{2}} e^{-n}.
$$

由 Wallis 公式，

$$
\lim\_{n \to \infty} \frac{(n!)^2 2^{2n}}{(2n)!} \cdot \frac{1}{\sqrt{n}} = \lim\_{n \to \infty} \frac{(2n)!!}{(2n - 1)!!} \cdot \frac{1}{\sqrt{n}} = \sqrt{\pi},
$$

将 $n!$ 和 $(2n)!$ 的表达式代入，有

$$
\sqrt{\pi} = \lim\_{n \to \infty} \frac{e^{2(1 - C\_n)} n^{2n+1} e^{-2n} \cdot 2^{2n}}{e^{1 - C\_{2n}} (2n)^{2n + \frac{1}{2}} e^{-2n} \cdot \sqrt{n}} = \frac{e^{1 - C}}{\sqrt{2}},
$$

这就得到 $n!$ 的如下表示

$$
n! = \sqrt{2\pi n} \left(\frac{n}{e}\right)^n e^{C - C\_n}, \quad \text{(Stirling)}
$$

其中

$$
1 < e^{C - C\_n} < e^{\frac{1}{12} \left(\frac{1}{n^2} + \frac{1}{n}\right)} < 1 + \frac{1}{12n} + \frac{1}{10n^2}, \quad \forall n > 1.
$$