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

将函数分解为基本函数的加权和,是研究函数十分重要的方法之一。比如上一篇中我们将函数分解为了单位冲激函数的加权和(对于离散时间信号是这样,连续信号则是卷积积分),就是将单位冲激函数作为了基本函数。再例如,如果将幂函数作为基本函数,那就是泰勒展开;而以三角函数作为基本函数,那就是本篇的重点——傅里叶展开了。
一、复指数信号的特殊性
为什么复指数函数可以作为基本函数呢?这主要是因为复指数信号在经过线性时不变系统输出后,只会产生幅值上的变化。也就是说,假设复指数信号 ,其中 是任意复数,那么:
幅值会产生 倍的变化,而且这个幅值还由 唯一确定!我们将这种经过系统作用之后的输出只是输入的常数倍的信号称为特征函数,而这个常数称为特征值。
证明. 为了证明以上结论,首先假设系统的单位冲激响应为 ,。那么系统对 的响应就是:
因此,若信号可以被分解为复指数信号的加权和:
显然输出可以直接确定为:
这就为研究系统带来了极大的便利,正如上篇文章中那般,只要确定单位冲激响应就能完全确定线性时不变系统对任意信号的响应,现在只要确定 也能完全确定线性时不变系统对某些信号的响应。
二、傅里叶级数的计算
在看到将函数分解为复指数函数带来便利之后,我们的下一个问题是:如何将一个函数分解为复指数函数的加权和?首先需要注意到复指数函数是周期信号,因此“只有”周期函数可以被展开。以下内容都事先默认假定要展开的函数是周期函数。
1.以三角函数作为基底
在高等数学中我们学习过傅里叶级数、三角展开,就是将一个函数变换为另一些函数的线性组合——而如果将函数视为向量的话,“另一些函数”的函数簇中的函数两两正交,那么就是将原函数在这些函数张成的无限维空间中正交分解罢了。傅里叶展开所选的正交基是三角函数系 ,这其中任取两不同函数,其积在 上的 范数都等于 ,也即两两正交:
注意到我们要求三角函数系中 为整数,这就要求 必须是这些三角函数的周期。这只是为了后续书写计算方便,之后会逐步放开限制。假设某周期为 的函数 可以傅里叶展开为:
直接对上式两边在 区间上积分,注意到上面提到的 ,那么:
因此得到 的计算式:
类似地,将 乘上 后做同样的操作,并注意三角函数系的性质:
于是:
将 乘上 后做同样的操作:
于是:
综上,我们得到了周期为 的函数 的傅里叶三角项级数的计算式:
那么,若 的周期不为 呢?不必重新计算,因为只是相当于做了时间伸缩变换。假设周期为 ,那么首先将 重新映射为周期 的函数,也就是只对 做变换 ,那么傅里叶级数还能照原样计算:
上式不甚好看,因此对其整体再重新做一次变换 ,于是:
或者,由于 ,进一步化为:
上式是教科书上常给的形式。似乎既不好看,也不好记。
2.以复指数函数作为基底
在信号系统中,为了克服三角项级数既不好看也不好记的缺点,我们实际更常使用复指数函数作为基底。从三角项级数当然可以直接推出指数项级数的计算式,但这里不妨还是重新推导一下。现在我们直接取消掉 的周期为 的限制,直接设置为 ,且 。
首先,函数的指数项展开不像三角项展开那样繁复,十分简洁:
这样复指数函数簇 的基波周期就是 。注意到:
其中 表示对 取共轭,实际就是简单地将 设为 即可。所以:
因此:
这与从三角项级数计算公式中推出的一样。同时,不像三角项级数中那样分别要算 的积分,这里一个式子即可涵盖所有情况。今后,傅里叶级数系数直指指数项级数 而非三角项级数 。傅里叶级数系数也称为频谱系数,因为这些系数描述了每一个谐波分量的大小。其中,称 称为直流分量。
3.傅里叶级数收敛的条件
在假定信号 能用复指数信号加权和 拟合的情况下,频谱系数就是拟合效果最好的选择。但所谓“拟合效果最好”,并不是指能在任意 时刻实现拟合函数的值无限趋近于实际值即 ,而是指两者的误差函数 在任意周期内的能量趋于零,即:
或者说,拟合函数和原信号在能量上并没有差距。这样,
以上是判定傅里叶级数收敛的其中一种方法。另一种方法是使用狄利克雷条件判断,而狄利克雷条件实际又包含三个条件、并要求这些条件同时成立:
- 在任意周期内,原信号的最大值和最小值的数量是有限的;
- 在任意区间内,原信号的不连续点的数量是有限的,并且在不连续点上的取值亦是有限的。
- 在任意周期内,原信号是绝对可积的,也就是:
4.吉布斯现象
吉布斯现象指的是,拟合函数会在原信号的不连续点附近出现高频的振荡,并且不论拟合函数的精度有多高,拟合函数的峰值相比原信号总有不超过约 的超量。
5.帕塞瓦尔定理
对于信号 及其傅里叶级数系数而言,下式成立:
其意义是,在任意周期中,原信号的平均功率是所有谐波分量的平均功率之和。
谐波分量的平均功率是:
三、傅里叶级数的性质
为了方便描述,以下用符号:
表示信号与其傅里叶级数系数的关系。
1.时间平移
假设信号 的时间发生平移 ,那么其傅里叶级数可以重新计算如下:
因此,若 ,那么:
2.时间反演
现在假设发生时间反演 ,仿照上述计算:
因此,若 ,那么:
由此可以引发对奇偶信号的讨论。若信号为奇信号即 ,由上述结论显然可得 。而对于偶信号 ,则有 。
3.时间伸缩
假设时间发生伸缩变换 ,其中要求 ,那么根据:
也就是说,傅里叶级数系数本身不发生改变,只是修改了谐波分量的基波频率,时间的伸缩倍率可以看做作用在基波频率上。
当然另一方面,当 时此伸缩变换就可以看做时间反演,而此时这个取负操作视作作用在 上:
4.信号共轭
假如将信号取共轭即 ,那么傅里叶级数如何变换?还是通过处理下式:
两边同时取共轭可得:
因此:
根据这个结论,我们可以对实信号进行讨论。因为实信号 ,显然可得 ,由此可得以下四式:
5.信号相加(线性)
假设两个具有相同周期的信号33 务必注意周期相同的条件是非常重要的,在这一节说明傅里叶级数的性质时,取两个信号总默认其周期相同:
而 ,那么其傅里叶级数系数 为:
因此 。
6.信号相乘(系数卷积)
同样还是取:
而 。这里我们不再直接计算 ,而是换一种思路——既然存在:
所以:
也就是说:
注意到上式其实就是一个卷积,因此也可以简写为 ——但需要注意,这样的写法实际将 看做一个离散序列,而非孤立的值。
7.信号卷积(系数相乘)
设 是下述两信号在一个周期上的卷积:
即 ,代入 可得:
注意到后半部分积分其实就是 的傅里叶级数系数计算式,只是有 的时间平移、并相差一个 的倍数关系。故利用上面已经得到的时间平移结论:
又恰好可以组成 的傅里叶级数系数计算式!从而得到以下结论:
8.信号微分
若对信号 求导,那么根据:
两边求导,立即可得:
因此:
9.信号积分
同样,若取 ,按同样的计算:
故:
10.信号奇偶
这里我们要求信号 是实信号,这样才能满足 的条件。将信号分解为偶部和奇部:
于是利用线性性质,对于偶部:
对于奇部:
四、傅里叶变换的计算
以上内容中,为了将信号分解为复指数信号的加权和,我们一直要求待分解信号是周期信号。并且我们也能发现,信号被分解之后,各分量之间是谐波关系,也即频率的取值是离散的、分立的、成倍的。现在我们不禁想问:非周期信号,是否也能展开为“傅里叶级数”?
1.非周期信号的傅里叶变换(狭义傅里叶变换)
答案是肯定的,但此时并不能说是将其展开为傅里叶级数,而是做傅里叶变换。注意到傅里叶级数系数的计算公式:
于是下式恒成立:
其中 是信号的周期。非周期信号为什么不可以看做周期为无穷大的周期信号呢?在此思想下,显然有 ,此时 可以视作连续变量 ,并将 重新记为 :
其中 的结果显然是关于新连续变量 的函数,或者一般认为是 的函数,记为 或 ,称为频谱函数。于是上式又能化为黎曼积分:
这样我们就将非周期信号 也表示成了复指数函数的加权和——但由于非周期信号的频谱连续,所以实际应写成积分的形式。必须要指出,虽然上式与下式44 一个是对非周期信号、一个是对周期信号:
看起来大有不同,但本质是相同的,无非一个连续、一个离散而已。而傅里叶级数系数计算式:
和傅里叶变换式:
虽然不能划等号,但本质也是相同的。其中的区别会在下一小节讲述。
2.周期信号的傅里叶变换(广义傅里叶变换)
上一小节中说明了对于非周期信号来说如何计算其傅里叶变换,但同样的计算方法却无法适用于周期信号。由于前半部分已经说明周期信号可以表示为成谐波的复指数信号的加权和,现在为了说明上一小节导出的计算方法不能适用于周期信号,只需分析基本复指数信号 即可。对其傅里叶变换:
这个积分无论如何都不收敛!这似乎导向了一个令人沮丧的结论,即所有周期信号都不能傅里叶变换;甚至含有周期信号的非周期信号也不能傅里叶变换。幸运的是,如果对奇异函数 即单位冲激函数做傅里叶变换,可以得到如下结果:
按照这个变换式来定义原本不收敛的傅里叶逆变换式:
也就是我们定义 ,而这等价于 ,相当于又定义了如下傅里叶变换:
这样我们就能利用此变换并结合傅里叶变换的性质,求出任意周期信号的傅里叶变换了。假定周期信号 具有傅里叶级数形式:
那么其傅里叶变换就是:
也就是说,周期函数的频谱函数实际就是按其傅里叶级数系数排列的冲激串,但在强度上还需要放大 倍。
五、傅里叶变换的性质
由于傅里叶变换式:
和傅里叶级数系数计算式:
之间相差并不多,只是一个系数 、参数 以及积分上下限不同,因此傅里叶级数系数的性质基本都可以挪用于此处。仿照之前傅里叶级数部分的写法,下述符号:
表示原信号 与其频谱函数 的对应关系。
1.时间平移
当发生时间平移 :
因此:
2.时间伸缩(频域伸缩)
当发生时间伸缩 :
因此:
3.频域平移
若频域发生平移 ,那么可以计算傅里叶逆变换:
因此:
4.信号共轭
考虑取信号的共轭 ,那么将傅里叶变换式两边同时取共轭:
因此 ,故:
根据这个结论,我们又可以展开对实信号的讨论。类似在傅里叶级数中的结论,实信号首先有共轭对称性,即 ,于是:
5.信号相加(线性)
对于 ,显然:
故:
值得注意的是,相比傅里叶级数部分,这里取消了 应有相同周期的限制——毕竟傅里叶变换的作用对象本就是非周期函数。
6.信号微分
将傅里叶逆变换式两边求导可得:
故:
7.信号积分
将傅里叶逆变换式两边求积可得:
故:
假设已知 ,并且 ,现在我们尝试求 的傅里叶变换。在这里值得注意的是:由于求导会使常数项归零,因此我们应当假设 中存在常数项;即,仅认为 是不全面的,而应是 。后面这个常数项会带来 的冲激。但因为前面一项即 已经使 的取值成为不可能,因此这个冲激可以略去不计。
当然,在奥本海姆的教科书上给出的傅里叶变换积分性质是:
这大抵是利用了 的性质之后,结合 和之后要讲的卷积性质导出的。
8.信号奇偶
类似于傅里叶级数中的讨论,我们还是要求信号为实信号以符合共轭对称性 。这样结合线性性质可得以下结论:
9.时域相乘(频域卷积)
傅里叶变换极其重要的性质之一,就是原信号相乘,那么其频谱函数则是原先两个频谱函数的卷积。设 ,那么:
现在令 ,所以:
从而可以看出 。故:
10.时域卷积(频域相乘)
傅里叶变换的另一重要性质,就是原信号卷积,那么其频谱函数则是原先两个频谱函数的乘积。设 ,那么:
故:
11.频域微分
对傅里叶变换式两边求导可得:
故:
当然,由于 ,上式我们更一般写为:
12.对偶性
对偶性指的是,傅里叶变换式和傅里叶逆变换式在形式上有极大的相似性:
这导致了:如果信号 在傅里叶变换之后得到频谱函数 ,再对频谱函数做变量替换、记 ,那么 在形式上又会与 有极大的相似性!典型的例子是 与单位冲激函数 的傅里叶变换对,可以参考这个例子琢磨一下“对偶性”究竟为何物:
当然您可能也已经注意到了:譬如时域相乘相当于频域卷积,反过来也成立——时域卷积相当于频域相乘;时域微分和频域微分导出的结果竟也差不多……这也是对偶性的体现之一。
六、基本傅里叶变换表
以下给出一些基本信号的傅里叶变换表:
表 1 基本傅里叶变换表
| 信号 | 傅里叶变换 | 条件或说明 |
| — | ||
| — | ||
| — | ||
| — | ||
| — | ||
| — | ||
| — | ||
| — | ||
| — | ||
| — | ||
| — | ||
| — |