数字信号处理巴特沃斯滤波南京工程学院你懂得.docx
- 文档编号:15392940
- 上传时间:2023-07-04
- 格式:DOCX
- 页数:13
- 大小:100.23KB
数字信号处理巴特沃斯滤波南京工程学院你懂得.docx
《数字信号处理巴特沃斯滤波南京工程学院你懂得.docx》由会员分享,可在线阅读,更多相关《数字信号处理巴特沃斯滤波南京工程学院你懂得.docx(13页珍藏版)》请在冰点文库上搜索。
数字信号处理巴特沃斯滤波南京工程学院你懂得
目录…………………………………………………………………………………1
课程设计任务书……………………………………………………………………2
语音信号滤波去噪——使用脉冲响应不变法设计的巴特沃斯滤波器…………3
1、引言………………………………………………………………………………3
2、设计原理…………………………………………………………………………4
3、设计步骤…………………………………………………………………………7
出现的问题及解决方法……………………………………………………………14
结束语………………………………………………………………………………15
参考文献……………………………………………………………………………15
课程设计任务书
课程名称
数字信号处理课程设计
时间
2010~2011学年第1学期14周
学生姓名
李晓伟
指导老师
金晓华
题目
语音信号滤波去噪——使用双线性变换法设计的切比雪夫II型滤波器
主要内容:
用麦克风采集一段语音信号,绘制波形并观察其频谱,给定相应技术指标,用双线性变换法设计的一个满足指标的切比雪夫II型IIR滤波器,(或用脉冲响应不变法设计的一个满足指标的巴特沃斯IIR滤波器)对该语音信号进行滤波去噪处理,比较滤波前后的波形和频谱并进行分析,根据结果和学过的理论得出合理的结论。
要求:
(1)滤波器指标必须符合工程实际。
(2)设计完后应检查其频率响应曲线是否满足指标。
(3)处理结果和分析结论应该一致,而且应符合理论。
(4)独立完成课程设计并按要求编写课程设计报告书。
应当提交的文件:
(1)课程设计论文。
(2)课程设计附件(主要是源程序)。
语音信号滤波去噪——使用脉冲响应不变法设计的巴特沃斯滤波器
摘要本课程设计主要运用麦克风采集一段语音信号,绘制波形并观察其频谱,给定相应技术指标,用脉冲响应不变法设计的一个满足指标的巴特沃斯IIR滤波器,对该语音信号进行滤波去噪处理,比较滤波前后的波形和频谱并进行分析,根据结果和学过的理论得出合理的结论。
关键词课程设计;滤波去噪;巴特沃斯滤波器;脉冲响应不变法;MATLAB
1引言
本课程设计主要利用麦克风采集一段8000Hz,8k的单声道语音信号,并绘制波形观察其频谱,再用MATLAB利用脉冲响应不变法设计一个巴特沃斯滤波器,将该语音信号进行滤波去噪处理。
1.1课程设计目的
《数字信号处理》课程设计是在学生完成数字信号处理和MATLAB的结合后的基本实验以后开设的。
本课程设计的目的是为了让学生综合数字信号处理和MATLAB并实现一个较为完整的小型滤波系统。
这一点与验证性的基本实验有本质性的区别。
开设课程设计环节的主要目的是通过系统设计、软件仿真、程序安排与调试、写实习报告等步骤,使学生初步掌握工程设计的具体步骤和方法,提高分析问题和解决问题的能力,提高实际应用水平。
1.2课程设计的要求
(1)滤波器指标必须符合工程设计。
(2)设计完后应检查其频率响应曲线是否满足指标。
(3)处理结果和分析结论应该一致,而且应符合理论。
(4)独立完成课程设计并按要求编写课程设计报告。
2设计原理
用麦克风采集一段语音信号,绘制波形并观察其频谱,给定相应技术指标,用脉冲响应不变法设计的一个满足指标的巴特沃斯IIR滤波器,对该语音信号进行滤波去噪处理,比较滤波前后的波形和频谱并进行分析。
2.1IIR滤波器
IIR滤波器设计方法有间接法和直接法,间接法是借助于模拟滤波器的设计方法进行的。
其设计步骤是:
先设计过渡模拟滤波器得到系统函数H(s),然后将H(s)按某种方法转换成数字滤波器的系统函数H(z)。
FIR滤波器比鞥采用间接法,常用的方法有窗函数法、频率采样发和切比雪夫等波纹逼近法。
对于线性相位滤波器,经常采用FIR滤波器。
对于数字高通、带通滤波器的设计,通用方法为双线性变换法。
可以借助于模拟滤波器的频率转换设计一个所需类型的过渡模拟滤波器,再经过双线性变换将其转换策划那个所需的数字滤波器。
具体设计步骤如下:
(1)确定所需类型数字滤波器的技术指标。
(2)将所需类型数字滤波器的边界频率转换成相应的模拟滤波器的边界频率,转换公式为Ω=2/Ttan(0.5ω)
(3)将相应类型的模拟滤波器技术指标转换成模拟低通滤波器技术指标。
(4)设计模拟低通滤波器。
(5)通过频率变换将模拟低通转换成相应类型的过渡模拟滤波器。
(6)采用双线性变换法将相应类型的过渡模拟滤波器转换成所需类型的数字滤波器。
我们知道,脉冲响应不变法的主要缺点是会产生频谱混叠现象,使数字滤波器的频响偏离模拟滤波器的频响特性。
为了克服之一缺点,可以采用双线性变换法。
下面我们总结一下利用模拟滤波器设计IIR数字低通滤波器的步骤:
(1)确定数字低通滤波器的技术指标:
通带边界频率、通带最大衰减,阻带截止频率、阻带最小衰减。
(2)将数字低通滤波器的技术指标转换成相应的模拟低通滤波器的技术指标。
(3)按照模拟低通滤波器的技术指标设计及过渡模拟低通滤波器。
(4)用双线性变换法,模拟滤波器系统函数转换成数字低通滤波器系统函数。
设IIR滤波器的输入序列为
,则IIR滤波器的输入序列
和输出序列
之间的关系可用下列方程式表示:
(2-1)
其中
和
是滤波器的系数,且
中至少有一个非零。
与上述差分方程对应的传递函数为:
(2-2)
由传递函数可以发现无限长单位冲激响应滤波器有如下特点:
(1)单位冲激响应
是无限长的。
(2)系统传递函数
在有限z平面上有极点存在。
(3)结构上存在着输出到输入的反馈,也就是结构上是递归型的。
设
是某个模拟滤波器的频率响应,那么低通滤波器在幅度平方响应上的技术指标给出为:
(2-3)
(2-4)
式中
是通带波纹参数,
是通带截止频率以rad/s(弧度每秒)计,A是阻带衰减参数,以及
是阻带截止频率以rad/s计。
(2-5)
(2-6)
参数
和A是分别与以dB计的参数
和
有关的,这些关系是
(2-7)
(2-8)
波纹
和
的绝对标尺是通过下式与
和A有关的:
利用幅度平方响应给出的模拟滤波器要求(2-1)和(2-2)式不包含任何相位信息。
现在,为了求s域的系统函数
,考虑
那么有
(2-12)
或者有
(2-9)
因此,幅度平方函数的零点和极点相对于
轴是以镜像对称方式分布的【2】。
接下来,我们根据语音信号的特点给出有关滤波器的技术指标:
低通滤波器的性能指标:
fp=1200Hz,fs=1200Hz,ft=8000Hz,As=20db,Ap=1dB
2.2巴特沃斯滤波器
MATLAB信号处理工具箱函数buttpbuttorbutter是巴特沃斯滤波器设计函数,其有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);
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所示:
图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频域幅度谱
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);
figure(3)
plot(w*8000*0.5/pi,abs(h));
Buttorworth波形如下图5:
图5IIR低通滤波器
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=1
Ap=1;%阻带波纹As=20
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');%求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);%求滤波后的信号
figure
(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(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滤波前后相位比较图可看出由于经过滤波,相位变得稀疏;经过MATLAB仿真,听滤波前后的声音,可以听出有明显的滤波效果。
因此利用脉冲响应不变法设计的巴特沃斯滤波器已经达到了设计的要求。
4.出现的问题及解决方法
在这次的课程设计中我们设计的滤波器对语音信号的滤噪处理或多或少会因为我们的疏忽或者实践能力不够还是会遇到问题的,当遇到问题时我们不但运用自己学到的知识还得到同学和老师的帮助将其解决,不断完善和修改。
设计中出现的问题及解决方法如下:
(1)设计的滤波器达不到要求,因为在采集语音信号时,录音机的音频格式应该是PCM8000kHz,8位,单声道。
若由非人发出的声音可能会被滤波器滤去,且声音文件最后在一秒左右。
(2)在运用Matlab设计滤波器时,当编辑完前面两条程序时无法放出声音,后来发现我们应当把采集的语音信号wav文件放到Matlab的work文件夹中。
(3)还要在滤波器性能曲线的wc处画一根竖线,这样更方便看出结果,其中wc处线的确定还需计算出wb/pi的值。
(4)所有的时间波形横坐标都要化为时间,滤波前后频谱的横坐标应是频率,这样在观察通带截止频率和阻带截止频率时更加精确,误差较小。
(5)要区分清楚ft、fs、fp所代表的意思。
5.结束语
本次经过一周的课程设计,让我们更熟悉了数字信号处理的基本知识和MATLAB的m语言,让我们把课上的理论知识运用到实际中去,让我们更近一步地巩固了课堂上所学的理论知识,并能很好地理解与掌握数字信号处理中的基本概念、基本原理、基本分析方法。
因为学过数字信号处理这门课,但这只是理论知识,通过实验我们才能真正理解其意义。
语音信号处理在以前的学习中并没有怎么接触,但是通过网上查找资料,我很快的理解了其基本原理。
经过这次的DSP课程设计,让我有机会将自己学到的理论知识运用到实际中,提高了自己的动手能力和思维能力。
在课程设计中发现自己的不足,所以在今后的学习和生活中我们要更加努力,学习好我们的专业知识并要能运用到实际。
参考文献
[1]程佩青数字信号处理清华大学出版社2007年2月
[2]从玉良数字信号处理原理及matlab实现电子工业出版社2005年7月
[3]刘波matlab信号处理电子工业出版社2006年1月
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字信号 处理 巴特沃斯 滤波 南京 工程学院 懂得
![提示](https://static.bingdoc.com/images/bang_tan.gif)