数字信号处理期末考试题及参考

数字信号处理期末考试题及参考数字信号处理作为电子信息、通信工程、自动化等专业的核心课程，其理论性与实践性均较强，对学生的综合分析和解决问题能力要求较高。本次期末考试旨在全面考察学生对课程核心知识点的掌握程度，以及运用这些知识分析和解决实际问题的能力。以下为本次期末考试的试题及参考解答，供同学们复习参考。一、填空题(每空3分，共24分)1.对于一个线性时不变（LTI）离散系统，其输入输出关系在时域可以用_______卷积_______运算来描述，在z域则可以用_______系统函数乘积_______来描述。2.序列x(n)的傅里叶变换X(e^jω)是ω的_______周期_______函数，其周期为_______2π_______。3.若序列x(n)是实序列且为偶对称，则其离散傅里叶变换X(k)的实部_______偶对称_______，虚部_______奇对称_______。4.用矩形窗设计FIR滤波器时，主要缺点是_______阻带衰减较小_______和_______过渡带较宽_______。为了改善阻带衰减特性，可以采用_______凯泽窗_______等其他窗函数。二、简答题(每题10分，共30分)1.简述离散傅里叶变换（DFT）的物理意义，并说明其与序列傅里叶变换（DTFT）的主要区别与联系。参考答案：离散傅里叶变换（DFT）的物理意义在于，它将有限长的时域离散序列变换到频域的离散点上，这些离散点对应于序列傅里叶变换（DTFT）在[0,2π)区间内均匀采样的结果。具体而言，对于长度为N的序列x(n)，其DFTX(k)(k=0,1,...,N-1)表示了x(n)的频谱在N个等间隔频率点ω_k=2πk/N(k=0,1,...,N-1)上的采样值。这使得数字信号的频谱分析得以在数字计算机上实现。区别：*DTFT是对任意（可无限长）离散序列进行的变换，其结果X(e^jω)是频率ω的连续周期函数，周期为2π。*DFT仅对有限长序列（或通过周期延拓将无限长序列转化为有限长周期序列）进行变换，其结果X(k)是频率的离散采样点，点数与序列长度（或延拓后的周期）N相同。联系：*DFT的结果X(k)是其对应的DTFTX(e^jω)在频率点ω_k=2πk/N(k=0,1,...,N-1)处的采样值。*对于有限长序列，DFT可以看作是对其DTFT进行等间隔采样而得到的离散频谱表示。2.什么是窗函数？在FIR滤波器设计中，窗函数法的基本原理是什么？参考答案：窗函数是一种在时域上具有特定形状，且在有限区间外取值为零的函数。在数字信号处理中，特别是FIR滤波器设计中，窗函数用于截取无限长的理想滤波器单位脉冲响应h_d(n)，以得到有限长的实际滤波器单位脉冲响应h(n)。窗函数法设计FIR滤波器的基本原理如下：1.根据设计指标（如通带截止频率、阻带截止频率、通带波纹、阻带衰减等），确定理想滤波器的频率响应H_d(e^jω)。2.对理想频率响应H_d(e^jω)进行傅里叶逆变换，得到理想的单位脉冲响应h_d(n)。通常，h_d(n)是无限长的，且可能是非因果的。3.选择合适的窗函数w(n)（如矩形窗、汉宁窗、汉明窗、布莱克曼窗、凯泽窗等）。4.将理想脉冲响应h_d(n)与窗函数w(n)在时域相乘，得到实际FIR滤波器的单位脉冲响应h(n)=h_d(n)*w(n)。这一步实现了对h_d(n)的截断，使其成为有限长。5.对h(n)进行傅里叶变换，得到实际滤波器的频率响应H(e^jω)，并检验其是否满足设计指标。若不满足，可调整窗函数类型或长度，重复上述过程。窗函数的类型和长度直接影响最终滤波器的性能，如过渡带宽度和阻带衰减。3.简述IIR滤波器和FIR滤波器的主要特点，并比较它们在设计方法和应用场景上的差异。参考答案：IIR滤波器（无限脉冲响应滤波器）的主要特点：*单位脉冲响应h(n)是无限长的。*结构上通常采用递归结构，即输出不仅与当前和过去的输入有关，还与过去的输出有关。*可以用较少的阶数实现较高的选择性，相位特性一般是非线性的（除了特定设计的线性相位IIR滤波器，但实现复杂）。*设计方法多基于模拟滤波器的原型（如巴特沃斯、切比雪夫、椭圆滤波器），通过脉冲响应不变法或双线性变换法等映射到数字域。*可能存在稳定性问题（极点需在单位圆内）。FIR滤波器（有限脉冲响应滤波器）的主要特点：*单位脉冲响应h(n)是有限长的。*结构上主要采用非递归结构，输出仅与当前和过去的输入有关，因此一定是稳定的。*可以容易地设计成线性相位特性（满足h(n)的偶对称或奇对称条件）。*为达到与IIR滤波器相同的选择性，通常需要更高的阶数。*设计方法主要有窗函数法、频率采样法、等波纹最佳逼近法等。*由于线性相位和稳定性，在数据通信、图像处理等对相位失真敏感的领域有广泛应用。应用场景差异：IIR滤波器适用于对相位要求不高，追求高选择性和低计算复杂度的场合；FIR滤波器适用于对线性相位要求严格，或对稳定性有绝对要求的场合，尽管可能需要更高的阶数和计算量。三、分析与计算题(共46分)1.(16分)已知某线性时不变离散系统的系统函数为H(z)=(1-0.5z^{-1})/(1-1.5z^{-1}+0.5z^{-2})。(1)求该系统的单位脉冲响应h(n)。(2)判断该系统是否稳定，并说明理由。(3)若输入x(n)=u(n)（单位阶跃序列），求系统的零状态响应y(n)。参考答案：(1)求单位脉冲响应h(n)：对H(z)进行部分分式展开。H(z)=(1-0.5z^{-1})/(1-1.5z^{-1}+0.5z^{-2})=(z-0.5)/(z^2-1.5z+0.5)分母因式分解：z^2-1.5z+0.5=(z-1)(z-0.5)故H(z)=(z-0.5)/[(z-1)(z-0.5)]=1/(z-1)，其中z=0.5是一个零点，与极点z=0.5抵消。化简后H(z)=1/(z-1)=z^{-1}/(1-z^{-1})查z变换表，1/(1-z^{-1})<->u(n)，根据时移性质，z^{-1}/(1-z^{-1})<->u(n-1)。因此，h(n)=u(n-1)。*(注：原H(z)存在零极点相消，简化了系统。)*(2)判断系统稳定性：系统函数H(z)的极点由分母决定。原分母为1-1.5z^{-1}+0.5z^{-2}，令其等于零：1-1.5z^{-1}+0.5z^{-2}=0→z^2-1.5z+0.5=0→解得z1=1,z2=0.5。但由于分子分母有(z-0.5)的公因式，系统的实际极点仅为z=1（二阶极点？不，此处是分子一次，分母二次，抵消后分母为一次，极点为z=1）。对于线性时不变离散系统，稳定的充要条件是系统函数H(z)的所有极点均位于单位圆内。此处，极点z=1位于单位圆上（非单位圆内）。因此，该系统为临界稳定系统，而非BIBO稳定系统。(3)求输入x(n)=u(n)时的零状态响应y(n)：方法一：时域卷积。y(n)=x(n)*h(n)=u(n)*u(n-1)。u(n)*u(n-1)=sum_{k=-infty}^{infty}u(k)u(n-1-k)。当n-1<0即n<1时，和为0。当n>=1时，k的取值范围是0<=k<=n-1，共有n项，每项为1。所以y(n)=nu(n-1)。方法二：z域相乘。Y(z)=X(z)H(z)。X(z)=Z[u(n)]=1/(1-z^{-1})，|z|>1。H(z)=1/(z-1)=z^{-1}/(1-z^{-1})，|z|>1。Y(z)=[1/(1-z^{-1})]*[z^{-1}/(1-z^{-1})]=z^{-1}/(1-z^{-1})^2。查z变换表，1/(1-z^{-1})^2<->(n+1)u(n)。根据时移性质，z^{-1}/(1-z^{-1})^2<->nu(n-1)。因此，y(n)=nu(n-1)。2.(15分)考虑长度为4的有限长序列x(n)=[1,2,3,4](n=0,1,2,3)。(1)利用DFT的定义计算X(k)=DFT[x(n)]，k=0,1,2,3。(2)若y(n)=x(n)*x(n)（4点圆周卷积），求y(n)。参考答案：(1)计算X(k)=DFT[x(n)]，N=4。X(k)=sum_{n=0}^{3}x(n)W_4^{kn},k=0,1,2,3，其中W_4=e^{-j2π/4}=-j。X(0)=x(0)W_4^0+x(1)W_4^0+x(2)W_4^0+x(3)W_4^0=1+2+3+4=10X(1)=x(0)W_4^0+x(1)W_4^1+x(2)W_4^2+x(3)W_4^3=1*1+2*(-j)+3*(-1)+4*(j)=(1-3)+(-2j+4j)=-2+2jX(2)=x(0)W_4^0+x(1)W_4^2+x(2)W_4^4+x(3)W_4^6=1*1+2*(-1)+3*(1)+4*(-1)（因为W_4^4=(W_4^4)=1,W_4^6=W_4^(4+2)=W_4^2=-1）=1-2+3-4=-2X(3)=x(0)W_4^0+x(1)W_4^3+x(2)W_4^6+x(3)W_4^9=1*1+2*(j)+3*(-1)+4*(-j)（W_4^3=j;W_4^6=W_4^2=-1;W_4^9=W_4^(8+1)=W_4^1=-j）=(1-3)+(2j-4j)=-2-2j所以X(k)=[10,-2+2j,-2,-2-2j](2)求y(n)=x(n)④x(n)（4点圆周卷积）。根据圆周卷积定理，时域圆周卷积的DFT等于频域DFT的乘积。即Y(k)=X(k)*X(k)，k=0,1,2,3。Y(0)=X(0)*X(0)=10*10=100Y(1)=X(1)*X(1)=(-2+2j)(-2+2j)=4-4j-4j+4j²=4-8j-4=-8jY(2)=X(2)*X(2)=(-2)*(-2)=4Y(3)=X(3)*X(3)=(-2-2j)(-2-2j)=4+4j+4j+4j²=4+8j-4=8j对Y(k)进行IDFT得到y(n)：y(n)=(1/4)sum_{k=0}^{3}Y(k)W_4^{-kn},n=0,1,2,3。y(0)=(1/4)[Y(0)W_4^0+Y(1)W_4^0+Y(2)W_4^0+Y(3)W_4^0]=(1/4)[100+(-8j)+4+8j]=(1/4)(104)=26y(1)=(1/4)[Y(0)W_4^0+Y(1)W_4^{-1}+Y(2)W_4^{-2}+Y(3)W_4^{-3}]=(1/4)[100*1+(-8j)(j)+4*(-1)+8j*(-j)]（W_4^{-1}=j,W_4^{-2}=-1,W_4^{-3}=-j）=(1/4)[100+(-8j²)-4+(-8j²)]=(1/4)[100+8-4+8]（因为j²=-1）=(1/4)(112)=28y(2)=(1/4)[Y(0)W_4^0+Y(1)W_4^{-2}+Y(2)W_4^{-4}+Y(3)W_4^{-6}]=(1/4)[100*1+(-8j)(-1)+4*(1)+8j*(-1)]（W_4^{-4}=1,W_4^{-6}=W_4^{-2}=-1）=(1/4)[100+8j+4-8j]=(1/4)(104)=26y(3)=(1/4)[Y(0)W_4^0+Y(1)W_4^{-3}+Y(2)W_4^{-6}+