MATLAB窗函数法实现FIR的高通,带通和低通滤波器的程序要点.docxVIP

发M火潺 MATLABI程设计报告 学 院：地球物理与石油资源学院 班 级： 测井（基）11001 姓 名： 大牛啊啊啊 学 号： 班内编号： 指导教师： 陈义群 完成日期： 2013 年6月3日 题目 FIR滤波器的窗函数设计法及性能比较 FIR滤波器简介 数字滤波器是一种用来过滤时间离散信号的数字系统， 通过对抽样数据进行数学处理来 达到频域滤波的目的。根据其单位冲激响应函数的时域特性可分为两类： 无限冲激响应（IIR） 滤波器和有限冲激响应（FIR）滤波器。与IIR滤波器相比，FIR滤波器的主要特点为： a.线性相位；b.非递归运算。 FIR滤波器的设计 FIR滤波器的设计方法主要有三种： a.窗函数设计法；b.频率抽样发；c.最小平法抽样法； 这里我主要讨论在 MATLAB环境下通过调用信号分析与处理工具箱的几类窗函数来设计滤 波器并分析与比较其性能。窗函数法设计 FIR滤波器的一般步骤如下： a.根据实际问题确定要设计的滤波器类型； b.根据给定的技术指标，确定期望滤波器的理想频率特性； c.求期望滤波器的单位脉冲响应； d.求数字滤波器的单位脉冲响应； e. 应用。 常用的窗函数有 (1) Hanning window 二 nw(n) = 二 n w(n) = [0.5 - 0.5 cos( )] RN (n) N -1 (2)Hamming window 2二 n w(n)=[0.54 —0.46 cos( )] RN (n) N -1 (3)Balckman window 2月n 4只n _ w(n) =[0.42 -0.5cos( ) 0.08cos( -)] RN (n) N -1 N -1 (9.2.2) (9.2.3) (9.2.4) (4)Kaiser window w(n)=[I0{ 1,1 —[2n/(N w(n)=[ I0{ 1,1 —[2n/(N -1)]2} ITTj ]RN(n) (9.2.5) 式中I0（x）是零阶Bessel函数，可定义为 oOI0(x) =1 oO I0(x) =1 八[ m =1 (x/2) m! m -]2 当x=0时与矩形窗一致；当x = 5.4414时与海明窗结果相同; 当x = 8.885时与布莱克曼窗结果相 同。 3.窗函数的选择标准 .较低的旁瓣幅度，尤其是第一旁瓣； .旁瓣幅度要下降得快，以利于增加阻带衰减; .主瓣宽度要窄，这样滤波器过渡带较窄。 .常用窗函数的参数 常用窗函数的几个参照效据 要有以下几种： hanning(N) hamming(N) blackman(N) kaiser(n+1,beta) kaiserord freqz filter6.实验具体步骤 要有以下几种： hanning(N) hamming(N) blackman(N) kaiser(n+1,beta) kaiserord freqz filter 6.实验具体步骤 hanning hamming blackman kaiser 窗函数的调用 窗函数的调用 窗函数的调用 窗函数的调用 计算kaiser 窗函数的相关参数 求取频率响应 对信号进行滤波的函数 窗类型 旁髓峰值(dB) 主解宽度 最小阻带波)t (UB) 矩形窗 -13 -21 三角窗 -25 加/N -25 汉宁窗 -31 8G 一 44 哈明窗 一 41 8G -53 布莱克曼窗 一 57 一74 5. FIR滤波器的MATLAB实现方式 在MATLAB信号分析与处理工具箱中提供了大量 FIR窗函数的设计函数，本次用到主 窗函数，给以相同本次实验分另 U通过调用 hanning ,hamming ,Blackman,kaiser 窗函数，给以相同 的技术参数，来设计低通，带通，高通滤波器，用上述窗函数的选择标准来比较各种窗 函数的优劣，并给以一个简谐波进行滤波处理，比较滤波前后的效果。达到综合比较的 效果。 二、源代码 .禾U用 hanning hamming blackman kaiser 窗，设计——个彳氐通 FIR function lowpassfilter clc; clear all; Fs=100;%采样频率 fp=20;%通带截止频率 fs=30;%阻带起始频率 wp=2*pi*fp/Fs;%将模拟通带截止频率转换为数字滤波器频率 ws=2*pi*fs/Fs;%将模拟阻带起始频率转换为数字滤波器频率 wn=(wp+ws)/2/pi;% 标准化的截止频率响应 Bt=ws-wp; N0=ceil(6.2*pi/Bt);% 滤波器长度 N=N0+mod(N0+1,2); window1=hanning(N);% 使用 hanning 窗函数 window2=hamming(N);%