2012~2013学年第1学期数字信号处理期末考试试卷学号:________________;姓名:__________________;成绩:___________一、填空计算题(每空1分,共25分,其中最后5个是判断题,填写“√”或“×”)1.设x[n]是一个如图所示的有限长序列,X(ejω)为其傅里叶变换,X(k)为其6点离散傅里叶变换,则可求得X(ej0)=____0____,X(ejπ)=____0____,X(0)=____0____,X(3)=___0_____,2()jXed=__22()20xn___;2.复指数序列ej0.5n(-∞n∞)的傅里叶变换DTFT_2(0.52)rr、ej0.5nu[n]的z变换为_0.511,11jzez_、2[]jmnNxne(0mN)的N点DFT_[]()XkNkm_;3.单位脉冲响应为h[n]=δ[n]-δ[n-1]的系统是____非时变___(时变、非时变)__因果____(因果、非因果)、__稳定___(稳定、不稳定)、___线性___(线性、非线性)系统;1z1z1z1z2481632)(ny1z)(nx图1某LTI系统的横截型结构4.某LTI系统的横截型结构如图1所示,该系统的单位脉冲响应为_2[()(6)]nunun__,系统函数为__61164()12zHzz___,该系统__不是__(是否)线性相位系统;5.FIR滤波器的窗函数设计法中,阻带衰减取决于____窗种类___,加特定形状窗口条件下,过渡带宽度取决于____窗口宽度____;6.一个时间连续的实信号xc(t),带宽限制在5KHz以下,即对于|Ω|≥2π(5000),Xc(jΩ)=0,以每秒10000个样本的采样率对信号xc(t)进行采样,得到一个长度为N=1000的序列x[n]=xc(nT)。x[n]的N点DFT记作X[k]。若已知X[400]=1+j,则X[600]=1-j,k=400对应Xc(jΩ)的连续频率是Ωk=24000rad/s,在该连续频率处Xc(jΩk)=441010j;7.任何信号通过线性时不变的离散时间系统不可能产生比输入信号本身更多的频率分量(√)8.离散时间系统的极点全部在Z平面的单位圆内,则系统一定是稳定的(×)9.因果线性时不变系统的其单位冲激响应未必是正半轴序列(×)10.线性常系数差分方程无论初始状态为何,总是代表线性时不变系统(×)11.线性时不变离散时间系统存在系统函数,则频率响应必存在且连续(×)二、(12分)某LTI因果系统用下面差分方程描述:y(n)=0.9y(n-1)+x(n)+0.9x(n-1)(a)求系统函数H(z)及单位脉冲响应h(n);(b)写出系统频率响应函数H(ejω)的表达式,说明该系统为低通滤波器还是高通滤波器?(c)该系统是否存在因果稳定的逆系统?解:(a)两边取z变换,得11()0.9()()0.9X()YzzYzXzzz,11()10.9H()()10.9YzzzXzz,111.8H()110.9zzz,0.9z。由于系统是因果系统,所以11()()1.80.9(1)0.9()0.90.9(1)nnnhnnununun.(b)收敛域包含单位圆。所以10.9()10.9jjjeHee,零点在0.9z,极点在0.9z,所以是低通滤波器。(c)系统1110.9H()10.9zzz,极点在0.9z,选择收敛域0.9z,则存在因果稳定的逆系统1110.9H()10.9zzz。三、(15分)在图3所示系统中,输入连续信号cxt的频谱jXc是带限的,即N时,jXc=0。离散时间系统jeH=10c其它。(a)为了使cyt=cxt,采样周期T最大可以取多少?(b)要使整个系统等效为低通滤波器,确定T的取值范围?(c)若给定采样频率1/T=20KHz,整个系统等效为截止频率为3kHz的理想低通滤波器,确定c及N的取值范围。C/DD/CtxcTjeHtycTjYejXe图3解:(a)等效模拟低通滤波器,现要求全通,则数字频率的最高频率限制在c之内,即NNcT,所以cNT。(b)等效模拟低通滤波器,则数字频率的最高频率大于c无混叠,或者有混叠,混叠大于c,即NNcT,且2NcT,化简得2ccNNT。(c)31230000.32010ccT,N的取值使采样后不混叠,且高于0.3,或混叠部分在0.3以上。所以NNcT,2()0.3NTT,0.320.3NT,即600034000N。四、(15分)已知序列[]4[]3[1]2[2][3]xnnnnn,其6点离散傅立叶变换(DFT)用X[k]表示。(a)若序列y[n]的长度为6,其6点离散傅立叶变换为Y[k]46[]kWXk,求y[n];(b)求[]xn﹡[]xn;(c)求[]xn④[]xn;解:(a)依据Y[k]46[]kWXk,y[]n是[]xn循环右移4位的结果,即6y[]((4))4[4]3[5]2[][1]nxnnnnn(b)()()16()24(1)25(2)20(3)10(4)4(5)(6)xnxnnnnnnnn(c)[]xn④[]xn26()28(1)26(2)20(3)nnnn五、(10分)采用Kaiser窗函数法设计一个广义线性相位的数字低通滤波器,经验公式如下0.40.1102(8.7)500.5842(21)0.07886(21)21500.021AAAAAA82.285AM要求性能指标为:0.4p,0.6s,通带纹波10.005,阻带纹波20.001。确定该滤波器的参数β、最小阶次及延迟;解:根据窗函数设计法的对称性,应当设20.001,截止频率()/20.5cps,0.2sp,1020log60AdB,所以5.653,最小阶次37M,延迟(1)/218M。六、(8分)研究一个如图所示长度为N的有限长序列x[n],实线表示序列在0和N-1之间取值的包络,1[]xn是x[n]后面补N个零的长度为2N的有限长序列。x[n]的N点DFT用Xk表示,1[]xn的2N点DFT用kX1表示,能否用kX1表示得出Xk,说明理由。解:解法一:由于DFT,Xk是DTFT在区间[0,2]的等间隔采样,所以2()()|jkNXkXe,122()()|jkNXkXe,所以1()(2)XkXk。解法二:211112200()()()NNknknNNnnXkxnWxnW,021kN。当2km为偶数时,1121200(2m)()()()NNmnmnNNnnXxnWxnWXm,01mN。所以1()(2)XkXk,01kN。七(15分)考虑两个实值有限长序列h[n]和x[n],0≤n≤58,若线性卷积为y[n]=x[n]*h[n],该线性卷积可用DFT进行计算,即分别计算出H[k]、X[k],然后通过IDFT计算出y[n]=IDFT{X[k]H[k]}。试问:(a)计算H[k]、X[k]的最小点数是多少?(b)若有复数基2-FFT程序可供使用,如何构造一序列z[n],通过一次调用该程序,并经简单计算得到H[k]和X[k],写出实现步骤。解:(a)两个序列线性卷积的长度为121117NN,则长度为117的圆周卷积可以计算线性卷积,所以计算DFT的最小点数是117.(b)步骤1:序列()xn,h()n均补零至长度71282N,满足121115NNN,即()0,1,,58()059,60,,127xnnxnn,h()0,1,,58h()059,60,,127nnnn步骤2:由于()xn,h()n的71282N点DFT分别为()Xk,H()k,所以构造序列z(n)=x(n)+jh(n),计算其71282N点DFTZ()k,得Z(k)=X(k)+jH(k)。步骤3:*()Z()(()(N))/2epXkkZkZk*Z()(()(N))/2opH(k)kZkZk