语音实验二信号的基音周期提取.docVIP

实验二语音信号的基音周期提取? 一、实验目的? 1、熟练运用MATLAB软件的运用，学习通过MATLAB软件编程来进行语音信号的基因周期提取。? 2、掌握语音信号的基音周期提取的方法，实现其中一种基频提取方法。? 3、学会用自相关法进行语音信号的基因检测。? 二、实验仪器设备及软件? HP?D538、MATLAB? 三、实验原理? 浊音信号的自相关函数在基因周期的整数倍位置上出现峰值，而清音的自相关函数没有明显的峰值出现。因此检测自相关函数是否有峰值就可以判断是清音还是浊音，而峰-峰值之间对应的就是基音周期。? 影响从自相关函数中正确提取基音周期的最主要原因是声道响应。当基音的周期性和共振峰的周期性混在一起时，被检测出来的峰值可能会偏离原来峰值的真实位置。另外，在某些浊音中，第一共振频率可能会等于或低于基音频率。此时，如果其幅度很高，它就可能在自相关函数中产生一个峰值，而该峰值又可以同基音频率的峰值相比拟。? 1、自相关函数?? 对于离散的语音信号x(n)，它的自相关函数定义为：? R(k)=Σx(n)x(n-k)，? 如果信号x(n))具有周期性，那么它的自相关函数也具有周期性，而且周期与信号x(n)?的周期性相同。自相关函数提供了一种获取周期信号周期的方法。在周期信号周期的整数倍上，它的自相关函数可以达到最大值，因此可以不考虑起始时间，而从自相关函数的第一个最大值的位置估计出信号的基音周期，这使自相关函数成为信号基音周期估计的一种工具。? ?2、短时自相关函数? 语音信号是非平稳的信号，所以对信号的处理都使用短时自相关函数。短时自相关函数是在信号的第N个样本点附近用短时窗截取一段信号，做自相关计算所得的结果? Rm(k)=Σx(n)x(n-k)?? 式中，n表示窗函数是从第n点开始加入。?? 3、算法? 通过对自相关基音检测原理的分析,考虑到检测准确度和检测速率2方面的因素,提出了算法实现方案,并对算法进行了Matlab编程实现。算法包含6个功能模块:带通滤波、取样、分帧、短时能量分析、相关运算、基音检测。? 语音信号???? 图2.1?基音检测框图? 四、实验步骤及程序? （1）实验步骤? 1、取一段录音作为音频样本。 ?2、对样本音频进行采样。 ?3、对采样后样本音频进行分帧。 ?4、对每一帧求短时自相关函数 ?5、算出对应周期。? 6、进行带通滤波去除野点。? （2）语音信号基音周期提取的程序流程图：? ??????? ? ???? 图2.2?语音信号基音周期提取的程序流程图 （3）实验源程序 [x,fs]=wavread(1.wav); [f0] = pitch_sift(x,0.38,fs); plot(f0); pitch_sift定义 % [f0] = pitch_sift(x,THRESHOLD,SAMP_FREQ,FRAME_SIZE,STEP_SIZE) % THRESHOLD is an optional argument between 0 and 1 for voiced-unvoiced % decision,defult is 0.38 % FRAME_SIZE is the frame size(non-overlapping) for the short time analysis.defult is 200 % STEP_SIZE is the step size between successive frames.defult is 100. % Pitch estimation by Autocorrelation method % % Levent Arslan August 4, 1996 function [f0,p_voice] = pitch_sift(x,THRESHOLD,SAMP_FREQ,FRAME_SIZE,STEP_SIZE); if isstr(x), x=readsp(x); end; if nargin 4, FRAME_SIZE = SAMP_FREQ*25/1000; end; %25ms if nargin 5, STEP_SIZE = 0.4*FRAME_SIZE; end; %10ms if nargin 3, SAMP_FREQ=16000; end% old is 8000 if nargin 2, THRESHOLD = 0.28; end%old is 0.38 norm = FRAME_SIZE/256; P=ceil(SAMP_FREQ/2000); M=P; b =[0.0003 0.0020 0.0051 0.0068 0.0051 0.0020 0