数字信号处理课程设计报告语音信号的数字滤波处理大学论文

信号处理综合设计

设计说明书

设计项目：语音信号的数字滤波处理

项目完成人：姓名：学号：

专业班级：电子信息工程二班

指导教师：

提交日期： 2017.1.7

机电工程学院电子信息工程

1.设计基本原理

1.1课题研究的背景

1)熟悉离散信号和系统的时域特性。

2)掌握数字信号处理的基本概念，基本理论和基本方法。

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

4)练习对MATLAB的使用，掌握MATLAB的程序设计方法。

5)学习利用MATLAB对语音信号进行频谱分析。

6)学习设计巴特沃斯滤波器以及采用hamming窗设计FIR数字低通，高通和带通滤波器的

方法。

1.2课题研究意义

通过本次课程设计，综合运用数字信号处理技术课程和其他有关先修课程的理论和生产实际知识去分析和解决具体问题，并使所学知识得到进一步巩固、深化和发展。初步培养学生对工程设计的独立工作能力，学习设计的一般方法。通过课程设计树立正确的设计思想，提高学生分析问题、解决问题的能力。进行设计基本技能的训练，如查阅设计资料和手册、程序的设计、调试，了解语音信号的产生、采集、能绘制语音信号的频率响应曲线及频谱图，学会用MATLAB对语音信号进行分析和处理，学习设计滤波器的一般方法，通过课程设计树立正确的设计思想，提高分析问题、解决问题的能力，掌握用滤波器去除语音信号噪声的方法，观察去噪前后的语音信号，进行设计基本技能的训练，如查阅设计资料和手册、程序的设计、调试等。

设计任务

2.1课题设计的任务

(1)

(2) 熟悉离散信号和离散系统的时域特性。 (3) 掌握序列快速傅里叶变换（FFT ）方法。 (4) 掌握滤波器的几种网络结构。

(5) 学会使用软件MATLAB ，掌握其程序设计方法。 (6) 掌握利用MATLAB 对语音信号进行频谱分析的方法。

2.2课题研究的内容

预习题部分

1、设计卷积运算的演示程序:

① 可输入任意两个序列x1(n)、x2(n)，指定x1(n)为自己的学号，例如x1(n)={2, 0, 1, 0, 5, 7, 0, 5, 0, 2, 0, 1}。x2(n)的内容和长度自选。例如x2(n)={0.5, 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.007}。

② 分别动态演示两个序列进行线性卷积x1(n)*x2(n)和圆周卷积x1(n)⊙x2(n)的过程；要求分别动态演示翻转、移位、乘积、求和的过程；

③ 圆周卷积默认使用两个序列中的最大长度，但卷积前可以指定卷积长度N 用以进行混叠分析；

④ 改变圆周卷积长度N ，根据实验结果分析两类卷积的关系；

⑤ 在计算机操作系统中选一段声音文件(XP 系统在“C:\WINDOWS\Media ”)，读取文件取10ms 的声音数据产生时域序列x1(n)，序列内容自定义。利用x2(n)={0.5, 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.007}。利用FFT 实现快速卷积，验证时域卷积定理，并与直接卷积进行效率对比（验证时采用Matlab 子函数）。

2、编写程序演示采样定理（时域采样、频谱周期延拓），同时演示采样频率小于2f c 时，产生的混叠效应：

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

00()sin()(),at a x t Ae t u t A n a n -=Ω==Ω==，，学号，A 为幅度因子，a

为衰减因子，0Ω为模拟角频率，其中n 为学号(例如，贺娜同学n=201)；

②要求输入采样频率f s（根据程序处理需要指定范围）后，在时域演示信号波形、

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

设计题部分

数字滤波器设计及其应用

①利用Windows下的录音机或其他软件，进行语音信号的采集（*.wav）；

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

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

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

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

幅频、相频图（设计3个IIR滤波器）

a．模拟滤波器类型：巴特沃斯滤波器（低通、带通、高通）

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

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

幅频、相频图（设计3个FIR滤波器）

a．滤波器类型：汉宁(Hanning)窗（低通、带通、高通）

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

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

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

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

回放语音信号。

设计提示

1．语音信号的采集

要求利用windows下的录音机（开始—程序—附件—娱乐—录音机，文件—属性—立即转换—8000Hz，8位，单声道）录制一段自己的话音，或者采用Windows自带的声音文件（默认为22050Hz），时间控制在几秒左右。然后在MATLAB软件平台下，利用函数wavread 对语音信号进行采样，记住采样频率和采样点数。通过wavread函数的使用，要求理解采样频率、采样位数等概念。

wavread函数调用格式：

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

[y,fs,nbits]=wavread(file)，采样值放在向量y中，fs表示采样频率（Hz），nbits