剑客
关注科技互联网

我理解的傅里叶变换

我理解的傅里叶变换

本科学的电子信息工程,用的最多的就是傅里叶变换。但大多数 EE 的学生只是会套用傅里叶变换的公式而已,对其数学本质却了解得比较少。我不敢说掌握了数学本质,只能说从某个角度有了一些理解。

之前知友@Heinrich 写过一篇傅里叶分析之掐死教程,这是一篇非常浅显易懂的傅里叶变换科普文,写得很好。而本文则侧重从数学的角度(准确地说是从线性空间的角度)介绍傅里叶变换,要求读者具备基本的高等数学和线性代数知识,目的是使读者对傅里叶变换有一个更好的理解。

本文分为五部分。第一部分介绍一些基本的数学概念,第二部分介绍傅里叶基和傅里叶级数,第三部分介绍傅里叶变换,第四部分回归 EE 的内容,第五部分总结。

一、基本概念

1.域

域是代数中一个非常重要的概念,是一种特殊的环(交换除环)。介绍域是为介绍线性空间做准备,而域的准确概念在本文中并不重要。因此为了方便,大家可以把域理解为一个 可以做加减乘除的集合 。也就是说,域由一个集合和四则运算构成,这个集合内的元素两两间可以做这四种运算(除了除以 0),结果依然在这个集合里。

常见的域有有理数域 我理解的傅里叶变换 、实数域 我理解的傅里叶变换 和复数域 我理解的傅里叶变换 ,它们分别是全体有理数、实数、复数关于数的加减乘除构成的。

2.线性空间

线性空间的概念是建立在域的基础上的。假设 我理解的傅里叶变换 是一个域,而 我理解的傅里叶变换 是一个集合。如果 我理解的傅里叶变换 内的元素之间可以做 加法 (也就是说两个元素做加法的结果仍然在集合内), 我理解的傅里叶变换我理解的傅里叶变换 可以做 数乘 (和域内的乘法不同,数乘是说 我理解的傅里叶变换 中的一个元素和 我理解的傅里叶变换 中的一个元素做数乘结果仍在 我理解的傅里叶变换 内。数乘可以理解为一个二元函数,它把 我理解的傅里叶变换 内的一个元素和 我理解的傅里叶变换 内的一个元素映射到 我理解的傅里叶变换 内的一个元素),这个 我理解的傅里叶变换 就叫做 我理解的傅里叶变换一个线性空间 。( 我理解的傅里叶变换 中的元素可以称为 向量

当然了,上边所说的加法和乘法还要满足几个性质,比如加法要构成阿贝尔群等等。但为了方便,这里不详细说明,只是举几个线性空间的例子。

比如对 我理解的傅里叶变换 (n 维欧式空间)是实数域 我理解的傅里叶变换 上的线性空间。因为空间中任意两个向量做加法(按照一般向量加法)或用一个实数乘一个向量,结果仍然是一个 n 维欧式空间里的向量。

而定义在闭区间 我理解的傅里叶变换 上的连续函数也构成一个 我理解的傅里叶变换 上的线性空间。因为任意两个连续函数的和仍然是原来区间上的连续函数,一个连续函数乘一个实数也是连续函数。

下面将介绍线性空间有关的一些概念,先介绍线性组合。

3.线性组合

我们知道线性空间内的元素可以做加法,线性空间中的元素可以和域上的元素做数乘。因此我们取 我理解的傅里叶变换 中的一些向量 我理解的傅里叶变换我理解的傅里叶变换 中的一些数 我理解的傅里叶变换 ,然后做运算

我理解的傅里叶变换

上式的 就叫做 我理解的傅里叶变换 的一个 线性组合

这样,给定一些 我理解的傅里叶变换 中的向量,就可以通过它们的线性组合得到更多向量。而我们比较关心能不能得到 0 向量(0 向量是每个线性空间中唯一满足下列条件的一个元素:对任意元素 我理解的傅里叶变换 )。显然当所有的 我理解的傅里叶变换 都取域中的 0 元素(和 0 向量的定义类似,因为域中的元素和线性空间中的元素均构成阿贝尔群)时结果为 0 。那么当至少有一个 我理解的傅里叶变换 不为 0 时,它们的线性组合能不能为 0 呢?

如果可以为 0 ,我们就称这些向量 线性相关 ,反之称为 线性无关 。我们更关心线性无关的情况,它是我们介绍基和维数必不可少的。

4.基与维数

我们在 我理解的傅里叶变换 中找到了一些线性无关的向量,这些向量可以通过线性组合构成无数其他向量。细心的你可能已经发现了:这些向量的所有线性组合构成的向量就构成了一个线性空间。这个空间中的元素一定属于 我理解的傅里叶变换 ,因此它叫做 我理解的傅里叶变换 的一个子空间。

如果构成的线性空间恰好是 我理解的傅里叶变换 ,我们称这些向量是 我理解的傅里叶变换一组基 ,而向量的个数叫做 我理解的傅里叶变换维数 。如果 我理解的傅里叶变换 中有无限多个线性无关的向量,它就是无限维的。

在这里我们不加证明地给出:任何非零线性空间均有基。

5.内积

之前我们介绍线性空间包括了两种运算:第一种是空间内向量的加法,第二种是空间内向量与域中元素的数乘。而内积则是空间内两向量的运算,或者我们可以理解为把空间中两个元素映射为域中一个元素的一个二元函数。我们规定这个函数应该满足一些性质。

首先我们规定,内积是双线性的。在介绍双线性之前,我们先讲讲线性。线性就是说一个一元函数把线性组合映射为线性组合,即 我理解的傅里叶变换 恒成立.

而双线性是针对二元函数的,它是说固定其中任何一个变量后,这个一元函数都是线性的。

除了双线性,内积还必须是 正定的 。即一个向量和它自己的内积必须是非负的,而非零向量和它自己的内积必须是正的。

容易验证,对于欧式空间中的向量点乘是满足这个性质的。而我们前边提到的闭区间上的连续函数,可以定义其内积为 我理解的傅里叶变换 .

最后,我们定义两个向量 正交 是指它们的内积为零。

6.赋范线性空间

我估计很多人能看到这已经准备选择狗带了,但是我还是要讲。

我保证这是最后一个概念了……而且我尽量简单粗暴。

首先我们介绍一个概念:范数。对于 我理解的傅里叶变换 上的线性空间,我们定义它的范数为

我理解的傅里叶变换

范数除了要满足正定性之外,还要满足三角不等式( 我理解的傅里叶变换 )、CS不等式( 我理解的傅里叶变换 )。

我理解的傅里叶变换 上的线性空间如果构造了这样的范数,就称它为 赋范线性空间。

我不知道有多少人有耐心读到这,如果上边的内容有疑问可以去看任何一本高等代数教材,各种概念并不需要理解得很具体,建议理解得不透彻的读者 利用欧式空间进行类比 。下边将开始介绍傅里叶基和傅里叶级数。

二、傅里叶基和傅里叶级数

1.傅里叶基

我们回忆一下刚才讲的正交的概念:两个向量的内积为 0 。如果一个空间有一组基两两正交,那么它就叫做一组 正交基

我们可以证明,定义在 我理解的傅里叶变换 上的所有平方可积的函数构成线性空间,规定其内积为

我理解的傅里叶变换

(为了简单起见,只考虑实值函数,复函数则其中一个取共轭即可)

我们不加证明地给出:这个空间内的每个向量都可以表示为其一组基的 无限线性组合 。这就是 Fourier 展开 。常取的基是

我理解的傅里叶变换我理解的傅里叶变换

容易证明这些基都是正交的,而且不能找到另外一个向量和它们全正交(这样的一组向量叫做 极大正交向量族 )。这样,我们就可以把任何一个 我理解的傅里叶变换 上的平方可积的函数通过这些基的线性组合表示出来,这组基就叫做 傅里叶基

2.傅里叶级数

傅里叶级数实际上就是把这个空间中的一个向量通过基的线性组合的方式写出来。我们已经知道了,两个正交的向量内积为 0 ,因此如果一个函数是由正交基的线性组合表示,我们可以很容易地求得它的系数。

我理解的傅里叶变换 (这就是傅里叶级数,实际上就是一个基的线性组合, 是 1 的系数, 和 我理解的傅里叶变换 分别是 和 我理解的傅里叶变换 的系数)。现在我们想计算每一个基的系数。为此,我们试着计算一下它和其中一项的内积,比如

我理解的傅里叶变换

不过注意到不同的向量都是正交的,所以非零项只有 我理解的傅里叶变换 .

好简单啊!也就是说, 我理解的傅里叶变换 .

这么说似乎还是不明显,我这么写:

我理解的傅里叶变换 .

类比一下欧几里得空间。比如说二维欧式空间 我理解的傅里叶变换 中,我取一组正交基 我理解的傅里叶变换 ,那么向量 我理解的傅里叶变换 写成这组基的线性组合时系数分别是多少呢?答案就是

我理解的傅里叶变换我理解的傅里叶变换 .

完全一样耶!

实际上从空间的角度,二者确实没有什么本质区别。

那么现在我们已经理解了,傅里叶基其实就是这个函数构成的空间中的一组正交基,而傅里叶级数就是把空间里的元素写成基的线性组合。但我们注意到这里的函数是 我理解的傅里叶变换 的,如果换一个区间,结果会如何呢?

三、傅里叶变换

我们把区间换成 我理解的傅里叶变换 ,这样相当于把每个函数“拉伸”了 我理解的傅里叶变换 ,这样每个基也拉伸了 我理解的傅里叶变换 倍,变成了

我理解的傅里叶变换我理解的傅里叶变换 , 我理解的傅里叶变换 .

而傅里叶级数就变成了

我理解的傅里叶变换我理解的傅里叶变换 .

下面我们考虑指数函数形式的傅里叶级数在 时的情形。

1.在区间变为 我理解的傅里叶变换 后,傅里叶系数变为 我理解的傅里叶变换

2.在 我理解的傅里叶变换 的傅里叶展开中,把傅里叶系数带入,得到

我理解的傅里叶变换

3.上式的积分式是关于变量 的,而外边的项 我理解的傅里叶变换 与 无关。因此可以写成

我理解的傅里叶变换

4.我们现在希望上式是黎曼和的形式,这样就可以写成定积分。为此取 我理解的傅里叶变换 ,这样上式变成

我理解的傅里叶变换

5.当 时, 我理解的傅里叶变换 ,上式作为黎曼和的形式可以写成定积分:

我理解的傅里叶变换

6.上式略作变形,得到

我理解的傅里叶变换

这样,我们从傅里叶级数出发,在 的情形下得到了一个式子。这个式子的本质仍然是一个求和式,只是因为极限所以写成了积分式。它把 我理解的傅里叶变换 写成了 我理解的傅里叶变换 的线性组合。

等等!

线性组合?!

确实如此啊。只是这里的变量 是取遍整个数轴的。相当于对实数轴上的每个点都对应了这个函数的一个“基”,而这个积分式就是这些基的“线性组合”。

括号里的那堆东西,也就是

我理解的傅里叶变换

就叫做这个函数的 傅里叶变换

傅里叶变换是这样一个函数,它在 处的函数值 我理解的傅里叶变换 表示函数 我理解的傅里叶变换 在 对应的基上的系数。至此我们就完成了傅里叶变换从空间角度的介绍。

补充:

前边我们已经定义了范数,而两个向量之间的距离可以用它们差的范数刻画。空间中有一向量序列 我理解的傅里叶变换柯西列 ,如果它满足对任意的 我理解的傅里叶变换 ,存在正整数 我理解的傅里叶变换 ,使得对任意的 我理解的傅里叶变换 ,有 我理解的傅里叶变换 。也就是说,在柯西列中,我们总可以去掉有限个元素,使得剩余的项两两之间的最大距离小于一个给定的正数。

柯西列都是收敛的。如果空间中任何一个柯西列都收敛到空间内部,我们就称这个空间是 完备的

完备的内积空间叫做 希尔伯特空间 。希尔伯特空间中的极大正交向量组称为 希尔伯特基 。傅里叶基是 Lebesgue 函数空间 我理解的傅里叶变换我理解的傅里叶变换 指平方可积函数构成的函数空间,它和平方可和序列构成的空间均为希尔伯特空间)上常取的一组基。

另外,从基的角度,先引出傅里叶级数,再通过极限去说明傅里叶变换更好一些。但其实傅里叶变换更加自然,而傅里叶级数只是紧致阿贝尔群上的傅里叶变换。

四、回归信号

我们再从信号分析的角度考虑上述结果。我们求傅里叶变换说的是把信号从“ 时间域 ”变换到“ 频率域 ”,而频率域对应值就代表“信号”在该点的“ 频率分量 ”。为什么我们可以说“频率分量”呢?就是因为每一个频率都代表了一个基,原信号可以写成这些基的线性组合,而每个基上的系数就代表了信号在这个基(频率)上的分量大小。

回忆我们高中时候学物理,做平抛、斜抛运动的题时,总是二话不说地水平竖直分解运动,为什么呢?因为对于平面内的运动,水平和竖直就是两组“基”(并不十分准确,理解即可)。那么为什么非得水平和竖直而不是斜 30° 角呢?这我们就要考虑重力了。重力只在竖直方向有作用,它在水平这组基下的分量为 0 。所以说,基的选取也很重要。

而信号的分解为什么要在三角函数或指数函数里进行分解呢?首先我们要说明,这两个分解的本质是一样的,因为三角函数的定义就是通过指数函数的(与欧拉公式无关,见 余翔关于欧拉公式的一个回答 )。一方面,三角函数具有非常好的性质,简单的三角函数非常容易分析,产生也很自然。另一方面,三角函数非常好地刻画了信号的频率,而这也是我们非常关心的。

从基的线性组合的角度,傅里叶变换的很多性质都是显然的,请读者自己试着证明。

五、总结

傅里叶变换是 EE 专业学生的“命根子”,几乎所有的信号处理都需要依靠傅里叶变换。毫无疑问,对傅里叶变换更深层的理解,有助于我们更好地学习和研究。本文的理解角度也是我认为最适合 EE 专业学生的理解方式。

然而本文只是一个概念性的介绍,如果要真正地理解,可以仔细阅读参考文献中列举的资料。如果文章有不妥之处,也请大家批评指正。(因为本文不是介绍知识,里边掺杂了自己的理解)

参考文献

[1] 张贤科. 高等代数学, (第二版)[M]. 清华大学出版社, 2004.

[2] Albert, Boggess, 赵树森. 小波与傅里叶分析基础[M]. 电子工业出版社, 2010.

分享到:更多 ()

评论 抢沙发

  • 昵称 (必填)
  • 邮箱 (必填)
  • 网址