学士学位论文—数字滤波在语音信号的应用张净丹.docVIP

数字滤波器在语音 信号的应用 姓 名： 学 号： 000000 专 业： 电气工程及自动化 授课教师： 00000 2016年11月4日—程序—附件—娱乐—录音机，文件—属性—立即转换—8000KHz，8位，单声道）或调用系统声音。然后在MATLAB软件平台下，利用函数audioread对语音信号进行采样，记住采样频率和采样点数。通过audioread函数的使用进一步理解采样频率、采样位数等概念。 2.1.2方法 本报告介绍了基于Matlab的对语音信号采集、处理及滤波器的设计，并使之实现的过程。理解与掌握课程中的基本概念、基本原理、基本分析方法，用Matlab进行数字语音信号处理，并阐述了课程设计的具体方法、步骤和内容。综合运用本课程的理论知识进行频谱分析以及滤波器设计，通过理论推导得出相应结论，并利用MATLAB作为工具进行实现，从而复习巩固课堂所学的理论知识，提高对所学知识的综合应用能力，并从实践上初步实现对数字信号的处理。 Matlab 语言是一种数据分析和处理功能十分强大的计算机应用软件 ,它可以将声音文件变换为离散的数据文件 , 然后利用其强大的矩阵运算能力处理数据 ,如数字滤波 、傅里叶变换 、时域和频域分析 、声音回放以及各种图的呈现等 , 信号处理是 Matlab 重要应用的领域之一 。本文是用 Matlab对含噪的的语音信号同时在时域和频域进行滤波处理和分析。整个设计思路的设计原理示意图如图所示： 2.1.3语音信号的读取 在matlab软件平台下，利用audioread函数对语音信号进行采样，记住采样频率和采样点数。通过使用audioread函数，理解采样频率、采样位数等概念。 audioread函数调用格式 y=audioread(file)，读取file所规定的音乐文件，返回采样值放在向量y中。 [y,fs,nbits]=audioread(file)，采样值放在向量y中，fs表示采样频率（hz），nbits表示采样位数。 y=audioread（file，N），读取前N点的采样值放在向量y中。 y=audioread（file，[N1,N2]），读取从N1到N2点的采样值放在向量y中。 对语音信号“China.wav”进行采样其程序如下： [y,fs,nbits]=audioread(China.wav) 可知该语音信号的采样频率fs = 8000Hz ,数据位nbits = 16Bit。 2.1.4语音信号的频谱分析 首先画出语音信号的时域波形，然后对语音信号进行频谱分析。在Matlab 中可以利用函数fft 对信号行快速傅里叶变换,得到信号的频谱特性，从而加深对频谱特性的理解。 其程序如下： [y,fs,nbits]=audioread (China.wav ); %把语音信号进行加载入Matlab 仿真软件平台中 sound(y,fs,nbits); %回放语音信号 p = audioplayer(x1,44100);play(p)%常用语音播放 n = length (y) ; %求出语音信号的长度 Y=fft(y,n); %快速傅里叶变换 subplot(2,1,1);plot(y);title(原始信号波形); subplot(2,1,2);plot(abs(Y));title(原始信号频谱) 当然我们也可以利用系统自带的声音进行模拟实验，调用系统声音程序如下： MATLAB也自带了很多声音： %鸟声 %火车 load train load chirp sound(y,Fs) sound(y,Fs) 2.1.5程序结构图 2.2实验结果 说明：本实验声音取自于系统，程序设计如下： %鸟声 load chirp sound(x1,Fs) %加入噪声后声音播放 load chirp t=0:1/13129:(size(y)-1)/13129; Au=0.03; d=[Au*cos(2*pi*10000*t)]; x2=y+d; p = audioplayer(x2,13129); play(p) %模拟滤波器 wp=0.25*pi; ws=0.3*pi; Rp=1; Rs=15; Fs=13129; Ts=1/Fs; wp1=2/Ts*tan(wp/2);