发明名称 一种自适应噪声带宽的载波环路跟踪方法
摘要 本发明提供了一种自适应噪声带宽的载波环路跟踪方法,依次经过载波剥离、伪码精确剥离、相干积分、同相支路相干累加值均方差检测和噪声带宽自适应调整、鉴相、环路滤波和反馈控制载波NCO,将同相通道相干累加值的归一化均方差作为判断门限自适应地调整载波跟踪环路噪声带宽,可以实时地、较好地平衡动态应力变化和环路噪声对环路收敛情况的影响,有效改善环路收敛时间和收敛精度性能。
申请公布号 CN103163534A 申请公布日期 2013.06.19
申请号 CN201310075841.4 申请日期 2013.03.11
申请人 西北工业大学 发明人 姚如贵;毕彦博;张兆林;王伶;徐少波
分类号 G01S19/29(2010.01)I 主分类号 G01S19/29(2010.01)I
代理机构 西北工业大学专利中心 61204 代理人 顾潮琪
主权项 1.一种自适应噪声带宽的载波环路跟踪方法,其特征在于包括下述步骤:1)将从跟踪前端处理系统输入的数字中频信号与载波NCO输出的经过映射后的两路本地载波分别混频,形成同相、正交双通道信号;2)同相、正交双通道信号分别输入码跟踪环前端,在本地超前、即时、滞后支路伪码相关器中做相关,输出为六路相关结果;所述的超前和滞后支路分别比即时支路超前、滞后半码片;3)对即时支路相关器输出的两路相关结果分别做1ms-10ms的相干积分累加,分别输出两路相干累加值,然后清除积分器的各寄存单元,接着进行下一时段的积分,如此重复不断。4)同相支路相干累加值均方差检测和噪声带宽自适应调整,详细过程如下:a)接收同相通道相干累加值I<sub>P</sub>并实时存储,将存储值进行N点采样后依次进行步骤b)中的归一化均方差计算和步骤c)中的门限判断,在此段时间内,将实时接收的I<sub>P</sub>值直接传输至载波鉴相器,用PLL环路中最佳噪声带宽作为环路参数参与计算,PLL环路中最佳噪声带宽:<maths num="0001"><![CDATA[<math><mrow><msub><mrow><mo>(</mo><msub><mi>B</mi><mi>n</mi></msub><mo>)</mo></mrow><mrow><mi>PLL</mi><mo>-</mo><mi>optimum</mi></mrow></msub><mo>=</mo><mroot><mfrac><msup><mrow><mo>(</mo><mrow><mo>(</mo><mn>2</mn><mi>n</mi><mo>)</mo></mrow><mfrac><mi>a</mi><mn>3</mn></mfrac><mi>&Delta;f</mi><mo>)</mo></mrow><mn>2</mn></msup><mrow><mfrac><mn>1</mn><mrow><mi>c</mi><mo>/</mo><msub><mi>n</mi><mn>0</mn></msub></mrow></mfrac><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>Tc</mi><mo>/</mo><msub><mi>n</mi><mn>0</mn></msub></mrow></mfrac><mo>)</mo></mrow><msup><mrow><mo>(</mo><mfrac><msub><mi>&lambda;</mi><mi>L</mi></msub><mrow><mn>2</mn><mi>&pi;</mi></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mn>7</mn></mroot></mrow></math>]]></maths>其中,a=B<sub>n</sub>ω,ω为环路无阻尼振荡频率,Δf为起始频差,c/n<sub>0</sub>为输入信号载噪比,T为相干积分时间,λ<sub>L</sub>为载波信号波长;b)设需要统计的I<sub>P</sub>采样值分别为{I<sub>pi</sub>},其中,i=0,1,...,N-1,根据I<sub>pi</sub>的正负号分别设定两个集合{I<sup>+</sup>}和{I<sup>-</sup>},满足下式<maths num="0002"><![CDATA[<math><mrow><msub><mi>I</mi><mi>pi</mi></msub><mo>&Element;</mo><mfenced open='{' close=''><mtable><mtr><mtd><msup><mi>I</mi><mo>+</mo></msup><mo>,</mo></mtd><mtd><msub><mi>I</mi><mi>pi</mi></msub><mo>&GreaterEqual;</mo><mn>0</mn></mtd></mtr><mtr><mtd><msup><mi>I</mi><mo>-</mo></msup><mo>,</mo></mtd><mtd><msub><mi>I</mi><mi>pi</mi></msub><mo>&lt;</mo><mn>0</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>进一步设<img file="FDA00002903081300013.GIF" wi="86" he="77" />和<img file="FDA00002903081300014.GIF" wi="92" he="90" />为当前统计N点中集合{I<sup>+</sup>}和{I<sup>-</sup>}中元素的个数,定义当前N点采样归一化均方差σ<sup>2</sup>为:<maths num="0003"><![CDATA[<math><mrow><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>=</mo><mfrac><mrow><mi>Var</mi><mrow><mo>(</mo><mo>{</mo><msup><mi>I</mi><mo>+</mo></msup><mo>}</mo><mo>)</mo></mrow><mo>&times;</mo><msub><mi>N</mi><msup><mi>I</mi><mo>+</mo></msup></msub><mo>+</mo><mi>Var</mi><mrow><mo>(</mo><mo>{</mo><msup><mi>I</mi><mo>-</mo></msup><mo>}</mo><mo>)</mo></mrow><mo>&times;</mo><msub><mi>N</mi><msup><mi>I</mi><mo>-</mo></msup></msub></mrow><mrow><mi>E</mi><mrow><mo>(</mo><mo>{</mo><msubsup><mi>I</mi><mi>pi</mi><mn>2</mn></msubsup><mo>}</mo><mo>)</mo></mrow><mo>&times;</mo><mi>N</mi></mrow></mfrac></mrow></math>]]></maths>其中,<maths num="0004"><![CDATA[<math><mrow><mi>Var</mi><mrow><mo>(</mo><mo>{</mo><msup><mi>I</mi><mo>+</mo></msup><mo>}</mo><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><msub><mi>N</mi><msup><mi>I</mi><mo>+</mo></msup></msub></mfrac><munderover><mi>&Sigma;</mi><mn>0</mn><mrow><msub><mi>N</mi><msup><mi>I</mi><mo>+</mo></msup></msub><mo>-</mo><mn>1</mn></mrow></munderover><msup><mrow><mo>(</mo><msubsup><mi>I</mi><mi>i</mi><mo>+</mo></msubsup><mo>-</mo><mi>E</mi><mrow><mo>(</mo><mo>{</mo><msup><mi>I</mi><mo>+</mo></msup><mo>}</mo><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><mi>E</mi><mrow><mo>(</mo><mo>{</mo><msup><mi>I</mi><mo>+</mo></msup><mo>}</mo><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><msub><mi>N</mi><msup><mi>I</mi><mo>+</mo></msup></msub></mfrac><munderover><mi>&Sigma;</mi><mn>0</mn><msub><mi>N</mi><msup><mi>I</mi><mo>+</mo></msup></msub></munderover><msubsup><mi>I</mi><mi>i</mi><mo>+</mo></msubsup></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><mi>Var</mi><mrow><mo>(</mo><mo>{</mo><msup><mi>I</mi><mo>-</mo></msup><mo>}</mo><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><msub><mi>N</mi><msup><mi>I</mi><mo>-</mo></msup></msub></mfrac><munderover><mi>&Sigma;</mi><mn>0</mn><mrow><msub><mi>N</mi><msup><mi>I</mi><mo>-</mo></msup></msub><mo>-</mo><mn>1</mn></mrow></munderover><msup><mrow><mo>(</mo><msubsup><mi>I</mi><mi>i</mi><mo>-</mo></msubsup><mo>-</mo><mi>E</mi><mrow><mo>(</mo><mo>{</mo><msup><mi>I</mi><mo>-</mo></msup><mo>}</mo><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><mi>E</mi><mrow><mo>(</mo><mo>{</mo><msup><mi>I</mi><mo>-</mo></msup><mo>}</mo><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><msub><mi>N</mi><msup><mi>I</mi><mo>-</mo></msup></msub></mfrac><munderover><mi>&Sigma;</mi><mn>0</mn><msub><mi>N</mi><msup><mi>I</mi><mo>-</mo></msup></msub></munderover><msubsup><mi>I</mi><mi>i</mi><mo>-</mo></msubsup></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><mi>E</mi><mrow><mo>(</mo><mo>{</mo><msubsup><mi>I</mi><mi>pi</mi><mn>2</mn></msubsup><mo>}</mo><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mn>0</mn><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>I</mi><mi>pi</mi></msub></mrow></math>]]></maths>c)将上一步骤中计算出的归一化均方差与设定的门限进行比较:若在门限范围内,则与捕获后初始跟踪时的状态相同,仍选用由最小环路跟踪误差计算出的最佳噪声带宽;若大于门限上限,则增大噪声带宽,以应对动态应力变化;若小于门限下限,则减小噪声带宽,使环路锁定更加稳定;5)将经过步骤4)处理的实时存储的同相、正交通道相干累加值输入载波鉴相器,其输出即为载波跟踪环路输入的数字中频信号与本地载波的相位误差或相位误差的函数;6)载波鉴相器输出的相位误差或相位误差的函数输入到环路滤波器,用步骤4)中得到的自适应调整后的噪声带宽参与计算环路滤波器参数;7)环路滤波器的输出反馈给载波NCO,以相应增大或减小载波NCO产生的本地载波频率,至此完成一次环路操作,后跳回至步骤1)重复处理,直到载波鉴相器输出的相位误差为零时,判定环路稳定收敛。
地址 710072 陕西省西安市友谊西路127号