数字信号处理第三西安科大出高西全丁玉美课后答案第精要.ppt

　　程序运行结果： 　　用矩形窗、 Hanning窗、 Hamming窗和Blackman窗设计的单位脉冲响应h(n)及其损耗函数曲线如题20*解图所示。 由图可见， 滤波器长度N固定时， 矩形窗设计的滤波器过渡带最窄， 阻带最小衰减也最小； blackman窗设计的滤波器过渡带最宽， 阻带最小衰减最大。  　　21*. 将要求改成设计线性相位高通FIR滤波器， 重作 题20。 　　解： 本题的设计程序除了在每个fir1函数的调用参数中加入滤波器类型参数“high”外， 与第20题的程序完全相同， 请读者修改并运行程序， 完成本题。 　　22*. 调用MATLAB工具箱函数remezord和remez设计线性相位低通FIR滤波器， 实现对模拟信号的采样序列x(n)的数字低通滤波处理。 指标要求： 采样频率为16 kHz； 通带截止频率为4.5 kHz， 通带最小衰减为1 dB； 阻带截止频率为6 kHz， 阻带最小衰减为75 dB。列出h(n)的序列数据， 并画出损耗函数曲线。 　　解： 本题设计程序ex722.m如下： 　　%ex722.m： 调用remezord和remez设计线性相位低通FIR滤波器 　　Fs=16000； f=［4500， 6000］； 　　　　%采样频率， 边界频率为模拟频率(Hz) 　　m=［1， 0］；  　　rp=1； rs=75； dat1=(10∧(rp/20)-1)/(10∧(rp/20)+1)； 　　 　　dat2=10∧(-rs/20)；  　　rip=［dat1， dat2］；  　　［M， fo， mo， w］=remezord(f， m， rip， Fs)； 　　 M=M+1； %边界频率为模拟频率(Hz)时必须 　　　　　　　　%加入采样频率Fs 　　hn=remez(M， fo， mo， w)  　　%以下为绘图部分（省略） 　　程序运行结果：  h（n） =［－0.0023 0.0026 0.0207 0.0131 　　　　　　 －0.0185 0.0032 0.0278 －0.0306　　　　　　 －0.0176 0.0705 －0.0402 －0.1075 0.2927 0.6227 0.2927 　－0.1075  －0.0402 0.0705 －0.0176 －0.0306 0.0278 0.0032 －0.0185 0.0131　　　　　　 0.0207 0.0026 －0.0023］　　单位脉冲响应h(n)及其损耗函数曲线如题22*解图所示。 题22*解图 　　23*. 调用MATLAB工具箱函数remezord和remez设计线性相位高通FIR滤波器， 实现对模拟信号的采样序列x(n)的数字高通滤波处理。 指标要求： 采样频率为16 kHz； 通带截止频率为5.5 kHz， 通带最小衰减为1dB； 过渡带宽度小于等于3.5 kHz， 阻带最小衰减为75 dB。 列出h(n)的序列数据， 并画出损耗函数曲线。 　　解： 滤波器的阻带截止频率fs=5500－3500=2000 Hz。 本题设计程序ex723.m如下： 　　%ex723.m： 调用remezord和remez设计线性相位高通FIR滤波器 　　Fs=16000； f=［2000， 5500］； %采样频率， 边界频率为模拟频率(Hz) 　　m=［0， 1］；  　　rp=1； rs=75； dat1=(10∧(rp/20)-1)/(10∧(rp/20)+1)； 　　 　　dat2=10∧(-rs/20)；  　　rip=［dat2， dat1］；  　　［M， fo， mo， w］=remezord(f， m， rip， Fs)； 　　　%边界频率为模拟频率(Hz)时必须加入采样频率Fs 　　hn=remez(M， fo， mo， w) 　　程序运行结果： 滤波器长度为N=M+1=11， 　　单位脉冲响应h(n)及其损耗函数曲线如题23*解图所示， 请读者运行程序查看h(n)的数据。 题2