[发明专利]一种音频信号中时变正弦干扰的自适应抑制方法

权利要求书

1.一种音频信号中时变正弦干扰的自适应抑制方法，其特征在于，包括以下步骤：

步骤一、输入原始音频信号，对原始音频信号进行AD转换，得到数字化的音频信号；

步骤二、使用自适应陷波算法计算音频信号中正弦干扰信号的初步频率；

步骤三、使用LMS梯度速降算法计算音频信号中正弦干扰信号的幅度、相位和精确频率；

步骤四、使用判决算法判定正弦干扰信号的可信度：如果可信度为真，则转到步骤五；如果可信度为假，则直接输出音频信号；

步骤五、从音频信号中减去计算出的正弦干扰信号，得到抑制干扰后的音频信号并输出。

2.根据权利要求1所述的一种音频信号中时变正弦干扰的自适应抑制方法，其特征在于，所述使用自适应陷波算法计算音频信号中正弦干扰信号的初步频率ω₀(n)，正弦干扰信号的初步频率ω₀(n)与初步频率相关系数k₀(n)的关系为：

|k₀(n)|＝|-cos(ω₀(n))|

初步频率相关系数k₀(n)的递推公式为：

k₀(n)＝-B(n)/A(n)

A(n)＝λA(n-1)+2w(n-1)·w(n-1)

B(n)＝λB(n-1)+w(n-1)·[w(n)+w(n-2)]

式中，遗忘因子λ的值大于等于0.95且小于等于0.99。

3.根据权利要求2所述的一种音频信号中时变正弦干扰的自适应抑制方法，其特征在于，如果初步频率相关系数k₀(n)的绝对值大于1，则使用下述公式对其进行非线性处理：

4.根据权利要求1所述的一种音频信号中时变正弦干扰的自适应抑制方法，其特征在于，使用LMS梯度速降算法计算音频信号中正弦干扰信号的幅度相位和精确频率采用以下递推公式：

P^(n+1)=P^(n)+2TsμAe(n)sinφ^(n)]]>

φ^(n+1)=φ^(n)+Tsω^(n)+2Tsμω′μδ′e(n)P^(n)cosφ^(n)]]>

ω^(n+1)=ω^(n)+2Tsμω′e(n)P^(n)cosφ^(n)]]>

式中，定义误差信号e(n)为：

e(n)=u(n)-η^(n)]]>

式中，μ_A，μ_ω，μ_δ分别为100，1000和0.02。