fir数字滤波器数字信号处理实验报告——用窗函数法设计fir数字滤波器iviweu.docVIP

PAGE fir数字滤波器 数字信号处理实验报告——用窗函数法设计FIR数字滤波器 实验四 用窗函数法设计 FIR数字滤波器 一、 实验目的 (1) 掌握用窗函数法设计FIR数字滤波器的原理和方法。 (2) 熟悉线性相位FIR数字滤波器特性。 (3) 了解各种窗函数对滤波特性的影响。 二、 实验原理、 滤波器的理想频率响应函数为Hd… 南京邮电大学 实验报告 系统频率响应和样本处理算法实现 加窗和离散傅氏变换 数字滤波器设计 课程名称：数字信号处理 姓 名： 学 号 开课时间 2011 /2012 学年， 第 2 学期 实验一 一．实验名称：采样、系统性质及滤波 二．实验目的和任务… 数字信号处理专业课程设计报告书 系部名称学生姓名专业名称班 级实习时间 电子与信息工程系 电子信息工程 2011年3月28日至2011年4月11日 ： ： ： ： ： 实验报告 题目一 一：实验题目 A sinusoid of frequency … 实验四 用窗函数法设计 FIR数字滤波器 一、 实验目的 (1) 掌握用窗函数法设计FIR数字滤波器的原理和方法。 (2) 熟悉线性相位FIR数字滤波器特性。 (3) 了解各种窗函数对滤波特性的影响。 二、 实验原理、 滤波器的理想频率响应函数为Hd(ejω)，则其对应的单位脉冲响应为 hd(n) = 12? ? ? ?? Hd(e j? )e j?n d? 窗函数设计法的基本原理是用有限长单位脉冲响应序列h(n)逼hd(n)。由于hd(n)往往是无限长序列，且是非因果的，所以用窗函数。w(n)将hd(n)截断，并进行加权处理： h(n) = hd(n) w(n) h(n)就作为实际设计的FIR数字滤波器的单位脉冲响应序列，其频率响应函数H(ejω)为 N?1 H(e) =? n?0 jω h(n)e ?j?n 用窗函数法设计的滤波器性能取决于窗函数w(n)的类型及窗口长度N的取值。设计过程中，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N。 三、 实验内容 1.MATALAB程序 任务一 N1=15; N2=33; b1=fir1(14,1/4,hanning(15)); b2=fir1(32,1/4,hanning(33)); [H1,W]=freqz(b1,1); H1_db=20*log10(abs(H1)); magH1=abs(H1); phaH1=angle(H1); [H2,W]=freqz(b2,1); H2_db=20*log10(abs(H2)); magH2=abs(H2); phaH2=angle(H2); figure(1); subplot(2,1,1);stem(b1); title(‘N=15时,汉宁窗得到的FIR滤波器脉冲响应’); subplot(2,1,2);stem(b2); title(‘N=33时,汉宁窗得到的FIR滤波器脉冲响应’); figure(2); subplot(3,1,1);plot(W/pi,magH1); title(‘N=15时汉宁窗得到的FIR滤波器的幅频特性’) subplot(3,1,2);plot(W/pi,H1_db); title(‘N=15时汉宁窗得到的FIR滤波器的对数幅频特性’) subplot(3,1,3);plot(W/pi,phaH1); title(‘N=15时汉宁窗得到的FIR滤波器的相频特性’) figure(3); subplot(3,1,1);plot(W/pi,magH2); title(‘N=33时汉宁窗得到的FIR滤波器的幅频特性’) subplot(3,1,2);plot(W/pi,H2_db); title(‘N=33时汉宁窗得到的FIR滤波器的对数幅频特性’) subplot(3,1,3);plot(W/pi,phaH2); title(‘N=33时汉宁窗得到的FIR滤波器的相频特性’) 任务二 N=33; b1=fir1(32,1/4,boxcar(33)); b2=fir1(32,1/4,hanning(33)); b3=fir1(32,1/4,bartlett(33)); b4=fir1(32,1/4,blackman(33)); [H1,W]=freqz(b1,1); H1_db=20*log10(abs(H1)); magH1=abs(H1); phaH1=angle(H1); [H2,W]=freqz(b2,1); H2_db=20*log10(abs(H2)); magH2=abs(H2); phaH2=angle(H2); [H3,W]=freqz(b3,1); H3_db=20*log10(abs