基于MATLAB_的语音信号分析和处理(20220301211143).pdf

《基于MATLAB_的语音信号分析和处理(20220301211143).pdf》由会员分享，可在线阅读，更多相关《基于MATLAB_的语音信号分析和处理(20220301211143).pdf（25页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、1 基于MATLAB 的语音信号分析和处理福建师范大学协和学院信息技术系电子信息科学与技术专124122006028 王祯飞指导老师黄小芬【摘要】本设计用微软录音机采集了一段语音，对其进行了时域分析，频谱分析,分析语音信号的特性。并应用matlab 平台对语音信号加入了不同的噪声，进一步用双线性变换法设计了一个的巴特沃思低通滤波器，基于巴特沃斯模拟滤波器设计数字带通滤波器和用窗函数法设计了FIR 低通滤波器，然后对加噪的语音信号进行滤波处理。最后对比滤波前后的语音信号的时域和频域特性，回放加噪语音信号和去噪语音信号，对比研究了巴特沃思IIR 滤波器和FIR 滤波器在信号处理尤其是信号滤波去噪方

2、面的特性和应用。论文从理论和实践上比较不同数字滤波器的滤波效果。【关键词】语音信号；频域特性;时域特性;滤波器2 1 绪论1.1 课题的研究意义语言是我们人类所特有的功能，它是传承和记载人类几千年文明史，没有语言就没有我们今天人类的文明。语音是语言最基本的表现形式，是相互传递信息最重要的手段，是人类最重要、最有效、最常用和最方便的交换信息的形式。语音信号处理属于信息科学的一个重要分支，大规模集成技术的高度发展和计算机技术的飞速前进，推动了这一技术的发展；它是研究用数字信号处理技术对语音信号进行处理的一门新兴学科，同时又是综合性的多学科领域和涉及面很广的交叉学科，因此我们进行语言信号处理具有时代

3、的意义。1.2 设计任务本设计先完成语音信号的采集，然后设计低通，高通，带通等滤波器对采集到的语音信号进行滤波处理，分析语音信号各频率段的特性。并对所采集的语音信号加入不同的干扰噪声，对加入噪声的信号进行频谱分析，针对受干扰语音信号的特点设计不同的滤波器，对加噪信号进行滤波，恢复原信号。把原始语音信号、加噪语音信号和滤波后的信号进行时域变换和频域变换，画出它们的时域波形和频域波形图，从视觉角度比较分析滤波的效果。也可将这3 类信号进行播放从听觉角度感受滤波的效果。3 2 基本原理2.1 语音信号概述语言是人类创造的，是人类区别于其他地球生命的本质特征之一。人类用语言交流的过程可以看成是一个复杂

4、的通信过程，为了获取便于分析和处理的语音信源，必须将在空气中传播的声波转变为包含语音信息并且记载着声波物理性质的模拟（或数字）电信号，即语音信号，因此语音信号就成为语音的表现形式或载体。语音学和数字信号处理的交叉结合便形成了语音信号处理。语音信号处理是建立在语音学和数字信号处理基础之上的，对语音信号模型进行分析、存储、传输、识别和合成等方面的一门综合性学科。2.2 数字滤波器原理2.2.1 数字滤波器的概念若滤波器的输入、输出都是离散时间信号，那么该滤波器的单位冲激响应h(n)也必然是离散的，这种滤波器称为数字滤波器（DF）。数字滤波器实质上是一种运算过程用来描述离散系统输入与输出关系的差分方

5、程的计算或卷积计算。数字滤波器的实质是用一有限精度算法实现的离散时间线性时不变系统，以完成对信号进行滤波处理的过程。它是数字信号处理的一个重要分支，具有稳定性好、精度高、灵活性强、体积小、质量轻等诸多优点。2.2.2 数字滤波器的分类数字滤波器根据不同的分类标准可以将滤波器分成不同的类别。（1）根据单位冲激响应h(n)的时间特性分类无限冲激响应（IIR）数字滤波器有限冲激响应（FIR）数字滤波器（2）根据实现方法和形式分类递归型数字滤波器非递归型数字滤波器快速卷积型（3）根据频率特性分类低通数字滤波器、高通数字滤波器、带通数字波器、带阻数字滤波器4 3 总体设计思想3.1 语音信号的采集Win

6、dows 自带了一个录音机程序（简称录音机），通过它可以驱动声卡采集、播放和简单处理语音信号。语音信号的采集可以用麦克风直接录制人的语音，也可以通过音频线或者MIDI 线将收音机、电视机或磁带录机中的语音信号采集到计算机中。在录音机中可以进行简单的声音处理，如加大或降低音量，加速或减速，声音的反转或添加回音效果等。加速或减速的改变可以完成变音功能，反转可以达到对声音文件保密的功能。本设计用是收音机接到电脑上用录音机录的，保存起来以备处理。3.2 语音信号处理工具的选择语音信号的进一步处理分析工作选用了matlab 平台。MATLAB 是一种科学计算软件，专门以矩阵的形式处理数据。MATLAB

7、将高性能的数值计算和可视化集成在一起，并提供了大量的内置函数，不断完善MATLAB 产品以提高产品自身的竞争能力Matlab 的数据分析和处理功能十分强大，运用它来进行语音信号的分析、处理和可视化相当便捷。在编程效率、程序可读性、可移植性和可扩充性上matlab远远优于其它的高级编程语言，而且编程易学、直观，代码非常符合人们的思维习惯。另外matlab为用户提供了丰富的 windows 图形界面设计方法，使用户能够在利用其强大的数值计算功能的同时可设计出友好的图形界面，它受到了越来越多的用户的欢迎Matlab 几乎可以在各种机型和操作系统上运行，所以在可移植性和可扩充性上，matlab 远优越

8、于其他的高级编程语言。Matlab 语言具有强大的数值计算能力和视图能力，其偏微分方程工具箱提供了有限元求解的一个强大而灵活的环境，并且有限元网格可做精细划分以满足要求。但是，和其他的高级语言相比，matlab 程序的执行速度较慢。在目前电脑处理速度不断提升的情况下，如果实时性要求不是非常高的情况下，使用matlab 开发就不存在此类问题了。3.3 数字滤波器的设计本设计在matlab 平台上，设计了IIR 和 FIR 滤波器各种类型的滤波器。用以对语音原始信号及加噪的语音信号进行滤波，重点对加噪信号进行了滤波。对于 IIR 滤波器，它的极点可以在单位圆内的任何位置，实现IIR 滤波器的阶次可

9、以较低，所用的存储单元较少，效率高，又由于IIR 数字滤波器能够保留一些模拟滤波器的优良特性，因此应用很广。设计 IIR 数字滤波器的常用方法是从模拟滤波器来设计数字IIR 滤波器。模拟滤波器设计已经有了一套相当成熟的方法，它不但有完整的设计公式，而且还有较为完整的图表供查询，因此充分利用这些已有的资源将会给数字滤波器的设计带来很大的方便。再将模拟滤波器转换为数字滤波器，其转换方法有冲激响应不变法和双线性Z 变换法。一般，当着眼于滤波器的时域瞬态响应时，采用脉冲响应不变法较好，而其它情况下，对于IIR 的设计，大多采用双线性变换法。对于 FIR 滤波器，它有精确、严格的线性相位特性，并且可以做

10、成既是因果的又是稳定的系统。所以FIR 滤波器的应用越来越广泛。IIR 与 FIR 滤波器的性能比较如下图3.1 所示。表 3.1 IIR 滤波器与 FIR 滤波器的性能比较5 IIR DF FIR DF（1）相位一般是非线性的（1）相位可以做到严格线性（2）不一定稳定（2）一定是稳定的（3）不能用FFT作快速卷积（3）信号通过系统可采用快速卷积（4）一定是递归结构（4）主要是非递归结构，也可含递归环节（5）对频率分量的选择性好（零极点可同时起作用）（5）选择性差（6）相同性能下阶次较低（6）相同性能下阶次高（7）有噪声反馈，噪声大（7）噪声小（8）运算误差大，有可能出现极限环振荡（8）运算误

11、差小，不会出现极限环振荡（9）设计有封闭形式的公式，一次完成（9）没有封闭形式的设计公式，须靠经验与反复调试（10）对计算手段的要求较低（10）一般需用计算机计算（11）主要用于设计分段常数的标准低通、高通、带通、带阻和全通滤波器（11）还可设计正交变换器、微分器、线性预测器、回波抵消器、均衡器、线性调频器等各种网络，适用范围广从上面的简单比较我们可以看到IIR 与 FIR 滤波器各有所长，所以在实际应用时应该从多方面考虑来加以选择。例如，从使用要求上来看，在对相位要求不敏感的场合，如语言通讯等，选用IIR 较为合适，这样可以充分发挥其经济高效的特点，而对于图像信号处理，数据传输等以波形携带信

12、息的系统，则对线性相位要求较高，如果有条件，采用FIR 滤波器较好，当然，在实际应用中应考虑经济上的要求以及计算工具的条件等多方面的因素。整体来看，IIR滤波器达到同样效果阶数少，延迟小，但是有稳定性问题，非线性相位；FIR 滤波器没有稳定性问题，线性相位，但阶数多，延迟大。而不论是IIR 滤波器还是FIR 滤波器的设计都包括三个步骤：(1)按照实际任务的要求，确定滤波器的性能指标。(2)用一个因果、稳定的离散线性时不变系统的系统函数逼近性能要求。根据系统的不同的要求可以考虑用 IIR 系统函数或FIR 系统函数去逼近。(3)利用有限精度算法实现系统函数，包括结构选择、字长选择等。6 4 语音

13、信号分析和滤波处理4.1 语音信号的采集把语音信号并保存为.wav 文件，长度小于30 秒，并对语言信号进行采样；录制的软件可以使用Windows自带的录音机，或者也可以使用其它专业的录音软件，录制时需要配备录音硬件（如麦克风），为了方便比较，需要在安静、无噪音、干扰小的环境下录。4.2 语音信号的读入与打开在 MATLAB 中，y,fs,bits=wavread(Blip,N1 N2);用于读取语音，采样值放在向量y 中，fs表示采样频率(Hz)，bits表示采样位数。N1 N2 表示读取的值从N1点到 N2点的值。sound(y);用于对声音的回放。向量y 则就代表了一个信号，也即一个复杂

14、的“函数表达式”，也可以说像处理一个信号的表达式一样处理这个声音信号。下面是语音信号在MATLAB 中的语言程序，它实现了语音的读入与打开，并绘出了语音信号的波形频谱图。x,fs,bits=wavread(wzf.wav);sound(x);X=fft(x,4096);magX=abs(X);angX=angle(X);subplot(221);plot(x);title(原始信号波形);subplot(222);plot(X);title(原始信号频谱);subplot(223);plot(magX);title(原始信号幅值);subplot(224);plot(angX);title(原

15、始信号相位);程序运行可以听到声音，得到的结果图4-1 所示：图 4-1 语音信号的读入与打开4.3 语言信号处理 4.3.1语音信号分析7 用 MATLAB 绘制出语音信号的时域波形图和原始语音信号的频率响应图和原始语音信号的FFT频谱图。程序设计如下：fs=22050;%语音信号采样频率为22050 x,fs,bits=wavread(wzf.wav);sound(x,fs,bits);%播放语音信号y1=fft(x,1024);%对信号做1024 点 FFT变换f=fs*(0:511)/1024;figure(1)plot(x)%做原始语音信号的时域波形图title(原始语音信号时域图)

16、;xlabel(time n);ylabel(fuzhi n);figure(2)freqz(x)%绘制原始语音信号的频率响应图title(频率响应图)figure(3)plot(f,abs(y1(1:512);title(原始语音信号频谱)xlabel(Hz);ylabel(fudu);如下图 4-2 所示：8 图 4-2 语言信号分析4.3.2 含噪语音信号的合成在 MATLAB 软件平台下，给原始的语音信号叠加上噪声，噪声类型分为如下几种：（1）单频噪色（正弦干扰）；（2）高斯随机噪声。绘出加噪声后的语音信号时域和频谱图，在视觉上与原始语音信号图形对比，也可通过Windows 播放软件从

17、听觉上进行对比，分析并体会含噪语音信号频谱和时域波形的改变。程序代码如下：fs=22050;%语音信号采样频率为22050 x,fs,bits=wavread(ding.wav);%读取语音信号的数据，赋给变量x%sound(x)%t=0:1/22050:(size(x)-1)/22050;y1=fft(x,1024);%对信号做 1024 点 FFT变换f=fs*(0:511)/1024;x1=rand(1,length(x);%产生一与x 长度一致的随机信号x2=x1+x;%t=0:(size(x)-1);%加入正弦噪音9%Au=0.3;%d=Au*sin(6*pi*5000*t);%x2

18、=x+d;sound(x2);figure(1)subplot(2,1,1)plot(x)%做原始语音信号的时域图形title(原语音信号时域图)subplot(2,1,2)plot(x2)%做原始语音信号的时域图形title(加高斯噪声后语音信号时域图)xlabel(time n);ylabel(fudu);y2=fft(x2,1024);figure(2)subplot(2,1,1)plot(abs(y1)title(原始语音信号频谱);xlabel(Hz);ylabel(fudu);subplot(2,1,2)plot(abs(y2)title(加噪语音信号频谱);xlabel(Hz);

19、ylabel(fudu);程序运行可以听到声音，得到的结果图4-4 和图 4-5 所示：10 图 4-4 高斯随机噪声图 4-5 正弦 sin 噪声的加噪语音信号11 4.4 滤波器的设计4.4.1 双线性变换法设计了巴特沃斯低通滤波器对加噪语音信号进行滤波对加入高斯随机噪声和正弦噪声的语音信号进行滤波。用双线性变换法设计了巴特沃斯数字低通IIR滤波器对两加噪语音信号进行滤波，并绘制了巴特沃斯低通滤波器的幅度图和两加噪语音信号滤波前后的时域图和频谱图。程序设计如下：x,fs,bits=wavread(wzf.wav);%sound(x)%随机噪声合成x2=rand(1,length(x);%产

20、生一与x 长度一致的随机信号y=x+x2;%加入正弦噪声%t=0:(size(x)-1);%Au=0.3;%d=Au*sin(2*pi*500*t);%y=x+d;wp=0.1*pi;ws=0.4*pi;Rp=1;Rs=15;Fs=22050;Ts=1/Fs;wp1=2/Ts*tan(wp/2);%将模拟指标转换成数字指标ws1=2/Ts*tan(ws/2);N,Wn=buttord(wp1,ws1,Rp,Rs,s);%选择滤波器的最小阶数Z,P,K=buttap(N);%创建 butterworth模拟滤波器Bap,Aap=zp2tf(Z,P,K);b,a=lp2lp(Bap,Aap,Wn)

21、;bz,az=bilinear(b,a,Fs);%用双线性变换法实现模拟滤波器到数字滤波器的转换H,W=freqz(bz,az);%绘制频率响应曲线figure(1)plot(W*Fs/(2*pi),abs(H)grid f1=filter(bz,az,y);figure(2)subplot(2,1,1)plot(t,y)%画出滤波前的时域图title(滤波前的时域波形);subplot(2,1,2)plot(t,f1);%画出滤波后的时域图title(滤波后的时域波形);sound(f1);%播放滤波后的信号F0=fft(f1,1024);f=fs*(0:511)/1024;figure(3

22、)12 y2=fft(y,1024);subplot(2,1,1);plot(f,abs(y2(1:512);%画出滤波前的频谱图title(滤波前的频谱)xlabel(Hz);ylabel(fuzhi);subplot(2,1,2)F1=plot(f,abs(F0(1:512);%画出滤波后的频谱图title(滤波后的频谱)xlabel(Hz);ylabel(fuzhi);程序运行可以播放滤波前面的语音信号，对比滤波前面的语音效果，得到的结果分别图4-6 和图 4-7 所示：13 图 4-6 高斯噪声滤波图 4-7 正弦噪声滤波4.4.2 利用双线性变换实现频率响应S域到 Z域的变换14 利

23、用双线性变换实现频率响应S 域到 Z 域的变换法 设计了巴特沃斯低通数字IIR 滤波器，对加入高斯随机噪声和正弦噪声的语音信号进行滤波，并绘制了两滤波器滤波前后的语音信号时域图和频谱图。程序设计如下：Ft=8000;Fp=1000;Fs=1200;wp=2*pi*Fp/Ft;ws=2*pi*Fs/Ft;fp=2*Ft*tan(wp/2);fs=2*Fs*tan(wp/2);n11,wn11=buttord(wp,ws,1,50,s);%求低通滤波器的阶数和截止频率b11,a11=butter(n11,wn11,s);%求 S 域的频率响应的参数num11,den11=bilinear(b11,

24、a11,0.5);%利用双线性变换实现频率响应S 域到 Z 域的变换x,fs,nbits=wavread(wzf.wav);n=length(x);%求出语音信号的长度t=0:(n-1);x2=rand(1,length(x);%产生一与x 长度一致的随机信号y=x+x2;%加入正弦噪声%t=0:(size(x)-1);%Au=0.03;%d=Au*sin(2*pi*500*t);%y=x+d;figure(1)f2=filter(num11,den11,y)subplot(2,1,1)plot(t,y)title(滤波前的加高斯噪声时域波形);subplot(2,1,2)plot(t,f2)

25、;%画出滤波后的时域图title(滤波后的时域波形);sound(f1);%播放滤波后的信号F0=fft(f1,1024);f=fs*(0:511)/1024;figure(2)y2=fft(y,1024);subplot(2,1,1);plot(f,abs(y2(1:512);%画出滤波前的频谱图title(滤波前加高斯噪声的频谱)xlabel(Hz);ylabel(fuzhi);subplot(2,1,2)F1=plot(f,abs(F0(1:512);%画出滤波后的频谱图title(滤波后的频谱)xlabel(Hz);ylabel(fuzhi);15 程序运行可以播放滤波前面的语音信号，

26、对比滤波前面的语音效果，得到的结果分别图4-8 和图 4-9 所示：图 4-8 高斯噪音滤波16 图 4-9 正弦噪音滤波4.4.3 基于巴特沃斯模拟滤波器设计数字带通滤波器设计了巴特沃斯带通数字IIR滤波器对加了高斯随机噪声和正弦噪声的语音信号进行滤波，并绘制了两滤波器滤波后的语音信号时域图和频谱图。程序设计如下：Wp=0.3*pi,0.7*pi;Ws=0.2*pi,0.8*pi;Ap=1;As=30;N,wn=buttord(Wp/pi,Ws/pi,Ap,As);%计算巴特沃斯滤波器阶次和截止频率b,a=butter(N,wn,bandpass);%频率变换法设计巴特沃斯带通滤波器%db,

27、mag,pha,grd,w=freqz_m(b,a);%数字滤波器响应%Plot(w/pi,mag);%Title(数字滤波器幅频响应|H(ejOmega)|)x,fs,nbits=wavread(wzf.wav);n=length(x);%求出语音信号的长度t=0:(size(x)-1);x2=rand(1,length(x);%产生一与x 长度一致的随机信号y=x+x2;%加入正弦噪声%n=length(x);%求出语音信号的长度%t=0:(n-1);%Au=0.03;%d=Au*sin(2*pi*500*t);%y=x+d;f=filter(bz,az,y);figure(1)freqz

28、(b,1,512)f2=filter(bz,az,y)figure(2)subplot(2,1,1)plot(t,y)title(滤波前的时域波形);17 subplot(2,1,2)plot(t,f2);title(滤波后的时域波形);sound(f2);%播放滤波后的语音信号F0=fft(f1,1024);f=fs*(0:511)/1024;figure(3)y2=fft(y,1024);subplot(2,1,1);plot(f,abs(y2(1:512);%画出滤波前的频谱图title(滤波前的频谱)xlabel(Hz);ylabel(fuzhi);subplot(2,1,2)F1=p

29、lot(f,abs(F0(1:512);%画出滤波后的频谱图title(滤波后的频谱)xlabel(Hz);ylabel(fuzhi);程序运行可以听到声音，得到的结果分别图4-10 和图 4-11 所示：18 图 4-10 高斯噪声滤波19 图 4-17 正弦噪声滤波4.4.4 基于窗函数法的FIR滤器波使用窗函数法,选用海明窗设计了数字FIR 低通滤波器对加了高斯随机噪声和正弦噪声的语音信号进行滤波，并绘制了两滤波器滤波后的语音信号时域图和频谱图。程序设计如下：fs=22050;x,fs,bits=wavread(wzf.wav);%sound(x)t=0:(size(x)-1);x2=r

30、and(1,length(x);%产生一与x 长度一致的随机信号y=x+x2;%加入正弦噪声t=0:(n-1);Au=0.03;d=Au*sin(2*pi*500*t);y=x+d;wp=0.25*pi;ws=0.3*pi;wdelta=ws-wp;N=ceil(6.6*pi/wdelta);%取整wn=(0.2+0.3)*pi/2;b=fir1(N,wn/pi,hamming(N+1);%选择窗函数，并归一化截止频率figure(1)freqz(b,1,512)f2=filter(bz,az,y)figure(2)subplot(2,1,1)plot(t,y)title(滤波前的时域波形);

31、subplot(2,1,2)plot(t,f2);title(滤波后的时域波形);sound(f2);%播放滤波后的语音信号20 F0=fft(f1,1024);f=fs*(0:511)/1024;figure(3)y2=fft(y,1024);subplot(2,1,1);plot(f,abs(y2(1:512);%画出滤波前的频谱图title(滤波前的频谱)xlabel(Hz);ylabel(fuzhi);subplot(2,1,2)F1=plot(f,abs(F0(1:512);%画出滤波后的频谱图title(滤波后的频谱)xlabel(Hz);ylabel(fuzhi);程序运行可以听

32、到声音，得到的结果分别图4-12 和图 4-13 所示：21 图 4-12 高斯噪声滤波图 4-13 正弦噪声滤波22 总结本设计采用了高效快捷的开发工具MATLAB，实现了语音信号的采集，对语音信号加噪声及设计滤波器滤除噪声的一系列工作。从频率响应图中可以看出:巴特沃斯滤波器具有单调下降的幅频特性，通带内是平滑的。海明窗设计的FIR 滤波器的频率特性几乎在任何频带上都比巴特沃斯滤波器的频率特性好，过渡带也比较小，只是海明窗设计的滤波器下降斜度较小对语言的过渡失真进行了补偿。论文初步完成了设计任务，由于本人能力有限，还存在许多不足的地方，比如滤波器的设计种类还比较单一，没有做更多的滤波效果比较

33、等。在以后的工作和学习中会更加努力来完善设计任务。23 致谢经过一个学期的努力和完善，系统终于可以正常的运行。在这里首先感谢我的毕业设计指导教师黄小芬老师，在这段时间一直给我的支持与鼓励。认真负责的监督我们毕业设计的进度，耐心的指导我们使我们能够按时的完成任务。还有特别感谢一下陈志标老师对我硬件部分提出的宝贵的建议，使得我硬件部分能够比较顺利的完成。还有感谢学校为我们提供的良好实验环境以及充足的实验设备，为我们的设计和调试提供了很大的方便。在这段时间学到了很多，虽然由于自身的不足没有能够为系统提出更好的解决方案。但这对我来说绝对是一个非常宝贵的历练。从中我切身体会到了理论和现实的差距，只有真正

34、动手去做才能发现问题。最后，感谢所有对我提供过帮助的人。谢谢。24 参考文献1(美)维纳?K?恩格尔,(美)约翰?G?普罗克斯.数字信号处理M.西安：西安交通大学出版社，2002.2 黄文梅，熊桂林，杨勇.信号分析与处理M.长沙：国防科技大学出版社，2000.6 陈怀琛.数字信号处理教程M.北京：电子工业出版社，2004.7 程佩青.数字信号处理教程（第二版)M.北京：清华大学出版社，2001.8 韩纪庆，张磊，郑铁然.语音信号处理M.北京：清华大学出版社，2004.9 易克初，田斌，付强.语音信号处理M.北京：国防工业出版社，2000.10 周辉，董正宏.数字信号处理基础及MATLAB 实现

35、 M.北京：北京希望电子出版社，2006.11 邹理和.语音信号处理 M.北京：国防工业出版社，1985.25 Matlab-based voice signal analysis and processingConcord University College Fujian Normal University Department of Information Technology Electronics and Information Science and Technology Professional 124122006028 WangZhenfei xiao feng Huang m

36、aster science【Abstract】This design using Microsoft recorder collected a speech,analyses the time-domain analysis,the spectrum analysis,the analysis of the characteristics of speech signal.Matlab and application of speech signal joined with different noise,further bilinear transformation method for t

37、he design of a filter based on butterworth simulation design digital filter bandpass filter and the design method of window functions with FIR low-pass filter,and then to the speech signal with noise filtered.Compares the speech signal after filtering time domain and frequency domain properties,play

38、back speech signal and noise and speech signal de-noising compared and studied butterworth beth IIR filter and FIR filters in signal processing,especially the signal filtering de-noising characteristics and applications.From the theory and practice of the filter is different digital filter.【Keywords】speech signal,Frequency characteristics,Temporal characteristics,Filter