基于MATLAB+GUI滤波器设计.docx
《基于MATLAB+GUI滤波器设计.docx》由会员分享,可在线阅读,更多相关《基于MATLAB+GUI滤波器设计.docx(61页珍藏版)》请在冰豆网上搜索。
![基于MATLAB+GUI滤波器设计.docx](https://file1.bdocx.com/fileroot1/2022-10/13/46000701-7eb4-4341-b968-22ff32b9d642/46000701-7eb4-4341-b968-22ff32b9d6421.gif)
基于MATLAB+GUI滤波器设计
MATLAB课程设计报告
基于MATLABGUI的“滤波器设计软件”设计
摘要
面对庞杂繁多的原始信号,如何提取所需信号、抑制不需要的信号这就需要使用滤波器。
滤波器的作用主要是选择所需频带的信号内容而抑制不需要的其他频带的信号内容。
数字滤波器因其精度高、可靠性好、灵活性大等优点,在语音信号处理、信号频谱估计、信号去噪、无线通信中的数字变频以及图像处理等工程实际应用中都很广泛。
根据其冲击响应函数的时域特性可将数字滤波器分为IIR(有限长冲击响应)和FIR(无限长冲击响应)。
作为强大的计算软件,MATLAB提供了编写图形用户界面的功能。
所谓图形用户界面,简称为GUI,是由各种图形对象,如图形窗口菜单按钮、文本框等构建的用户界面。
MATALB可以创建图形用户界面GUI(GraphicalUserInterface),它是用户和计算机之间交流的工具。
MATLAB将所有GUl支持的用户控件都集成在这个环境中并提供界面外观、属性和行为响应方式的设置方法,随着版本的提高,这种能力还会不断加强。
而且具有强大的绘图功能,可以轻松的获得更高质量的曲线图。
关键词:
MATLABGUIIIR滤波器FIR滤波器
4.3基于MatlabGUI的数字滤波器设计实现16
5设计总结和心得33
附录
1设计任务
1设计目的
1.巩固所学习过的关于数字信号处理的滤波器的设计
2.学习使用MATLABGUI设计工具应用平面
2设计要求
要求:
“滤波器设计软件”设计
1、要有人机交互界面。
2、滤波器设计的类型、参数均有多个值可供选择。
3、当用户选择(或输入)滤波器的参数时,能够及时显示设计好的滤波器的频谱。
4、该软件可设计FIR或IIR型滤波器至少一种。
2MATLABGUI的简介
MATLAB是美国MathWorks公司出品的商业数学软件,用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境,主要包括MATLAB和Simulink两大部分。
MATLAB是矩阵实验室(MatrixLaboratory)的简称,和Mathematica、Maple并称为三大数学软件。
它在数学类科技应用软件中在数值计算方面首屈一指。
MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等,主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。
MATLAB的基本数据单位是矩阵,它的指令表达式与数学、工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完相同的事情简捷得多,并且mathwork也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。
在新的版本中也加入了对C,FORTRAN,C++,JAVA的支持。
可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用。
MATALB可以创建图形用户界面GUI(GraphicalUserInterface),它是用户和计算机之间交流的工具。
MATLAB将所有GUl支持的用户控件都集成在这个
环境中并提供界面外观、属性和行为响应方式的设置方法,随着版本的提高,这种能力还会不断加强。
而且具有强大的绘图功能,可以轻松的获得更高质量的曲线图。
3滤波器设计原理
3.1滤波器概述
随着信息时代和数字世界的到来,数字信号处理已成为当今一门极其重要的学科和技术领域。
数字信号处理在通信、语音、图像,自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。
在数字信号处理中,数字滤波器占有极其重要的地位。
现代数字滤波器可以用软件或设计专用的数字处理硬件两种方式来实现,用软件来实现数字滤波器优点是随着滤波器参数的改变,很容易改变滤波器的性能。
根据数字滤波器单脉冲响应的时域特性可将数字滤波器分为两种,即IIR(InfiniteImpulseResponse)无限长脉冲响应数字滤波器和FIR(FiniteImpulseResponse)有限长脉冲响应数字滤波器。
从功能上分类,可分为低通、高通、带通、带阻滤波器。
3.2IIR数字滤波器
3.2.1IIR数字滤波器设计原理
滤波器的设计质上是寻找一个既能物理实现,又能满足给定频率特性指标要求的系统传输函数。
IIR滤波器一般采用递归型的结构,系统的输入与输出服从N阶差分方程:
相应的传输函数为:
设计IIR数字滤波器就是要确定传输函数中的系数、或零极点增益、、A,使滤波器的频率特性满足给定的性能指标要求。
设计原理主要包括两个方面:
一是根据设计指标,先设计出相应的模拟滤波器再通过脉冲响应不变法或双线性变换法转换成对应的数字滤波器;二是选择一种优准则,如最小均方准则,再在
,先最误差此准则下求出滤波器传输函数的系数。
根据设计理论,在MATLAB环境下设计IIR数字滤波器主要有四种方法:
一是典型设计法;二是完全设计法;三是最优设计法;四是工具设计法。
由于完全设计法程序简单,我们在这里利用完全设计法设计滤波器。
3.2.2IIR滤波器设计思想
IIR滤波器设计思想是:
利用已有的模拟滤波器设计理论,首先根据设计指标设计一个合适的模拟滤波器,然后再通过脉冲响应不变法或双线性变换法,完成从模拟到数字的变换。
常用的模拟滤波器有巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Ellipse)滤波器、贝塞尔(Bessel)滤波器等,这些滤波器各有特点,供不同设计要求选用。
滤波器的模拟数字变换,通常是复变函数的映射变换,也必须满足一定的要求
由于数字滤波器传输函数只与频域的相对值有关,故在设计时可先将滤波器设计指标进行归一化处理,设采样频率为Fs,归一化频率的计算公式是:
利用完全设计法设计数字滤波器的步骤:
(1)将设计指标归一化处理。
(2)根据归一化频率,确定最小阶数N和频率参数Wn。
可供选用的阶数选择函数有:
buttord,cheblord,cheb2ord,ellipord等。
(3)运用最小阶数N设计模拟低通滤波器原型。
根据最小阶数直接设计模拟低通滤波器原型,用到的函数有:
butter,chebyl,cheby2,ellip和bessel。
如[B,A]=butter(N,Wn,'type')设计'type'型巴特沃斯(Butterworth)滤波器filter。
N为滤波器阶数,Wc为截止频率,type决定滤波器类型,type=high,设计高通IIR滤波器,ftype=stop,设计带阻IIR滤波器。
(4)再用freqz函数验证设计结果。
3.2.3IIR滤波器设计编程实现
例如选择设计IIR的Butterworth低通滤波器,其Fs=22050Hz,Fp1=3400Hz,
Fs1=5000Hz,Rp=2dB,Rs=20dB
程序和效果图(图2)如下:
Fs=22050;Fp1=3400;Fs1=5000;Rp=3;Rs=20;%设计指标
wp1=2*Fp1/Fs;ws1=2*Fs1/Fs;%求归一化频率
%确定butterworth的最小阶数N和频率参数Wn
[n,Wn]=buttord(wp1,ws1,Rp,Rs);
[B,A]=butter(N,Wn);%确定传递函数的分子、分母系数
[h,f]=freqz(b,a,Nn,Fs_value);%生成频率响应参数
plot(f,20*log(abs(h)))%画幅频响应图
plot(f,angle(h));%画相频响应图
%[N,Wn]=buttord(Wp,Ws,Rp,Rs)确定butterworth的N和Wn
%[N,Wn]=cheblord((Wp,Ws,Rp,Rs)确定Chebyshev滤波器的N和Wn
%[N,Wn]=cheb2ord(Wp,Ws,Rp,Rs)确定Chebyshev2滤波器的N和Wn
%[N,Wn]=ellipord(Wp,Ws,Rp,Rs)确定椭圆(Ellipse)滤波器的N和Wn
%[B,A]=butter(N,Wn,'type')设计'type'型巴特沃斯(Butterworth)滤波器filter.
%[B,A]=cheby1(N,R,Wn,'type')设计'type'型切比雪夫Ⅰ滤波器filter.
%[B,A]=cheby2(N,R,Wn,'type')设计'type'型切比雪夫Ⅱ滤波器filter.
%[B,A]=ellip(N,Rp,Rs,Wn,'type')设计'type'型椭圆filter.
图2Butterworth低通滤波器
3.3FIR数字滤波器
3.3.1FIR数字滤波器设计原理
根据数字滤波器冲激响应的时域特征,可将数字滤波器分为2种,即无限长冲激响应滤波器(IIRDF)和有限长冲激响应滤波器(FIRDF)。
FIRDF具有突出的优点:
系统总是稳定的、易于实现线性相位、允许设计多通带(或多阻带)滤波器。
因此FIRDF在数字信号处理中得到广泛的应用。
但与IIRDF相比,在满足同样的阻带衰减的情况下需要较高的阶数。
滤波器阶数越高将占用更多的DSP运算时间。
因此,对FIRDF的设计目标是在满足指标要求的情况下尽量减少滤波器的阶数。
数字滤波器可以理解为是一个计算程序或算法,将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列,并在转化过程中,使信号按预定的形式变化。
FIRDF的冲激响应h(k)是有限长的M阶FIRDF系统函数可表示为
滤波器的输出
它的设计问题实质上是确定能满足所要求的转移序列或脉冲响应的常数问题,设计方法主要有窗函数法、频率采样法和等波纹最佳逼近法等。
若要逼近的理想滤波器的频率响应为常用的有巴特沃思滤波器、切比雪夫型滤波器、椭圆滤波器和巴塞尔滤波器。
3.3.2FIR滤波器设计思想
在这里我们利用窗函数法设计FIR滤波器。
窗函数法的基本想法是选取某一种合适的理想频率选择性滤波器(这种滤波器总是有一个非因果,无限长的脉冲响应),然后将它的脉冲响应截断(或加窗)以得到一个线形相位和因果的FIR滤波器。
因此,这种方法的重点在于选择某种恰当的窗函数和一种合适的理想滤波器。
窗函数法又称傅立叶级数法,一般是先给定所要求的滤波器的频率响应Hd(ejw),要求设计一个FIR滤波器的频率响应H(ejw)来逼近Hd(ejw).设计是在时域进行的,首先由傅立叶变换导出hd(n),因此hd(n)一定是无限长的序列,而h(n)是有限长的,即要用有限长的h(n)来逼近无限长的hd(n),最有效的方法是截断hd(n),或者说用一个有限长的窗口函数w(n)来截取hd(n),即:
h(n)=hd(n)w(n),因而窗函数的形状及长度的选择就很关键了.在MATLAB中常用的窗函数有矩形窗,Hanning窗,Hamming窗,Blackman窗,Kaiser窗等,这些窗函数各有优缺点,各有利弊,默认值是Hamming窗。
利用完全设计法设计数字滤波器的步骤:
(1)将设计指标归一化处理。
(2)根据归一化频率,选择函数buttord确定最小阶数N和频率参数Wn。
(3)确定窗口值。
Windows指定窗函数类型,默认为Hamming窗;可选Han