有噪声语音信号处理的理论概述

目录

1.语音信号的采样

2.IIR数字滤波器概述

3.FIR数字滤波器概述

3.数字滤波器概述

---

       数字滤波器, 是数字信号处理中及其重要的一部分。随着信息时代和数字技术的发展，受到人们越来越多的重视。数字滤波器可以通过数值运算实现滤波，所以数字滤波器处理精度高、稳定、体积小、重量轻、灵活不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊功能。数字滤波器种类很多，根据其实现的网络结构或者其冲激响应函数的时域特性，可分为两种，即有限冲激响应( FIR，Finite Impulse Response)滤波器和无限冲激响应( IIR，Infinite Impulse Response)滤波器。

1.语音信号的采样

       在进行模拟/数字信号的转换过程中，当采样频率fs.max大于信号中，最高频率fmax的2倍时，即：fs.max>=2fmax,则采样之后的数字信号完整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的5～10倍；采样定理又称奈奎斯特定理。 理想低通信道的最高大码元传输速率=2W*log2 N (其中W是理想低通信道的带宽,N是电平强度)

      在MATLAB软件平台下，利用wavread函数对语音信号进行采样，记住采样频率和采样点数，Wavread函数调用格式：

y=wavread（file）%读取file所规定的wav文件，返回采样值放在向量y中。

[y,fs,bits]=wavread(file) %采样值放在向量y中，fs表示采样频率（hz），bits表示采样位数。

y=wavread（file，N）%读取前N点的采样值放在向量y中。

y=wavread（file，[N1,N2]）%读取从N1到N2点的采样值放在向量y中。

对语音信号Orisound.wav进行采样其程序如下：

[y,fs,nbits]=wavread (‘OriSound’); %把语音信号加载入MATLAB仿真软件平台中

       画出语音信号的时域波形，再对语音信号进行频谱分析。MATLAB提供了快速傅里叶变换算法FFT计算DFT的函数fft,其调用格式如下：

Xk=fft(xn,N)        

       参数xn为被变换的时域序列向量，N是DFT变换区间长度，当N大于xn的长度时，fft函数自动在xn后面补零。，当N小于xn的长度时，fft函数计算xn的前N个元素，忽略其后面的元素。在本次设计中，我们利用fft对语音信号进行快速傅里叶变换，就可以得到信号的频谱特性。程序如下：

[y,fs,nbits]=wavread(‘OriSound’);       %语音信号的采集

sound(y,fs,nbits);                     %语音信号的播放

n=length(y) ;

Y=fft(y,n);                          %快速傅里叶变换

figure;

subplot(2,1,1);

plot(y);

title(‘原始信号波形’,’fontweight’,’bold’);

axis([ 78000 80000 -1 1]);

grid;

subplot(2,1,2);

plot(abs(Y));

title(‘原始信号频谱’,’fontweight’,’bold’);

axis([ 0 150000 0 4000]);

grid;

2.IIR数字滤波器概述

       IIR(Infinite Impulse Response)数字滤波器，又名“无限脉冲响应数字滤波器”，或“递归滤波器”。递归滤波器，也就是IIR数字滤波器，顾名思义，具有反馈，一般认为具有无限的脉冲响应。IIR滤波器有以下几个特点：

（1）封闭函数：IIR数字滤波器的系统函数可以写成封闭函数的形式。

（2）IIR数字滤波器采用递归型结构：IIR数字滤波器采用递归型结构，即结构上带有反馈环路。IIR滤波器运算结构通常由延时、乘以系数和相加等基本运算组成，可以组合成直接型、正准型、级联型、并联型四种结构形式，都具有反馈回路。由于运算中的舍入处理，使误差不断累积，有时会产生微弱的寄生振荡。

（3）借助成熟的模拟滤波器的成果：IIR数字滤波器在设计上可以借助成熟的模拟滤波器的成果，如巴特沃斯、契比雪夫和椭圆滤波器等，有现成的设计数据或图表可查，其设计工作量比较小，对计算工具的要求不高。在设计一个IIR数字滤波器时，我们根据指标先写出模拟滤波器的公式，再通过一定的变换，将模拟滤波器的公式转换成数字滤波器的公式。

（4）需加相位校准网络：IIR数字滤波器的相位特性不好控制，对相位要求较高时，需加相位校准网络。

3.FIR数字滤波器概述

        FIR数字滤波器（finite impulse response filter）又名“有限脉冲响应数字滤波器”，这类滤波器对于脉冲输入信号的响应最终趋向于0，因此而得名。有限脉冲响应滤波器（FIR filter）的优点：

（1） 脉冲响应（impulse response）为有限长：造成当输入数位讯号为有限长的时候，输出数位讯号也为有限长。

（2） 比无限脉冲响应滤波器（IIR filter）较容易最佳化（optimize）。

（3） 线性相位（linear phase）：造成h(n),是偶对称（even）或奇对称（odd）且有限长。

（4）一定是稳定的（stable）：因为Z转换（Z transform）后所有的极点（pole）都在单位圆内

        低通滤波器：对于不同滤波器而言，每个频率的信号的减弱程度不同。当使用在音频应用时，它有时被称为高频剪切滤波器, 或高音消除滤波器。低通滤波器概念有许多不同的形式，其中包括电子线路（如音频设备中使用的hiss 滤波器、平滑数据的数字算法、音障（acoustic barriers）、图像模糊处理等等，这两个工具都通过剔除短期波动、保留长期发展趋势提供了信号的平滑形式。低通滤波器在信号处理中的作用等同于其它领域如金融领域中移动平均数（moving average)所起的作用；低通滤波器有很多种，其中，最通用的就是巴特沃斯滤波器和切比雪夫滤波器。

       高通滤波器：去掉信号中不必要的低频成分，去掉低频干扰的滤波器。在电力系统中，谐波补偿时用高通滤波器滤除某次及其以上的各次谐波。高通滤波器是指车载功放中能够让中、高频信号通过而不让低频信号通过的电路，其作用是滤去音频信号中的低音成分，增强中音和高音成分以驱动扬声器的中音和高音单元。此外高通滤波器常常和低通滤波器成对出现，不论哪一种，都是为了把一定的声音频率送到应该去的单元。

       带通滤波器：能通过某一频率范围内的频率分量、但将其他范围的频率分量衰减到极低水平的滤波器，与带阻滤波器的概念相对。一个模拟带通滤波器的例子是电阻-电感-电容电路(RLC circuit)。这些滤波器也可以用低通滤波器同高通滤波器组合来产生。

3.数字滤波器概述

       用自己设计的各滤波器分别对加噪的语音信号进行滤波，在Matlab中，FIR滤波器利用函数fftfilt对信号进行滤波，IIR滤波器利用函数filter对信号进行滤波。

       函数fftfilt用的是重叠相加法实现线性卷积的计算。调用：y=fftfilter(h,x,M)。其中，h是系统单位冲击响应向量；x是输入序列向量；y是系统的输出序列向量；M是有用户选择的输入序列的分段长度，缺省时，默认的输入向量的重长度M=512。

       函数filter的调用格式：yn=filter(B，A.xn)，它是按照直线型结构实现对xn的滤波。其中xn是输入信号向量，yn输出信号向量。由于所设计的滤波器有六个，故我们需要进行六次滤波，然后比较得出哪种滤波器的效果最好。