数字信号处理(第三版)课后习题全程序.ppt

　　(2) 对低通滤波器单位脉冲响应乘以cosπn可得到高通滤波器单位脉冲响应：  　　· 矩形窗： · 升余弦窗： · 改进升余弦窗： · 布莱克曼窗： 题12图 　　12． 利用窗函数（哈明窗）法设计一数字微分器， 逼近题12图所示的理想微分器特性， 并绘出其幅频特性。 　　解: (1) 由于连续信号存在微分， 而时域离散信号和数字信号的微分不存在， 因而本题要求设计的数字微分器是指用数字滤波器近似实现模拟微分器， 即用数字差分滤波器近似模拟微分器。 下面先推导理想差分器的频率响应函数。 　　设模拟微分器的输入和输出分别为x(t)和y(t)， 即 令x(t)=ejΩt， 则 　　　　　　　　y(t)=jkΩeΩt=jkΩx(t) 对上式两边采样（时域离散化）， 得到 其中ω=ΩT。 将x(nT)和y(nT)分别作为数字微分器的输入和 输出序列， 并用Hd(ejω)表示数字理想微分器的频率响应函数， 则 即 根据题12图所给出的理想特性可知 所以应取k=T， 所以 　　　　　　　　　Hd(ejω)=jω 取群延时τ=(N－1)/2， 则逼近频率响应函数应为 Hd(ejω)=jωe－jωτ=ωe－j(ωτ－π/2) 　　设FIR滤波器h(n)长度为N， 一般取τ=(N－1)/2。 加窗后得到 我们知道， 微分器的幅度响应随频率增大线性上升， 当频率ω=π时达到最大值， 所以只有N为偶数的情况4才能满足全频带微分器的时域和频域要求。 因为N是偶数， τ=N/2－1/2=正整数－1/2， 上式中第一项为0， 所以 ① 　　①式就是用窗函数法设计的FIR数字微分器的单位脉冲响应的通用表达式， 且具有奇对称特性h(n)= －h(N－1－n)。 选定滤波器长度N和窗函数类型， 就可以直接按①式得到设计结果。 当然， 也可以用频率采样法和等波纹最佳逼近法设计。 　本题要求的哈明窗函数： ② 将②式代入①式得到h(n)的表达式： ③ 　　(2) 对3种不同的长度N=20，40和41，用MATLAB计算单位脉冲响应h(n)和幅频特性函数，并绘图的程序ex712.m如下： 　　%ex712.m： 用哈明窗设计线性相位FIR微分器 　　clear all； close all；  　　N1=20； n=0： N1－1； tou=(N1－1)/2；  　　h1n=sin((n－tou)*pi)./(pi*(n-tou).∧2).*(hamming(N1))′； 　　N2=40； n=0： N2－1； tou=(N2－1)/2；  　　h2n=sin((n－tou)*pi)./(pi*(n-tou).∧2).*(hamming(N2))′； 　　N3=41； n=0： N3－1； tou=(N3－1)/2；  　　h3n=sin((n－tou)*pi)./(pi*(n－tou).∧2).*(hamming(N3))′; 　　h3n((N3－1)/2+1)=0； 　　　　 %因为该点分母为零， 无定义， 所以赋值0 　　%以下为绘图部分（省略） 　　程序运行结果即数字微分器的单位脉冲响应和幅频特性函数曲线如题12解图所示。 由图可见， 当滤波器长度N为偶数时， 逼近效果好。 但N=奇数时（本程序中N=41）， 逼近误差很大。 这一结论与教材给出的理论一致（对第二类线性相位滤波器， N=奇数时不能实现高通滤波特性）。 题12解图 　　也可以采用调用等波纹最佳逼近法设计函数remez来设计FIR数字微分器的方法。  hn=remez(N－1， f， m， ′defferentiator′) 设计N－1阶FIR数字微分器， 返回的单位脉冲响应向量hn具有奇对称特性。 在大多数工程实际中， 仅要求在频率区间0≤ω≤ωp上逼近理想微分器的频率响应特性， 而在区间ωpω≤π上频率响应特性不作要求， 或要求为零。 对微分器设计， 在区间ωpω≤π上频率响应特性要求为零时， 调用参数f=［0, ωp/π, (ωp+B)/π, 1］， m=［0,ωp/π, 0, 0］， 其中B为过渡带宽度（即无关区）， ωp不能太靠近π， B也不能太小， 否则设计可能失败。 调用等波纹最佳逼近法设计函数remez设计本题要求的FIR数字微分器的程序ex712b.m如下： 　　% ex712b.m： 调用remez函数设计FIR微分器　　Wp=0.9； B=0.09； 　　　　%设置微分器边界频率（关于π归一化）　　N=40； f=［0，wp，wp