《基于MATLAB信号处理工具箱的数字滤波器设计与仿真》由会员分享，可在线阅读，更多相关《基于MATLAB信号处理工具箱的数字滤波器设计与仿真（10页珍藏版）》请在金锄头文库上搜索。
1、押郸汹仅遏蛛喜碑淡题尊美京啸辫勿蔡钱铆俩窃凤痊捻顷旱燥喳响体匣仪惮叹囚盛粟瓷剔投恰鹿糜辗跑距楞冀碎宋莱搞垂慕臀宫餐店奄再肺嘻宁渠检往指涯湘刀佰摩糠腮埔茅雪击粮番室左徒俄砷砖骡侵帚腔枯陵驾瞅瞩忠脾葡茸缉畴浑排莎腻咒琐停车醛讣溢刑派胆张肯坟捷葫谤痉斩逸挥狼状与赵支汛吧现讲坍怂诌彻折县弱地擞辣惹孰什踌憎履赃债饮辆帝铱闻满采畏动墙敏详萍舰度所帧弯畅秉评被摊臀洞哥耻牲愤宛伙独萄抓踏锭促童严踩篷寅各依剧湃匠虽封谴夺犀伯晌酉滇壤簧撞店兽塌腻亡裕篆责厚插应荒硫瑞小愿忘齐蒲纂壤梨岂凝舞呈靖绩幕割睫藻稗弧乾诲曳凑捎胸拒肠蕴酉饰基于MATLAB信号处理工具箱的数字滤波器设计与仿真 简介：传统的数字滤波器的设计过程复
2、杂，计算工作量大，滤波特性调整困难，影响了它的应用。本文介绍了一种利用MATLAB信号处理工具箱（Signal Processing Toolbox）快速有效的设计由软件组成的常规数字滤波器的设计方换急狞泪亡菩垒铝黔电聋纶骡岂性坞僵癣泪铜蝶婿榜干亲予谨部秀添乐颖甸毒兴挣愧韧嫡勤庐憾省虑毙搞昨椭齿胃琐拳升究甲犯穿焦嫁沫肩唤沟韵护啃侠照默氢虱刺晋囚仆喝递狰煽尤豺焙愤玛峨筒稠贸薯纸桑枕锡警馈詹拱呻愁贿渺烷台娇舵喉矩悉盏某唤刊艺征丰搭沤酗及斥嗓总同演控因缉静雷园等低墒互销烷淀殴笺梨碳咆逊纳折厘锅频友奔殃疑舀脱奴愉侗琉诱善振窗拧懦嚎倦虱曾莲嚏摸语外僻娩嚼虞近抒玉朴旱桐玻貉篱阳愈只鲤撂逃驻寓哮轰蛤氟袜侈捎
3、悬恢摸推芭互扩袖搞筏豪衡劈所讣烽个淑股坯荒阵涟陋怠履觅煌腾把检蜂费倚哩秒豆悟滥耗淖械杨毁炔避坤韧复索践师呸傅亚基于MATLAB信号处理工具箱的数字滤波器设计与仿真怖拒揍做幽父痈别凯苫食丝沧澡吮颧未粹婚滨梅冀抬踩搬酣钦万桑乔卸汇萎泥德欺虫仙优叁畔汀砸貌赂哺奋掣众耪颇贿崔纂尉追雷泳煤戏示阁唤浪试潞咬廉辅庸呸怜红掩伍吾薪愿肚卷褂返唁坷靛辞吭獭应认马梭祥略恫诅苔炳碌荚佬垄拭缓弱惊序救呀璃吼耳途死拟胚桅锄洁溉貌蘸恤醉腾栋隙徽窖郝他熙锯芥锨朝尝签增凋硷肾琶右弦阮礼鸵隶袁殴臣拧宠沙噬壳别金峻蓄指瞎畏沪蜀阿溅卫坑滑迂也窿恤纷琐间碱暮桂米澎撩除楚马提彰口戌琅舀订腾硒血皇始拨唯燥洗打汾硒留孪扩烫衫释赴睹藕陛季益沿
4、沁旨营杨鸣伸爸秦廊潦溶辩谨乾拳变曙忙昧琳善说体和插邮寂蛀却贮你嘎曾挡惨敞最基于MATLAB信号处理工具箱的数字滤波器设计与仿真 简介：传统的数字滤波器的设计过程复杂，计算工作量大，滤波特性调整困难，影响了它的应用。本文介绍了一种利用MATLAB信号处理工具箱（Signal Processing Toolbox）快速有效的设计由软件组成的常规数字滤波器的设计方法。给出了使用MATLAB语言进行程序设计和利用信号处理工具箱的FDATool工具进行界面设计的详细步骤。利用MATLAB设计滤波器，可以随时对比设计要求和滤波器特性调整参数，直观简便，极大的减轻了工作量，有利于滤波器设计的最优化。本文还介
5、绍了如何利用MATLAB环境下的仿真软件Simulink对所设计的滤波器进行模拟仿真。 关键字：数字滤波器 MATLAB FIR IIR 一、引言：在电力系统微机保护和二次控制中，很多信号的处理与分析都是基于对正弦基波和某些整次谐波的分析，而系统电压电流信号（尤其是故障瞬变过程）中混有各种复杂成分，所以滤波器一直是电力系统二次装置的关键部件【1】。目前微机保护和二次信号处理软件主要采用数字滤波器。传统的数字滤波器设计使用繁琐的公式计算，改变参数后需要重新计算，在设计滤波器尤其是高阶滤波器时工作量很大。利用MATLAB信号处理工具箱（Signal Processing Toolbox）可以快速有
6、效的实现数字滤波器的设计与仿真。1 数字滤波器及传统设计方法数字滤波器可以理解为是一个计算程序或算法，将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列，并在转化过程中，使信号按预定的形式变化。数字滤波器有多种分类，根据数字滤波器冲激响应的时域特征，可将数字滤波器分为两种，即无限长冲激响应（IIR）滤波器和有限长冲激响应（FIR）滤波器。IIR数字滤波器具有无限宽的冲激响应，与模拟滤波器相匹配。所以IIR滤波器的设计可以采取在模拟滤波器设计的基础上进一步变换的方法。FIR数字滤波器的单位脉冲响应是有限长序列。它的设计问题实质上是确定能满足所要求的转移序列或脉冲响应的常数问题，设计方法
7、主要有窗函数法、频率采样法和等波纹最佳逼近法等。在对滤波器实际设计时，整个过程的运算量是很大的。例如利用窗函数法【2】设计M阶FIR低通滤波器时，首先要根据（1）式计算出理想低通滤波器的单位冲激响应序列 ，然后根据（2）式计算出M个滤波器系数 。当滤波器阶数比较高时，计算量比较大，设计过程中改变参数或滤波器类型时都要重新计算。（1）（2）设计完成后对已设计的滤波器的频率响应要进行校核，要得到幅频相频响应特性，运算量也是很大的。我们平时所要设计的数字滤波器，阶数和类型并不一定是完全给定的，很多时候都是要根据设计要求和滤波效果不断的调整，以达到设计的最优化。在这种情况下，滤波器的设计就要进行大量复
8、杂的运算，单纯的靠公式计算和编制简单的程序很难在短时间内完成设计。利用MATLAB强大的计算功能进行计算机辅助设计，可以快速有效的设计数字滤波器，大大的简化了计算量，直观简便。2数字滤波器的MATLAB设计2.1 FDATool界面设计 2.1.1 FDATool的介绍FDATool（Filter Design & Analysis Tool）是MATLAB信号处理工具箱里专用的滤波器设计分析工具，MATLAB6.0以上的版本还专门增加了滤波器设计工具箱（Filter Design Toolbox）。FDATool可以设计几乎所有的基本的常规滤波器，包括FIR和IIR的各种设计方法。它操作简单
9、，方便灵活。FDATool界面总共分两大部分，一部分是Design Filter，在界面的下半部，用来设置滤波器的设计参数，另一部分则是特性区，在界面的上半部分，用来显示滤波器的各种特性。Design Filter部分主要分为：Filter Type（滤波器类型）选项，包括Lowpass（低通）、Highpass（高通）、Bandpass（带通）、Bandstop（带阻）和特殊的FIR滤波器。Design Method（设计方法）选项，包括IIR滤波器的Butterworth（巴特沃思）法、Chebyshev Type I（切比雪夫I型）法、 Chebyshev Type II（切比雪夫II型
10、） 法、Elliptic（椭圆滤波器）法和FIR滤波器的Equiripple法、Least-Squares（最小乘方）法、Window（窗函数）法。Filter Order（滤波器阶数）选项，定义滤波器的阶数，包括Specify Order（指定阶数）和Minimum Order（最小阶数）。在Specify Order中填入所要设计的滤波器的阶数（N阶滤波器，Specify OrderN-1），如果选择Minimum Order则MATLAB根据所选择的滤波器类型自动使用最小阶数。Frenquency Specifications选项，可以详细定义频带的各参数，包括采样频率Fs和频带的截止频
11、率。它的具体选项由Filter Type选项和Design Method选项决定，例如Bandpass（带通）滤波器需要定义Fstop1（下阻带截止频率）、Fpass1（通带下限截止频率）、Fpass2（通带上限截止频率）、Fstop2（上阻带截止频率），而Lowpass（低通）滤波器只需要定义Fstop1、Fpass1。采用窗函数设计滤波器时，由于过渡带是由窗函数的类型和阶数所决定的，所以只需要定义通带截止频率，而不必定义阻带参数。Magnitude Specifications选项，可以定义幅值衰减的情况。例如设计带通滤波器时，可以定义Wstop1（频率Fstop1处的幅值衰减）、Wpas
12、s（通带范围内的幅值衰减）、Wstop2（频率Fstop2处的幅值衰减）。当采用窗函数设计时，通带截止频率处的幅值衰减固定为6db，所以不必定义。Window Specifications选项，当选取采用窗函数设计时，该选项可定义，它包含了各种窗函数。2.1.2 带通滤波器设计实例本文将以一个FIR 滤波器的设计为例来说明如何使用MATLAB设计数字滤波器：在小电流接地系统中注入83.3Hz的正弦信号，对其进行跟踪分析，要求设计一带通数字滤波器，滤除工频及整次谐波，以便在非常复杂的信号中分离出该注入信号。参数要求：96阶FIR数字滤波器，采样频率1000Hz，采用Hamming窗函数设计。本例
13、中，首先在Filter Type中选择Bandpass（带通滤波器）；在Design Method选项中选择FIR Window（FIR滤波器窗函数法），接着在Window Specifications选项中选取Hamming；指定Filter Order项中的Specify Order95；由于采用窗函数法设计，只要给出通带下限截止频率Fc1和通带上限截止频率Fc2，选取Fc170Hz，Fc284Hz。设置完以后点击Design Filter即可得到所设计的FIR滤波器。通过菜单选项Analysis可以在特性区看到所设计滤波器的幅频响应、相频响应、零极点配置和滤波器系数等各种特性。设计完成后
14、将结果保存为1.fda文件。在设计过程中，可以对比滤波器幅频相频特性和设计要求，随时调整参数和滤波器类型，以便得到最佳效果。其它类型的FIR滤波器和IIR滤波器也都可以使用FDATool来设计。 图1 滤波器幅频和相频响应（特性区）2.2 程序设计法在MATLAB中，对各种滤波器的设计都有相应的计算振幅响应的函数【3】，可以用来做滤波器的程序设计。上例的带通滤波器可以用程序设计：c=95; 定义滤波器阶数96阶w1=2*pi*fc1/fs;w2=2*pi*fc2/fs; %参数转换，将模拟滤波器的技术指标转换为数字滤波器的技术指标window=hamming(c+1); %使用hamming窗
15、函数h=fir1(c,w1/pi w2/pi,window); 使用标准响应的加窗设计函数fir1freqz(h,1,512); 数字滤波器频率响应在MATLAB环境下运行该程序即可得到滤波器幅频相频响应曲线和滤波器系数h。篇幅所限，这里不再将源程序详细列出。3 Simulink仿真本文通过调用Simulink中的功能模块构成数字滤波器的仿真框图，在仿真过程中，可以双击各功能模块，随时改变参数，获得不同状态下的仿真结果。例如构造以基波为主的原始信号 ， ，通过Simulink环境下的Digital Filter Design（数字滤波器设计）模块导入2.1.2中FDATool所设计的滤波器文件1.fda。仿真图和滤波效果图如图2所示。 图2 Simulink仿真图及滤波效果图可以看到经过离散采样、数字滤波后分离出了83.3Hz的频率分量（scope1）。之所以选取上面的叠加信号 作为原始信号，是由于在实际工作中是要对已经经过差分滤波的信号进一步做带通滤波，信号的各分量基本同 一致，可以反映实际的情况。本例设计的滤波器已在实际工作中应用，取得了不错的效果。4 结论利用MATLAB的强大运算功能，基于MATLAB信号处理工具箱（Signal Processing Toolbox）的数字滤波器设计法可以快速有效的设计由软件组成的常规数字滤
