在这一篇
Blog 中,博主打算记录一下周期信号傅里叶级数的一些重要性质。在阅读本文之前,请明确一件事情:对于周期信号我们讲的是傅里叶级数展开。对于非周期信号我们讲傅里叶变换。
一、什么样的周期信号才能够做傅里叶展开?
当时傅里叶的论文中有这样一句话:“所有的连续时间信号都能够表示成成谐波关系的复指数信号的加权和”。在当时引来了像拉普拉斯等人的强烈反对。傅里叶认为:对于周期方波信号,只要我取的正弦信号足够多,那么我就一定能够完美地拟合出来。如下图所示:
但实际上并不是这样,无论我们的正弦信号取了多少,在原方波的间断点处不可避免地会出现震荡和超量。超量的幅度并不会随着我们所选取的正弦波的数量增多而减少,选取的正弦波的数量增多只会使得超量的震荡频率更大,并且朝着间断点处压缩。这就是所谓的 “吉布斯现象”。
那么,到底什么样的连续时间信号才能够展开成傅里叶级数呢?我们下面给出三个条件(其中条件三是狄里克雷第一定理)。连续时间的周期信号只要满足三者之一,就可以展开成傅里叶级数:
【条件一】:信号全部连续
【条件二】:信号在一个周期内能量有限:
T01∫T0∣x(t)∣2dt<∞
【条件三】:信号在一个周期内绝对可积:
T01∫T0∣x(t)∣<∞
二、周期信号傅里叶级数的重要性质
2.1 线性
首先我们给出定义:
连续时间信号傅里叶级数的线性性:若信号
x1(t)的傅里叶级数是
ak,信号
x2(t) 的傅里叶级数是
bk,那么,信号
Ax1(t)+Bx2(t) 的傅里叶级数就是
Aak+Bbk。
首先,我们先知道求傅里叶系数的过程就是一个积分的过程,如果我们用
∫ 表示计算傅里叶系数(当然这样的写法不够准确,但是为了说明问题暂且先这样用)。
所以,我们如果对信号
Ax1(t)+Bx2(t) 求傅里叶系数:
∫Ax1(t)+Bx2(t),就可以这样变换:
∫Ax1(t)+Bx2(t)=∫Ax1(t)+∫Bx2(t)=A∫x1(t)+B∫x2(t)
而我们知道:
ak=∫x1(t);
bk=∫x2(t),所以信号
Ax1(t)+Bx2(t) 的傅里叶系数就是:
Aak+Bbk
2.2 时移特性
我们同样先看定义:
连续时间信号傅里叶级数的时移特性:假设
x(t) 的傅里叶系数是
ak,那么
x(t−t0) 的傅里叶系数就是
ake−jkω0t0,也即是说,时移并不会改变傅里叶级数的幅度,改变的是傅里叶级数的相位。
首先根据定义,我们可以得到:
ak=T01∫T0x(t)e−jkω0tdt
下面我们对信号
x(t−t0) 求傅里叶级数:因为
x(t−t0) 只是时移,所以周期并不会改变,仍为
T0
T01∫x(t−t0)e−jkω0tdt=e−jkω0t0T01∫x(t−t0)e−jkω0(t−t0)d(t−t0)=e−jkω0t0T01∫x(t)e−jkω0tdt
注意:在最后一行我们做了变量代换:令
t=t−t0
2.3 尺度变换
这个稍微难理解一点点,我们先推导,再给出定义:
首先对于信号
x(t) ,我们有:
ak=T01∫T0x(t)e−jkω0tdt
而我们知道,尺度变换相当于对信号做拉伸或者压缩,自然会改变原信号的周期。例如,如果信号
x(t) 的周期是 T ,那么信号
x(at) 的周期就是
T1=aT。那么角频率就是
ω1=aω0。好,下面我们计算
x(at) 的傅里叶系数,根据定义,有:
T11∫T1x(at)e−jkω1tdt=T0a∫aTx(at)e−jkaω0tdt
所以,我们对于计算带时移的周期信号的傅里叶级数,我们首先要做的就是明确这个信号和原本信号周期、角频率之间的关系。然后才能根据定义求解。
那么,下面我们给出定义:对于带时移的周期信号
x(at),其傅里叶级数表示为:
bk=T0a∫aT0x(at)e−jkaω0tdt(1)
更进一步讲,如果我们使用
τ 表示
at,那么,因为我们的积分范围就是在一个周期
aT0 内,因此,我们知道,
t 的取值范围是:
0≤t≤aT0(当然这个范围并不唯一)。那么
τ 的取值范围就是
0≤τ≤T0 下面我们就把
τ 带入(1) 式:
bk=T0a∫T0x(τ)ejkω0τd(aτ) =a1T0a∫T0x(τ)ejkω0τdτ=T01∫T0x(τ)ejkω0τdτ=ak
2.4 反转
这里我们直接给出结论:对于周期信号
x(t),其傅里叶系数是
ak。那么其反转信号
x(−t) 的傅里叶系数就是
a−k,即也是对于
ak 的反转。进一步讲,若
x(t) 是偶信号,那么
ak 也是偶的。如果
x(t)是奇信号,那么
a−k=−ak
2.5 时域相乘等价于频域卷积
我们先给出结论:若
x(t),y(t) 的傅里叶级数分别是
ak,bk,那么有:
x(t)y(t) 的傅里叶级数就是:
ak∗bk
我们证明一下,直接用求傅里叶系数的公式即可:
T1∫Tx(t)y(t)e−jω0ktdt=T1∫Tl=−∞∑+∞akejω0lty(t)e−jω0ktdt=l=−∞∑+∞akT1∫Ty(t)ejω0(l−k)tdt=l=−∞∑+∞akbl−k=ak∗bk
2.6 周期卷积定理(和2.5 对偶)
对于连续时间的、周期相同的周期信号
x(t),y(t),我们在一个周期内的卷积,就等于它们对应的傅里叶级数相乘,再乘上周期。表述为:
∫Tx(τ)y(t−τ)dτ=Takbk
2.7 共轭以及共轭对称性
首先我们看看第一个结论:假设周期信号
x(t) 的傅里叶级数是
ak,那么,如果取
x(t) 的共轭:
x∗(t),那么这个共轭信号
x∗(t) 的傅里叶系数就应该要对
ak 取共轭,并且进行反转,即:
a−k∗
用数学语言表述即为:
x(t)
FSak x∗(t)
FSa−k∗
更进一步:如果
x(t) 是实信号,即
x(t)=x∗(t),那么应有:
ak=a−k∗,不过我们更常用的是:
ak∗=a−k
下面的几个推导大家需要跟上,否则可能会造成混乱:
【1】若
x(t) 是实信号,而且是偶函数,即
x(t)=x(−t)。那么因为
x(−t) 的傅里叶级数应该是
a−k。所以应有:
ak=a−k,同时,我们根据
x(t) 是实信号,又可以得到:
ak=a−k∗。综上,我们得到了一串等式:
ak=a−k=ak∗。前面一个等号说明
ak 也是偶的、后面的等号说明
ak 是实的。 即若
x(t) 是实偶函数,那么其频谱将会是实偶函数
【2】若
x(t) 是实信号,同时又是奇函数
x(t)=−x(−t),因为
−x(−t) 的傅里叶级数是
−a−k。所以我们也可以得到一串等式:
ak=−a−k=−ak∗。前一个等号表示
ak 是奇函数、后一个等号表示
ak 是纯虚数,因此,若
x(t) 是实奇函数,那么其频谱
ak 将会是纯虚奇函数。
2.8 微分性
我们先给出结论:若
x(t) 的傅里叶级数是
ak,那么
∂t∂x(t) 的傅里叶系数就是:
akjω0t
下面给出证明:我们先从
x(t) 的傅里叶展开表达式入手。由于:
x(t)=k=−∞∑+∞akejkω0t
下面我们直接对等式两边求导,得:
∂t∂x(t)=k=−∞∑+∞(akjkω0) ejkω0t
因此,新的傅里叶系数就是:
akjω0t
利用微分型,我们一样可以计算出连续时间周期矩形信号的频谱。这里不详细展开,但是提几个突破口:
- 对周期矩形信号求导,将会得到有上有下的单位冲激函数。
- 周期为 T 的单位冲激函数串的傅里叶系数的幅值都是
ak=T1
- 有上有下的单位冲激函数是可以写成单位冲激函数串右移
T1 和左移
T1 的差值。再利用傅里叶系数时移的特点,即可求出周期矩形信号的频谱。
三、帕斯瓦尔定理
这个定理的表示简单,但是证明是比较困难的。下面我们直接给出:
对于一个周期信号的平均功率的计算,可以把它的所有傅里叶系数的平方加起来。
T1∫T∣x(t)∣2dt=k=−∞∑+∞∣ak∣2
同时,我们也引入周期信号功率谱的概念:周期信号的
∣ak∣2 随着
kω0 变化的情况称为功率谱。 对应地,非周期信号还有功率谱密度,我们以后再来介绍。