1、傅立叶级数和傅立叶变换解析

1、傅立叶法国数学家及物理学家。 最早使用定积分符号，改进符号法则及根数判别方法。 傅立叶级数（三角级数）创始人。法国数学家、物理学家。1768年3月21日生于欧塞尔，1830年5月16日卒于巴黎。9岁父母双亡，被当地教堂收养。12岁由一主教送入地方军事 学校读书。17岁 （1785）回乡教数学，1794到巴 黎，成为高等师范学校的首批 学员， 次年到巴黎综合工科学校执教。1798年随拿破仑远征埃及时任军中文书 和埃及研究院秘书，1801年回国后任伊泽尔 省地方长官。1817年当选为科学院 院 士，1822年任该院终身秘书，后又任法兰西学院终身秘书和理工科大学校务 委员会主席。主要 贡献是在研究热

2、的传播时创立了一套数学理论。1807年向巴黎科 学院呈交 热的传播论文， 推导 出着名的热传导方程 ，并在求解该方程时 发现解函数可以由三角函数构成的级数形式表示， 从而提出任一函数都可以展成 三角函数的无穷级数。1822年在代表作热的分析理论中解 决了热在非均匀加热的 固体 中分布传播问题，成为分析学在物理中应用的最早例证之一，对19世纪数学和 理论物理学的发展产生深远影响。傅立叶级数（即三角级数）、傅立叶分析等理 论 均由此创始。 其他贡献有： 最早使用定积分符号， 改进了代数方 程符号法则 的证法和实根个数 的判别法等。3-2 信号在正交函数集中的分解为了形象地说明信号的分解，首先我们复

3、习矢量的分解。一、矢量的分解（1）矢量的一维分解：用一个标准矢量Ai乘以一个标量&得到的新矢量Ai，去近似近似矢量 A , 并要求误差尽可能小，应该取多少？下图通过几何方法表示了Cl的确定方法。从几何或者解析角度，都可以得到使误差最小的系数为:A1AC1二A1A1其中的C1称为矢量 A 和 A 的相似系数。如果C1=0（或A A1A A二0），则表明A和A A1相垂直（又称为正交）。（2）矢量的二维分解用两个标准矢量A Ai、A A2的线性组合ciA Aic2A A2，去近似近似矢量 A，并要求误差尽可能小，C1、C2各应该取多少？下图通过几何方法表示了 &、C2的确 定方法。

4、在上图表示的情况下，G G、C C2的取值都同时与A1、A2有关，计算公式 可能比较复杂。如果标准向量A A1、A A2相互垂直（正交），计算就很简 单了：1容易得到此时的系数计算公式为:A1Ac1:AiAi此时每一个系数只与其相关的标准矢量有关， 系数计算公式与一维情况下的 计算公式相似。上图中表示的是用两个矢量表示一个二维的矢量，误差为零。如果用两 个矢量表示一个二维以上的矢量，误差就不一定等于0了。但是可以证 明，在这种系数情况下误差最小。显然，如果知道了标准矢量A Ai、A A2和相应的系数 G G、Q Q，就可以确定任意矢量。这实际上就是我们在平面几何中见到的笛卡尔坐标系。（3）矢量

5、的多维分解：上面二维的情况可以推广到任意维，可以将矢量表示成为一系列标准矢量 （基）的线性组合：nA ClA1C2A2CnAn二CiAii二显然，如果知道了标准矢量A Ai和响应的系数C Ci,就可以确定任意矢量。 如果矢量A Ai两两正交，可以证明相应的最佳系数的计算公式为：如果标准矢量基A Ai的长度都为1，则A AiA Ai=1，上面的公式可以简化为:Ci= AiAC21（4）标准矢量基的几个限制条件:为了便于计算系数，实际使用的标准正交矢量集最好满足以下几个条件：1）归一化：标准矢量Ai的模等于12）正交化：标准矢量两两正交3）完备性：可以不失真地组合出任意矢量其中归一化和正交化是为了

6、计算系数时比较方便；而完备性则是为了保证可 以完整、没有误差地表示任意矢量，使这种分解更有实用性。二、信号的分解与矢量分解相似，我们也可以推导出信号分解。1、单个标准信号下的分解：在时间区间魚1）内，用qfi（t）近似任意函数f（t），并使误差进可能小。（这里假设所有函数都是实数函数）误差：；(t)二f(t) Gf1(t)如何衡量函数误差的大小?-可以采用方均误差：1t2；二t；(t)dt t2 711取什么值的时侯何时误差最小？或者何时系数最佳?f (t) f (t)dt G二t；f1(t) f1(t)dtt1G也称为函数f（t）和f1（t）的相似系数。最佳系数的证明:误差：；(t)二f(t

7、) Gfjt)方均误差:最佳系数:ti21；(t)：t2-tl1t?-ti 1-t2为了求使；2(t)最小的ci，将上式对ci求偏导并令其为零，可以得到:-t2t22 n厂2f f(t)fi(t)dt +2ci (fi(t)dt=O由此可得:f f(t)fi(t)dtTiG = -12(fi(t)fi(t)dtf (t) f|(t)dt = 0如果Ci(或ti)，则称f(t)和fi(t)正交-这个正交的含义与矢量中的正交类似。如果f(t)和fi(t)是复函数，则方均误差的定义应该改为：2it22it2*孑化)=卜(t)欽t)(t)dtttitittiti相应的最佳系数计算公式为：t2*tf(t

8、)fi(t)dtTiG =*tfi(t)fi(t)dtti2、多个标准信号下的分解：将信号表示为多个标准信号的线性组合：nf (t) “ fi(t) C2f2(t) . Cnfn(t)Cjfj(t)i Tt22r(t)dtt1f f(t) - cJi(t)dt:f2(t)dt - 2ci：f (t) fi(t)dtCi2Jf(t)dt2t2这里的同样难以确定。但是如果标准函数fi之间两两正交, 证明：f f(t)fi*(t)dtLit2*fi(t)fi(t)dt则可以我们实际上在高等数学等前期课程中已经见到过几个这样的标准信号集了。例如：泰勒级数使用的是：1,x,.,.?在本章中将要用到的标准

9、函数集为三角函数集：1,cost,sint,cos2t,sin 2t,., coskt,sin kt,.”3、对标准信号集的要求：与矢量分解中的情况一样，这里对于用于分解函数的标准函数集也有以下的 要求：3）完备性：可以用其线性组合表示任意信号正交性标准函数集的首要条件。只有在这种情况下系数才可以用上 美的公式计算，而且可以保证方均误差最小。其他两个条件都会受 到实际应用的限制，可能难于达到。完备正交函数集一般都包含无穷多个函数，例如： 三角函数集， 沃 尔什函数集等。但在实际应用中不可能用无穷多个，只可能用有限个函数，只能近 似表示任意函数。函数与矢量的运算与分解有很大的相似性，很多函数分解

10、中的概念（例如正 交等）也是从矢量运算中引用过来的。这里用一个表格作比较：矢量函数加法A1+ A2f1(t)+ f2(t)标乘c A2c f1(t)乘法A1A2=A1A2cosrf1(t) f2*(t)dtt1正交A1A2= 0占*fjt) t (t)dt = 0t11）归2）正交化:t1fi(t) fj*(t)dt =0f(t) f：(t)dt =1t2归一A= 1t2 *f(t) f (t)dt=1t1误差名=A1- A2哄t)= f1(tf2(t)误差代价函数忖221t22名2(t)=- /(t)dtt tt1L2L1系数A A1c1=A1A1f2f(t)f1*(t)dtt15 一t2*

11、! f1(t)f1(t)dtt1 3 -3 信号表示为傅利叶级数傅利叶级数是最常用的一种正交函数集。它在工程中有很广泛的用途。一、三角函数形式的傅利叶级数1、 三角正交函数集Hcos tsin 4, cos2t,sin2 4,.,cosk 4,sin2兀2兀12 =- =-其中：Tt2 _ti或将正交函数集表示为：Zos(nt),s in (nt) n =0,1,2,.；可以证明该函数集满足正交性：函数集中的函数两两相正交。t2cos( nt) si n(mt)dt =0t1t2f cos(nOt)cos(m0t)dt =0m式nL2f sin(nOt)sin(m0t)dt=O丫1J2、 任意

12、信号在三角函数集中的分解可以将任意函数f(t)在这个三角函数集中展开(表示成该正交函数集函数 的线性组合)：f(t)=玄ycosjt a2cos2 4 . ancosnt .dsinxt b2sin2t . bnsinn-.4 .-ho二a0 ancosni】t bnsinnn 4其中的系数可以根据前面的公式计算出:这个公式中的an的表达不太方便。为此将分解式改写:-bef(t)二色亠 _iancos nit bns inn 4 2n二则系数为：anf(t)cos(nt)dtt2 7屯2t2bnf (t)sin( n4)dtt2 7 %通过这种分解，可以将信号可以表示成为直流信号和一系列正弦信

13、号之 和。3、任意信号在仅余弦三角函数集中的分解在原来的信号分解公式a乂f (t)0亠 一i.ancos nt bnsin nt2心中，利用三角函数公式，令K$ 22叫=-arctan二An=n +0a.则可以将上式表达成:an =2f(t)cos(nt)dttt2f (t)cos(ni】t)dt n = 0t1f cos2(n Ct) dtt1t2t2jf(t)dtt2j f (t)sin(nQt)dt1bn=t2sin2(nt)dtt1t2f (t) sin(nCt)dtt1f (t)一 AnCOS n：n2n 4它可以看成是下列正交信号集：匕os(nt)n -0,1,2,./的平移后的线

14、性组合。从系数计算公式可以看出，如果f(t)是实数信号，贝an和A是n的偶函数；bl和是n的奇函数。上面的分解等式的左右两边的函数是否相等，没有误差？或者,是否随着n趋向于无穷大，等式右边的函数收敛于左边的函数？- Direchlet证明，只要满足下面三个条件，等式就一定收敛：Cf (t) dt =1)f(t)绝对可积，即：I2)f(t)在区间内有有限个间断点；3)f(t)在区间内有有限个极值点。这个条件被称为Direchlet条件。实际信号大都满足这个条件，所以都可以 这样分解。这个分解等式中，等号右边是多个周期为T的函数的和，它仍然是周期为T的函数。显然，如果f(t)本身也是一个周期为 T

15、 的函数，则如果它可以在一个周 期内用上面的公式分解，则它同时也可以在整个时间区间内分解。这种分解可以用在两个场合：1)研究任意函数在(t1,t2)区间内的分解2)研究周期为T的函数在整个时间区间内的分解。本课程中讨论的主要是后一种情况。如果f(t)周期为T的函数,为了方便讨论，一般函数的主值区间取T门 22.丿在函数的分解中:a0兀称为信号的直流分量;aiCOSF、aisin沁或A cos（t 1）称为信号的基波分量；Kcos zt、Ksi rm川或A1cos（n4n）称为信号的n次谐波分量;一般情况下，n无法计算到无穷大，只能取有限。这时，这种正交展开是有 误差的。n越大，误差越小。下面通

16、过一个实例进一步讨论傅里叶级数的一些特性。例：求方波f（t）=*的傅利叶1-1十级数。0 ct2Tct cT 2/anbn.an2- bn2cos(申n)=：n2sin(n) =an *bnAn二an- jbn=2f f (t)cos(n0t)dt j2f f(t)sin(n0t)dt T JT叫f (t) Cos(rrt) - j sin(ni】t) dt f (t)eEdt两种推导过程得到的答案应该相同。对比两个系数计算公式，可以得到:_ An _ Aej _an - jbncn2 2 2这个等式反映了Cn与A、n或an、bn之间的关系3-4 周期性信号的频谱周期性函数可以在傅利叶级数中展

17、开。 如果给定了各个频率分量的幅度 和相位，就可以确定信号。频谱是信号的一种图形表示方法，它将信号各个频率分量上的系数关系 用图形的方法表示出来。它可以说明信号的特性，而且可以给信号的变 换和处理计算带来很多方便之处。频谱图有两个组成部分：振幅频谱：表示信号含有的各个频率分量的幅度。 其横坐标为频率, 纵坐标各个对应频率分量的幅度。相位频谱：表示信号含有的各个频率分量的相位。 其横坐标为频率, 纵坐标各个对应频率分量的相位。频谱图有两种形式：1、如果用正弦函数展开式形式的傅里叶级数，则相应的表达式为：bof (t)二西、AnCOs nti：n2n则振幅为：/ 22An_,an 0相位为：二ar

18、cta n按照这种定义做出的频谱，因为只有n（或-0）时才有意义，做出的图只有n亠0的一边，所以又被称为单边频谱。例：周期性方波的单边频谱。4bn = *| n兀an =1所以:22.-4当n为偶数An =、anbn= n二0当n为奇数，广bn -=当n为偶数屮n= _arcta n=丿2an0当n为奇数由此可以作出其频谱图2、如果用复数正弦函数展开式形式的傅里叶级数，则相应的表达式为：f （t）=送CnejnnR的傅里叶级数表达式，贝U：振幅为Cn相位为an gG）。按照这种定义做出的频谱在n大于和小于零的两边都有意义，做出的图 又被称为双边频谱。当n为偶数当n为奇数由于对于实数信号而言，其

19、频谱具有对称性，所以一般情况下对于双边频谱也只要作出门一（或-0）部分就可以了。这样一来的频谱与单边 的频谱就有些相似，但是含义不同。在频谱形状上，两者的相位频谱相 同，但是振幅频谱的幅度大小是单边谱的一半。单边频谱在物理概念上容易理解，但是双边频谱对于后续的处理带来很 大的好处在后面的内容中，频谱往往都是用双边频谱。单边频谱中，对于n=（或，=）点上的幅度频谱，有一些与其它频率 点上的不同之处：1）如果认为幅度频谱表示的是是信号在各个频率上的信号分量幅度的A0大小，则信号真正的直流分量应该为2,频谱在=上的分量的 大小应该减半。2）如果认为幅度频谱表示的是A1随频率变化的规律，则幅度频谱不用

20、 变化。周期性信号的频谱有下面三个特点：1.离散性：它有不连续的线条组成；2.谐波性：线条只出现在基波频率的整数倍点上；3收敛性：实际信号的幅频特性总是随频率趋向无穷大而趋向于零。例：周期性方波脉冲的频谱：kT : t kT2 2其它根据上面的公式可以画出信号的频谱。 该例中信号的振幅频谱和相位频谱可 以合二为一f(t)AiCn4A .- sinn0T 2AT Tz宀a n = T2A.nM 1-sin -nOTI 2丿A弋L Tn ASac Tn = 0nJ2Sax壬其中：x，称为抽样函数。根据周期性方波的频谱，我们可以得到关于信号特性的几个一般性结论：1、T增加一一Sa()函数不变一一频谱

21、的包络不变，收敛性不变。但是：1) 谱线幅度降低；2)谱线密度加大。信号周期加大，对振幅的收敛性没有影响，但会使谱线密度增加。当T趋向无穷大时，信号成为非周期信号，这时，谱线幅度降低为无穷 小,谱线密度加大，信号分量出现在所有频率上。2、下降一一Sa()尺度扩大一一收敛性变差，但是谱线间隔不变。信号时间宽度变小，将使信号能量向高频扩散，信号的频带增加。3、 信号的频带：由于信号的频谱的收敛性，一般可以在一个信号分量主要集中的频率区间内 研究信号的特性，而忽略信号其它部分的分量。响应的频率区间就是信号的频带。信号的频带有很多种定义方法：1)以信号最大幅度的1/10为限，其它部分忽略不计；2)以信

22、号振幅频谱中的第一个过零点为限，零点以外部分忽略不计;3)以包含信号总能量的90%处为限，其余部分忽略不计；4、信号的边沿对信号频带的影响 信号的边沿变化越快，信号的频带越宽 例：三角脉冲函数的频谱：AT汽4 Jj(t_rec_p.m)3-5非周期性信号的频谱非周期性信号可以看成周期信号在周期趋向无穷大时的极限一、从周期信号到非周期信号从傅利叶级数到傅利叶变换根据周期信号傅利叶级数展开公式，其各个频率分量的幅度为:2兀0T，此时:1)频谱间隔趋进无穷小，信号在各个频率点上都有信号分量一一频率取值变成连续的1t2= f(t)eCn2)在每一个频率点上的频率分量大小趋向零。其中第二点给计算带来了麻

23、烦，所以无法用傅利叶级数表示非周期信号。这时，为了消除系数公式中趋向无穷小的部分，定义：ct?F(n二TCn( = 2n) = f(t)entdtQ11这时上式可以得到一个非零的值。令T一；心，则0，而成为一个连续的变量，假设其表示为连续的变量-，则可以得到傅利叶变换公式：F(j ) f(t)edt因为该式有“单位频带内信号幅度”的量纲，所以被称为“频谱密度函 数”。它表示信号在该频率点上的分量的相对大小， 而信号在此频率点 上的实际分量分量大小为零。与傅利叶级数一样，如果f(t)是实数函数，F(j)的幅度是的偶函数，F( j* )的相位是的奇函数。、傅利叶反变换一一怎样用F(j )计算f(t

24、)jni tF( j n； ；. )jn i tf (t)二lim % Cnelim、eTT_jpcTn n 亠11i tF(j )ejd 2二这个公式实际上也表示了将信号分解为一系列复数三角函数的子信号之和(积分)。这个公式也可以表达成为一个在物理上更容易理解的实数三角函数形式：三、正反傅利叶变换由此可以得到正反傅利叶变换公式为：FT. F3)=Ff(t)=J；f(t)e*dtf (t) = F亠忙(j叭鼻亠J：F (j灼)ej%ccIFT:2二-f(t)和F ()之间是一一对应的，根据其中的一个可以确定另外一个。 可以认为，它们包含了相同的信息，只不过自变量不同，它们是相同信 号的不同表达

25、形式。正变换将以时间为自变量的函数变成了以频率为变量的函数，将信号从时域变换到了频域。所以建立在这种变换上的系统 分析方法称为变换域法。这种变换通常经过积分计算得出，所以也称为 积分变换。傅利叶变换所牵涉的两个函数都是连续函数， 所以它完成的是从连续函 数到连续函数的变换；而傅利叶级数则是完成从连续函数到离散函数的 变换。傅利叶变换存在的条件依然是Direchlet条件，只不过这时考虑的时间这里，在频域中我们用j作自变量，目的是为后面引入拉普拉斯变换打下伏笔四、非周期信号的频谱这里同样可以用图的形式，在变换域中表示信号。响应的频谱图称为信号的 幅频特性曲线和相频特性曲线。五、傅利叶变换的另外几

26、种形式：IFT:这种形式上正反傅利叶变换形式上比较对称。但是使用时并不方便。2、一些文献上也可以见到另一种形式的傅利叶变换公式:FT. F(2)=Ff(t) = 2H：f(t)edt f(t) =F亠9(浮)=(乍(浮)e血d国IFT : -6 常用信号的 F.T常用信号的FT见P125129表。现在将一些结论列举如1、将频域中的自变量从变成，则:FT:-beF(j2f)二J(t)e_2ftdt1亠.f(t)F(j2：f)ejftd 2二fIFT: 或:2:二F(j2：f)ej2ftdfF(f)_：f(t)ejdtf(t) = _.：；F(f)ej2ftdfFT:区间为以上两个信号的FT只在:

27、0时存在4、门函数：G (tKSa(f)名(t)兀石() +5、阶跃信号：厂6、直流：12二)阶跃信号和直流信号并不满足绝对可积条件，严格地说不 存在傅里叶变换。但是通过引入冲激函数，也可以找到其 傅里叶变换的表达式，从而也可以用傅里叶变换的方法进 行分析；3-7 周期性信号的傅利叶变换周期信号只是一个相对的概念。如果忽略其周期性，它应该也可以被看成是 非周期信号处理，进行傅里叶变换。但周期信号是功率信号，不满足绝对可积条 件。但是通过引入冲激函数，一样可以找到傅里叶变换。1、复正弦信号的傅里叶变换：el 2二c)根据这个变换以及后面要证明的傅里叶变换的线性特性，可以推导出：cos(ct) I

28、 菽(:c)、G；-：c)丨2、周期性信号的傅里叶变换周期性信号可以展开成傅里叶级数：f (t)= _討沁n = ：:由此可以得到周期性信号的傅里叶变换为下:1、2、3、冲激函数：川八1单边指数信号:双边指数信号:e；(t)，-boF（ j J = 2二Cn、.（，-n）n二 3可见，周期性信号的傅里叶变换是一系列间隔均匀的冲激序列。3、脉冲信号f（t）（FT为F（j）按照周期T进行周期化后信号 匸的FT:（这里假设周期化后各个脉冲没有重叠）f（t）周期化后可以表示成为傅利叶级数：命=Cn川1n =:所以：-bofT（t） i 2二vCn（，-n1）1n二二:其中：Cn二丄匚仃“也二丄,（t）

29、ejdt二F（jn）T=T所以：fT（th -F（jn）、（ 一n）】Tnjoo-boF （jn1）、（ -n1）1n =.：:通过查表，可以很方便地得到：1）非周期信号的FT2）周期信号的FT3）周期信号的傅利叶级数对照傅利叶级数和傅利叶变换的定义，可以得到：3-8傅利叶变换的性质1、 线性特性：a （t） b f-（t），a F!（j ） b F-（j ）2、 延时特性：f（t to），F（）ej t03、 移频特性：f （t）ej cFj（ 一c）1移频特性与延时特性互成对偶。推论：1d.nf(t)cos( ct)Fj(门亠 Lc) h：；Fj (；-：；%) P24、尺度变换：f(at

30、)1 2F j-a 1-信号的宽度沿时间轴压缩-倍，信号的频率宽度B沿频 率轴扩展-倍。脉冲信号的宽度和频带宽度B的乘积等 于常数。数据传输中总希望信号的脉冲宽度尽可能小，占用的信号频带同时也尽可能小。但从该性质可以看出，信号脉冲宽度 的频带宽度是一对矛盾。5、奇偶虚实性假设：F(j )二.；：f(t)ej也= f (t)cos，t dt j f (t) sin t dt二R( ) - jX ()二F(j)其中：Rd) =Uf(t)cotdt，为Fj)的实部;X(沪：f(t)sintdt，为F()的虚部;F(冋，为F(冋的幅度;1-f(-t)F(-)b、f*( F*(-j )c、f*( F*(

31、j )2如果信号f(t)是实数信号，贝U:aR(j是的偶函数；X()是的奇函数; 或：F(j F*(-j )()=-rct-nlR)丿,为F(2)的相角;b、T（j如 是的偶函数；申）是的奇函数;3） 如果f（t）是实偶函数，则F（）也是实偶函数;如果f（t）是实奇函数，则F（j）是虚奇函数；4） 思考：如果f（t）是虚函数，情况怎样？5） 对称特性如果f（t）i F（j ），则:F（jt） 2计（- ）7、微分特性&积分特性t1.f( ) -F (0)、( )F(j )-j 1F(j ) j 9、频域微积分d-jtf （t）F（j ）do二f（0）、（t）jp，F（r）dt一-（1）

32、灼:;（t） j- f（t）_F（r）d门it）或：、n0(t) jj f(t).；F(j)d.dndnF(j )如果dtf(t)存在并且满足Direchlet条件，则:dn推广:dtnf(t)，j 、(j )lim如果F（0） =0，或畀F（ ） 存tJ(.)d 10、卷积定理fl（t）* f2（t），Fi（j ）F2（j ）1fi（t）f2（t）Fi（j .）*F2（j .）2二利用这十个性质，结合傅利叶变换表，可以求解很多工程 上的信号的傅利叶变换。-9 能量频谱与功率频谱能量频谱和功率频谱从能量或功率的角度研究信号在各个频率分量上的能 量或功率，以频谱的形式表达出。这种频谱对确定性信号

33、意义不大， 对于随机信 号有很大意义。但为了方便讨论，这里我们从确定性信号的角度进行研究。一、周期性信号的功率谱周期性信号的能量无穷大，无法从能量上进行研究。但是它的功率有限, 可以从功率上进行研究。1、周期性信号的功率谱：将周期性信号在各个频率上的分量的功率大小用图的方法表示出。横坐标: 频率；纵坐标：信号分量的功率。对于单边功率谱，在每个不等于零（非直流）频率上子信号功率1212AA2。直流信号的功率为4对于双边功率谱，在每个频率点上子信号功率2=Cn功率谱只有大小（幅度），没有相位。2、Parseval定理：周期信号的功率等于该信号在完备正交函数集中分解后各个子信号功 率的和。二、 能量

34、信号（脉冲信号）的能量谱1、能量谱1） 能量信号的功率为零，能量为有限，只可以从能量角度研 究其分布；2） 信号在各个频谱上的实际分量大小为无穷小， 只能用能量 密度谱描述G（）描述单位频带内的信号能量。信号总能量：:2: *W f (t)dt f(t)f (t)dtT，(j). jtLdtd_(r )F*(r )d -He21址2fJF()血=-j0 F(j砒do由此定义单位频带内信号的(2)如果信号是实数信号，则还可以得到其单边能量谱为:“单位频带”指的是什么频率:a般情况下指角频率b、也可以用一般频率f(单位Hz)2由此可以得到双边能量谱：G(f) =F(j2f和单边能量谱：G(f2F(

35、j2Tf)能量谱同样只有大小(幅度)，没有相位2、Rayleigh定理：1 W f(t)dt&12n1 :.?Wf(t) dtj0即：信号在时域和频域的能量相等脉冲信号的脉冲宽度和频带宽度对于一般的信号，可以通过其频谱密度函数或功率谱函数定义其频带宽度, 其定义方法与3.4节中讨论的相似F(j )ej td dt*F(j ) f (t)ejt.dtdjt.dtd(1)双边能量谱为:1G()二2F(jo)G( ) =1JtF(j )2此时:咨.2W= J*F(j2兀f). df2I 2j) d(0f)df- 2F(jw) d国=2 J F( j2斗)df12二11、 脉冲宽度：脉冲的绝大

36、部分能量集中的时间区间to .2(f(t) dt=W2、 频带宽度：脉冲的绝大部分能量集中的频率区间丄F(念)2血=HWn 03、 对于一种脉冲而言，B =常数Matlab 函数1、产生信号sincos sinesquarechirp(t,fO,t1,f1)x = square(t)产生周期为2 n的方波信号x = square(t,duty)，占空比为duty的方波x = sawtooth(t)产生周期为2 n的锯齿波或三角波x = sawtooth(t,width) 2、FFT和ifft,fftshift,ifftshiftNX(k)=$/)唏 47Nx(j) = (1/N)乞X吨丿k =

37、 l这里Y = fft(X)使用快速傅立叶算法计算，返回序列X的DFT(离散傅立叶变换) 如果X是矩阵，则按列求傅立叶变换。Y = fft(X, n)返回X的n点DFT。如果X的长度比n小，则自动在X末尾补零；如果X长度比n大, 操作。如果X是矩阵，则对其列长度进行自动调整。Y = fft(X,dim)Y = fft(X, n, dim)对X的第dim维进行DFTy= ifft(X)*y= ifft(X, n)y= ifft(X,dim)y= ifft(X,n ,dim)举例：1、一个频率为50Hz的正弦信号求其FFT。不截断和截断的情况f=50; fs=10*f; t=0:1/fs:1; x

38、=si n( 2*pi*f*t); y=fft(x); y2=fft(x,256); figure;subplot(3,1,1) plot(t,x); subplot(3,1,2) plot(abs(y); subplot(3,1,3)plot(abs(y2);则进行截短-130002001000100200300400500600100501001502002503000.20.30.40.50.60.70.80.90.150若画图改成用stem函数0-100.10.20.30.40.50.60.70.80.91frrrr1.显 4J. H J.* J *工上*-k- J 4 J丄J.*血乩

39、出扯鼻丿rtf30405060708090100110Lk 1ipi朴豐呐.J;.二.Mi) 111lib.806040200501001502002502、分析截短对fft结果的影响-1300对连续的单一频率周期信号按采样频率1,-；.：采样，截取长度N分别选N =20和N=16，观察其DFT结果的幅度谱。解此时离散序列二一二匕二二二二，即k=8。用MATLAB计算并作图，函数fft用于计算离散傅里叶变换DFT程序如下：k=8;n1= 0:1:19;xa1=si n(2*pi* n1/k);subplot(2,2,1)plot( n1,xa1)xlabel(t/T);ylabel(x( n)

40、;xk1=fft(xa1);xk1=abs(xk1);subplot(2,2,2)stem( n1,xk1)xlabel(k);ylabel(X(k);n 2=0:1:15;xa2=s in( 2*pi* n2/k);subplot(2,2,3)plot( n2,xa2)xlabel(t/T);ylabel(x( n);xk2=fft(xa2);xk2=abs(xk2);subplot(2,2,4)stem( n2,xk2) xlabel(k);ylabel(X(k);图2 1不同载駅长度的正弦信号蜃其DFT结果计算结果示于图2.1,(a)和(b)分别是N=20时的截取信号和DFT结果，由 于截取了两个半周期，频谱出现泄漏；(c)和(d)分别是N=16时的截取信号和DFT结果，由于截取了两个整周期，得到单一谱线的频谱。上述频谱的误差主要 是由于时域中对信号的非整周期截断产生的频谱泄漏。3、分析FFT所求结果的对称性例：序列x1=1 2 3 4 5 6 7 ,x2=1 2 3 4 5 6 7 8,分别求其fft，看结果如何x1=1 2 3 4 5 6 7 ;x2=1 2 3 4 5 6 7 8; y1=fft(x1); y2=fft(x2);z仁fftshift(yl);