发明名称 一种正弦信号频率估计的相位匹配方法
摘要 本发明涉及信号处理领域,特别是正弦信号的频率估计方法。本发明的适用对象为正弦信号频率估计,包括以下步骤:首先,利用正弦信号的线性预测性质、自相关计算同频信号的自相关序列,消除初相位不同对正弦信号信息融合的影响;然后,对自相关序列进行时域平均处理,提高信号的信噪比,为改善信号频率估计精度创造有利条件;最后,根据自相关序列的特点生成零初相正弦序列,并利用两者的相关系数构造误差函数,使误差函数最小获得信号频率估计值。本发明涉及的正弦信号的频率估计方法实现简单、抗噪性能好、精度高,消除了多段信号信息融合对初相位的依赖性,提高了信号的频率估计精度。
申请公布号 CN104142425A 申请公布日期 2014.11.12
申请号 CN201410342480.X 申请日期 2014.07.13
申请人 中国人民解放军后勤工程学院 发明人 涂亚庆;沈艳林
分类号 G01R23/02(2006.01)I 主分类号 G01R23/02(2006.01)I
代理机构 代理人
主权项 一种正弦信号频率估计的相位匹配方法,其特征在于:适用对象为M段正弦信号的频率估计。该方法包括以下步骤:第一步:利用<img file="FSA0000106406310000011.GIF" wi="568" 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个值,<img file="FSA0000106406310000012.GIF" wi="653" he="59" />第二步:利用<img file="FSA0000106406310000013.GIF" wi="212" he="113" />对自相关序列u<sub>m</sub>进行时域平均处理,获得时域平均信号序列λ,提高信号的信噪比。第三步:利用下式获得信号频率的粗略估计值ω<sub>1</sub>,<maths num="0001" id="cmaths0001"><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><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>[</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>]</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>]]></math><img file="FSA0000106406310000014.GIF" wi="589" he="265" /></maths>第四步:根据时域平均信号序列λ的特点和频率粗略估计值ω<sub>1</sub>生成零初相正弦序列{cos(kω<sub>1</sub>)}、{sin(kω<sub>1</sub>)},代入下式获得信号频率估计值<img file="FSA0000106406310000016.GIF" wi="68" he="55" />。<maths num="0002" id="cmaths0002"><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><mo>(</mo><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><msub><mi>k&omega;</mi><mn>1</mn></msub><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>sin</mi><mrow><mo>(</mo><msub><mi>k&omega;</mi><mn>1</mn></msub><mo>)</mo></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><mo>[</mo><mi>k </mi><mi>cos</mi><mrow><mo>(</mo><msub><mi>k&omega;</mi><mn>1</mn></msub><mo>)</mo></mrow><mi>sin</mi><mrow><mo>(</mo><msub><mi>k&omega;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>]</mo><mo>}</mo></mrow><mrow><msup><mrow><mo>(</mo><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>sin</mi><mrow><mo>(</mo><msub><mi>k&omega;</mi><mn>1</mn></msub><mo>)</mo></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>[</mo><mi>k </mi><mi>sin</mi><mrow><mo>(</mo><msub><mi>k&omega;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>]</mo></mrow><mn>2</mn></msup></mrow></mfrac></mrow>]]></math><img file="FSA0000106406310000015.GIF" wi="1360" he="280" /></maths>
地址 401311 重庆市沙坪坝区大学城中国人民解放军后勤工程学院