CN102629869B - 基于卡尔曼滤波和最小二乘算法的数字延迟锁定环 - Google Patents
基于卡尔曼滤波和最小二乘算法的数字延迟锁定环 Download PDFInfo
- Publication number
- CN102629869B CN102629869B CN201210114242.4A CN201210114242A CN102629869B CN 102629869 B CN102629869 B CN 102629869B CN 201210114242 A CN201210114242 A CN 201210114242A CN 102629869 B CN102629869 B CN 102629869B
- Authority
- CN
- China
- Prior art keywords
- msub
- time delay
- multipath
- mrow
- direct signal
- 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
- 238000001914 filtration Methods 0.000 title claims abstract description 20
- 238000004891 communication Methods 0.000 claims abstract description 5
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 claims abstract description 4
- 230000008030 elimination Effects 0.000 claims description 9
- 238000003379 elimination reaction Methods 0.000 claims description 9
- 238000000034 method Methods 0.000 claims description 3
- 230000002596 correlated effect Effects 0.000 description 6
- 230000000875 corresponding effect Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000001934 delay Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000005311 autocorrelation function Methods 0.000 description 1
- 235000009508 confectionery Nutrition 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000002054 transplantation Methods 0.000 description 1
Landscapes
- Noise Elimination (AREA)
Abstract
本发明为一种基于卡尔曼滤波和最小二乘算法的数字延迟锁定环,包括可编程门阵列单元和数字信号处理单元;其中可编程门阵列单元包括相关器、控制信号产生模块、通信接口,数字信号处理单元包括卡尔曼滤波器、多径估计模块、多径重构模块、多径消除模块、伪码鉴相器、环路滤波器、决策模块、数字控制振荡器;本发明解决了同时估计直接信号和多径信号时间延迟的延迟锁定环对噪声信号敏感、不能抑制时间延迟较短的多径干扰造成的跟踪误差的问题,实现了一种利用卡尔曼滤波对相关输出中的噪声进行滤波、利用最小二乘算法对滤波后的相关输出进行参数估计的数字延迟锁定环。
Description
技术领域
本发明涉及一种延迟锁定环,特别是涉及一种基于卡尔曼滤波和最小二乘算法的数字延迟锁定环。
背景技术
延迟锁定环作为通信***中的伪码跟踪环路,是定位***中接收机的主要跟踪环路之一,延迟锁定环对接收到的伪码信号的跟踪精度直接决定定位定位***的定位精度。多径干扰是高精度定位***的主要误差源,随着人们对定位精度要求的提高,多径误差消除技术受到广泛的关注。随着接收机的数字化,基于数据处理的多径误差消除技术成为研究的热点。目前已有的基于数据处理的延迟锁定环主要有只估计直接信号时间延迟的延迟锁定环、同时估计直接信号和多径信号时间延迟的延迟锁定环。
只估计直接信号时间延迟的延迟锁定环根据参数估计模块估计直接信号的时间延迟来确定环路的同步点,具有计算量小、结构简单、硬件实现容易的特点,是对传统延迟锁定环的鉴相器进行改进后的环路。但是,只估计直接信号时间延迟的延迟锁定环通常对噪声敏感、不能抑制时间延迟较短的多径的影响、不能抑制载波多径误差的缺点,无法满足噪声环境下对精度要求较高的定位***的要求。
同时估计直接信号和多径信号时间延迟的延迟锁定环由参数估计模块估计直接信号和多径信号的时间延迟和幅度参数,根据估计的多径参数重构多径相关输出,并从接收信号的相关输出中减去多径相关输出,得到直接信号的相关输出,再根据直接信号的相关输出进行直接信号时间延迟的估计。同时估计直接信号和多径信号时间延迟的延迟锁定环不仅可以抑制伪码多径误差而且可以抑制载波多径误差,可以有效提高定位***的定位精度。但是,对噪声比较敏感、不能抑制时间延迟较短的多径带来的误差仍然是同时估计直接信号和多径信号时间延迟的延迟锁定环的瓶颈问题。
发明内容
本发明的目的是为了解决同时估计直接信号和多径信号时间延迟的延迟锁定环对噪声信号敏感、不能抑制时间延迟较短的多径干扰造成的跟踪误差的问题,实现了一种利用卡尔曼滤波对相关输出中的噪声进行滤波、利用最小二乘算法对滤波后的相关输出进行参数估计的数字延迟锁定环,即:基于卡尔曼滤波和最小二乘算法的数字延迟锁定环。
为了解决上述技术问题,本发明采用的技术方案是:
基于卡尔曼滤波和最小二乘算法的数字延迟锁定环,包括可编程门阵列单元和数字信号处理单元;其中可编程门阵列单元包括相关器、控制信号产生模块、通信接口,数字信号处理单元包括卡尔曼滤波器、多径估计模块、多径重构模块、多径消除模块、伪码鉴相器、环路滤波器、决策模块、数字控制振荡器;数字信号处理单元读出可编程门阵列单元中相关器的相关输出值,并将相关输出值送给卡尔曼滤波器进行滤波;滤波后的相关输出输入给基于最小二乘的多径估计模块完成对直接信号和多径信号的时间延迟和幅度的估计,其中直接信号时间延迟估计的结果记为第一直接信号时间延迟估计值;参数估计结果作为多径重构模块的输入以实现多径相关输出的重构,多径重构模块将输出的多径相关输出值输入给多径消除模块;多径消除模块从接收到的卡尔曼滤波器的相关输出中减去重构的多径相关输出,所得到的直接信号相关输出为伪码鉴相器的输入,环路滤波器对伪码鉴相器的输出进行滤波后计算出第二直接信号时间延迟估计值;决策模块根据第一直接信号时间延迟估计值与上一时刻环路确定的直接信号时间延迟的偏差、第二直接信号时间延迟估计值与上一时刻环路确定的直接信号时间延迟的偏差分别赋予第一直接信号时间延迟估计值和第二直接信号时间延迟估计值不同的权值,第一直接信号时间延迟估计值和第二直接信号时间延迟估计值的加权和为决策判决模块的输出,即当前时刻估计的时间延迟,该输出通过数字控制振荡器输出到可编程门阵列单元。
所述的相关器支路为26路,其中即时支路相关器1个,早码支路相关器10个,晚码支路相关器15个。
所述的环路滤波器采用一阶环路。
与现有技术相比,本发明的优点是:
1)本发明通过可编程门阵列和数字信号处理器来实现数字延迟锁定环,具有计算速度快、易于算法移植的特点。
2)本发明通过卡尔曼滤波器实现了对相关输出中的高斯噪声的滤波,消除了噪声对相关输出的影响,从而克服了噪声对多径估计精度的影响,提高了噪声环境下多径消除的性能。
3)本发明中采用最小二乘算法进行直接信号和多径信号的时间延迟和幅度估计,可以区分时间延迟较短的多径和直接信号。
4)本发明中的决策判决模块考虑了前馈和反馈结果对直接信号时间延迟的影响,即保持了环路的稳定性又提高了对直接信号时间延迟的估计精度。
附图说明
图1为基于卡尔曼滤波和最小二乘算法的数字延迟锁定环的框图;
图2为卡尔曼滤波和最小二乘算法实现的结构图。
具体实施方式
下面结合附图对本发明作进一步详细说明。
本发明利用接收到的伪码信号与本地产生的伪码信号的相关输出值进行卡尔曼滤波、多径估计、多径重构、多径消除,根据消除多径和噪声影响的直接信号相关输出值进行伪码鉴相,利用时间延迟偏差与鉴相输出之间的关系可以检测出本地时间延迟估计偏差,最后根据决策模块计算的直接信号时间延迟调整本地码,使之与接收到的伪码信号同步。
本发明的工作原理:
接收到的数字伪码信号r(k)为:
式中:c(k)——伪随机码;A0——直接伪码幅度;τ0——直接伪码的传输时间延迟;M——多径信号个数;Ai——第i个多径信号幅度;τi——第i个多径信号相对于直接信号的传输时间延迟;n(k)——白高斯噪声。
图1中,多径重构模块根据估计的多径参数和可以重构出即时支的多径相关输出
多径消除模块从即时支路相关输出的的卡尔曼滤波结果y中减去可得直接信号的相关输出结果为:
同理,可以得到其他支路的直接信号相关输出结果所有支路的直接信号相关输出向量记为 图1中,采用的鉴相器的鉴相函数为:
其中,分别为直接信号的早、晚码相关输出, 环路滤波器对相关输出结果进行滤波之后,可以根据(16)确定相关输出与本地估计时间偏差ε之间的关系,并根据确定k时刻的直接信号的时间延迟τ0,k,1。同时将最小二乘参数估计模块估计的直接信号时间延迟记为τ0,k,2,记
图1中,在决策模块,根据
τ0,k=ω1τ0,k,1+ω2τ0,k,2 (17)
确定k时刻的直接信号时间延迟τ0,k,并根据τ0,k控制数字振荡器调整其输出,将输出结果用于控制信号产生模块产生k时刻的本地信号。通信接口实现数字信号处理器与可编程门阵列之间的通信,将本地产生的信号转换为易于用可编程门阵列实现与接收信号相关的信号。
图2中,当估计的直接信号的时间延迟为时,本地产生的即时伪码信号为r(k)与相关之后得到即时相关支路的相关输出y为:
其中:
R(·)——自相关函数;
n1——即时支路相关输出噪声。
与即时支路间距为d的相关支路的相关输出yd为:
其中,n1,d为与即时支路相关间距为d的相关输出噪声。当d>0时,相关支路为早码支路,反之为晚码支路。本发明中相关器的支路数为26,其中早码支路数为10,晚码支路数为15,即时支路数为1,相邻相关支路的间距为Ts=Tc/10,Tc为一个码片的持续时间。容易证明n1,d也为白高斯噪声,k时刻的相关器输出观测方程可以表示为
yd,k=βyd,k-1+nd,k-1, (4)
β相关输出的系数。
卡尔曼滤波的目的是滤除相关输出中的噪声nd,因此可以建立相关输出yd,k的状态空间模型:
状态模型:yd,k=βyd,k-1+nd,k-1, (5)
观测模型:zd,k=Ayd,k-1+vd,k, (6)
由于多径和相关间距不发生改变时,相关输出也不发生改变,所以β=1。A是根据相关器输出得到的状态转移矩阵,相关输出在估计时间内不发生改变,所以A=1。用卡尔曼滤波器对每一路相关输出进行滤波,过程为
状态预测:
协方差预测:
增益计算:
预测更新:
协方差更新:
其中Qk是***方差,Rk是测量方差,Kk是卡尔曼滤波增益。
图2中,滤波后的26路相关输出结果作为最小二乘估计模块的输入。定义相关输出向量为y=[y1,…,y26]T∈R26×1,y是相关输出向量,其中y(1),y(2),…,y(10)分别为第1,第2,…,第10路早码相关器的滤波后输出结果,且第1路相关器与即时支路相关器的间距为10Ts,第2路相关器与即时支路相关器的间距为9Ts,以此类推,第10路相关器与即时支路相关器的间距为Ts;y(11)为即时支路的相关输出,对应的本地估计时间延迟为y(12),…,y(26)分别为第1路至第15路晚码支路的相关输出,且与即时支路的相关间距分别为-Ts至-15Ts。根据相关输出与信号参数A0,Ai以及τ0,τi之间的关系,可以将理想的相关输出结果表示为
y=Gh, (12)
其中G是相关输出矩阵,元素为gij=R((i-j)Ts/Tc),i,j=1,…,26,h是26×1维向量。对于最小二乘算法,向量h可以通过最小化||y-Gh||2的2范数来计算。h的最小二乘解表示为:
hLS=(GHG)-1GHy. (13)
所以求解多径参数的问题转化为寻找hLS中非零元素及其对应的索引1的问题。如果在hLS的第1个位置存在直接信号或者多径信号,那么对应的实际时间延迟为Loc+[1-(11)],此时时间延迟的单位为Ts;否则,不存在多径信号。其中hLS中最大元素对应的位置为直接信号时间延迟对应的位置。
Claims (3)
1.基于卡尔曼滤波和最小二乘算法的数字延迟锁定环,包括可编程门阵列单元和数字信号处理单元;其特征在于:可编程门阵列单元包括相关器、控制信号产生模块、通信接口,数字信号处理单元包括卡尔曼滤波器、多径估计模块、多径重构模块、多径消除模块、伪码鉴相器、环路滤波器、决策模块、数字控制振荡器;数字信号处理单元读出可编程门阵列单元中相关器的相关输出值,并将相关输出值送给卡尔曼滤波器进行滤波;滤波后的相关输出输入给基于最小二乘的多径估计模块完成对直接信号和多径信号的时间延迟和幅度的估计,其中直接信号时间延迟估计的结果记为第一直接信号时间延迟估计值;参数估计结果作为多径重构模块的输入以实现多径相关输出的重构,多径重构模块将输出的多径相关输出值输入给多径消除模块;多径消除模块从接收到的卡尔曼滤波器的相关输出中减去重构的多径相关输出,所得到的直接信号相关输出为伪码鉴相器的输入,环路滤波器对伪码鉴相器的输出进行滤波后计算出第二直接信号时间延迟估计值;
决策模块根据第一直接信号时间延迟估计值与上一时刻环路确定的直接信号时间延迟的偏差、第二直接信号时间延迟估计值与上一时刻环路确定的直接信号时间延迟的偏差分别赋予第一直接信号时间延迟估计值和第二直接信号时间延迟估计值不同的权值,第一直接信号时间延迟估计值和第二直接信号时间延迟估计值的加权和为决策判决模块的输出,即当前时刻估计的时间延迟,该输出通过数字控制振荡器输出到可编程门阵列单元;
其中,在决策模块,根据τ0,k=ω1τ0,k,1+ω2τ0,k,2确定k时刻的直接信号时间延迟τ0,k,并根据τ0,k控制数字振荡器调整其输出,将输出结果用于控制信号产生模块产生k时刻的本地信号,
其中,τ0,k,1为k时刻根据鉴相器输出计算得到的直接信号时间延迟,记为第一直接信号时间延迟;τ0,k,2为根据最小二乘参数估计模块估计的直接信号时间延迟,记为第二直接信号时间延迟; 且 为k-1时刻直接信号时间延迟的估计值。
2.如权利要求1所述的基于卡尔曼滤波和最小二乘算法的数字延迟锁定环,其特征在于:所述的相关器支路为26路,其中即时支路相关器1个,早码支路相关器10个,晚码支路相关器15个。
3.如权利要求1或2所述的基于卡尔曼滤波和最小二乘算法的数字延迟锁定环,其特征在于:所述的环路滤波器采用一阶环路。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210114242.4A CN102629869B (zh) | 2012-04-18 | 2012-04-18 | 基于卡尔曼滤波和最小二乘算法的数字延迟锁定环 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210114242.4A CN102629869B (zh) | 2012-04-18 | 2012-04-18 | 基于卡尔曼滤波和最小二乘算法的数字延迟锁定环 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102629869A CN102629869A (zh) | 2012-08-08 |
CN102629869B true CN102629869B (zh) | 2015-08-26 |
Family
ID=46588041
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210114242.4A Active CN102629869B (zh) | 2012-04-18 | 2012-04-18 | 基于卡尔曼滤波和最小二乘算法的数字延迟锁定环 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102629869B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103760576B (zh) * | 2013-12-02 | 2016-03-02 | 北京理工大学 | 一种基于相关峰的可变间距采样多径检测方法 |
US9876581B2 (en) * | 2016-02-10 | 2018-01-23 | Inphi Corporation | Circuit for multi-path interference mitigation in an optical communication system |
CN106291604B (zh) * | 2016-08-02 | 2018-06-26 | 桂林电子科技大学 | 卫星导航信号接收机的改进码跟踪方法及环路 |
CN106842240B (zh) * | 2017-03-17 | 2019-04-30 | 太原理工大学 | 基于最小误差熵和ε等级差分进化的多径估计器 |
CN107682235A (zh) * | 2017-11-24 | 2018-02-09 | 重庆邮电大学 | 一种基于卡尔曼滤波的传输时延预测数据调度方法 |
DE102018212369A1 (de) * | 2018-07-25 | 2020-01-30 | Robert Bosch Gmbh | Verfahren zur Überwachung einer Energieversorgung in einem Kraftfahrzeug |
CN110505171B (zh) * | 2019-07-12 | 2022-03-25 | 四川安迪科技实业有限公司 | 一种无数据辅助频偏估计方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5398034A (en) * | 1993-03-29 | 1995-03-14 | Stanford Telecommunications, Inc. | Vector delay lock loop processing of radiolocation transmitter signals |
CN1311575A (zh) * | 2000-02-28 | 2001-09-05 | 三菱电机株式会社 | 副载波频率信号解调装置 |
CN101084451A (zh) * | 2004-11-24 | 2007-12-05 | SiRF技术公司 | 基准振荡器频率校正*** |
-
2012
- 2012-04-18 CN CN201210114242.4A patent/CN102629869B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5398034A (en) * | 1993-03-29 | 1995-03-14 | Stanford Telecommunications, Inc. | Vector delay lock loop processing of radiolocation transmitter signals |
CN1311575A (zh) * | 2000-02-28 | 2001-09-05 | 三菱电机株式会社 | 副载波频率信号解调装置 |
CN101084451A (zh) * | 2004-11-24 | 2007-12-05 | SiRF技术公司 | 基准振荡器频率校正*** |
Also Published As
Publication number | Publication date |
---|---|
CN102629869A (zh) | 2012-08-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102629869B (zh) | 基于卡尔曼滤波和最小二乘算法的数字延迟锁定环 | |
CN113238261B (zh) | 低轨卫星扩频通信体制信号捕获跟踪*** | |
KR101173925B1 (ko) | 시간 지연 측정 | |
CN101826889B (zh) | 一种二进制偏移载波信号跟踪环路 | |
CN105917585B (zh) | 用于接收复合信号的方法和接收器 | |
CN103592662B (zh) | 一种gps信号接收机的载波跟踪方法及环路 | |
CN102435999B (zh) | Gps接收机基带模块和gps信号捕获及跟踪方法 | |
Borio et al. | Doppler measurement accuracy in standard and high-sensitivity global navigation satellite system receivers | |
CN110138698B (zh) | 高阶调制线性混合信号频偏初相联合估计方法及装置 | |
CN102183770A (zh) | 一种抗多径干扰的gps伪随机码跟踪环路及其抗多径干扰方法 | |
JP2018529107A (ja) | 飛行時間検出のシステム及び方法 | |
CN103338024B (zh) | 天线组阵中时延的互补卡尔曼滤波装置与方法 | |
CN111399006B (zh) | 一种高灵敏度gnss载波跟踪环路优化方法 | |
CN103439719B (zh) | 一种gps l2c信号跟踪方法 | |
CN112910533B (zh) | 一种具有并行结构的宽带信号组阵*** | |
Morton et al. | GNSS receiver signal tracking | |
CN103873107B (zh) | 通信设备 | |
Vila-Valls et al. | An interactive multiple model approach for robust GNSS carrier phase tracking under scintillation conditions | |
Cortés et al. | Low-complexity adaptive direct-state Kalman filter for robust GNSS carrier tracking | |
CN115865127B (zh) | 一种直扩信号的参数估计和解调方法 | |
CN114578387A (zh) | 多径干扰消除方法及装置、估计方法及装置、接收机 | |
Jianfeng et al. | Low C/N0 carrier tracking loop based on optimal estimation algorithm in GPS software receivers | |
CN113050131A (zh) | 一种基于预处理fft和栅栏效应校正的捕获方法 | |
Chang | Anti-Multipath Filter with Multiple Correlators in GNSS Receviers | |
Sokolova et al. | Doppler rate measurements in standard and high sensitivity GPS receivers: Theoretical analysis and comparison |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C53 | Correction of patent for invention or patent application | ||
CB03 | Change of inventor or designer information |
Inventor after: Gan Minggang Inventor after: Chen Jie Inventor after: Cheng Lan Inventor before: Chen Jie Inventor before: Gan Minggang Inventor before: Cheng Lan |
|
COR | Change of bibliographic data |
Free format text: CORRECT: INVENTOR; FROM: CHEN JIE GAN MINGGANG CHENG LAN TO: GAN MINGGANG CHEN JIE CHENG LAN |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |