《数字信号处理》试题答案

《数字信号处理》期末考试A卷1一、填空题（本大题共7小题，每小题1分，共7分）1.序列x(n)=sin(0.3n+0.25)，该序列的周期N为20。2.序列x(n)存在傅里叶变换的充分条件是nnx)(。3.用DFT对序列进行谱分析时，对序列截断引起主谱线向附近展宽的现象称为频谱泄露。4.全通滤波器的极点和零点是互为共轭倒易关系。5.对12点长序列x(n)做DIF-基2FFT计算，其运算流图中每级的蝶形个数是8个。6.设计IIR滤波器的脉冲响应不变法，不适合设计高通、带阻滤波器。7.用频率采样法设计FIR数字滤波器，为了提高阻带衰减，可在频响间断点处内插一个或几个采样点。二、判断改错题，正确打“”，错误打“”，并改错。（本大题共4小题，每小题2分，共8分）8.周期序列的傅里叶级数仍是周期离散的。（）9.DIT-基2FFT分解的基本方法是将序列x(n)按n值前后对半分为2个序列。（）的奇偶10.序列x(n)的N点DFT为X(k)，则序列x*(n)的DFT变换为X*(Nk)。（），且X(N)=X(0)11.因果稳定的LTI时域离散系统，其系统函数所有零点都必须在单位圆内。（）极点三、计算题（本大题共6小题，共42分）12.已知序列()(1)2(3),()2()(2)xnnnhnnn，试计算循环卷积()()()ynxnhn，且循环卷积区间长度L=4。（6分）解：求x(n)和h(n)的DFT：304)()(nknWnxkXkkWW3442《数字信号处理》期末考试A卷2304)()(nknWnhkHkW24-2求X(k)与H(k)的乘积：)()()(kXkHkYkW24-2kkWW3442kkkkkkkkkW343求Y(k)的反变换得：)3(3nny13.若序列x(n)波形如下，且x(n)的FT变换为X(ej)，不直接求X(ej)，完成下列运算：求(1)X(ej)=？（2）2()?jXed（8分）解：（1）∵nnnnjjnxnjnnxenxeX)(]sin)[cos()()(∴3121111)()(42njnxeX（2）由帕斯维尔定理，有njnxdeX22)(2)(2341102111211nx(n)《数字信号处理》期末考试A卷3∴18)141111(2)(2)(2)(42222nnjnxnxdeX14.用微处理机对实序列作谱分析，要求谱分辨率F100Hz，信号最高频率为4kHz，试确定以下各参数：(1)最小记录时间TPmin；(2)最大取样间隔Tmax；(3)最少采样点数Nmin；(4)若信号频带不变，采用基2FFT做谱分析，求使谱分辨率提高1倍的N值。（8分）解：(1)已知最大谱分辨率F=100Hz，所以sFTP01.010011min(2)msffTs125.0400021211maxminmax(3)801025.101.04maxminminTTNP(4)频带不变，则取样率不变，分辨率提高1倍，则要求最小记录时间扩大1倍，此时有：1601025.102.04maxminminTTNP采用基2FFT做谱分析，由于128=27160256=28,所以N应取256。15.若某计算机实现复数乘法平均需要4s，实现复数加法平均需要1s，分别用直接计算和基2FFT计算N=1024点复数信号x(n)的N点DFTX(k)时，若只考虑乘法和加法运算，所需要的运算时间各是多少？（6分）解：直接计算所需运算量：复乘法次数N2次，复加法次数为N(N-1)次；基2FFT所需运算量：复乘法次数NN2log2次，复加法次数为NN2log次。因此，直接计算所需运算时间为：sNNNtD241856.510104755210419430410110231024104)1024(101)1(10466662662基2FFT计算所需运算量：《数字信号处理》期末考试A卷4sNNNNtF030720.01010240102048010110102410410512101log104log26666626216.已知某模拟系统22()276aHsss，采样间隔T=2s，试用脉冲响应不变法将()aHs转换为数字滤波器，求该数字滤波器系统函数H(z)。（8分）解：由因式分解可得：22232)23)(2(16722)(2sssssssHa即：23;221ss，因为T=2，所以有：421113141223122101134.90681.0063.0)0498.0)(0183.0()0315.0(20498.020183.020498.0140183.01414141414)(zzzzzzzzzzzzzezezezezH17.已知序列x(n)长度为N，X(k)=DFT[x(n)]，(),01()0,21xnnNynNnNY(k)为2N点y(n)的DFT值，试用X(k)表示Y(k)。（6分）解：由DFT定义有：整数2,)2()()()()(1021021202kkXWnxWnxWnykYNnnkNNnnkNNnnkN《数字信号处理》期末考试A卷5五、画结构图题（本大题共2小题，每小题8分，共16分）18.已知IIR系统为：51()(1)(2)()3(1)(2)66ynynynxnxnxn－－＋2，试：(1)确定系统的零、极点；（2）画出该系统的直接型网络结构图（II型）；（3）画出系统的级联型网络结构图。解：（1）由差分方程有：)2(2)1(3)()2(61)1(65)(nxnxnxnynyny方程两端求Z变换：)(2)(3)()(61)(65)(2121zXzzXzzXzYzzYzzY系统函数为：)31)(21()1)(2(61652361651231)()()(222121zzzzzzzzzzzzzXzYzH系统函数零点z1=2，z2=1；系统函数极点：31,2121pp（2）系统直接Ⅱ型结构流图如下图所示：系统函数可以写作：1111311121121)31)(21()1)(2()(zzzzzzzzzH由此可以画出系统的一种级联型结构流图如下图所示：x(n)y(n)-3z-15/6z-12-1/6x(n)y(n)z-1-21/31/2-1z-1《数字信号处理》期末考试A卷619.已知FIR滤波器系统函数为12456()10.81.51.50.8Hzzzzzz，求：(1)该滤波器的h(n)及阶数N；(2)判断系统是否具有线性相位，若是线性相位，求出相位函数和画出线性相位网络结构图；若不是线性相位，画出直接型网络结构图。解：（1）对H(z)求Z反变换可得：)6()5(8.0)4(5.1)2(5.1)1(8.0)()(nnnnnnnh滤波器阶数N=7（2）因为h(n)=-h(N-n-1)，而且h(n)为实数，所以该FIR滤波器是第二类线性相位。相位函数为：32)21(2)(N滤波器线性相位结构流图如下：六、证明题（本大题共1小题，每小题7分，共7分）20.若x(n)为实序列，且x(n)=x(Nn)，已知X(k)＝DFT[x(n)]N，证明：X(k)为实偶对称序列，即证明X(k)=X(Nk)且X(k)为实序列。解1：∵x(n)可写作x(n)=xep(n)+xop(n)，由x(n)=x(Nn)知，xop(n)=0，即x(n)=xep(n)；从而，由DFT的共轭对称性，DFT[xep(n)]=Re[X(k)]，即X(k)为实序列。由X(k)为实序列，即有X(k)=X(Nk)=X(Nk)。解2：因为：102)()]([NnnkjNenxnxDFT10]2sin2)[cos(NnnkNjnkNnxk=0，1，…，N-1由于x（n）是关于N的实偶序列，而nkN2sin是关于N的奇序列，所以有：02sin)(10NnnkNnx亦即：102cos)()]([NnnkNnxnxDFT为实序列；z-1z-1z-1x(n)-1-1-1z-1z-1z-10.8-1.5y(n)《数字信号处理》期末考试A卷7又有：)(2cos)(]2sin2sin2cos2)[cos()(2cos)()(101010kXnkNnxkNnkNnnxkNnNnxkNXNnNnNn七、设计题（本大题共2小题，每小题10分，共20分）21.用窗函数法设计一个FIR线性相位低通滤波器，要求阻带最小衰减为40dB，过渡带宽度为400Hz，采样频率为4kHz，确定设计要求滤波器的窗函数类型和长度；并构造逼近的频率响应函数()jdHe。附：常用窗函数的基本参数：窗函数类型旁瓣峰值an/dB过渡带宽度Bt阻带最小衰减as/dB矩形窗131.8/N21三角窗256.1/N25汉宁窗316.2/N44哈明窗416.6/N53布莱克曼窗5711/N74解：由题知2.040004002tB由阻带最小衰减为40dB和满足指标条件下N最小的原则，滤波器窗函数应选择汉宁窗，由附表可求得滤波器单位脉冲响应长度：312.02.6N设该FIR低通滤波器的通带截止频率为c，由滤波器阶数N=31知，该滤波器群延时为：1521N容易构造此滤波器所逼近的理想频率响应函数为：ccjjdeeH0)15（22.用双线性变换法设计IIR数字低通滤波器，要求通带频率低于0.25时，容许幅度差在2dB之内，频率在0.6到之间阻带衰减大于15dB，试采用巴特沃斯原型滤波器进行设计，采样间隔T=1s。附：已知归一化低通巴特沃斯滤波器的传输函数为：N＝1时，；N＝2时，；11aHpp2121aHppp《数字信号处理》期末考试A卷8N＝3时，巴特沃斯模拟滤波器阶数N计算公式为：；0.10.1101101spaspak；解：（1）将数字滤波器指标转换为相应的模拟滤波器指标（含预畸变）sradtgtgTpp/8284.0225.0222，dBp2sradtgtgTss/7528.226.0222,dBs15（2）计算模拟滤波器阶数N3230.38284.07528.2pssp2363.71101101.01.0psaaspk6481.15215.08595.03230.3lg2363.7lglglgspspkN取N=2（3）由给出的归一化原型巴特沃思低通滤波器传递函数可得所设计的模拟滤波器归一化传递函数为：（4）求出3dB频率：sradsNscs/1702.14251.07528.24251.0)110(211.0sradpNpcp/9472.01435.18284.01435.1)110(211.0由教材知（p158）：若选取c=1.1702rad/s，则阻带指标刚好满足要求，通带指标有富余，若选取c=0.9472rad/s，则通带指标刚好满足要求，阻带指标有富余。因为巴特沃斯曲线是单调的，为计算简单，取sradc/1