信号系统(三)——连续时间傅里叶

将函数分解为基本函数的加权和,是研究函数十分重要的方法之一。比如上一篇中我们将函数分解为了单位冲激函数的加权和(对于离散时间信号是这样,连续信号则是卷积积分),就是将单位冲激函数作为了基本函数。再例如,如果将幂函数作为基本函数,那就是泰勒展开;而以三角函数作为基本函数,那就是本篇的重点——傅里叶展开了。

注 1.
或者,由于三角函数可以进一步分解为复指数函数,傅里叶展开也可以认为是以复指数函数作为基本函数。以下,对复指数信号的分析,也可视作对三角函数的分析。

一、复指数信号的特殊性

为什么复指数函数可以作为基本函数呢?这主要是因为复指数信号在经过线性时不变系统输出后,只会产生幅值上的变化。也就是说,假设复指数信号 e𝑠𝑡 ,其中 𝑠 是任意复数,那么:

e𝑠𝑡→𝐻(𝑠)e𝑠𝑡

幅值会产生 𝐻(𝑠) 倍的变化,而且这个幅值还由 𝑠 唯一确定!我们将这种经过系统作用之后的输出只是输入的常数倍的信号称为特征函数,而这个常数称为特征值。

证明. 为了证明以上结论,首先假设系统的单位冲激响应为 ℎ(𝑡) ,。那么系统对 𝑥(𝑡) 的响应就是:

𝑦(𝑡)=𝑥(𝑡)∗ℎ(𝑡)=∫−∞+∞𝑥(𝜏)ℎ(𝑡−𝜏)𝑑𝜏=∫−∞+∞𝑥(𝑡−𝜏)ℎ(𝜏)𝑑𝜏=∫−∞+∞e𝑠(𝑡−𝜏)ℎ(𝜏)𝑑𝜏=e𝑠𝑡∫−∞+∞e−𝑠𝜏ℎ(𝜏)𝑑𝜏
记 𝐻(𝑠)=∫−∞+∞e−𝑠𝜏ℎ(𝜏)𝑑𝜏 ,假定这个积分收敛,那么上述结论得证。并且结论对离散时间信号亦成立。
□

因此,若信号可以被分解为复指数信号的加权和:

𝑥(𝑡)=𝑎0e𝑠0𝑡+𝑎1e𝑠1𝑡+𝑎2e𝑠2𝑡+…=∑𝑎𝑖e𝑠𝑖𝑡

显然输出可以直接确定为:

𝑦(𝑡)=𝐻(𝑠0)𝑎0e𝑠0𝑡+𝐻(𝑠1)𝑎1e𝑠1𝑡+𝐻(𝑠2)𝑎2e𝑠2𝑡+…=∑𝐻(𝑠𝑖)𝑎𝑖e𝑠𝑖𝑡

这就为研究系统带来了极大的便利,正如上篇文章中那般,只要确定单位冲激响应就能完全确定线性时不变系统对任意信号的响应,现在只要确定 𝐻(𝑠) 也能完全确定线性时不变系统对某些信号的响应。


二、傅里叶级数的计算

在看到将函数分解为复指数函数带来便利之后,我们的下一个问题是:如何将一个函数分解为复指数函数的加权和?首先需要注意到复指数函数是周期信号,因此“只有”周期函数可以被展开。以下内容都事先默认假定要展开的函数是周期函数。

1.以三角函数作为基底

在高等数学中我们学习过傅里叶级数、三角展开,就是将一个函数变换为另一些函数的线性组合——而如果将函数视为向量的话,“另一些函数”的函数簇中的函数两两正交,那么就是将原函数在这些函数张成的无限维空间中正交分解罢了。傅里叶展开所选的正交基是三角函数系 [1,cos𝑘𝑥,sin𝑘𝑥],𝑘=1,2,… ,这其中任取两不同函数,其积在 [−𝜋,𝜋] 上的 2− 范数都等于 0 ,也即两两正交:

∫−𝜋𝜋cos𝑘𝑥=∫−𝜋𝜋sin𝑘𝑥𝑑𝑥=0∫−𝜋𝜋cos𝑘𝑥⋅sin𝑛𝑥𝑑𝑥=0,𝑘≠𝑛∫−𝜋𝜋cos𝑘𝑥⋅cos𝑛𝑥𝑑𝑥=0,𝑘≠𝑛∫−𝜋𝜋sin𝑘𝑥⋅sin𝑛𝑥𝑑𝑥=0,𝑘≠𝑛

注意到我们要求三角函数系中 𝑘 为整数,这就要求 2𝜋 必须是这些三角函数的周期。这只是为了后续书写计算方便,之后会逐步放开限制。假设某周期为 2𝜋 的函数 𝑥(𝑡) 可以傅里叶展开为:

𝑥(𝑡)=𝑎02+∑𝑘=1+∞𝑎𝑘cos𝑘𝑡+𝑏𝑘sin𝑘𝑡

直接对上式两边在 [−𝜋,+𝜋] 区间上积分,注意到上面提到的 ∫−𝜋𝜋cos𝑘𝑥=∫−𝜋𝜋sin𝑘𝑥𝑑𝑥=0 ,那么:

∫−𝜋𝜋𝑥(𝑡)𝑑𝑡=∫−𝜋𝜋𝑎02𝑑𝑡+∑𝑘=1+∞𝑎𝑘∫−𝜋𝜋cos𝑘𝑡𝑑𝑡+𝑏𝑘∫−𝜋𝜋sin𝑘𝑡𝑑𝑡=∫−𝜋𝜋𝑎02𝑑𝑡=𝑎0𝜋

因此得到 𝑎0 的计算式:

𝑎0=1𝜋∫−𝜋𝜋𝑥(𝑡)𝑑𝑡

类似地,将 𝑥(𝑡) 乘上 cos𝑛𝑡 后做同样的操作,并注意三角函数系的性质:

∫−𝜋𝜋𝑥(𝑡)cos𝑛𝑡𝑑𝑡=∫−𝜋𝜋𝑎02cos𝑛𝑡𝑑𝑡+∑𝑘=1+∞𝑎𝑘∫−𝜋𝜋cos𝑘𝑡⋅cos𝑛𝑡𝑑𝑡+𝑏𝑘∫−𝜋𝜋sin𝑘𝑡⋅cos𝑛𝑡𝑑𝑡=𝑎𝑛∫−𝜋𝜋cos2𝑛𝑡𝑑𝑡=𝑎𝑛𝜋

于是:

𝑎𝑛=1𝜋∫−𝜋𝜋𝑥(𝑡)cos𝑛𝑡𝑑𝑡

将 𝑥(𝑡) 乘上 sin𝑛𝑡 后做同样的操作:

∫−𝜋𝜋𝑥(𝑡)sin𝑛𝑡𝑑𝑡=∫−𝜋𝜋𝑎02sin𝑛𝑡𝑑𝑡+∑𝑘=1+∞𝑎𝑘∫−𝜋𝜋cos𝑘𝑡⋅sin𝑛𝑡𝑑𝑡+𝑏𝑘∫−𝜋𝜋sin𝑘𝑡⋅sin𝑛𝑡𝑑𝑡=𝑏𝑛∫−𝜋𝜋sin2𝑛𝑡𝑑𝑡=𝑏𝑛𝜋

于是:

𝑏𝑛=1𝜋∫−𝜋𝜋𝑥(𝑡)sin𝑛𝑡𝑑𝑡

综上,我们得到了周期为 2𝜋 的函数 𝑥(𝑡) 的傅里叶三角项级数的计算式:

{𝑎𝑛=1𝜋∫−𝜋𝜋𝑥(𝑡)cos𝑛𝑡𝑑𝑡𝑏𝑛=1𝜋∫−𝜋𝜋𝑥(𝑡)sin𝑛𝑡𝑑𝑡

那么,若 𝑥(𝑡) 的周期不为 2𝜋 呢?不必重新计算,因为只是相当于做了时间伸缩变换。假设周期为 𝑇=2𝜋𝜔0 ,那么首先将 𝑥(𝑡) 重新映射为周期 2𝜋 的函数,也就是只对 𝑥(𝑡) 做变换 𝑡→𝑡𝜔0 ,那么傅里叶级数还能照原样计算:

{𝑎𝑛=1𝜋∫−𝜋𝜋𝑥(𝑡𝜔0)cos𝑛𝑡𝑑𝑡𝑏𝑛=1𝜋∫−𝜋𝜋𝑥(𝑡𝜔0)sin𝑛𝑡𝑑𝑡

上式不甚好看,因此对其整体再重新做一次变换 𝑡→𝜔0𝑡 ,于是:

{𝑎𝑛=𝜔0𝜋∫−𝜋𝜔0𝜋𝜔0𝑥(𝑡)cos𝑛𝜔0𝑡𝑑𝑡𝑏𝑛=𝜔0𝜋∫−𝜋𝜔0𝜋𝜔0𝑥(𝑡)sin𝑛𝜔0𝑡𝑑𝑡

或者,由于 𝑇=2𝜋𝜔0 ,进一步化为:

{𝑎𝑛=2𝑇∫−𝑇2𝑇2𝑥(𝑡)cos𝑛𝜔0𝑡𝑑𝑡𝑏𝑛=2𝑇∫−𝑇2𝑇2𝑥(𝑡)sin𝑛𝜔0𝑡𝑑𝑡

上式是教科书上常给的形式。似乎既不好看,也不好记。

2.以复指数函数作为基底

在信号系统中,为了克服三角项级数既不好看也不好记的缺点,我们实际更常使用复指数函数作为基底。从三角项级数当然可以直接推出指数项级数的计算式,但这里不妨还是重新推导一下。现在我们直接取消掉 𝑥(𝑡) 的周期为 2𝜋 的限制,直接设置为 2𝑇 ,且 𝑇=2𝜋𝜔0 。

首先,函数的指数项展开不像三角项展开那样繁复,十分简洁:

𝑥(𝑡)=∑𝑘=−∞+∞𝑐𝑘ej𝑘𝜔0𝑡

这样复指数函数簇 [ej𝑘𝜔0𝑡],𝑘=0,±1,±2,… 的基波周期就是 𝑇=2𝜋𝜔0 。注意到:

其中 表示对 ej𝑛𝜔0𝑡 取共轭,实际就是简单地将 j 设为 −j 即可。所以:

∫𝑇𝑥(𝑡)e−j𝑛𝜔0𝑡𝑑𝑡=∑𝑘=−∞+∞𝑐𝑘∫𝑇ej𝑘𝜔0𝑡⋅e−j𝑛𝜔0𝑡𝑑𝑡=∑𝑘=−∞+∞𝑐𝑘∫𝑇ej(𝑘−𝑛)𝜔0𝑡𝑑𝑡=𝑐𝑛𝑇

因此:

𝑐𝑛=1𝑇∫𝑇𝑥(𝑡)e−j𝑛𝜔0𝑡𝑑𝑡

这与从三角项级数计算公式中推出的一样。同时,不像三角项级数中那样分别要算 𝑥(𝑡)cos𝑛𝑡,𝑥(𝑡)sin𝑛𝑡 的积分,这里一个式子即可涵盖所有情况。今后,傅里叶级数系数直指指数项级数 𝑐𝑛 而非三角项级数 𝑎𝑛,𝑏𝑛 。傅里叶级数系数也称为频谱系数,因为这些系数描述了每一个谐波分量的大小。其中,称 𝑐0 称为直流分量。

3.傅里叶级数收敛的条件

在假定信号 𝑥(𝑡) 能用复指数信号加权和 𝑥∗(𝑡)=∑𝑐𝑘ej𝑘𝜔0𝑡 拟合的情况下,频谱系数就是拟合效果最好的选择。但所谓“拟合效果最好”,并不是指能在任意 𝑡=𝑡0 时刻实现拟合函数的值无限趋近于实际值即 𝑥∗(𝑡0)→𝑥(𝑡0) ,而是指两者的误差函数 𝑒(𝑡)=𝑥(𝑡)−𝑥∗(𝑡) 在任意周期内的能量趋于零,即:

∫𝑇|𝑒(𝑡)|2𝑑𝑡=0

或者说,拟合函数和原信号在能量上并没有差距。这样,

定理 1.(傅里叶级数收敛的条件)
只要保证原信号在任意周期内的能量收敛到有限值,那么傅里叶级数一定是收敛的。 11 也就是只要保证 ∫𝑇|𝑥(𝑡)|2𝑑𝑡<∞ 即可。

以上是判定傅里叶级数收敛的其中一种方法。另一种方法是使用狄利克雷条件判断,而狄利克雷条件实际又包含三个条件、并要求这些条件同时成立:

定理 2.(狄利克雷条件)
  1. 在任意周期内,原信号的最大值和最小值的数量是有限的;
  2. 在任意区间内,原信号的不连续点的数量是有限的,并且在不连续点上的取值亦是有限的。
  3. 在任意周期内,原信号是绝对可积的,也就是:
∫𝑇|𝑥(𝑡)|𝑑𝑡<∞

4.吉布斯现象

吉布斯现象指的是,拟合函数会在原信号的不连续点附近出现高频的振荡,并且不论拟合函数的精度有多高,拟合函数的峰值相比原信号总有不超过约 9% 的超量。

注 2.(吉布斯现象的条件)

需要注意,出现吉布斯现象的条件是原信号存在不连续点,而不是存在不光滑点22 复习:连续的定义是左极限和右极限相等。例如,方波信号:

𝑥(𝑡)={−1, −𝑇2≤𝑡≤01, 0≤𝑡≤𝑇2

会在 𝑡=0 处出现吉布斯现象,而三角波信号:

𝑥(𝑡)=|𝑡|,−𝑇≤𝑡≤𝑇
则不会出现吉布斯现象。

5.帕塞瓦尔定理

定理 3.(帕塞瓦尔定理)

对于信号 𝑥(𝑡) 及其傅里叶级数系数而言,下式成立:

1𝑇∫𝑇|𝑥(𝑡)|2𝑑𝑡=∑𝑘=−∞+∞|𝑐𝑘|2

其意义是,在任意周期中,原信号的平均功率是所有谐波分量的平均功率之和。

注 3.

谐波分量的平均功率是:

𝑃𝑘=1𝑇∫𝑇|𝑐𝑘ej𝑘𝜔0𝑡|2𝑑𝑡=1𝑇∫𝑇|𝑐𝑘|2𝑑𝑡=|𝑐𝑘|2

三、傅里叶级数的性质

为了方便描述,以下用符号:

𝑥(𝑡)→ℱ︀𝒮︀𝑐𝑘

表示信号与其傅里叶级数系数的关系。

1.时间平移

假设信号 𝑥(𝑡) 的时间发生平移 𝑡→𝑡−𝑡0 ,那么其傅里叶级数可以重新计算如下:

𝑐𝑛′=1𝑇∫𝑇𝑥(𝑡−𝑡0)e−j𝑛𝜔0𝑡𝑑𝑡=1𝑇∫𝑇𝑥(𝑡)e−j𝑛𝜔0(𝑡+𝑡0)𝑑𝑡=e−j𝑛𝜔0𝑡01𝑇∫𝑇𝑥(𝑡)e−j𝑛𝜔0𝑡𝑑𝑡=e−j𝑛𝜔0𝑡0𝑐𝑛

因此,若 𝑥(𝑡)→ℱ︀𝒮︀𝑐𝑘 ,那么:

𝑥(𝑡−𝑡0)→ℱ︀𝒮︀e−j𝑛𝜔0𝑡0𝑐𝑘

2.时间反演

现在假设发生时间反演 𝑡→−𝑡 ,仿照上述计算:

𝑐𝑛′=1𝑇∫𝑇𝑥(−𝑡)e−j𝑛𝜔0𝑡𝑑𝑡=−1𝑇∫−𝑇𝑥(𝑡)ej𝑛𝜔0𝑡𝑑𝑡=1𝑇∫𝑇𝑥(𝑡)e−j(−𝑛)𝜔0𝑡𝑑𝑡=𝑐−𝑛

因此,若 𝑥(𝑡)→ℱ︀𝒮︀𝑐𝑘 ,那么:

𝑥(−𝑡)→ℱ︀𝒮︀𝑐−𝑘

由此可以引发对奇偶信号的讨论。若信号为奇信号即 𝑥(−𝑡)=−𝑥(𝑡) ,由上述结论显然可得 𝑐−𝑘=−𝑐𝑘 。而对于偶信号 𝑥(−𝑡)=𝑥(𝑡) ,则有 𝑐−𝑘=𝑐𝑘 。

3.时间伸缩

假设时间发生伸缩变换 𝑡→𝛼𝑡 ,其中要求 𝛼>0 ,那么根据:

𝑥(𝑡)=∑𝑘=−∞+∞𝑐𝑘ej𝑘𝜔0𝑡
⇒𝑥(𝛼𝑡)=∑𝑘=−∞+∞𝑐𝑘ej𝑘𝜔0⋅𝛼𝑡=∑𝑘=−∞+∞𝑐𝑘ej𝑘(𝛼𝜔0)𝑡

也就是说,傅里叶级数系数本身不发生改变,只是修改了谐波分量的基波频率,时间的伸缩倍率可以看做作用在基波频率上。

注 4.(由时间伸缩到时间反演)

当然另一方面,当 𝛼=−1 时此伸缩变换就可以看做时间反演,而此时这个取负操作视作作用在 𝑘 上:

𝑥(−𝑡)=∑𝑘=−∞+∞𝑐𝑘ej𝑘𝜔0−𝑡̇=∑𝑘=−∞+∞𝑐𝑘ej(−𝑘)𝜔0𝑡=∑𝑘=−∞+∞𝑐−𝑘ej𝑘𝜔0𝑡
这就为上一小节 “2.时间反演” 做了一个更简短的证明。同时也允许了 𝛼<0 的情况。基于此,之后对时间反演不再单独讨论,因为这可以直接体现在对时间伸缩变换中。

4.信号共轭

假如将信号取共轭即 ,那么傅里叶级数如何变换?还是通过处理下式:

𝑥(𝑡)=∑𝑘=−∞+∞𝑐𝑘ej𝑘𝜔0𝑡

两边同时取共轭可得:

因此:

根据这个结论,我们可以对实信号进行讨论。因为实信号 ,显然可得 ,由此可得以下四式:

{ℛ︀ℯ︀{𝑐𝑘}=ℛ︀ℯ︀{𝑐−𝑘}ℐ︀𝓂︀{𝑐𝑘}=−ℐ︀𝓂︀{𝑐−𝑘}|𝑐𝑘|=|𝑐−𝑘|arg(𝑐𝑘)=−arg(𝑐−𝑘)

5.信号相加(线性)

假设两个具有相同周期的信号33 务必注意周期相同的条件是非常重要的,在这一节说明傅里叶级数的性质时,取两个信号总默认其周期相同:

{𝑥(𝑡)→ℱ︀𝒮︀𝑎𝑘𝑦(𝑡)→ℱ︀𝒮︀𝑏𝑘

而 𝑧(𝑡)=𝛼𝑥(𝑡)+𝛽𝑦(𝑡) ,那么其傅里叶级数系数 𝑐𝑛 为:

𝑐𝑛=1𝑇∫𝑇(𝛼𝑥(𝑡)+𝛽𝑦(𝑡))e−j𝑛𝜔0𝑡𝑑𝑡=𝛼1𝑇∫𝑇𝑥(𝑡)e−j𝑛𝜔0𝑡𝑑𝑡+𝛽1𝑇∫𝑇𝑦(𝑡)e−j𝑛𝜔0𝑡𝑑𝑡=𝛼𝑎𝑛+𝛽𝑏𝑛

因此 𝑧(𝑡)→ℱ︀𝒮︀𝛼𝑎𝑘+𝛽𝑏𝑘 。

6.信号相乘(系数卷积)

同样还是取:

{𝑥(𝑡)→ℱ︀𝒮︀𝑎𝑘𝑦(𝑡)→ℱ︀𝒮︀𝑏𝑘

而 𝑧(𝑡)=𝑥(𝑡)𝑦(𝑡) 。这里我们不再直接计算 𝑐𝑛=1𝑇∫𝑇𝑥(𝑡)𝑦(𝑡)e−j𝑛𝜔0𝑡𝑑𝑡 ,而是换一种思路——既然存在:

{𝑥(𝑡)=∑𝑘=−∞+∞𝑎𝑘ej𝑘𝜔0𝑡𝑦(𝑡)=∑𝑙=−∞+∞𝑏𝑙ej𝑙𝜔0𝑡

所以:

𝑧(𝑡)=∑𝑘=−∞+∞𝑎𝑘ej𝑘𝜔0𝑡⋅∑𝑙=−∞+∞𝑏𝑙ej𝑙𝜔0𝑡=∑𝑘=−∞+∞∑𝑙=−∞+∞𝑎𝑘𝑏𝑙ej(𝑘+𝑙)𝜔0𝑡=∑𝑛=−∞+∞{∑𝑘=−∞+∞𝑎𝑘𝑏𝑛−𝑘}ej𝑛𝜔0𝑡

也就是说:

𝑐𝑛=∑𝑘=−∞+∞𝑎𝑘𝑏𝑛−𝑘

注意到上式其实就是一个卷积,因此也可以简写为 𝑐𝑛=𝑎𝑛∗𝑏𝑛 ——但需要注意,这样的写法实际将 𝑎𝑛,𝑏𝑛 看做一个离散序列,而非孤立的值。

7.信号卷积(系数相乘)

设 𝑧(𝑡) 是下述两信号在一个周期上的卷积:

{𝑥(𝑡)→ℱ︀𝒮︀𝑎𝑘𝑦(𝑡)→ℱ︀𝒮︀𝑏𝑘

即 𝑧(𝑡)=∫𝑇𝑥(𝜏)𝑦(𝑡−𝜏)𝑑𝜏 ,代入 𝑐𝑛=1𝑇∫𝑇𝑧(𝑡)e−j𝑛𝜔0𝑡𝑑𝑡 可得:

𝑐𝑛=1𝑇∫𝑇{∫𝑇𝑥(𝜏)𝑦(𝑡−𝜏)𝑑𝜏}e−j𝑛𝜔0𝑡𝑑𝑡=1𝑇∫𝑇∫𝑇𝑥(𝜏)𝑦(𝑡−𝜏)e−j𝑛𝜔0𝑡𝑑𝜏𝑑𝑡=1𝑇∫𝑇𝑥(𝜏)𝑑𝜏∫𝑇𝑦(𝑡−𝜏)e−j𝑛𝜔0𝑡𝑑𝑡

注意到后半部分积分其实就是 𝑦(𝑡) 的傅里叶级数系数计算式,只是有 𝜏 的时间平移、并相差一个 𝑇 的倍数关系。故利用上面已经得到的时间平移结论:

𝑐𝑛=𝑇𝑏𝑛𝑇∫𝑇𝑥(𝜏)e−j𝑛𝜔0𝜏𝑑𝜏=𝑇2𝑎𝑛𝑏𝑛𝑇=𝑇𝑎𝑛𝑏𝑛

又恰好可以组成 𝑥(𝑡) 的傅里叶级数系数计算式!从而得到以下结论:

∫𝑇𝑥(𝜏)𝑦(𝑡−𝜏)𝑑𝜏→ℱ︀𝒮︀𝑇𝑎𝑛𝑏𝑛

8.信号微分

若对信号 𝑥(𝑡)→ℱ︀𝒮︀𝑐𝑘 求导,那么根据:

𝑥(𝑡)=∑𝑘=−∞+∞𝑐𝑘ej𝑘𝜔0𝑡

两边求导,立即可得:

𝑑𝑥(𝑡)𝑑𝑡=∑𝑘=−∞+∞𝑗𝑘𝜔0𝑐𝑘ej𝑘𝜔0𝑡

因此:

𝑥(𝑡)→ℱ︀𝒮︀𝑐𝑘⇒𝑑𝑥(𝑡)𝑑𝑡→ℱ︀𝒮︀𝑗𝑘𝜔0𝑐𝑘

9.信号积分

同样,若取 𝑧(𝑡)=∫−∞𝑡𝑥(𝑡)𝑑𝑡 ,按同样的计算:

∫−∞𝑡𝑥(𝑡)𝑑𝑡=∑𝑘=−∞+∞𝑐𝑘𝑗𝑘𝜔0ej𝑘𝜔0𝑡

故:

𝑥(𝑡)→ℱ︀𝒮︀𝑐𝑘⇒∫−∞𝑡𝑥(𝑡)𝑑𝑡→ℱ︀𝒮︀𝑐𝑘𝑗𝑘𝜔0

10.信号奇偶

这里我们要求信号 𝑥(𝑡) 是实信号,这样才能满足 的条件。将信号分解为偶部和奇部:

{ℰ︀𝓋︀{𝑥(𝑡)}=12(𝑥(𝑡)+𝑥(−𝑡))𝒪︀𝒹︀{𝑥(𝑡)}=12(𝑥(𝑡)−𝑥(−𝑡))

于是利用线性性质,对于偶部:

ℰ︀𝓋︀{𝑥(𝑡)}→ℱ︀𝒮︀12(𝑐𝑘+𝑐−𝑘)=ℛ︀ℯ︀{𝑐𝑘}

对于奇部:

𝒪︀𝒹︀{𝑥(𝑡)}→ℱ︀𝒮︀12(𝑐𝑘−𝑐−𝑘)=jℐ︀𝓂︀{𝑐𝑘}

四、傅里叶变换的计算

以上内容中,为了将信号分解为复指数信号的加权和,我们一直要求待分解信号是周期信号。并且我们也能发现,信号被分解之后,各分量之间是谐波关系,也即频率的取值是离散的、分立的、成倍的。现在我们不禁想问:非周期信号,是否也能展开为“傅里叶级数”?

1.非周期信号的傅里叶变换(狭义傅里叶变换)

答案是肯定的,但此时并不能说是将其展开为傅里叶级数,而是做傅里叶变换。注意到傅里叶级数系数的计算公式:

𝑐𝑛=1𝑇∫𝑇𝑥(𝑡)e−j𝑛𝜔0𝑡𝑑𝑡

于是下式恒成立:

𝑥(𝑡)=1𝑇∑𝑘=−∞+∞∫𝑇𝑥(𝑡)e−j𝑘𝜔0𝑡𝑑𝑡⋅ej𝑘𝜔0𝑡=12𝜋∑𝑘=−∞+∞∫𝑇𝑥(𝑡)e−j𝑘𝜔0𝑡𝑑𝑡⋅ej𝑘𝜔0𝑡⋅𝜔0

其中 2𝑇 是信号的周期。非周期信号为什么不可以看做周期为无穷大的周期信号呢?在此思想下,显然有 𝑇→∞,𝜔0→0 ,此时 𝑘𝜔0 可以视作连续变量 𝜔 ,并将 𝜔0 重新记为 𝑑𝜔 :

𝑥(𝑡)=12𝜋∑𝑘=−∞+∞∫−∞+∞𝑥(𝑡)e−j𝜔𝑡𝑑𝑡⋅ej⋅𝑘𝑑𝜔⋅𝑡𝑑𝜔

其中 ∫−∞+∞𝑥(𝑡)e−j𝜔𝑡𝑑𝑡 的结果显然是关于新连续变量 𝜔 的函数,或者一般认为是 j𝜔 的函数,记为 𝑋(j𝜔) 或 ℱ︀{𝑥(𝑡)} ,称为频谱函数。于是上式又能化为黎曼积分:

定义 1.(非周期信号的傅里叶逆变换式/综合公式)
𝑥(𝑡)=12𝜋∫−∞+∞𝑋(j𝜔)⋅ej𝜔𝑡𝑑𝜔
该式可记为 ℱ︀−1{𝑋(j𝜔)} 。

这样我们就将非周期信号 𝑥(𝑡) 也表示成了复指数函数的加权和——但由于非周期信号的频谱连续,所以实际应写成积分的形式。必须要指出,虽然上式与下式44 一个是对非周期信号、一个是对周期信号:

𝑥(𝑡)=∑𝑘=−∞+∞𝑐𝑘ej𝑘𝜔0𝑡

看起来大有不同,但本质是相同的,无非一个连续、一个离散而已。而傅里叶级数系数计算式:

𝑐𝑛=1𝑇∫𝑇𝑥(𝑡)e−j𝑛𝜔0𝑡𝑑𝑡

和傅里叶变换式:

定义 2.(傅里叶变换式)
𝑋(j𝜔)=ℱ︀{𝑥(𝑡)}=∫−∞+∞𝑥(𝑡)e−j𝜔𝑡𝑑𝑡

虽然不能划等号,但本质也是相同的。其中的区别会在下一小节讲述。

2.周期信号的傅里叶变换(广义傅里叶变换)

上一小节中说明了对于非周期信号来说如何计算其傅里叶变换,但同样的计算方法却无法适用于周期信号。由于前半部分已经说明周期信号可以表示为成谐波的复指数信号的加权和,现在为了说明上一小节导出的计算方法不能适用于周期信号,只需分析基本复指数信号 𝑥(𝑡)=ej𝜔0𝑡 即可。对其傅里叶变换:

𝑋(j𝜔)=∫−∞+∞ej𝜔0𝑡⋅e−j𝜔𝑡𝑑𝑡=∫−∞+∞ej(𝜔0−𝜔)𝑡𝑑𝑡

这个积分无论如何都不收敛!这似乎导向了一个令人沮丧的结论,即所有周期信号都不能傅里叶变换;甚至含有周期信号的非周期信号也不能傅里叶变换。幸运的是,如果对奇异函数 𝛿(𝑡) 即单位冲激函数做傅里叶变换,可以得到如下结果:

ℱ︀{𝛿(𝑡)}=∫−∞+∞𝛿(𝑡)⋅e−j𝜔𝑡𝑑𝑡=∫−∞+∞𝛿(𝑡)⋅e0𝑑𝑡=1

按照这个变换式来定义原本不收敛的傅里叶逆变换式:

ℱ︀−1{1}=𝛿(𝑡)=12𝜋∫−∞+∞1⋅ej𝜔𝑡𝑑𝜔

也就是我们定义 ∫−∞+∞ej𝜔𝑡𝑑𝜔=2𝜋𝛿(𝑡) ,而这等价于 ∫−∞+∞e−j𝜔𝑡𝑑𝑡=2𝜋𝛿(𝜔) ,相当于又定义了如下傅里叶变换:

ℱ︀{1}=2𝜋𝛿(𝜔)

这样我们就能利用此变换并结合傅里叶变换的性质,求出任意周期信号的傅里叶变换了。假定周期信号 𝑥(𝑡) 具有傅里叶级数形式:

𝑥(𝑡)=∑𝑘=−∞+∞𝑐𝑘ej𝑘𝜔0𝑡

那么其傅里叶变换就是:

定理 4.(周期信号的傅里叶变换)
𝑋(j𝜔)=2𝜋∑𝑘=−∞+∞𝑐𝑘𝛿(𝜔−𝑘𝜔0)

也就是说,周期函数的频谱函数实际就是按其傅里叶级数系数排列的冲激串,但在强度上还需要放大 2𝜋 倍。

注 5.
在未定义 1 的傅里叶变换时、无法求出周期信号的频谱函数 𝑋(j𝜔) 的本质原因,就是因为周期信号的频谱是离散的、而非周期信号的频谱则是连续的。在未引入单位冲激函数时,连续函数 𝑋(j𝜔) 天然地不具备表达离散函数的能力,直到我们定义 ℱ︀{1}=2𝜋𝛿(𝜔) 从而引入单位冲激函数,便可以解决这个问题。

五、傅里叶变换的性质

由于傅里叶变换式:

𝑋(j𝜔)=ℱ︀{𝑥(𝑡)}=∫−∞+∞𝑥(𝑡)e−j𝜔𝑡𝑑𝑡

和傅里叶级数系数计算式:

𝑐𝑛=1𝑇∫𝑇𝑥(𝑡)e−j𝑛𝜔0𝑡𝑑𝑡

之间相差并不多,只是一个系数 1𝑇 、参数 𝜔,𝑛𝜔0 以及积分上下限不同,因此傅里叶级数系数的性质基本都可以挪用于此处。仿照之前傅里叶级数部分的写法,下述符号:

𝑥(𝑡)→ℱ︀𝑋(j𝜔)

表示原信号 𝑥(𝑡) 与其频谱函数 𝑋(j𝜔) 的对应关系。

1.时间平移

当发生时间平移 𝑡→𝑡−𝑡0 :

𝑋′(j𝜔)=∫−∞+∞𝑥(𝑡−𝑡0)e−j𝜔𝑡𝑑𝑡=∫−∞+∞𝑥(𝑡)e−j𝜔(𝑡+𝑡0)𝑑𝑡=e−j𝜔𝑡0∫−∞+∞𝑥(𝑡)e−j𝜔𝑡𝑑𝑡=e−j𝜔𝑡0𝑋(j𝜔)

因此:

𝑥(𝑡)→ℱ︀𝑋(j𝜔)⇒𝑥(𝑡−𝑡0)→ℱ︀e−j𝜔𝑡0𝑋(j𝜔)

2.时间伸缩(频域伸缩)

当发生时间伸缩 𝑡→𝛼𝑡 :

𝑋′(j𝜔)=∫−∞+∞𝑥(𝛼𝑡)e−j𝜔𝑡𝑑𝑡=1|𝛼|∫−∞+∞𝑥(𝑡)e−j𝜔𝑡𝛼𝑑𝑡=1|𝛼|∫−∞+∞𝑥(𝑡)e−j𝜔𝛼𝑡𝑑𝑡=1|𝛼|𝑋(j𝜔𝛼)

因此:

𝑥(𝑡)→ℱ︀𝑋(j𝜔)⇒𝑥(𝛼𝑡)→ℱ︀1|𝛼|𝑋(j𝜔𝛼)⇒𝑥(−𝑡)→ℱ︀𝑋(−j𝜔)

3.频域平移

若频域发生平移 𝜔→𝜔−𝜔0 ,那么可以计算傅里叶逆变换:

𝑥′(𝑡)=12𝜋∫−∞+∞𝑋(𝑗(𝜔−𝜔0))ej𝜔𝑡𝑑𝜔=12𝜋∫−∞+∞𝑋(j𝜔)ej(𝜔+𝜔0)𝑡𝑑𝜔=ej𝜔0𝑡12𝜋∫−∞+∞𝑋(j𝜔)ej𝜔𝑡𝑑𝜔=ej𝜔0𝑡𝑥(𝑡)

因此:

𝑥(𝑡)→ℱ︀𝑋(j𝜔)⇒ej𝜔0𝑡𝑥(𝑡)→ℱ︀𝑋(𝑗(𝜔−𝜔0))

4.信号共轭

考虑取信号的共轭 ,那么将傅里叶变换式两边同时取共轭:

因此 ,故:

根据这个结论,我们又可以展开对实信号的讨论。类似在傅里叶级数中的结论,实信号首先有共轭对称性,即 ,于是:

{ℛ︀ℯ︀{𝑋(j𝜔)}=ℛ︀ℯ︀{𝑋(−j𝜔)}ℐ︀𝓂︀{𝑋(j𝜔)}=−ℐ︀𝓂︀{𝑋(−j𝜔)}

5.信号相加(线性)

对于 𝑧(𝑡)=𝛼𝑥(𝑡)+𝛽𝑦(𝑡) ,显然:

𝑍(j𝜔)=∫−∞+∞(𝛼𝑥(𝑡)+𝛽𝑦(𝑡))e−j𝜔𝑡𝑑𝑡=𝛼∫−∞+∞𝑥(𝑡)e−j𝜔𝑡𝑑𝑡+𝛽∫−∞+∞𝑦(𝑡)e−j𝜔𝑡𝑑𝑡=𝛼𝑋(j𝜔)+𝛽𝑌(j𝜔)

故:

𝛼𝑥(𝑡)+𝛽𝑦(𝑡)→ℱ︀𝛼𝑋(j𝜔)+𝛽𝑌(j𝜔)

值得注意的是,相比傅里叶级数部分,这里取消了 𝑥(𝑡),𝑦(𝑡) 应有相同周期的限制——毕竟傅里叶变换的作用对象本就是非周期函数。

6.信号微分

将傅里叶逆变换式两边求导可得:

𝑥(𝑡)=12𝜋∫−∞+∞𝑋(j𝜔)ej𝜔𝑡𝑑𝜔
⇒𝑑𝑥(𝑡)𝑑𝑡=12𝜋𝑑𝑑𝑡∫−∞+∞𝑋(j𝜔)ej𝜔𝑡𝑑𝜔=12𝜋∫−∞+∞𝑋(j𝜔)⋅𝑑(ej𝜔𝑡)𝑑𝑡𝑑𝜔=12𝜋∫−∞+∞j𝜔𝑋(j𝜔)⋅ej𝜔𝑡𝑑𝜔

故:

𝑑𝑥(𝑡)𝑑𝑡→ℱ︀j𝜔𝑋(j𝜔)

7.信号积分

将傅里叶逆变换式两边求积可得:

∫−∞𝑡𝑥(𝑡)𝑑𝑡=12𝜋∫−∞𝑡∫−∞+∞𝑋(j𝜔)ej𝜔𝑡𝑑𝜔𝑑𝑡=12𝜋∫−∞𝑡ej𝜔𝑡𝑑𝑡∫−∞+∞𝑋(j𝜔)𝑑𝜔=12𝜋∫−∞+∞1j𝜔𝑋(j𝜔)ej𝜔𝑡𝑑𝜔

故:

∫−∞𝑡𝑥(𝑡)𝑑𝑡→ℱ︀1j𝜔𝑋(j𝜔)
注 6.(与教科书的区别)

假设已知 𝑥(𝑡)→ℱ︀𝑋(j𝜔) ,并且 𝑥(𝑡)=𝑦′(𝑡) ,现在我们尝试求 𝑦(𝑡) 的傅里叶变换。在这里值得注意的是:由于求导会使常数项归零,因此我们应当假设 𝑦(𝑡) 中存在常数项;即,仅认为 𝑦(𝑡)=∫−∞𝑡𝑥(𝑡)𝑑𝑡 是不全面的,而应是 𝑦(𝑡)=∫−∞𝑡𝑥(𝑡)𝑑𝑡+𝐶 。后面这个常数项会带来 2𝜋𝐶𝛿(𝜔) 的冲激。但因为前面一项即 1j𝜔𝑋(j𝜔) 已经使 𝜔=0 的取值成为不可能,因此这个冲激可以略去不计。

当然,在奥本海姆的教科书上给出的傅里叶变换积分性质是:

∫−∞𝑡𝑥(𝑡)𝑑𝑡→ℱ︀1j𝜔𝑋(j𝜔)+𝜋𝑋(0)𝛿(𝜔)

这大抵是利用了 ∫−∞𝑡𝑥(𝑡)𝑑𝑡=𝑥(𝑡)∗𝑢(𝑡) 的性质之后,结合 ℱ︀{𝑢(𝑡)}=1j𝜔+𝜋𝛿(𝜔) 和之后要讲的卷积性质导出的。

8.信号奇偶

类似于傅里叶级数中的讨论,我们还是要求信号为实信号以符合共轭对称性 。这样结合线性性质可得以下结论:

{ℰ︀𝓋︀{𝑥(𝑡)}→ℱ︀ℛ︀ℯ︀{𝑋(j𝜔)}𝒪︀𝒹︀{𝑥(𝑡)}→ℱ︀jℐ︀𝓂︀{𝑋(j𝜔)}

9.时域相乘(频域卷积)

傅里叶变换极其重要的性质之一,就是原信号相乘,那么其频谱函数则是原先两个频谱函数的卷积。设 𝑧(𝑡)=𝑥(𝑡)𝑦(𝑡) ,那么:

𝑧(𝑡)=14𝜋2∫−∞+∞𝑋(j𝜑)ej𝜑𝑡𝑑𝜑⋅∫−∞+∞𝑌(j𝜓)ej𝜓𝑡𝑑𝜓=14𝜋2∫−∞+∞∫−∞+∞𝑋(j𝜑)𝑌(j𝜓)ej(𝜑+𝜓)𝑡𝑑𝜑𝑑𝜓

现在令 𝜔=𝜑+𝜓 ,所以:

𝑧(𝑡)=14𝜋2∫−∞+∞∫−∞+∞𝑋(j𝜑)𝑌(𝑗(𝜔−𝜑))ej𝜔𝑡𝑑𝜑𝑑𝜔=12𝜋∫−∞+∞{12𝜋∫−∞+∞𝑋(j𝜑)𝑌(j(𝜔−𝜑))𝑑𝜑}ej𝜔𝑡𝑑𝜔

从而可以看出 𝑍(j𝜔)=12𝜋∫−∞+∞𝑋(j𝜑)𝑌(j(𝜔−𝜑))𝑑𝜑=12𝜋𝑋(j𝜑)∗𝑌(j𝜑) 。故:

𝑥(𝑡)𝑦(𝑡)→ℱ︀12𝜋𝑋(j𝜔)∗𝑌(j𝜔)

10.时域卷积(频域相乘)

傅里叶变换的另一重要性质,就是原信号卷积,那么其频谱函数则是原先两个频谱函数的乘积。设 𝑧(𝑡)=𝑥(𝑡)∗𝑦(𝑡) ,那么:

𝑍(j𝜔)=∫−∞+∞{∫−∞+∞𝑥(𝜏)𝑦(𝑡−𝜏)𝑑𝜏}e−j𝜔𝑡𝑑𝑡=∫−∞+∞𝑥(𝜏)𝑑𝜏∫−∞+∞𝑦(𝑡−𝜏)e−j𝜔𝑡𝑑𝑡=𝑌(j𝜔)∫−∞+∞𝑥(𝜏)e−j𝜔𝜏𝑑𝜏=𝑋(j𝜔)𝑌(j𝜔)

故:

𝑥(𝑡)∗𝑦(𝑡)→ℱ︀𝑋(j𝜔)𝑌(j𝜔)

11.频域微分

对傅里叶变换式两边求导可得:

𝑋(j𝜔)=∫−∞+∞𝑥(𝑡)e−j𝜔𝑡𝑑𝑡
⇒𝑑𝑋(j𝜔)𝑑𝜔=𝑑𝑑𝜔∫−∞+∞𝑥(𝑡)e−j𝜔𝑡𝑑𝑡=∫−∞+∞𝑥(𝑡)𝑑(e−j𝜔𝑡)𝑑𝜔𝑑𝑡=∫−∞+∞−j𝑡𝑥(𝑡)e−j𝜔𝑡𝑑𝑡

故:

−j𝑡𝑥(𝑡)→ℱ︀𝑑𝑋(j𝜔)𝑑𝜔

当然,由于 −1j=j ,上式我们更一般写为:

𝑡𝑥(𝑡)→ℱ︀j𝑑𝑋(j𝜔)𝑑𝜔

12.对偶性

对偶性指的是,傅里叶变换式和傅里叶逆变换式在形式上有极大的相似性:

{𝑋(j𝜔)=ℱ︀{𝑥(𝑡)}=∫−∞+∞𝑥(𝑡)e−j𝜔𝑡𝑑𝑡𝑥(𝑡)=ℱ︀−1{𝑋(j𝜔)}=12𝜋∫−∞+∞𝑋(j𝜔)⋅ej𝜔𝑡𝑑𝜔

这导致了:如果信号 𝑥(𝑡) 在傅里叶变换之后得到频谱函数 𝑋(j𝜔) ,再对频谱函数做变量替换、记 𝑦(𝑡)=𝑋(𝑗𝑡) ,那么 𝑌(j𝜔) 在形式上又会与 𝑥(𝑡) 有极大的相似性!典型的例子是 1 与单位冲激函数 𝛿(𝑡) 的傅里叶变换对,可以参考这个例子琢磨一下“对偶性”究竟为何物:

{1→ℱ︀2𝜋𝛿(𝜔)𝛿(𝑡)→ℱ︀1

当然您可能也已经注意到了:譬如时域相乘相当于频域卷积,反过来也成立——时域卷积相当于频域相乘;时域微分和频域微分导出的结果竟也差不多……这也是对偶性的体现之一。


六、基本傅里叶变换表

以下给出一些基本信号的傅里叶变换表:

表 1 基本傅里叶变换表

信号𝑥(𝑡)傅里叶变换𝑋(𝑗𝜔)=ℱ︀{𝑥(𝑡)}条件或说明
12𝜋𝛿(𝜔)—
ej𝑘𝜔0𝑡2𝜋𝛿(𝜔−𝑘𝜔0)—
∑𝑘=−∞+∞𝑎𝑘ej𝑘𝜔0𝑡2𝜋∑𝑘=−∞+∞𝑎𝑘𝛿(𝜔−𝑘𝜔0)—
𝛿(𝑡)1—
𝛿(𝑡−𝑡0)e−j𝜔𝑡0—
𝑢(𝑡)1j𝜔+𝜋𝛿(𝜔)—
∑𝑘=−∞+∞𝛿(𝑡−𝑘𝑇)2𝜋𝑇∑𝑘=−∞+∞𝛿(𝜔−2𝑘𝜋𝑇)—
cos(𝜔0𝑡)𝜋(𝛿(𝜔−𝜔0)+𝛿(𝜔+𝜔0))—
sin(𝜔0𝑡)𝜋j(𝛿(𝜔−𝜔0)−𝛿(𝜔+𝜔0))—
e−𝑎𝑡𝑢(𝑡)1𝑎+j𝜔ℛ︀ℯ︀{𝑎}>0
𝑡e−𝑎𝑡𝑢(𝑡)1(𝑎+j𝜔)2ℛ︀ℯ︀{𝑎}>0
𝑡𝑛−1(𝑛−1)!e−𝑎𝑡𝑢(𝑡)1(𝑎+j𝜔)𝑛ℛ︀ℯ︀{𝑎}>0
sin(𝑊𝑡)𝜋𝑡{1, |𝜔|<𝑊0, |𝜔|>𝑊—
{1, |𝑡|<𝑇10, |𝑡|>𝑇12sin(𝜔𝑇1)𝜔—
{1, |𝑡|<𝑇10, 𝑇1<|𝑡|≤𝑇2𝑥(𝑡+𝑇)=𝑥(𝑡)∑𝑘=−∞+∞2sin(𝑘𝜔0𝑇1)𝑘𝛿(𝜔−𝑘𝜔0)—