第 12 章 Fourier 分析
前面各章把函数在空间区域上积分;本章换一条思路——把「足够好的」函数写成三角函数(或复指数)的叠加,在频率域里分析信号与周期现象。Fourier 级数处理周期函数;平方平均收敛给出 $L^2$ 意义下「最佳逼近」;Dirichlet 定理回答何时级数在每一点收敛;Fourier 积分与变换把周期推广到全直线,是偏微分方程、信号处理、量子力学的共同语言。
本章主线:正交性 $\Rightarrow$ 系数公式 $\Rightarrow$ 级数展开 $\Rightarrow$ 三种收敛(一致 / 逐点 / 平方平均)$\Rightarrow$ 非周期情形的 Fourier 变换。计算题要分清:周期延拓方式、奇偶性、半区间展开类型;理论题要分清:$L^2$ 收敛条件最宽,Dirichlet 条件保证逐点收敛到 $\tfrac12[f_++f_-]$。
12.1 Fourier 级数
12.1.1 三角函数系的正交性
在区间 $[-\pi,\pi]$ 上,考虑函数系 $$1,\ \cos x,\ \cos 2x,\ \ldots,\ \cos nx,\ \ldots;\qquad \sin x,\ \sin 2x,\ \ldots,\ \sin nx,\ \ldots.$$ 对任意 $m,n\in\mathbb{N}$($m,n\ge 1$),有 $$\int_{-\pi}^{\pi}\cos mx\,\cos nx\,dx= \begin{cases} \pi, & m=n,\\ 0, & m\neq n, \end{cases}\qquad \int_{-\pi}^{\pi}\sin mx\,\sin nx\,dx= \begin{cases} \pi, & m=n,\\ 0, & m\neq n, \end{cases}$$ $$\int_{-\pi}^{\pi}\sin mx\,\cos nx\,dx=0,\qquad \int_{-\pi}^{\pi}\cos nx\,dx=0,\ \int_{-\pi}^{\pi}\sin nx\,dx=0.$$ 此外 $\displaystyle\int_{-\pi}^{\pi}1\cdot 1\,dx=2\pi$。
验证一例:$m\neq n$ 时, $$\int_{-\pi}^{\pi}\cos mx\cos nx\,dx=\frac12\int_{-\pi}^{\pi}\big[\cos(m+n)x+\cos(m-n)x\big]\,dx=0,$$ 因 $m\pm n\neq 0$,$\cos kx$ 在完整周期上积分为零。$m=n$ 时 $\cos^2 nx=\tfrac12(1+\cos 2nx)$,得 $\pi$。$\sin mx\cos nx$ 可写成 $\tfrac12[\sin(m+n)x+\sin(m-n)x]$,两项在 $[-\pi,\pi]$ 上亦为零。
内积与正交:对平方可积函数 $f,g$,定义 $$\langle f,g\rangle=\frac{1}{\pi}\int_{-\pi}^{\pi} f(x)g(x)\,dx.$$ 则上述系在加权意义下两两正交;$\{\cos nx\}_{n\ge0}$、$\{\sin nx\}_{n\ge1}$ 各自构成正交系($\cos 0=1$ 单独计)。
正交性的计算技巧
利用 $\cos mx\cos nx=\tfrac12[\cos(m+n)x+\cos(m-n)x]$ 等积化和差,把乘积化为单项三角函数在完整周期上的积分。$\int_{-\pi}^{\pi}\cos kx\,dx=0$($k\neq0$)来自「正负半周抵消」;$\int_{-\pi}^{\pi}\cos^2 kx\,dx=\pi$ 可用倍角公式 $\cos^2 kx=\tfrac12(1+\cos 2kx)$。
归一化写法:令 $\phi_0=1/\sqrt2$,$\phi_n=\cos nx$,$\psi_n=\sin nx$($n\ge1$),则 $\langle\phi_m,\phi_n\rangle=\delta_{mn}$,$\langle\psi_m,\psi_n\rangle=\delta_{mn}$,$\langle\phi_m,\psi_n\rangle=0$。Fourier 系数可写 $a_n/\sqrt2=\langle f,\phi_n\rangle$($n\ge0$),$b_n=\langle f,\psi_n\rangle$,与有限维正交基坐标一致。
由正交性导出系数公式:设 $f$ 在 $[-\pi,\pi]$ 上可展开为 $$f(x)\sim\frac{a_0}{2}+\sum_{n=1}^{\infty}\big(a_n\cos nx+b_n\sin nx\big).$$ 形式上把 $f$ 与级数等同,两边同乘 $\cos kx$ 并在 $[-\pi,\pi]$ 上积分(逐项积分在收敛性未证前只是形式操作,但导出正确的候选系数): $$\int_{-\pi}^{\pi}f(x)\cos kx\,dx=\frac{a_0}{2}\underbrace{\int_{-\pi}^{\pi}\cos kx\,dx}_{=0}+\sum_{n=1}^{\infty}a_n\underbrace{\int_{-\pi}^{\pi}\cos nx\cos kx\,dx}_{=\pi\delta_{nk}}+\sum_{n=1}^{\infty}b_n\underbrace{\int_{-\pi}^{\pi}\sin nx\cos kx\,dx}_{=0}=\pi a_k.$$ 故 $a_k=\dfrac{1}{\pi}\int_{-\pi}^{\pi} f(x)\cos kx\,dx$($k\ge 1$)。$k=0$ 时与 $1$ 作内积: $$\int_{-\pi}^{\pi}f(x)\,dx=\frac{a_0}{2}\cdot 2\pi=a_0\pi\quad\Rightarrow\quad a_0=\frac{1}{\pi}\int_{-\pi}^{\pi} f(x)\,dx.$$ 与 $\sin kx$ 作内积得 $$b_k=\frac{1}{\pi}\int_{-\pi}^{\pi} f(x)\sin kx\,dx,\quad k\ge 1.$$ 上述公式对任意可积 $f$ 均有意义;它们是否给出 $f$ 的「正确展开」取决于收敛性(§12.2、§12.3)。
计算技巧:若 $f$ 为偶函数,则 $b_n=0$,只需算 $a_n$;若 $f$ 为奇函数,则 $a_n=0$(含 $a_0$),只需算 $b_n$。若 $f(\pi-x)=f(x)$,则仅含 $\cos 2nx,\sin 2nx$ 等偶次谐波;$f(\pi-x)=-f(x)$ 则仅含奇次谐波。
[要点]
- 正交性 ⇒ 系数是「在各基函数方向上的投影」。
- $a_n,b_n$ 由积分公式唯一确定,与 $f$ 是否收敛到级数无关。
- 计算时先检验 $f$ 的奇偶性,常可化简积分。
12.1.2 周期函数的 Fourier 级数
设 $f$ 以 $2\pi$ 为周期:$f(x+2\pi)=f(x)$。若 $f$ 在 $[-\pi,\pi]$ 上 Riemann 可积(或 Lebesgue 可积),其 Fourier 级数 为 $$S_f(x)=\frac{a_0}{2}+\sum_{n=1}^{\infty}\big(a_n\cos nx+b_n\sin nx\big),$$ 系数由 §12.1.1 的公式给出。记号 $\sim$ 表示「形式上的 Fourier 展开」,收敛性需另证。
Riemann–Lebesgue 引理(陈述):$f$ 可积 ⇒ $a_n,b_n\to 0$($n\to\infty$)。直观上,高频振荡项 $\cos nx,\sin nx$ 与 $f$ 的乘积积分相互抵消。该引理也是 Dirichlet 证明中估计余项的关键工具。
因 $f$ 周期,积分区间可平移,例如 $$a_n=\frac{1}{\pi}\int_{-\pi}^{\pi} f(x)\cos nx\,dx=\frac{1}{\pi}\int_{0}^{2\pi} f(x)\cos nx\,dx.$$
一般周期 $2L$:设 $f(x+2L)=f(x)$,令 $t=\dfrac{\pi x}{L}$,则 $F(t)=f\big(\tfrac{Lt}{\pi}\big)$ 以 $2\pi$ 为周期。换元得 $$f(x)\sim\frac{a_0}{2}+\sum_{n=1}^{\infty}\Big(a_n\cos\frac{n\pi x}{L}+b_n\sin\frac{n\pi x}{L}\Big),$$ 其中 $$a_n=\frac{1}{L}\int_{-L}^{L} f(x)\cos\frac{n\pi x}{L}\,dx,\qquad b_n=\frac{1}{L}\int_{-L}^{L} f(x)\sin\frac{n\pi x}{L}\,dx.$$
换元验证:令 $x=Lt/\pi$,则 $f(x)=F(t)$ 以 $2\pi$ 为周期,标准公式给出 $a_n=\dfrac{1}{\pi}\int_{-\pi}^{\pi}F(t)\cos nt\,dt$,换回 $x$ 即得上式。周期为 $T$ 时,角频率 $\omega_n=2\pi n/T$,公式中 $\pi/L=\pi/(T/2)=2\pi/T$ 与物理频率对应。
例 12.1 求 $f(x)=x$($-\pi $f$ 为奇函数,故 $a_n=0$。$b_n=\dfrac{1}{\pi}\int_{-\pi}^{\pi} x\sin nx\,dx$。分部积分:
$$\int_{-\pi}^{\pi} x\sin nx\,dx=\Big[-\frac{x\cos nx}{n}\Big]_{-\pi}^{\pi}+\frac{1}{n}\int_{-\pi}^{\pi}\cos nx\,dx
=\frac{(-1)^{n+1}2\pi}{n}.$$
故 $b_n=\dfrac{2(-1)^{n+1}}{n}$,级数为
$$x\sim\sum_{n=1}^{\infty}\frac{2(-1)^{n+1}}{n}\sin nx,\quad x\in(-\pi,\pi).$$
在 $x=\pm\pi$ 及整数倍处,级数收敛到 $0$(跳跃间断点的中点),见 §12.3。 例 12.2 矩形波:$f(x)=1$($0 例 12.2' 方波 $f(x)=x$($-\pi\le x\le\pi$),周期 $2\pi$。$f$ 奇,$b_n=\dfrac{2(-1)^{n+1}}{n}$,即
$$x\sim 2\sum_{n=1}^{\infty}\frac{(-1)^{n+1}}{n}\sin nx.$$
与例 12.1 差一因子 $2$ 是因为区间端点约定不同;在 $(-\pi,\pi)$ 内结论一致。 函数 $f$ 仅定义在 $[0,L]$(或 $[-L,L]$)上时,没有天然的周期;要写成三角级数,须先延拓成周期函数,再取 Fourier 级数。延拓方式决定级数形式。 偶延拓(余弦级数):定义
$$F_{\mathrm{even}}(x)=\begin{cases} f(x), & 0\le x\le L,\\ f(-x), & -L\le x<0, \end{cases}$$
再以 $2L$ 为周期延拓。$F_{\mathrm{even}}$ 为偶函数,$b_n=0$,得余弦级数
$$f(x)\sim\frac{a_0}{2}+\sum_{n=1}^{\infty}a_n\cos\frac{n\pi x}{L},\quad x\in[0,L],$$
$$a_0=\frac{2}{L}\int_0^L f(x)\,dx,\qquad a_n=\frac{2}{L}\int_0^L f(x)\cos\frac{n\pi x}{L}\,dx.$$ 奇延拓(正弦级数):定义
$$F_{\mathrm{odd}}(x)=\begin{cases} f(x), & 0\le x\le L,\\ -f(-x), & -L\le x<0, \end{cases}$$
周期 $2L$ 延拓。得正弦级数
$$f(x)\sim\sum_{n=1}^{\infty}b_n\sin\frac{n\pi x}{L},\quad
b_n=\frac{2}{L}\int_0^L f(x)\sin\frac{n\pi x}{L}\,dx.$$ 半区间展开:在 $[0,\pi]$ 上常取 $L=\pi$。同一 $f$ 的余弦展开与正弦展开一般不同——它们对应不同的周期延拓,分别适合 Neumann 型与 Dirichlet 型边界条件的偏微分方程问题。 偶延拓的几何意义:把 $[0,L]$ 上的图像关于 $x=0$ 镜像到 $[-L,0]$,再以 $2L$ 为周期平铺。图像在 $x=0,L,\pm L,\ldots$ 处连续(若 $f$ 在端点连续),适合要求边界斜率为零(绝热边界)的扩散问题。 奇延拓的几何意义:镜像时取负值,图像过原点(若 $f(0)=0$)并在 $x=\pm L$ 处可能出现跳跃。适合 Dirichlet 边界 $u=0$ 的弦振动、热传导问题。 系数公式为何是 $\dfrac{2}{L}\int_0^L$:在 $[0,L]$ 上 $\{\cos\tfrac{n\pi x}{L}\}$ 与 $\{\sin\tfrac{n\pi x}{L}\}$ 关于权 $\tfrac{2}{L}\int_0^L$ 正交归一;偶(奇)延拓后全区间 $[-L,L]$ 上的标准 Fourier 系数化简即得上式。 例 12.3 在 $[0,\pi]$ 上 $f(x)=x(\pi-x)$。 余弦展开:$a_0=\dfrac{2}{\pi}\int_0^{\pi}x(\pi-x)\,dx=\dfrac{2}{\pi}\Big[\dfrac{\pi x^2}{2}-\dfrac{x^3}{3}\Big]_0^{\pi}=\dfrac{\pi^2}{3}$。
$$a_n=\frac{2}{\pi}\int_0^{\pi}x(\pi-x)\cos nx\,dx.$$
分部积分两次(或利用 $x(\pi-x)$ 关于 $x=\pi/2$ 对称),得 $a_n=-\dfrac{2}{n^2}[1+(-1)^n]$:$n$ 偶时 $a_n=0$;$n=2k+1$ 时 $a_{2k+1}=-\dfrac{4}{(2k+1)^2}$。故
$$x(\pi-x)\sim\frac{\pi^2}{6}-\frac{4}{\pi}\sum_{k=0}^{\infty}\frac{\cos(2k+1)x}{(2k+1)^2},\quad x\in[0,\pi].$$ 正弦展开:$b_n=\dfrac{2}{\pi}\int_0^{\pi}x(\pi-x)\sin nx\,dx$。设 $I_n=\int_0^{\pi}x(\pi-x)\sin nx\,dx$,分部积分得 $I_n=\dfrac{2\pi}{n^2}[1-(-1)^n]$,故 $b_n=\dfrac{4\big(1-(-1)^n\big)}{\pi n^3}$:$n$ 偶时 $b_n=0$;$n=2k+1$ 时 $b_{2k+1}=\dfrac{8}{\pi(2k+1)^3}$。 例 12.3' $f(x)=x^2$($-\pi\le x\le\pi$),周期 $2\pi$。$f$ 偶,$b_n=0$。$a_0=\dfrac{2}{\pi}\int_0^{\pi}x^2\,dx=\dfrac{2\pi^2}{3}$;$a_n=\dfrac{2}{\pi}\int_0^{\pi}x^2\cos nx\,dx=\dfrac{4(-1)^n}{n^2}$。故
$$x^2\sim\frac{\pi^2}{3}+4\sum_{n=1}^{\infty}\frac{(-1)^n}{n^2}\cos nx.$$
取 $x=\pi$ 得 $\pi^2\sim\dfrac{\pi^2}{3}+4\sum\dfrac{1}{n^2}$,即 $\displaystyle\sum_{n=1}^{\infty}\frac{1}{n^2}=\frac{\pi^2}{6}$(Euler 问题;亦可用 Parseval 从 $x$ 的展开得到)。 [要点] Euler 公式 $e^{inx}=\cos nx+i\sin nx$ 把三角系化为指数系 $\{e^{inx}\}_{n\in\mathbb{Z}}$。在 $[-\pi,\pi]$ 上
$$\int_{-\pi}^{\pi} e^{inx}\,e^{-imx}\,dx=2\pi\delta_{nm}.$$
设
$$f(x)\sim\sum_{n=-\infty}^{\infty}c_n e^{inx},\qquad
c_n=\frac{1}{2\pi}\int_{-\pi}^{\pi} f(x)e^{-inx}\,dx.$$
与实系数的关系:
$$c_0=\frac{a_0}{2},\quad c_n=\frac{a_n-ib_n}{2}\ (n>0),\quad c_{-n}=\overline{c_n}=\frac{a_n+ib_n}{2}.$$
若 $f$ 为实值函数,则 $c_{-n}=\overline{c_n}$(Hermitian 对称)。 一般周期 $2L$:
$$f(x)\sim\sum_{n=-\infty}^{\infty}c_n e^{in\pi x/L},\qquad
c_n=\frac{1}{2L}\int_{-L}^{L} f(x)e^{-in\pi x/L}\,dx.$$ 实部与虚部:$f$ 实值时 $c_{-n}=\overline{c_n}$;$a_n=c_n+c_{-n}=2\operatorname{Re}c_n$,$b_n=i(c_{-n}-c_n)=-2\operatorname{Im}c_n$($n>0$)。计算时可先求 $c_n$,再还原 $a_n,b_n$,往往比分部积分更省事。 例 12.4 $f(x)=e^x$($-\pi 复形式计算 $f(x)=x$($-\pi,\pi$):$c_n=\dfrac{1}{2\pi}\int_{-\pi}^{\pi}x e^{-inx}\,dx$。$n=0$ 时 $c_0=0$;$n\neq 0$ 时分部积分得 $c_n=\dfrac{i(-1)^n}{n}$,故
$$x\sim\sum_{n\neq 0}\frac{i(-1)^n}{n}e^{inx}=2\sum_{n=1}^{\infty}\frac{(-1)^{n+1}}{n}\sin nx,$$
与例 12.1 一致。 为何引入复形式 复指数 $e^{inx}$ 在求导、卷积、偏微分方程(如热方程分离变量)中运算更简洁;物理上 $|c_n|^2$ 直接给出第 $n$ 个频率分量的能量(Parseval,§12.2)。对称性 $c_{-n}=\overline{c_n}$ 保证实函数的级数仍为实值。 设 $f,g\in L^2[-\pi,\pi]$(平方可积),内积
$$\langle f,g\rangle=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)g(x)\,dx,\qquad
\|f\|_{L^2}=\sqrt{\langle f,f\rangle}=\sqrt{\frac{1}{\pi}\int_{-\pi}^{\pi}|f(x)|^2\,dx}.$$
这是 Hilbert 空间 $L^2[-\pi,\pi]$ 上的标准内积;三角函数系在该内积下正交。 Fourier 部分和
$$S_N(f)(x)=\frac{a_0}{2}+\sum_{n=1}^{N}\big(a_n\cos nx+b_n\sin nx\big)$$
的系数 $a_n,b_n$ 恰为 $f$ 在 $\cos nx,\sin nx$ 方向上的「坐标」。记 $e_0=1/\sqrt2$,$e_n=\cos nx$($n\ge1$),$\tilde{e}_n=\sin nx$,则
$$S_N(f)=\sum_{n=0}^{N}\langle f,e_n\rangle e_n+\sum_{n=1}^{N}\langle f,\tilde{e}_n\rangle\tilde{e}_n$$
(在适当归一化下即正交投影)。 定理(最佳逼近):在 $L^2$ 范数下,$S_N(f)$ 使 $\|f-S_N(f)\|_{L^2}$ 达到最小;即对任意次数 $\le N$ 的三角多项式 $T_N$,有
$$\|f-S_N(f)\|_{L^2}\le\|f-T_N\|_{L^2}.$$
证明要点:$f-S_N(f)\perp T_N-S_N(f)$(正交补),由勾股定理 $\|f-T_N\|^2=\|f-S_N\|^2+\|T_N-S_N\|^2\ge\|f-S_N\|^2$。 记 $f_n=\cos nx$($n\ge0$,$f_0=1/\sqrt2$ 归一化亦可)、$g_n=\sin nx$($n\ge1$)。正交归一系 $\{e_n\}$ 下,部分和系数满足
$$\sum_{n=0}^{N}|c_n|^2\le\|f\|_{L^2}^2\quad\text{(Bessel 不等式)},$$
其中 $c_n$ 为 Fourier 系数(复形式或实形式等价)。等号对一切 $N$ 成立当且仅当 $f\in\mathrm{span}\{e_n\}_{n\le N}$;对无穷级数,等号成立即 Parseval。当 $N\to\infty$ 时,若 $f\in L^2$,则
$$\|f-S_N(f)\|_{L^2}\to 0\quad\text{(平方平均收敛)},$$
且 Parseval 等式
$$\|f\|_{L^2}^2=\frac{a_0^2}{2}+\sum_{n=1}^{\infty}(a_n^2+b_n^2)
=\sum_{n=-\infty}^{\infty}|c_n|^2.$$
实形式下 $\tfrac{a_0^2}{2}+\sum(a_n^2+b_n^2)$ 表示信号总能量在频率上的分配。 证明 Bessel 不等式(概要):设 $r_N=f-S_N(f)$。$S_N(f)\perp r_N$(正交投影),故
$$\|f\|_{L^2}^2=\|S_N(f)\|_{L^2}^2+\|r_N\|_{L^2}^2\ge\|S_N(f)\|_{L^2}^2=\sum_{|n|\le N}|c_n|^2.$$
令 $N\to\infty$ 即得 Bessel;平方平均收敛成立时范数极限给出 Parseval。 实系数形式的 Parseval:由 $a_n=c_n+\overline{c_{-n}}$、$b_n=i(c_{-n}-c_n)$($n>0$)及 $c_{-n}=\overline{c_n}$ 可验证
$$\|f\|_{L^2}^2=\frac{1}{\pi}\int_{-\pi}^{\pi}|f|^2=\frac{a_0^2}{2}+\sum_{n=1}^{\infty}(a_n^2+b_n^2).$$ 三种收敛的比较: 一致 $\Rightarrow$ 逐点(在适当条件下)$\Rightarrow$ 平方平均;反向一般不成立。 为何 $L^2$ 收敛最重要 工程与概率中更关心能量 $\int|f|^2$ 而非逐点值;$L^2$ 收敛保证截断级数的能量误差趋于零。数值计算中 $N$ 项部分和在 $L^2$ 意义下已是「最优」$N$ 阶三角逼近,Parseval 给出可检验的误差上界 $\|r_N\|_{L^2}^2=\|f\|_{L^2}^2-\sum_{|n|\le N}|c_n|^2$。 例 12.5 对例 12.1 的 $f(x)=x$($-\pi,\pi$),Parseval 验证
$$\frac{1}{\pi}\int_{-\pi}^{\pi}x^2\,dx=\frac{2\pi^2}{3}=\sum_{n=1}^{\infty}\frac{4}{n^2}=4\cdot\frac{\pi^2}{6}=\frac{2\pi^2}{3}.$$ 例 12.5' 对 $f(x)=x^2$(例 12.3'),$\|f\|_{L^2}^2=\dfrac{1}{\pi}\int_{-\pi}^{\pi}x^4\,dx=\dfrac{2\pi^4}{5}$。又 $\dfrac{a_0^2}{2}+\sum a_n^2=\dfrac{2\pi^4}{9}+16\sum\dfrac{1}{n^4}$,得 $\displaystyle\sum_{n=1}^{\infty}\frac{1}{n^4}=\frac{\pi^4}{90}$。 [要点] Dirichlet 条件:$f$ 在 $[-\pi,\pi]$ 上分段光滑(分段 $C^1$,有限个第一类间断点),作 $2\pi$ 周期延拓。 Dirichlet 定理:在任意点 $x_0$,Fourier 级数 $S_f(x_0)$ 收敛到
$$\frac{f(x_0+)+f(x_0-)}{2}.$$
特别地,若 $f$ 在 $x_0$ 连续,则 $S_f(x_0)=f(x_0)$;在跳跃间断处收敛到左右极限的平均(Gibbs 现象:间断点附近部分和有过冲,但不影响积分意义下的收敛)。 证明思路(概要): 一致收敛的补充:若 $f$ 连续且以 $2\pi$ 为周期的导数 $f'$ 分段连续,则 Fourier 级数在 $\mathbb{R}$ 上一致收敛于 $f$(M 判别法:$|a_n|,|b_n|=O(1/n)$,$\sum 1/n^2$ 收敛 ⇒ 系数可和)。 Gibbs 现象 在跳跃间断附近,$S_N(f)$ 出现过冲(约 9% 跃变量),随 $N$ 增大过冲区变窄但不消失。这不与 Dirichlet 定理矛盾——定理说的是逐点极限在间断点等于平均值,而非一致收敛。 定理:若 $f\in L^2[-\pi,\pi]$,则 $\|f-S_N(f)\|_{L^2}\to 0$。 证明思路: Riesz–Fischer 定理(陈述):对任意满足 $\sum(|a_n|^2+|b_n|^2)<\infty$ 的系数序列,存在 $f\in L^2[-\pi,\pi]$ 以之为 Fourier 系数。即 Fourier 映射是 $L^2$ 到 $\ell^2$ 的等距同构(在 Parseval 意义下)。 对比:Dirichlet 定理要分段光滑(点态);$L^2$ 收敛只需平方可积,但不要求逐点极限等于 $f$。存在连续函数其 Fourier 级数在某点发散(Du Bois-Reymond 型反例,超出本课程范围);但 $L^2$ 理论保证「几乎处处」意义下 Cesàro 平均收敛(Fejér 定理)。 Fejér 定理(补充):对连续周期 $f$,Fejér 和 $\sigma_N(f)=\dfrac{1}{N+1}\sum_{k=0}^{N}S_k(f)$ 一致收敛于 $f$。说明三角级数在 Cesàro 意义下表现良好,而 Gibbs 现象只影响 $S_N$ 而非 $\sigma_N$。 例 12.6 $f(x)=\begin{cases}1,&0 例 12.6' $f(x)=\begin{cases}x,&0\le x<\pi,\\x-2\pi,&-\pi\le x<0,\end{cases}$ 周期 $2\pi$(锯齿波)。$f$ 连续但 $f'(\pi-)=1,\ f'(\pi+)=1$ 在 $x=\pi$ 处导数跳跃。Fourier 系数 $b_n=\dfrac{2(-1)^n}{n}$,级数在 $\mathbb{R}$ 上收敛到 $f$ 的周期延拓;在 $x=\pi+2k\pi$ 处 Dirichlet 值仍为 $\pi$(连续点)。 当 $f$ 定义在全直线 $\mathbb{R}$ 上且非周期时,可视为周期 $2L$ 的极限 $L\to\infty$。记 $\omega_n=\dfrac{n\pi}{L}$,$\Delta\omega=\dfrac{\pi}{L}$,复系数
$$c_n=\frac{1}{2L}\int_{-L}^{L}f(x)e^{-i\omega_n x}\,dx\approx\frac{\Delta\omega}{2\pi}\int_{-\infty}^{\infty}f(x)e^{-i\omega x}\,dx.$$
形式极限给出 Fourier 变换
$$\hat{f}(\omega)=\mathcal{F}[f](\omega)=\int_{-\infty}^{\infty}f(x)e^{-i\omega x}\,dx,$$
逆变换
$$f(x)=\mathcal{F}^{-1}[\hat{f}](x)=\frac{1}{2\pi}\int_{-\infty}^{\infty}\hat{f}(\omega)e^{i\omega x}\,d\omega.$$
(常数 $\dfrac{1}{2\pi}$ 的放置依约定不同;本书采用上式。) Fourier 积分公式:在适当条件下(如 $f$ 绝对可积且分段光滑),
$$f(x)=\frac{1}{2\pi}\int_{-\infty}^{\infty}\Big[\int_{-\infty}^{\infty}f(t)e^{-i\omega t}\,dt\Big]e^{i\omega x}\,d\omega.$$ 形式推导:周期 $2L$ 时 $f(x)=\sum_{n=-\infty}^{\infty}c_n e^{in\pi x/L}$,$c_n=\dfrac{1}{2L}\int_{-L}^{L}f(t)e^{-in\pi t/L}\,dt$。令 $\omega_n=n\pi/L$,$\Delta\omega=\pi/L$,则
$$f(x)\approx\sum_n c_n e^{i\omega_n x}=\sum_n\Big[\frac{\Delta\omega}{2\pi}\int_{-L}^{L}f(t)e^{-i\omega_n t}\,dt\Big]e^{i\omega_n x}\to\frac{1}{2\pi}\int_{-\infty}^{\infty}\hat{f}(\omega)e^{i\omega x}\,d\omega.$$ 收敛性:若 $f\in L^1(\mathbb{R})$ 且在 $x$ 处满足 Dirichlet 型条件,则上述反演公式在 $x$ 成立;若 $f\in L^2(\mathbb{R})$,则
$$\lim_{R\to\infty}\frac{1}{2\pi}\int_{-R}^{R}\hat{f}(\omega)e^{i\omega x}\,d\omega=f(x)\quad\text{($L^2$ 意义)}.$$ 设 $\mathcal{F}[f]=\hat{f}$,在函数足够正则时: 卷积定理的证明梗概:$\mathcal{F}[f*g](\omega)=\int\int f(x-t)g(t)e^{-i\omega x}\,dt\,dx$,换元 $u=x-t$ 得 $\int f(u)e^{-i\omega u}\,du\cdot\int g(t)e^{-i\omega t}\,dt=\hat{f}(\omega)\hat{g}(\omega)$。 Plancherel 定理:对 $f\in L^2(\mathbb{R})$,$\mathcal{F}$ 可延拓为 $L^2$ 上的等距同构(差一个 $2\pi$ 因子),即 Parseval 等式对 $L^2$ 函数成立;这是 Fourier 变换在 Hilbert 空间框架下的严格表述。 Parseval(连续情形):
$$\int_{-\infty}^{\infty}|f(x)|^2\,dx=\frac{1}{2\pi}\int_{-\infty}^{\infty}|\hat{f}(\omega)|^2\,d\omega.$$ 例 12.7 矩形脉冲 $f(x)=\begin{cases}1,&|x|a,\end{cases}$
$$\hat{f}(\omega)=\int_{-a}^{a}e^{-i\omega x}\,dx=\frac{2\sin a\omega}{\omega}\quad(\omega\neq 0),\qquad \hat{f}(0)=2a.$$
$\hat{f}$ 为 $\operatorname{sinc}$ 型;$\omega\to\infty$ 时振荡衰减,与 Riemann–Lebesgue 引理一致。物理上,窄脉冲(小 $a$)⇒ 频谱展宽。 例 12.8 高斯函数 $f(x)=e^{-x^2/2}$。利用 $\int_{-\infty}^{\infty}e^{-ax^2+bx}\,dx=\sqrt{\pi/a}\,e^{b^2/(4a)}$ 得
$$\hat{f}(\omega)=\sqrt{2\pi}\,e^{-\omega^2/2},$$
即高斯的 Fourier 变换仍是高斯(相差常数与尺度)。 例 12.9 用微分性质求 $\mathcal{F}[e^{-|x|}]$。令 $g(x)=e^{-|x|}$,则 $-g''+g=2\delta(x)$(分布意义),变换得 $(\omega^2+1)\hat{g}(\omega)=2$,故 $\hat{g}(\omega)=\dfrac{2}{\omega^2+1}$(Lorentz 型谱)。 例 12.10 求 $\mathcal{F}[\mathbf{1}_{[-a,a]}(x)]$ 与 $\mathcal{F}[\sin\omega_0 x]$。前者即例 12.7;后者 $\mathcal{F}[\sin\omega_0 x]=\dfrac{i}{2}\big[\delta(\omega+\omega_0)-\delta(\omega-\omega_0)\big]$(分布意义),体现单频信号的线谱。 偏微分方程:一维热方程 $u_t=ku_{xx}$ 在 $[-\pi,\pi]$ 上分离变量 $u=T(t)X(x)$ 得 $X''+\lambda X=0$,本征函数为 $\sin nx,\cos nx$,解展开为 Fourier 级数;无穷长 rod 上用 Fourier 变换把 PDE 化为 $\hat{u}_t=-k\omega^2\hat{u}$,形式解 $\hat{u}(\omega,t)=\hat{u}_0(\omega)e^{-k\omega^2 t}$。 信号处理:$|\hat{f}(\omega)|$ 为频谱幅度;卷积定理使滤波(乘窗)对应时域卷积,是数字信号处理的基础。 Heisenberg 不确定性(直观):在 $L^2$ 意义下,$f$ 与其 Fourier 变换不能同时高度局域——$\|xf\|\cdot\|\omega\hat{f}\|$ 有下界(量子力学中位置与动量的对偶关系与此同构)。严格表述需分布理论,此处仅作背景联系。 $2\pi$ 约定的提醒 不同教材在 $\mathcal{F}$、$\mathcal{F}^{-1}$ 前放置 $2\pi$ 的方式不同(有的在正变换、有的在逆变换、有的各放 $\sqrt{2\pi}$)。本书采用 $\hat{f}(\omega)=\int f e^{-i\omega x}\,dx$,$f(x)=\dfrac{1}{2\pi}\int\hat{f}(\omega)e^{i\omega x}\,d\omega$。换书查表时务必核对 Parseval 式中的系数。 [要点] 12.1 求 $f(x)=|x|$($-\pi 12.2 将 $f(x)=1$($0 12.3 设 $f(x)=x^2$($-\pi\le x\le\pi$),周期延拓。求 Fourier 级数,并用 Parseval 等式计算 $\displaystyle\sum_{n=1}^{\infty}\frac{1}{n^2}$。 12.4 证明:若 $f\in L^2[-\pi,\pi]$ 且 $a_n=b_n=0$ 对所有 $n\ge1$,则 $f(x)=\dfrac{a_0}{2}$ 几乎处处成立(提示:$S_N(f)\to f$ 于 $L^2$,部分和为零)。 12.5 求 $\mathcal{F}[f]$,其中 $f(x)=e^{-\alpha|x|}$($\alpha>0$)。验证 Parseval 等式 $\displaystyle\int_{-\infty}^{\infty}|f|^2=\frac{1}{2\pi}\int_{-\infty}^{\infty}|\hat{f}|^2$。 提示概要:12.1 注意 $|x|$ 为偶函数;12.2 正弦展开在 $(0,\pi)$ 内和为 $1$,余弦展开和为 $1$ 但延拓不同;12.3 见例 12.3';12.4 用 $L^2$ 完备性与正交系;12.5 化为 $\dfrac{2\alpha}{\alpha^2+\omega^2}$,左侧积分为 $\dfrac{2}{\alpha}$,右侧用 $\int_{-\infty}^{\infty}(\alpha^2+\omega^2)^{-2}\,d\omega=\dfrac{\pi}{\alpha^3}$。 (详解可在后续修订中补入。)12.1.3 有限区间上的展开与延拓
12.1.4 复数形式的 Fourier 级数
12.2 平方平均收敛
12.2.1 $L^2$ 内积与正交投影
12.2.2 Bessel 不等式与 Parseval 等式
收敛类型
定义
典型条件
一致收敛
$\sup_x\|S_N-f\|\to 0$
$f$ 分段 $C^1$ 且端点等化(较强)
逐点收敛
$S_N(x_0)\to f(x_0)$
Dirichlet 条件
平方平均
$\|S_N-f\|_{L^2}\to 0$
$f\in L^2$(最宽)
12.3 收敛性定理
12.3.1 Dirichlet 定理(陈述与证明思路)
12.3.2 平方平均收敛的证明思路
12.4 Fourier 积分与 Fourier 变换
12.4.1 从周期到非周期
12.4.2 基本性质与常用变换
性质
公式
线性
$\mathcal{F}[\alpha f+\beta g]=\alpha\hat{f}+\beta\hat{g}$
平移
$\mathcal{F}[f(x-x_0)](\omega)=e^{-i\omega x_0}\hat{f}(\omega)$
调制
$\mathcal{F}[e^{i\omega_0 x}f(x)](\omega)=\hat{f}(\omega-\omega_0)$
微分
$\mathcal{F}[f'](\omega)=i\omega\hat{f}(\omega)$
卷积
$\mathcal{F}[f*g]=\hat{f}\cdot\hat{g}$,$(f*g)(x)=\int f(x-t)g(t)\,dt$
12.4.3 应用简述
本章综合习题(选)