语音数字滤波系统的设计

《语音数字滤波系统的设计》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、课程设计课程名称数字信号处理题目名称语音数字滤波系统的设计学生学院自动化专业班级电子（2）学号3110001576学生姓名W指导教师徐迎晖2013年6月6日目录—，实验g白勺2二实验环境2三实验过程五、实验总€5%录音文件放在E盘目录下%播放原始声音%取声频右声道长度%制造频率81OOHZ的噪声Q/o加噪声后的音频Q/o播放加噪声后的音频输出第一个图，原始声音x语音数字滤波系统的设计-、实验目的语音数字信号在实际的采集过程中必然会伴随着背景噪声，下面针对语音信号处理过程中高频噪声干扰信号应用低通消除噪声

2、，可使带有噪声的原始语言信号得到明显改善。二、实验环境MATLAB三、实验过程(1)先对一段录音文件进行取样，(2)然后人为加噪声(3)再根据指标设计滤波器滤掉噪声四、实验数据源程序及实验结果如下：[x，fs，nbit]=wavread(’e:sound.wav’)sound(x，fs)x=x(:，2);N=length(x);n=1:N;c=3*sin(2*pi*8100*n./fs);x1=x+c’；sound（x1，fs）;subplot（3,2,1）;stem（n，x）;titleC原始声音信

3、号TxIabelCiT)X=fft(x);magX=abs(X(1:(N+1)/2+1));%傅立叶链换k=0:1:(N+1)/2;w=2*pi/N*k;subplot(3,2,2〉;stem(w/pi，magX)第二个，原始声音的幅频ltitle（•原始声音信号幅频图》;xlabel（•以pi为单位的频率•）subplot(3,2,3)；stem(n,x1);Q/o第三个，加噪声的信号titleC加入噪声的声音信号TxIabelCn’)X1=fft(x1);magX1=abs(X1(1:(N+1)/2

4、+1));subplot(3,2,4);stem(w/pi,magX1);%第四个，加噪声幅频ItitleC加入噪声信号幅频图》;xlabelC以pi为单位的频率Wp=0.1*pi;Ws=0.2*pi;Rp=1;As=15;T=1；%用双线性变换设计巴特沃兹0megaP=(2/T)*tan(Wp/2);0megaS=(2/T)*tan(Ws/2);ep=sqrt(10A(Rp/10)-1);Ripple=sqrt(1/(1+ep*ep));Attn=1/(10A(As/20));[N，wn]=butto

5、rd(OmegaP，OmegaS，Rp，As，’s’);[cs，ds]=butter(N，wn，’s’)[b，a]=bHinear(cs，ds，T)y=filter(b，a，x1);Y1=fft(y);magY1=abs(Y1(1:(N+1)/2+1));%滤波后的声音信号k=0:1:(N+1)/2;w=2*pi/N*k;subplot(3,2,5);stem(n，y);title(’滤波后的声音信号》；subplot(3,2,6);stem(w/pj，magY1)title(’滤波后的幅频图％xlab

6、elC以pi为单位的频率’)sound(y，fs);原始声音信号012345678n加入哇声的声音信号x10151001235678910n1I00.10.20.30.40.50.6以pi为单位的頻率迪波后的翻图0.70.80.91波进后的声音值号x1001231IIIII1911111403]2010000.10.20.30.40.50.6以pi为单位的頻率0.70.80.91五、实验总结—开始，得出的X是双声道，对称的数组，跟噪声C无法相加，所以取x的单声道，再经过转置后解决了这个问题。然后关于双线

7、性变换法设计巴特沃兹滤波器也费了不少功夫。就是第六幅图片，滤波后的声音幅频特性不知是不是这样。感觉噪声所在的频率为所设频率的一半左右，导致所给的设计指标参数不太对应。最后得出的sound(y,fs)基本跟原始声音一致。