第八章-傅里叶变换课件

1、1第八章 傅里叶变换28.1 傅里叶积分8.1.1 傅里叶积分的概念在工程计算中, 无论是电学还是力学, 经常要和随时间而变的周期函数fT(t)打交道. 例如:t3具有性质fT(t+T)=fT(t), 其中T称作周期, 而1/T代表单位时间振动的次数, 单位时间通常取秒, 即每秒重复多少次, 单位是赫兹(Herz, 或Hz).t最常用的一种周期函数是三角函数fT(t)=Asin(wt+j) 其中w=2p/T4而Asin(wt+j)又可以看作是两个周期函数人们发现, 所有的工程中使用的周期函数都可以用一系列的三角函数的线性组合来逼近.方波sinwt和coswt的线性组合Asin(wt+j)=as

2、inwt+bcoswt54个正弦波的逼近100个正弦波的逼近6研究周期函数实际上只须研究其中的一个周期内的情况即可, 通常研究在闭区间-T/2,T/2内函数变化的情况. 并非理论上的所有周期函数都可以用傅里叶级数逼近, 而是要满足狄利克雷(Dirichlet)条件, 即在区间-T/2,T/2上这两个条件实际上就是要保证函数是可积函数.1, 连续或只有有限个第一类间断点;2, 只有有限个极值点;7第二类间断点第一类间断点第一类间断点和第二类间断点的区别:8不满足狄氏条件的例: 而在工程上所应用的函数, 尤其是物理量的变化函数, 全部满足狄氏条件. 实际上不连续函数严格上讲都是不存在的, 但经常用

3、不连续函数来近似一些函数, 使得思维简单一些.9 在区间-T/2,T/2上满足狄氏条件的函数的全体也构成一个集合, 这个集合在通常的函数加法和数乘运算上也构成一个线性空间V, 此空间的向量就是函数, 线性空间的一切理论在此空间上仍然成立. 更进一步地也可以在此线性空间V上定义内积运算, 这样就可以建立元素(即函数)的长度(范数), 及函数间角度, 及正交的概念. 两个函数f和g的内积定义为:10一个函数f(t)的长度为11而在区间-T/2,T/2上的三角函数系1, coswt, sinwt, cos 2wt, sin 2wt, ., cos nwt, sin nwt, .是两两正交的, 其中w

4、=2p/T, 这是因为cos nwt和sin nwt都可以看作是复指数函数ejnwt的线性组合. 当nm时,12这是因为13由此不难验证14而1, coswt, sinwt, ., cos nwt, sin nwt, .的函数的长度计算如下:15因此, 任何满足狄氏条件的周期函数fT(t), 可表示为三角级数的形式如下: 16为求an, 计算fT(t), cosnwt, 即17同理, 为求bn, 计算fT(t), sin nwt, 即18最后可得:19而利用三角函数的指数形式可将级数表示为:20如令wn=nw (n=0,1,2,.)21给定fT(t), cn的计算如下:2223如图所示:1-1

5、otf(t)1例 定义方波函数为24现以f(t)为基础构造一周期为T的周期函数fT(t), 令T=4, 则1-12T=4f4(t)t-225则26sinc函数介绍27sinc(x)xsinc函数的图形:28w前面计算出291-17T=8f8(t)t现在将周期扩大一倍, 令T=8, 以f(t)为基础构造一周期为8的周期函数f8(t)30则31w则在T=8时,32w如果再将周期增加一倍, 令T=16, 可计算出33一般地, 对于周期T34当周期T越来越大时, 各个频率的正弦波的频率间隔越来越小, 而它们的强度在各个频率的轮廓则总是sinc函数的形状, 因此, 如果将方波函数f(t)看作是周期无穷大

6、的周期函数, 则它也可以看作是由无穷多个无穷小的正弦波构成, 将那个频率上的轮廓即sinc函数的形状看作是f(t)的各个频率成份上的分布, 称作f(t)的傅里叶变换.35对任何一个非周期函数f(t)都可以看成是由某个周期函数fT(t)当T时转化而来的. 作周期为T的函数fT(t), 使其在-T/2,T/2之内等于f(t), 在-T/2,T/2之外按周期T延拓到整个数轴上, 则T越大, fT(t)与f(t)相等的范围也越大, 这就说明当T时, 周期函数fT(t)便可转化为f(t), 即有36Otf(t)OtfT1(t)OtfT2(t)3738O w1 w2 w3 wn-1wnw如图3940此公式

7、称为函数f(t)的傅里叶积分公式, 简称傅氏积分公式,41定义8.1 称广义积分为傅里叶积分.其中积分变量t取实值且从-到+,为实值参数.可看成复变函数的朗级数加以推广得到.428.2.2 傅里叶积分的物理意义-频谱周期为T的函数fT(t)的傅里叶级数展开式：其中傅里叶系数1. 非正弦周期的频谱序列43 从物理的观点来看，式（1.2)说明fT(t)可表示为频率2n/T的诸振动的叠加.频率为2n/T的第n次振动的振幅2rn和相位n.44 所有出现的诸振动的振幅和相位的全体在物理上由fT(t)所描写的自然现象的频谱. 由于cn的下标n取离散值，所反映的诸振动振幅随频率变化的图形是不连续的状态，故称

8、为离散谱.例8.3 求周期性矩形脉冲函数的频谱序列.45解：经计算可得 故所反映频率为2n/T的第n次振动的振幅2rn和相位n分别为462. 非周期函数的频谱函数对定义在区间（- ,+ )上的非周期函数f(t),可看成周期为+ ，将(1.3)代入（1.2),令T + ,得47若记我们可以讲,由f(t)通过公式（1.4)得到频谱函数F(),反之借助频谱函数又将f(t)作为角频率为的诸振动F( )ei t d /2的迭加形式（1.5)给出.频谱函数F()的模|F()|通常称作f(t)的振幅频谱，这个频谱的图形是连续的，称为连续谱.4849tf(t)|F()|503.傅里叶积分的物理意义频谱函数F(

9、)恰好反映前述频谱序列的和.518.1 若f(t)在(-, +)上满足条件: 1, f(t)在任一有限区间上满足狄氏条件; 2, f(t)在无限区间(-, +)上绝对可积, 则有8.1.3 傅里叶积分定理52例8.5 求矩形脉冲函数的傅里叶积分，傅里叶积公式.解：此函数显然满足傅里叶积分定理条件.故傅里叶积分为53故傅里叶积分公式为由傅里叶积分定理可得到548.2傅里叶变换8.2.1 傅氏变换的定义我们知道, 若函数f(t)满足傅氏积分定理的条件, 则在f(t)的连续点处, 有55(8.6)式叫做f(t)的傅氏变换式, (8.8)式为F(w)的傅式逆变换式.f(t)与F(w)可相互转换,可记为

10、F(w)=Lf(t) 和 f(t)=L-1F(w)56还可以将f(t)放在左端, F(w)放在右端, 中间用双向箭头连接: f(t) F(w) (8.6)式右端的积分运算, 叫做f(t)的傅氏变换, 同样, (8.8)式右端的积分运算, 叫做F(w)的傅氏逆变换. F(w)称作f(t)的象函数, f(t)称作F(w)的象原函数.可以说象函数F(w)和象原函数f(t)构成了一个傅氏变换对.57tf(t)58这就是指数衰减函数的傅氏变换.根据(8.6)式, 有59根据(8.8)式, 有6061可见钟形函数的傅氏变换也是钟形函数因此有如果令b=1/2, 就有62求钟形脉冲函数的积分表达式, 根据(1

11、.8)式63例3 解积分方程解：给函数f(x)在区间(- ,0)上补充定义，使f(x)在区间(- , +)上成为偶函数，则6465这一讲介绍傅氏变换的几个重要性质, 为了叙述方便起见, 假定在这些性质中, 凡是需要求傅氏变换的函数都满足傅氏积分定理中的条件, 在证明这些性质时, 不再重述这些条件.8.2.2 傅里叶变换的性质66这个性质的作用是很显然的, 它表明了函数线性组合的傅氏变换等于各函数傅氏变换的线性组合. 它的证明只需根据定义就可推出.设F1(w)=F f1(t), F2(w)=F f2(t), a,b是常数, 则F af1(t)+bf2(t)=aF1(w)+bF2(w) (1.12

12、)1.线性性质67同样, 傅氏逆变换亦具有类似的线性性质, L -1aF1(w)+bF2(w)=af1(t)+bf2(t) （1.13)例1求函数的傅氏逆变换.解：因为6869证 由傅氏变换的定义, 可知2. 位移性质70例2 求函数的傅氏逆变换.解：因为71例3 证明证：因为723.微分性质如果f(t)在(-, +)上连续或只有有限个可去间断点, 且当|t|+时, f(t)0, 则F f (t)=iwF f(t).(1.16)证 由傅氏变换的定义, 并利用分部积分可得73推论 L-1 f(n)(t)=(iw)n L-1 f(t).(1.17)同样, 我们还能得到象函数的导数公式, 设F f(

13、t)=F(w), 则74例4 求函数的傅氏变换解：设则75764. 积分性质77的解, 其中t+, a,b,c均为常数.例2 求微分积分方程根据傅氏变换的微分性质和积分性质, 且记F x(t)=X(w), F h(t)=H(w).解：在方程两边取傅氏变换, 可得 78运用傅氏变换的线性性质, 微分性质以及积分性质, 可以把线性常系数微分方程转化为代数方程, 通过解代数方程与求傅氏逆变换, 就可以得到此微分方程的解.另外, 傅氏变换还是求解数学物理方程的方法之一798081性质小结: 若F f(t)=F(w), F g(t)=G(w)8283实际上, 只要记住下面四个傅里叶变换, 则所有的傅里叶

14、变换都无须从公式直接推导而从傅里叶变换的性质就可导出.846.卷积与卷积定理称为函数f1(t)与f2(t)的卷积, 记为f1(t)*f2(t)若已知函数f1(t), f2(t), 则积分定义8.385卷积的图示f1(t)f2(t)tOf2(-t)OtOttf2(t-t)86在积分中, 令u=t-t, 则t=t-u, du=-dt, 则即卷积满足交换律.87下证卷积满足结合律, 即f1(t)*f2(t)*f3(t)=f1(t)*f2(t)*f3(t)为此, 令则88交换二重积分的次序, 得令v=t-u, 则u=t-v, 89例1 证明证 根据卷积的定义f1(t)*f2(t)+f3(t)=f1(t

15、)*f2(t)+f1(t)*f3(t)90例2 若求f1(t)*f2(t)f1(t)1OttOf2(t-t)1t91由卷积的定义有tO1-e-t192定理1.2(卷积定理 ) 假定f1(t), f2(t)都满足傅氏积分定理中的条件, 则f1(t) * f2(t) F1(w)F2(w)如 f1(t) F1(w) f2(t)F2(w)以及93证: 按傅氏变换的定义, 有94在物理和工程技术中, 常常会碰到单位脉冲函数. 因为有许多物理现象具有脉冲性质, 如在电学中, 要研究线性电路受具有脉冲性质的电势作用后产生的电流; 在力学中, 要研究机械系统受冲击力作用后的运动情况等. 研究此类问题就会产生我

16、们要介绍的单位脉冲函数.8.3 函数及其傅氏变换95例8.21 在原来电流为零的电路中, 某时间为t=0时刻进入一单位电量的脉冲, 现在要确定电路上的电流i(t). 8.3 函数的概念解：设q(t)表示电路中的电荷函数, 则1. 概念的引入96此外,电路在t=0以后到任何时刻的总电量为97例8.2 对某静止的单位质量物体施以瞬时外力F(t),使其速度v(t)突然增加一个单位. 现在对外力F(t)给以表示. 解：由牛顿第二运动定律知：98而物体在t=0的瞬时动量的增量为99这表明在通常意义下的函数类中找不到一个函数能够表示这样的电流强度. 为了确定这样的电流强度, 引进一称为狄拉克(Dirac)

17、的函数, 简单记成d-函数. 有了这种函数, 对于许多集中于一点或一瞬时的量, 例如点电荷, 点热源, 集中于一点的质量及脉冲技术中的非常窄的脉冲等, 就能够象处理连续分布的量那样, 以统一的方式加以解决.1002. 函数的定义定义1.4 满足以下两个条件的函数称为 函数101的函数称为 (t-t0)函数定义8.5 满足以下两个条件102当趋于0时的极限函数(t)称为 函数定义8.6 函数序列103当趋于0时的极限函数(t- t0)称为 (t- t0)函数定义8.7 函数序列104工程上将d函数称为单位脉冲函数, 可将d函数用一个长度等于1的有向线段表示, 这个线段的长度表示d-函数的积分值,

18、 称为d函数的强度.tOd(t)1tOd(t- t0)1t01058.3.2 d函数的性质1. 筛选性质1061071082. d函数为偶函数将 d函数的定义利用矩形脉冲形式109例8. 25 证明证明：利用筛选性质1101113. 相似性质设a为实常数,则4. d函数是单位阶跃函数的导数1121138.3.3 d函数的傅氏变换114tOd(t)1wOF(w)1可见, 单位脉冲函数d(t)与常数1构成了一傅氏变换对. 同理, d(t-t0)和 亦构成了一个傅氏变换对.115在物理学和工程技术中, 有许多重要函数不满足傅氏积分定理中的绝对可积条件, 即不满足条件例如常数, 符号函数, 单位阶跃函数以及正, 余弦函数等, 然而它们的广义傅氏变换也是存在的, 利用单位脉冲函数及其傅氏变换就可以求出它们的傅