打爆傅里叶变换的狗头:为什么离散谱线会变成连续频谱?

建议先看我前面写的傅里叶文章后,再来读这个,我这里不在详细的回顾之前的推导。
https://mp.weixin.qq.com/s/5C5uVWAxcGw-hxEDbZGDSw
https://mp.weixin.qq.com/s/P_yuV0OPA2citU_bMubtWg

傅里叶级数里明明是一根根离散谱线,为什么到了傅里叶变换里,就突然变成了一条连续曲线?
傅里叶级数的系数是 Cn,傅里叶变换里为什么不继续叫 C(ω),而要重新定义一个 F(ω)

这篇文章想从傅里叶级数的公式出发并计算 Cn,为了是回忆一下傅里叶级数,然后一步一步推到傅里叶变换,并重点解释中间最容易混乱的几个地方。

不要看到计算就害怕,因为我写的真的很清楚了。

核心思想可以先放在这里:

傅里叶变换不是凭空出现的,它可以看作是傅里叶级数在周期 T 时的极限形式。


一、傅里叶级数系数 Cn

对于一个周期为 T 的连续周期信号,其复指数形式的傅里叶级数为:

f(t)=n=Cnejnω0t

其中:

ω0=2πT

ω0 称为基波角频率,也叫基频角频率

注意:

  • ω0 的单位是 rad/s;
  • f0=1T 称为基频,单位是 Hz;
  • “基波”通常指频率为 ω0 的那一项信号分量,而不是 ω0 这个数本身。

频域中的谱线只会出现在基波角频率的整数倍位置:

ω=nω0,n=0,±1,±2,

这些离散频率点上的谱线共同组成了周期信号的频谱。


二、Cn 的物理意义

Cn 表示频率 ω=nω0 处那根谱线的复振幅。

也就是说:

  • C0:表示 ω=0 处谱线的复振幅,也就是直流分量
  • C1:表示 ω=ω0 处谱线的复振幅
  • C1:表示 ω=ω0 处谱线的复振幅
  • C2:表示 ω=2ω0 处谱线的复振幅
  • C2:表示 ω=2ω0 处谱线的复振幅

一般地:

Cnω=nω0 处谱线的复振幅

因此:

  • 谱线的位置由 nω0 决定;
  • 谱线的复振幅由 Cn 决定;
  • 谱线的高度通常画为 |Cn|
  • 谱线的相位通常画为 Cn

image.png|679


三、如何计算 Cn

傅里叶级数系数由下面的积分给出:

Cn=1TT/2T/2f(t)ejnω0tdt

这个积分的作用是:用 ejnω0t 去“试探”信号 f(t) 中是否含有频率为 nω0 的成分。

哲学:你不试探一下你自己,你就不知道自己有多菜。🤣

这里的 n 可以取:

n=0,±1,±2,

也就是说,n 表示第几条谱线,或者说第几阶谐波。

如果计算结果:

Cn0

说明频率 ω=nω0 处有谱线。

如果计算结果:

Cn=0

说明频率 ω=nω0 处没有谱线。


四、计算具体谱线

如果题目要求计算:

n=0,±1,±2

对应的谱线,那么只需要分别计算:

C0, C1, C1, C2, C2

1. 直流分量

C0=1TT/2T/2f(t)dt

对应频率:

ω=0

2. 正一阶谐波

C1=1TT/2T/2f(t)ejω0tdt

对应频率:

ω=ω0

3. 负一阶谐波

C1=1TT/2T/2f(t)ejω0tdt

对应频率:

ω=ω0

4. 正二阶谐波

C2=1TT/2T/2f(t)ej2ω0tdt

对应频率:

ω=2ω0

5. 负二阶谐波

C2=1TT/2T/2f(t)ej2ω0tdt

对应频率:

ω=2ω0


五、从 Cn 到频谱

傅里叶级数的频谱由一根根离散谱线组成。第 n 条谱线位于

ω=nω0

其复振幅为

Cn=|Cn|ejCn

其中 |Cn| 是幅度谱,Cn 是相位谱。计算频谱时,只需在离散频率点 ω=nω0 上计算对应的 Cn

Cn=1TT/2T/2f(t)ejnω0t,dt


六、下面傅里叶变换开始了:大家坐好了,准备go

现在我们把 Cn 代回傅里叶级数:

f(t)=n=[1TT/2T/2f(τ)ejnω0τdτ]ejnω0t

因为:

ω0=2πT

所以:

1T=ω02π

于是上式可以改写为:

f(t)=n=[ω02πT/2T/2f(τ)ejnω0τdτ]ejnω0t

12π 提出来,并把 ω0 放到求和项末尾:

f(t)=12πn=[T/2T/2f(τ)ejnω0τdτ]ejnω0tω0

这个形式非常关键,因为它已经很像一个积分了:

n=()ω0

看上面的公式就可以分析出:当频率间隔 ω0 变得越来越小时,它就会变成连续积分。


七、当周期 T 时,发生了什么?

现在进入最关键的一步。

让周期:

T

这一步可以理解为:

把周期信号的周期拉得越来越长,使相邻两次重复之间的距离越来越远,最后得到一个非周期信号。

例如一个单次脉冲、一次雷声、一个瞬间衰减的波形,都可以看作是周期无限大时的极限情形。

T 时,公式中的几个量会同时发生变化。

1. 积分区间从有限变成无限

原本傅里叶级数只在一个周期内积分:

T/2T/2

T 时,积分区间变成整条时间轴:

2. 谱线间隔变成无穷小

因为:

ω0=2πT

所以当 T 时:

ω00

也就是说,原本相邻谱线之间的间隔越来越小。

在极限情况下,可以把离散频率间隔 ω0 看成微分量:

ω0dω

3. 离散频率变成连续频率

原来的频率点是:

nω0

它只能取一串离散值。

但当 ω00 时,这些频率点变得无限密集,最终铺满整条频率轴,于是:

nω0ω

其中 ω 成为连续频率变量。

4. 求和变成积分

原来的求和形式是:

n=()ω0

ω0dω,并且离散频率点变成连续频率后,它就变成了积分:

()dω


八、傅里叶变换自然出现了

回到刚才的公式:

f(t)=12πn=[T/2T/2f(τ)ejnω0τdτ]ejnω0tω0

T 时,有:

T/2T/2nω0ωω0dωn=()ω0()dω

于是得到:

f(t)=12π[f(τ)ejωτdτ]ejωtdω

这时,我们把中括号里的部分单独定义为 F(ω)

F(ω)=f(τ)ejωτdτ

于是得到傅里叶变换对:

F(ω)=f(t)ejωtdtf(t)=12πF(ω)ejωtdω

这就是傅里叶变换。


九、为什么不直接让 Cn 变成 C(ω)

很多人学到这里会问:

傅里叶级数的系数是 Cn,现在频率连续了,为什么不直接写成 C(ω)

原因是:

T 时,单根谱线的复振幅 Cn 会趋近于 0。

这不是记号问题,而是本质问题。

回看傅里叶级数系数:

Cn=1TT/2T/2f(τ)ejnω0τdτ

这里有一个关键因子:

1T

当:

T

时:

1T0

对于非周期有限能量信号来说,信号只在有限时间内有明显变化,而周期被拉到无限大。于是每一根离散谱线分到的复振幅越来越小。最后,单根谱线的复振幅趋近于 0。所以如果我们强行令:

CnC(ω)

那么这个 C(ω) 通常会变成处处为 0,无法描述频域结构。

这就像你把一杯水摊到一条无限长的线上,每一个孤立点分到的水量都是 0。但这并不代表整条线上没有水,而是水被连续地分布开了。


十、傅里叶变换其实对应 TCn 的极限

为了得到一个有限的、有意义的频域函数,我们不能直接看 Cn,而要看:

TCn

由傅里叶级数系数公式:

Cn=1TT/2T/2f(τ)ejnω0τdτ

两边乘以 T,得到:

TCn=T/2T/2f(τ)ejnω0τdτ

T 时:

TCnf(τ)ejωτdτ

于是:

F(ω)=limTTCn

其中:

ω=nω0

这句话非常重要:

傅里叶变换 F(ω) 不是傅里叶级数系数 Cn 的简单连续化,而是 TCnT 时的极限。

换句话说:

Cn1TF(nω0)

又因为:

ω0=2πT

所以:

1T=ω02π

因此也可以写成:

Cnω02πF(nω0)

这正好解释了傅里叶逆变换中为什么会出现 12πdω


十一、连续频谱的纵坐标到底是什么意思?

傅里叶级数中,Cn 是一根谱线的复振幅。

但傅里叶变换中,F(ω) 不是某个孤立频率点上的真实振幅。

更准确地说:

在角频率形式下,F(ω) 可以理解为相对于角频率 ω 的频谱密度。

在连续频谱中,单个频率点没有宽度,因此单个频率点本身对应的实际贡献为 0。

真正有意义的是一个很小频率区间内的贡献。

傅里叶逆变换为:

f(t)=12πF(ω)ejωtdω

其中一个很小频带 [ω,ω+dω] 对信号的贡献大约是:

12πF(ω)ejωtdω

所以连续频谱里真正参与叠加的不是孤立的 F(ω),而是带有频率宽度的微元项:

F(ω)dω

在角频率 ω 的写法下,如果 f(t) 的单位是 V,则:

F(ω)

的单位是:

Vs

因为:

F(ω)=f(t)ejωtdt

积分中 dt 的单位是 s。

dω 的单位是 rad/s,rad 通常视为无量纲,所以 dω 的量纲可看作 1/s。因此:

F(ω)dω

的单位回到 V。

再乘上无量纲的 12πejωt,最后得到的仍然是原信号的单位 V。

如果使用普通频率 f,单位是 Hz,那么傅里叶变换常写成:

X(f)=x(t)ej2πftdtx(t)=X(f)ej2πftdf

这时 X(f) 可以理解为相对于 Hz 的频谱密度。

如果 x(t) 的单位是 V,那么:

X(f)

的单位是:

Vs

由于:

1 Hz=1/s

所以也可以写成:

VHz

注意:角频率 ω 和普通频率 f 只差一个 2π

ω=2πf

所以不同教材里单位和 2π 的位置会略有不同,本质是同一件事。


十二、一个生活类比:速度和位移

“单点频率上的贡献为 0”这件事很反直觉。

可以用速度来类比。

假设一辆车正在高速路上行驶,速度表显示:

120 km/h

现在问一个问题:

在正午 12 点整这个没有时间宽度的瞬间,汽车前进了多少公里?

答案是:

0

因为一个孤立瞬间没有时间长度。

再快的速度,如果乘以:

dt=0

得到的位移也是 0。

真正有意义的是一小段时间内的位移:

ds=v(t)dt

速度 v(t) 本身不是位移,它是“位移对时间的密度”。

傅里叶变换中的 F(ω) 也是类似的。

F(ω) 本身不是某个孤立频点上的真实振幅,它更像是“振幅对频率的密度”。

真正参与合成的是:

F(ω)dω

再乘上基函数:

ejωt

最后对所有频率积分,才还原出原信号。


十三、傅里叶级数和傅里叶变换的公式对齐

现在把两套公式放在一起看,会更清楚。

1. 傅里叶级数:周期信号,离散频谱

分析公式:

Cn=1TT/2T/2f(t)ejnω0t,dt

合成公式:

f(t)=n=Cnejnω0t

其中:

ω0=2πT

特点是:频率只出现在 ω=nω0 上,频谱是一根根离散谱线,Cn 是对应谱线的复振幅。周期信号通常讨论平均功率,而不是总能量。

2. 傅里叶变换:非周期信号,连续频谱

分析公式:

F(ω)=f(t)ejωt,dt

合成公式:

f(t)=12πF(ω)ejωt,dω

特点是:频率 ω 连续变化,频谱是连续曲线,F(ω) 可以理解为相对于角频率的频谱密度。非周期有限能量信号通常讨论能量分布。

3. 两者之间的关键关系

当周期 T 越来越大时,

ω0=2πT0

所以离散频率点

ω=nω0

会越来越密,最后趋向连续频率变量 ω

傅里叶变换和傅里叶级数系数的关系为:

F(ω)=limTTCn,ω=nω0

也可以近似写成:

Cn1TF(nω0)

又因为:

1T=ω02π

所以:

Cnω02πF(nω0)

代回傅里叶级数合成公式:

f(t)=n=Cnejnω0t

得到:

f(t)n=ω02πF(nω0)ejnω0t

ω00 时,上式变成积分:

f(t)=12πF(ω)ejωt,dω

这就是傅里叶逆变换。

十四、最核心的逻辑链条

整件事可以压缩成一条逻辑链:

Tω0=2πT谱线越来越密

同时:

Cn1T单根谱线越来越矮

最后当 T 时:

ω0dω,nω0ω()ω0(),dωCn0,TCnF(ω)

所以,傅里叶变换不是“把 Cn 改名成 F(ω)”,而是在“谱线无限变密、单根无限变矮”的极限中,提取出一个有限的频谱密度函数 F(ω)

我写完,你读完,那就撒花:🎉🎉🎉🎉🎉🎉🎉🎉🎉。

暂无评论

发送评论 编辑评论


				
|´・ω・)ノ
ヾ(≧∇≦*)ゝ
(☆ω☆)
(╯‵□′)╯︵┴─┴
 ̄﹃ ̄
(/ω\)
∠( ᐛ 」∠)_
(๑•̀ㅁ•́ฅ)
→_→
୧(๑•̀⌄•́๑)૭
٩(ˊᗜˋ*)و
(ノ°ο°)ノ
(´இ皿இ`)
⌇●﹏●⌇
(ฅ´ω`ฅ)
(╯°A°)╯︵○○○
φ( ̄∇ ̄o)
ヾ(´・ ・`。)ノ"
( ง ᵒ̌皿ᵒ̌)ง⁼³₌₃
(ó﹏ò。)
Σ(っ °Д °;)っ
( ,,´・ω・)ノ"(´っω・`。)
╮(╯▽╰)╭
o(*////▽////*)q
>﹏<
( ๑´•ω•) "(ㆆᴗㆆ)
😂
😀
😅
😊
🙂
🙃
😌
😍
😘
😜
😝
😏
😒
🙄
😳
😡
😔
😫
😱
😭
💩
👻
🙌
🖕
👍
👫
👬
👭
🌚
🌝
🙈
💊
😶
🙏
🍦
🍉
😣
Source: github.com/k4yt3x/flowerhd
颜文字
Emoji
小恐龙
花!
上一篇