欢迎来到天天文库
浏览记录
ID:14325022
大小:82.00 KB
页数:4页
时间:2018-07-27
《基于matlab的fir滤波器设计与滤波(含matlab代码)》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、基于MATLAB的FIR滤波器设计与滤波摘 要介绍了利用MATLAB信号处理工具箱进FIR滤波器设计的三种方法:程序设计法、FDATool设计法和SPTool设计法,给出了详细的设计步骤,并将设计的滤波器应用到一个混和正弦波信号,以验证滤波器的性能。关键词 MATLAB,数字滤波器,有限冲激响应,窗函数,仿真1 前言 数字滤波器是一种用来过滤时间离散信号的数字系统,通过对抽样数据进行数学处理来达到频域滤波的目的。根据其单位冲激响应函数的时域特性可分为两类:无限冲激响应(IIR)滤波器和有限冲激响应(FIR)滤波器。与II
2、R滤波器相比,FIR的实现是非递归的,总是稳定的;更重要的是,FIR滤波器在满足幅频响应要求的同时,可以获得严格的线性相位特性。因此,它在高保真的信号处理,如数字音频、图像处理、数据传输、生物医学等领域得到广泛应用。2 FIR滤波器的窗函数设计法 FIR滤波器的设计方法有许多种,如窗函数设计法、频率采样设计法和最优化设计法等。窗函数设计法的基本原理是用一定宽度窗函数截取无限脉冲响应序列获得有限长的脉冲响应序列,主要设计步骤为: (1) 通过傅里叶逆变换获得理想滤波器的单位脉冲响应hd(n)。 (2) 由性能指标确
3、定窗函数W(n)和窗口长度N。 (3) 求得实际滤波器的单位脉冲响应h(n),h(n)即为所设计FIR滤波器系数向量b(n)。 (4) 检验滤波器性能。 本文将针对一个含有5Hz、15Hz和30Hz的混和正弦波信号,设计一个FIR带通滤波器,给出利用MATLAB实现的三种方法:程序设计法、FDATool设计法和SPTool设计法。参数要求:采样频率fs=100Hz,通带下限截止频率fc1=10Hz,通带上限截止频率fc2=20Hz,过渡带宽6Hz,通阻带波动0.01,采用凯塞窗设计。3 程序设计法 M
4、ATLAB信号处理工具箱提供了各种窗函数、滤波器设计函数和滤波器实现函数。本文的带通滤波器设计及滤波程序如下: [n,Wn,beta,ftype]=kaiserord([7131723],[010],[0.010.010.01],100); %得出滤波器的阶数n=38,beta=3.4 w1=2*fc1/fs;w2=2*fc2/fs;%将模拟滤波器的技术指标转换为数字滤波器的技术指标 window=kaiser(n+1,beta);%使用kaiser窗函数 b=fir1(n,[w1w2],window);
5、使用标准频率响应的加窗设计函数fir1 freqz(b,1,512);%数字滤波器频率响应 t=(0:100)/Fs; s=sin(2*pi*t*5)+sin(2*pi*t*15)+sin(2*pi*t*30);%混和正弦波信号 sf=filter(b,1,s);%对信号s进行滤波 程序执行的结果如图1所示: (1)滤波器幅频特性和相频特性 (2)滤波前后的波形图1滤波器特性和滤波效果图4 FDATool设计法 FDATool(FilterDesign&AnalysisTool
6、)是MATLAB信号处理工具箱专用的滤波器设计分析工具,操作简单、灵活,可以采用多种方法设计FIR和IIR滤波器。在MATLAB命令窗口输入FDATool后回车就会弹出FDATool界面。4.1带通滤波器设计 已知滤波器的阶数n=38,beta=3.4。本例中,首先在FilterType中选择Bandpass;在DesignMethod选项中选择FIRWindow,接着在Window选项中选取Kaiser,Beta值为3.4;指定FilterOrder项中的Specifyorder为38;采样频率Fs=100Hz,截止频
7、率Fc1=10Hz,Fc2=20Hz。设置完以后点击窗口下方的DesignFilter,在窗口上方就会看到所设计滤波器的幅频响应,通过菜单选项Analysis还可以看到滤波器的相频响应、组延迟、脉冲响应、阶跃响应、零极点配置等。设计完成后将结果保存为kaiser15.fda文件。4.2 Simulink仿真 在Simulink环境下,将滤波器文件kaiser15.fda导入DigitalFilterDesign模块,输入信号为s(t)=sin(10πt)+sin(30πt)+sin(60πt),生成的仿真图和滤波效果如图
8、2所示。 (1)Simulink仿真图 (2)滤波前后的离散波形图2Simulink仿真图和滤波效果图5 SPTool设计法 SPTool是MATLAB信号处理工具箱中自带的交互式图形用户界面工具,它包含了信号处理工具箱中的大部分函数,可以方便快捷地完
此文档下载收益归作者所有