CN103163534B - 一种自适应噪声带宽的载波环路跟踪方法 - Google Patents
一种自适应噪声带宽的载波环路跟踪方法 Download PDFInfo
- Publication number
- CN103163534B CN103163534B CN201310075841.4A CN201310075841A CN103163534B CN 103163534 B CN103163534 B CN 103163534B CN 201310075841 A CN201310075841 A CN 201310075841A CN 103163534 B CN103163534 B CN 103163534B
- Authority
- CN
- China
- Prior art keywords
- loop
- mrow
- carrier
- msup
- noise bandwidth
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 18
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 12
- 230000001427 coherent effect Effects 0.000 claims abstract description 28
- 230000010354 integration Effects 0.000 claims abstract description 13
- 238000005070 sampling Methods 0.000 claims description 11
- 238000009825 accumulation Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000007423 decrease Effects 0.000 claims description 4
- 239000000969 carrier Substances 0.000 claims description 3
- 230000002596 correlated effect Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 3
- 230000009191 jumping Effects 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 230000010355 oscillation Effects 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 2
- 238000001914 filtration Methods 0.000 abstract description 2
- 230000001133 acceleration Effects 0.000 description 3
- 230000000875 corresponding effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000011897 real-time detection Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
本发明提供了一种自适应噪声带宽的载波环路跟踪方法,依次经过载波剥离、伪码精确剥离、相干积分、同相支路相干累加值均方差检测和噪声带宽自适应调整、鉴相、环路滤波和反馈控制载波NCO,将同相通道相干累加值的归一化均方差作为判断门限自适应地调整载波跟踪环路噪声带宽,可以实时地、较好地平衡动态应力变化和环路噪声对环路收敛情况的影响,有效改善环路收敛时间和收敛精度性能。
Description
技术领域
本发明涉及一种自适应噪声带宽载波环路跟踪方法,属于卫星导航领域。
背景技术
全球导航卫星系统(GNSS)能够全天候的在全球范围内提供精确的PVT(位置、速度、时间)信息,具有十分广阔的用途和发展前途。GNSS接收机是接收卫星信号的设备,其利用接收到的卫星信号获得必要的导航和定位信息及观测量,并经过数据处理最终实现导航和定位。在一次完整的导航定位过程中,首先通过捕获操作得到卫星信号的多普勒频移以及码相位的初始值,然后,进入到载波跟踪环路,进一步确定载波频率和扩频码的码相位,同时对卫星信号进行连续准确的跟踪,进而解调出所需的导航数据,并且通过提取星历、伪距计算、卫星位置解算、GNSS接收机位置解算等操作得到最终的位置以及速度信息。
目前,关于载波跟踪环路的研究有很多,但结论大多都是在载波跟踪环路噪声带宽不随输入信号信噪比改变而变化的前提下得到的。在此前提下,提高跟踪灵敏度最直接的方法就是提高输入信号信噪比,即增长相干累加时间。但相干时间的增长又会导致载波跟踪环路应对高动态环境的能力减弱,甚至出现载波跟踪环路失锁的现象。
文献1“Rapid Parallel GPS Signal Acquisition[ION GPS,Salt Lake City,UT,September19-22,2000]”中设计的锁相环采用了两套不同的环路参数,分别对应两个噪声带宽不同的环路滤波器,当接收机信号通道进入跟踪阶段后,先让锁相环采用带宽为110Hz的牵引滤波器,并让其运行30ms,使相位跟踪误差迅速减小,然后再切换到带宽为15Hz的跟踪滤波器。但上述方案只是用跟踪时间作为状态切换条件,不能随着噪声变换情况和动态应力强弱不同而实时、自适应地对环路进行调整,从而无法工作在高动态的环境之下。
发明内容
为了克服现有技术的不足,本发明提供一种自适应噪声带宽载波跟踪环路方法,以同相支路相干累加值IP的统计方差作为判定门限,自适应地调整载波跟踪环路噪声带宽,改善跟踪载波跟踪环路在高动态环境下的跟踪性能。
本发明解决其技术问题所采用的技术方案包括以下步骤:
1)将从跟踪前端处理系统输入的数字中频信号与载波NCO输出的经过映射后的两路本地载波分别混频,形成同相、正交双通道信号;
2)同相、正交双通道信号分别输入码跟踪环前端,在本地超前、即时、滞后支路伪码相关器中做相关,输出为六路相关结果;所述的超前和滞后支路分别比即时支路超前、滞后半码片;
3)对即时支路相关器输出的两路相关结果分别做1ms-10ms的相干积分累加,分别输出两路相干累加值,然后清除积分器的各寄存单元,接着进行下一时段的积分,如此重复不断。
4)同相支路相干累加值均方差检测和噪声带宽自适应调整,详细过程如下:
a)接收同相通道相干累加值IP并实时存储,将存储值进行N点采样后依次进行步骤b)中的归一化均方差计算和步骤c)中的门限判断,在此段时间内,将实时接收的IP值直接传输至载波鉴相器,用PLL环路中最佳噪声带宽作为环路参数参与计算,PLL环路中最佳噪声带宽:
其中,a=Bnω,ω为环路无阻尼振荡频率,Δf为起始频差,c/n0为输入信号载噪比,T为相干积分时间,λL为载波信号波长;
b)设需要统计的IP采样值分别为{Ipi},其中,i=0,1,...,N-1,根据Ipi的正负号分别设定两个集合{I+}和{I-},满足下式
进一步设和为当前统计N点中集合{I+}和{I-}中元素的个数,定义当前N点采样归一化均方差σ2为:
其中,
c)将上一步骤中计算出的归一化均方差与设定的门限进行比较:若在门限范围内,则与捕获后初始跟踪时的状态相同,仍选用由最小环路跟踪误差计算出的最佳噪声带宽;若大于门限上限,则增大噪声带宽,以应对动态应力变化;若小于门限下限,则减小噪声带宽,使环路锁定更加稳定;
5)将经过步骤4)处理的实时存储的同相、正交通道相干累加值输入载波鉴相器,其输出即为载波跟踪环路输入的数字中频信号与本地载波的相位误差或相位误差的函数;
6)载波鉴相器输出的相位误差或相位误差的函数输入到环路滤波器,用步骤4)中得到的自适应调整后的噪声带宽参与计算环路滤波器参数;
7)环路滤波器的输出反馈给载波NCO,以相应增大或减小载波NCO产生的本地载波频率,至此完成一次环路操作,后跳回至步骤1)重复处理,直到载波鉴相器输出的相位误差为零时,判定环路稳定收敛。
本发明的有益效果是:本发明将同相通道相干累加值的归一化均方差作为判断门限自适应地调整载波跟踪环路噪声带宽,可以实时地、较好地平衡动态应力变化和环路噪声对环路收敛情况的影响,有效改善环路收敛时间和收敛精度性能。
附图说明
图1是自适应噪声带宽载波跟踪环路方法流程图;
图2是不同噪声带宽方案IP值归一化均方差对比图;
图3是自适应噪声带宽方案环路调整情况图;
图4是跟踪不同阶段IP值归一化均方差变化情况图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
本发明提出一种自适应噪声带宽的载波跟踪环路方法,如附图1所示,包括以下步骤:
1)载波剥离:将从跟踪前端处理系统输入的数字中频信号与载波NCO(numericalcontrolled oscillator)输出的经过sin/cos映射后的两路本地载波分别混频,形成同相、正交双通道信号。
2)伪码精确剥离:同相、正交双通道信号分别输入码跟踪环前端,在本地超前、即时、滞后支路(超前和滞后支路分别比即时支路超前、滞后半码片)伪码相关器中做相关,输出为六路相关结果。为简化模型,此处只取即时支路相关器及其输出的两路相关结果进行分析,并认为相关后可以精确剥离伪码。
3)相干积分:对即时支路相关器输出的两路相关结果分别做一定时间(通常为1ms-10ms)的相干积分累加,分别输出两路相干累加值,然后清除积分器的各寄存单元,接着进行下一时段的积分,如此重复不断。
4)同相支路相干累加值均方差检测和噪声带宽自适应调整算法为本发明的关键部分,详细过程如下:
a)接收同相通道相干累加值IP并实时存储,将存储值进行N(取值范围根据系统实时处理精度选定)点采样后依次进行步骤b)中的归一化均方差计算和步骤c)中的门限判断。因为采样和归一化均方差计算均不能实时进行,所以在此段时间内,将实时接收的IP值直接传输至载波鉴相器,用文献2中提出的由最小环路跟踪误差计算出的初始最佳噪声带宽作为环路参数参与计算。文献2“动态环境下基于FLL/PLL相结合的载波跟踪[弹箭与制导学报,第24卷第1期,171-173]”提出环路跟踪误差的均方差最小时对应有最佳环路带宽,并求得:
其中,(Bn)PLLoptimum是PLL环路中最佳噪声带宽,a=Bnω,ω为环路无阻尼振荡频率,Δf为起始频差,c/n0为输入信号载噪比,T为相干积分时间,λL为载波信号波长。
b)归一化均方差计算:假设需要统计的IP采样值的点数为N点,采样值分别为{Ipi},其中,i=0,1,...,N-1。根据Ipi的正负号分别设定两个集合{I+}和{I-},满足下式
进一步假设和为当前统计N点中集合{I+}和{I-}中元素的个数,显然,由此,定义当前N点采样归一化均方差σ2为:
其中,
c)门限判断:将上一步骤中计算出的归一化均方差与本步骤中系统已设定的门限(门限取值范围根据系统实时处理精度选定)进行比较:若在门限范围内,则与捕获后初始跟踪时的状态相同,仍选用由最小环路跟踪误差计算出的最佳噪声带宽(如公式1所示);若大于门限上限,则增大噪声带宽,以应对动态应力变化;若小于门限下限,则减小噪声带宽,使环路锁定更加稳定。噪声带宽的每次调整量视系统实时处理精度而定。实际应用中,复杂噪声背景下,也可在本步骤门限判断前增设门限值训练模块,根据噪声变化设定不同门限等级,通过实时检测和反馈训练出最适合当前噪声情况的门限值。
5)鉴相:将经过步骤4)处理的实时存储的同相、正交通道相干累加值输入载波鉴相器,其输出即为载波跟踪环路输入的数字中频信号与本地载波的相位误差(或相位误差的函数)。
6)环路滤波:载波鉴相器输出的相位误差(或相位误差的函数)输入到环路滤波器,用步骤4)中得到的自适应调整后的噪声带宽参与计算环路滤波器参数。
7)反馈控制载波NCO:环路滤波器的输出反馈给载波NCO,以相应增大或减小载波NCO产生的本地载波频率,至此完成一次环路操作,后跳回至步骤1)重复处理,直到载波鉴相器输出的相位误差为零时,判定环路稳定收敛。
本实施例中采用实测GPS信号,输入的数字中频信号为15.58MHz,采样频率为62MHz,经捕获后的多普勒频率偏移控制在±250Hz之内,突发动态应力加速度的变化为2g,g是重力加速度,载波NCO本地正余弦查找表存储深度为512,数据位宽为8位,相干积分时间为1ms,门限范围下限设为1.0e-4,上限设为1.0e-3,需要统计的IP采样值的点数N为25。
附图1中,I1、Q1信号分别为输入的数字中频信号与载波NCO产生的本地载波混频形成的同相、正交双通道信号;I2、Q2为I1、Q1信号与码跟踪环路的即时支路(P支路)相关器做相关运算输出的结果;对两路相关结果I2、Q2分别做相干积分累加和清零操作,对应输出两路相干累加值IP、QP。
附图2所示为按照本实施例参数设置,分别采用固定噪声带宽10Hz、60Hz、固定最佳噪声带宽25Hz和基于IP值归一化均方差门限的自适应噪声带宽时,IP值归一化均方差变化情况。如图可见,随着噪声带宽减小,环路收敛时的IP值归一化均方差逐渐减小;数据采样值5N处(对应跟踪时间125ms处),IP值归一化均方差小于门限下限1.0e-4,该处将噪声带宽由初始25Hz下调至20Hz;采样值20N处(对应跟踪时间500ms)检测到IP值归一化均方差已小于初始门限下限的0.1倍,将噪声带宽再次下调。附图3所示为与附图2相同参数设置下,自适应噪声带宽条件时相应的环路调整情况。图中可以看到随着两次噪声带宽的调整,频率偏移量逐渐趋于稳定,使之更适应当前环路噪声状况,很好地改善了环路收敛精度。
下表所示为突发动态应力加速度为2g(其他参数设置与实施例相同)时,分别选用四种噪声带宽(固定噪声带宽10Hz、固定最佳噪声带宽25Hz、固定噪声带宽60Hz和基于IP值归一化均方差门限的自适应噪声带宽),环路首次收敛时间、环路稳定收敛后IP值均方差和突发动态应力后环路重新收敛情况的对比。分析数据表明,自适应噪声带宽方案较好地平衡了收敛时间和收敛精度,整体改善了环路性能。
表1不同噪声带宽设置下的环路性能对比
Claims (1)
1.一种自适应噪声带宽的载波环路跟踪方法,其特征在于包括下述步骤:
1)将从跟踪前端处理系统输入的数字中频信号与载波NCO输出的经过映射后的两路本地载波分别混频,形成同相、正交双通道信号;
2)同相、正交双通道信号分别输入码跟踪环前端,在本地超前、即时、滞后支路伪码相关器中做相关,输出为六路相关结果;所述的超前和滞后支路分别比即时支路超前、滞后半码片;
3)对即时支路相关器输出的两路相关结果分别做1ms-10ms的相干积分累加,分别输出两路相干累加值,然后清除积分器的各寄存单元,接着进行下一时段的积分,如此重复不断;
4)同相支路相干累加值均方差检测和噪声带宽自适应调整,详细过程如下:
a)接收同相通道相干累加值IP并实时存储,将存储值进行N点采样后依次进行步骤b)中的归一化均方差计算和步骤c)中的门限判断,在此段时间内,将实时接收的IP值直接传输至载波鉴相器,用PLL环路中最佳噪声带宽作为环路参数参与计算,PLL环路中最佳噪声带宽:
其中,a=Bn/ω,ω为环路无阻尼振荡频率,Δf为起始频差,c/n0为输入信号载噪比,T为相干积分时间,λL为载波信号波长;
b)设需要统计的IP采样值分别为{Ipi},其中,i=0,1,...,N-1,根据Ipi的正负号分别设定两个集合{I+}和{I-},满足下式
进一步设和为当前统计N点中集合{I+}和{I-}中元素的个数,定义当前N点采样归一化均方差σ2为:
其中,
c)将上一步骤中计算出的归一化均方差与设定的门限进行比较:若在门限范围内,则与捕获后初始跟踪时的状态相同,仍选用由最小环路跟踪误差计算出的最佳噪声带宽;若大于门限上限,则增大噪声带宽,以应对动态应力变化;若小于门限下限,则减小噪声带宽,使环路锁定更加稳定;
5)将经过步骤4)处理的实时存储的同相、正交通道相干累加值输入载波鉴相器,其输出即为载波跟踪环路输入的数字中频信号与本地载波的相位误差或相位误差的函数;
6)载波鉴相器输出的相位误差或相位误差的函数输入到环路滤波器,用步骤4)中得到的自适应调整后的噪声带宽参与计算环路滤波器参数;
7)环路滤波器的输出反馈给载波NCO,以相应增大或减小载波NCO产生的本地载波频率,至此完成一次环路操作,后跳回至步骤1)重复处理,直到载波鉴相器输出的相位误差为零时,判定环路稳定收敛。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310075841.4A CN103163534B (zh) | 2013-03-11 | 2013-03-11 | 一种自适应噪声带宽的载波环路跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310075841.4A CN103163534B (zh) | 2013-03-11 | 2013-03-11 | 一种自适应噪声带宽的载波环路跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103163534A CN103163534A (zh) | 2013-06-19 |
CN103163534B true CN103163534B (zh) | 2014-09-10 |
Family
ID=48586776
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310075841.4A Active CN103163534B (zh) | 2013-03-11 | 2013-03-11 | 一种自适应噪声带宽的载波环路跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103163534B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103414493B (zh) * | 2013-08-02 | 2015-05-27 | 北京航空航天大学 | 一种通用的非相干直接序列扩频信号跟踪方法 |
CN104283552B (zh) * | 2014-05-06 | 2017-07-28 | 浙江大学 | 一种用于载波提取的锁相环环路带宽动态调整的方法 |
CN104215981B (zh) * | 2014-08-28 | 2017-03-15 | 四川九洲电器集团有限责任公司 | 一种接收机高动态环境下自适应跟踪方法 |
CN105572702A (zh) * | 2015-12-23 | 2016-05-11 | 航天恒星科技有限公司 | 一种滑动窗口环路跟踪方法及装置 |
CN105954777B (zh) * | 2016-04-22 | 2019-05-14 | 北京大学 | 基于相位估计与补偿的载波跟踪方法及其实现装置 |
CN106019324A (zh) * | 2016-05-16 | 2016-10-12 | 河北科技大学 | 高动态卫星导航接收机环路跟踪器及其跟踪方法 |
CN106842253B (zh) * | 2017-03-03 | 2019-04-19 | 中国电子科技集团公司第五十四研究所 | 一种快速收敛的自适应伪码延迟锁定环路 |
CN108931792A (zh) * | 2017-05-22 | 2018-12-04 | 上海华测导航技术股份有限公司 | 一种环路参数动态加权的方法 |
CN108600134B (zh) * | 2018-04-23 | 2020-01-17 | 北京邮电大学 | 一种载波跟踪方法及装置 |
CN109917154B (zh) * | 2019-03-13 | 2021-04-06 | 东南大学 | 一种硅微谐振式加速度计带宽测试系统及方法 |
CN109743075B (zh) * | 2019-03-14 | 2020-12-04 | 哈尔滨工程大学 | 一种三环联动的非同源扩频码跟踪环路同步装置及方法 |
CN111123331B (zh) * | 2019-10-23 | 2023-07-25 | 湖北三江航天险峰电子信息有限公司 | 一种北斗导航伪距监测方法及系统 |
CN112305566B (zh) * | 2020-09-08 | 2023-06-23 | 南京低功耗芯片技术研究院有限公司 | 基于频率补偿的gnss信号跟踪环路失锁检测方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101581775B (zh) * | 2009-06-18 | 2011-06-15 | 哈尔滨工业大学 | 基于四维ukf的高动态gnss载波的开环补偿跟踪方法 |
CN102445698B (zh) * | 2011-11-07 | 2013-06-05 | 东南大学 | 一种gps接收机卫星跟踪失锁的判定方法 |
CN102621561B (zh) * | 2012-04-16 | 2013-05-15 | 南京航空航天大学 | 一种卫星导航接收机环路自调节方法 |
-
2013
- 2013-03-11 CN CN201310075841.4A patent/CN103163534B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN103163534A (zh) | 2013-06-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103163534B (zh) | 一种自适应噪声带宽的载波环路跟踪方法 | |
CN104253774B (zh) | 一种高动态环境下多普勒频偏估计系统和方法 | |
CN101776752B (zh) | 机群链路高动态信号的精密跟踪与测量方法 | |
CN101266292B (zh) | 一种gnss反射信号的频域处理装置和方法 | |
CN102621561B (zh) | 一种卫星导航接收机环路自调节方法 | |
CN102435999B (zh) | Gps接收机基带模块和gps信号捕获及跟踪方法 | |
CN103558612B (zh) | 锁相环与副载波环联合鉴相跟踪环路 | |
CN103954977B (zh) | 一种gnss欺骗干扰感知方法和系统 | |
CN104570016A (zh) | 一种高动态运动载体的北斗信号捕获跟踪接收方法 | |
CN104280746A (zh) | 一种惯性辅助gps的深组合半实物仿真系统 | |
US20140347218A1 (en) | Apparatus and methods for determining status of a tracking loop | |
CN104536016A (zh) | 一种gnss新体制信号捕获装置及方法 | |
CN103926603A (zh) | Gnss接收机极弱信号的跟踪方法 | |
CN103592662A (zh) | 一种gps信号接收机的载波跟踪方法及环路 | |
CN103116038A (zh) | 一种利用卫星接收机载波跟踪环路测量加速度的方法 | |
CN102253396A (zh) | 一种高动态gps载波环跟踪方法 | |
CN105093243A (zh) | 一种基于随机共振算法的gnss载波环路跟踪方法 | |
CN102801671A (zh) | 一种自适应调节参数的载波跟踪装置 | |
CN106526625A (zh) | 一种基于能量的频率鉴别方法和装置 | |
CN102854516A (zh) | 一种gnss接收机中的载噪比估计方法及其系统 | |
Zhang et al. | Tracking GPS signals under ionosphere scintillation conditions | |
Jian-Feng et al. | Adaptively robust phase lock loop for low C/N carrier tracking in a GPS software receiver | |
Chen et al. | High dynamic GPS signal tracking based on UKF and carrier aiding technology | |
CN104199063B (zh) | 一种基于叉积算法的盲鉴频器处理方法 | |
Morton et al. | GNSS receiver signal tracking |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |