基于matlab的语音信号处理(4)

grid;

subplot(2,2,3); plot(y);

title('滤波前信号的波形','fontweight','bold');

grid;

subplot(2,2,4);

plot(z11);

title('滤波后的信号波形','fontweight','bold');

grid;

由程序运行结果，如图4-3可知，所涉及的滤波器符合要求。

图4-3验证所涉及的滤波器

4.5滤波

用自己设计的各滤波器分别对加噪的语音信号进行滤波，在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输出信号向量。

由于所设计的滤波器有六个，故我们需要进行六次滤波，然后比较得出哪种滤波器的效果最好。

对原始语音信号具体程序及运行结果如附录II比较滤波前后语音信号的波形及频谱