傅里叶变换与傅里叶级数_第1页
傅里叶变换与傅里叶级数_第2页
傅里叶变换与傅里叶级数_第3页
傅里叶变换与傅里叶级数_第4页
傅里叶变换与傅里叶级数_第5页
已阅读5页,还剩14页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、重温傅里叶笔记篇 本文记录的大多是基础的公式,还有一些我认为比较重要的有参考价值的说明。(如果对这些公式已经很熟悉,可以直接看第三部分:总结性说明)重温傅里叶笔记篇一、傅里叶级数$ 关于三角函数系的正交性: 三角函数系包括: 1, cos x, sinx , cos2x, sin 2x, cos nx, sinnx, “正交性”是说,三角函数系中的任何一项与另一项的乘积,在 (-, ) 区间内的积分为0。(任何两相的积总可以展成两个频率为整数倍基频的正余弦函数之和或差,而这两个展开后的正余弦在(-, )上积分都为0)。 不同频率(但都是整数倍基频)的两个正弦函数之积,在(-, )上积分恒为0。

2、 同频率的两个正弦函数之积,只有在这两个正弦的相位正交时,其在(-, )上积分才是0。 三角函数系中除“1”以外的任何一项的平方,在(-, )上的积分恒为,“1”在这个区间上的积分为2。 $ 上公式 ! 当周期为2时:式(1):上式成立的条件是f(x)满足狄立克雷充分条件:1.在任意有限区间内连续,或只有有限多个第一类间断点;2.任意的有限区间,都可被分成有限多个单调区间(另一种说法是:任意有限区间内只有有限多个极值点,其实是一样的) 式(1)第一行中的a0/2 就是f(x)的周期平均值,而且第一行的式子只对f(x)是连续函数的情况成立;如果f(x)不连续,则应表示成“(1/2) f(x-0)

3、+f(x+0)”,即f(x)左右极限的算术平均。下面的类似情况都是这样,之后就不再专门说明,这些大家应该都懂。 第三、四行中,n的取值都是:1,2,3,4,n,(都为正,且不包含0)。 当周期为2L时(这也是最一般的情形): 式(2): 第一行中的a0/2 就是f(x)的周期平均值; 第三、四行中,n的取值都是:1,2,3,4,n,(都为正,且不包含0)。$ 傅里叶级数的复数表达方式 同样设周期为2L。根据欧拉公式,正余弦函数都可以用复指数表示出来。这样上面式(2)中的第一行: 可以表示为:令:cn与c-n互为共轭。这样式(4)变为:由式(5)和式(2)中对 a0 b0an bn c0 cn

4、c-n的定义,可以发现cn可统一表达为:将傅里叶级数用复数表示后,就是式(6)和式(7)这样简洁的形式。简单分析:若f(x)为偶(或奇)函数,则所有的bn(或an)将为0,此时的cn将变为实数(或纯虚数),且an(或bn)是转换后所得的cn的2(或2i)倍,而c-n与cn相等(或纯虚共轭)。二、复变函数中的傅里叶变换$ 先上公式: 定理:若f(t)在(-,+)上绝对可积,即 f(t)的绝对值在(-,+)上收敛,则F()在(-,+)上存在且连续(F()的连续性在复变函数的教科书中一般都有证明)。F()是实变复值函数,即变量是在实数区间(-,+)定义,而函数值F()却在复数空间。 式(9)的条件是

5、:f(t)在(-,+)上绝对可积,并在任一有限区间 满足狄立克雷充分条件。$ 若f(t)为偶函数,则F()将为纯实数,且同为偶函数; 若f(t)为奇函数,则F()将为纯虚数,且同为奇函数; 而对任意f(t),F() 与F(-)始终共轭,这意味着 |F()| 与 | F(-)| 恒相等,即F()的绝对值是偶函数。$ 由于要求f(t)绝对可积,所以对于周期函数一般是不能用傅里叶变换的,只能用傅里叶级数分析。(周期函数往往不能收敛)。 三、总结性说明 周期函数可以看成由很多频率是原函数频率整数倍的正余弦波叠加而成,每个频率的波都有各自的振幅和相位,必须将所有频率的振幅和相位同时记录才能准确表达原函数

6、。但从上面的公式来看,我们好像从没涉及到相位?其实不然,从式(2)来看,我们将每个频率的波分成了一个正弦分量和一个余弦分量,同时记录了这两个分量的振幅an、bn其实就已经包含了这个频率的波的相位信息;而对于式(6a),每个频率的波被分成了正负两个频率的复数“波”,这种方式其实比正余弦形式更加直观,因为复振幅cn恰好同时记录了这个频率的振幅和相位,它的物理意义很明显:cn的幅值 |cn| 即为该频率的振幅(准确的说是振幅的一半),而其辐角恰好就是相位(准确的说是反相的相位,c-n的辐角才恰好代表该频率波分量的相位)。 傅里叶变换针对的是非周期函数,或者说,周期为无穷的函数。它是傅里叶级数的一个特

7、例(好吧,我曾经一直以为傅里叶级数是傅里叶变换的一个特例,正好相反,刚前几天才想通透)。当傅里叶级数的周期L趋于无穷时,自然就变成了上面的傅里叶变换。这种关系从二者的表达式中大概能看出点端倪,但是也不是特别明显,毕竟它们的表达形式差别还挺大。如果不把傅里叶级数表达成复数形式,那就更加难看出二者之间的联系了,这也是为什么本文中详细列出了复数形式的傅里叶级数。傅里叶变换要求f(t)在(-,+)上绝对可积,其实可以理解成“傅里叶级数要求函数在一个周期内的积分必须收敛”。在深入篇中,我再好好说说二者是如何联系的。重温傅里叶-深入篇1-傅里叶级数与傅里叶变换的关系以及频谱图的介绍 在读本文前,请先大致浏

8、览一下笔记篇里的东西,下面使用的符号及其意义都跟笔记篇里是一致的。笔记篇里记录的大都是基础的公式,教科书上都可以找到。(抱歉,刚发现有点小错误:在式(6-4)和式(11)里,积分项中的“dx”都应改为“d”,由于改图不太好改,就只在这里说明了。请读者看的时候注意)为了下面叙述方便,我先做几点约定和说明: 本文中提到的傅里叶级数都是复数形式的级数,下标n都是负无穷到正无穷; 对于笔记篇里经常出现的“ n/L ”,它可以看成一个角频率,用表示。(角频率与频率(通常用f表示)之间的关系是:=2f)。(参见笔记篇中的式(3)、(4)、(6)等); 进一步,我将“/L”称为“角基频”, 这样的话“ n/

9、L ”就是n倍角基频。当周期为2时,角基频恰好为1; 一定别搞混:cn代表的不是角频率为n的波分量的振幅,而是角频率为n倍角基频的波分量的振幅; 对于周期函数,除了角频率为整数倍(包括负整数倍)角基频的波分量振幅可以不为0外,角频率为其他值的波分量振幅都是0。(下面介绍频谱图时会再提到此事); *对于周期L等于无穷大的函数(非周期函数),其角基频为/L = 0 ,这样实数范围内的所有角频率都可以看成整数倍角基频了,因此非周期函数在所有的角频率处都有波分量!(就是说,频谱图由离散变得连续了)。什么,那不乱套了?如果所有的角频率都有波分量而且每个波分量都有一个不为0的振幅,那级数怎么可能收敛?还好

10、,每个cn的表达式中都有一个 1/2L 的系数,这样周期无穷大时,所有的振幅cn也都变成“0”了,所以不会乱套,但是这么多0加一块应该还是0,怎么能凑出原来的f(x)呢?这就像对一个函数积分一样,函数在任意一个点处的积分都是0(好吧我知道这说法不科学,但是方便理解),但对一个区间积分,这么多0加起来就成了一个有限值。好了,不乱说了,越说越乱,本文就从这里开始,看完下面的几段大家就能清楚的知道是怎么一回事了。 为了方便大家翻阅,我先将一会儿涉及到的几个公式重新贴一遍在这里。这些公式及公式的标号都与笔记篇中相同。周期级数公式如式(6)和式(7)那样,我们现在要做的是,搞明白为什么周期L趋于无穷时,

11、就会有式(9)和式(8)的结果。好,现在我们对式(6)和式(7)进行第一步加工:将式中的“ n/L ”用角频率n来表示,代表n倍角基频。这样,会产生下面的新式子:对比式(7-1)和式(8),发现他们右边的积分式主体部分形式几乎是一样的,只是上下限和系数不同。好吧,为了更直观的对比,我再创造一个符号,Fn,将它定义如下: Fn= cn 2L 这样我们就可以彻底抛弃cn 这个碍眼的符号了,全部用Fn代替。然后重写式(6)和式(7):再拿式(7-2)和式(8)对比,会发现很让人兴奋的结果,他们的形式几乎一样!但是式(6-2)和式(9)貌似差别还不小,他们的系数一个是(1/2L),一个(1/2)。好吧

12、,接着来,我们再创造一个符号,定义如下: = (/ L) (其实就是角基频的大小)利用它来再次加工式(6):(式(7-2)不变,但还是一块列了出来)重新对比式(6-3)和式(9),发现形式已经很相近了,只不过一个是积分一个是和式等一下!和式?再仔细看看看式(6-3),发现这时它很像一个函数积分的和式展开式!那我们现在来构造两个函数吧:F* ()和* (),构造方法如下: F*() = Fn 当 ( n - 1/2 ) ( n + 1/2 ) 时;* () = n 当 ( n - 1/2 ) ( n + 1/2 ) 时;这是两个分段跳跃函数,它们都以为自变量,并每隔,函数值变化一次。好吧,数字太

13、不直观,我把F*()的函数图象大致画出来方便大家理解: 上面这个阶梯状的东西就是F* ()的函数图象。* ()的图像也是类似的阶梯状,而且它的更简单,是一个从负无穷到正无穷逐步升高的形状(每次升高一个角基频的大小)。 这里有必要说明一下,以免误导大家:Fn 一般都是复数,只有在f(x)本身是偶函数时才是实数,因此函数F*的值也应为复数。也就是说,将F*的函数图象画成图1那样的实数形式其实是不合理的。我这样做只是为了方便大家理解(6-3)中的和式是如何变成积分式的。好了,有了这两个函数,我们再来仔细看看式(6-3),不难看出,这个和式其实就是函数F*在(-,+)上的积分(面积)!这次我们再进一步

14、,将上面两个式子中的Fn和n也都换掉,使其变成*和F*这两个函数之间的关系式(离成功不远了): 这就是转换后的结果。笔记篇中的式(6b)与式(7),跟现在推出的式(6-4)与式(7-4),是完全等价的,因为后面的两个就是根据前两个换算来的,只不过借助了F*()和* ()这两个新构造的函数而已。 表达的意义一样,适用范围也一样(都适用于周期函数),但形式却大变! 这时再回头看看式(9)和式(8),我们终于可以松口气了,形式完全一样!好了,现在我们再看看看周期L趋于无穷时会发生什么。如果直接分析笔记篇中的式(6b)与式(7),我们会很失望,因为L趋于无穷时,它们都“退化”了,很难直接地从这两个式子

15、中得到有用的信息(如果用这两个式子,我们所能得到的“直观”结果就是:cn 全变0了,所以f(x)是0。显然这是错的)。但我们后来创造出来的式(6-4)与式(7-4),适应环境的能力就很强了。1. 首先,L趋于无穷时,会变得越来越小直至变成0(是什么?忘了?前面有, = (/ L);2. 同时,对于* () = n,由于其实就是角基频,而相邻的两个n差就是一个角基频,根据1可知,L趋于无穷时,* ()就由阶梯跳跃变得连续了,这时* () =。3. 同时,两个 相邻的Fn,他们的差别也越来越小直至变成0,(Fn =cn 2L ,从cn的表达式可以看出,L趋于无穷时cn本身就是一个与(1/L)同阶的

16、无穷小量,那相邻的cn之间的差值就是比(1/L)更高阶的无穷小量,因此相邻的Fn之间的差值就趋于0了)。OK完结,多么简单,可是以前就没想到,刚现在才开窍。 数字游戏玩完之后,我们再好好理解一下式(8)(9)中的F()。从我们刚才的证明过程中,可以看到 Fn = cn2L ,在笔记篇中我说过,cn其实就代表某个频率波分量的振幅和相位,而Fn与cn是成正比的,它的值同样可以表征一个波分量的振幅和相位。F()与Fn有相同的意义,因此F()的分布其实就代表了各角频率波分量的分布。具体的说: |F()| 的分布正比地体现了各个角频率波分量的振幅分布。(别忘了F()是复数) F()的辐角体现了各个角频率

17、波分量的相位分布。 我们平时所说的“频谱图”,其实指的就是| F()|的函数图象,它始终是偶函数(这个就是实数了,因为我们取的是F()的幅值而不是F()本身)。对于满足傅里叶变换条件的非周期函数,他们的频谱图一般都是连续的;而对于周期函数,他们的频谱则都是离散的点,只在整数倍角基频的位置有非零的频谱点存在。根据频谱图可以很容易判断该原函数是周期函数还是非周期的(看频谱图是否连续就行了),而且对于周期函数,可以从频谱图读出周期大小(相邻的离散点之间的横轴间距就是角基频,这个角频率对应的周期就是原函数的周期)。那怎样读出每个频率的振幅呢?| F()| 与振幅成正比,要想读出某个频率波分量的实际振幅

18、,只需让 |F()| 乘以相邻离散点的横轴间距再除以即可。其实就是让 |F()| 除以原函数周期值的一半(即L),参考一下我们上面说到的Fn和cn之间的关系式以及我在笔记篇中提到的“|cn|的幅值是实际振幅的一般”,就可以轻松得到得到这个结论。对于非周期函数来说,其频谱图已趋于连续,相邻“离散点”的横轴间距就是一个无穷小量,而 |F()| 是有限值,那么每个频率波分量的实际振幅就都是0了。所以对于非周期函数,说“|F()| 代表了振幅密度的大小”比说“ |F()| 代表了振幅的大小”更贴切一点。在某个宽度为的区间内(频带),对这个“密度”进行积分,(其实还要再除以的)就能得到这个宽度为的频带中所有频率产生的振幅之和(虽然大家的振幅都是趋于0,但无数个加一块就有非零值了)。怎么理解呢?先把这个连续频谱图想象成一个由很多离散点组成的离散频谱图,只不过相邻离散点之间的横轴间距特别小(用d表示吧,方便我叙述),其实相当于先把这个非周期函数想象成了一个周期很长的周期函数(周期越长,相邻离散点的横轴间距/L 越小),然后用周期函数那一套

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论