M09 · 第 3 章 · 第二编 傅里叶变换

傅里叶变换与卷积

固定连续傅里叶变换的归一化约定,说明逆变换的适用条件,推导平移、缩放、微分和卷积性质,并用双边指数函数与矩形脉冲完成可复算的时域—频域对应。

报告页面错误
预备知识傅里叶级数、收敛与 Gibbs 现象反常积分复数与复平面

本章目标

  1. 在一套固定的 $2\pi$ 归一化约定下计算变换并写出逆变换。
  2. 区分绝对可积、平方可积与逐点反演条件,不把形式积分自动当作点态等式。
  3. 由变量代换推导平移、调制和缩放性质,并正确保留相位与绝对值因子。
  4. 在可用 Fubini 定理的条件下证明卷积定理,并说明乘积对应的频域卷积常数。
  5. 直接复算双边指数、矩形脉冲和三角函数形卷积的变换对。
页面阅读位置0% · 仅保存在此浏览器
章节未开始
本册完成进度0/6 章 · 0%
本页目录

连续频谱需要先固定约定

傅里叶级数用整数频率表示周期函数。对非周期函数,频率不再局限于离散整数,求和相应变成对连续频率的积分。不同教材会把 2π2\pi 分配在正变换、逆变换或指数中;每种约定都可以成立,但同一推导中不能混用。

本章始终采用角频率 ξ\xi

傅里叶变换与逆变换约定

fL1(R)f\in L^1(\mathbb R),定义

f^(ξ)=F[f](ξ)=f(x)eiξxdx.\widehat f(\xi)=\mathcal F[f]\,(\xi) =\int_{-\infty}^{\infty}f(x)e^{-i\xi x}\,\mathrm dx.

在逆变换成立的意义下,

f(x)=F1[f^](x)=12πf^(ξ)eiξxdξ.f(x)=\mathcal F^{-1}[\widehat f]\,(x) =\frac1{2\pi}\int_{-\infty}^{\infty} \widehat f(\xi)e^{i\xi x}\,\mathrm d\xi.

指数中的负号和逆变换前的 1/(2π)1/(2\pi) 是一套整体约定。若改用普通频率 ν\nu 与指数 e2πiνxe^{-2\pi i\nu x},缩放、卷积和微分公式中的 2π2\pi 位置会一起改变。

fL1f\in L^1 时,f^\widehat f 有界、连续,并在 ξ|\xi|\to\infty 时趋于零;这并不保证 f^L1\widehat f\in L^1。一个简洁的充分条件是 ff 属于 Schwartz 空间:它及所有阶导数都比任意多项式倒数衰减得更快,此时正反变换都可逐点进行。若 f,f^L1f,\widehat f\in L^1,逆积分定义一个连续函数,并与 ff 几乎处处相等;若 ff 取其连续代表,等式处处成立。

对分段光滑且满足适当衰减的函数,截断逆积分在跳跃处通常趋于左右极限的平均,而不是任意指定的单点值。对 L2L^2 函数,Plancherel 定理把变换延拓为均方意义的算子,并给出

f(x)2dx=12πf^(ξ)2dξ.\int_{-\infty}^{\infty}|f(x)|^2\,\mathrm dx =\frac1{2\pi}\int_{-\infty}^{\infty}|\widehat f(\xi)|^2\,\mathrm d\xi.

这里的相等是 L2L^2 意义,不自动提供每一点的反演值。

可直接复算的双边指数变换

例 1:双边指数衰减的有理频谱

a>0a>0f(x)=eaxf(x)=e^{-a|x|}。把积分按零点分开,并利用偶性:

f^(ξ)=20eaxcos(ξx)dx.\widehat f(\xi) =2\int_0^\infty e^{-ax}\cos(\xi x)\,\mathrm dx.

0e(aiξ)xdx=1aiξ\int_0^\infty e^{-(a-i\xi)x}\,\mathrm dx =\frac1{a-i\xi}

取实部,得到

f^(ξ)=2aa2+ξ2.\widehat f(\xi)=\frac{2a}{a^2+\xi^2}.

两端都绝对可积,所以可用逆变换核对:

eax=12π2aa2+ξ2eiξxdξ.e^{-a|x|} =\frac1{2\pi}\int_{-\infty}^{\infty} \frac{2a}{a^2+\xi^2}e^{i\xi x}\,\mathrm d\xi.

x=0x=0 时,这等价于 (a2+ξ2)1dξ=π/a\int_{-\infty}^{\infty}(a^2+\xi^2)^{-1}\mathrm d\xi=\pi/a。频谱宽度随 aa 增大而变宽:时域衰减更快、形状更集中,对应更广的频率范围。

平移、调制与缩放

这些规则都可从定义中的变量代换得到。令 g(x)=f(xx0)g(x)=f(x-x_0),取 u=xx0u=x-x_0,则

g^(ξ)=eiξx0f^(ξ).\widehat g(\xi)=e^{-i\xi x_0}\widehat f(\xi).

时域平移只增加线性相位,不改变频谱模长。反过来,乘以复指数会平移频谱:

F[eiξ0xf(x)](ξ)=f^(ξξ0).\mathcal F[e^{i\xi_0x}f(x)]\,(\xi) =\widehat f(\xi-\xi_0).

g(x)=f(ax)g(x)=f(ax)a0a\ne0,变量代换必须处理积分方向,得到

g^(ξ)=1af^ ⁣(ξa).\widehat g(\xi)=\frac1{|a|}\widehat f\!\left(\frac\xi a\right).

因子是 1/a1/|a|,不是 1/a1/a。当 a>1|a|>1 时,函数在时域被压窄,频谱横向展开并降低幅值,以保持积分尺度一致。

ff 绝对连续,且 f,fL1f,f'\in L^1 并满足无穷远边界项消失,分部积分给出

F[f](ξ)=iξf^(ξ).\mathcal F[f']\,(\xi)=i\xi\widehat f(\xi).

同样,在 xf(x)L1xf(x)\in L^1 等条件下,

ddξf^(ξ)=iF[xf(x)](ξ).\frac{\mathrm d}{\mathrm d\xi}\widehat f(\xi) =-i\mathcal F[xf(x)]\,(\xi).

微分强化高频、乘以位置对应频谱求导;若边界项不消失或导数含跳跃产生的广义函数,上述经典公式必须改用分布框架解释。

卷积把邻域贡献叠加起来

卷积

对可积函数 f,gf,g,定义

(fg)(x)=f(t)g(xt)dt.(f*g)(x)=\int_{-\infty}^{\infty}f(t)g(x-t)\,\mathrm dt.

f,gL1(R)f,g\in L^1(\mathbb R) 时,卷积几乎处处有定义,且 fg1f1g1\lVert f*g\rVert_1\le\lVert f\rVert_1\lVert g\rVert_1

g(xt)g(x-t) 表示把核翻转后移动到 xx,再与 f(t)f(t) 逐点相乘并积分。若 gg 是归一化、集中在零附近的非负核,卷积可视为局部加权平均;若核含正负振荡,它也可检测边缘或选取频段,不能一概理解为“平滑”。

f,gL1f,g\in L^1,被积函数的绝对积分有限,可以用 Fubini 定理交换积分次序:

F[fg](ξ)=RRf(t)g(xt)eiξxdtdx=Rf(t)eiξtdtRg(u)eiξudu=f^(ξ)g^(ξ).\begin{aligned} \mathcal F[f*g]\,(\xi) &=\int_{\mathbb R}\int_{\mathbb R} f(t)g(x-t)e^{-i\xi x}\,\mathrm dt\,\mathrm dx\\ &=\int_{\mathbb R}f(t)e^{-i\xi t}\,\mathrm dt \int_{\mathbb R}g(u)e^{-i\xi u}\,\mathrm du\\ &=\widehat f(\xi)\widehat g(\xi). \end{aligned}

这就是卷积定理。与本章约定配套,时域乘积对应

F[fg](ξ)=12π(f^g^)(ξ),\mathcal F[fg]\,(\xi) =\frac1{2\pi}(\widehat f*\widehat g)(\xi),

其中的 1/(2π)1/(2\pi) 不可遗漏。

例 2:两个矩形脉冲卷积成三角函数

r(x)=1[1/2,1/2](x).r(x)=\mathbf1_{[-1/2,\,1/2]}(x).

它的变换可直接积分:

r^(ξ)=1/21/2eiξxdx=2sin(ξ/2)ξ,\widehat r(\xi) =\int_{-1/2}^{1/2}e^{-i\xi x}\,\mathrm dx =\frac{2\sin(\xi/2)}{\xi},

ξ=0\xi=0 处用极限补成一。卷积 (rr)(x)(r*r)(x) 等于区间 [1/2,1/2][-1/2,1/2][x1/2,x+1/2][x-1/2,x+1/2] 的重叠长度,因此

(rr)(x)={1x,x1,0,x>1.(r*r)(x)= \begin{cases} 1-|x|,& |x|\le1,\\ 0,& |x|>1. \end{cases}

卷积定理立即给出三角函数的频谱

F[rr](ξ)=(2sin(ξ/2)ξ)2.\mathcal F[r*r]\,(\xi) =\left(\frac{2\sin(\xi/2)}{\xi}\right)^2.

可在 ξ=0\xi=0 核对两边:左侧等于三角函数面积一,右侧极限也为一。这个零频检查常能发现归一化或尺度错误。

卷积定理怎样求解线性方程

常系数线性微分算子在频域中变成关于 ξ\xi 的乘子。例如形式上对

u+a2u=f,-u''+a^2u=f,

变换后得到

(ξ2+a2)u^(ξ)=f^(ξ),u^(ξ)=f^(ξ)ξ2+a2.(\xi^2+a^2)\widehat u(\xi)=\widehat f(\xi), \qquad \widehat u(\xi)=\frac{\widehat f(\xi)}{\xi^2+a^2}.

GG 的变换为 (ξ2+a2)1(\xi^2+a^2)^{-1},则 u=Gfu=G*f。这把“求逆微分算子”转化为“与 Green 核卷积”。但分母若在实频率上为零,就涉及可解性、边界条件或分布意义,不能简单逐点相除。

在信号处理中,线性平移不变系统的输出常写成 y=hxy=h*x,频域关系为 y^=h^x^\widehat y=\widehat h\,\widehat xh^\widehat h 是系统对不同频率的增益与相位响应。由有限采样计算离散卷积时,还需区分线性卷积、循环卷积和边界填充;连续定理本身不替代离散实现的索引检查。

常见误区

写出逆积分就已经证明逐点反演

变换存在不代表逆积分绝对收敛。必须说明使用 Schwartz、L1L^1 加变换可积、分段光滑的截断极限,还是 L2L^2 均方反演。不同条件给出的等式意义不同。

时域平移会把频谱图整体平移

时域平移乘上线性相位,频谱模长不动;时域调制才使频谱平移。混淆二者会错误判断滤波器的通带位置。

函数压缩后的频谱只需替换自变量

f(ax)f(ax) 的变换除了把 ξ\xi 替换为 ξ/a\xi/a,还有 1/a1/|a|。遗漏该因子会使零频值、积分和能量尺度不一致。

频谱可积不是所有可积函数自动拥有的性质

矩形脉冲 rL1r\in L^1,但其频谱 2sin(ξ/2)/ξ2\sin(\xi/2)/\xi 的绝对积分发散。因此不能仅凭 rL1r\in L^1 直接套用“fff^\widehat f 都可积”的绝对收敛反演条件;矩形脉冲的跳跃点需用截断反演或 L2L^2 理论解释。

练习

练习

f(x)=1[2,2](x)f(x)=\mathbf1_{[-2,2]}(x)。按本章约定求 f^(ξ)\widehat f(\xi),并写出 ξ=0\xi=0 处的连续补值。

查看解答

直接积分得 f^(ξ)=22eiξxdx=2sin(2ξ)/ξ\widehat f(\xi)=\int_{-2}^{2}e^{-i\xi x}\mathrm dx=2\sin(2\xi)/\xi。当 ξ0\xi\to0 时,sin(2ξ)2ξ\sin(2\xi)\sim2\xi,所以补值为四;这也等于 ff 的积分,即区间长度。

练习

已知 f^(ξ)=F(ξ)\widehat f(\xi)=F(\xi),求 g(x)=f(3(x2))g(x)=f(3(x-2)) 的变换。分别指出频谱模长的缩放和相位因子。

查看解答

先令 h(x)=f(3x)h(x)=f(3x),则 h^(ξ)=13F(ξ/3)\widehat h(\xi)=\frac13F(\xi/3);而 g(x)=h(x2)g(x)=h(x-2),所以

g^(ξ)=e2iξ13F(ξ/3).\widehat g(\xi)=e^{-2i\xi}\frac13F(\xi/3).

频谱横向扩大三倍、幅值乘 1/31/3,平移二只增加模长为一的相位因子。

练习

h(x)=ex1[0,)(x)h(x)=e^{-x}\mathbf1_{[0,\infty)}(x)。先求 h^\widehat h,再直接计算 hhh*h,并核对卷积定理。

查看解答

h^(ξ)=0e(1+iξ)xdx=(1+iξ)1\widehat h(\xi)=\int_0^\infty e^{-(1+i\xi)x}\mathrm dx=(1+i\xi)^{-1}。当 x<0x<0 时卷积为零;当 x0x\ge0 时,只有 0tx0\le t\le x 有贡献,故

(hh)(x)=0xete(xt)dt=xex.(h*h)(x)=\int_0^x e^{-t}e^{-(x-t)}\mathrm dt=xe^{-x}.

另一方面,卷积定理给出 F[hh](ξ)=h^(ξ)2=(1+iξ)2\mathcal F[h*h]\,(\xi)=\widehat h(\xi)^2=(1+i\xi)^{-2},与直接对 xex1[0,)xe^{-x}\mathbf1_{[0,\infty)} 积分所得结果一致。

练习

给定频谱 F(ξ)=1[Ω,Ω](ξ)F(\xi)=\mathbf1_{[-\Omega,\Omega]}(\xi),其中 Ω>0\Omega>0。形式上求逆变换,补出 x=0x=0 的值,并说明为何结果不宜仅用 L1L^1 绝对反演来解释。

查看解答

由逆变换

f(x)=12πΩΩeiξxdξ=sin(Ωx)πx,f(x)=\frac1{2\pi}\int_{-\Omega}^{\Omega}e^{i\xi x}\mathrm d\xi =\frac{\sin(\Omega x)}{\pi x},

x=0x=0 处的极限为 Ω/π\Omega/\pi。该函数按 1/x1/|x| 振荡衰减,通常不属于 L1(R)L^1(\mathbb R),但属于 L2(R)L^2(\mathbb R);因此这组变换对应更自然地由 L2L^2 的 Plancherel 理论或适当的广义反演理解。

知识关系

参考资源

课程 · 2024

Topics in Fourier Analysis

Daniel Stroock

用于核对正交展开、Fourier 变换、反演条件和分布理论的对象边界。

打开官方来源

MIT OpenCourseWare 的 Fourier Analysis 专题课程系统覆盖 Fourier 级数、L1/L2L^1/L^2 变换与温和分布,适合核对本章的归一化、反演条件和卷积定理。

课程 · 2011

MIT 18.03SC Differential Equations

Arthur Mattuck, Haynes Miller

为傅里叶方法、模态分解及其微分方程应用提供连续课程背景。

打开官方来源

MIT 微分方程课程展示 Fourier 方法怎样把线性方程和边值问题分解为独立模式,可用于继续学习频域求解。

课程 · 2016

MIT 8.03SC Physics III: Vibrations and Waves

Yen-Jie Lee

连接波动方程的数学解、边界条件、驻波和物理观测。

打开官方来源

MIT 振动与波课程提供频谱、叠加和波动传播的物理背景,适合把变换性质与可观测信号联系起来。

书籍 · 2016

Calculus Volume 2

Gilbert Strang, Edwin Herman

用于核对数列与级数的收敛条件、积分技巧、反常积分和幂级数例题。

打开官方来源

OpenStax Calculus Volume 2 在此仅作为广义积分、分部积分和无穷级数等先修计算的复习来源;高级 Fourier 反演与 L2L^2 理论以专门课程资源为准。

后续学习

下一章将引入广义函数,使 Dirac 冲激、常数函数和纯频率指数也能进入统一的变换规则,并进一步讨论采样与 Poisson 求和。随后在偏微分方程中,卷积核会具体表现为热核、Green 函数或传播算子。