2018年西北工业大学电子信息学院816数字信号处理之数字信号处理考研强化五套模拟题
● 摘要
一、计算题
1． 对于一个稳定的离散时间线性时不变系统, 其输入
和输出
的关系为
确定其单位取样响应.
【答案】对差分方程两边取Z 变换, 得
则
其零点为边序列的和, 即
2． 利用窗函数（哈明窗）法设计一数字微分器, 逼近如图所示的理想微分器特性, 并绘出其幅频特性
.
极点为
根据系统是稳定的要求, 其收敛域要包含单位圆, 因此单位取样响应为一个左边序列和一个右
图
【答案】（1）由于连续信号存在微分, 而时域离散信号和数字信号的微分不存在, 因而本题要求设计的数字微分器是指用数字滤波器近似实现模拟微分器, 即用数字差分滤波器近似模拟微分
器. 下面先推导理想差分器的频率响应函数.
设模拟微分器的输入和输出分别为x （t ）和y （t ）, 即
令
则
对上式两边采样（时域离散化）, 得到
其中
. 将
和y （nT ）分别作为数字微分器的输入和输出序列, 并用
表示数
字理想微分器的频率响应函数, 则
即
根据图所给出的理想特性可知
所以应取k=T, 所以
取群延时
, 则逼近频率响应函数应为
设FIR 滤波器h （n ）长度为N , 一般取
. 加窗后得到
我们知道, 微分器的幅度响应随频率增大线性上升, 当频率上式中第一项为0, 所以
①式就是用窗函数法设计的FIR 数字微分器的单位脉冲响应的通用表达式, 且具有奇对称特性h （n ）=-h（N-l-n ）. 选定滤波器长度N 和窗函数类型, 就可以直接按①式得到设计结果. 当然,
时达到最大值, 所以只有N 为
,
偶数的情况4才能满足全频带微分器的时域和频域要求. 因为N 是偶数,
也可以用频率采样法和等波纹最佳逼近法设计.
本题要求的哈明窗函数:
将式代入
式得到h （n ）的表达式:
]
（2）对3种不同的长度N=20, 40和41, 用MATLAB 计算单位脉冲响应h （n ）和幅频特性函数, 并绘图的程序
如下:
程序运行结果即数字微分器的单位脉冲响应和幅频特性函数曲线如图所示. 由图可见, 当滤波器长度N 为偶数时, 逼近效果好. 但N=奇数时（本程序中N=41）, 逼近误差很大. （对第二类线性相位滤波器, N=奇数时不能实现高通滤波特性）
.