CN102023299B - 一种抑制多径的相关器方法 - Google Patents
一种抑制多径的相关器方法 Download PDFInfo
- Publication number
- CN102023299B CN102023299B CN2010105299011A CN201010529901A CN102023299B CN 102023299 B CN102023299 B CN 102023299B CN 2010105299011 A CN2010105299011 A CN 2010105299011A CN 201010529901 A CN201010529901 A CN 201010529901A CN 102023299 B CN102023299 B CN 102023299B
- Authority
- CN
- China
- Prior art keywords
- multipath
- correlator
- strobe
- sign indicating
- indicating number
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
Description
技术领域
本发明涉及GPS***中的相关器(correlator),尤其是一种抑制多径的相关器方法,属于GPS通信技术领域。
背景技术
GPS卫星发射的信号为将伪码(C/A码)和导航数据码调制在载波上所得到的信号。现代通用的数字GPS接收机主要由GPS天线、变频器、基带信号通道和数据输出等组成,完成信号的接收、捕获、跟踪、解调和定位解算等功能。基本的工作流程为:GPS天线接收所有可见GPS卫星信号,并将微弱电磁波转化为微弱电流,然后通过前置低噪声放大器对GPS信号电流进行放大,接着通过变频转换将高频信号下变频为中频数字信号后送入各个接收的数字基带信道进行处理。这个处理过程为:首先捕获模块通过对所有不确定的C/A码的时延和所有不确定的多普勒频移完成二维搜索,从而找到可视卫星、近似的C/A码延迟和多普勒频移,然后将结果输入到跟踪模块,而跟踪模块主要完成C/A码、相位和多普勒频移的跟踪,从而获得这些参数的精确估计值。采用相位和多普勒频移的估计来复制本地载波频率与接收到的数字信号的载波频率一致,从而去掉载波,得到C/A码和导航数据码,此为载波跟踪环路;利用估计的码延迟通过相关器的运算来得到超前码E、即时码P和滞后码L,从而复制本地C/A码和接收到的数字信号中的C/A码一致,进而去掉C/A码,此为码跟踪环路,其和载波跟踪环路统称为跟踪环路。然后得到导航数据,导航数据中包含了卫星的位置以及卫星信号的发射时间,结合本地时间我们可以得到卫星信号从卫星到接收机的传播时间,乘以光速后得到伪距(这个距离里面有诸如大气时延、接收机时钟误差以及接收机的测量误差等等所引起的距离上的偏差,故称其为伪距。),至此我们得到了卫星位置和伪距,我们至少要4颗卫星的数据,就可以联立解方程组来得到接收机所在的位置,此为解算过程。
实际的GPS测量中,接收机天线接收的信号包括直接来自卫星方向的直达信号和直达信号经过一次或多次在接收天线四周的反射面反射后得到的信号,即多径信号,多径信号和直达信号的传播路径不同,因此有不同的幅度、相位和时延,从而导致合成信号的波形发生扭曲,从而影响定位精度,产生定位误差。对多径抑制方法的研究通常从空域、环路和数据后处理三个方面进行,空域抑制方法着眼于尽量减少接收机天线对多径信号的接收,如调整天线的摆放位置;环路抑制方法则是从环路内部的算法上考虑,或者估计出各个多径的参数并将其从合成信号中除去,或者由改进的相关器设置和码环相关器来直接获得多径被抑制后的结果;数据后处理方法则是通过一些滤波方法来减少多径的影响,如卡尔曼滤波、小波变换等等。它们各有利弊,一个完善的多径抑制策略往往是考虑设计目标后对这些方法的综合运用。在接收机内部环路的数字信号处理过程中采取多种措施来抑制多径效应显得相当现实和有效。
而跟踪环路抑制方法中一种代表性的为Strobe相关器,这是一种软件实现方法。它是一种双Δ相关器,它利用两组早迟码E2、E1、P、L1、L2所构成的一种相关器,它可以看成由一个早迟间隔为d的窄相关相关器和一个早迟间隔为2d的较宽相关器构成,由码环的相关器的输出结果,我们可以得到代入Strobe公式:得到相关结果,此即Strobe相关器。
Strobe相关器对于长多径的抑制效果很好,但是在码相位偏移较小时,误差依然较大,即现有的Strobe对于短多径的抑制效果并不明显。
发明内容
针对现有Strobe相关器对于迟延小于d的短多径无能为力的缺点,提出一种抑制多径的相关器方法。这种方法细化对于同相情况下多径的分类建模,在现有Strobe相关器的基础上增加了一个修正因子,使其在更加广阔的延迟范围内具有多径抑制作用。
本发明采取的技术方案是:一种抑制多径的相关器设计方法,以跟踪环路抑制方法的Strobe相关器为基础,环路中的三个复制码:第一超前码、即时码、第一滞后码,即E1、P、L1,它们两两之间的间隔为1/2码片(即1/2chip),而第二超前码、即时码、第二滞后码,即E2、P、L2,它们两两之间的间隔为1码片,经过载波剥离的输入信号I分别与它们相乘,相乘的结果经过累加和存储器后得到输出值信号功率的第二超前分量、第一超前分量、即时分量、第一滞后分量和第二滞后分量,即将代入Strobe公式:获得相关结果,其特征在于:
式中:Dcorr-in-phase表示本发明带有修正因子的Strobe相关器表达式;Dstrobe表示现有Strobe相关器表达式;Dδ-m-phase表示本发明的修正项。
与现有Strobe相关器相比,带有修正因子的Strobe相关器有如下优点:
1)带有修正因子的Strobe相关器更加深入的发掘了多径信号的特征,进一步对同相多径下的误差进行修正,提供了相关器在短多径下的抑制性能。
2)同相多径情况下,增加的修正因子的Strobe相关器使得相关结果不依赖于对噪声的估计,提供更加精确的码相位估计,从而得到更加精准的定位,这个性能对于城市环境中的定位意义很大。
3)本发明环路的多径抑制方法,同样利用了5个相关器的输出值,和现有Strobe相关器相比,它在提高性能的同时,并没有增加额外的硬件消耗。
附图说明
图1是GPS***中多径信号产生的示意图;
图2是同相多径信号和反相多径信号与直达信号叠加后的示意图;
图3是Strobe相关器中相关三角的示意图;
图4是现有Strobe相关器的实现过程,包括预处理过程和现有Strobe相关器;
图5是存在一条延迟为0.15chip的多径信号的情况下,此时得到的直达信号、同相多径信号和合成信号的示意图,图中表明了合成后的信号分为五段;
图7是本发明相关器的实现过程,包括预处理过程和本发明相关器,这里只画出了同相之路I部分的操作,正交之路Q的操作是一样的。
具体实施方式
参看图1,卫星101发出信号,一部分为直达信号102,直接到达接收天线103;而另一部分则经过地平面或者其他建筑物105反射后到达天线,此即多径信号104,多径信号干扰了定位的精度。
参看图2,根据多径信号的相位和直达信号的关系,可以将多径信号分为同相信号和反相信号两种,图中给出了两种信号以及他们分别和直达合成后的信号的示意图。
参看图3,相关器输出的相关三角如图,而不同的相关间隔可以得到不同的超前分量和滞后分量,图中d1=d,d2=2d,d为一个C/A码码片的间隔,理想情况下,分别得到的E2、E1、P、L1、L2在相关三角上位置如图;若存在多径,他们的位置就会发生变动,就不是两边对称的位置。
现有Strobe相关器,参看图4,包括:输入信号401、本地振荡器402、伪随机噪声码发生器403、乘法器404、累加&存储器405、现有Strobe相关器406。其中401~405为预处理部分,将得到的结果输入到现有Strobe相关器406,从而得到相关结果。
参看图5,当多径延迟Δτm<d时,锁定峰值P位于实际峰值的右边,受多径干扰,不在同一条直线上,但是由于窄相关间隔d较小,始终处于同一条直线上,此时就产生了鉴相误差因此只要能对该误差进行正确估计,并且作为修正项加到现有Strobe相关器表达式上,就能使其辨认出Δτm<d时的误差。
同相多径干扰下,为方便起见,设直达信号的幅值为a,多径信号的幅值为b,延迟为Δτm,则可得到直达信号和多径信号的表达式分别为:
根据在各个区间段直达信号和多径信号的叠加关系,合成信号分为五段,如图5所示:
①:y=a(x+1)
②:y=a(x+1)+b(x+1-Δτm)
③:y=-a(x-1)+b(x+1-Δτm)
④:y=-a(x-1)-b(x-1-Δτm)
⑤:y=-b(x-1-Δτm)
P=PLOS+Pmultipath
=-a(δ-1)+b(δ+1-Δτm)=(-a+b)δ+a+(1-Δτm)b
由此可得到:
其中的2a是未知的常数,需用等价的表达式进行替换。由于存在短多径的情形下IP和在同一直线上,有:
这里,Dcorr-in-phase表示同相情形下的带有修正因子的Strobe相关器表达式;Dstrobe表示现有Strobe相关器表达式;Dδ-in-phase表示同相多径情形下的修正项。
由此可见,带有修正因子的Strobe相关器的输出并不依赖于对噪声n的估计,因此有更好的噪声抑制性能,可以提供更精确的码相位估计。
参看图7,输入信号701、本地振荡器702、伪随机噪声码发生器703、乘法器704、累加&存储器705、模式鉴别器706、现有Strobe相关器707(与图4结构相同)、带有修正因子的Strobe相关器708。环路中的三个复制码E1、P、L1两两之间的间隔为1/2码片(即1/2chip),而E2、P、L2两两之间的间隔为1码片,经过载波剥离的输入信号I701分别与它们相乘,相乘的结果再经过累加&存储器705后得到输出值此时,将得到的值输入到模式鉴别器706中,若满足条件1:
成立,则选择现有Strobe相关器:
Claims (1)
1.一种抑制多径的相关器方法,以跟踪环路抑制方法的Strobe相关器为基础,环路中的三个复制码:第一超前码、即时码、第一滞后码,即E1、P、L1,它们两两之间的间隔为1/2个C/A码码片,第二超前码、即时码、第二滞后码,即E2、P、L2,它们两两之间的间隔为1个C/A码码片,经过载波剥离的输入信号I分别与它们相乘,相乘的结果经过累加和存储器后得到输出值信号功率的第二超前分量、第一超前分量、即时分量、第一滞后分量和第二滞后分量,即 Ip、 将 Ip、 代入Strobe相关公式: 获得相关结果,其特征在于:
式中:Dcorr-in-phase表示带有修正因子的Strobe相关器表达式;Dstrobe表示现有Strobe相关器表达式;Dδ-in-phase表示修正因子,d为一个C/A码码片的长度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010105299011A CN102023299B (zh) | 2010-11-03 | 2010-11-03 | 一种抑制多径的相关器方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010105299011A CN102023299B (zh) | 2010-11-03 | 2010-11-03 | 一种抑制多径的相关器方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102023299A CN102023299A (zh) | 2011-04-20 |
CN102023299B true CN102023299B (zh) | 2012-11-28 |
Family
ID=43864831
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010105299011A Expired - Fee Related CN102023299B (zh) | 2010-11-03 | 2010-11-03 | 一种抑制多径的相关器方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102023299B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183770A (zh) * | 2011-05-20 | 2011-09-14 | 哈尔滨工程大学 | 一种抗多径干扰的gps伪随机码跟踪环路及其抗多径干扰方法 |
CN102338878B (zh) * | 2011-07-19 | 2013-11-06 | 电子科技大学 | 一种多径抑制boc码跟踪方法以及码跟踪环 |
CN109828293B (zh) * | 2017-11-23 | 2020-09-11 | 北京紫光展锐通信技术有限公司 | 一种信号处理方法及其相关设备 |
CN109143276B (zh) * | 2018-02-24 | 2022-09-30 | 上海华测导航技术股份有限公司 | 一种基于cusum算法的组合抗多径鉴相器的设计方法 |
CN108776326B (zh) * | 2018-04-13 | 2022-04-05 | 北京邮电大学 | 一种多径效应抑制方法、装置及设备 |
CN108594273B (zh) * | 2018-05-03 | 2020-03-20 | 中国人民解放军国防科技大学 | 一种基于准时码参考波形设计的载波相位多径抑制方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1219040A (zh) * | 1998-10-28 | 1999-06-09 | 香港大学 | 快速修正误差的切普时钟恢复装置 |
CN101762803A (zh) * | 2008-12-24 | 2010-06-30 | 复旦大学 | 一种应用于全球导航卫星***接收机的差分相关器 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000091954A (ja) * | 1998-09-16 | 2000-03-31 | Oki Electric Ind Co Ltd | 相関演算回路 |
-
2010
- 2010-11-03 CN CN2010105299011A patent/CN102023299B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1219040A (zh) * | 1998-10-28 | 1999-06-09 | 香港大学 | 快速修正误差的切普时钟恢复装置 |
CN101762803A (zh) * | 2008-12-24 | 2010-06-30 | 复旦大学 | 一种应用于全球导航卫星***接收机的差分相关器 |
Non-Patent Citations (3)
Title |
---|
JP特开2000-91954A 2000.03.31 |
一种Strobe相关器及其多径抑制性能研究;纪元法等;《宇航学报》;20070930;第28卷(第5期);第1094-1099页 * |
纪元法等.一种Strobe相关器及其多径抑制性能研究.《宇航学报》.2007,第28卷(第5期),第1094-1099页. |
Also Published As
Publication number | Publication date |
---|---|
CN102023299A (zh) | 2011-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101132191B (zh) | 一种gnss接收机基带信号处理方法 | |
CN102023299B (zh) | 一种抑制多径的相关器方法 | |
CN101646955B (zh) | 时延测量 | |
CN103645483B (zh) | 一种弱信号环境下北斗信号捕获方法 | |
CN102435999B (zh) | Gps接收机基带模块和gps信号捕获及跟踪方法 | |
CN106291614B (zh) | 用于在多路径环境中跟踪卫星无线电导航信号的装置 | |
CN102183770A (zh) | 一种抗多径干扰的gps伪随机码跟踪环路及其抗多径干扰方法 | |
CN107255822A (zh) | 多径环境下gnss接收机信号参数估计方法 | |
CN106842263A (zh) | 基于导航信号空时处理与矢量跟踪相结合的卫星导航方法 | |
CN101520505B (zh) | 全球卫星导航***弱信号下的自适应相干累积联合捕获方法 | |
CN104536016A (zh) | 一种gnss新体制信号捕获装置及方法 | |
CN101216549B (zh) | 中短波扩频导航***距离差观测量提取方法 | |
CN104614740A (zh) | 一种导航信号数据导频联合跟踪方法及装置 | |
CN102565825B (zh) | 接收信号可靠度判定装置、方法及码相位误差算出方法 | |
CN106918822B (zh) | 计算用于解析副载波跟踪模糊度的非模糊鉴别器的gnss接收器 | |
Yang et al. | Post-correlation semi-coherent integration for high-dynamic and weak GPS signal acquisition | |
CN105372678B (zh) | 一种正弦boc调制信号的无模糊跟踪方法 | |
CN103760578B (zh) | 一种gnss卫星导航信号的无模糊跟踪方法 | |
KR101467320B1 (ko) | 균등 세분된 부분상관함수에 기초한 tmboc(6,1,4/33) 신호를 위한 비모호 상관함수 생성 방법, tmboc 신호 추적 장치 및 이를 이용한 위성 항법 신호 수신 시스템 | |
CN102680985B (zh) | 在全球导航卫星***中抑制多径信号的鉴相方法 | |
Khan et al. | Acquisition strategies of GNSS receiver | |
CN201716419U (zh) | 机载多普勒/延迟映射接收机 | |
CN102508266A (zh) | 一种抑制多径干扰及噪声的滤波器方法 | |
CN105527632A (zh) | 全球卫星导航***弱信号处理方法及装置 | |
CN103439715A (zh) | 基于二维跟踪结构的boc信号抗多径方法 |
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 | ||
CB02 | Change of applicant information |
Address after: 211300 Gaochun County of Jiangsu Province Economic Development Zone Branch Center Building Room 405 Applicant after: Southeast University Address before: 210096 Jiangsu city Nanjing Province four pailou No. 2 Applicant before: Southeast University |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121128 Termination date: 20151103 |
|
EXPY | Termination of patent right or utility model |