1.基于倒谱及贝叶斯决策的周期信号基频提取与匹配方法,其特征在于包括如下步骤:步骤1、根据T种不同周期的周期信号的基频,得到其组成的基频集合S={μ1,μ2,...,μT};
步骤2、获得未知基频的周期信号数据s(m),m=1,2,...,Ns,周期信号数据s(m)的基频μx属于基频集合S,Ns为数据采样点数;
步骤3、按照帧长Lf、帧移Lo对周期信号数据s(m)进行分帧,得到包含Nf个数据帧的数据帧集合:步骤4、对每一帧数据xi(n),i=1,2,...,Nf,n=1,2,...,Lf,使用倒谱法提取出基频估计值ffi:单帧信号xi(n)帧长Lf,采样频率Fs,针对其倒谱法基频提取过程如下:
4.1、使用快速傅里叶变换(FFT)得到xi(n)的频谱X(ω):X(ω)=FFT[x(n)],ω=0,Fs/Lf,...,(Lf-1)*Fs/Lf;
4.2、依据频谱X(ω)计算得到幅值谱Xa(ω):Xa(ω)=|X(ω)|,ω=0,Fs/Lf,...,(Lf-1)*Fs/Lf;
4.3、依据幅值谱Xa(ω)计算得到幅值对数谱Xl(ω):Xl(ω)=log[Xa(ω)],ω=0,Fs/Lf,...,(Lf-1)*Fs/Lf;
4.4、将幅值对数谱Xl(ω)进行快速傅里叶逆变换(IFFT)得到倒谱
4.5、根据基频检索频率的范围[fmin,fmax],计算得到倒谱索引范围[nl,nr]:nl=Fs/fmax,nr=Fs/fmin;
4.6、在倒谱序列 的索引范围[nl,nr]中,找到最大值 及对应的索引:
4.7、根据ntop计算得到第i帧数据的基频的估计值ffi:ffi=Fs/ntop
最终将所有帧的基频估计值组成基频估计值序列{ffi},i=1,2,...,Nf;
步骤5、求基频估计值序列{ffi}的均值:步骤6、依据贝叶斯决策得到准确的基频μx,即: