频率失配条件下的自适应无偏相位差估计方法与流程

技术特征：

1.频率失配条件下的自适应无偏相位差估计方法，适用对象为频率已知的两路同频余弦信号；

该方法包括以下步骤：

第一步：利用获得的信号频率值ω，构建频率信号cosωn和sinωn，算法步长μ＝0.01；

第二步：计算第一路余弦信号x₁(n)的估计值s₁’(n)＝a₁’(n)cosωn-b₁’(n)sinωn，其中a₁’(n)，b₁’(n)分别为x₁(n)中参数a₁和b₁的估计值；

第三步：计算第二路余弦信号x₂(n)的估计值s₂’(n)＝a₂’(n)cosωn-b₂’(n)sinωn，其中a₂’(n)，b₂’(n)分别为x₂(n)中参数a₂和b₂的估计值；

第三步：构建误差信号e₁(n)＝x₁(n)-s₁’(n)和e₂(n)＝x₂(n)-s₂’(n)；

第四步：利用自适应迭代计算方式分别估计a₁’(n+1)，b₁’(n+1)，a₂’(n+1)和b₂’(n+1)，迭代计算方法为：

a₁’(n+1)＝a₁’(n)+μe₁(n)cosωn；

b₁’(n+1)＝b₁’(n)-μe₁(n)sinωn；

a₂’(n+1)＝a₂’(n)+μe₂(n)cosωn；

b₂’(n+1)＝b₂’(n)-μe₂(n)sinωn；

第五步：根据计算的a₁’(n+1)，b₁’(n+1)，a₂’(n+1)和b₂’(n+1)值，计算频率已知的两路同频余弦信号的相位差，则相位差为：

$\mathrm{\Δ}\mathrm{\θ}(n+1)=arctan\left(\frac{{a_2}^{\′}(n+1){b_1}^{\′}(n+1)-{a_1}^{\′}(n+1){b_2}^{\′}(n+1)}{{a_1}^{\′}(n+1){a_2}^{\′}(n+1)+{b_1}^{\′}(n+1){b_2}^{\′}(n+1)}\right).$