1.基于EMD重构的相关时延估计方法,其特征是,包括如下步骤:
1)采集两路目标信号,获得两路目标信号的采样序列y1(n)、y2(n),并提取两路目标信号信道的前导背景噪声序列b1(n)、b2(n),因为两路目标信号到达两接收传感器时存在时延差σ,因此采样序列y1(n)与y2(n)之间也存在时延σ;
2)对采样序列y1(n)进行离散傅里叶变换得到Y1(ω),将Y1(ω)的幅值取对数得到再对 进行傅里叶逆变换得到采样序列y1(n)的倒谱序列 通过同样的方式得到目标信号信道前导背景噪声序列b1(n)的倒谱序列
3)将 进行滤波,分别得到 的声道响应 与 的声道响应
对 分别进行傅里叶变换得到 的频谱包络线 的频谱
包络线 对 做差消除对等噪声,并除以瞬时背景噪声能量 获
得谱差曲线D1(ω),再对谱差曲线D1(ω)进行归一化得到归一化谱差曲线
4)用门限法获取 幅值高于门限T1的频段区域,假设满足条件的频段区域有k个,对应目标信号中所含语音信号成分主导频段记为:[ωp11,ωq11]、......、[ωp1k,ωq1k];
5)对采样序列y1(n)进行EMD分解,依次得到m个基本模式分量,记为:h11(n),h12(n),……,h1m(n),和一个分解后的余量r1(n),这样y1(n)就可以表示为基本模式分量和余项的和,即
6)通过傅里叶变换计算并获取每组分量h11(n)......h1m(n)的功率谱分布曲线H11(ω)......H1m(ω);计算每组分量的功率谱分布曲线在频段[ωp11,ωq11]、......、[ωp1k,ωq1k]内的幅值累加与全局幅值累加的比值η1,如果对应分量计算得到的η1大于预设门限值TH,则该分量作为一路目标信号组成重构信号C1(n)的分量之一,否则将该分量丢弃;
7)将另一路采样序列y2(n)按照步骤2)-步骤6)进行同样的操作,获取另一路目标信号重构信号C2(n),两路目标信号重构后的信号分别表示为C1(n)和C2(n),将C1(n)和C2(n)进行二次相关得到C1(n)、C2(n)的二次相关序列 对二次相关序列 的峰值进行检测可得到时延估计值
2.根据权利要求1所述的基于EMD重构的相关时延估计方法,其特征是,所述采样序列y1(n)、y2(n)是两个传感器同步采集的时间采样序列,采样频率为Fs=8000Hz。
3.根据权利要求1所述的基于EMD重构的相关时延估计方法,其特征是,所述的门限T1为0.7-0.9。
4.根据权利要求1所述的基于EMD重构的相关时延估计方法,其特征是,所述的门限T1为0.8。
5.根据权利要求1所述的基于EMD重构的相关时延估计方法,其特征是:所述的预设门限值TH为0.5。
6.根据权利要求1所述的基于EMD重构的相关时延估计方法,其特征是:所述的谱差曲线D1(ω)为: