[发明专利]基于鲁棒自适应盲系统辨识的多通道时延估计方法

摘要

本发明公开了一种基于鲁棒自适应盲系统辨识的多通道时延估计方法，初始化各传声器对应通道的自适应滤波器系数，然后各传声器分别对声信号进行采集并分帧保存，利用各传声器在同一时间段采集的信号帧进行多通道时延估计：首先由各传声器采集的时域信号与其前一帧时域信号、前一帧时域信号对应的自适应滤波器系数以及Cauchy估计器的一阶和二阶导数计算自适应滤波器系数的更新量，利用该更新量对前一次迭代获得的自适应滤波器系数进行更新，得到当前帧对应的自适应滤波器系数，然后根据计算得到的自适应滤波器系数进行时延估计。本发明使用Cauchy估计器定义新的频域代价函数，提升时延估计的性能。

主权项

1.一种基于鲁棒自适应盲系统辨识的多通道时延估计方法，其特征在于，包括以下步骤：S1：初始化自适应滤波器系数其中为长度为L的系数列向量，k＝1,2,…,M，设置自适应步长μf、遗忘因子λ、Cauchy估计器的参数c；S2：M个传声器分别对声源信号进行持续采集，每个传声器将采集的时域信号以每L个样本为一帧进行保存，第i个传声器的第m帧时域信号记为xi,L(m)＝[xi(mL-L) xi(mL-L+1)… xi(mL-1)]T，其中i＝1,2,…,M，m＝1,2,…；当第1帧样本采集并保存完毕后，每当完成一帧样本采集，则进入步骤S3；S3：根据每个传声器的第m帧、第m+1帧时域信号及其对应的自适应滤波器系数迭代计算得到第m+1帧时域信号对应的自适应滤波器系数其具体步骤包括：S3.1：计算自适应滤波器系数向量的傅里叶变换向量 h ^ ‾ k 10 ( m ) = F 2 L × 2 L h ^ k T ( m ) 0 1 × L T ]]>其中，k＝1,2,…,M，F2L×2L表示维度为2L×2L的傅里叶矩阵，01×L表示长度为L的所有元素为0的行向量；S3.2：计算传声器采集的时域信号的傅里叶变换向量xi,2L(m)：xi,2L(m)＝F2L×2Lxi,2L(m)其中，i＝1,2,…,M，xi,2L(m)＝[xi(mL-L) xi(mL-L+1) … xi(mL+L-1)]T；S3.3：按照以下步骤依次计算得到时域先验误差信号的傅里叶变换向量eik(m)： e ‾ i k ( m ) = F 2 L × 2 L - 1 e ‾ ‾ i k