第四章 连续系统的频域分析 时域分析:f(t) yf(t)=h(t)*f(t) ( 分解 ( 基本信号((t)→ LTI →h(t) 频域分析: f(t) yej(t =h(t)* H(j()Fej(t ( 分解 ( 基本信号 sin(t → LTI → H(j()ej(t ej(t H(j():系统的频域响应函数,是信号角频率(的函数,与t无关. 主要内容: 一、信号的分解为正交函数。 二、周期信号的频域分析(付里叶级数(求和),频谱的特点。 信号 三、非周期信号的频域分析(付里叶变换(积分),性质。 分析 四、LTI系统的频域分析:频域响应H(j();y(j()= H(j()?F(j(). (系统分析) 五、抽样定理:连续信号(离散信号. §4.1 信号分解为正交函数 一、正交: 两个函数满足 φ1(t) φ2(t)dt=0, 称φi(t),φj(t)在区间( t1 ,t2)正交。 二、正交函数集:几个函数 φi(t) φi(t)dt= 0 当i≠j; Ki 当i=j. 三、完备正交函数集:在{φ1(t)… φn(t)}之外, 不存在((t)满足 ( (t) φi(t)dt= 0 (i=1,2,…n). 例、三角函数集:{1,cos(t,cos2(t,… ,cosm(t,…,sin(t, sin2(t,…sin(n(t),…}区间:(t0,t0+T),t=2π/(为周期. 满足: cosm(tcosn(tdt= 0 m≠n T/2 m=n≠0 T m=n=0 sin(m(t)sin(n(t)dt= 0 m≠n T/2 m=n≠0 sin(m(t)cos(n(t)dt= 0. 所有的m和n. 结论:三角函数集是完备正交集。 推导: cosm(tcosn(tdt =(1/2) [cos(m+n) (t+cos(m-n) (t]dt =(1/2)sin(m+n)(t +(1/2)sin(m-n)(t =(1/2)[sin(m+n) ((t0+T)-sin(m+n)(t0] +(1/2)[sin(m-n) ((t0+T)-sin(m-n)(t0] =0 当m≠n时. m=n≠0,原式=(1/2) [ cos(m+n)(t+1]dt=(1/2)?t =T/2 m=n=0 , 原式=(1/2) [1+1]dt=T. 4、复函数的正交函数集: 几个复函数集{φi(t)}, φi(t) φi* (t)dt= 0 i≠j ki i=j 例:复函数集{ ejnΩt}(n=0,±1,±2…) 区间(t0,t0+T),T=2π/(为周期。 满足 ejm (t(ejn(t)*dt= ej(m-n)(t dt =[1/(j(m-n)Ω)] ej(m-n)(t dt =0 m≠n = 1dt=T m=n. 结论:{ ejnΩt}是完备正交集。(n=0,±1,±2…) 二、信号分解为正交函数集。 1、分解: 二维 A=c1vx + c2yy { vx ,v}y二维正交矢量集 三维 A= c1vx +c2vy +c3vz { vx ,vy,vz }三维正交矢量集 n维:{φ1(t)… φn(t)}在( t1 ,t2)构成正交函数集。 f(t)≈c1φ1 (t)+ c2φ2(t)+…cnφn(t)+(t)= cjφj(t) 任意一个函数可以用这几个正交函数的线性组合来近似。 2、系数cj的选择。 方均误差定义:=[1/(t2-t1)] [f(t)- cjφj(t)]2dt 使 最小,对第i个系数ci来说,应使/ci =0. ∴cj= [ f(t) φj(t)dt]/ ( [φj(t)]2dt) =(1/Kj) f(t) φj(t)dt 最佳近似条件下的方均误差: =[1/(t2-t1)]( [f(t)]2 dt - cj2Kj). ∵ ≥0,n(, (; ∴n(∞, (0. 则 [f(t)]2 dt= cj2Kj(称帕斯瓦尔方程。 f(t)= cjφj(t). 即函数f(t)在区间( t1 ,t2)可分解为无穷多项正交函数之和。 §4.2付里叶级数 一、付里叶级数:(三角形式) f(t)=(a0/2)·1+a1cos(t+a2cos2(t+…+b 1sin(t+b 2sin2(t+… = a0/2+ ancos(n(t)+ b nsin(n(t). 积分区间:t0 t0+T, 0T, -T/2T/2 Ki= (cos(n(t))2 dt=T/2. an=(2/T)  f(t)cos(n(t)dt bn=(2/T)f(t)sin(n(t)dt 形制:a-n=an是偶函数 b-n=-bn时奇函数 (其中n=0,1,2…). 2、三角形式二:同频率项合并。 f(t)=a0/2+A1cos((t+φ1)+A2cos(2(t+φ2) +… = a0/2+ Ancos(n(t+φn). A0=a0 an= bn =-arctg(bn / an). 由性质可知:a0= A0 an=Ancosφn bn= bn sinφn 3、物理意义;同周期信号可分解为各次谐波之和。 f(t)= a0/2+A1cos((t+φ1)+A2cos(2(t+φ2) +…+Ancos((t+φn)+… 例4.2-1 f(t)为方波,分解为付里叶级数。 周期:T 频率:1/T 角频率:(=2(/T. 区间:(-T/2,T/2) (1)f(t)= a0/2+ ancos(n(t)+ bnsin(n(t) an=(2/T) f(t)cos(nΩt)dt =0 bn=(2/T) f(t)sin(nΩt)dt= 0 n=2,4,6… . 4/(n() n=1,3,5… ∴f(t)=(4/()[sin(t+(1/3)sin(3(t)+…+ (1/n)sin(n(t)+ …] 结论:方波只含有1,3,5等奇次谐波分量,无直流分量。 (2)方均误差(有限项逼近) =[1/(t2-t1)][ f2(t)dt- c2jKj] =(1/T)[ 1dt-(T/2) (bj)2]=1-(1/2) (bj)2 只取基波:=1-(1/2)(4/()2=0.189. 取基三次谐波:=1-(1/2)[(4/()2+(4/3()2=0.0994. 基“+”3,”+”5次: =1-(1/2)[(4/()2+(4/3()2+(4/5()2]=0.0669 (3)方波分解的特点 1、它包含的基波分量越多,越接近方波,其均方误差越小。 2、当合成波所含基波次数n(∞,在间断点仍有约9(偏差,在间断点出尖峰下的面积非常小以致趋近于零。 二、奇偶函数的付里叶系数的特点: 1、为偶函数:f(-t)=f(t),关于纵坐标对称。 an=(2/T) f(t)cos(n(t)dt =(2/T)f(t) cos(n(t)dt +(2/T)f(t) cos(n(t)dt ∴an=(4/T)f(t) cos(n(t)dt bn=(2/T) f(t)sin(n(t)dt+(2/T)f(t)sin(n(t)dt ∴ bn= 0. 当f(t)为偶函数时 an=(4/T)f(t) cos(n(t)dt An= |an| bn/ an=0 bn= 0 (n= m( arctgbn/an角度为0,( 2、f(t)为奇函数。F(-t)=-f(t),波形关于原点对称。 当f(t)为奇函数时: an=0 An= |bn| bn=(4/T)f(t)sin(nΩt)dt (n= (2m+1)(/2. bn/an(∞. ∴ 奇函数只有正弦项。 任意函数 f(t)=fod(t)+fev(t) ( fod(t)=(f(t)-f(-t))/2. f(-t)= fod(-t)+fev(-t)= -fod(t)+fev(t) fev(t)=(f(t)+f(-t))/2. 3、f(t)为奇谐函数。(半波对称函数) f(t)=- f(t(T/2),移动T/2后,关于横轴对称。 付里叶级数只含奇次谐波,不含偶次谐波。 a0= a2= a4= a6=( b0= b2= b4=(=0 例4.2-2 把锯齿波信号展为付里叶级数。 解:     方法1:f(t)=t/T既不是偶函数也不是奇函数, 直接在[0,T]区间上求an ,bn . 方法二:把分为奇偶两部分。 fev(t)=(1/2)[f(t)-f(-t)]=(1/2)[t/T+(-t+T)/T]=1/2. fod(t)=(1/2)[f(t)+f(-t)]=(1/2)[t/T-(-t+T)/T]=t/T-1/2=(t-T/2)/T. 奇函数部分分解为:an bn =(4/T)[t/T-1/2]sin(n(t)dt =(4/T2)[sin(n(t)-n(cos(n(t)]/(n()2 +(2/T)[cos(n(t)]/(n()]=-1/n(. n=1,2,3… ∴f(t)= fev(t)+fod(t)=1/2+ bn sin(n(t) =1/2-(1/()[sin(t+(1/2)sin(2(t)+(1/3) sin(3(t)+…]. 锯齿波含直流分量和各次谐波分量。 三、周期信号分解为指数形付里叶级数。 1、定义式:(由三角形式推导) f(t)=A0/2+ Ancos(n(t+φn) = A0/2+ (An/2)[ej(nΩt+φn)+e -j(nΩt+φn)]。 ∴ f(t)= FnejnΩt 2、确定付里叶系数Fn Fn=(1/2) Anejφn+(1/2)[Ancosφn)+jAnsinφn]=(1/2)(an-jbn) =(1/2)(2/T) f(t)cos(nΩt)dt -j(1/2)(2/T) f(t)sin(nΩt)dt =(1/T)f(t)[cos(nΩt)-jsin(nΩt)]dt ∴ Fn=(1/T) f(t)e-jnΩtdt. n=0,±1,±2… 3、物理意义:周期信号可分解为许多不同频率(n()的虚指数信号(ejnΩt)之和。 每个分量的大小用Fn来表示,分为幅度和相位。 各三角函数型和指数型付里叶级数及其系数,以及各系数间的关系见表4-1。 §4.3 周期信号的频谱 一、频谱的概念: 频谱分为 ( 幅度频谱:以频率ω(或角频率()为横坐标,An/|Fn|为纵坐标。 ( 相位频谱:以频率ω(或角频率()为横坐标,φn 为纵坐标。 f(t)=A0/2+ Ancos(nΩt+φn) A0为直流分量幅度;An为n次谐波的振幅;φn为n次谐波的初相角。  周期信号的频谱是离散的。 结论:正如波形是信号在时域的表示一样,频谱则是信号在频域的表示。 描述了一个信号的频谱就等于描述了这个信号。 信号分解:从已知信号绘制其频谱图。 合成:根据其频谱图反过来和成原有的信号。 波形f(t) 频谱 Fn与An比较: An:每条谱线代表一个完整的谐波分量的幅度,物理意义明确。 Fn:从数学上将cosn(t分成ejnΩt和e-jnΩt,有负频率,没有物理意义。变化趋势一致都可进行信号的频谱分析。|Fn|=(1/2)An. 3、周期信号频谱的特点:离散性;谐波性(是基波频率的整数倍)。 二、周期矩形脉冲的频谱。 f(t)幅度为1,脉冲宽度为(;周期为T. 1、求频谱:f(t)=  nejnΩt Fn=(1/T)f(t) e-jnΩtdt=(1/T) e-jnΩtdt =((/T)[sin(n(()/(n()]=((/T)[sin(n((/2)/(n((/2)] = ((/T)Sa(n((/2) 或(=2(/T. n=((/T)[sin(n2((/2T)/(n2((/2T)] =((/T)Sa(n((/2). N=0,(1,(2( (1) ( f(t)= ((/T)Sa(n((/2) ejnΩt 是指数形式的付里叶级数展开式。 由(1)式画出矩形脉冲信号频谱图。设T=4( Fn=((/4()Sa(n((/4()=(1/4)[sin(n((/4)/(n(/4)] = sin(n()/(n() n=0,(1,(2( n=0 F0 =1/4=0.25 [∵Sa(x)=1,当x(0时] n=1 F1= sin((/4)/(=0.225.  n=2 F2= sin((/2)/2(=0.16 n=3 F3= sin(3(/4)/3(=0.075. n=4 F4= sin(()/4(=0. n=5 F5= sin(5(/4)/5(=-0.045. n=6 F6= sin(3(/2)/6(=-0.053. n=7 F7= sin(7(/4)/7(=-0.032. n=8 F8= sin(2()/8(=0. 特点:1、是离散的,仅含有(=n(的各分量。(n取整数)。 2、谱线间隔为(((=2(/T) T( 间隔小,密 T( 间隔大,疏 3、第一零点在2(/(处,与(有关 (( 主瓣宽 (( 主瓣窄。 2、脉冲宽度与频谱的关系:(( 直流分量F0=(/T( 频带宽度(F=1/(( ∴ 保持第一零点内能量不变。 脉冲宽度(() 频谱幅度(F0=(/T) 第一零点 (F=1/( (=T/4 F0=1/4 2(/(=8(/T 4/T (=T/8 F0=1/8 2(/(=16(/T 8/T (=T/16 F0=1/16 2(/(=32(/T 16/T   3、周期与频谱的关系。谱线间隔保持第一零点内能量不变 F0=(/4 (=2(/T T=4( F0=1/4 (=2(/4( T=8( F0=1/8 (=2(/8( T=16( F0=1/16 (=2(/16(   T(( ,频谱趋于一个脉冲。 三、周期信号的功率 p= (1/T)f2(t)dt=(1/T)[ A0/2+ Ancos(nΩt+φn)]2dt = (A0/2)2 + (An)2/2. P=(1/T)f2(t)dt= |F0|2+2 |Fn|2 =|F-n|2+|F0|2+ |Fn|2=|Fn|2 例4.3-1 T=1,(=0.2 解:p= (1/T)f2(t)dt=(1) 2dt=0.2 Fn = ((/T)Sa(n((/T)=0.2·Sa(0.2(() n=0,1,2,3,4,5. 确定第一个零点:2(/(= 2(/0.2=10(, (= 2(/T= 2(,n= 10(/2(=5. P10(=(0.2)2+2(0.2)2[Sa2(0.2()+Sa2(0.4()+Sa2(0.6()+Sa2(0.8()+ Sa2(()]=0.1806。 P10(/p=0.1806/0.2=90.3(. 4、=H(jn() n =(1/)Sa(0.2n()e-jarctg0.5n y(t)=  ejn(t= [Sa(0.2n()/] e-jarctg0.5n . 输出波形与时域分析相同。 §4.4 非周期信号的频谱(付里叶变换 信号分析; ( 周期信号:可展开为付里叶级数,频谱n是离散的,求和形式,满足狄里赫利条件。 非周期信号:存在付里叶变换,频谱密度F(j()是连续的,积分形式, |f(t)|dt<( 一、付里叶变换。 由周期信号非周期信号,推导出付里叶变换的定义。 1、频谱密度函数 定义:F(j()=n/(1/T)= n·T称为频谱密度函数。 n/f表示单位频率的频谱,类似于单位体积的质量,定义为物体的密度。 T((,即为非周期。 2、付里叶变换的定义: 周期信号 n·T=f(t) e-jnΩtdt (1) f(t)= n·T·ejnΩt·(1/T) = [f(t)e-jnΩtdt]·ejnΩt·(/2(. (2) 非周期信号: F(j()= n·T def f(t) e-j(t dt f(t)= [f(t)e-j(tdt]ej(td(/2(def=(1/2()f(j()ej(td( F(j()=[f(t)] f(t)=[F(j()] f(t)( F(j() F(j()与n一样,也是一复函数,讨论时可分开写为: F(j()=|F(j()|ej((()=R(()+jX(() =|F(j()|cos [((()]+j|F(j()|sin [((()]. 3、复里叶变换的物理意义 三角形式: f(t)= (1/2() F(j()ej(td(=(1/2() |F(j()|ej[(t + ((()]d( =(1/2()|F(j()|cos[(t+(((t)]d( +j(1/2()|F(j()|sin [(t+(((t)]d( =(1/()|F(j()|cos [(t+(((t)]d( 定义:非周期信号可看作是由不同频率的各余弦“分量”组成,它包含了频率从零到无限大的一切频率“分量”。 三要素:1、它包含了频率从零到无限大的一切频率“分量”,且是连续的。 2、各分量的振幅为:(1/()|F(j()|d( 它是无穷小量。 3、相位为((()。 4、付里叶变换的条件: 充分条件: f(t)在无限区间内绝对可积,即f(t) e-j(t dt<(,但这并非必要条件。 在引入((t)函数后,可将条件放宽,使许多不满足绝对可积条件的函数也能进行付里叶变换。 例4.4-3:双边指数函数(((0,衰减) e-(|t| ( (2()/((2+(2) 实函数。 例4.4-4: f2(t)= - e-(t t<0 e-(t t>0 ((>0)满足绝对可积条件。 F2(j()=- e(t·e-j(tdt+ e-(t·e-j(tdt =-1/((-j()+1/((+j()=-j2(/((2((2)。 F2(j()= R2(()+jX2(() 5、典型常用信号的付里叶变换。 ①门函数g( (t),幅度为1,宽度为(.     F(j()=  f(t)e-j(tdt=1 e-j(tdt = (e-j((((-e-j(((()/(-j()=2[sin(((/2)]/(=(·Sa((()/2 零点幅值:F(0)=( 第一零点位置在((/2=(,(=2(/(处。 信号的宽度(F=1/( , (( (F( ②单边指数函数f(t)= e-(t·((t),((>0),满足绝对可积条件。  F(j()=  f(t)e-j(tdt= e-(t·e-j(tdt =[ -1/((+j()]·e-((+j()t =[0-1]/[-( (+j()] = 1/((+j().(>0. 复函数 ( F(j()(=1/ . 偶函数 ((()=-arctg((/(). 奇函数 F(j()=  f(t)e-j(tdt=  f(t)cos((t)dt-jf(t)sin((t) dt. 特点:若f(t)是t的偶函数 (F(j()是的实函数 若f(t)是t的奇函数 ( F(j()是的虚函数。 若f(t)非奇非偶 ( F(j()为复函数,用幅度和相位才能表示。 二、奇异函数的付里叶变换。 1、((t)的频谱 ①由定义[((t)]=  ((t)e-j(tdt=1 其频谱密度在-(<(<(区间处处相等。 ②由极限概念[(1/()g((t)]= (1/()(((Sa(((/2) [((t)]= Sa(((/2)=1. 2、单位直流信号的频谱: f(t)=1, -(<t<(不满足绝对可积条件,不能用定义。 只有极限概念得到:引入(()函数后,条件放宽。 双边指数: f1(t)=e-((t( ((0 f1(t)=1.   F1(j()=2(/((2((2) ((0 2(/((2((2)= 0 ((0 2/( (=0 是一个以(为位自变量的冲激函数,强度有冲激函数定义求出。 强度:2(/((2((2)d(=2(/[1(((/()2]d((/() =2arctg((/() = 2[arctg(()-arctg(-()]=2[(/2-(-(/2)]=2(. ( [1]=2((((). 3、符号函数 sgn(t)= -1 t<0 0 t=0 1 t>0 不满足决度可积条件,不能用定义。用极限 f2(t)= -e-(t t<0 ((0 Sgn(t)= -1 t<0 e(t t>0 1 t>0  ((0  F2(j()=-2j(/((2((2) -2j(/((2((2)= -2j/(=2/j( ((0 0 (=0 F2(j()是的奇函数,在 (=0处时值为0.  4、((t)的频谱 v(t)不满足绝对可积条件,不能用定义求. v(t)=1/2+(1/2)Sgn(t).  [( (t)]=(1/2)[1]+(1/2)[Sgn(t)]=(((()(1/j(. 5、((t)的频谱 [(((t)]=(((t)e-j(tdt= j(. 同理:[((n)(t)]=(j()(n). § 4.5 付里叶变换的性质 连续时间信号有两种描述方法: ( 时域描述 f(t) ( 频域描述 F(j() 一、线性 a1f1(t)+ a2f2(t) ( a1F1(j()+ a2F2(j() 利用该性质,可将所求信号表示成已知频谱信号的线性组合,用间接方式求出频谱函数。 二、奇偶性 大前提,f(t)是实函数;f(t)与F(j()奇偶虚实关系: 推导:F(j()=f(t)e-j(tdt=f(t)cos(tdt-jf(t)sin(tdt F(j()=R(()+jX(()=(F(j()(ej((() (F(j()(= ,( (()=arctg[X(()]/[R(()]. 1、实部是偶函数R(()=R(-(),虚部是奇函数X(()=-X(-(). 模是偶函数,(F(j()(=(F(-j()(,相角是奇函数( (()=-( (-()。 2、若f(t)是偶函数f(t)=f(-t), 则X(()=0, F(j()=R(()是实函数,也是偶函数 若f(t)是奇函数f(t)=-f(-t), 则R(()=0, F(j()=jX(()是虚函数,也是奇函数 偶:∵f(t)sin(t是的奇函数,虚部积分为0,∴只有实部。 奇:∵f(t)cos(t是的奇函数,实部积分为0,∴只有虚部。 3、f(-t)( F(-j()= F((j()。 推导:[f(-t)] = f(t)e-j(tdt((-t f(()ej((d(-() = f(()e-j(-()(d(()=F(-j(). F(-j()= R(-()+jX(-()=R(()-jX(()=F((j(). 三、对称性 若f(t) ( F(j() 则F(jt) (2(f(-() 当f(t)为偶函数时: f(t) ( R((),R(t) (2(f((). 推出:F(j()=f(t)e-j(tdt f(t)=(1/2()  F(j()ej(td( 令t=-t,f(-t)= (1/2()  F(j()e-j(td( 令t=(/ (=t ,f(-()= (1/2()  F(jt)e-j(tdt ( (2() f(-()=  F(jt)e-j(tdt. 定义:时间函数F(jt) [与F(j()形式相同]的付里叶变换是(2()f((). 例:((t)( 1 1((2()((()利用对称性,可以很方便地求出一些函数的付里叶变换。 例:4.5-1 Sa(t)=sint/t. 门函数g((t) (( Sa(((/2)令(/2=1.则(=2. ∴(1/2)(g((t) (2((1/2)(Sa(()=Sa((). ( 由对称性知:Sa(t) ( 2(((1/2)(g2(()= (g2(() ( 例:f(t)=t ∵ (((t) ( j(. jt(2(( (((-()=-2(( (((() ∴ t(+j2(( ((((). 例:f(t)=1/t. 已知:sgn(t)(2/j( 则2/(jt)( 2(sgn((). ∴1/t (j(sgn(-()=-j(sgn((). 四、尺度变换(时域展缩) 若 f(t) ( F(j() 则f(at) ( (1/(a() F(j(/a)       结论:信号的等效脉冲宽度与占有的等效宽度成反比。若言压缩信号的持续时间,则不得不以占宽频带作代价。在通信中,通信速度与占用频带宽度是一对矛盾。通信系统的设计便是寻找矛盾的合理解决方案。 五、时移特性 若f(t) ( F(j() 则f(t(t0) ( e(j(t0F(j()= (F(j()( ej[(((((j(t0] 即时域中信号延时t0 (频域中所有频率“分量”相应落后一相位( t0,而幅度不变。 既有时移又有尺度变换 若f(t) ( F(j() 则f(at+b) ( (1/(a() e-j(b/a)( F(j(/a) b=0 尺度变换 a=1,时移。 例4.5-3 已知f1(t)= g2(t), F(j()= (( Sa(((/2)= 2Sa(()=zsin(/(.   解:(1) f2(t)= f1(t+1)+f1(t-1) F2(j()= ej( F1(j()- e-j( F1(j() = (ej( - e-j()2sin(/ (=4j(sin2(/ (. (2) f3(t)= f2(2t) F3(j() = (1/2) F2(j(/2) = (1/2) 4j(sin2((/2)/((/2)= j4(sin2((/2)/(. 例4.5-5 已知f(t) ( F(j() 求f(3-2t) ( ej4t的付里叶变换。 解:f(t) ( F(j() 时移:f(t+3) (ej3( F(j() 尺度变换 a=-2, f(-2t+3) ((1/(-2() ej3(/(-2) F(j(/(-2))= (1/2) e-j(3/2)( F(-j(/2) 频移特性: ej4t f(3-2t) ( (1/2) e-j(3/2)((-4) F(-j((-4)/2). 六、频移特性(调制特性) 若f(t) ( F(j() (0为常数 f(t)e(j(0t ( F[j(((0)]. 证:[f(t) ej(0t]= f(t) ej(0t(e-j(tdt=f(t)(e-j((-(0)tdt= F[j((-(0)]. 应用:频域搬移技术在通信系统中得到广泛应用,诸调频、同步解调、变频等过程都是在频谱搬移的基础上完成的。频谱搬移的实现原理是将信号乘以载频信号cos(0t或sin(0t。 应用例子:y(t)= g((t) cos(0t y(t)= (1/2)g((t) (e-j(0t+(1/2)g((t)(ej(0t Y(j()=((/2)Sa[((+(0)(/2)+((/2)Sa[(-(0)(/2]. 表示:已调信号的频谱是将原频谱一分为二,分别向左和向右搬移(0,在搬移中幅度谱的形状并未改变。  应用:收音机,分成各个波段,将声音调制在不同的频段上。 七、卷积定理 1、时域卷积定理:f1(t)( f2(t)( F1(j()F2(j() 证:[ f1(t)( f2(t)]= [ f1(()f2(t-()d(](e-j(tdt = f1(()[ f2(t-()(e-j(t dt]d( = f1(()(e-j(t F2(j()d( =F2(j() f1(()(e-j(td(= F2(j()F1(j() 2、频域卷积定理:f1(t)f2(t) ((1/2() F1(j()(F2(j() =(1/2() F1(j()F2(j(-j()d(. 例:4.5-7 r(t)=tv(t)的频谱函数 解:1、t的频谱 t(j2(((((() 2、((t)的频谱 ((t) (((((()+1/(j() 3、由频域卷积定理 [t((t)]= (1/2()[ j2(((((()]([((((()+1/(j()] = j(((((()(((()+(((()((1/() = j(((((()+( (()((1/()(= j(((((()+ j(((((()-1/(2 ∴ t((t) (j(((((()-1/(2 可推出:(t((-2/(2 八、时域微积分 1、微分:f(t) ( F(j() f((t) ( j(F(j(). f(n)(t) ( (j()n F(j() 证:f((t)=f(t)*(((t) [ f((t)]= [ f(t)] ([(((t)]= F(j()(j(. 2、积分: f(-1)(t) ( (F(0) (((()+(1/j()F(j() f(-1)(t) ( (1/j()F(j() 当F(0)=  f2(t)dt=0时 证:f(-1)(t)= f(-1)(t)* ((t)= f (t)* ((-1)(t)= f(t)*((t). [ f(-1)(t)] = [ f(t)] ([((t)]= (F(0)(((()+(1/j() = (F(j()(((()+(1/j() F(j()=(F(0)(((()+(1/j() F(j() 如果F(0)=0即F(0)=  f2(t)dt=0则f(-1)(t) ( (1/j()F(j() 例:4.5-8 仍求f((t)的频谱函数。 解:f((t)( f(((t)( f((((t).冲激函数变换 F(j()   f((t)= f(-2)(t) 为双重积分。 f(t)=(2/()((t+(/2)- (4/()((t)+ (2/()((t-(/2) F(j()=(2/() ej((/2- (4/()+ (2/() e-j((/2 =(4/()[cos(((/2)-1]=-8sin2(((/4)/(. 要用积分公式限判断F(0)=  f((t)dt=0, 因为双重积分,判断F(0)=  f( (t)dt=0, F(0)=  f(((t)dt=0. F((j()=(1/j()2 F(j()=8sin2(((/4)/((()2= (2/()Sa2(((/4). 九、频域微积分 1、微分:f(t) ( F(j() -jtf (t) ( F((j(). (-jt)n f(t) ( F(n)(j() 2、积分: (f(0)(((t)-(1/jt)f(t)( F(-1)(j() ∴ -(1/jt)f (t) ( F(-1) (j() 当f(0)= 0时. 例:4.5-10. 求t((t) 的付氏变换。 解: ((t) (((((()+1/(j() -jt ((t) (d[((((()+1/(j()]/d(=(((((()-1/(j(2). ∴t((t) (j(((((()+1/(2 例:4.5-11 求 Sa(t)=sint/t的付氏变换。 f(t)=sint=(1/2j)(ejt-e-jt) ( (1/2j)[2((((-1)-((( +1)] =j([((( +1)-(((-1)]=F(j() ∵f(0)=0 ∴sint/(-jt) ( F(-1)(j()=j([(((+1)-(((-1)]d( = j([((( +1)-(((-1)] ∴ sint/t(([((( +1)-(((-1)]= (g2(()。 十、能量谱和功率谱: 1、能量谱 信号的能量:E= f2(t)dt= f(t)([(1/2() F(j()(ej(t d(]dt =(1/2() (F(j()(2 d(. 能量谱: ξ(()=(F(j()(2 单位频率的信号能量,能量密度函数。 E= f2(t)dt =(1/2()ξ(()d(. 2、功率谱 信号的功率:P=  (1/T)f2(t)dt =(1/2()[((F(j()(2)/T]d(. 功率谱: ((()=((F(j()(2)/T单位频率的信号能量,功率密度函数。 P= (1/2()((()d(. § 4.6 周期信号的付里叶变换 信号分析; ( 周期信号:可展开为付里叶级数(求和),频谱Fn是离散的,满足狄里赫利条件。 非周期信号:存在付里叶变换(积分),频谱密度F(j()是连续的, |f(t)|dt<(,可放宽。 目的:统一分析方法。 一、正余弦函数的付里叶变换(典型的周期信号) 1(2(((() ej(0t(2((((-(0) e-j(0t(2((((+(0) cos(0t=(1/2)(ej(0t+e-j(0t)(([(((-(0)+ (((+(0)] sin(0t=(1/2j)(ej(0t–e-j(0t)(([(((+(0)- (((-(0)] 二、一般周期函数的付里叶变换:fT(t),周期为T. 方1、先求Fn,再求F(j()=[ fT(t)] fT展开成付里叶级数:fT(t)= FnejnΩt Fn=(1/T)f(t)e-jnΩtdt. F(j()=[ fT(t)]= [ FnejnΩt]=(2() Fn(((-n() 含义:有无穷多个冲击函数组成,位置在n(处,强度2(Fn. 例4.6-1  解:先求 Fn=((/T)Sa(n((/2) F(j()=[ PT(t)]=(2()((/T)Sa(n((/2) (((-n() = [2sin(n((/2)/n](((-n() 比较 fn 与 F(j()图形相位似,含义不同。 Fn是虚指数分量的幅度和相位;F(j()是频谱密度。 例4.6-2   求(T(t)= ((t-nT)的付里叶变换 解:先求Fn=(1/T) (T(t)e-jnΩtdt=(1/T) ((t)e-jnΩtdt=1/T 再求F(j()=[(T(t)]=(2()Fn(((-n()=(2(/T)(((-n() = ( (((-n()=(((() ∴ (T(t)( (((((). 方2、先求第一周期函数f0(t) ( F0(j(),再求F(j().   先求第一周期函数f0(t) ( F0(j(), fT(t)= f0(t)* ((t-nT)= f0(t)*(T(t) 再求 F(j()=[fT(t)]= [f0(t)]˙[(T(t)] = F0(j()((((()= (  F0(j()(((-n() =(  F0(jn()(((-n() 例4.6-3 P0(t)=g((t) 的付里叶变换.  解:先求[P0(t)]= (Sa(((/2) 再求: [PT(t)]= ( (Sa(n((/2) (((-n() =(2()((/T) Sa(n((/2) (((-n() 三、付里叶系数与付里叶变换 比较 [fT(t)]= (2()Fn (((-n() [fT(t)]= (2(/T)F0(jn()(((-n(). 关系:Fn=(1/T) F0(jn()=(1/T) F0(j() 表明:付里叶变换中的许多性质,定理也可用于付里叶级数,并提供一种求Fn的方法。 例4.6-4:将fT(t)展开成指数形式付里叶级数。  解:∵f1(t)( F1(j()=(T/2)S2a((T/4) f0(t)( F1(j()e-j(T/2)(=(T/2)S2a((T/4) e-j(T/2)(= F0(j() ∴Fn=(1/T) F0(j()=(1/2) S2a(n(T/4)e-j(T/2)n( = (1/2) S2a(n(/2)e-jn( fT(t)= (1/2) S2a(n(/2)e-jn(e-jn(t § 4.7 LTI系统的频域分析 一、频域响应: f(t)( LTI系统 ( y(t)=yf(t) ej(t ( LTI系统 ( H(j() ej(t f(t)=(1/2()F(j()ej(td(( LTI系统 (y(t)= (1/2() H(j()F(j()ej(td( 1、频域响应的定义;当输入为f(t)= ej(t 时 系统的零状态响应y(t)=h(t)*f(t)=h(t)* ej(t = h(()ej((t-()d(=h(()e-j((d(ej(t= H(j()ej(t 定义:H(j()=h(()e-j(td(为频域响应函数/系统函数; 关系:h(t)( H(j() 用来描述系统的特性;表示H(j()=(H(j()(ej((() h(t)描述时域特性 ; H(j()描述频域特性 ; 2、频域分析的基础方法:当激励为任意信号f(t)时 f(t)( F(j() 系统函数h(t)( H(j() y(t)( Y(j()= H(j()F(j() 时域分析 f(t)*h(t) = y(t)    频域分析 的关系: H(j()F(j()=Y(j() 非周期信号激励的系统响应: Y(j()= H(j()F(j() y(t)=[ Y(j()] 周期信号激励的系统响应: Y(j()= H(j()F(j()= H(j() (2()Fn(((-n() =(2() Fn H(jn()(((-n() y(t)= [ Y(j()]=[(2()  H(jn()Fn(((-n()] = Fn H(jn()ejnΩt=YnejnΩt ∴ yn= H(jn()Fn y(t)= ynejnΩt 例4.7-1:已知H(j(),f(t)=2+4cos(5t)+4cos(10t)求y(t).  解:f(t)是周期信号,基波(=5rad/s (1)用付里叶级数法 yn= H(jn()Fn= F(j()Fn. f(t)=2+4cos((t)+4cos(2(t) =2+2ejΩt+2e -jΩt+2ej2Ωt+2e-j2Ωt=2ejnΩt H(j()= |H(j()|ej( (() H(jn()= |H(jn()|ej( (n() yn= H(jn()Fn y0= H(0)F0 =1(ej0(2=2 y1= H(j()F1=(1/2)(e-j((/2)(2= e-j((/2) y2= H(j2()F2=0 y-1= H(-j()F-1=(1/2)(ej((/2)(2= ej((/2) y-2= H(-j2()F-2=0 y(t)= ynejnΩt =y-1e-jnΩt+y0+y1ejnΩt=e-j(5t-(/2)+2+ej(5t-(/2) =2+2cos(5t+900). (2)用付里叶变换法 F(j()=4(((()+4([(((+()+(((-()]+4([(((+2()+(((-2()] = 4((((-n(). Y(j()=H(j()4((((-n()=4(  H(jn()(((-n() = 4([0.5 ej(((+()+((()+0.5 e-j(((-()]. y(t)= [ Y(j()]= e-j(n (-)+2+ ej(n (-)=2+2cos((t-900). 例:4.7-4 已知f(t)=[sin(2t)]/t,s(t)=cos(3t), 系统H(j()= 1 |(|<3 0 |(|>3, 求y(t).    H(j() ( y(t) ( ( ( ( 调制信号 ((t)载波信号 已调信号 传输   解:从频域角度考虑系统在传输过程中的影响。 ①求 X(j(). x(t)=f(t)(s(t) X(j()=(1/2()F(j()*S(j() ∵f(t)=2Sa(2t) ( F(j()=(g4(() s(t)=cos(3t) ( S(j()=([(((+3)+(((-3)] ∴X(j()=(1/2()(g4(()*([(((+3)+(((-3)]=((/2)[g4((+3)+g4((-3)]. ②系统 H(j()=g6(() ③ Y(j()= X(j()( H(j()=((/2)[g2((+2)+g2((-2)] = (1/2()(( g2(()([(((+2)+(((-2)] ④ y(t)=Sa(t)(cos(2t)=(sint/t)(cos2t. 分析:传输中有失真。 3、求频域响应的方法: 由定义:H(j()=h(t)e-j(tdt 由微分方程:y((t)+a1y((t)+a0y(t)=f(t) 两边取付氏变换:(j()2Y(j()+ a1(j()Y(j()+a0Y(j()=F(j() ∴ H(j()=Y(j()/X(j()=1/[ (j()2+ a1(j()+a0]. 向量法(简单电路) 由电路分析可知,电路在正弦稳态作用下可应用向量法分析,此时电路有阻抗特征或导纳特征,且其阻抗和导纳均为的复函数。 H(j()=Y(j()/X(j()= 例4.7-2 已知:y((t)+2y(t)=f(t),f(t)=e-t((t).求y(t). 解:方程取付氏变换: j(Y(j()+2Y(j()=F(j(). H(j()= Y(j()/ F(j()=1/(j(+2). f(t) =e-t((t)( 1/(j(+1)= F(j(). Y(j()= H(j() F(j()=1/[(j(+2)(j(+1)]= 1/(j(+1)- 1/(j(+2). ∴y(t)= [Y(j()]= e-t((t)- e-2t((t). 例:4.7-3已知RC电路. u s(t)= ((t).求u c(t). 解:向量法:H(j()=Uc(j()/Us(j()=(1/j(c)/[R+(1/j(c)] =[1/(RC)]/[j(+(1/RC)]=(/((+j(). 激励:u s(t)= ((t)( Us(j()=(((()+1/(j(). Uc(j()= H(j()Us(j()=[(/((+j()][(((()+1/(j()] = [((/((+j()]((()+ (/[((+j()(j()]. = (((()+1/(j()- 1/((+j(). u c(t)= [Uc(j()]=1/2+(1/2)Sgn(t)-e-(t((t)=(1- e-(t)((t). 二、无失真传输 引起信号失真的根本原因是系统函数。包括两个方面: 幅度失真:(H(j()(对信号频率分量加权不同,使信号幅度产生相应变化。 相位:((() 对信号频率分量附加了不同相移,而是合成波形也产生了变化。 1、定义:y(t)=Kf(t-td). 2、频域条件:H(j()应满足的条件 Y(j()=Ke-j(td F(j()  H(j()=Ke-j(td (H(j()(=K全部频带内,幅频特性为一常数;读不同频谱分量,加 权相同。 ((() =-(td 向频特性为原点的直线,是线性的,斜率为- td。 3、时域条件:h(t)= [H(j()]=[ Ke-j(td]=K((t-td). 应为冲激函数的K倍,延时td。(以上为理想条件。) 三、理想低通滤波器 1、定义:( H(j()( = 1 ((( < (c 0 ((( > (c ((() =-(td 系统函数:H(j()= e-j(td ((( < (c 0 ((( > (c = e-j(td g2(td((). 宽度为2(c,幅度为1的门函数。 2、冲激响应 h(t)= [e-j(td g2(()]=((c/()Sa[(c(t- td)] = ((c/()(sin[(c(t- td)]/ (c(t- td)) 特点: 中心:td 峰值:(c/( 零点间隔(/(c t<0时,h(t)(0,为非因果系统。 3、阶跃响应 g(t)(G(j(). G(j()= H(j()[v(t)] = [(((()+(1/j()] e-j(td ((( < (c 0 ((( > (c g(t)= (1/2()G(j()ej(td( = (1/2()[(((()+(1/j()] e-j(td ej(td( = 1/2+(1/2() e-j((t-td)/(j()d( = 1/2 + (1/()Si[(c(t- td)].其中Si(x)=(sin()/(d( ①阶跃响应的斜率: dg(t)/dt=h(t)= (c/(。 ②信号的上升时间tr,定义为t=td处斜率的倒数。 Tr = (/(c =(/2(fc =0.5/fc=0.5/8. 通带愈宽,上升时间愈短,波形愈陡直,上升时间与通带宽度成反比。 ③基波斯现象:间断点处过冲。 第一个极大值发生在t=td+(/(c处 gmax =1/2+(1/()Si[(c(t-td)]=1/2+(1/()Si(()=1.0895.幅度为9%. 结论:与(c无关,不可能通过增大通带宽度来改变。 四、实际低通滤波器: 1、传输函数H(j() H(j()= UR(j()/ US(j()=1/[1-(2LC+j((L/R)] =1/[1-((/(c )2+j((/(c )] (c =,R= 幅度特性:H(j()= 1/ 相频特性; ((()= -arctg(((/(c )/[1-((/(c )2]). H(±(c )=1/ ((±(c)= ±(/2.  与理想滤波器有相似之处,(-(c ,(c )阶数越高,逼近程度越好。 2、冲激响应;h(t)=[ (cesin()]((t). 阶跃响应: g(t)=[1-esin(+(/4)]((t). 五、网络可实现的条件 1、时域条件: h(t)=0 t<0 g(t)=0 t>0. 2、频域条件:  ( H(j()(d(<( 无失真,不满足. 且 [(ln ( H(j()((]/(1+(2)d(<( 佩利——维钠准则或定理. §4.8 取样定理 一、信号的取样:连续信号 离散信号,以便用计算机处理。 1、抽样的模型: f(t)((fs(t) ( ((t) fs(t)=f(t) ((t) Ts:抽样周期; fs=1/Ts 抽样频率;(s=2(fs抽样角频率. 冲激抽样 2、冲激抽样 时域(频域变化。  S(t)= (T(t)= (Ts(t) 2([ (Ts(t)]= ((((()=(s ((s(()=(s  (((-n(s). 表示:有无限个冲激序列组成,间隔(s,强度(s f(t)( F(j()带限信号. fs(t)=f(t)(Ts(t). Fs(j()=(1/2()F(j()*(s  (((-n(s) =((s/2() F(j()* (((-n(s)=(1/Ts) F[j((-n(s)]. 表示:是原信号频谱F(j()的周期性重复,重复周期:(s强度:1/Ts. 频谱混叠想象: 当(s ( 2(m或Ts(1/(2fm) 时,无混叠现象,可从fm(t)中恢复f(t)。 当(s (2(m或时,有混叠,无法分开,不能恢复f(t). 3、矩形脉冲抽样,时域(频域  令s(t)是幅度为1,宽度为(的门函数序列. S(t)=PT(t)= g((t-nTs). 由例4.6-1 P(j()=(2((/T) Sa(n((/2)(((-n() = (2((/Ts) Sa(n(s(/2)(((-n(s) 表明:它是冲激函数序列,周期为(s,包络线是(2((/Ts)Sa(((/2). f(t)( F(j() 频带受限信号,当(( ( ( (m时,不为0. fs(t)=f(t)PT(t). Fs(j()= (1/2()F(j()*P(j() =((/Ts)F(j()*[Sa(n(s(/2)(((-n(s)] = ((/Ts) [Sa(n(s(/2)F(j()*(((-n(s)] = ((/Ts) [Sa(n(s(/2)F(j((-n(s)].。 表示:有F(j()无限个频移组成,间隔(s,包络线((/Ts)Sa(((/2). 混叠现象: 若(s( 2(m可恢复,不混叠。 若(s ( 2(m不可恢复,混叠。g(t) t td tr 0 0.5 1.0895 -0.0895  三、时域抽样定理: 两个条件:1、信号f(t)必须是带限信号,在F(j()时( ( (>(m为零。 抽样频率; 2、(s( 2(m;抽样间隔不能太长Ts(1/(2fm). 满足时,可从fm(t)中恢复f(t),通过理想低通滤波器((c=(s/2)。 最低抽样频率fs= 2fm 奈奎斯特频率。 最大抽样间隔Ts=1/(2fm) 奈奎斯特间隔。 二、信号的恢复:离散信号fs(t)(连续信号f(t). 理想低通滤波器 1、理想低通滤波器的选择。    (c: (m( (c ((s/2. 幅度:Ts 相位: ((()=0,td=0. H(j()= Ts ((( <(c 0 ((( >(c 2、理想低通的h(t) h(t)=  = . 取(c =(s/2, Ts=2(/(s=(/(c ( h(t)=Ts((1/Ts)(Sa((st/2). 表示:f(t)可展开为正交抽样函数[Sa函数]的无穷级数 系数为f(nTs),各点的抽样值。 在采样点,f(nTs)与f(t)完全相等,采样点间由Sa函数的无穷波形叠加而成,无失真。 四、频域抽样定理:对称性。 频域( 时域     Fs(j()= F(j()((s(() fs(t)= f(t)* [((s(()]. ∵((t-nTs) ( (s(((-n(s)= (s((s((). ∴(1/(s)((t-nTs) (((s((). fs(t)= f(t)* (1/(s)((t-nTs)= (1/(s)f(t)*((t-nTs) = (1/(s)f(t-nTs).表示时域信号以为Ts周期重复,幅度为1/(s . 混叠现象: 若选Ts<2tm [1/Ts= fs<1/(2tm)],则fs(t)无混叠,可用选通信号无失真恢复。 若选Ts>2tm,会产生混叠,无法恢复。 频域取样定理:一个在时域区间((-tm ,tm )以外为零的有限时间信号f(t)的频谱函数F(j(),可惟一地由其在均匀频率间隔fs( fs<1/(2tm))上的样点值F(jn(s)确定。 2个条件: 有限时间信号f(t),(-tm ,tm )间有值。 取样频率fs<1/(2tm),可恢复。