发明名称 一种正弦信号频率估计的相位匹配方法
摘要 本发明涉及信号处理领域,特别是正弦信号的频率估计方法。本发明的适用对象为正弦信号频率估计,包括以下步骤:首先,利用正弦信号的线性预测性质、自相关计算同频信号的自相关序列,消除初相位不同对正弦信号信息融合的影响;然后,对自相关序列进行时域平均处理,提高信号的信噪比,为改善信号频率估计精度创造有利条件;最后,根据自相关序列的特点生成零初相正弦序列,并利用两者的相关系数构造误差函数,使误差函数最小获得信号频率估计值。本发明涉及的正弦信号的频率估计方法实现简单、抗噪性能好、精度高,消除了多段信号信息融合对初相位的依赖性,提高了信号的频率估计精度。
申请公布号 CN104142425B 申请公布日期 2017.05.03
申请号 CN201410342480.X 申请日期 2014.07.13
申请人 中国人民解放军后勤工程学院 发明人 涂亚庆;沈艳林
分类号 G01R23/02(2006.01)I 主分类号 G01R23/02(2006.01)I
代理机构 代理人
主权项 一种正弦信号频率估计的相位匹配方法,其特征在于:适用对象为M段正弦信号的频率估计;该方法包括以下步骤:第一步:利用<img file="FSB0000160326640000011.GIF" wi="551" he="123" />计算同频信号的自相关序列u<sub>m</sub>,消除初相位不同对同频信号信息融合的影响;其中x<sub>m,n</sub>表示第m段信号第n个点的值,N<sub>m</sub>表示第m段信号长度,m∈[1,M],u<sub>m</sub>表示对应第m段信号的自相关序列,u<sub>m,k</sub>表示自相关序列u<sub>m</sub>的第k个值,k=1,2…p,p≤[(min(N<sub>m</sub>)‑1)/2];第二步:利用<img file="FSB0000160326640000012.GIF" wi="206" he="113" />对自相关序列u<sub>m</sub>进行时域平均处理,获得时域平均信号序列λ,提高信号的信噪比;第三步:利用下式获得信号频率的粗略估计值ω<sub>1</sub>,<maths num="0001"><math><![CDATA[<mrow><msub><mi>&omega;</mi><mn>1</mn></msub><mo>=</mo><msup><mi>cos</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>{</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>2</mn></mrow><mrow><mi>p</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>&lambda;</mi><mi>k</mi></msub><mo>&lsqb;</mo><msub><mi>&lambda;</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>&lambda;</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>&rsqb;</mo></mrow><mrow><mn>2</mn><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>2</mn></mrow><mrow><mi>p</mi><mo>-</mo><mn>1</mn></mrow></munderover><msubsup><mi>&lambda;</mi><mi>k</mi><mn>2</mn></msubsup></mrow></mfrac><mo>}</mo></mrow><mo>;</mo></mrow>]]></math><img file="FSB0000160326640000013.GIF" wi="592" he="263" /></maths>第四步:根据时域平均信号序列λ的特点和频率粗略估计值ω<sub>1</sub>生成零初相正弦序列{cos(kω<sub>1</sub>)}、{sin(kω<sub>1</sub>)},代入下式获得信号频率估计值<img file="FSB0000160326640000015.GIF" wi="42" he="47" />,<maths num="0002"><math><![CDATA[<mrow><mover><mi>&omega;</mi><mo>^</mo></mover><mo>=</mo><msub><mi>&omega;</mi><mn>1</mn></msub><mo>+</mo><mfrac><mrow><mo>{</mo><mrow><mrow><mo>(</mo><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msub><mi>&lambda;</mi><mi>k</mi></msub><mi>cos</mi><mrow><mo>(</mo><mrow><msub><mi>k&omega;</mi><mn>1</mn></msub></mrow><mo>)</mo></mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msub><mi>&lambda;</mi><mi>k</mi></msub><mi>k</mi><mi> </mi><mi>sin</mi><mrow><mo>(</mo><mrow><msub><mi>k&omega;</mi><mn>1</mn></msub></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msubsup><mi>&lambda;</mi><mi>k</mi><mn>2</mn></msubsup><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><mrow><mo>&lsqb;</mo><mrow><mi>k</mi><mi> </mi><mi>cos</mi><mrow><mo>(</mo><mrow><msub><mi>k&omega;</mi><mn>1</mn></msub></mrow><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><mrow><msub><mi>k&omega;</mi><mn>1</mn></msub></mrow><mo>)</mo></mrow></mrow><mo>&rsqb;</mo></mrow></mrow><mo>}</mo></mrow><mrow><msup><mrow><mo>(</mo><mrow><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msub><mi>&lambda;</mi><mi>k</mi></msub><mi>k</mi><mi> </mi><mi>sin</mi><mrow><mo>(</mo><mrow><msub><mi>k&omega;</mi><mn>1</mn></msub></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msubsup><mi>&lambda;</mi><mi>k</mi><mn>2</mn></msubsup><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msup><mrow><mo>&lsqb;</mo><mrow><mi>k</mi><mi> </mi><mi>sin</mi><mrow><mo>(</mo><mrow><msub><mi>k&omega;</mi><mn>1</mn></msub></mrow><mo>)</mo></mrow></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>.</mo></mrow>]]></math><img file="FSB0000160326640000014.GIF" wi="1343" he="280" /></maths>
地址 401311 重庆市沙坪坝区大学城中国人民解放军后勤工程学院