刚学会高数，xsinx/(x^2+16) 从负无穷到正无穷的积分怎么算？

题目

求解 ${\int }_{-\infty }^{+\infty }\frac{x\sin x}{x^2+16}\mathrm{d}x$

歪解

因为 $y=x$ 是奇函数，$y=\sin x$ 是奇函数，$y=\frac{1}{x^2+16}$ 是偶函数，奇函数乘上奇函数乘上偶函数为奇函数，奇函数在关于 0 对称的范围内积分等于 0，所以 ${\int }_{-\infty }^{+\infty }\frac{x\sin x}{x^2+16}\mathrm{d}x=0$，完毕（这不是很简单？）。

答案

好做完了，然后上网搜一下答案：
$$\begin{array}{rl}& x^2+16=0\\ & x^2=-16\\ & x=\pm \sqrt{-16}=\pm 4i\\ & J={\oint }_C\frac{z{e}^{iz}}{z^2+16}\mathrm{d}z=2\pi i\cdot \frac{z{e}^{iz}}{2z}{|}_{z=4i}=\pi i\cdot e^{i\cdot 4i}=i\pi e^{-4}\\ & \therefore {\int }_{-\infty }^{\infty }\frac{x\sin x}{x^2+16}\mathrm{d}x=\mathrm{Im}(J)=\pi e^{-4}\end{array}$$

啊？怎么不是 0？看到答案后脑子里冒出一个又一个问题：

• 为什么要求分母满足 $x^2+16=0$ 的解？
• 为什么一个一元函数的积分可以转化成曲线积分？曲线积分不是在平面上才有的吗？
• $z$ 跟 $x$ 是什么关系？$x\sin x=z{e}^{iz}$ 吗？
• $2\pi i$ 是从哪里来的？
• $\frac{z{e}^{iz}}{2z}$ 分母的那个 $2z$ 又是什么东西？
• 为什么只取 $z=4i$ 的值代入计算，$z=-4i$ 直接扔掉了？
• 最后那个 $\mathrm{Im}(J)$ 里 $\mathrm{Im}$ 又是什么函数？

总之是没有一步是看得懂的，这真的是刚学完高数能做的积分吗？
答案不等于 $0$ 说明这不是一个奇函数，先画个图看一看：

怎么是个偶函数？噢对不起，奇函数乘上奇函数是一个偶函数，大错特错嘞。上网搜搜，解这题需要用到什么知识点，搜到了 B 站视频：合集·高等数学（物理类专业用-四川大学编），看起来是要用到留数定理，开学！

以下为看懂本题答案需要的所有前置知识（除高数以外）（不，其实是我个人的学习笔记）。

本文中略的部分可能属于以下几种情况之一：

• 太简单
• 高数书的旧知识
• 跟本题无关
• 证明太繁琐，仅给出思路

留数定理，启动！

1.1 复数

1.1.1 复数域

复数：形如 $z=x+iy$ 的数，其中 $x,y\in R$，$i$ 为虚数单位（$i^2=-1$），称 $x$ 为 $z$ 的实部，记 $x=\mathrm{Re}(z)$，$y$ 为 $z$ 的虚部，记 $y=\mathrm{Im}(z)$，基本运算规则略（仅为说明答案中 $\mathrm{Im}$ 的含义）。

1.1.2 复平面

如下图建立平面直角坐标系，横坐标为实轴，纵坐标为虚轴，则 $z=x+iy$ 可以在平面上唯一确定一点，从 $(0,0)$ 指向 $(x,y)$ 的有向线段可表示为一条向量，借鉴向量的含义，复数的模表示为 $r=|z|$，幅角表示为 $\mathrm{Arg}(z)=\theta +2k\pi (k\in Z)$，复数的幅角不唯一，一般将 $(-\pi ,\pi ]$ 或 $[0,2\pi )$ 的幅角称为 $\mathrm{Arg}(z)$ 的主值，记为 $\arg (z)=\theta$，该平面称为复平面。

由图可得：$z=x+iy=r(\cos \theta +i\sin \theta )=r{e}^{i\theta }$，表示模长为 $r$ 的向量从实轴开始逆时针旋转 $\theta$ 角得到的向量所代表的复数为 $z=x+iy$。

特别地，$e^{i\pi }=-1⟺e^{i\pi }+1=0$，$0$、$1$（来自算术），$\pi$（来自几何），$e$（来自分析学），$i$（来自代数）这五个数以及加法、乘法、指数运算这 3 种重要的数学运算巧妙的结合在一起，优美，太优美了！

其中 $\cos \theta +i\sin \theta =e^{i\theta }$ 为欧拉公式，证略；由复数的指数形式推导模长、幅角计算公式略。

1.1.3 复数的乘幂与方根（略）

感兴趣可以学习视频：【官方双语】欧拉公式与初等群论，从另一个角度理解复数域的运算（本节与答案无关）

1.2 复变函数

1.2.1 区域与若尔当曲线（略）

本节是对复数在复平面上的拓扑概念的严格数学定义，后面看到这些概念其实可以自然理解。

1.2.2 复变函数的概念

概念、定义内容略。

过渡内容：对于复数 $z=x+iy$，$f(z)$ 总能表示为 $u(x,y)+iv(x,y)$ 的形式，其中 $u(x,y),v(x,y)$ 均为实函数，后面对于复变函数的研究总可以从二元函数角度出发进行。

1.2.3 极限与连续（略）

本节使用 $\epsilon \delta$ 语言严格定义复变函数的极限与连续，与实函数类似。

利用 $f(z)=u(x,y)+iv(x,y)$ 可将复变函数的极限、连续性问题转化为实二元函数的类似问题，相关公式、定理略。

注：由于复数不能比较大小，复变函数的连续性与有界性表现为函数值模的连续性与有界性。

2.1 解析函数的概念与柯西·黎曼条件

2.1.1 导数与微分

导数定义：设 $w=f(z)$ 在区域 $D$ 内有定义。若点 $z\in D$ 有 ${\lim }_{\Delta z\to 0}\frac{\Delta w}{\Delta z}={\lim }_{\Delta z\to 0}\frac{f(z+\Delta z)-f(z)}{\Delta z}$ 存在，则称 $f(z)$ 在 $z$ 点可导，此极限称为 $f(z)$ 在 $z$ 点的导数，记为 $f^{\prime }(z)$。

微分定义：如果 $w=f(z)$ 在点 $z\in D$ 处的改变量 $\Delta w=f(z+\Delta z)-f(z)$ 可写成 $\Delta w=A(z)\cdot \Delta z+\rho (\Delta z)$ 的形式，则称 $f(z)$ 在 $z$ 点可微。线性部分 $A(z)\cdot \Delta z$ 称为 $f(z)$ 在 $z$ 处的微分，记作 $dw=df=A(z)dz$。
其中 $\rho (\Delta z)$ 满足 ${\lim }_{\Delta z\to 0}\frac{\rho (\Delta z)}{\Delta z}=0$，即 $\rho (\Delta z)$ 是关于 $\Delta z$ 的复变函数，且它是 $\Delta z$ 的高阶无穷小。

• 注 1：复变函数可导等价于可微，证略。
• 注 2：在实函数中，$\Delta x$ 指分别从 $-\infty$ 或者 $+\infty$ 两个方向趋于 $x$ 做差，在复变函数中 $\Delta z$ 指从复平面上任意方向趋于 $z$ 做差。
• 注 3：复变函数对应的每个实二元函数 $u(x,y),v(x,y)$ 都可微，该复变函数不一定可微。反例：$f(z)=\bar{z}$，由于与答案无关，证略。

求导公式略。

2.1.2 柯西·黎曼条件

若复变函数 $f(z)=u(x,y)+iv(x,y)$ 可微，则满足 $\frac{\partial u}{\partial x}=\frac{\partial v}{\partial y}$ 且 $\frac{\partial u}{\partial y}=-\frac{\partial v}{x}$（柯西·黎曼条件，也称 C-R 条件）。

证明思路：将 $u,v$ 代入导数的定义：$f^{\prime }(z)={\lim }_{\Delta z\to 0}\frac{f(z+\Delta z)-f(z)}{\Delta z}$，分别从实轴正方向与虚轴正方向趋近 $0$ 得到的极限应该相等，得到 $f(z)=\frac{\partial u}{\partial x}+i\frac{\partial v}{\partial x}=\frac{\partial v}{\partial y}-i\frac{\partial u}{\partial y}$。

$f(z)=u(x,y)+iv(x,y)$ 在区域 $D$ 内有定义，$f(z)$ 在点 $z=x+iy$ 处可微的充要条件为：$u(x,y),v(x,y)$ 在 $(x,y)$ 处可微，并且 $u_x=v_y,u_y=-v_x$（$u_x$ 表示 $\frac{\partial u}{\partial x}$，其他符号同理）

证明思路：

• 必要性：由 $f(z)$ 可微得 $\Delta f(z)=A(z)\cdot \Delta z+\rho (\Delta z)$，设 $A(z)=a+bi$，另外可将 $\Delta f(z)=f(z+\Delta z)-f(z)$ 拆解为关于 $u(x,y),v(x,y)$ 的函数，拆解后实部虚部分别相等，可得
  $$\begin{array}{c}\{\begin{array}{l}u(x+\Delta x,y+\Delta y)-u(x,y)=a\Delta x-b\Delta y+\mathrm{Re}(\rho (\Delta z))\\ v(x+\Delta x,y+\Delta y)-v(x,y)=b\Delta x+a\Delta y+\mathrm{Im}(\rho (\Delta z))\end{array}\end{array}$$
  只需要证明 $\mathrm{Re}(\rho (\Delta z))$ 与 $\mathrm{Im}(\rho (\Delta z))$ 均为 $\sqrt{x^2+y^2}$ 的高阶无穷小，即可得到 $a=u_x=v_y,b=-u_y=v_x$，证略。
• 充分性：由 $u,v$ 在 $(x,y)$ 处可微可得式 $(1)$，其中需要分别将 $\mathrm{Re}(\rho (\Delta z))$ 与 $\mathrm{Im}(\rho (\Delta z))$ 改为 ${\eta }_1,{\eta }_2$，表示为两个二元函数微分后的高阶无穷小量，最后需要证明 ${\eta }_1+i{\eta }_2$ 是 $\Delta z$ 的高阶无穷小，即可说明 $f(z)$ 在 $z=x+iy$ 处可微，证略。

2.1.3 解析函数

解析函数定义：$f(z)$ 在点 $z_0$ 的某领域内有定义，且处处可微，则称 $z_0$ 为 $f(z)$ 的解析点，或称 $f(z)$ 在 $z_0$ 处解析。若 $f(z)$ 在区域 $D$ 内点点解析，则称 $f(z)$ 在区域 $D$ 内解析，称 $f(z)$ 是 $D$ 内的解析函数。

奇点的定义：若 $f(z)$ 在点 $z_0$ 不解析，但在 $z_0$ 的任意邻域内总有 $f(z)$ 的解析点，则称 $z_0$ 是 $f(z)$ 的奇点。

注：解析函数有非常好的光滑性，解析函数在解析点上无限次可微！（证略或者后面再证明）

2.2 解析函数与调和函数的关系（略）

本节内容如题，主要通过解析函数的 C-R 条件 + 可微性质推导其与调和函数的关系，由于与答案无关，略。

2.3 初等解析函数

2.3.1 初等单值函数（重要）

• 幂函数相关：
  • 正整数次的幂函数（$f(z)=z^n$）在复平面上解析；
  • 多项式函数（$f(z)={\sum }_{i=0}^n{a}_i{z}^i$）在复平面上解析；
  • 有理函数（$f(z)=\frac{{\sum }_{i=0}^n{a}_i{z}^i}{{\sum }_{j=0}^m{b}_j{z}^j}$）在所有满足 ${\sum }_{j=0}^m{b}_j{z}^j=0$ 的点上为奇点，在复平面上除奇点外均解析，2.1.2 节的充要条件可证，证略；
• 指数函数（$f(z)=e^z$）在复平面上解析，证略；
• 三角函数 $\cos z=\frac{e^{iz}+e^{-iz}}{2}$，$\sin z=\frac{e^{iz}-e^{-iz}}{2i}$（欧拉公式可得）在复平面上解析，证略；
• 相关公式与当 $z$ 为纯实数时的形式相同（如 $(z^n{)}^{\prime }=n{z}^{n-1},(e^z{)}^{\prime }=e^z,(\sin z{)}^{\prime }=\cos z$），而领一些与实函数不同的性质由于与答案无关，略；证略；

到这里我们终于可以看懂答案的第一步是在做什么了：由于被积函数 $\frac{x\sin x}{x^2+16}$ 中分子在复平面上解析，分母是一个有理函数，因此求分母满足 $x^2+16=0$ 的解，即找出被积函数在复平面上的奇点（这里变量为 $x$ 或者 $z$ 并无区别），之后围绕奇点做一些计算。
搞懂了第一步在算奇点，但还不知道为什么要算这个奇点，所以只能算是向成功迈出了 0.5 步吧（终于有进展了，泪目/(ㄒoㄒ)/~~）

2.3.2 初等多值函数

复数幅角的多值性导致复数函数的多值性（如函数 $f(z)=z^{\frac{1}{3}}$ 当自变量 $z=-1$ 时，因变量可以为 $e^{i\frac{\pi +2k\pi }{3}},k\in Z$），对于多值函数无法讨论在某一点上的解析、可微等性质，因此本节引入单值分支的概念限制值域、支割线的概念限制定义域，将多值函数限制在单值函数中，再进行后续内容的学习。

3.1 复积分的概念与性质

3.1.1 定义

设 $C$ 是一条以 $z_0$ 为起点，$z^{\prime }$ 为终点的有向曲线，函数 $f(z)$ 在 $C$ 上有定义。沿着 $C$ 的正向依次取 $z_0,z_1,\cdots ,z_n=z^{\prime }$，考虑 ${\sum }_{i=1}^{n}f({\xi }_i)\Delta z_i$，其中 $\Delta z_i=z_i-z_{i-1}$，${\xi }_i$ 为弧段 $\stackrel{⌢}{z_{i-1}{z}_i}$ 上的一点。

设 λ = max ⁡ 1 ≤ i ≤ n { ∣ Δ z i ∣ } \lambda=\max_{1\leq i\leq n}\{|\Delta z_i|\} λ=max1≤i≤n​{∣Δzi​∣}，若极限 ${\lim }_{\lambda \to 0}{\sum }_{i=1}^{n}f({\xi }_i)\Delta z_i$ 存在，且与划分方式、${\xi }_i$ 选取方式无关，则称 $f(z)$ 沿 $C$ 可积。

定义 ${\int }_{C}f(z)\mathrm{d}z={\lim }_{\lambda \to 0}{\sum }_{i=1}^{n}f({\xi }_i)\Delta z_i$ 为 $f(z)$ 沿 $C$ 的积分，$C$ 称为积分路径。

以上定义方式严谨但不好理解，这里做一些简单的解释以帮助理解这个积分的含义：
假设有一个人在复平面上从 $z_0=(x_0,y_0)$ 点沿着曲线 $C$ 的路径方向走到 $z^{\prime }=(x^{\prime },y^{\prime })$，在这个过程中他用大小与方向不断变化着的力 $F$ 推着一个箱子前进，而这个力 $F$ 是关于他所经过路径上的点 $z$ 的函数 $F=f(z)$，我们知道力在路径上的做功为 $W=Fs$，当 $F$ 为恒力且 $s$ 为线段时可简化为直接相乘，而变力作用在曲线路径的情况下，需要将 $s$ 细分为无穷多小段 $s_i$，在其中的每一段 $F$ 都是一个关于箱子所在位置的函数，即 $W=Fs={\lim }_{s_i\to 0}{\sum }_{i=1}^{n}f(z_i)s_i={\int }_{C}f(z)\mathrm{d}z$，在这里积分的结果就是变力 $F$ 沿着曲线路径 $C$ 所做的功。
注意：这里用一个复数 $z=x+iy$ 表示在平面上的位置，而用另一个复数 $f(z)$ 表示力的大小与方向（向量）。

注 1：简单闭曲线的正向（下图可帮助理解）：

• 面朝正向，左手指向内部，右手指向外部（二维平面适用）
• 外法向为 $x$ 轴方向，正向为 $y$ 轴方向（可拓展至高维空间）

注 2：积分路径不同，积分值可能不同，例：积分 ${\int }_C\mathrm{Re}(z)\mathrm{d}z$，沿着以下两条路径积分值不同

• $C$ 为链接 $0$ 与 $1+i$ 的直线段，积分值为 $\frac{1+i}{2}$
• $C$ 为链接 $0$ 与 $1$，$1$ 与 $1+i$ 的折线段，积分值为 $\frac{1}{2}+i$

TODO: 补充一张图，参考 013 节板书

若能将 $z$ 写成关于 $t$ 的参数方程 $z=z(t)$，则 ${\int }_{C}f(z)\mathrm{d}z={\int }_{\alpha }^{\beta }f(z(t))z^{\prime }(t)\mathrm{d}t$，推导略。

★★★ 计算以下积分，曲线 $C$ 为以 $a$ 为中心，$\rho (\rho >0)$ 为半径的圆周（非常非常重要）：

$$I={\oint }_C\frac{\mathrm{d}z}{(z-a{)}^n}=\begin{array}{c}\{\begin{array}{ll}2\pi i,& \text{当}n=1\text{时}\\ 0,& \text{当}n\ne 1\text{且 n 为整数时}\end{array}\end{array}$$

考虑 $z\in C$，设 $z(\theta )=a+\rho e^{i\theta },\theta \in [0,2\pi )$，则 $\mathrm{d}z=i\rho e^{i\theta }\mathrm{d}\theta$

$$\begin{array}{c}\begin{array}{rl}I& ={\oint }_C\frac{\mathrm{d}z}{(z-a{)}^n}\\ & ={\int }_0^{2\pi }\frac{i\rho e^{i\theta }\mathrm{d}\theta }{(a+\rho e^{i\theta }-a{)}^n}\\ & ={\int }_0^{2\pi }\frac{i\rho e^{i\theta }}{{\rho }^n{e}^{in\theta }}\mathrm{d}\theta \\ & =i{\rho }^{1-n}{e}^{i\theta (1-n)}{|}_0^{2\pi }\end{array}\end{array}$$

当 $n=1$ 时，$I=2\pi i$，当 $n\ne 1$ 且 $n$ 为整数时：

$$\begin{array}{c}\begin{array}{rl}I& =i{\rho }^{1-n}(e^{2\pi i(1-n)}-1)\\ & =i{\rho }^{1-n}(\cos (2\pi (1-n))+i\sin (2\pi (1-n))-1)\\ & =i{\rho }^{1-n}(1+i\times 0-1)\\ & =0\end{array}\end{array}$$

到这里回顾一下答案，可以发现已经出现了答案中 ${\oint }_{C}f(z)\mathrm{d}z=2\pi ig(z){|}_{z=4i}$ 这一行的部分关键内容。我们大胆猜测，这是在 $f(z)$ 的其中一个奇点 $z=4i$ 附近找了一条半径任意小的圆形积分路径，对这条路径做积分，但我们仍不清楚这么做的原因是什么。

3.1.2 性质

• ${\int }_C\mathrm{d}z=z_n-z_1$，$z_1,z_n$ 分别为曲线 $C$ 的起点与终点；
• 线性：${\int }_C(a_1{f}_1(z)+a_2{f}_2(z))\mathrm{d}z=a_1{\int }_C{f}_1(z)\mathrm{d}z+a_2{\int }_C{f}_2(z)\mathrm{d}z$，其中 $a_1,a_2$ 为常数；
• 区域可加性：${\int }_{C}f(z)\mathrm{d}z={\int }_{C_1}f(z)\mathrm{d}z+{\int }_{C_2}f(z)\mathrm{d}z$，其中曲线 $C$ 由 $C_1,C_2$ 首尾连接；
• ${\int }_{C}f(z)\mathrm{d}z={\int }_{C^{-}}f(z)\mathrm{d}z$，其中曲线 $C^{-}$ 为与 $C$ 方向相反的同一曲线；
• $|{\int }_{C}f(z)\mathrm{d}z|\le {\int }_C|f(z)||\mathrm{d}z|={\int }_C|f(z)|\mathrm{d}s$，其中 $\mathrm{d}s$ 仅表示为曲线的曲线的一小段长度微元；
• 若 $|f(z)|\le M$，曲线 $C$ 的长度 ${\int }_C\mathrm{d}s\le \rho$，则 $|{\int }_{C}f(z)\mathrm{d}z|\le \rho M$

3.2 柯西积分定理

3.2.1 柯西积分定理