傅里叶变换的性质PPT课件.ppt

《傅里叶变换的性质PPT课件.ppt》由会员分享，可在线阅读，更多相关《傅里叶变换的性质PPT课件.ppt（101页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第第4 4章章 连续连续系系统统的的频频域分析域分析 信号的正交分解与傅里叶信号的正交分解与傅里叶级级数数 信号的信号的频谱频谱傅里叶傅里叶变换变换的性的性质质线线性非性非时变时变系系统统的的频频域分析域分析傅里叶傅里叶变换计变换计算机模算机模拟举拟举例例4.1 信号的正交分解与傅里叶信号的正交分解与傅里叶级级数数 4.1.1信号的正交分解数学上给定条件下的函数可展开为由某种基本函数形式所构成的一组多项式,例如函数的泰勒级数展开式。信号是随时间变化的函数,在一定条件下也可展开成这样一组多项式。这就是信号的分解,用式（41）描述:(i,n为整数)（41）当上述函数集中任意两个函数i(t),j(t

2、)之间,在区间例如,三角函数集1,cost,cos2t,cosmt,sint,sin2t,sinnt,在区间（t0,t0+）(式中T=2/)组成正交函数集,而且是完备的正交函数集。这是因为（ki为与之有关的常量）(42)(43)即三角函数集满足正交性式（42）,因而是正交函数集。其完备性这里不去讨论。对于调幅信号（=5）f(t)=A(1+Bcos)cos(44)利用三角公式2coscos=cos(-)+cos(+)可写为f(t)=Acost+ABcos(-)t+ABcos(+)t(45)式(45)即是信号f(t)在三角函数集上的正交分解。图4.1中绘出了有关信号的波形。图4.1调幅信号及其频谱

3、4.1.2傅里叶级数19世纪初叶,法国数学家吉傅里叶证明:任何正常的周期为T的函数f(t)都可分解为无限个正弦和余弦函数的代数和。即通常称（46）式为傅里叶级数。如果已知f(t),则可通过式(47)、(48)和(49)分别求出an,bn,c的值。(46)(47)(48)(49)根据三角函数的运算法则,式(46)还可写成式(410)。(410)(411)(413)(412)式(46)还可写为如下形式式中,An=A-n,n=-n。最后,由欧拉公式,上式可写为(414)(415)对于式（410）,(414),同式(46)一样,也是傅里叶级数,只是形式不同而已。式(46)和(410)称为三角函数式傅里

4、叶级数,式（414）称为复指数形式的傅里叶级数。由于式（414）的数学表示更为简洁,故在后续章节中,这一式子用得更多。4.1.3信号的傅里叶级数正交分解由于傅里叶级数具有正交性及完备性,故任何周期信号均可正交分解成傅里叶级数。这种分解,在对信号进行分析时将会表现出很大的优势。例41试将图4.2所示的方波信号f(t)展开为傅里叶级数。图4.2方波信号的傅里叶级数解我们将信号按式(46)分解成傅里叶级数,并按式(47)、(48)、(49)分别计算an,bn及c。4.2 信号的信号的频谱频谱4.2.1信号频谱上一节我们指出,信号可分解为傅里叶级数,即信号可由系列复数指数函数加权之和构成。一般我们称这

5、里的复数指数函数ejnt为n次谐波,在该函数上所加的权为谐波的振幅,n为谐波的角频率,可以说所有的信号均是由系列角频率不同的谐波叠加而成的(角频率可简称为频率）。4.2.2周期信号的频谱以周期性矩形脉冲为例,说明周期信号频谱的特点。设有一幅度为1,脉冲宽度为的周期性矩形脉冲,其周期为T,如图4.3所示。根据式（46）,可求得其傅里叶系数图4.3矩形脉冲考虑到=2/T,上式也可以写为根据式（416）可写出该周期性矩形脉冲的指数形式傅里叶级数展开式为图4.4画出了T=4的周期性矩形脉冲的频谱。由于Fn为实数,相位n=0,故而没有单独画出其相位频谱。图4.4周期矩形脉冲的频谱（T=4）1.频谱的物理

6、意义前面讲过,任何信号均由多次谐波叠加而成,我们通过仪器观察谐波时,只有由三角函数所描述的谐波Akcos(kt+k)才能被观察到,而复指数谐波ckejkt是通过数学方法由前者构造而成,它不能直接被观察得到。两者的关系为即有(418)2.频带宽度从周期矩形脉冲频谱可以看出,谱线有无限多条。矩形脉冲信号的频带宽度或称信号的带宽,用符号f表示,即3.周期信号的功率了解周期信号功率在各次谐波中的分布情况,是信号频谱的一个重要应用。分析信号的功率关系,一般都将信号f(t)看作电压或电流,而考察其在1电阻上所消耗的平均功率,即（419）(420)将f(t)表示成傅里叶级数并代入上式可得(420)（421）

7、4.2.3非周期信号的频谱非周期信号可视为周期足够长的周期信号来处理。因此,我们可以从周期信号的频谱分析来推测非周期信号的频谱。重写周期信号的频谱函数如下:(422)(423)(424)(425)现将信号f(t)的傅里叶级数展开式重写如下将式(425)代入式(426)中,同时将求和号改为积分号,n改为,则有(426)(427)式(424)和(427)是非常重要的一对式子,重写如下,并称前式为f(t)的傅里叶变换,后式为函数F()的傅里叶逆变换,F()称为f(t)的频谱函数,f(t)称为F()的原函数。(428)(429)式(428)和(429)可简记为(431)(432)4.2.4常见信号的频

8、谱分析举例例42求冲激信号(t)的频谱。解由频谱函数的定义式（428）有(434)(435)图4.5冲激信号及其频谱例43求矩形脉冲信号g(t)的频谱。图4.6矩形脉冲信号及其频谱解矩形脉冲信号g(t)是一个如图4.6（a）所示的门函数。其定义为(436)g(t)的傅里叶变换为(437)(438)(439)例44求单边指数信号的频谱。解单边指数信号是指(441)(440)图4.7单边指数信号及其频谱例45求双边指数信号的频谱。解双边指数信号是指(442)从频谱函数的定义式出发(443)图4.8双边指数信号及其频谱例46求单位直流信号的频谱。解幅度为1的单位直流信号可表示为f(t)=1,-t0)

9、(4-51)4.3 傅里叶傅里叶变换变换的性的性质质为了方便起见,我们将傅里叶变换式重写如下(4-52)(4-53)(4-54)4.3.1线性某一域内函数的值作线性变换时与之对应的另一域中的象函数的值也作等比例的线性变换。例48求单位阶跃函数u(t)的频谱函数。解单位阶跃函数u(t)可看作是幅度为1/2的直流信号与幅度为1/2的符号函数sgn(t)之和,即4.3.2奇偶虚实性实际存在的信号都是实信号,虚信号是我们为了数学运算上的方便而引入的。现在研究时间函数f(t)与其频谱F()之间的奇偶虚实关系。先来看f(t)为实函数的情况。此时傅里叶变换可写为（460）式中,频谱函数的实部和虚部分别为(4

10、61)(462)频谱函数的模和相角分别为(463)(464)由此可看出,此时F()是虚函数且是的奇函数。对于f(t)为虚函数的情况,分析方法同上,结论相反。上述讨论的结果如下:例49利用奇偶虚实性求图4.11单边指数信号f(t)=2e-tu(t)的频谱。图4.11单边指数信号及其频谱解从波形图（a）上可见,单边指数信号f(t)是非偶非奇函数,但可分解为如图（b）,（c）所示的偶函数和奇函数两部分,见下式。f(t)=2e-tu(t)=fe(t)+fo(t)其中4.3.3对称性傅立叶变换可用（4-52）表示图4.12抽样函数Sa()及其频谱4.3.4尺度变换将时间函数f(t)中的t换成at(a为常

11、量),考察与之对应的频谱函数。现在来求时间函数f(t)尺度变换后的频谱函数。设f(at)=f(t),则有（4-72）例411已知求g(2t)的频谱函数。解根据傅里叶变换的尺度变换性质,g(2t)的频谱函数为图4.13尺度变换4.3.5时移特性在时间函数f(t)中,当时间t变为t+t0时,就会引起相应的频谱函数的变换,称为时移特性。这里t0为实常量。设f(t+t0)=f1(t),时移特性可作如下推导:因此,时移特性可表示如下。若且t0为常数,则(474)例412求移位冲激函数(t-t0)的频谱函数。解由于已知冲激函数(t)的频谱函数为1,求移位冲激函数(t-t0)的频谱函数，此时可利用傅里叶变换

12、的时移特性式(474)。(475)4.3.6频移特性频移特性与时移特性对称,此时所考虑的是频谱函数F()中,频率变为+0,相应的时间函数怎样随之而变。这里0为实常量。因此,时移特性可简写如下且0为实常数,则例413求高频脉冲信号p(t)=g(t)cos0t(477)的频谱函数。解由于(476)故有根据频移特性有图4.14频移特性4.3.7卷积定理现在我们讨论卷积在傅里叶变换中的规律。1.时域卷积性质设有两个时间函数f1(t)和f2(t),它们分别对应的频谱函数为F1()和F2()。两个函数卷积的傅里叶变换为(4-78)(4-79)这就是时域卷积定律，可简记为（4-80）2.频域卷积性质同时域卷

13、积定律一样，我们也可以证明频域卷积定律。在这里略去证明，只写出结论。(4-81)(4-82)例414求图4.15所示梯形脉冲的傅里叶变换。图4.15梯形脉冲的傅里叶变换解梯形脉冲可看作是两个不同宽度的矩形脉冲f1(t)与f2(t)的卷积，如图4.15所示。f(t)=f1(t)*f2(t)而矩形脉冲的傅里叶变换已在例43中求出,具体来说图4.16半波正弦脉冲图4.17三角形脉冲及其一、二街导的波形表42傅里叶变换的性质4.4 线线性非性非时变时变系系统统的的频频域分析域分析4.4.1频域分析在第二章线性非时变系统的时域分析中,我们已经指出线性非时变系统的零状态响应yf(t)是激励f(t)与冲击响

14、应h(t)的卷积积分。即(4112)(4113)现在我们设h(t)、f(t)和yf(t)各自对应的傅里叶变换式为H()、F()和Yf(),即(4114)(4115)(4117)(4116)图4.18频域分析示意图接下来我们讨论一下H()。很显然在前面的讨论中,H()已有了两个方面的含义:一是前面我们在式（4114）中定义的H()是与冲激响应h(t)对应的频谱函数;二是通过式（4117）所表达的H()的含义是零状态响应频谱函数Yf()与激励函数F()的比值,即(4118)另外,考虑到式（4115）在看式（4-116）和（4-117）（4-119）（4-200）例420如图4.19所示,试分析单位

15、阶跃信号u(t)通过RC高通网络传输后的波形。图4.19解显然,当输入信号uS(t)为复指数信号ejt时,如图有则按H()的定义有对于单位阶跃信号u(t)而言,此时最后一步考虑了冲激函数的取样性质。因此4.4.2无失真传输系统的频域分析无失真传输系统是指这样一个系统,它的输出信号与输入信号相比只有幅度的大小和出现时间的先后不同,而没有波形上的变化。系统的数学模型具有如下的形式:yf(t)=Kf(t+td)(4121)设输入信号f(t)的谱频函数为F(),输出信号yf(t)的谱频函数为Yf()。根据傅里叶变换的时移特性对上式进行傅里叶变换后可得输出信号频谱和输入频谱之间的关系:Y()=KejtF

16、()(4122)由上式可见,为使信号传输无失真,系统频率函数应为H()=Ke-jtd从中可看出系统频率函数的模和相位为(4123)图4.20无失真传输系统的幅频和相频特性其实，此时的系统输出信号yf(t)既是系统的冲激响应h(t),既有（2-124）（4-125）4.4.3理想低通滤波器的频域分析理想低通滤波器是指频率特性为式(4126)所限制的系统，即由图4.21中可看出H()可看作是在频域中宽度为2C,幅度为1的门函数,可写为(4126)(4127)图4.21理想低通滤波器的频率特性图4.22理想低通滤波器对单位冲激信号的响应波形则输出信号的频谱为由此可求出理想底通滤波器对单位阶跃信号的影

17、响（4-128）图4.23理想低通滤波器对单位阶跃信号的响应波以上讨论了理想低通滤波器对单位冲激信号和单位阶跃信号的响应,这里我们还需要注意以下几点:(1)由响应的波形图可见,响应的时间比激励滞后,延迟时间为td。(2)阶跃信号的响应不像阶跃信号那样陡直,而是倾斜的，这说明输出信号的建立需要一定的时间。一般以阶跃响应中幅度由0到1作为计算建立时间的标准。查Six正弦函数积分表可知响应建立时间为(3)由响应的波形图可见,输出信号在输入信号建立之前和后都有,向延伸且振荡。由此,早在t=0时刻以前在无信号输入的情况下就已有信号输出,这显然违背了自然界的因果律。这是因为理想低通滤波器是根据式（4127

18、）设计的,过于理想化,现实中不可能实现。4.5 傅里叶傅里叶变换计变换计算机模算机模拟举拟举例例运用计算机对傅里叶变换进行模拟的主要任务就是在设定系统频率函数的基础上根据不同的激励信号来模拟与之相应的响应。整个过程一般分为以下几个步骤:（1）输入系统频率函数;（2）输入激励信号;（3）按一定的算法完成激励信号的傅里叶变换得到激励信号频谱;（4）根据系统频率函数和激励信号频谱按一定的算法计算系统的响应频谱;（5）按一定的算法完成频率响应的傅里叶逆变换得到系统响应;（6）输出系统响应及其频谱。(4129)(4130)(4131)我们在每个分区内用梯形面积代替实际积分值。当分区相当多时,所有这些梯形

19、面积的和可逼近实际的积分值。计算公式如下:当时运算结束。Ri,j为积分近似解。龙贝格积分的计算次序是依上述算法实现的相应C语言算法如下:floatrbg(float(*f)(),floata,floatb)floath,sum,R101101,err;inti,j,k,n;h=(b-a);R11=h*(f(a)+f(b)/2;n=100;for(i=1;i=n;i+)sum=0;for(k=1;k=2*i-2;k+)sum=sum+f(a+(k-0.5)*h)h=h/2;Ri1=0.5*(Ri-11+h*sum);for(j=2;j=i;j+)err=(Rij-1-Ri-1j-1)/(powe

20、r(4,j-1)-1);Rij=Rij-1+err;err=(Rij-1-Ri-1j-1)/(power(4,j-1)-1);Rij=Rij-1+err;if(abs(err)E)thenreturnRijreturn0例421计算抽样函数的傅里叶变换。解式(4132)抽样函数的傅里叶变换可写为(4132)这样我们可利用龙贝格算法求出该函数的数值解,程序如下:floatt;floatD1(floatw)returnsin(t)*cos(w*t)/t;floatS1(floatw)returnsin(t)*sin(w*t)/t;main()inti,a,b;a=LEFTVALUE;b=RIGHTVALUE;dt=STEP;t=STARTVALUE;for(i=1;in;i+)Di=rbg(D1,a,b);Si=rbg(S1,a,b);t=t+dt;