语音信号的数字滤波处理——巴特沃思、Blackman窗滤波器设计.docxVIP

第1页共55页

1绪论

1.1设计目的

通过本次数字信号处理课程设计，综合运用数字信号处理技术课程和其他有关课程的理论和生产实际知识去分析和解决具体问题，并使所学知识得到进一步巩固、深化和发展。学习设计滤波器的一般方法，通过课程设计树立正确的设计思想，提高分析问题、解决问题的能力。进行设计基本技能的训练，如查阅设

计资料和手册、程序的设计、调试等。

1.2设计要求

(1)掌握数字信号处理的基本概念，基本理论和基本方法。(2)熟悉离散信号和系统的时域特性。

(3)掌握序列快速傅里叶变换方法。

(4)学会MATLAB的使用，掌握MATLAB的程序设计方法。

(5)掌握利用MATLAB对语音信号进行频谱分析。

(6)掌握滤波器的网络结构。

1.3设计内容

1.3.1预习题部分

(1)设计卷积运算的演示程序：

可输入任意两个序列x1(n)、x2(n),指定x1(n)为自己的学号，例如

x1(n)={2,0,0,7,5,7,1,7,0,2,2,7}。x2(n)的内容和长度自选。例如x2(n)={1,2.43,6.17,12.93,22.17,32.25,40.88,45.87,45.87,40.88,32.25,22.17,

12.93,6.17,2.43,1.0000}。

(2)编写程序演示采样定理(时域采样、频谱周期延拓),同时演示采样频率小于

2fc时，产生的混叠效应：

①对下面连续信号进行采样：

x。()=Aesin(Q?I)u(1),A=n,a=Q?=n√2π,n=学号，A为幅度因子，a为衰

减因子，Q,为模拟角频率，其中n为学号(例如，王墨同学n=23)

第2页共55页

②要求输入采样频率fs(根据程序处理需要指定范围)后，在时域演示信号波

形、采样脉冲及采样后信号；在频域演示不同采样频率下对应信号的频谱

1.3.2设计题部分

(1)利用Windows下的录音机或其他软件，选择Windows系统的“叮·●

(ding.wav),并对该信号进行采样；

(2)语音信号的频谱分析，画出采样后语音信号的时域波形和频谱图；

(3)产生噪声信号加到语音信号中，得到被污染的语音信号，并回放语音信号；

(4)污染信号的频谱分析，画出被污染的语音信号时域波形和频谱；

(5)根据有关的频谱特性，采用间接法设计IIR数字滤波器，并画出相应滤波器

的幅频图(设计3个IIR滤波器)

a.模拟滤波器类型：巴特沃思滤波器(低通、带通、高通)

b.总体要求：MATLAB原程序+仿真波形+技术指标

(6)根据有关的频谱特性，采用直接法设计FIR数字滤波器，并画出相应滤波器

的幅频图(设计3个FIR滤波器)

a.滤波器类型：Blackman窗(低通、带通、高通)

b.总体要求：MATLAB原程序+仿真波形+技术指标+窗函数

(7)用自己设计的这些滤波器分别对被不同噪声污染的信号进行滤波；

(8)分析得到信号的频谱，画出滤波后信号的时域波形和频谱，并对滤波前后的

信号进行对比，分析信号的变化；

(9)回放语音信号。

2课程设计方案及相应原理

2.1课程设计方案

选择windows自带的“ding.wav”声音作为语音信号，产生噪声信号加到语音信号中，得到被污染的语音信号，经过频谱分析，分别设计巴特沃思滤波器以及Blackman窗滤波器(低通、带通、高通)对被污染的语音信号滤波，滤掉

噪音信号，得到符合要求的语音信号，方案设计流程图如图2.1、图2.2所示。

第3页共55页

仿真开始

仿真开始

噪声信号

语音信号(dingwav)

语音信号+噪声信号

巴特沃思滤波器滤波

(低通、带通、高通)

滤波后的信号

仿真结束

时域波形

幅频图

图2.1总体方案设计流程图(巴特沃思)

仿

仿真开始

噪声信号

语音信号(ding.wav)

语音信号+噪声信号

Blackman窗滤波器滤波

(低通、带通、高通)

滤波后的信号

仿真结束

时域波形

幅频图

图2.2总体方案设计流程图(Blackman窗)

2.2课程设计原理

第4页共55页

设计数字滤波器的任务就是寻求一个因果稳定的线性时不变系统，并使系统函数H(z)具有指定的频率特性。数字滤波器从实现的网络结构或者从单位冲激响应分类，可以分成无限长单位冲激响应(IIR)数字滤波器和有限长单位冲激

响应(FIR)数字滤波器。

本次数字信号处理课程设计我们小组选用巴特沃思数字滤波器(IIR)以及

Bl