1数字信号处理第0章绪论1.数字信号处理是利用计算机或专用处理设备,以数字形式对信号进行采集、变换、滤波、估值、增强、压缩、识别等处理,以得到符合人们需要的信号形式。2.DSP系统构成输入抗混叠滤波A/DDSP芯片D/A平滑滤波输出输入信号首先进行带限滤波和抽样,然后进行A/D(AnalogtoDigital)变换将信号变换成数字比特流。根据奈奎斯特抽样定理,为保证信息不丢失,抽样频率至少必须是输入带限信号最高频率的2倍。DSP芯片的输入是A/D变换后得到的以抽样形式表示的数字信号。3.信号的形式(1)连续信号在连续的时间范围内有定义的信号。连续--时间连续。2(2)离散信号在一些离散的瞬间才有定义的信号。离散--时间离散。4.数字信号处理主要包括如下几个部分(1)离散时间信号与系统的基本理论、信号的频谱分析(2)离散傅立叶变换、快速傅立叶变换(3)数字滤波器的设计第一章离散时间信号一、典型离散信号定义1.离散时间信号与数字信号时间为离散变量的信号称作离散时间信号;而时间和幅值都离散化的信号称作为数字信号。32.序列离散时间信号-时间上不连续上的一个序列。通常定义为一个序列值的集合{x(n)},n为整型数,x(n)表示序列中第n个样值,{·}表示全部样本值的集合。离散时间信号可以是通过采样得到的采样序列x(n)=xa(nT),也可以不是采样信号得到。二.常用离散信号1.单位抽样序列(也称单位冲激序列))(n0,00,1)(nnnδ(n):在n=0时取值为142.单位阶跃序列)(nu,0,00,1)(nnnu3.矩形序列,其它nNnnRN,010,1)(4.实指数序列,)()(nuanxn,a为实数55.正弦型序列)sin()(nAnx式中,ω为数字域频率,单位为弧度。15On1100sinnωt0sin16.复指数序列njenx)(0)(7.周期序列如果对所有n存在一个最小的正整数N,使下面等式成立:)()(Nnxnx,则称x(n)为周期序列,最小周期为N。三.序列的运算(一)基本运算离散时间信号(即序列)的基本运算包括移位、反折、求积、乘积、差分运算和尺度变换,下面分别介绍。1.序列的移位设某一序列为()xn,若0m,则()xnm表示序列()xn整体右移了6m个样点形成的新序列,也称()xnm是()xn的m个样点的延迟。此时()xnm表示序列()xn整体左移了m个样点形成的新序列,也称()xnm是()xn的m个样点的超前。例如,()xn如图(a)所示,则(2)xn和(2)xn分别如图(b)和图(c)所示。序列的移位2.序列的反折设某一序列为()xn,则()xn是以0n为对称轴将序列()xn水平翻转,()xn称为序列()xn的反折。若()xn如图(a)所示,则()xn如图(b)所示。序列的反折3.序列的求和()xn与()yn两个序列之和是指两个序列同序号(即n相同)的序列值逐项对就相加构成一个新的序列()zn,表示为()()()znxnyn【例1】已知1(),1()20,1nnxnn,1(),0()21,0nnynnn,求()()xnyn7解:根据序列求和定义,得11(),02()()()2,11,1nnznxnynnnn()xn、()yn和()()xnyn的图形分别如图(a)、(b)和(c)所示。序列的求和4.序列的乘积()xn与()yn两个序列的乘积是指两个序列同序号(即n相同)的序列值逐项对就相乘构成一个新的序列()zn,表示为()()()znxnyn【例2】已知1(),1()20,1nnxnn,1(),0()21,0nnynnn求()xn与()yn的乘积()()xnyn。解:根据序列乘积的定义可得1(),0()()()40,0nnznxnynn另外,序列还可以与标量相乘,定义为()()ynxn8序列()xn与标量相乘相当于()()ynxn中每一个样值是相同序号(即n相同)的()xn样值的倍。5.尺度变换(1)抽取序列)(mnx是)(nx序列每隔m点取一点形成的,相当于时间轴n压缩了m倍。(2)插值序列mnx是)(nx序列相邻抽样点间补(m-1)个零值点,表示零值插值。(二)其它计算1.累加(等效积分)nkkxny)()(2.差分运算(1)前向差分)()1()(nxnxnx(2)后向差分)1()()(nxnxnx93.卷积mmnhmxnhnxny)()()()()(设x(n)和h(n)两序列的长度分别是N和M,线性卷积后的序列长度为(N+M-1)。四.时域离散系统1.线性系统:满足叠加原理的系统。判定公式:则:例:)792sin()()(nnxny是线性系统。证:)792sin()()(11nnxny)792sin()()(22nnxny)792sin()]()([)]()([22112211nnxanxanxanxaT)792sin()]()([)()(22112211nnxanxanyanya由上式得:)()()]()([22112211nyanyanxanxaT2.时不变系统(移不变系统):判定公式:若[()]()TxnYn,则[()]()TxnkYnk例:bnaxny)()(10证:bnnaxnnxT)()]([00bnnaxnnxT)()]([00bnnaxnny)()(00由上式得:)]([)(00nnxTnny3.线性移不变系统设系统的输入序列为x(n),它可以表示为单位取样序列的移位加权和,即:mxnxmnm那么,系统对应的输出为:[]mynTxnTxmnm如果该系统是一线性移不变系统,根据其线性则有:mmynTxmnmxmTnm又根据移不变性和h(n)定义,则有:冲激响应:()[()]hnmTnm所以此时系统输出为:()()*()ynxnhn,()()()YjXjHj,()()()YzXzHz4.稳定系统:有界的输入产生的输出也有界的系统,即:若|()|xn,则|()|yn线性移不变系统是稳定系统的充要条件:|()|nhn5.因果系统:0n时刻的输出0()yn只由0n时刻之前的输入0(),xnnn决定线性移不变系统是因果系统的充要条件:()0,0hnn(1)对于连续时间系统:t=t1的输出y(t1)只取决于t≤t1的输入x(t≤t1)时,则此系统为因果系统,11(2)对于离散时间系统:n=n1的输出y(n1)只取决于n≤n1的输入x(n≤n1)时,则此系统为因果系统,例:(1)y(t)=x(sin(t))不是因果系统,因为y(-π)=x(0),表明y(t)在一段时间内可能取决于未来的x(t)。(2)y(t)=x(t)cos(t+1)是因果系统,cos(t+1)是时变函数,相当于一个已知的函数波形,所以x(t)的当前值影响了y(t)的当前值。6.稳定因果系统:同时满足上述两个条件的系统。线性移不变系统是因果稳定系统的充要条件:|()|nhn,()0,0hnn或:H(z)的极点在单位园内H(z)的收敛域满足:||,1xxzRR第二章采样定理连续时间信号的处理往往是通过对其采样得到的离散时间序列的处理来完成的。1.实际抽样与理想抽样12xa(t)ot(a)(b)xa(t))(ˆatxTp(t)tttt(c)(e)(d)(f)s(t)xp(t))(ˆatxooooT1T)()()(ˆtstxtxaa)()()(ˆnTtnTxtxana2.理想采样信号的频谱)()()(ˆtstxtxaadtetsjSdtetxjXtjtjaa)()()()()()(21)(ˆjSjXjXaakaaTjkjXTjX)2(1)(ˆ130Ωc-ΩcXa(jΩ)P(jΩ)-ΩsΩsΩΩ0Xa(jΩ)Ω0Xa(jΩ)ΩΩcΩs(a)(b)(c)(d)^^2s0-ΩsΩs-Ωsδ2s2s由上图可知:当Ωs≥2Ωc时,采样信号的频谱不产生重叠,通过理想滤波器,信号可以无失真地复现出来,这种情况称过采样;如果Ωs2Ωc,采样信号的频谱出现重叠,信号产生失真,称不足采样,这种频谱的重叠称混淆。3.抽样的恢复2||02||)(ssTjG140Xa(jΩ)Ω^G(jΩ)xa(t)ya(t)0G(jΩ)Ω-π/Tπ/T0Xa(jΩ)Ω(a)(b)(c)(d)^4.奈奎斯特取样定理(1)对连续信号进行等间隔采样形成采样信号,采样信号的频谱是原连续信号的频谱以采样频率为周期进行周期性的延拓形成的。(2)设连续信号xa(t)属带限信号,最高截止频率为Ωc,如果采样角频率Ωs≥2Ωc(fs≤1/2fc),那么让采样信号x^a(t)通过一个增益为T,截止频率为Ωs/2的理想低通滤波器,可以唯一地恢复出原连续信号xa(t)。否则Ωs2Ωc会造成采样信号中的频谱混叠现象,不可能无失真地恢复原连续信号。ttya0cos)((3)不足采样情况如果Ωs2Ω0,其中Ωs是采样频率,Ω0是输入信号频率。例如Ω0=6π,Ωs=10π可以得出Ωs=5/3Ω0Ωs-Ω0=2/3Ω015ttysa)cos()(0于是,复现信号为)(tya=cosΩ0t2/3=cos6πt-00Ya(j)Ya(j)-(s-0)ππππ混叠无混叠(e)(d)0>T0<Toos-0(4)折叠频率Ts2为防混叠滤波器,低通滤波器的截止频率为Ωs/2。第二章信号与信号频谱一.信号的定义及种类信号的概念广泛出现于各领域中。这里所说的均指电信号,一般可表示为一个或多个变量的函数。按照信号随时间变化的特点,可分为确定信号与随机信号连续时间信号与离散时间信号周期信号与非周期信号16其它分类如:奇信号与偶信号,调制信号与载波信号,能量有限信号与功率有限信号……二.信号相关分析原理1.信号的互能量与互能谱(1)信号的能量:指信号f(t)的归一化能量,即信号的电压(电流)加在1欧电阻上所消耗的能量。dttfE2|)(|(2)能量谱与功率谱能量谱:dFdttfE22)(21)(其中|F(ω)|2表明了信号能量在频域的分布情况,所以被称为能量谱密度,简称能谱。记作:2)()(FW功率谱:设)(0tfT是)(tf的截短函数220000)()(TTTtttftf则f(t)的功率谱密度函数为02)(lim)(00TFSTTdSP)(21(3)两信号的互能量两信号x(t)、y(t)之和的能量为:17信号的互能量为:dttytxExy)()(2(4)广义瑞利公式、互能谱广义瑞利公式:若信号x(t)和y(t)为实函数,其频谱密度分别为)()(YX和,则dYXdttytxyx)()(21)()(),(互能谱:)()()(YXWxyWxy(ω)称为信号x(t)、y(t)的互能谱密度,简称互能谱。2.信号的自相关函数(1)定义:为了定量地确定信号x(t)与时移副本x(t-τ)的差别或相似程度,通常用自相关函数:dttxtxRx)()()((2)特点:自相关函数是偶函数)()(RR当τ=0时,自相关函数等于信号的能量xxEdttxR