傅里叶变换

来自testwiki
跳转到导航 跳转到搜索

Template:NoteTA

Template:傅里叶变换傅里叶变换Template:Lang-frTemplate:Lang-en,缩写:FT)是一种线性变换,通常定义为一种积分变换。其基本思想是一个函数可以用(可数或不可数,可数的情况对应于傅里叶级数)无穷多个周期函数的线性组合来逼近,从而这些组合系数在保有原函数的几乎全部信息的同时,还直接地反映了该函数的“頻域特征”。

傅里叶变换将函数的时域(红色)与频域(蓝色)相关联。频谱中的不同成分频率在频域中以峰值形式表示。

因其基本思想首先由法国学者约瑟夫·傅里叶系统地提出,所以以其名字来命名以示纪念。在现代数学理论中,傅里叶积分变换可以得到各种推广,并在分析学中有广泛应用,构成了調和分析这一数学领域。

经过傅里叶变换生成的函数 f^ 称作原函数 f 的傅里叶变换,应用意义上称作频谱。在特定情況下,傅里叶变换是可逆的,即将 f^ 通过逆变换可以得到其原函数 f。通常情况下,f 是一个实函数,而 f^ 则是一个复数值函数,其函数值作为复数可同时表示振幅相位

定义

对于不同种类的函数,有一些不同版本的傅里叶变换的定义。对于定义在欧几里得空间 d 上的函数,可给出通常的连续傅里叶变换;对于定义在 d环面 𝕋d 上的函数,或者说周期函数,就给出傅里叶级数。进行傅里叶变换的函数的定义域可以推广到拓扑群,如局部紧交换群

若「傅里叶变换」一词不加任何限定语,则往往是指所谓「连续傅里叶变换」,下文我们默认讨论连续傅里叶变换。其他的常见变种列于傅里叶变换的其他变种一节。

傅里叶积分变换

对于一个多实变量的复函数 f:d ,其傅里叶变换的结果(记作 f^:d(f) )最常見的定义方式是下面的傅里叶积分

(f)(𝐤)=f^(𝐤)=df(𝐱) e2πi𝐤𝐱 d𝐱,

称为傅里叶积分变换,其中粗体的 𝐱,𝐤dd 维实向量,点乘号 则表示欧几里得内积。值得一提的是,不同的作者可能在定义中对积分号前的系数 1 和指数上的系数 2πi 进行调整,不同领域有不同的惯用约定[1],参见变换参数的常见约定

常可定义另一积分变换 1

1(f)(𝐱)=fˇ(𝐱)=df(𝐤) e2πi𝐱𝐤 d𝐤,

它称作傅里叶积分逆变换,满足 1=id

若某函数的傅里叶积分不收敛,则这一版本的傅里叶变换对其就无法定义;即便傅里叶积分收敛,所得的函数的积分逆变换也可能不收敛;或者这两个变换非互逆关系。所以须了解什么样的函数是可变换的,而且满足不同假設的函数的傅里叶变换可能性质不同。

另外,可行上述变换的函数太少,如各种多项式函数都无法用上面的积分定义,这些情况对于信号频域分析等直接应用而言也十分重要,从而有必要进行推广,这就是下面几小节的主题。

收敛性

对于一个 d 上的连续函数 f 而言,要使其黎曼积分收敛,主要需限制其在趋向无穷远时的衰减速度。为使前述傅里叶积分收敛,可以考虑使得 |xd+ϵ||f(x)| 有界的连续函数 f ,其中 ϵ>0 。容易验证这些函数在通常加法和数乘下构成一个向量空间,记作 ϵ 。这些函数都是黎曼可積[2],并且乘上模为1的指数函数所得函数仍在 ϵ 中,也就是说其傅里叶积分也是收敛的。

速降函数

ϵ 上的傅里叶变换已具备许多良好的性质。然而,由于未必有 fMϵf^Mϵ ,所以 :ϵϵ 可能不具可逆性。一些分析表明, f^ 趋向无穷远时的衰减行为与 f 的连续性与可微性有关:为使 f^ 更快地衰减, f 应具有更好的光滑性,反过来也一样[2]。这启发我们研究所谓速降函数,其任意阶导数的衰减速度都快于任意负幂次函数,其构成的向量空间称为速降函数空间,记作 𝒮 。速降函数的傅里叶变换仍是速降函数,那么上面的积分逆变换可以定义在整个 𝒮 上。可以证明该积分确实是逆变换,于是傅里叶积分变换是 𝒮 上的一个自同构

勒贝格积分与勒贝格可积函数

上述积分变换是基于黎曼积分的。然而相比之下,使用基于测度勒貝格積分来定义傅里叶积分变换有许多的优势。下文提到傅里叶积分变换时都默认是勒贝格积分意义上的。

d 上全体勒贝格可积的函数构成的向量空间记作 L1(d) 或此处省略地记为 L1 ,其是Lp空间的一种。对于勒贝格可积函数 fL1 而言,由于 |f(x)e2πi𝐤𝐱|=|f(x)||e2πi𝐤𝐱|=|f(x)|f(x)e2πi𝐤𝐱 也是勒贝格可积的,也就是说傅里叶积分变换在 L1 上有定义。

一般来说可积函数 fL1 的傅里叶变换 f^ 未必是可积的。但对于其中可积的 f^L1 可以证明[3],积分逆变换

df^(𝐤) e2πi𝐱𝐤 d𝐤

对于 𝐱d 幾乎處處收敛于 f(𝐱)

平方可积函数

对于勒贝格积分定义的傅里叶积分变换而言,它仍是速降函数空间 𝒮 上的自同构。更具体地说,这是一个連續线性等距自同构(可查阅连续线性算子以了解可能导出的其他性质;另,这里所涉及的范数L2范数)。

一个重要的事实是, 𝒮平方可積函數空间 L2稠密,而 L2 是一个希尔伯特空间。这意味着 𝒮 上的任一连续线性算子可以Template:LeL2 上且保持连续性,且这样的扩张唯一。也就是说,若可找到速降函数序列 fn 收敛于一个平方可积函数 f ,那么 f 的傅里叶变换可定义为序列 fn^极限,而不产生收敛性与对 fn 选择的依赖问题。

更进一步地,这个扩张还可以保持等距同构性质,在 L2 上满足这一性质意味着傅里叶变换是一个幺正算子,这就是普朗歇尔定理的内容。

这一变换并不对于所有的平方可积函数都具有一个积分定义式,对于补集 L2L1 中的函数则需要通过极限来定义。故不再称其为积分变换,而是称为傅里叶-普朗歇尔变换或简单称为傅里叶变换普朗歇尔变换

缓增分布

Template:Main如前面所提到的,非零多项式函数不是可积或平方可积的,从而无法使用上面的方法来定义傅里叶变换。这一点可以通过考虑Template:Le来解决。

分布(也称为是一种广义函数)是测试函数空间上的连续线性泛函(即,该空间的連續對偶空間的成员)。而缓增分布则是以速降函数为测试函数的情况。Template:注释

缓增分布 f 的傅里叶变换 f^ 定义为[4]

f^:𝒮:φf(φ^).

也就是说,若记 𝒮 上的傅里叶变换为 0𝒮 的连续对偶空间为 𝒮 ,则 0共轭算子 0*𝒮 上的限制就是缓增分布的傅里叶变换。

这个变换也是可逆的,并且在某种意义上0 的扩张,也就是说可以用它来变换的函数比 𝒮 中的更多。为说明这一点,下面在一维情况下举些例子。

例子

对于实轴上的速降函数 g𝒮() ,其有这样唯一一个由积分定义的连续线性泛函 g 与之对应(这一点实际上只需局部可积函数):

g:φφ(x)g(x) dx,

现在对其做傅里叶变换,得到

g^(φ)=φ^(x)g(x) dx=φ(x)g^(x)dx,

其中第二个等号源于 0 的如下性质(由富比尼定理易证):

dψ^(𝐱)ϕ(𝐱)d𝐱=dψ(𝐱)ϕ^(𝐱)d𝐱.

由此可看出线性泛函 g^ 以前述的方式唯一对应于函数 g^ (唯一性在这样意义上理解:对应相同线性泛函的函数几乎处处相等)。在这个意义上它对于 g𝒮 是与 0 的定义相重合的。

然而并非所有连续线性泛函都有这样的积分表示,如所谓求值泛函。0处的求值泛函作用于每个函数时,都给出该函数在0处的值:

δ0:𝒮:φφ(0),

它正是所谓狄拉克δ函数,其傅里叶变换满足

δ0^(φ)=δ0(φ^)=φ^(0)=φ(x) dx,

而这正是常值函数 1 如前述方式对应的线性泛函 1 。也就是说在这个意义上,狄拉克δ“函数”的傅里叶变换是 1 。从频谱意义上理解,这意味着常值函数的频谱集中在零频率处(或者说周期无穷大) exp(ikx)|k=0=1 的情况。

变换参数的常见约定

如前面提到的,傅里叶积分变换中有可调整的参数,对它们的调整不会造成变换性质的显著变化,而仅仅是对变换结果的定义域或值域进行了放缩。通过显式地引入参数 a,b ,傅里叶积分变换的通式可以写为[1]

f^(𝐤):=(|b|(2π)1a)d/2df(𝐱)eib𝐤𝐱d𝐱,

相应的满足 1=id=1 的逆变换定义为

fˇ(𝐱):=(|b|(2π)1+a)d/2df(𝐤)eib𝐤𝐱d𝐤.

特殊参数选择的特性
编号 选择 特性
1 a=0 变换与逆变换有相同的前置因子,且变换是幺正
2 (2π)1a=|b| 正变换没有前置因子
3 (2π)1+a=|b| 逆变换没有前置因子
4 |b|=1 变换的自变量 k 对应于角频率
5 |b|=2π 变换的自变量 k 对应于频率
6 b<0 exp(ix)贡献正频率部分的谱
7 b>0 exp(ix)贡献正频率部分的谱

一些选择组合因同时满足上面多个特性或在特定领域中自然出现而变得常见,列在下表。

一些常见选择[1]
编号 a b 满足特性 场景
1 0 1 1、4、7 现代物理
2 1 -1 2、4、6 纯数学、系统工程
3 -1 1 3、4、7 传统物理
4 0 2π 1、2、3、5、6 信号处理

本节中使用的约定是 (a,b)=(0,2π)

应用

傅里叶变换在医学数据科学物理学声学光学结构力学量子力学数论组合数学概率论统计学信号处理密码学大氣科學海洋学通讯金融等领域都有着广泛的应用。

数学

Template:Expand section

物理学

在处理具有波动方程背景的函数时,频域的信息处理起来通常更为方便,且信号的滤波等频域操作在器件方面有简单的实现。

除了力学振动、振荡电路等有明显波动方程背景的问题外,也有一些其他情况使得傅里叶变换在理论中自然地出现,如:

基本性质

下面性质的更直观的写法可参见常用傅里叶变换表

线性性质

傅里叶变换是线性映射。也就是说对于 α,β,f,gDom(), 也存在 (αf+βg)=α(f)+β(g).

平移性质

可定义函数间的映射 τα 满足

τα(f)(x)=f(xα),

称为平移算子。其可以推广到缓增分布上,定义为共轭算子 τ~α=τα*,也就是说对于

ϕ𝒮,φ𝒮,ϕ(ταφ)=τ~αf(φ),

其中省略了表示平移算子作用于函数所需的括号。下文不再区分函数与分布的平移,采用相同记号。

傅里叶变换与平移算子满足如下关系:

ταf^(k)=exp(2πiαk)f^(k),

反过来,对于函数 g(x):=exp(2πiαx)f(x) ,也有 g^(k)=ταf^(k)

也就是说平移与相移相关联。

放缩性质

f(ax)    1|a|f^(ξa), a0

a=1 的情况即给出所谓反射性质

导数关系

傅里叶变换与导数算子满足如下关系:

α(f)^(k)=exp(2πik)αf^(k),

其中 α 是高阶导数,多元情况则一般化为多重指标

反过来,对于函数 g(x):=(2πix)αf(x) ,也有 g^(k)=α(f^)(k)

也就是说求导在乘上频率是相关联的。

这里的导数算子也可以是缓增分布的导数算子,同样由共轭算子定义为 ~α=(1)α(α)*

卷积定理

Template:Further

若函数 f,g 有傅里叶变换,且存在卷积 f*g:xf(xξ)g(ξ)dξ ,则该卷积的傅里叶变换即 f^g^

也可以推广地定义分布 ϕ 与测试函数 φ 的卷积,这时同样有 ϕ*φ^=φ^ϕ^

互相关定理

Template:Expand English

实虚部与奇偶性

Template:Expand English

帕塞瓦尔定理

Template:Main L2 上的傅里叶变换是等距映射。或者说,若函数 f(x) 平方可積,则+|f(x)|2dx=+|f^(ω)|2dω

黎曼-勒贝格定理

Template:Main

本征函数

Template:Expand English

傅里叶变换的其他变种

傅里叶级数

Template:Main 连续形式的傅里叶变换其实是傅里叶级数(Fourier series)的推广,因为积分其实是一种极限形式的求和算子而已。对于周期函数,其傅里叶级数是存在的:

f(x)=n=Fneinx,

其中Fn为复振幅。对于实值函数,函数的傅里叶级数可以写成:

f(x)=a02+n=1[ancos(nx)+bnsin(nx)]

其中anbn频率分量的振幅。

傅里叶分析最初是研究周期性现象,即傅里叶级数的,后来通过傅里叶变换将其推广到了非周期性现象。理解这种推广过程的一种方式是将非周期性现象视为周期性现象的一个特例,即其周期为无限长。

离散时间傅里叶变换

Template:Main

离散傅里叶变换是离散时间傅里叶变换(DTFT)的特例(有时作为后者的近似)。DTFT在时域上离散,在频域上则是周期的。DTFT可以被看作是傅里叶级数的逆轉換。

离散傅里叶变换

Template:Main 为了在科学计算和数字信号处理等领域使用计算机进行傅里叶变换,必须将函数xn定义在离散点而非连续域内,且须满足有限性周期性条件。这种情况下,使用离散傅里叶变换,将函数xn表示为下面的求和形式:

Xk=n=0N1xnei2πNknk=0,,N1

其中Xk是傅里叶振幅。直接使用这个公式计算的计算复杂度𝒪(n2),而快速傅里叶变换(FFT)可以将复杂度改进为𝒪(nlogn)。计算复杂度的降低以及数字电路计算能力的发展使得DFT成为在信号处理领域十分实用且重要的方法。

在阿贝尔群上的统一描述

Template:Main 以上的傅里叶变换都可以被统一描述为任意局部紧致阿贝尔群上的傅里叶变换。这一问题属于调和分析的范畴。在调和分析中,一个变换从一个群变换到它的对偶群。此外,将傅里叶变换与卷积相联系的卷积定理在调和分析中也有类似的结论。傅里叶变换的广义理论基础参见庞特里亚金对偶性中的介绍。

时频分析变换

Template:Main 小波变换Chirplet变换分数傅里叶变换的都是为了得到时间信号的频率信息。同时解析频率和时间的能力在数学上受不确定性原理的限制。

傅里叶变换家族

Template:Main

下表列出了傅里叶变换家族的成员。容易发现,函数在时(频)域的离散对应于其像函数在频(时)域的周期性。反之连续则意味着在对应域的信号的非周期性。下表给出详细的情形:

变换 时间 频率
连续傅里叶变换 連續,非週期性 連續,非週期性
傅里叶级数 連續,週期性 離散,非週期性
离散时间傅里叶变换 離散,非週期性 連續,週期性
离散傅里叶变换 離散,週期性 離散,週期性

傅里叶-斯蒂尔杰斯变换

Template:Expand English

常用傅里叶变换表

下面的表记录了一些封闭形式的傅立叶变换。对于函数f(x), g(x)h(x),它们的傅立叶变换分别表示为f^, g^h^。只包含了三种最常见的形式。注意条目105给出了一个函数的傅里叶变换与其原函数,这可以看作是傅里叶变换及其逆变换的关系。

函数关系

下表列出的常用的傅里叶变换对可以在Template:HarvtxtTemplate:Harvtxt中找到。

函数 傅立叶变换
么正,普通的频率
傅立叶变换
么正,角频率
傅立叶变换
非么正,角频率
注释
f(x) f^(ξ)=

f(x)e2πixξdx

f^(ω)=

12πf(x)eiωxdx

f^(ν)=

f(x)eiνxdx

基本定义
101 af(x)+bg(x) af^(ξ)+bg^(ξ) af^(ω)+bg^(ω) af^(ν)+bg^(ν) 线性性质
102 f(xa) e2πiaξf^(ξ) eiaωf^(ω) eiaνf^(ν) 时域平移
103 e2πiaxf(x) f^(ξa) f^(ω2πa) f^(ν2πa) 频域平移,变换102的频域对应
104 f(ax) 1|a|f^(ξa) 1|a|f^(ωa) 1|a|f^(νa) 在时域中定标。如果|a|值较大,则f(ax)会收缩到原点附近,而1|a|f^(ωa)会扩散并变得扁平。当|a|趋向无穷时,f(ax)成为狄拉克δ函数
105 f^(x) f(ξ) f(ω) 2πf(ν) 傅里叶变换的二元性性质。这里f^的计算需要运用与傅里叶变换那一列同样的方法。通过交换变量xξων得到。
106 dnf(x)dxn (2πiξ)nf^(ξ) (iω)nf^(ω) (iν)nf^(ν) 傅里叶变换的微分性质
107 xnf(x) (i2π)ndnf^(ξ)dξn indnf^(ω)dωn indnf^(ν)dνn 变换106的频域对应
108 (f*g)(x) f^(ξ)g^(ξ) 2πf^(ω)g^(ω) f^(ν)g^(ν) 记号f*g表示fg的卷积—这就是卷积定理
109 f(x)g(x) (f^*g^)(ξ) (f^*g^)(ω)2π 12π(f^*g^)(ν) 变换108的频域对应。
110 f(x)是实变函数 f^(ξ)=f^(ξ) f^(ω)=f^(ω) f^(ν)=f^(ν) 埃尔米特对称。z表示复共轭
111 f(x)是实偶函数 f^(ω), f^(ξ)f^(ν)都是实偶函数
112 f(x)是实奇函数 f^(ω), f^(ξ)f^(ν)都是奇函数
113 f(x) f^(ξ) f^(ω) f^(ν) 复共轭,110的一般化

平方可积函数

时域信号 角频率表示的
傅里叶变换
弧频率表示的
傅里叶变换
注释
g(t)

12πG(ω)eiωtdω
G(ω)

12πg(t)eiωtdt
G(f)

g(t)ei2πftdt
10 rect(at) 12πa2sinc(ω2πa) 1|a|sinc(fa) 矩形脉冲和归一化的sinc函数
11 sinc(at) 12πa2rect(ω2πa) 1|a|rect(fa) 变换10的频域对应。矩形函数是理想的低通滤波器,sinc函数是这类滤波器对反因果冲击的响应。
12 sinc2(at) 12πa2tri(ω2πa) 1|a|tri(fa) tri三角形函数
13 tri(at) 12πa2sinc2(ω2πa) 1|a|sinc2(fa) 变换12的频域对应
14 eαt2 12αeω24α παe(πf)2α 高斯函数exp(αt2)的傅里叶变换是其本身;只有当Re(α)>0时,该函数可积的
15 eiat2=eαt2|α=ia 12aei(ω24aπ4) πaei(π2f2aπ4) 光学领域应用较多
16 cos(at2) 12acos(ω24aπ4) πacos(π2f2aπ4)
17 sin(at2) 12asin(ω24aπ4) πasin(π2f2aπ4)
18 ea|t| 2πaa2+ω2 2aa2+4π2f2 a>0
19 1|t| 1|ω| 1|f| 变换本身就是一个公式
20 J0(t) 2πrect(ω2)1ω2 2rect(πf)14π2f2 J0(t)0阶第一类贝塞尔函数
21 Jn(t) 2π(i)nTn(ω)rect(ω2)1ω2 2(i)nTn(2πf)rect(πf)14π2f2 上一个变换的推广形式; Tn (t)第一类切比雪夫多项式
22 Jn(t)t 2πin(i)nUn1(ω)

   1ω2rect(ω2)

2in(i)nUn1(2πf)

   14π2f2rect(πf)

Un (t)第二类切比雪夫多项式

分布

时域信号 角频率表示的
傅里叶变换
弧频率表示的
傅里叶变换
注释
g(t)

12πG(ω)eiωtdω
G(ω)

12πg(t)eiωtdt
G(f)

g(t)ei2πftdt
基本定义
23 1 2πδ(ω) δ(f) δ(ω)代表狄拉克δ函数分布.这个变换展示了狄拉克δ函数的重要性:该函数是常函数的傅立叶变换
24 δ(t) 12π 1 变换23的频域对应
25 eiat 2πδ(ωa) δ(fa2π) 由变换103和23得到
26 cos(at) 2πδ(ωa)+δ(ω+a)2 δ(fa2π)+δ(f+a2π)2 由变换101和25得到,应用了欧拉公式cos(at)=(eiat+eiat)/2.
27 sin(at) 2πδ(ωa)δ(ω+a)2i δ(fa2π)δ(f+a2π)2i 由变换101和25得到
28 tn in2πδ(n)(ω) (i2π)nδ(n)(f) 这里, n是一个自然数. δ(n)(ω)是狄拉克δ函数分布的n阶微分。这个变换是根据变换107和24得到的。将此变换与101结合使用,我们可以变换所有多項式函数。
29 1t iπ2sgn(ω) iπsgn(f) 此处sgn(ω)符号函数;注意此变换与变换107和24是一致的.
30 1tn iπ2(iω)n1(n1)!sgn(ω) iπ(i2πf)n1(n1)!sgn(f) 变换29的推广
31 sgn(t) 2π1i ω 1iπf 变换29的频域对应
32 u(t) π2(1iπω+δ(ω)) 12(1iπf+δ(f)) 此处u(t)单位阶跃函数;此变换根据变换101和31得到.
33 eatu(t) 12π(a+iω) 1a+i2πf u(t)单位阶跃函数,且a>0.
34 n=δ(tnT) 2πTk=δ(ωk2πT) 1Tk=δ(fkT) Template:Le——有助于解释或理解从连续到离散时间的转变.

二元函数

时域信号 傅立叶变换
单一,普通频率
傅立叶变换
么正,角频率
傅立叶变换
非么正,角频率
400 f(x,y) f^(ξx,ξy)=
f(x,y)e2πi(ξxx+ξyy)dxdy
f^(ωx,ωy)=
12πf(x,y)ei(ωxx+ωyy)dxdy
f^(νx,νy)=
f(x,y)ei(νxx+νyy)dxdy
401 eπ(a2x2+b2y2) 1|ab|eπ(ξx2/a2+ξy2/b2) 12π|ab|e(ωx2/a2+ωy2/b2)4π 1|ab|e(νx2/a2+νy2/b2)4π
402 circ(x2+y2) J1(2πξx2+ξy2)ξx2+ξy2 J1(ωx2+ωy2)ωx2+ωy2 2πJ1(νx2+νy2)νx2+νy2
注释

400: 变量ξxξyωxωyνxνy为实数。二重积分是对整个平面积分。

401: 这两个函数都是高斯函数,而且可能不具有单位体积。

402: 此圆有单位半径,如果把 circ(t)认作阶梯函数 u(1t); Airy分布用 J1一阶第一类贝塞尔函数)表达。Template:Harv

三元函数

时域信号 角频率表示的
傅里叶变换
弧频率表示的
傅里叶变换
注释
circ(x2+y2+z2) 2πsin[ω]2πfrcos[ω]ω3 4πsin[2πfr]2πfrcos[2πfr](2πfr)3 此球有单位半径;fr是频率矢量的量值{fx,fy,fz}.

参见

Template:Div col

Template:Div col end

參考資料

引用

Template:Reflist

注释

Template:NoteFoot

来源

Template:Refbegin

Template:Refend

外部連結

Template:Authority control