太原理工大学数字信号处理实验三用窗函数法设计FIR滤波器.doc

实验四 用窗函数法设计FIR滤波器 一、实验目的 熟悉FIR滤波器设计的基本方法。 掌握用户窗函数设计FIR数字滤波器的原理及方法，熟悉相应的计算机高级语言编程。 熟悉线性相位FIR滤波器的幅频特性和相位特性。 了解各种不同窗函数对滤波器性能的响应。 二、实验原理和方法 FIR滤波器的设计 目前FIR滤波器的设计方法主要有三种：窗函数法、频率取样法和切比雪夫等波纹逼近的最优化设计方法。本实验中的窗函数法比较简单，可应用现成的窗函数公式，在技术指标要求不高的时候是比较灵活方便的。它是从时域出发，用一个窗函数截取一个理想的得到h(n)，以有限长序列h(n)近似理想的：如果从频域出发，用理想的在单位圆上等角度取样得到H（k），根据h(k)得到H(z)将逼近理想的这就是频率取样法。 窗函数设计法 首先要对滤波器提出性能指标。一般是给定一个理想的频率响应，使所设计的FIR滤波器的频率响应去逼近所要求的理想的滤波器的响应窗函数设计的任务在于寻找一个可实现（有限长单位脉冲响应）的传递函数。 三、实验内容及步骤 1. 用Hanning窗设计线性带通滤波器（N=15，N=45） Window=Hanning(16); h=fir1(15,[0.3,0.5],Window) freqz(h,1) h = Columns 1 through 8 -0.0028 -0.0048 0.0350 0.0701 -0.0436 -0.1972 -0.0753 0.2184 Columns 9 through 16 0.2184 -0.0753 -0.1972 -0.0436 0.0701 0.0350 -0.0048 -0.0028 Window=Hanning(46); h=fir1(45,[0.3,0.5],Window) Window=Hanning(46); h=fir1(45,[0.3,0.5],Window) freqz(h,1) h = Columns 1 through 8 -0.0001 -0.0001 0.0002 -0.0003 0.0005 0.0039 0.0021 -0.0085 Columns 9 through 16 -0.0112 0.0050 0.0161 0.0040 -0.0044 0.0055 -0.0074 -0.0459 Columns 17 through 24 -0.0220 0.0800 0.1025 -0.0471 -0.1740 -0.0586 0.1599 0.1599 Columns 25 through 32 -0.0586 -0.1740 -0.0471 0.1025 0.0800 -0.0220 -0.0459 -0.0074 Columns 33 through 40 0.0055 -0.0044 0.0040 0.0161 0.0050 -0.0112 -0.0085 0.0021 Columns 41 through 46 0.0039 0.0005 -0.0003 0.0002 -0.0001 -0.0001 2. 用Rectangle窗设计线性带通滤波器（N=15，N=45） Window=boxcar(16); b=fir1(15,[0.3,0.5],Window) freqz(b,1) b = Columns 1 through 8 -0.0518 -0.0233 0.0799 0.0976 -0.0433 -0.1555 -0.0514 0.1392 Columns 9 through 16 0.1392 -0.0514 -0.1555 -0.0433 0.0976 0.0799 -0.0233 -0.0518 Window=boxcar(46); b=fir1(45,[0.3,0.5],Window) freqz(b,1) b = Columns 1 through 8 -0.0210 -0.0044 0.0041 -0.0043 0.0051 0.0270 0.0112 -0.0345 Columns 9 through 16 -0.0368 0.0136 0.0378