1、指导教师:摘要 1一、基本原理 21.1 IIR数字滤波器设计原理1.2 FIR数字滤波器设计原理 3二、滤波器详细设计 52.1 信号源2.2 低通滤波器 62.3 带通滤波器 92.3 高通滤波器 11总结 13参考文献 14附录 15摘要滤波器广泛应用于电子电路中,具有平滑信号、选频、去噪和使交流信号变为直流信号的功能。可分为低通滤波器、高通滤波器、带通滤波器和带阻滤波器,也可分为模拟滤波器和数字滤波器,数字滤波器分为IIR数字滤波器和FIR数字滤波器。设计滤波器之前先根据实际需要确定滤波器类型和指标,通过MATLAB分别设计低通、带通和高通三种滤波器,首先通过三种滤波器对时域信号滤波,
2、然后对信号频域变换后滤波,最后对信号进行时域搬移后滤波。通过MATLAB得出的频谱图对信号的三种形式通过三种滤波器后的结果进行分析,通过分析比较各滤波器的优点和缺点。关键词:MATLAB IIR滤波器 FIR滤波器 频域变换 时域搬移一、基本原理数字滤波器是一个离散系统,其系统函数一般可表示为的有理多项式形式,即当ai;i=1,2,N都为0时,上式描述的系统称为有限脉冲响应数字滤波器,简称FIR(Finite-Impulse Response)数字滤波器。当系数ai;i=1,2,N中至少有一个是非0时,上式描述的系统称为无限脉冲响应数字滤波器,简称IIR(Infinite-Impulse Re
3、sponse)数字滤波器。对于IIR数字滤波器,一般满足MN,这时将系统称为N阶IIR数字滤波器。对于FIR数字滤波器,系统函数中的有理多项式的最高次幂M就是其阶数。在设计数字滤波器时,首先根据具体的滤波成分,确定待设计数字滤波器的技术指标,再求解出满足设计指标的离散系统的系统函数H(z)。数字滤波器的技术指标通常用通带截频,通带最大衰减Ap,阻带截频,阻带最小衰减AS表示。1.1 IIR数字滤波器设计原理IIR滤波器的设计就是根据给定的数字滤波器技术指标,确定上式中滤波器的阶数N和系数ai,bj。在满足技术指标的条件下,滤波器的阶数应尽可能低。应为滤波器阶数越低,实现滤波器的成本就越低。由于
4、模拟滤波器的设计技术已经非常成熟,且可得闭合形式的解,因此在设计IIR滤波器时,一般是通过模拟滤波器来设计滤波器。设计方法是先将数字滤波器技术指标转换为对应的模拟滤波器技术指标,然后设计满足技术指标的模拟滤波器H(s),再将模拟滤波器H(s)转换为对应的数字滤波器H(z)。因此,在IIR数字滤波器设计中,模拟滤波器设计是基础,模拟滤波器到数字滤波器的转换是核心。而模拟滤波器的设计都是通过设计模拟低通滤波器来实现,比较常见的模拟低通滤波器有巴特沃思和切比雪夫等。将模拟滤波器变换为数字滤波器的主要方法有脉冲响应不变法和双线性变换法。图1-1 IIR数字滤波器的设计过程1.2 FIR数字滤波器设计原
5、理表1-1 四种线性相位FIR滤波器的特性类型IIIIIIIV阶数M偶数奇数Hk的对称性偶对称奇对称0.5M0.5A()关于=0的对称性=的对称性)的周期24A(0)任意)可适用的滤波器类型LP,HP,BP,BS等LP,BP微分器,Hilbert(希尔伯特)变换器变换器,HP窗函数法设计FIR滤波器的基本思想是在时域逼近理想滤波器的单位脉冲相应。首先根据待逼近的理想滤波器的频率响应,由IDTFT求出理想滤波器的单位脉冲响应,再将无限长的加窗截断得到有限长序列。为了获得线性相位FIR滤波器,在窗函数法设计FIR滤波器的过程中,需要将线性相位因子加入理想滤波器的频率响应,具体步骤如下:(1) 根据
6、所需设计的滤波器,确定线性相位滤波器的类型(I型,II型,III型,IV型)。(2) 确定理想滤波器的幅度函数(3) 确定理想滤波器的相位对I型和II型线性相位FIR滤波器,对III型和IV型线性相位FIR滤波器(4) 计算(5) 截断得其中wk是长度N=M+1的矩形窗。二、滤波器详细设计2.1 信号源产生一个连续信号,包含低频、中频、高频分量,对其进行采样,进行频谱分析分别设计三种高通、低通、带通滤波器对信号进行滤波处理,观察滤波后信号的频谱。合成信号频谱程序如下:t=(1:100)/2000;x=sin(2*pi*t*20)+sin(2*pi*t*200)+sin(2*pi*t*500);
7、plot(t,x);n=1:100;t=n/2000;X=fft(x,512);w=(0:255)/256*1000;plot(w,abs(X(1:256);xlabel(Hz);ylabel(频率响应幅度title(合成信号频谱图图2-1 合成信号波形图图2-2 合成信号频谱图2.2 低通滤波器根据数字滤波器的设计原理,首先将数字域的指标转化为模拟域的指标设计模拟低通滤波器,然后应用双线性变换法将模拟滤波器转化为数字滤波器将指标转换成归一化模拟低通滤波器的指标,通过归一化的模拟低通滤波器阶数N和3dB截止频率的计算,将模拟域频率变换成模拟低滤波器H(s),并用双线性变换法将H(s)转换成数字
8、低通滤波器H(z),由此得到低通数字滤波器。IIR低通滤波器的设计程序为:Ft=8000;Fp=1000;Fs=1200;As=100;Ap=1;wp=2*pi* Fp/Ft;ws=2*pi*Fs/Ft;fp=2*Fp*tan(wp/2);fs=2*Fs*tan(ws/2);n11,wn11=buttord(wp,ws,1,50,sb11,a11=butter(n11,wn11,num11,den11=bilinear(b11,a11,0.5);h,w=freqz(num11,den11);plot(w*8000*0.5/pi,abs(h);低通滤波器图2-3 IIR低通滤波器读入信号对信号傅
9、里叶变换进行频谱分析,利用filter函数滤波,对得到的信号傅里叶变换进行频谱分析,与滤波前进行比较。图2-4 滤波前信号波形图2-5滤波后信号波形图2-6滤波后信号的频谱图2-7滤波前信号的频谱2.3 带通滤波器IIR带通滤波器的设计程序为:Fp1=1200;Fp2=3000;Fs1=1000;Fs2=3200;wp1=tan(pi*Fp1/Ft);wp2=tan(pi*Fp2/Ft);ws1=tan(pi*Fs1/Ft);ws2=tan(pi*Fs2/Ft);w=wp1*wp2/ws2;bw=wp2-wp1;wp=1;ws=(wp1*wp2-w.2)/(bw*w);n12,wn12=but
10、tord(wp,ws,1, 50,b12,a12=butter(n12,wn12,num2,den2=lp2bp(b12,a12,sqrt(wp1*wp2),bw);num12,den12=bilinear(num2,den2,0.5);h,w=freqz(num12,den12);带通滤波器图2-8 IIR带通滤波器图2-9滤波后的信号波形图2-10滤波后信号的频谱读入信号对信号傅里叶变换进行频谱分析,利用filter函数滤波,对得到的信号傅里叶变换进行频谱分析,与滤波前进行比较.并在滤波前后用sound函数播放语音信号,观察滤波前后的变化。2.3 高通滤波器IIR高通滤波器的设计程序为:Fp=4000;Fs=3500;wp1=tan(pi*Fp/Ft);ws1=tan(pi*Fs/Ft);ws=wp1*wp/ws1;n13,wn13=cheb1ord(wp,ws,1,50,b13,a13=cheby1(n13,1,wn13,num,den=lp2hp(b13,a13,wn13);num13,den13=bilinear(num,den,0.5);h,w=freqz(num13,den13);plot(w*21000*0.5/pi,abs(h);高通滤波器图2-11 IIR高通滤波器图2-12滤波后的信号波形图2-13滤波后的信号频谱
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1