设计IIR,FIR滤波器对语音信号进行滤波

题目设计IIR，FIR滤波器对语音信号进行滤波

专业

班级

学号

姓名

1、主要内容

1)掌握IIR,FIR滤波器设计方法

2)掌握使用MATLAB设计滤波器的方法，并能够编程对采样的语音信号进

行滤波。

3)掌握利用MATLAB对信号时频域进行分析的方法。

2、基本要求

用MATLAB编程实现下面的要求：

1)得到语音信号

2)分析语音信号的频谱

3)根据语音信号的频谱设计滤波器的技术指标

4)根据滤波器的技术指标设计FIR，或者IIR滤波器。

5)对语音信号进行滤波，分析滤波结果

3、主要参考文献

1)高西全，丁玉美。数字信号处理（第三版）。西安电子科技大学出版社

出版社。2010

2)张德丰编著。MATLAB数字信号处理与应用。清华大学出版社。2010

3)（美)英格尔(Ingle,V.K.) 等著。数字信号处理(MATLAB版)。西安交

通大学出版社。2008

4)丛玉良, 王志宏编著。数字信号处理原理及其MATLAB实现。电子工业

出版社。2009

完成期限：

指导教师签名：

课程负责人签名：

程序代码：

[z1,fs,bite]=wavread('D:\S.H.E.wav',[400000,800000]);%对语音信号进行采样 sound(z1,fs,bite); %播放原始语音信号

z=z1(:,1); %由双通道选单通道

subplot(2,2,1);plot(z);title('时域波形'); %信号的时域波形

y1=fft(z,1024); %信号的傅里叶变换

subplot(2,2,2);plot(abs(y1));title('幅频特性'); %信号的频谱特性性

wp=0.018*pi;ws=0.131*pi; %通带便捷频率，阻带截止频率 Rp=1;As=30; %通带衰减，阻带衰减

wc=(ws+wp)/2; %通带截止频率

tr_width=ws-wp; %过渡带宽度

M=ceil(6.2*pi/tr_width)+1;N=0:1:M-1; %h 【n 】的长度

alpha=(M-1)/2;

n=0:1:(M-1);

m=n-alpha+eps;

hd=sin(wc*m)./(pi*m); %频率响应hd 【n 】

w_hanning=(hanning(M))'; %窗函数

h=hd.*w_hanning; %h[n]

x=filter(h,[1],z); %滤波

subplot(2,2,3);plot(x);title('滤波后时域波形');

x1=fft(x,1024); %求滤波后信号的傅里叶变换 subplot(2,2,4);plot(abs(x1));title('滤波后幅频特性');%滤波后信号幅频特性

wavwrite(x,fs,bite,'d:\滤波.wav')

sound(x,fs,bite); %滤波后的语音信号

024

6x 10

5-1-0.5

0.5

1

时域波形

05001000150005

10幅频特性024

6x 105-1-0.5

0.5

1

滤波后时域波形

050010001500

05

10滤波后幅频特性