最新傅里叶级数和傅里叶变换.docx
《最新傅里叶级数和傅里叶变换.docx》由会员分享,可在线阅读,更多相关《最新傅里叶级数和傅里叶变换.docx(95页珍藏版)》请在冰豆网上搜索。
最新傅里叶级数和傅里叶变换
傅里叶级数和傅里叶变换
第九章傅里叶级数和傅里叶变换
在自然界中广泛地存在各种各样的周期性运动(即相隔一定时间间隔往复循返的过程)。
例如,日月星球的运动,海洋潮汐的运动,电磁波与声波的运动,工厂里机器部件的往复运动,时钟摆的摆动以及人体心脏的跳动等等,都是周期性运动。
为了描述周期性的运动过程,数学上是借助某类函数来描述的。
当然这类函数也要体现出周期性。
这类函数称为周期函数。
在前面几章中,为了研究函数的性质,常常采用分析表示法,将这些函数在某区域展开成幂级数的形式,如泰勒级数或罗朗级数。
但是,这种幂级数形式的展开式是体现不出周期性来的,那么,对于周期性函数应采取怎样的分析表示法呢?
这就是本章要讨论的内容。
9.1周期函数和傅里叶级数
9.1.1周期函数
凡满足以下关系式:
«SkipRecordIf...»(T为常数)(9.1.1)
的函数,都称为周期函数。
周期的定义
(1)满足式(9.1.1)的T值中的最小正数,即为该函数的周期;
(2)一个常数以任何正数为周期。
9.1.2基本三角函数系
按某一规律确定的函数序列称为函数系。
如下形式的函数系:
1,«SkipRecordIf...»,«SkipRecordIf...»,«SkipRecordIf...»,«SkipRecordIf...»,…,«SkipRecordIf...»,«SkipRecordIf...»,…(9.1.2)
称为基本三角函数系。
所有这些函数具有各自的周期,例如«SkipRecordIf...»和«SkipRecordIf...»的周期为«SkipRecordIf...»,但它们的共有周期为«SkipRecordIf...»(即所有周期的最小公倍数)。
通常这个周期命名为函数系的周期。
所以式(9.1.2)的三角函数系的周期为«SkipRecordIf...»。
如果我们将基本三角函数系中的函数,任意取n个组合,则我们可以得到一个较复杂的函数。
例如图9.1(a)是两个函数的组合«SkipRecordIf...»;图9.1(b)是三个函数的组合«SkipRecordIf...»。
如果我们取跟多的函数组合,甚至全体的组合,将会得到更复杂的函数或我们期望的函数。
9.1.3傅里叶级数
现在我们讨论上述问题的逆问题。
即如果给定一个周期为«SkipRecordIf...»的任意周期函数«SkipRecordIf...»
«SkipRecordIf...»(9.1.3)
我们能否将它表示成简单的三角函数(有限个或无限个)之和呢?
即能否将«SkipRecordIf...»分解成如下形式:
«SkipRecordIf...»(9.1.4)
如果能实现这种分解,那么对许多复杂的函数就可以通过简单的三角函数来研究其性质了。
上述问题的回答是肯定的,并称式(9.1.4)为函数«SkipRecordIf...»的傅里叶级数(有时我们也称它为狭义傅里叶展开式)。
若函数«SkipRecordIf...»按非三角函数系«SkipRecordIf...»…)进行展开,所得的级数称为广义傅里叶展开式。
因为这一问题,最早由工程师J.Fourier提出来的。
他在1807年12月21日向权威的法国科学院宣告:
任意的周期函数«SkipRecordIf...»都能展开成正弦及余弦的无穷级数,即形式(9.1.4)的级数。
他的宣告震怒了整个科学院。
当时许多杰出的院士,包括著名的法国数学家拉格朗日等,都认为他的结果是荒谬的。
因为那时它在数学上没有得到严格的证明。
然而,现在数学家已经把“傅里叶级数理论”发展到相当高的水平,已有许多专著,并建立这级数收敛的非常明确的条件。
这结果被认为是20世纪所发现的最重要的数学理论之一。
以上这种分解不仅具有严格的数学基础,而且还具有真实的物理背景。
我们也可以用实验来证明这种分解过程。
例如将矩形脉冲或半波整流电路的输出波形输入到一个选频放大器中,再将选频放大器的输出端接到示波器上,调整选频放大器的频率,就可以看到示波器上出现各种不同频率的正弦波,这说明矩形脉冲或半波整流电路的输出波形可以看作许多不同频率的正弦波的叠加。
特别是近几年来配备有数字电子计算机的专用仪器相应问世(如频率分析仪、快速傅里叶变换处理机、信号处理机等),想这种分解过程在很短的时间内就可以完成。
将一个周期为«SkipRecordIf...»的任意周期函数按基本三角函数系展开,首先需要解决如下两个问题:
(1)在什么条件下«SkipRecordIf...»才能按基本三角函数系展开?
(2)如何确定展开式中系数?
也就是说,如果«SkipRecordIf...»可以展开成下式:
«SkipRecordIf...»
(其中,«SkipRecordIf...»前的«SkipRecordIf...»系数,是为了使今后系数公式更为对称而引人的)那么«SkipRecordIf...»…及«SkipRecordIf...»…这些系数应如何确定?
为了解决这个问题,还需要介绍如下几个概念。
§9.2完备正交函数系
由于我们前两章的知识,我们可以在任意区间«SkipRecordIf...»上建立一组完备正交函数系,首先从如下标准的S-L本征值问题(其中«SkipRecordIf...»)出发,
«SkipRecordIf...»«SkipRecordIf...»(9.2.1)
以满足方程的解«SkipRecordIf...»代人边界条件
«SkipRecordIf...»
即
«SkipRecordIf...»«SkipRecordIf...»,«SkipRecordIf...»,«SkipRecordIf...»,…)
«SkipRecordIf...»(9.2.2)
由于周期性边界条件有两个线性独立的解(复数解的实部和虚部)。
因此,我们这样来排序,令
«SkipRecordIf...»
«SkipRecordIf...»
«SkipRecordIf...»,«SkipRecordIf...»,«SkipRecordIf...»
(9.2.3)
根据定理(8.2.1),上述{«SkipRecordIf...»}在区间«SkipRecordIf...»上构成完备正交系,并且任何一个在«SkipRecordIf...»上具有连续的一阶导数和分段连续的二阶导数,又适合本征值问题的边界条件的函数«SkipRecordIf...»,都可按{«SkipRecordIf...»}展开成在«SkipRecordIf...»上绝对且一致收敛的级数,即
«SkipRecordIf...»
其中
«SkipRecordIf...»«SkipRecordIf...»(9.2.4)
由于
«SkipRecordIf...»
«SkipRecordIf...»(«SkipRecordIf...»=1,2,3,…)
也可表示成
«SkipRecordIf...»(9.2.5)
其中
«SkipRecordIf...»
«SkipRecordIf...»(9.2.6)
«SkipRecordIf...»
如果令«SkipRecordIf...»=0,«SkipRecordIf...»=2«SkipRecordIf...»,即在区间«SkipRecordIf...»上下列函数系«SkipRecordIf...»:
«SkipRecordIf...»(«SkipRecordIf...»=1,2,3…)
(9.2.7)
构成完备正交函数系。
如果令«SkipRecordIf...»=-«SkipRecordIf...»,«SkipRecordIf...»=«SkipRecordIf...»,即在区间«SkipRecordIf...»上下列函数系{«SkipRecordIf...»}:
«SkipRecordIf...»(«SkipRecordIf...»=1,2,3…)(9.2.8)
构成完备正交函数系。
这函数系就是通常的基本三角函数系(9.1.2)。
展开式(9.2.5)中的本征函数系排序是按区间«SkipRecordIf...»上本征函数的零点个数来排序的,但为了照顾习惯和历史,我们仍按通俗的方式来排序。
于是,«SkipRecordIf...»按基本三角函数系展开式为
«SkipRecordIf...»(9.2.9)
其中系数
«SkipRecordIf...»(«SkipRecordIf...»=0,1,2,…)(9.2.10a)
«SkipRecordIf...»(«SkipRecordIf...»=1,2,3,…)(9.2.10b)
利用帕塞瓦尔等式(7.3.8),有
(ⅰ)若«SkipRecordIf...»,则«SkipRecordIf...»(9.2.11)
(ⅱ)若«SkipRecordIf...»;«SkipRecordIf...»,则
«SkipRecordIf...»(9.2.12)
应用到«SkipRecordIf...»按基本三角函数系展开式(9.2.9)上,并注意«SkipRecordIf...»
(ⅰ)«SkipRecordIf...»(9.2.13)
(ⅱ)«SkipRecordIf...»(9.1.14)
综上所述,对于任意函数«SkipRecordIf...»可按函数系{«SkipRecordIf...»}进行展开的条件是
(1)函数«SkipRecordIf...»在«SkipRecordIf...»上是绝对可积的连续函数(今后这条件可放宽,由狄利克雷条件代替);
(2)函数系必须在«SkipRecordIf...»上是完备正交系。
以上两个条件才能保证«SkipRecordIf...»展开如下形式:
«SkipRecordIf...»(9.2.15)
其中系数由下式确定:
«SkipRecordIf...»(«SkipRecordIf...»=0,1,2,…)(9.2.16)
对于基本三角函数系(9.1.2),正是因为她在区间«SkipRecordIf...»上构成正交完备系,才能让绝对可积的连续函数«SkipRecordIf...»展开成(9.2.9)形式的傅里叶级数,其中系数由公式(9.2.10)确定。
如果用«SkipRecordIf...»来表示函数«SkipRecordIf...»与其展开式的前«SkipRecordIf...»+1项之部分和的均方偏差,即
«SkipRecordIf...»(9.2.17)
要使«SkipRecordIf...»取极小值,对所有«SkipRecordIf...»下式必须成立:
«SkipRecordIf...»(«SkipRecordIf...»=0,1,2,…,n)
«SkipRecordIf...»(«SkipRecordIf...»=1,2,…,n)
这就导致求系数«SkipRecordIf...»及«SkipRecordIf...»的公式,即傅里叶系数公式(9.2.10a)及(9.2.10b),并且这些公式与我们取多少项的部分和求均方偏差无关。
在复杂波形的分析中(海洋潮汐、地震、乐调等),也许更方便的是将傅里叶级数写成如下形式:
«SkipRecordIf...»(9.2.20)
其中系数«SkipRecordIf...»及«SkipRecordIf...»与上述傅里叶级数系数«SkipRecordIf...»及«SkipRecordIf...»的关系为
«SkipRecordIf...»
«SkipRecordIf...»(«SkipRecordIf...»=1,2,3…)(9.2.21)
系数«SkipRecordIf...»又称功率谱,它明确地与«SkipRecordIf...»有关,并且在相角
改变之后,并不变化。
§9.3傅里叶级数的性质
9.3.1收敛性
定理9.3.1傅里叶级数的收敛准则——狄利克雷(Dirichlet)定理
若
(1)«SkipRecordIf...»在«SkipRecordIf...»«SkipRecordIf...»上或者连续,或者只有有限个间断点,在间断处函数的左、右极限都存在;
(2)«SkipRecordIf...»在«SkipRecordIf...»上只有有限个极大值点与极小值点;
(3)«SkipRecordIf...»在«SkipRecordIf...»外是周期函数,其周期为2«SkipRecordIf...»,则级数
«SkipRecordIf...»(9.3.1)
证明
«SkipRecordIf...»«SkipRecordIf...»
=«SkipRecordIf...»
=«SkipRecordIf...»
=«SkipRecordIf...»
因为
«SkipRecordIf...»
及
«SkipRecordIf...»
所以
«SkipRecordIf...»证毕
例9.3.1试将锯齿波«SkipRecordIf...»在区间«SkipRecordIf...»上«SkipRecordIf...»展开为傅里叶级数。
解如图9.2所示,我们要将«SkipRecordIf...»在«SkipRecordIf...»之外视作是2«SkipRecordIf...»的周期函数,按傅里叶级数公式(9.2.10a)及(9.2.10b)有
«SkipRecordIf...»(«SkipRecordIf...»=0,1,2,…)
及
«SkipRecordIf...»
=«SkipRecordIf...»(«SkipRecordIf...»=1,2,3,…)
因此,所求级数为
«SkipRecordIf...»(9.3.2)
由于«SkipRecordIf...»=0是«SkipRecordIf...»的连续点,所以上式两边可划等号。
事实上,也正是如此,可代入数字验证。
而«SkipRecordIf...»=«SkipRecordIf...»是«SkipRecordIf...»间断点,由图9.2可知
«SkipRecordIf...»
按收敛准则,«SkipRecordIf...»傅里叶级数在间断点处应收敛到
«SkipRecordIf...»
事实上,以«SkipRecordIf...»=«SkipRecordIf...»代入级数(9.3.2),得级数和为零。
必须注意,狄利克雷定理中加在«SkipRecordIf...»上的条件
(1)和
(2)是充分的,但不是必要的。
在实际中这些条件通常是满足的,目前还不知道傅里叶级数收敛的必要且充分的条件是什么。
值得注意的是,单从«SkipRecordIf...»的连续性考虑还不能保证傅里叶级数收敛。
9.3.2积分
定理9.3.2如果«SkipRecordIf...»在区间«SkipRecordIf...»上分段连续,其傅里叶级数为
«SkipRecordIf...»
则
F«SkipRecordIf...»(9.3.3)
证明
«SkipRecordIf...»(9.3.4)
利用公式(9.3.2)以及公式(9.2.14),得
«SkipRecordIf...»(9.3.5)
上式代入式(9.3.4),即得所证。
如果原级数中«SkipRecordIf...»,只要用«SkipRecordIf...»代替公式(9.3.4)中的«SkipRecordIf...»即可。
9.3.3微分
定理9.3.3若«SkipRecordIf...»在«SkipRecordIf...»上连续,又«SkipRecordIf...»绝对可积,则有
«SkipRecordIf...»
«SkipRecordIf...»(9.3.6)
其中«SkipRecordIf...»。
利用求系数公式(9.2.10)及分部积分,可以证明
«SkipRecordIf...»(«SkipRecordIf...»=0,1,2,…)
«SkipRecordIf...»(«SkipRecordIf...»=1,2,3,…)
如果«SkipRecordIf...»,则«SkipRecordIf...»的傅里叶级数可通过对«SkipRecordIf...»的傅里叶级数进行逐项求导而得,即
«SkipRecordIf...»(9.3.7)
微分与积分大不相同,例如考虑下列函数(锯齿波):
«SkipRecordIf...»«SkipRecordIf...»
的傅里叶级数为
«SkipRecordIf...»(9.3.7)
对上式逐项微分得
«SkipRecordIf...»
于是得到不收敛的级数
其次,再考虑三角波
«SkipRecordIf...»«SkipRecordIf...»
它的傅里叶级数
«SkipRecordIf...»
是一个收敛得相当快的级数,且在«SkipRecordIf...»上一致收敛。
对上式逐项微分得
«SkipRecordIf...»
上式正是方波
«SkipRecordIf...»
的傅里叶级数。
事实上,三角波得导数正数方波。
从上面的例子可知,与积分相反,微分之后每一个系数前却添加了一个增长因子«SkipRecordIf...»,这就降低了收敛程度。
所以上面第一个例子微分后得一发散级数。
事实上,第一个例子中的级数在«SkipRecordIf...»区间上一致收敛。
一般来说,微分使级数的收敛程度降低。
有时将可以逐项微分的条件表示成如下形式:
(9.3.8)
此外,函数的光滑程度可以从该函数的傅里叶级数的系数上反映出来。
一般而言,一个满足狄利克雷条件的周期函数。
其傅里叶级数中的系数«SkipRecordIf...»和«SkipRecordIf...»随着«SkipRecordIf...»趋向于无穷大时,他们至少应与«SkipRecordIf...»(其中«SkipRecordIf...»为与«SkipRecordIf...»无关的常数)一样快的趋向于零。
如果函数包含一个或几个间断点,那么不是«SkipRecordIf...»就是«SkipRecordIf...»,一般情况是二者都不能比«SkipRecordIf...»更快的趋向于零。
如果函数以及它的前(«SkipRecordIf...»-1)阶导数满足狄利克雷条件,而且处处连续,那么随着«SkipRecordIf...»趋向于无穷大,«SkipRecordIf...»的傅里叶级数的系数«SkipRecordIf...»和«SkipRecordIf...»至少应与«SkipRecordIf...»一样快趋向于零。
如果«SkipRecordIf...»的«SkipRecordIf...»阶导数不处处连续,那么不是«SkipRecordIf...»就是«SkipRecordIf...»,一般情况是二者都不能比«SkipRecordIf...»更快地趋向于零。
因此,函数愈光滑,其傅里叶级数的系数收敛得越快,反之,只要考虑某函数的傅里叶级数的系数的收敛快慢程度,就可以判断该函数的光滑程度。
9.3.4吉布斯现象
在间断点处,傅里叶级数呈现奇特的现象。
作为例子,我们考虑方波(图9.3)
«SkipRecordIf...»
其傅里叶级数为
«SkipRecordIf...»(9.3.9)
除去间断点«SkipRecordIf...»(«SkipRecordIf...»=0,«SkipRecordIf...»,…)之外,上式右边级数之和与«SkipRecordIf...»一致,设级数前N项
的部分和«SkipRecordIf...»,它与«SkipRecordIf...»之差命名为误差项«SkipRecordIf...»,即
«SkipRecordIf...»=«SkipRecordIf...»-«SkipRecordIf...»
因为
«SkipRecordIf...»
«SkipRecordIf...»=«SkipRecordIf...»
又
«SkipRecordIf...»=«SkipRecordIf...»
其中也利用了公式
«SkipRecordIf...»
因此,可将«SkipRecordIf...»表达为积分形式
«SkipRecordIf...»=«SkipRecordIf...»
在间断点«SkipRecordIf...»,
«SkipRecordIf...»
由分析可知,误差项在间断点处取极大值,特别是N→∞,«SkipRecordIf...»为间断点),而«SkipRecordIf...»趋于非零值,即误差总是存在的。
在间断点«SkipRecordIf...»的傅里叶级数之和越过这函数«SkipRecordIf...»的值,超出幅度约为18%,这种现象称为吉布斯现象。
这是由于J.W.Gibbs(1839~1903)第一个向自然杂志投稿,发表了他的研究成果,解释了这一现象。
因此后来就将部分和在间断点附近的异常行为称为吉布斯现象。
9.3.5偶函数和奇函数
如果周期函数为奇函数或偶函数,则傅里叶级数具有如下鲜明的特征:
(1)奇函数
奇函数«SkipRecordIf...»的特征是
«SkipRecordIf...»=«SkipRecordIf...»(9.3.10)
因此,由系数公式(9.2.10),得
«SkipRecordIf...»(«SkipRecordIf...»=0,1,2,…)
及
«SkipRecordIf...»(«SkipRecordIf...»=1,2,3,…)
所以,奇函数«SkipRecordIf...»的傅里叶级数为
«SkipRecordIf...»=«SkipRecordIf...»(9.3.11)
由正弦函数性质可知,在«SkipRecordIf...»=0与«SkipRecordIf...»=«SkipRecordIf...»两端,级数收敛到零,于是
«SkipRecordIf...»=«SkipRecordIf...»=0
综上所述,奇函数的