数字信号处理 巴特沃斯滤波 南京工程学院 你懂得

1、 数字信号处理课程设计 第页目录目录1课程设计任务书2语音信号滤波去噪使用脉冲响应不变法设计的巴特沃斯滤波器31、引言32、设计原理43、设计步骤7出现的问题及解决方法14结束语15参考文献15课程设计任务书 课程名称数字信号处理课程设计时间20102011学年第1学期14周学生姓名李晓伟指导老师金晓华题 目语音信号滤波去噪使用双线性变换法设计的切比雪夫II型滤波器主要内容： 用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用双线性变换法设计的一个满足指标的切比雪夫II型IIR滤波器，（或用脉冲响应不变法设计的一个满足指标的巴特沃斯IIR滤波器）对该语音信号进行滤波去噪处理，

2、比较滤波前后的波形和频谱并进行分析，根据结果和学过的理论得出合理的结论。要求：（1）滤波器指标必须符合工程实际。（2）设计完后应检查其频率响应曲线是否满足指标。（3）处理结果和分析结论应该一致，而且应符合理论。（4）独立完成课程设计并按要求编写课程设计报告书。应当提交的文件：（1）课程设计论文。（2）课程设计附件（主要是源程序）。语音信号滤波去噪使用脉冲响应不变法设计的巴特沃斯滤波器 摘要 本课程设计主要运用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用脉冲响应不变法设计的一个满足指标的巴特沃斯IIR滤波器，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析，

3、根据结果和学过的理论得出合理的结论。 关键词 课程设计；滤波去噪；巴特沃斯滤波器；脉冲响应不变法；MATLAB1 引 言本课程设计主要利用麦克风采集一段8000Hz，8k的单声道语音信号，并绘制波形观察其频谱，再用MATLAB利用脉冲响应不变法设计一个巴特沃斯滤波器，将该语音信号进行滤波去噪处理。1.1 课程设计目的数字信号处理课程设计是在学生完成数字信号处理和MATLAB的结合后的基本实验以后开设的。本课程设计的目的是为了让学生综合数字信号处理和MATLAB并实现一个较为完整的小型滤波系统。这一点与验证性的基本实验有本质性的区别。开设课程设计环节的主要目的是通过系统设计、软件仿真、程序安排与

4、调试、写实习报告等步骤，使学生初步掌握工程设计的具体步骤和方法，提高分析问题和解决问题的能力，提高实际应用水平。1.2 课程设计的要求（1）滤波器指标必须符合工程设计。（2）设计完后应检查其频率响应曲线是否满足指标。（3）处理结果和分析结论应该一致，而且应符合理论。（4）独立完成课程设计并按要求编写课程设计报告。2 设计原理用麦克风采集一段语音信号，绘制波形并观察其频谱，给定相应技术指标，用脉冲响应不变法设计的一个满足指标的巴特沃斯IIR滤波器，对该语音信号进行滤波去噪处理，比较滤波前后的波形和频谱并进行分析。2.1 IIR滤波器IIR滤波器设计方法有间接法和直接法，间接法是借助于模拟滤波器的

5、设计方法进行的。其设计步骤是：先设计过渡模拟滤波器得到系统函数H（s），然后将H（s）按某种方法转换成数字滤波器的系统函数H(z)。FIR滤波器比鞥采用间接法，常用的方法有窗函数法、频率采样发和切比雪夫等波纹逼近法。对于线性相位滤波器，经常采用FIR滤波器。 对于数字高通、带通滤波器的设计，通用方法为双线性变换法。可以借助于模拟滤波器的频率转换设计一个所需类型的过渡模拟滤波器，再经过双线性变换将其转换策划那个所需的数字滤波器。具体设计步骤如下：（1）确定所需类型数字滤波器的技术指标。（2）将所需类型数字滤波器的边界频率转换成相应的模拟滤波器的边界频率，转换公式为=2/T tan(0.5)(3)

6、将相应类型的模拟滤波器技术指标转换成模拟低通滤波器技术指标。（4）设计模拟低通滤波器。（5）通过频率变换将模拟低通转换成相应类型的过渡模拟滤波器。（6）采用双线性变换法将相应类型的过渡模拟滤波器转换成所需类型的数字滤波器。我们知道，脉冲响应不变法的主要缺点是会产生频谱混叠现象，使数字滤波器的频响偏离模拟滤波器的频响特性。为了克服之一缺点，可以采用双线性变换法。下面我们总结一下利用模拟滤波器设计IIR数字低通滤波器的步骤：（1）确定数字低通滤波器的技术指标：通带边界频率、通带最大衰减，阻带截止频率、阻带最小衰减。（2）将数字低通滤波器的技术指标转换成相应的模拟低通滤波器的技术指标。（3）按照模拟

7、低通滤波器的技术指标设计及过渡模拟低通滤波器。（4）用双线性变换法，模拟滤波器系统函数转换成数字低通滤波器系统函数。设IIR滤波器的输入序列为x(n)，则IIR滤波器的输入序列x(n)和输出序列y(n)之间的关系可用下列方程式表示：y(n)=i=0Mbixn-i+j=1Najy(n-j) (2-1)其中bi和aj是滤波器的系数，且aj中至少有一个非零。与上述差分方程对应的传递函数为：Hz=Y(z)X(z)=b0+b1z-1+bMz-M1-a1z-1-aNz-N (2-2)由传递函数可以发现无限长单位冲激响应滤波器有如下特点：（1） 单位冲激响应h(n)是无限长的。（2）系统传递函数H(z)在有

8、限z平面上有极点存在。（3）结构上存在着输出到输入的反馈，也就是结构上是递归型的。设Ha（j）是某个模拟滤波器的频率响应，那么低通滤波器在幅度平方响应上的技术指标给出为：11+2|Ha（j）|21，|p （2-3）0|Ha（j）|21A2，s| （2-4）式中是通带波纹参数，p是通带截止频率以rad/s（弧度每秒）计，A是阻带衰减参数，以及s是阻带截止频率以rad/s计。|Haj|2=11+2，=p （2-5）|Haj|2=1A2，=p （2-6）参数和A是分别与以dB计的参数Rp和As有关的，这些关系是Rp=-10log1011+2=10Rp10-1 （2-7）As=-10log101A2A

9、=10As20 （2-8）波纹1和2的绝对标尺是通过下式与和A有关的：1-11+1=11+2=211-1 21+1=1AA=1+12 利用幅度平方响应给出的模拟滤波器要求（2-1）和（2-2）式不包含任何相位信息。现在，为了求s域的系统函数Ha（s），考虑Haj=Ha（s）|s=j 那么有|Haj|2=HajHa*j=HajHa-j=HasHa-s|s=j （2-12）或者有HasHa-s=|Haj|2|=sj （2-9）因此，幅度平方函数的零点和极点相对于j轴是以镜像对称方式分布的【2】。接下来，我们根据语音信号的特点给出有关滤波器的技术指标：低通滤波器的性能指标：fp=1200Hz，fs=

10、1200Hz，ft=8000Hz, As=20db ,Ap=1dB2.2 巴特沃斯滤波器MATLAB信号处理工具箱函数buttp buttor butter是巴特沃斯滤波器设计函数，其有5种调用格式，本课程设计中用到的是N,wc=butter(N,wc,Rp,As,s),该格式用于计算巴特沃斯模拟滤波器的阶数N和3dB截止频率wc。函数butter,cheby1和ellip设计IIR滤波器时都是默认的双线性变换法，所以在设计滤波器时只需要代入相应的实现函数即可。下面我们设计的IIR数字滤波器的主要程序：wp=2*pi*Fp/Ft;ws=2*pi*Fs/Ft;fp=2*Ft*tan(wp/2);

11、fs=2*Ft*tan(ws/2);n,wn=buttord(wp,ws,Ap,As,'s'); b,a=butter(n,wn,'s'); num,den=bilinear(b,a,1); h,w=freqz(num,den); 3.设计步骤3.1设计流程图语音信号滤波去噪使用脉冲不变响应法设计的巴特沃斯滤波器的设计流程如图2所示：开始语音信号的采集（wavread函数），画时域图快速傅里叶变换，并且画频谱图设定滤波器性能指标，通带截止频率fb=1100，阻带截止频率fc=1200，通带波纹Ap=1，阻带波纹As=20脉冲响应不变法法设计巴特沃斯滤波器验证并进

12、行频谱分析设计好的滤波器进行滤波处理比较滤波前后语音信号的波形及频谱回放语音信号结束图1 脉冲响应不变法巴特沃斯滤波器对语音信号去噪流程图3.2语言信号的采集 图3是基于PC机的语音信号采集过程:通过计算机录音系统录制1秒： 图2 语音采集过程采集程序：y,fs,nbits=wavread('2.wav');%载入波形N=length(y)%计算输入波形长度Y=fft(y,N);%进行快速傅里叶变换sound(y);%播放声音figure(4);plot(y);figure(5);plot(abs(Y); 图3 录音过程截图3.3语音信号的频谱分析图4.1 时域波形图4.2 频

13、域幅度谱3.4滤波器设计将数字滤波器的设计指标设为通带截止频率fb=1100HZ,阻带频率fc=1200HZ,通带波纹Ap=1dB，阻带波纹As=20dB,要求确定H(z)。设计程序如下：Fp=1200;Ft=8000;As=20;Ap=1;wp=2*pi*Fp/Ft;ws=2*pi*Fs/Ft;fp=2*Ft*tan(wp/2);fs=2*Ft*tan(ws/2);n,wn=buttord(wp,ws,Ap,As,'s'); b,a=butter(n,wn,'s'); num,den=bilinear(b,a,1); h,w=freqz(num,den); f

14、igure(3)plot(w*8000*0.5/pi,abs(h); Buttorworth波形如下图5：图5 IIR低通滤波器3.5 信号滤波过程y,fs,nbits=wavread('2.wav');%载入波形N=length(y)%计算输入波形长度Y=fft(y,N);%进行快速傅里叶变换sound(y);%播放声音figure(4);plot(y);figure(5);plot(abs(Y);Fp=1200;%阻带截止频率Fs=1100;%通带截止频率Ft=8000;%采集频率As=20;%通带波纹Ap=1Ap=1;%阻带波纹As=20wp=2*pi*Fp/Ft;ws=

15、2*pi*Fs/Ft;fp=2*Ft*tan(wp/2);fs=2*Ft*tan(ws/2);n,wn=buttord(wp,ws,Ap,As,'s'); %求低通滤波器的阶数和截止频率b,a=butter(n,wn,'s'); %求S域的频率响应的参数 num,den=bilinear(b,a,1); %双线性变换实现S域到Z域的变换 h,w=freqz(num,den); %根据参数求出频率响应figure(3)plot(w*8000*0.5/pi,abs(h);z=filter(num,den,y);sound(z);m=fft(z); %求滤波后的信号f

16、igure(1)subplot(2,2,3);plot(abs(m),'r');title('滤波后信号的频谱');grid;subplot(2,2,4);plot(z,'b');title('滤波后的信号波形');grid;subplot(2,2,2);plot(y,'b');title('滤波前信号的波形');grid;subplot(2,2,1);plot(abs(Y),'r');title('滤波前信号的频谱');grid;figure(2);p=angle(

17、m);q=angle(Y);subplot(2,1,1);plot(q,'b');title('滤波前相位');grid;subplot(2,1,2);plot(p,'b');title('滤波后相位');grid; 所得波形为图6.1 滤波前后的时域对比波形和滤波前后频域频谱图图6.2滤波前后相位比较图3.5 结果分析 由图6.1中滤波前后波形比较可看出，经过滤波后的波形比原波形的振幅有所减小，去除了很多由于噪声所产生的干扰；从滤波前后的频谱比较可以看出经过滤波后除了原本的声音外，中间由于噪声产生的频谱波形已经滤除；由图6.2滤

18、波前后相位比较图可看出由于经过滤波，相位变得稀疏；经过MATLAB仿真，听滤波前后的声音，可以听出有明显的滤波效果。因此利用脉冲响应不变法设计的巴特沃斯滤波器已经达到了设计的要求。 4.出现的问题及解决方法 在这次的课程设计中我们设计的滤波器对语音信号的滤噪处理或多或少会因为我们的疏忽或者实践能力不够还是会遇到问题的，当遇到问题时我们不但运用自己学到的知识还得到同学和老师的帮助将其解决，不断完善和修改。设计中出现的问题及解决方法如下：（1）设计的滤波器达不到要求，因为在采集语音信号时，录音机的音频格式应该是PCM 8000kHz，8位，单声道。若由非人发出的声音可能会被滤波器滤去，且声音文件最后在一秒左右。（2）在运用Matlab设计滤波器时，当编辑完前面两条程序时