2016年郑州大学产业技术研究院数字信号处理复试笔试最后押题五套卷
● 摘要
一、综合题
1. 利用窗函数(哈明窗)法设计一数字微分器,逼近如图所示的理想微分器特性,并绘出其幅频特性。
图
【答案】(1)由于连续信号存在微分,而时域离散信号和数字信号的微分不存在,因而本题要求设计的数字微分器是指用数字滤波器近似实现模拟微分器,即用数字差分滤波器近似模拟微分器。下面先推导理想差分器的频率 响应函数。
,即
设模拟微分器的输入和输出分别为x (t )和y (t )
令
则
,得到
对上式两边采样(时域离散化)
其中
将
和
分别作为数字微分器的输入和输出序列,并用
表示数字理
想微分器的频率响应函数,则
即
根据图所给出的理想特性可知
所以应取取群延时
所以
则逼近频率响应函数应为
设FIR
滤波器
长度为N ,一般取
加窗后得到
我们知道,微分器的幅度响应随频率增大线性上升,当频率式中第一项为0, 所以
时达到最大值,所以只有N 为偶
上
数的情况4才能满足全频带微分器的时域和频域要求。因为N 是偶数,
式就是用窗函数法设计的FIR 数字微分器的单位脉冲响应的通用表达式,
且具有奇对称特性
选定滤波器长度N 和窗函数类型,就可以直接按式得到设计结果。当然,
也可以用频率采样法和等波纹最佳逼近法设计。 本题要求的哈明窗函数:
将
式代入
式得到
的表达式:
(2)对3种不同的长度N=20, 40和41,用MA TLAB 计算单位脉冲响应h (n )和幅频特性函数,并绘图的程序
如下:
程序运行结果即数字微分器的单位脉冲响应和幅频特性函数曲线如图所示。由图可见,当滤波器,逼近误差很大。这一结论与教长度N 为偶数时,逼近效果好。但N=奇数时(本程序中N=41)
材给出的理论一致(对 第二类线性相位滤波器,N=奇数时不能实现高通滤波特性)。
图
也可以采用调用等波纹最佳逼近法设计函
数
设计
率响应特性,
而在区间
来设计FIR 数字微分器的方法
。阶FIR 数字微分器,返回的单位脉冲响应向
上逼近理想微分器的频
量hn 具有奇对称特性。在大多数工程实际中,仅要求在频率区间
上频率响应特性要求为零时,
调用参数
,其中B 为过渡带宽度(即无关区)波纹最佳逼近法设计函数
不能太靠近
上频率响应特性不作要求,或要求为零。对微分器设计,
在区间
B 也不能太小,否则设计可能失败。调用等
如下:
设计本题要求的FIR 数字微分器的程序
调用remez 函数设计FIR 微分器
设置微分器边界频率(关于归一化)
调用remez 函数设计HR 微分器
以下为绘图部分(省略) 请读者运行该程序,观察设计效果。
2. 完成下面各题: (1)设
(2)
(3)并求出
求求的离散傅里叶级
以及收敛域;
,并定性画出幅频特性曲线;
画出
的波形,
将x (n )以5为周期进行周期性延拓,形成周期序列
相关内容
相关标签