连续系统的频域分析.ppt

《连续系统的频域分析.ppt》由会员分享，可在线阅读，更多相关《连续系统的频域分析.ppt（99页珍藏版）》请在装配图网上搜索。

1、第四章 连续系统的频域分析,4.1 信号分解为正交函数,一、矢量正交与正交分解,时域分析，以冲激函数为基本信号，任意输入信号可分解为一系列冲激函数；而yf(t) = h(t)*f(t)。 本章将以正弦信号和虚指数信号ejt为基本信号，任意输入信号可分解为一系列不同频率的正弦信号或虚指数信号之和。 这里用于系统分析的独立变量是频率。故称为频域分析。,矢量Vx = ( vx1, vx2, vx3)与Vy = ( vy1, vy2, vy3)正交的定义： 其内积为0。即,由两两正交的矢量组成的矢量集合-称为正交矢量集,如三维空间中，以矢量 vx=（2，0，0）、vy=（0，2，0）、vz=（0，0，

2、2) 所组成的集合就是一个正交矢量集。,例如对于一个三维空间的矢量A =(2，5，8)，可以用一个三维正交矢量集 vx，vy，vz分量的线性组合表示。即 A= vx+ 2.5 vy+ 4 vz 矢量空间正交分解的概念可推广到信号空间，在信号空间找到若干个相互正交的信号作为基本信号，使得信号空间中任意信号均可表示成它们的线性组合。,二、信号正交与正交函数集,1. 定义：,定义在(t1，t2)区间的两个函数 1(t)和 2(t),若满足,(两函数的内积为0),则称 1(t)和 2(t) 在区间(t1，t2)内正交。,2. 正交函数集：,若n个函数 1(t)， 2(t)， n(t)构成一个函数集，当

3、这些函数在区间(t1，t2)内满足,则称此函数集为在区间(t1，t2)的正交函数集。,3. 完备正交函数集：,如果在正交函数集1(t)， 2(t)， n(t)之外，不存在函数(t)(0）满足,则称此函数集为完备正交函数集。,例如：三角函数集1，cos(nt)，sin(nt)，n=1,2, 和虚指数函数集ejnt，n=0，1，2，是两组典型的在区间(t0，t0+T)(T=2/)上的完备正交函数集。,( i =1，2，n),例如：,三、信号的正交分解,设有n个函数 1(t)， 2(t)， n(t)在区间(t1，t2)构成一个正交函数空间。将任一函数f(t)用这n个正交函数的线性组合来近似，可表示为

4、 f(t)C11+ C22+ Cnn,如何选择各系数Cj使f(t)与近似函数之间误差在区间(t1，t2)内为最小。,通常使误差的方均值(称为均方误差)最小，非平均误差最小。均方误差为,为使上式最小,展开上式中的被积函数，并求导。上式中只有两项不为0，写为,即,所以系数,代入，得最小均方误差（推导过程见教材）,在用正交函数去近似f(t)时，所取得项数越多，即n越大，则均方误差越小。当n时（为完备正交函数集），均方误差为零。此时有,上式称为(Parseval)巴塞瓦尔公式，表明：在区间(t1,t2) f(t)所含能量恒等于f(t)在完备正交函数集中分解的各正交分量能量的总和。,当n 时，均方误差为

5、0， 函数f(t)可分解为无穷多项正交函数之和,4.2 傅里叶级数,一、傅里叶级数的三角形式,设周期信号f(t)，其周期为T，角频率=2/T，当满足狄里赫利(Dirichlet)条件时，它可分解为如下三角级数 称为f(t)的傅里叶级数,系数an , bn称为傅里叶系数,可见， an 是n的偶函数， bn是n的奇函数。,式中，A0 = a0,上式表明，周期信号可分解为直流和许多余弦分量。 其中， A0/2为直流分量； A1cos(t+1)称为基波或一次谐波，它的角频率与原周期信号相同； A2cos(2t+2)称为二次谐波，它的频率是基波的2倍； 一般而言，Ancos(nt+n)称为n次谐波。,可

6、见An是n的偶函数， n是n的奇函数。 an = Ancosn， bn = Ansin n，n=1,2,将上式同频率项合并，可写为,例41 试将下图所示的方波信号f(t)展开为傅里叶级数。,解：傅里叶系数, an,方波的组成,二、波形的对称性与谐波特性,1 .f(t)为偶函数对称纵坐标,bn =0，展开为余弦级数。,2 .f(t)为奇函数对称于原点,an =0，展开为正弦级数。,实际上，任意函数f(t)都可分解为奇函数和偶函数两部分，即 f(t) = fod(t) + fev(t) 由于f(-t) = fod(-t) + fev(-t) = -fod(t) + fev(t) 所以,3 .f(t

7、)为奇谐函数f(t) = f(tT/2),此时 其傅里叶级数中只含奇次谐波分量，而不含偶次谐波分量即 a0=a2=b2=b4=0,三、傅里叶级数的指数形式,三角形式的傅里叶级数，含义比较明确，但运算常感不便，因而经常采用指数形式的傅里叶级数。可从三角形式推出：利用 cosx=(ejx + ejx)/2,上式中第三项的n用n代换，A n=An， n= n， 则上式写为,令A0=A0ej0ej0t ，0=0,所以,令复数,称其为复傅里叶系数，简称傅里叶系数。,n = 0, 1, 2，,表明：任意周期信号f(t)可分解为许多不同频率的虚指数信号之和。 F0 = A0/2为直流分量。,直流和n次谐波分

8、量在1电阻上消耗的平均功率之和。 n0时， |Fn| = An/2。,周期信号一般是功率信号，其平均功率为,4.3 周期信号的频谱及特点,一、信号频谱的概念,从广义上说，信号的某种特征量随信号频率变化的关系，称为信号的频谱，所画出的图形称为信号的频谱图。 周期信号的频谱是指周期信号中各次谐波幅值、相位随频率的变化关系，即 将An和n的关系分别画在以为横轴的平面上得到的两个图，分别称为振幅频谱图和相位频谱图。因为n0，所以称这种频谱为单边谱。 也可画|Fn|和n的关系，称为双边谱。若Fn为实数，也可直接画Fn 。,例：周期信号 f(t) = 试求该周期信号的基波周期T，基波角频率，画出它的单边频

9、谱图。,解 首先应用三角公式改写f(t)的表达式，即,显然1是该信号的直流分量。,的周期T1 = 8,的周期T2 = 6,所以f(t)的周期T = 24，基波角频率=2/T = /12,是f(t)的/4/12 =3次谐波分量；,是f(t)的/3/12 =4次谐波分量；,画出f(t)的单边振幅频谱图、相位频谱图如图,二、周期信号频谱的特点,举例：有一幅度为1，脉冲宽度为的周期矩形脉冲，其周期为T，如图所示。求频谱。,令Sa(x)=sin(x)/x (取样函数）, n = 0 ,1，2，,Fn为实数，可直接画成一个频谱图。设T = 4画图。,零点为,特点： (1)周期信号的频谱具有谐波(离散)性。

10、谱线位置是基频的整数倍；(2)一般具有收敛性。总趋势减小。,谱线的结构与波形参数的关系(图见教材P131P132)：,(a) T一定，变小，此时（谱线间隔）不变。两零点之间的谱线数目：1/=（2/）/(2/T)=T/ 增多。 (b) 一定，T增大，间隔减小，频谱变密。幅度减小。 如果周期T无限增长（这时就成为非周期信号），那么，谱线间隔将趋近于零，周期信号的离散频谱就过渡到非周期信号的连续频谱。各频率分量的幅度也趋近于无穷小。,直流和n次谐波分量在1电阻上消耗的平均功率之和。 n0时， |Fn| = An/2。,周期信号一般是功率信号，其在1电阻上消耗的 平均功率，称为归一化平均功率。,三、周

11、期信号的功率,傅里叶级数展开式代入,例： 试计算图所示信号在频谱第一零点及以内各分量的功率所占总功率的百分比,解：由图可求得信号f(t)的功率,将f(t)展开为指数型傅里叶级数,频谱图如右图，频谱的第一个零点在n=5，=10,在频谱第一个零点内的各分量的功率和为,将|Fn|带入，得,即频谱第一个零点以内各分量的功率占总功率的90.3%,4.4 非周期信号的频谱傅里叶变换,一、傅里叶变换,非周期信号f(t)可看成是周期T时的周期信号。 前已指出当周期T趋近于无穷大时，谱线间隔趋近于无穷小，从而信号的频谱变为连续频谱。各频率分量的幅度也趋近于无穷小，不过，这些无穷小量之间仍有差别。 为了描述非周期

12、信号的频谱特性，引入频谱密度的概念。令,(单位频率上的频谱）,称F(j)为频谱密度函数。,考虑到：T，无穷小，记为d； n （由离散量变为连续量），而,同时， ,于是，,傅里叶变换式“-”,傅里叶反变换式,F(j)称为f(t)的傅里叶变换或频谱密度函数，简称频谱。 f(t)称为F(j)的傅里叶反变换或原函数。,根据傅里叶级数,也可简记为,F(j) = F f(t) f(t) = F 1F(j) 或 f(t) F(j),F(j)一般是复函数，写为 F(j) = | F(j)|e j () = R() + jX(),说明 (1)前面推导并未遵循严格的数学步骤。可证明，函数f(t)的傅里叶变换存在的

13、充分条件: 在无限区间里面f(T)绝对可积,(2)用下列关系还可方便计算一些积分,二、常用函数的傅里叶变换,单边指数函数f(t) = et(t)， 0实数,振幅频谱 相位频谱,2. 双边指数函数f(t) = et , 0,3. 门函数(矩形脉冲),4. 冲激函数(t)、(t),均匀谱或者白色频谱,5. 常数1,有一些函数不满足绝对可积这一充分条件，如1，(t) 等，但傅里叶变换却存在。直接用定义式不好求解。 可构造一函数序列fn(t)逼近f (t) ，即,而fn(t)满足绝对可积条件，并且fn(t)的傅里叶变换所形成的序列Fn(j)是极限收敛的。则可定义f(t)的傅里叶变换F (j)为,这样定

14、义的傅里叶变换也称为广义傅里叶变换。,构造 f(t)=e-t ， 0,所以,又,因此， 12(),另一种求法： (t)1代入反变换定义式，有,将t，t-,再根据傅里叶变换定义式，得,6. 符号函数,7. 阶跃函数(t),归纳记忆：,1. F 变换对,2. 常用函数 F 变换对：,(t),(t),e -t (t),g(t),sgn (t),e |t|,1,1,2(),4.5 傅里叶变换的性质,一、线性(Linear Property),If f1(t) F1(j)， f2(t) F2(j) then,Proof: F a f1(t) + b f2(t),= a F1(j) + b F2(j) ,

15、a f1(t) + b f2(t) a F1(j) + b F2(j) ,For example F(j) = ?,Ans: f (t) = f1(t) g2(t),f1(t) = 1 2(),g2(t) 2Sa(), F(j) = 2() - 2Sa(),-,二 奇偶性,其中，f(t)是t的实函数,得到：1、,2、3、,三、对称性质(Symmetrical Property),If f (t) F(j) then,Proof:,（1）,in (1) t ，t then,（2）,in (2) - then, F(j t) 2f () end,F( jt ) 2f (),For example,

16、 F(j) = ?,Ans:,if =1,例： 求取样函数Sa（t）=sin t /t的频谱函数,解： 已知，宽度为，幅度为1的门函数，g(t)的频谱函数为,取/2=1, =2，幅度为1/2,即,根据对称性， g2(t)是偶函数,四、时移性质(Timeshifting Property),If f (t) F(j) then,式中 “t0” 是一个实常数,Proof: F f (t t0 ) ,例: 图f(t) F(j) = ?,Ans: f1(t) = g6(t - 5) , f2(t) = g2(t - 5),g6(t - 5) ,g2(t - 5) , F(j) =,+,五、尺度变换性质

17、(Scaling Transform Property),If f (t) F(j) then,where “a” is a nonzero real constant.,Proof:,F f (a t ) =,For a 0 ,F f (a t ) ,for a 0 ,F f (a t ) ,That is ,f (a t ) ,Also,letting a = -1,f (- t ) F( -j),例 1,已知 f (t)F( j), 求 f (at b) ?,Ans: f (t b),e -jb F( j),f (at b) ,or,f (at) ,f (at b) =,例 ：2,f(t

18、) = F(j) = ?,解,利用已知单边函数的福利叶变换,利用尺度变化性质， a = -1,得,六、频移性质(Frequency Shifting Property),如 f (t) F(j) 那么,Proof:,式中 “0” 为实常数,F e j0t f(t),= F j(-0) end,For example 1,f(t) = ej3t F(j) = ?,Ans: 1 2() ej3t 1 2(-3),For example 2,f(t) = cos0t F(j) = ?,Ans:,F(j) = (+0)+ (-0),For example 3,Given that f(t) F(j),

19、那么信号 f(t) cos0t ?,七、卷积性质(Convolution Property),卷积积分 在时域里面：,If f1(t) F1(j)， f2(t) F2(j) Then f1(t)*f2(t) F1(j)F2(j),卷积积分在频域里面：,If f1(t) F1(j)， f2(t) F2(j),Then f1(t) f2(t) F1(j)*F2(j),Proof:,F f1(t)*f2(t) =,利用时移性质：,得到,F f1(t)*f2(t) =,= F1(j)F2(j),For example,Ans:,利用对称性：,七、时域的微分和积分 (Differentiation an

20、d Integration in time domain),If f (t) F(j) then,Proof:,f(n)(t) = (n)(t)*f(t) (j )n F(j) f(-1)(t)= (t)*f(t) ,f(t)= 1/t2 ?,For example 1,Ans:,例2,Determine f (t) F (j),Ans:,f ”(t) = (t+2) 2 (t) + (t 2),F2(j)= F f ”(t) = e j2 2 + e j2= 2cos(2) 2,F (j) =,Notice:,d(t)/dt = (t) 1,(t) 1/(j),If f (t) F(j) t

21、hen,(jt)n f (t) F(n)(j),where,例 1,f (t) = t(t) F (j)=?,解:,九、频域微分和积分,Notice: t(t) =(t) * (t) ,Its wrong. Because ()() and (1/j)() is not defined.,例 2,求右式的值,Ans:,十、能量谱和功率谱,|F(j)|2 单位频率上的频谱 (能量密度谱）Js,上式称为帕斯瓦尔方程或能量等式。,定义能量频谱E(w),为单位频率的信号能量,能量谱的定义,例题：,求信号的能量谱,解：由,得到,类似能量谱,定义能量频谱E(w),为单位频率的信号能量,功率谱的定义,P(w

22、)是w的偶函数，由频谱函数的模量决定，单位频谱的能量，单位 Ws,4.6 周期信号的傅里叶变换,一、正、余弦的傅里叶变换,12() 由频移特性得 e j 0 t 2(0 ) e j 0 t 2(+0 ) cos(0t)=(e j 0 t + e j 0 t) (0 ) +(+0 ) sin(0t)= (e j 0 t - e j 0 t)/(2j) j(+0 ) ( 0 ),二、一般周期信号的傅里叶变换,例1：周期为T的单位冲激周期函数T(t)=,解：,(1),例2：周期信号如图，求其傅里叶变换。,解：周期信号f(t)也可看作一时限非周期信号f0(t)的周期拓展。即,f(t) = T(t)*

23、f0(t),F(j) = () F0(j),F(j) =,本题 f0(t) = g2(t),(2),(2)式与上页(1)式比较，得,这也给出求周期信号傅里叶级数的另一种方法。,4.7 LTI系统的频域分析,傅里叶分析是将任意信号分解为无穷多项不同频率的虚指数函数之和。,对周期信号：,对非周期信号：,其基本信号为 ej t,一、基本信号ej t作用于LTI系统的响应,说明：频域分析中，信号的定义域为(，)，而t= 总可认为系统的状态为0，因此本章的响应指零状态响应，常写为y(t)。,设LTI系统的冲激响应为h(t)，当激励是角频率的基本信号ej t时，其响应,而上式积分 正好是h(t)的傅里叶变

24、换，记为H(j )，常称为系统的频率响应函数。,y(t) = H(j ) ej t,H(j )反映了响应y(t)的幅度和相位。,y(t) = h(t)* ej t,二、一般信号f(t)作用于LTI系统的响应,ej t,H(j ) ej t,F(j ) ej t d ,F(j )H(j ) ej t d ,可加性,f(t),y(t) =F 1F(j )H(j ) ,Y(j ) = F(j )H(j ),齐次性,频率响应H(j)可定义为系统零状态响应的傅里叶变换Y(j)与激励f(t)的傅里叶变换F(j)之比，即,H(j)称为幅频特性（或幅频响应）；()称为相频特性（或相频响应）。H(j)是的偶函数

25、，()是的奇函数。,频域分析法步骤：,傅里叶变换法,对周期信号还可用傅里叶级数法。,周期信号,若,则可推导出,例：某LTI系统的H(j)和()如图， 若f(t)= 2 + 4cos(5t) + 4cos(10t)，求系统的响应。,解法一：用傅里叶变换法,F(j) = 4() + 4(5) + (+5) + 4(10) + (+10),Y(j) = F(j)H(j) = 4() H(0) + 4(5) H(j5) + (+5) H(-j5) + 4(10) H(j10) + (+10) H(-j10) ,H(j)=H(j)ej(),= 4() + 4-j0.5(5) + j0.5(+ 5) ,y

26、(t) = F-1Y(j) = 2 + 2sin(5t),解法二：用三角傅里叶级数,f(t)的基波角频率=5rad/s,f(t)= 2 + 4cos(t) + 4cos(2t) Fn=4 n=0,1,2,由图得到， H(0) =1， H(j) = 0.5e-j0.5， H(j2) = 0,y(t) = 2 + 40.5cos(t 0.5) = 2 + 2sin(5t),三、频率响应H(j)的求法,1. H(j) = F h(t),2. H(j) = Y(j)/F(j) 由微分方程求，对微分方程两边取傅里叶变换。 由电路直接求出。,例1：某系统的微分方程为 y(t) + 2y(t) = f(t)

27、 求f(t) = e-t(t)时的响应y(t)。,解：微分方程两边取傅里叶变换,jY(j) + 2Y(j) = F(j),f(t) = e-t(t),Y(j) = H(j)F(j),y(t) = (e-t e-2t )(t),例2：如图电路，R=1，C=1F，以uC(t)为输出，求其h(t)。,解：画电路频域模型,h(t)= e-t (t),四、无失真传输与滤波,系统对于信号的作用大体可分为两类：一类是信号的传输，一类是滤波。传输要求信号尽量不失真，而滤波则滤去或削弱不需要有的成分，必然伴随着失真。,1、无失真传输,（1）定义：信号无失真传输是指系统的输出信号与输入信号相比，只有幅度的大小和出

28、现时间的先后不同，而没有波形上的变化。即 输入信号为f(t)，经过无失真传输后，输出信号应为 y(t) = K f(ttd) 其频谱关系为 Y(j)=Ke jtdF(j),系统要实现无失真传输，对系统h(t)，H(j)的要求是： (a)对h(t)的要求： h(t)=K(t td) (b)对H(j)的要求： H(j)=Y(j)/F(j)=Ke-jtd 即 H(j)=K ,()= td,上述是信号无失真传输的理想条件。当传输有限带宽的信号是，只要在信号占有频带范围内，系统的幅频、相频特性满足以上条件即可。,(2)无失真传输条件：,例：系统的幅频特性|H(j)|和相频特性如图(a)(b)所示，则下列

29、信号通过该系统时，不产生失真的是,(A) f(t) = cos(t) + cos(8t) (B) f(t) = sin(2t) + sin(4t) (C) f(t) = sin(2t) sin(4t) (D) f(t) = cos2(4t),2、理想低通滤波器,具有如图所示幅频、相频特性的系统称为理想低通滤波器。c称为截止角频率。 理想低通滤波器的频率响应可写为：,(1)冲激响应,h(t)= -1g 2 c()e-jtd =,可见，它实际上是不可实现的非因果系统。,(2)阶跃响应,g(t)=h(t)*(t)=,经推导，可得,称为正弦积分,特点：有明显失真，只要c，则必有振荡，其过冲比稳态值高约

30、9%。这一由频率截断效应引起的振荡现象称为吉布斯现象。,gmax=0.5+Si()/=1.0895,3、物理可实现系统的条件,就时域特性而言，一个物理可实现的系统，其冲激响应在t0时必须为0，即 h(t)=0 ,t0 即 响应不应在激励作用之前出现。 就频域特性来说，佩利（Paley)和维纳（Wiener)证明了物理可实现的幅频特性必须满足,并且,称为佩利-维纳准则。（必要条件） 从该准则可看出，对于物理可实现系统，其幅频特性可在某些孤立频率点上为0，但不能在某个有限频带内为0。,4.8 取样定理,取样定理论述了在一定条件下，一个连续信号完全可以用离散样本值表示。这些样本值包含了该连续信号的全

31、部信息，利用这些样本值可以恢复原信号。可以说，取样定理在连续信号与离散信号之间架起了一座桥梁。为其互为转换提供了理论依据。,一、信号的取样,所谓“取样”就是利用取样脉冲序列s(t)从连续信号f(t)中“抽取”一系列离散样本值的过程。 这样得到的离散信号称为取样信号。,如图一连续信号f(t),用取样脉冲序列s(t)(开关函数）进行取样，取样间隔为TS，fS =1/TS称为取样频率。,得取样信号 fS(t) = f(t)s(t),取样信号fS(t)的频谱函数为 FS(j)=(1/2)F(j)*S(j),冲激取样,若s(t)是周期为Ts的冲激函数序列Ts(t),则称为冲激取样。,如果f(t) 是带限

32、信号 即f(t)的频谱只在区间（- m，m)为有限值，而其余区间为0 。,设f(t)F(j)，取样信号fS(t)的频谱函数,FS(j)= (1/2)F(j)* S s(),S =2/TS,s(t)=Ts(t) S s(),=,*,=,上面在画取样信号fS(t)的频谱时，设定S 2m ,这时其频谱不发生混叠，因此能设法(如利用低通滤波器)，从FS(j)中取出F(j)，即从fS(t)中恢复原信号f(t)。否则将发生混叠，而无法恢复原信号。,二、时域取样定理,当S 2m 时，将取样信号通过下面的低通滤波器,其截止角频率C取m C S -m 。即可恢复原信号。,由于 fs(t)= f(t)s(t) =

33、 f(t),H(j) h(t) =,为方便，选C = 0.5S ,则TsC /=1,所以,根据f(t)=fS(t)*h(t) ，有,只要已知各取样值f(nTs),就出唯一地确定出原信号f(t)。,时域取样定理： 一个频谱在区间（-m，m)以外为0的带限信号f(t),可唯一地由其在均匀间隔Ts Ts1/(2fm) 上的样值点f(nTs)确定。,注意：为恢复原信号，必须满足两个条件：（1）f(t)必须是带限信号；（2）取样频率不能太低，必须fs2fm，或者说，取样间隔不能太大，必须Ts1/(2fm);否则将发生混叠。,通常把最低允许的取样频率fs=2fm称为奈奎斯特（Nyquist)频率，把最大允

34、许的取样间隔Ts=1/(2fm)称为奈奎斯特间隔。,频域取样定理： 根据时域与频域的对偶性，可推出频域取样定理。P191 一个在时域区间（-tm,tm)以外为0的时限信号f(t)的频谱函数F(j)，可唯一地由其在均匀频率间隔fsfs1/(2tm)上的样值点F(jns)确定。,4.9 序列的傅里叶变换,一、周期序列的离散傅里叶变换级数 类似连续信号可以分解成为一系列虚指数 之和 周期序列也可以分解为一系列虚指数 之和 即可写成,注：周期序列的傅里叶级数展开式仅为有限项（N项）,求解系数，得到周期序列的离散傅里叶级数（DFS）,请结合连续信号傅里叶级数记忆,令 ，得到求离散傅里叶系数的表达式,二、

35、非周期序列的离散时间傅里叶变换（DTFT）,当N 时，周期信号变为非周期信号,当N 时， 并趋向于无穷小 n ， 的周期变为2,得到,例：求单位取样序列的DTFT,求方波序列的DTFT,4.10 离散傅里叶变换及其性质,一、离散傅里叶变换（DFT）,利用周期序列离散傅里叶级数的概念推导出有限长序列的傅里叶变换,将有限长序列看作该序列在定义域上的无限延展序列的“主值区间”,周期序列的离散傅里叶级数变换对的两个求和公式都只限于“主值区间” 设有限长序列的长度为N，f(k)离散傅里叶变换和其逆变换的定义分别为,例： 对方波求离散傅里叶变换（设N=10）,解：P195,二、离散傅里叶变换的性质,1、线性,2、对称性,3、时移特性,4、频移特性,5、时域循环卷积定理,6、频域循环卷积定理,