南京理工大学数字信号处理DSP复习VIP

1、复习目录1 一些基础11.1 卷积、环形卷积11.2 xcs/xca、xev/xod21.3 周期、抽样率、抽样定理、0与的关系21.4 线性、时不变、因果、稳定31.5 时域频域信号对应关系41.6 全通系统极点关于单位圆镜像对称42 三种变换42.1 dtft42.1.1 定义42.1.2 性质52.1.3 常用52.2 z变换62.2.1 定义62.2.2 性质82.2.3 常用82.3 dft92.3.1 定义92.3.2 性质112.3.3 fft算法：蝶形图、指标123 两套系统133.1 差分方程y（n）时域描述133.1.1 建立/解差分方程133.2 系统函数h（z）频域描述

2、143.2.1 iir滤波器信号流图结构、低通iir设计、脉冲响应/双线性不变法143.2.2 fir滤波器信号流图结构、线性相位系统、窗函数法、四种类型181 一些基础1.1 卷积、环形卷积卷积：y（n）=x（n）*h（n）eg:卷积 xn=1,2,3,4,5hn=1,1,1,1xn*hn=1,3,6,10,14,12,9,5 表格对角线斜向相加,长度为5+4-1=8h【n】 x【n】12345112345112345112345112345 环形卷积：1）设序列h(n)和x(n)的长度分别为n和m。h(n)与x(n)的l点循环卷积定义为式中，l称为循环卷积区间长度，lmaxn，m。2) 循

3、环卷积矩阵eg:环形卷积 （1）计算下面给出的两个长度为4的序列h(n)与x(n)的4点和8点循环卷积。 hn=1,1,1,1与xn=1,2,3,4的循环卷积：注意，hx位置在前，且循环卷积的长度与循环卷积区间长度l相同 解： 4点循环卷积矩阵形式为 8点循环卷积矩阵形式为矩阵计算：行x列，8点第一行x第一列：1x1+0x1+0x1+0x1+0x0+4x0+3x0+2x0=1 第二行x第一列：2x1+1x1+0x1+0x1+0x0+0x0+4x0+3x0=3（2）xn=3n,hn=2+(-2)n,0n3.求6点圆周卷积xn=1,3,9,27 hn=3,0,6,610027933-1593100

4、2790993100276332793100-636027931009300279310108 1.2 xcs/xca、xev/xodxn：对象序列 x*：共轭运算实不变，虚取反xn = xcs + xca = xev + xodxcs共轭对称序列：xcs满足xn=x*-n若xn无虚部，则xn为偶序列xevxcs=12（xn+x*-n） xev=12（xn+x-n）xca共轭反对称序列：xca满足xn=-x*-n若xn无虚部，则xn为奇序列xodxca=12（xn-x*-n） xod=12（xn-x-n）eg:共轭对称，共轭反对称，偶序列，奇序列 xn= -1+j3 , 2-j7 , 4-j5

5、 , 3+5j , -2-j , -2n2，求xcs、xca、xev、xodxcs=xca=xev=xod= 1.3 周期、抽样率、抽样定理、0与的关系周期：若正弦序列x(n)= sin(n)=sin(30n/120)是周期的,则周期n=2，抽样率：r=ftft抽样定理：一个限带模拟信号xa(t)，若其频谱的最高频率为f0，对它进行等间隔抽样而得x(n)，抽样周期为t，或抽样频率为fs=1/t；只有在抽样频率fs2f0时，才可由xa(t)准确恢复x(n)。0与的关系：0=0t，0：离散时间信号归一化角频率，0：连续时间型号角频率，t：抽样周期eg:抽样相关 (1)the sequence xn

6、=cos(4n) is obtained by sampling an analog signal x(t)=cos(0t) at a sampling rate of 1000 samples/sec.two possible value of 0 are 250，1750that could have resulted in the sequence xn.41000=250，(-4+2)1000=1750(2)有一连续信号xa(t)=cos(2ft+)（1）求出xa(t)的周期。（2）用采样间隔t=0.02s对xa(t)进行采样，试写出采样信号 的表达式。（3）求出对应 的时域离散信号(

7、序列) x(n)，并求出x(n)的周期。解：（1） 周期为 （2）（3）x(n)的数字频率=0.8，故 ，因而周期n=5，所以x(n)=cos(0.8n+/2) 1.4 线性、时不变、因果、稳定eg:线性、时不变、因果、稳定 a discrete-time system is described by the equation yn=nsin(xn+4)=22nsinxn+22ncos(xn)线性liner：判定条件：t(ax1n+bx2n) = at(x1n) +bt(x2n)t(ax1n+bx2n)= 22nsinax1n+bx2n+22ncos(ax1n+bx2n)= 22nsinax1

8、n)cos(bx2n+cosax1n)sin(bx2n+22ncosax1ncosbx2n-sinax1nsinbx2nother hand:at (x1n)= 22ansinx1n+22ancos(x1n) , bt(x2n)= 22bnsinx2n+22bncos(x2n)at(x1n) +bt(x2n)= 22ansinx1n+22ancos(x1n)+ 22bnsinx2n+22bncos(x2n)at(x1n) +bt(x2n) t(ax1n+bx2n),so not liner.时不变shift-invariant：判定条件：y(n-n0) = t(xn-n0)yn=t(xn) ,

9、 if the imput is xn-n0,t(xn-n0) = nsin(xn-n0+4)but ,y(n-n0) = (n-n0)sin(xn-n0+4)y(n-n0) t(xn-n0)，so not shift-invariant因果causal：判定条件：h(n)0 (n0) 或 roc: r1 |rx-|。5)对于反因果序列xn, 其 z变换的收敛域roc由其离原点最近的极点确定，其形式为|z|rx+|。eg:z变换 1）求z正换（1）；（2）；（3） 解：（1）（2）（3）2）求z逆变换(1) ，收敛域roc 为,则 的z反变换为( )。说明：本题要求掌握序列的时域特性域z变换收

10、敛域之间的对应关系。具体说，有限长序列的z变换的roc是怎样的，右边序列的z变换的roc是怎样的，因果序列的z变换的roc是怎样的，左边序列的z变换的roc是怎样的，反因果序列的z变换的roc是怎样的。典型序列的z变换表达式是否记住了？这两个典型z变换对，对求z变换或逆z变换非常重要。(2) ，试求与x(z)对应的所有可能的序列x(n)。解：同一个z变换函数，收敛域不同，对应的序列也不同。本题没有给定收敛域，所以必须先确定收敛域。x(z)有两个极点：z1=0.5,z2=2，因为收敛域总是以极点为边界，所以收敛域有以下三种情况：|z|0.50.50.6.求xn.xz=3zz2+0.3z-0.18

11、=1z-0.3+2z+0.6=1z11-0.3z-1+21+0.6z-1，时移性质gn-n0z变换z-n0g(z)xn=0.3n-1+2-0.6n-1n-1 2.2.2 性质名称性质描述收敛域共轭g*nz变换g*(z*)rg时间反转g-nz变换g(1z)1rg线性gn+hnz变换g(z)+h(z)rgrh时移gn-n0z变换z-n0g(z)rg，可能除了点z=0或z=与指数相乘ngnz变换g(z)|rgg(z)的差分ngnz变换-zdg(z)dzrg，可能除了点z=0或z=卷积gnhnz变换g(z)h(z)rgrh调制gn+hnz变换g(z)+h(z)rgrh2.2.3 常用序列nz变换1 |

12、z|0nz变换11-z-1 |z|1nnz变换11-z-1 |z|n-n-1z变换11-z-1 |z|nnnz变换z-1(1-z-1)2 |z|(n+1)nnz变换1(1-z-1)2 |z|2.3 dft2.3.1 定义xk=n=0n-1xnwnknxk=1nn=0n-1xnwn-kn其中wn=e-j2n , 0kn-1，周期性：wn(n+n)k=wnnk=wn(k+n)n，对称性：wnn+n2=-wnneg:dft (1)某dft的表达式是，则变换后数字频域上相邻两个频率样点之间的间隔是。(2)某序列dft的表达式是，由此可看出，该序列的时域长度是 n ，变换后数字频域上相邻两个频率样点之间

13、隔是。（3）若x (n)= 3，2，1，2，1，2 ，0n5，1) 求序列x(n)的6点dft，x (k)=？2) 若，试确定6点序列g(n)=?3) 若y(n) =x(n)x(n)，求y(n)=？解：1） 2） 3） （4）let the sequence xn=2,4,-2,1,3,1,-2,4,0n7，and its 8-point dft is xk,a)x0=n=07xn=11b)x4=x0-x1+x2-x3=-9c)k=07xk=8x0=16d) k=07|xk|2=8xn8点圆周卷积xnn=0=440e) k=07e-j(3k4)xk=8x8-3=8(5)compute the

14、4-point dft xn of sequence xn1,2,0,1xk=n=03xnw4nk=w40k+2w4k+0w42k+w43k，wnk=e-j2nk=cos(2nk)-jsin(2nk),欧拉公式：eix=cosx+isinxx0=w40+2w40+0w40+w40=1+2+0+1=4 x1=w40+2w41+0w42+w43=1+2(cos241-jsin241)+0+cos241-jsin243=1-2j+0+j=1-jx2=w40+2w42+0w44+w46=1+2(cos242-jsin242)+0+cos246-jsin246=1-2+0-1=-2x3=w40+2w43

15、+0w46+w49=1+2(cos243-jsin243)+0+cos243-jsin249=1-0+0+j=1+2j+0-j=1-j(5)compute the 4-point dft xn of sequence xn3,2,0,2xk=n=03xnw4nk=3w40k+2w4k+0w42k+2w43kx0=3+2+0+2=7x1=3-j+0+j=3x2=3-2+0-2=-1x3=3+j+0-j=3 2.3.2 性质名称描述线性gn+hndftgk+hk圆周时移gn-n0ndftwnkn0gk圆周频移wn-k0ngndftgk-k0n对偶gndftng-kn反转g-nndftg-kn与共轭

16、反转相加 xn+x*ndftxk+x*-kn=2xcsk共轭对称xk=x*n-keg:dft的性质 (1)给定一16-点实序列x（n）, 其 16-点 dft 记为x(k), 已知x(13)= 2 + j3,则x*(3) = ( 2 + j3 )。(2)xn+3ndftwn-3kxk x2ndft1nxk圆周卷积xk x-nndftx-kn e-j8nnxndftxk+4n xn+x*ndftxk+x*-kn=2xcsk(3)consider a length-n sequence xn whith an n-point dft xk,n is an even number.1.yn=(-1)

17、nxn,determine the n-point dft yk (0kn-1)yn=-1nxn=ejxn=ej2nnn2xnyk=xn , 0kn-12.wn=xn,n is even（偶）0,n is odd(奇)，determine the n-point dft wk (0nn-1)wn=12xn+(-1)nxnwk=12xk+xn , xn3.zn=xn+xn+n2, determine the n2-point dft zk (0nn2-1)zn=x2k , 0nn2-1 2.3.3 fft算法：蝶形图、指标计算n=2m个点复加次数复乘次数fftnlog2nn2log2ndftn(

18、n-1)n2p4443 两套系统3.1 差分方程y（n）时域描述3.1.1 建立/解差分方程对lti系统1.将差分方程转化为h(z),注意roc。2.求解h(z),逆z变换得到h(n)eg:解差分方程 (1)a causal lti system is described by the following difference equationyn+0.2yn-1-0.15yn-2 = 2xn1)determine the impulse response hn of this system2)if the imput is xn= n-.05n-1,determine the output

19、sequence yn解：hz=21+0.2z-1-0.15z-2=2（1-0.3z-1）（1+0.5z-1),|z|0.5hz=341-0.3z-1+341+0.5z-1，hn=145-0.5n+30.3nnyn=145-0.5n+30.3nn-185-0.5n-1+30.3n-1n-1（2）yn-yn-1-yn-2=xn-11.h(z) and rochz=z-11-z-1-z-2=z(z-1-52)(z-1+52)roc:|z|1+522.inpulse responsehz=z-11-z-1-z-2=z-1(1-1-52z-1)(1-1+52z-1)=55(11-1+52z-1-11-

20、1-52z-1)nnz变换11-z-1 |z|hn=551+52n-1-52nn3.stable?roc包含单位圆n-n-1z变换11-z-1hn=551+52n-n-1-1-52nn(3)yn=xn-0.81yn-21.determine the system function h(z)2.sketch hej and ()解：复数z=a+bi,模为a2+b2，幅角w=arctanba,然后用（a,b）的象限确定w的值1.hz=11+0.81z-2=1(1-0.9z-1)(1+0.9z-1) ,(z0.9) (注意题目中式子的形式)2. hz=11+0.81z-2带几个关键点，z=ej=0，

21、2，32,算出z=a+bj，带入h(z)计算 hej () 3.2 系统函数h（z）频域描述3.2.1 iir滤波器信号流图结构、低通iir设计、脉冲响应/双线性不变法iir : hz=p0+p1z-1+p2z-2+p3z-31+d1z-1+d2z-2+d3z-3所有零点都在单位圆内叫最小相位，反之叫最大相位信号流图结构：直接型、直接、级联型、并联型eg:由h(z)画结构流图 hz=1+2z-1+2z-2+z-31-2z-1+z-3=1+-61-z-1+6+4z-11-z-1-z-2(并联)=1+z-11-z-11+z-1+z-21-z-1-z-2(级联)直接 直接级联 并联 低通iir滤波器

22、设计p348通带截止频率p 阻带截止频率s 3db通带截止频率c 通带最大衰减p=-20lg1-pdb 阻带最小衰减s=-20lg(s) 通带波纹峰值p=11+2(?) 阻带波纹峰值s=1a(?) 1)求阶数：模拟频带边界频率p=tan(p2) s=tan(s2) 过度比 k=ps 分辨率k1=a2-1滤波器阶数n=lg(1k1)lg(1k)2)n阶低通巴特沃兹haj2=11+(c)2nhajp2=11+2hajs2=1a2脉冲响应不变法1), 将ha(s)进行逆拉氏变换，得到： 式中，是单位阶跃函数。2)对ha(t)进行等间隔采样，采样间隔为t，得到：3)对上式进行变换，得到数字滤波器的系统

23、函数，即优点：时域逼近良好，=t；（为模拟频率）缺点：容易产生混叠失真，只适用于带限滤波器；双线性不变法双线性变换法：优点：设计运算简单；避免了频谱的混叠效应，适合各种类型滤波器；缺点：=2ttan2，会产生非线性频率失真。eg:脉冲响应/双线性不变法 (1)图示是由rc组成的模拟滤波器1.写出传输函数 ；2.选用一种合适的转换方法将 转换成数字滤波器 ，设采样间隔为t；解：1.由图可知该滤波器为模拟高通滤波器，电阻s,电容1sc，电感sl2.采用双线性变换法，因为选用脉冲响应不变法，会在高频处发生频率混叠现象。(2)已知模拟滤波器传输函数为 ，设 ，用脉冲响应不变法和双线性变换法将 转换为数

24、字滤波器系统函数 。解：用脉冲响应不变法（令）将转换为数字滤波器系统函数。用双线性变换法将转换为数字滤波器系统函数。 (3)has=1s2+s+1 ,design a corresponding digital filter by bilinear transform respectively(assume t=2)1.hz=ha(s)s=2t1-z-11+z-1=1+z-123+z-2 2. 复数z=a+bi,模为a2+b2，幅角w=arctanba,然后用（a,b）的象限确定w的值带几个关键点，z=ej=0，2，32,算出z=a+bj，带入h(z)计算 hej () 3.2.2 fir滤波

25、器信号流图结构、线性相位系统、窗函数法、四种类型fir : hz=k=0nhkz-k 信号流图结构：直接型、级联型、线性相位型：直接型 级联型yn=h0xn+ h1xn-1+ h2xn-2+ h3xn-3+ h4xn-4 hz=h0k=1k(1+1kz-1+2kz-2)线性相位型(奇数7) 线性相位型(偶数8)h(z)=h0(1+z-6)+ h1x(z-1+z-5)+ h2(z-2+z-4)+ h3z-3 h(z)=h0(1+z-7)+ h1x(z-1+z-6)+ h2(z-2+z-5)+ h3(z-3+z-4)信线性相位eg:线性相位fir系统与其流图 (1)h1n=1,2,3,2,1and

26、 h2n=-2,1,0,1,21.determine the impulse response of overall system h1nh2n=-2,-5,-8,-6,0,6,8,5,22.show that oberall system is a generalized linear phase system,and plot the liner-phase structure flow graph.hej=h1ejh2ej =(1+2e-j+3e-j2+2e-j3+e-j4)(-2-e-j+e-j3+e-j4) =2e-j2cos2+2cos+3-2je-j2(2sin2+sin) =-4je-j4cos2+2cos+3(2sin2+sin)so it is a generalized linear system3.why cannot realize lowpass and highpass filter?when =0 or , hej=0(2)hz=2+z-11+2z-11