CN105510932B - 一种码环鉴别器及短多径抑制方法 - Google Patents
一种码环鉴别器及短多径抑制方法 Download PDFInfo
- Publication number
- CN105510932B CN105510932B CN201510996958.5A CN201510996958A CN105510932B CN 105510932 B CN105510932 B CN 105510932B CN 201510996958 A CN201510996958 A CN 201510996958A CN 105510932 B CN105510932 B CN 105510932B
- Authority
- CN
- China
- Prior art keywords
- branch
- pseudo code
- pseudo
- lagging
- leading
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/22—Multipath-related issues
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明实施例提供了这一种码环鉴别器及短多径抑制方法,其中,码环鉴别器包括:载波剥离运算器、伪码发生器、伪码环路滤波器、伪码鉴相器、十个积分清零器和十个卷积运算器;十个卷积运算器中的每一卷积运算器的输入端均分别连接至载波剥离运算器的输出端和伪码发生器的输出端,十个卷积运算器输出端一一对应连接至十个积分清零器的输入端;十个积分清零器的输出端连接至伪码鉴相器的输入端,伪码鉴相器的输出端连接至伪码环路滤波器的输入端,伪码环路滤波器输出端连接至伪码发生器。应用本发明实施例,可以对短多径信号进行精确地码相位跟踪,从而有效地抑制短多径效应,获得高精度的室内定位服务。
Description
技术领域
本发明涉及室内定位技术领域,特别是涉及一种码环鉴别器及短多径抑制方法。
背景技术
目前,卫星定位***已被广泛的应用,主要的应用场景分为室外环境场景和室内环境场景。
其中,由于室外环境中的定位信号的多径传播方式大多是长多径传播和中多径传播,而目前提出的抑制多径的技术对抑制长多径和中多径非常有效,因而在室外环境下能够获得比较准确的定位精度。而由于在室内环境中定位信号的多径传播方式大多是近距离的、短时延的短多径传播,但是现有的抑制多径传播技术并不能有效的抑制短多径效应,这样使得定位技术在室内环境下的定位精度并不高。
因此,如何设计一种新的码环鉴别器及应用于该新的码环鉴别器的短多径抑制算法,以实现对短多径信号进行精确地码相位跟踪,从而有效地抑制短多径效应,提高室内定位服务的定位精度,成为亟待解决的技术问题。
发明内容
本发明实施例的目的在于提供一种码环鉴别器及短多径抑制方法,实现了对短多径信号进行精确地码相位跟踪,有效地抑制短多径效应,提高了室内定位服务的定位精度。具体技术方案如下:
第一方面,本发明实施例提供了一种码环鉴别器,所述码环鉴别器包括:载波剥离运算器、伪码发生器、伪码环路滤波器、伪码鉴相器、十个积分清零器和十个卷积运算器;
所述十个卷积运算器中的每一卷积运算器的输入端均分别连接至所述载波剥离运算器的输出端和所述伪码发生器的输出端,所述十个卷积运算器输出端一一对应连接至所述十个积分清零器的输入端;
所述十个积分清零器的输出端连接至所述伪码鉴相器的输入端,所述伪码鉴相器的输出端连接至所述伪码环路滤波器的输入端,所述伪码环路滤波器输出端连接至所述伪码发生器。
可选地,所述十个卷积运算器中的每一卷积运算器的输入端均分别连接至所述载波剥离运算器的输出端和所述伪码发生器的输出端,包括:
所述十个卷积运算器中的五个卷积运算器的输入端均分别连接至所述载波剥离运算器的同相支路输出端,且所述伪码发生器的第一超前支路、第二超前支路、即时支路、第一滞后支路和第二滞后支路的输出端分别与所述五个卷积运算器的输入端唯一对应连接;
所述十个卷积运算器中的剩余五个卷积运算器的输入端均分别连接至所述载波剥离运算器的正交支路输出端,且所述伪码发生器的所述第一超前支路、所述第二超前支路、所述即时支路、所述第一滞后支路和所述第二滞后支路的输出端分别与所述剩余五个卷积运算器的输入端唯一对应连接。
可选地,所述伪码发生器的所述第一超前支路超前即时支路d个码片,所述第二超前支路超前所述第一超前支路d个码片;
所述伪码发生器的所述第一滞后支路滞后即时支路d个码片,所述第二滞后支路滞后所述第一滞后支路d个码片,其中,0.2≤d≤0.3。
第二方面,本发明实施例提供了一种短多径抑制方法,应用于上述任一项所述的码环鉴别器,所述方法包括:
将所述码环鉴别器中的伪码发生器产生的与第一超前支路对应的第一超前支路伪码序列、与第二超前支路对应的第二超前支路伪码序列、与第一滞后支路对应的第一滞后支路伪码序列和与第二滞后支路对应的第二滞后支路伪码序列,分别与在载波剥离运算器中进行载波剥离后的中频数字信号在八路并行的卷积运算器中进行解扩,以得到解扩后的八路信号,其中,所述中频数字信号在经载波剥离运算器进行载波剥离后得到在同相支路进行载波剥离后的第一剥离信号和在正交支路进行载波剥离后的第二剥离信号;
在所述八个积分清零器中对所述解扩后的八路信号进行相关运算,以得到解扩时与同相支路所相关的第一超前支路、第二超前支路、第一滞后支路和第二滞后支路所对应的第一类幅值,和解扩时与正交支路所相关的第一超前支路、第二超前支路、第一滞后支路和第二滞后支路所对应的第二类幅值;
基于所得到的第一类幅值和第二类幅值,计算所述伪码发生器中的各个支路经过积分清零器后的自相关幅值;
基于所述伪码发生器中的各个支路的自相关幅值和预设的非相干超前混合滞后公式,计算由伪码发生器产生的即时支路伪码序列与接收码的相位差,其中,所述第一剥离信号和所述第二剥离信号中携带有接收码;
根据所述相位差,调节所述伪码发生器产生的即时伪码序列的相位,以使即时支路伪码序列与接收码在相位上保持一致。
可选地,计算伪码发生器中的任一支路的自相关幅值的预设自相关幅值计算公式具体为:
其中,所述IE是支路所对应的第一类幅值,所述QE是所述支路所对应的第二类幅值,所述E为所述支路经过积分清零器计算后得到的自相关幅值。
可选地,所述预先构建的非相干超前混合滞后公式可为:
其中,所述δ为由伪码发生器产生的即时伪码序列与接收码之间的相位差,所述t为跟踪时间,b是程序调试获取的经验参数,所述E1为第二超前支路的自相关幅值、E为第一超前支路的自相关幅值,L1为第二滞后支路的自相关幅值、L为第一滞后支路的自相关幅值。
本发明实施例提供的一种码环鉴别器及短多径抑制方法,将在主流码环鉴别器中的六个积分清零器增加至十个,并在伪码发生器中对应增加一对第二超前支路和一对第二滞后支路;该种码环鉴别器再通过结合本申请所提出的非相干超前混合滞后算法,实现了对短多径信号进行精确地码相位跟踪,有效地抑制短多径效应,从而提高了室内定位服务的定位精度。当然,实施本发明的任一产品或方法必不一定需要同时达到以上所述的所有优点。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种码环鉴别器的原理图;
图2为本发明实施例提供的一种短多径抑制方法的示意流程图;
图3为本发明实施例的直达相关函数与同相情况下与一条多径相关函数合成的相关函数的示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为了解决现有技术问题,本发明实施例提供了一种码环鉴别器及应用于该码环鉴别器的短多径抑制方法。
首先,本发明实施例提供了一种码环鉴别器,如图1所示,该码环鉴别器包括:
载波剥离运算器11、伪码发生器12、伪码环路滤波器13、伪码鉴相器14、十个积分清零器15a-15j和十个卷积运算器16a-16j;
该十个卷积运算器16a-16j中的每一卷积运算器的输入端均分别连接至该载波剥离运算器11的输出端和该伪码发生器12的输出端,该十个卷积运算器16a-16j输出端一一对应连接至该十个积分清零器15a-15j的输入端;
该十个积分清零器15a-15j的输出端连接至该伪码鉴相器14的输入端,该伪码鉴相器14的输出端连接至该伪码环路滤波器13的输入端,该伪码环路滤波器13输出端连接至该伪码发生器12。
其中,值得注意的是,该十个卷积运算器16a-16j中的每个卷积运算器的输入端均连接至载波剥离运算器11的输出端,且十个卷积运算器16a-16j中的每个卷积运算器的输入端均和伪码发生器14的输出端相连接;另外,该十个卷积运算器16a-16j输出端一一对应连接至该十个积分清零器15a-15j的输入端所表达的具体含义为:卷积运算器15a的输出端连接至积分清零器16a的输入端,卷积运算器15b的输出端连接至积分清零器16b的输入端,卷积运算器15c的输出端连接至积分清零器16c的输入端,……,卷积运算器15j的输出端连接至积分清零器16j的输入端。
此外,该载波剥离运算器11用于将射频前端处理得到的数字中频信号sIF,I和sIF,Q进行载波剥离,从而得到在相关支路i上的第一剥离信号i(n)和在正交支路q上的第二剥离信号q(n);该伪码发生器12用于产生与第一超前支路E对应的第一超前支路伪码序列、与第二超前支路E1对应的第二超前支路伪码序列、与即时支路P对应的即时支路伪码序列、与第一滞后支路L对应的第一滞后支路伪码序列和与第二滞后支路L1对应的第二滞后支路伪码序列;该十个卷积运算器16a-16j中卷积运算器16a至16e用于实现各个支路的伪码序列与第一剥离信号i(n)进行载波相关,依次得到载波相关函数iE1(n)、iE(n)、iP(n)、iL(n)和iL1(n);卷积运算器中16f至16j用于实现各个支路的伪码序列与第二剥离信号q(n)进行载波相关,依次得到载波相关函数qL1(n)、qL(n)、qP(n)、qE(n)和qE1(n);该十个积分清零器15a-15j中的积分清零器15a至15e用于跟卷积运算器16a至16e一一对应连接,且积分清零器在清空上一次积分值后对卷积运算器16a至16e的输出值进行相干积分计算,依次得到如图1所示的相干积分值IE1、IE、IP、IL和IL1;积分清零器15f至15j用于跟卷积运算器16f至16j一一对应连接,且积分清零器在清空上一次积分值后对卷积运算器16f至16j的输出值进行相干积分计算,依次得到相干积分值QL1、QL、QP、QE和QE1,以实现各个支路伪码序列与接收码相关;该伪码鉴相器14用于求解各个支路伪码序列对应的自相关幅值且用于鉴定即时伪码序列与接收码的相位差;该伪码环路滤波器用于抑制噪声及高频分量,并且控制着环路相位校正的速度与精度。
本发明实施例提供的一种码环鉴别器及短多径抑制方法,将在主流码环鉴别器中的六个积分清零器增加至十个,并在伪码发生器中对应增加一对第二超前支路和一对第二滞后支路;该种码环鉴别器再通过结合本申请所提出的非相干超前混合滞后算法,实现了对短多径信号进行精确地码相位跟踪,有效地抑制短多径效应,从而提高了室内定位服务的定位精度。
更进一步地,在本发明实施例中,该十个卷积运算器中的每一卷积运算器的输入端均分别连接至该载波剥离运算器的输出端和该伪码发生器的输出端,包括:
该十个卷积运算器中的五个卷积运算器的输入端均分别连接至该载波剥离运算器的同相支路输出端,且该伪码发生器的第一超前支路、第二超前支路、即时支路、第一滞后支路和第二滞后支路的输出端分别与该五个卷积运算器的输入端唯一对应连接;
该十个卷积运算器中的剩余五个卷积运算器的输入端均分别连接至该载波剥离运算器的正交支路输出端,且该伪码发生器的该第一超前支路、该第二超前支路、该即时支路、该第一滞后支路和该第二滞后支路的输出端分别与该剩余五个卷积运算器的输入端唯一对应连接。
其中,值得注意的是,该十个卷积运算器16a-16j中的任意五个卷积运算器的输入端均连接至载波剥离运算器11的同相支路i的输出端,且该任意五个卷积运算器的输入端一一和伪码发生器14的第一超前支路E、第二超前支路E1、即时支路P、第一滞后支路L和第二滞后支路L1的输出端相连接;另外,该十个卷积运算器16a-16j中的剩余五个卷积运算器的输入端均连接至载波剥离运算器11的正交支路q的输出端,且该剩余五个卷积运算器的输入端一一和伪码发生器14的第一超前支路E、第二超前支路E1、即时支路P、第一滞后支路L和第二滞后支路L1的输出端相连接。
更进一步地,在本发明实施例中,该伪码发生器的该第一超前支路超前即时支路d个码片,该第二超前支路超前该第一超前支路d个码片;
该伪码发生器的该第一滞后支路滞后即时支路d个码片,该第二滞后支路滞后该第一滞后支路d个码片,其中,0.2≤d≤0.3。
可以理解的是,在该种实现方式中,在现有主流六个相关器,即六个积分清零器的基础上,即存在第一超前支路、即时支路、第一滞后支路,且第一超前支路超前即时支路d个码片,第一滞后滞后滞后即时支路d个码片的基础上,增加第二超前支路和第二滞后支路,且设置第二超前支路超前第一超前支路d个码片,第二滞后支路超前第一滞后支路d个码片,且设置码片的大小为0.2≤d≤0.3,可以将伪码环路滤波器的带宽设计得更窄,进而减少环路噪声。
相应于上述码环鉴别器,本发明实施例还提供了一种适用于上述码环鉴别器的短多径抑制方法,如图2所示,该适用于上述码环鉴别器的短多径抑制方法可以包括如下步骤:
S101:将该码环鉴别器中的伪码发生器产生的与第一超前支路对应的第一超前支路伪码序列、与第二超前支路对应的第二超前支路伪码序列、与第一滞后支路对应的第一滞后支路伪码序列和与第二滞后支路对应的第二滞后支路伪码序列,分别与在载波剥离运算器中进行载波剥离后的中频数字信号在八路并行的卷积运算器中进行解扩,以得到解扩后的八路信号,其中,该中频数字信号在经载波剥离运算器进行载波剥离后得到在同相支路进行载波剥离后的第一剥离信号和在正交支路进行载波剥离后的第二剥离信号;
S102:在八个积分清零器中对该解扩后的八路信号进行相关运算,以得到解扩时与同相支路所相关的第一超前支路、第二超前支路、第一滞后支路和第二滞后支路所对应的第一类幅值,和解扩时与正交支路所相关的第一超前支路、第二超前支路、第一滞后支路和第二滞后支路所对应的第二类幅值;
S103:基于所得到的第一类幅值和第二类幅值,计算该伪码发生器中的各个支路经过伪码鉴相器后的自相关幅值;
S104:基于该伪码发生器中的各个支路的自相关幅值和预设的非相干超前混合滞后公式,计算由伪码发生器产生的即时支路伪码序列与接收码的相位差,其中,所述第一剥离信号和所述第二剥离信号中携带有接收码;
S105:根据该相位差,调节该伪码发生器产生的即时伪码序列的相位,以使即时支路伪码序列与接收码在相位上保持一致。
其中,伪码发生器在第一超前支路产生第一超前支路伪码序列、在第二超前支路产生第二超前支路伪码序列、在第一滞后支路产生第一滞后支路伪码序列及在第二滞后支路产生第二滞后支路伪码序列;该四路伪码序列分别与同相支路、正交支路进行载波相关,得到对应的八路信号;将每路信号在一个唯一对应的积分清零器中进行相干积分,得到每路信号所对应的相干积分,再将四路伪码序列中的每路伪码序列在同相支路和正交支路上的相干结果进行平方后再开方,得到各个伪码序列的自相关幅值;将得到的各个伪码序列的自相关幅值相应代入预设的非相干超前混合滞后公式,计算此时即时伪码序列和接收码的相位差;根据相位差调整伪码发生器产生的伪码序列的相位,以实现即时伪码序列与即时码实现同步,从而实现准确的获取卫星***发送的数据信息,从而实现准确定位。
值得强调的是,具体的定位算法是现有技术,在实现接收信号与发射信号同步的前提下实现精确的定位算法,在此不做详述。
本发明实施例提供的一种新的码环鉴别器及短多径抑制方法,将在主流码环鉴别器中的六个积分清零器增加至十个,并在伪码发生器中对应增加一对第二超前支路和一对第二滞后支路;该种新的码环鉴别器再通过结合本申请所提出的非相干超前混合滞后算法,实现了对短多径信号进行精确地码相位跟踪,有效地抑制短多径效应,从而提高了室内定位服务的定位精度。
更进一步地,计算伪码发生器中的任一支路的自相关幅值的预设自相关幅值计算公式具体为:
其中,该IE是支路所对应的第一类幅值,该QE是该支路所对应的第二类幅值,该E为该支路经过伪码鉴相器后得到的自相关幅值。
可以理解的是,在该种实现方式中,通过将四路伪码序列中的每路伪码序列在同相支路和正交支路上的相干结果进行平方后再开方,得到各个伪码序列的自相关幅值的方法,避免了过度依赖载波环的工作形式和状态,即在载波环的锁相环没有达到稳定状态,或者载波环采用的是锁频环的形式时,则相位跟踪误差则无法确定,此时需要消除载波环相位跟踪误差对码环相关积分值的影响,则利用对码环鉴别器进行改良,将同相支路与正交支路上的相干积分结果进行平方后再相加再开方,由此获取超前、即时和滞后支路上的伪码序列的自相关幅值。
更进一步地,该预先构建的非相干超前混合滞后公式可为:
其中,所述δ为由伪码发生器产生的即时伪码序列与接收码之间的相位差,所述t为跟踪时间,b是程序调试获取的经验参数,所述E1为第二超前支路的自相关幅值、E为第一超前支路的自相关幅值,L1为第二滞后支路的自相关幅值、L为第一滞后支路的自相关幅值。
其中,在主流的具有六个积分清零器的跟踪环路的基础上,增加四个积分清零器,并利用非相干超前混合滞后法求解即时伪码序列和接收码的相位差,从而精确地将即时伪码序列的相位调整至与接收码的相位一致,从而使得能够有效地抑制短多径效用,取得精确的室内定位结果。
可以理解的是,在该种实现方式中,举例而言,中频数字信号在经过载波剥离运算器进行剥离载波后得到在同相支路上的第一剥离载波信号i(n)和正交支路上的第二载波信号q(n),且第一剥离载波信号i(n)和第二载波信号q(n)分别对应与本地产生的即时伪码序列在卷积运算器16c和卷积运算器16h中进行解扩,并且解扩后输出的结果分别为ip(n)与qp(n),并且具体表达式为:
ip(n)=AD(n)R(τp)cos(ωe(n)t(n)+θe)
qp(n)=AD(n)R(τp)sin(ωe(n)t(n)+θe)
其中,该A为输入中频信号的幅值,该D(n)是数据比特电平值,该R(·)代表最大值为1的即时伪码序列的自相关函数,τp为即时伪码序列与接收码之间的相位差,ωe和θe分别为角频率差和初相位误差。
基于上述解扩结果,再对该解扩结果在即时伪码序列对应的积分清零器15c和积分清零器15h中分别进行相干积分,则得到的第一剥离载波信号i(n)和第二剥离载波信号q(n)所对应的相干积分结果如下:
其中,fe为频率误差,Tcoh为相干积分时间,为载波环输出的载波相位差。
值得强调的是,超前与滞后支路的相干积分计算方法与上述即时支路的相干积分计算方法一致,差别仅在于在于τE与τL分别与τp相差d个码片,即在本申请的非相干超前混合滞后算法只需要计算超前与滞后支路的相干积分的情况下,只需要按照上述即时支路的相干积分计算方法计算即可。
当载波环采用锁相环形式并且锁相环工作在初始状态时,可按照码环鉴别器相干点积功率法计算即时伪码序列与接收码的相位误差,其中码环鉴别器相干点积功率计算公式为:
δ=(IE-IL)×IP
但是,通过研究各个支路的相干积分值计算原理,可知相干积分值的大小依赖于载波环相位跟踪误差因此,为消除此依赖关系提出利用对码环鉴别器进行改良,将同相支路与正交支路上的相干积分结果进行平方相加后再开方,且代入各个支路得到第一超前支路、第二超前支路、第一滞后支路、第二滞后支路和即时支路所对应的相干积分值所对应的自相关幅值分别为:
E=AR(τE)sinc(feTcoh)
E1=AR(τE1)sinc(feTcoh)
P=AR(τP)sinc(feTcoh)
L=AR(τL)sinc(feTcoh)
L1=AR(τL1)sinc(feTcoh)
值得强调的是,该将相干积分结果进行平方相加后再开方的方法为现有技术,在此不做详述。其中,该E为第一超前支路所对应的自相关幅值、该E1为第二超前支路所对应的自相关幅值、该L为第一滞后支路所对应的自相关幅值、该L1为第二滞后支路对应的自相关幅值;其中,本申请的非相干超前混合滞后计算公式并不需要计算即时支路所对应的自相关幅值P,在此只是对应举例说明。
其中,上述自相关幅值计算式的物理意义是同相支路与正交支路的能量之和,码环工作不依赖于载波环的载波相位,且主流的六个积分清零器对应的码环鉴别器的相位差的计算方法为:
从同相支路与正交支路的能量出发,得到码环鉴别器的非相干超前减滞后幅值法归一化计算公式:
且可以将超前支路与滞后支路的非相干积分功率直接相减,即得到非相干超前减滞后功率法归一化计算公式为:
值得强调的是,上述非相干超前减滞后幅值法归一化算法和非相干超前减滞后功率法归一化算法为现有技术,在此不做详述。
本发明结合主流的六个积分清零器对应的码环鉴别器的相位差的非相干超前减滞后功率法归一化计算方法对相位差进行计算,具体地:
采用非相干超前混合滞后公式并结合非相干超前减滞后功率法归一化计算公式对即时伪码序列和接收码的相位差进行计算,即给两个算法对应添加一个指数型参数权重,其中,给非相干超前减滞后功率法归一化计算公式添加的权重为(e-bt),给非相干超前混合滞后公式添加权重为(1-e-bt),即在环路开始工作阶段,主要通过非相干超前减滞后功率法归一化计算公式对该相位差进行计算,提高码环对码相位的跟踪速度;在环路工作稳定阶段,主要通过非相干超前混合滞后公式对该相位差进行计算,提高码环对码相位的跟踪精度。
其中,采用非相干超前混合滞后公式并结合非相干超前减滞后功率法归一化计算公式对即时伪码序列和接收码的相位差进行计算的具体公式为:
下面,结合图3对本发明实施例的码环***利用非相干超前混合滞后算法能够对短多径信号进行有效抑制进行验证。具体地,图3示出了直达相关函数、一条多径相关函数及由该直达相关函数和该多径相关函数合成的相关函数的示意图。
为方便说明,由于实际情况下多径信号对定位结果的影响是线性叠加的,所以对一条多径的研究就能很好地体现算法的效果,因此,下面只针对仅存在一条短多径的情况进行说明。
将直达信号能量与多径信号能量的比值定义为SMR=20lg|1/a|,且令a=Am/A,其中,该A为直达信号幅度,Am为多径信号幅度;并且由于实际中直达信号幅度大于多径信号幅度,所以a的绝对值小于1,且当a为正时代表直达信号与多径信号同相,a为负时代表直达信号与多径信号反向。
假设多径信号的延迟为τm,码片持续时间为常数t,且由于伪码序列的自相关函数主峰具有对称性,则可得到只存在直达信号的相关函数S和多径信号的相关函数Sm的表达式分别为:
并由于本发明主要是解决对短多径信号进行有效抑制问题,故将τm设置为小于0.25个码片,并在仿真时将τm具体设定为τm=0.1,且直达信号能量与多径信号能量的比值SMR设置为14dB。则具体可得如图3所示的直达信号的相关函数与多径信号的相关函数,以及由该直达相关函数和该多径相关函数合成的相关函数。
根据直达信号的相关函数与多径信号的相关函数在各个区间段的线性叠加关系,合成5段相关函数,表达式如下:
其中,1+t,-1≤t<-1+τ对应图3所示的①段;1+t+a(1+t-τ),-1+τ≤t<0对应图3所示的②段;1-t+a(1+t-τ),0≤t<τ对应图3所示的③段;1-t+a(1-t+τ),τ≤t<1对应图3所示的④段;a(1-t+τ),1≤t≤1+τ对应图3所示的⑤段。
设定积分清零器间距为d=0.25个码片,则可计算得到E1路,E路,L路和L1路对应的自相关幅值,并将t=-0.5,t=-0.25分别代入②段所对应的计算公式,将t=0.25,t=0.5分别代入④段所对应的计算公式,即可分别得到:E1=0.5+a(0.5-τ)、E=0.75+a(0.75-τ)、L=0.75+a(0.75+τ)、L1=0.5+a(0.5+τ)。
并将上述值分别对应代入非相干超前混合滞后计算公式,得到在短多径干扰下的相位误差:
公式化简得:
而在设定积分清零器间距为d=0.25,且由于此时计算的是受多径干扰下的误差,因此代入受干扰的E值,则得到利用非相干超前减滞后功率法在短多径干扰下多径误差估计计算公式为:
由上述公式可知当τm小于0.25个码片时,δm<δ,即可知单独利用传统的非相干超前减滞后功率法计算相位差时存在较大的误差,因此当有短多径干扰存在时,利用非相干超前混合滞后鉴相法计算可以有效的减少相位差。具体地,非相干超前混合滞后鉴相算法中的第二超前支路E1路和第一超前支路E路将同时增加a(1+t-τ),第二滞后支路L1路和第一滞后支路L路将同时增加a(1-t+τ),当超前支路混合滞后支路进行码环鉴相时则可有效抑制短多径干扰,使定位接收机在复杂室内环境中提供更高精度定位服务。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
本说明书中的各个实施例均采用相关的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于***实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
以上所述仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等,均包含在本发明的保护范围内。
Claims (5)
1.一种码环鉴别器,其特征在于,所述码环鉴别器包括:载波剥离运算器、伪码发生器、伪码环路滤波器、伪码鉴相器、十个积分清零器和十个卷积运算器;
所述十个卷积运算器中的每一卷积运算器的输入端均分别连接至所述载波剥离运算器的输出端和所述伪码发生器的输出端,所述十个卷积运算器输出端一一对应连接至所述十个积分清零器的输入端;
所述十个积分清零器的输出端连接至所述伪码鉴相器的输入端,所述伪码鉴相器的输出端连接至所述伪码环路滤波器的输入端,所述伪码环路滤波器输出端连接至所述伪码发生器;
其中,所述伪码鉴相器用于:求解所述伪码发生器产生的各个支路伪码序列经过相应的积分清零器后对应的自相关幅值;将得到的各个伪码序列对应的自相关幅值代入预设的非相干超前混合滞后公式,计算即时伪码序列和接收码的相位差;
其中,所述非相干超前混合滞后公式为:
<mrow>
<mi>&delta;</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>b</mi>
<mi>t</mi>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<msup>
<mi>E</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msup>
<mi>E</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>b</mi>
<mi>t</mi>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mi>E</mi>
<mn>1</mn>
<mo>&times;</mo>
<mi>L</mi>
<mo>-</mo>
<mi>E</mi>
<mo>&times;</mo>
<mi>L</mi>
<mn>1</mn>
</mrow>
<mrow>
<mi>E</mi>
<mn>1</mn>
<mo>&times;</mo>
<mi>L</mi>
<mo>+</mo>
<mi>E</mi>
<mo>&times;</mo>
<mi>L</mi>
<mn>1</mn>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
其中,所述δ为由伪码发生器产生的即时伪码序列与接收码之间的相位差,所述t为跟踪时间,b是程序调试获取的经验参数,所述E1为所述伪码发生器产生的第二超前支路伪码序列的自相关幅值、E为所述伪码发生器产生的第一超前支路伪码序列对应的自相关幅值,L1为所述伪码发生器产生的第二滞后支路伪码序列的自相关幅值、L为所述伪码发生器产生的第一滞后支路伪码序列的自相关幅值。
2.根据权利要求1所述的码环鉴别器,其特征在于,所述十个卷积运算器中的每一卷积运算器的输入端均分别连接至所述载波剥离运算器的输出端和所述伪码发生器的输出端,包括:
所述十个卷积运算器中的五个卷积运算器的输入端均分别连接至所述载波剥离运算器的同相支路输出端,且所述伪码发生器的第一超前支路、第二超前支路、即时支路、第一滞后支路和第二滞后支路的输出端分别与所述五个卷积运算器的输入端唯一对应连接;
所述十个卷积运算器中的剩余五个卷积运算器的输入端均分别连接至所述载波剥离运算器的正交支路输出端,且所述伪码发生器的所述第一超前支路、所述第二超前支路、所述即时支路、所述第一滞后支路和所述第二滞后支路的输出端分别与所述剩余五个卷积运算器的输入端唯一对应连接。
3.根据权利要求2所述的码环鉴别器,其特征在于,
所述伪码发生器的所述第一超前支路超前即时支路d个码片,所述第二超前支路超前所述第一超前支路d个码片;
所述伪码发生器的所述第一滞后支路滞后即时支路d个码片,所述第二滞后支路滞后所述第一滞后支路d个码片,其中,0.2≤d≤0.3。
4.一种短多径抑制方法,其特征在于,应用于权利要求1-3中任一项所述的码环鉴别器,所述方法包括:
将所述码环鉴别器中的伪码发生器产生的与第一超前支路对应的第一超前支路伪码序列、与第二超前支路对应的第二超前支路伪码序列、与第一滞后支路对应的第一滞后支路伪码序列和与第二滞后支路对应的第二滞后支路伪码序列,分别与在载波剥离运算器中进行载波剥离后的中频数字信号在八路并行的卷积运算器中进行解扩,以得到解扩后的八路信号,其中,所述中频数字信号在经载波剥离运算器进行载波剥离后得到在同相支路进行载波剥离后的第一剥离信号和在正交支路进行载波剥离后的第二剥离信号;
在所述八个积分清零器中对所述解扩后的八路信号进行相关运算,以得到解扩时与同相支路所相关的第一超前支路、第二超前支路、第一滞后支路和第二滞后支路所对应的第一类幅值,和解扩时与正交支路所相关的第一超前支路、第二超前支路、第一滞后支路和第二滞后支路所对应的第二类幅值;
基于所得到的第一类幅值和第二类幅值,计算所述伪码发生器中的各个支路经过积分清零器后的自相关幅值;
基于所述伪码发生器中的各个支路的自相关幅值和预设的非相干超前混合滞后公式,计算由伪码发生器产生的即时支路伪码序列与接收码的相位差,其中,所述第一剥离信号和所述第二剥离信号中携带有接收码;
根据所述相位差,调节所述伪码发生器产生的即时伪码序列的相位,以使即时支路伪码序列与接收码在相位上保持一致;
所述预设的非相干超前混合滞后公式为:
<mrow>
<mi>&delta;</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>b</mi>
<mi>t</mi>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<msup>
<mi>E</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msup>
<mi>E</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>L</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>b</mi>
<mi>t</mi>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mi>E</mi>
<mn>1</mn>
<mo>&times;</mo>
<mi>L</mi>
<mo>-</mo>
<mi>E</mi>
<mo>&times;</mo>
<mi>L</mi>
<mn>1</mn>
</mrow>
<mrow>
<mi>E</mi>
<mn>1</mn>
<mo>&times;</mo>
<mi>L</mi>
<mo>+</mo>
<mi>E</mi>
<mo>&times;</mo>
<mi>L</mi>
<mn>1</mn>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
其中,所述δ为由伪码发生器产生的即时伪码序列与接收码之间的相位差,所述t为跟踪时间,b是程序调试获取的经验参数,所述E1为第二超前支路的自相关幅值、E为第一超前支路的自相关幅值,L1为第二滞后支路的自相关幅值、L为第一滞后支路的自相关幅值。
5.根据权利要求4所述的方法,其特征在于,计算伪码发生器中的任一支路的自相关幅值的预设自相关幅值计算公式具体为:
<mrow>
<mi>E</mi>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<msub>
<mi>I</mi>
<mi>E</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<msub>
<mi>Q</mi>
<mi>E</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>,</mo>
</mrow>
其中,所述IE是支路所对应的第一类幅值,所述QE是所述支路所对应的第二类幅值,所述E为所述支路经过积分清零器计算后得到的自相关幅值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510996958.5A CN105510932B (zh) | 2015-12-25 | 2015-12-25 | 一种码环鉴别器及短多径抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510996958.5A CN105510932B (zh) | 2015-12-25 | 2015-12-25 | 一种码环鉴别器及短多径抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105510932A CN105510932A (zh) | 2016-04-20 |
CN105510932B true CN105510932B (zh) | 2018-05-11 |
Family
ID=55719041
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510996958.5A Active CN105510932B (zh) | 2015-12-25 | 2015-12-25 | 一种码环鉴别器及短多径抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105510932B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107255822A (zh) * | 2017-06-01 | 2017-10-17 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 多径环境下gnss接收机信号参数估计方法 |
CN109743075B (zh) * | 2019-03-14 | 2020-12-04 | 哈尔滨工程大学 | 一种三环联动的非同源扩频码跟踪环路同步装置及方法 |
CN114793190B (zh) * | 2022-03-09 | 2023-10-31 | 北京邮电大学 | 一种tc-ofdm定位***的多径效应抑制方法和装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101776752A (zh) * | 2010-01-29 | 2010-07-14 | 中国科学院空间科学与应用研究中心 | 机群链路高动态信号的精密跟踪与测量方法 |
CN102013972A (zh) * | 2010-11-10 | 2011-04-13 | 北京航天自动控制研究所 | 一种载波误锁修正方法 |
CN102183770A (zh) * | 2011-05-20 | 2011-09-14 | 哈尔滨工程大学 | 一种抗多径干扰的gps伪随机码跟踪环路及其抗多径干扰方法 |
CN102375149A (zh) * | 2010-08-18 | 2012-03-14 | 中国科学院微电子研究所 | 多相关支路混合间距延迟锁定环及其应用方法 |
CN102801442A (zh) * | 2012-07-25 | 2012-11-28 | 北京邮电大学 | 一种测距的方法、装置和接收机 |
-
2015
- 2015-12-25 CN CN201510996958.5A patent/CN105510932B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101776752A (zh) * | 2010-01-29 | 2010-07-14 | 中国科学院空间科学与应用研究中心 | 机群链路高动态信号的精密跟踪与测量方法 |
CN102375149A (zh) * | 2010-08-18 | 2012-03-14 | 中国科学院微电子研究所 | 多相关支路混合间距延迟锁定环及其应用方法 |
CN102013972A (zh) * | 2010-11-10 | 2011-04-13 | 北京航天自动控制研究所 | 一种载波误锁修正方法 |
CN102183770A (zh) * | 2011-05-20 | 2011-09-14 | 哈尔滨工程大学 | 一种抗多径干扰的gps伪随机码跟踪环路及其抗多径干扰方法 |
CN102801442A (zh) * | 2012-07-25 | 2012-11-28 | 北京邮电大学 | 一种测距的方法、装置和接收机 |
Non-Patent Citations (2)
Title |
---|
GNSS软件接收机码环鉴别器性能评估与改进研究;孟骞 等;《第五届中国卫星导航学术年会论文集-S7北斗/GNSS用户终端技术》;20140521;第1-7页 * |
一种改进的BOC信号码跟踪环路设计方法;杨文津 等;《无线电工程》;20140215(第2期);第21-24页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105510932A (zh) | 2016-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7822105B2 (en) | Cross-correlation removal of carrier wave jamming signals | |
US7116704B2 (en) | Strong signal cancellation to enhance processing of weak spread spectrum signal | |
JP3964621B2 (ja) | 不均一で繋がりのない複数のサンプル・セグメントの相関分析を用いる信号検出器 | |
US8170082B2 (en) | Cross-talk mitigation in global navigation satellite systems | |
TW444475B (en) | Mobile station assisted timing synchronization in a CDMA communication system | |
EP1916535B1 (en) | GNSS receiver with cross-correlation rejection | |
CN102338878B (zh) | 一种多径抑制boc码跟踪方法以及码跟踪环 | |
CN109407121B (zh) | 一种可配置的导航信号兼容捕获跟踪装置 | |
JP2002540434A5 (ja) | 不均一で繋がりのない複数のサンプル・セグメントの相関分析を用いる信号検出器 | |
CN105510932B (zh) | 一种码环鉴别器及短多径抑制方法 | |
RU2658625C1 (ru) | Способ формирования сигнала c расширенным спектром, устройство формирования сигнала, способ приема сигнала и приемное устройство | |
JP2011524979A (ja) | 符号空間探索における多重相関処理 | |
Shamaei et al. | Pseudorange and multipath analysis of positioning with LTE secondary synchronization signals | |
JP2014228536A (ja) | 追跡ループのステータスを決定するための装置および方法 | |
CN105372678B (zh) | 一种正弦boc调制信号的无模糊跟踪方法 | |
US20110032974A1 (en) | Carrier frequency synchronization detection circuit and correlation calculator | |
JP5765902B2 (ja) | 衛星信号のコード追尾装置 | |
AU2002305881A1 (en) | System and method of estimation of time of arrival | |
JP2895399B2 (ja) | 同期追従方法 | |
JP2004219076A (ja) | Gps受信機 | |
JP5933559B2 (ja) | 信号捕捉装置及び方法 | |
Schreiber et al. | Software defined radio based receiver for TDOA positioning system | |
KR20180061625A (ko) | AltBOC(15,10) 신호를 추적하는 방법 및 AltBOC(15,10) 신호 추적 장치 | |
US9236903B2 (en) | Multi-path detection | |
JP2008252162A (ja) | 拡散信号受信装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20211229 Address after: No. 2056, 16th floor, scientific research building, No. 13 Jiuxianqiao Road, Chaoyang District, Beijing 100015 Patentee after: Beijing Duwei Technology Co.,Ltd. Address before: 100876 Beijing city Haidian District Xitucheng Road No. 10 Patentee before: Beijing University of Posts and Telecommunications |