CN109239744B - 一种基于复信号相位的快速抗比特反转快速捕获方法 - Google Patents

一种基于复信号相位的快速抗比特反转快速捕获方法 Download PDF

Info

Publication number
CN109239744B
CN109239744B CN201811221034.8A CN201811221034A CN109239744B CN 109239744 B CN109239744 B CN 109239744B CN 201811221034 A CN201811221034 A CN 201811221034A CN 109239744 B CN109239744 B CN 109239744B
Authority
CN
China
Prior art keywords
code phase
code
representing
bit
phase
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
Application number
CN201811221034.8A
Other languages
English (en)
Other versions
CN109239744A (zh
Inventor
吴超
孙闽红
刘二小
滕旭阳
简志华
钟华
苏明坤
乔磊
仇兆炀
宋慧娜
汪立新
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Dianzi University
Original Assignee
Hangzhou Dianzi University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hangzhou Dianzi University filed Critical Hangzhou Dianzi University
Priority to CN201811221034.8A priority Critical patent/CN109239744B/zh
Publication of CN109239744A publication Critical patent/CN109239744A/zh
Application granted granted Critical
Publication of CN109239744B publication Critical patent/CN109239744B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/24—Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/30—Acquisition or tracking or demodulation of signals transmitted by the system code related

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明主要公开了一种基于复信号相位的快速抗比特反转快速捕获方法，通过对接收信号进行下变频处理，并对本地相关信号和接收信号进行分块，将得到矩的阵进行码相位搜索处理，得到待检测变量，再通过相干积分检测，得到码相位搜索处理的相干积分值，对接收信号进行码相位处理得到码相位处理待检测变量，得出的估算的比特数据反转位置，再进行相干积分检测，最后得到待求参数码相位，最终获得估计的捕获参数码相位和频率，从而实现了GNSS信号的快速捕获。本发明方法可以快速精确估计码相位，从而实现比特符号反转情况下捕获参数估计。

Description

一种基于复信号相位的快速抗比特反转快速捕获方法
技术领域
本发明属于通信技术领域，具体涉及一种基于复信号相位的快速抗比特反转快速捕获方法。
背景技术
随着科技水平的进步，无线通信技术和全球卫星定位***(GPS)技术越来越多地应用于日常生活的方方面面开始。其中，全球导航卫星***(Global NavigationSatellite System，缩写为GNSS)，其定位是利用一组卫星的伪距、星历、卫星发射时间等观测量来的，同时还必须知道用户钟差。卫星导航接收机捕获并跟踪多个GNSS卫星的信号，然后解调其中调制的导航数据。卫星导航接收机利用测距码计算GNSS卫星与用户的相对距离，利用导航数据中的星历数据解算卫星位置和时间模型，进而计算出用户的位置。卫星导航定位技术目前已基本取代了地基无线电导航、传统大地测量和天文测量导航定位技术，并推动了大地测量与导航定位领域的全新发展。因此，GNSS***是国家安全和社会经济发展的基础设施。
在检测接收GNSS信号(即卫星信号捕获)过程中，待估计参量频率和码相位会受到数据反转的影响(比如GPS L1 C/A码接收信号，其中L1为载波信号，C/A码为GPS卫星发出的一种伪随机码用于粗测距和捕获GPS卫星)，导致检测峰值降低，从而使得估计参数不准确影响检测速度。
为了提升检测效率，当代主要捕获方法采用基于快速傅里叶变换方法估计捕获参数。但是捕获速度有待进一步提升，而且待检测信号数据比特反转的随机性，需要制定更有效的捕获策略。为了进一步提升捕获速度，降低捕获算法复杂度，需要进一步研究受到比特数据符号反转情况下的捕获参数快速捕获问题。
发明内容
本发明的目的是为了解决在数据比特反转情况下提升捕获方法速度的问题，提供一种快速且精准的估计GNSS信号的码相位和频率，从而实现基于复信号相位的快速抗比特反转快速捕获方法。
为了达到上述目的，本发明所采取的具体技术方案为：
一种基于复信号相位的快速抗比特反转快速捕获方法，包括以下步骤：
S1.对接收信号R(n)进行下变频处理，并对本地相关信号和接收信号进行分块，得到分块后的矩阵
Figure GDA0002620616840000021
和
Figure GDA0002620616840000022
其中，i表示第i次积分时间，D表示有数据比特影响时本地伪随机码、ND表示没有数据比特影响时本地伪随机码，M、N表示分块矩阵参量；
S2.将S1得到矩阵进行码相位搜索处理，得到待检测变量
Figure GDA0002620616840000023
和
Figure GDA0002620616840000024
再通过相干积分检测，得到码相位搜索处理的相干积分值
Figure GDA0002620616840000025
其中i0为估算的比特数据反转位置，fd,k0为估计频率，dr为码相位搜索处理估算码相位参数值，i表示第i次积分时间，D表示有数据比特影响时本地伪随机码、ND表示没有数据比特影响时本地伪随机码，fd,k表示搜索频率，dc表示搜索的码片偏移；
S3.对接收信号进行码相位处理得到码相位处理待检测变量
Figure GDA0002620616840000027
根据S1得出的估算的比特数据反转位置i0，对
Figure GDA0002620616840000028
进行相干积分检测，最后得到待求参数码相位，最终获得估计的捕获参数码相位和频率，从而实现了GNSS信号的快速捕获。
进一步的，接收信号R(n)为：
Figure GDA0002620616840000026
其中，q为接收信号的码相位，fd为多普勒频率，C()为伪随机码，D()为调制数据，fs为采样频率，p(n)表示实部和虚部均服从均值为0、方差为σ2的高斯噪声，n表示采样点，n＝0,1,…,
接收信号R(n)进行下变频处理，得到下变频信号为：
Figure GDA0002620616840000031
其中，下变频信号f(n)为：
Figure GDA0002620616840000032
其中，fd,k为搜索频率，fs为采样频率，n表示采样点。
两种本地相关信号分别为：
LD(n)＝(-1)k-1C(n)；
LND(n)＝C(n)；
其中，n＝iT+k0，T为PN码周期采样点数，k0＝0,...,T-1，C()为伪随机码,i表示第i次积分时间，LD(n)和LND(n)分别表示有数据比特影响时本地伪随机码、没有数据比特影响时本地伪随机码；
将本地相关信号和下变频信号进行分块(以下公式涉及的所有字母符号等需解释，前面若有出现，此处也需重新解释)：
Figure GDA0002620616840000033
Figure GDA0002620616840000041
Figure GDA0002620616840000042
其中，NL＝2(i-1)MN+dc，i＝1,...,Nc，Nc为相干积分时间，NR＝(i-1)MN，M、N表示分块矩阵参量，i表示第i次积分时间。
进一步的，所述将S1得到矩阵进行码相位搜索处理，具体方法如下：
Figure GDA0002620616840000043
Figure GDA0002620616840000044
其中，mod((c2+n),2N2)表示(c2+n)除以2N2的余数，i表示第i次积分时间，D表示有数据比特影响时本地伪随机码、ND表示没有数据比特影响时本地伪随机码，fd,k表示搜索频率，dc表示搜索的码片偏移，当dc-1<0，dc-1＝dc-1+M；如果dc+1＞M，dc+1＝dc+1-M，通过相干积分检测：
Figure GDA0002620616840000045
Figure GDA0002620616840000046
Figure GDA0002620616840000051
其中，i1表示数据比特符号可能的翻转位置，D表示有数据比特影响时本地伪随机码、ND表示没有数据比特影响时本地伪随机码，fd,k表示搜索频率，dc表示搜索的码片偏移，dr表示码相位搜索处理预测的码偏移量，i0表示预测数据比特翻转位置，fd,k0为估计的信号频率。
进一步的，所述对接收信号进行码相位处理的码相位估计处理，具体方法如下：
Figure GDA0002620616840000052
Figure GDA0002620616840000053
其中，D表示有数据比特影响时本地伪随机码、ND表示没有数据比特影响时本地伪随机码，
Figure GDA0002620616840000054
表示在码相位处理的码相位估计时没有比特翻转时第iN次积分时间码片偏移为d2时接收信号积分值，
Figure GDA0002620616840000055
表示在码相位处理的码相位估计时有比特翻转时第iN次积分时间码片偏移为d2时接收信号积分值，iN表示没有比特翻转时第iN次积分时间，iD＝i0+1表示有数据比特翻转时第iD积分时间，M、N表示分块矩阵参量，i0表示预测数据比特翻转位置，N＝M2N2，M2、N2表示码相位预测参量，
再进行相干积分检测：
Figure GDA0002620616840000056
Figure GDA0002620616840000057
Figure GDA0002620616840000061
其中，dr2表示码相位处理的码相位预测参量，
Figure GDA0002620616840000062
码相位处理的码相位预测时的相干积分后的检测量，D表示有数据比特影响时本地伪随机码、ND表示没有数据比特影响时本地伪随机码，i0表示预测数据比特翻转位置，d2表示码相位处理的码相位预测时码片偏移，根据PN码相关性质，可得当d2＝dr2时取得最大值为：
Figure GDA0002620616840000063
其中，cr2、dr2表示码相位处理的码相位预测参量，M表示分块矩阵参量，M2、N2表示码相位预测参量，M2、N2表示码相位处理的码相位预测分块矩阵参量。令X2(dr2)相位角为
Figure GDA0002620616840000064
则
Figure GDA0002620616840000065
其中，cr2表示码相位处理的码相位预测参量，N2表示码相位预测参量。最后，接收的估计码相位可以表示为
q＝(cr2M2+dr2)M+dr；
其中，cr2、dr2表示码相位处理的码相位预测参量，dr表示码相位搜索处理的码相位预测参量，M表示码相位搜索处理的码相位预测分块矩阵参量，M2表示码相位处理的码相位预测分块矩阵参量。
采用以上技术方案，本发明所具有的有益效果是：本发明的比特符号反转可以在积分区间内部，采用两步基于变换信号复相位估算信号码相位快速估计方法，将码相位分为不同区间段估计码相位。基于上述两步估计码相位方法可以快速精确估计码相位，从而实现比特符号反转情况下捕获参数估计。
附图说明
图1为本发明一种基于复信号相位的快速抗比特反转快速捕获方法的流程图。
图2为现有技术方法与本发明基于复信号相位的快速抗比特反转快速捕获方法对比流程图。
图3为现有技术方法与本发明基于复信号相位的快速抗比特反转快速捕获方法的复数加法计算量对比图。
图4为现有技术方法与本发明基于复信号相位的快速抗比特反转快速捕获方法的复数乘法计算量对比图。
具体实施方式
如图1所示，一种基于复信号相位的快速抗比特反转快速捕获方法，其特征在于，包括以下步骤：
S1.对接收信号R(n)进行下变频处理，并对本地相关信号和接收信号进行分块，得到分块后的矩阵
Figure GDA0002620616840000071
和
Figure GDA0002620616840000072
其中，i表示第i次积分时间，D表示有数据比特影响时本地伪随机码、ND表示没有数据比特影响时本地伪随机码，M、N表示分块矩阵参量。
接收信号R(n)为：
Figure GDA0002620616840000073
其中，q为接收信号的码相位，fd为多普勒频率，C()为伪随机码，D()为调制数据，fs为采样频率，p(n)表示实部和虚部均服从均值为0、方差为σ2的高斯噪声，n表示采样点，n＝0,1,…,
接收信号R(n)进行下变频处理，得到下变频信号为：
Figure GDA0002620616840000074
其中，下变频信号f(n)为：
Figure GDA0002620616840000075
其中，fd,k为搜索频率，fs为采样频率，n表示采样点。
两种本地相关信号分别为：
LD(n)＝(-1)k-1C(n)；
LND(n)＝C(n)；
其中，n＝iT+k0，T为PN码周期采样点数，k0＝0,...,T-1，C()为伪随机码,i表示第i次积分时间，LD(n)和LND(n)分别表示有数据比特影响时本地伪随机码、没有数据比特影响时本地伪随机码；
将本地相关信号和下变频信号进行分块：
Figure GDA0002620616840000081
Figure GDA0002620616840000082
Figure GDA0002620616840000083
其中，NL＝2(i-1)MN+dc，i＝1,...,Nc，Nc为相干积分时间，NR＝(i-1)MN，M、N表示分块矩阵参量，i表示第i次积分时间。
S2.将S1得到矩阵进行码相位搜索处理，得到待检测变量
Figure GDA0002620616840000091
和
Figure GDA0002620616840000092
再通过相干积分检测，得到码相位搜索处理的相干积分值
Figure GDA0002620616840000093
其中i0为估算的比特数据反转位置，
Figure GDA0002620616840000094
为估计频率，dr为码相位搜索处理估算码相位参数值，i表示第i次积分时间，D表示有数据比特影响时本地伪随机码、ND表示没有数据比特影响时本地伪随机码，fd,k表示搜索频率，dc表示搜索的码片偏移。
所述将S1得到矩阵进行码相位搜索处理，具体方法如下：
Figure GDA0002620616840000095
Figure GDA0002620616840000096
其中，M、N表示分块矩阵参量，i表示第i次积分时间，D表示有数据比特影响时本地伪随机码、ND表示没有数据比特影响时本地伪随机码，fd,k表示搜索频率，dc表示搜索的码片偏移，当dc-1<0，dc-1＝dc-1+M；如果dc+1＞M，dc+1＝dc+1-M，通过相干积分检测：
Figure GDA0002620616840000097
Figure GDA0002620616840000098
Figure GDA0002620616840000099
其中，i1表示数据比特符号可能的翻转位置，D表示有数据比特影响时本地伪随机码、ND表示没有数据比特影响时本地伪随机码，fd,k表示搜索频率，dc表示搜索的码片偏移，dr表示码相位搜索处理预测的码偏移量，i0表示预测数据比特翻转位置，fd,k0为估计的信号频率。
S3.对接收信号进行码相位处理得到码相位处理待检测变量
Figure GDA0002620616840000101
根据S1得出的估算的比特数据反转位置i0，对
Figure GDA0002620616840000102
进行相干积分检测，最后得到待求参数码相位，最终获得估计的捕获参数码相位和频率，从而实现了GNSS信号的快速捕获。
所述对接收信号进行码相位处理的码相位估计处理，具体方法如下：
Figure GDA0002620616840000103
Figure GDA0002620616840000104
其中，D表示有数据比特影响时本地伪随机码、ND表示没有数据比特影响时本地伪随机码，
Figure GDA0002620616840000105
表示在码相位处理的码相位估计时没有比特翻转时第iN次积分时间码片偏移为d2时接收信号积分值，
Figure GDA0002620616840000106
表示在码相位处理的码相位估计时有比特翻转时第iN次积分时间码片偏移为d2时接收信号积分值，iN表示没有比特翻转时第iN次积分时间，iD＝i0+1表示有数据比特翻转时第iD积分时间，M、N表示分块矩阵参量，i0表示预测数据比特翻转位置，N＝M2N2，M2、N2表示码相位预测参量，
再进行相干积分检测：
Figure GDA0002620616840000107
Figure GDA0002620616840000108
Figure GDA0002620616840000111
其中，dr2表示码相位处理的码相位预测参量，
Figure GDA0002620616840000112
码相位处理的码相位预测时的相干积分后的检测量，D表示有数据比特影响时本地伪随机码、ND表示没有数据比特影响时本地伪随机码，i0表示预测数据比特翻转位置，d2表示码相位处理的码相位预测时码片偏移，根据PN码相关性质，可得当d2＝dr2时取得最大值为：
Figure GDA0002620616840000113
其中，cr2、dr2表示码相位处理的码相位预测参量，M表示分块矩阵参量，M2、N2表示码相位预测参量，M2、N2表示码相位处理的码相位预测分块矩阵参量。令X2(dr2)相位角为
Figure GDA0002620616840000114
则
Figure GDA0002620616840000115
其中，cr2表示码相位处理的码相位预测参量，N2表示码相位预测参量。最后，接收的估计码相位可以表示为
q＝(cr2M2+dr2)M+dr；
其中，cr2、dr2表示码相位处理的码相位预测参量，dr表示码相位搜索处理的码相位预测参量，M表示码相位搜索处理的码相位预测分块矩阵参量，M2表示码相位处理的码相位预测分块矩阵参量。
如图2所示，其中，C(n)为本地信号，NF为单位积分时间内采样点数。R0(m)为没有数据调制一个伪码周期得到积分值，R1(m)为有数据调制一个伪码周期得到积分值，m为码片偏移。假设第i个伪码周期没有数据调制时积分值为R0，i(m)，假设第i个伪码周期没有数据调制时积分值为R1，i(m)，则上图积累结果为：
Figure GDA0002620616840000121
Figure GDA0002620616840000122
其中K为积累时间，k＝1，…，K。根据上述积累峰值可以得出最大值，最大值对应码片偏移即为接收信号码相位估计值。
表1是根据现有技术和提出方法流出，得出现有技术方法与本发明基于复信号相位的快速抗比特反转快速捕获方法的复加计算量和复乘计算量对比：
Figure GDA0002620616840000123
表1
其中Nc表示积分时间，NF单位积分时间内采样点数。其中M、N、M2和N2如公式步骤S2及S3表示的参数，表示第一、二步分块参数。
图3是现有技术方法与本发明基于复信号相位的快速抗比特反转快速捕获方法的复数加法计算量对比图。
图4是现有技术方法与本发明基于复信号相位的快速抗比特反转快速捕获方法的复数乘法计算量对比图。
从图3、4，我们可以发现本发明方法有采用将估计码相位转换到估计复信号相位上比现有基于快速傅里叶变换消耗更加少的复乘计算量和复加计算量，所以可以更加快速估计捕获参数，达到快速捕获的目的。
其中，Nc＝1～11ms，NF＝1023，M＝N＝32，M2＝4，N2＝8，P1＝P2＝8。从上图可以看出提出方法比当代方法可以极大地减小计算量。
其中，FFT是FFT(Fast Fourier Transformation)是离散傅氏变换(DFT)的快速算法。
本文中伪随机码是指(Pseudo-Noise，PN码)具有与白噪声类似的自相关性质的0和1所构成的编码序列。
本领域的技术人员应该理解，本发明的具体实施方式仅用于解释本发明的原理，而并不限制本发明。凡是依据本发明中的设计精神所做出的等效变化或修饰，均应落入本发明的保护范围。

Claims (4)

1.一种基于复信号相位的快速抗比特反转快速捕获方法，其特征在于，包括以下步骤：
S1.对接收信号R(n)进行下变频处理，并对本地相关信号和接收信号进行分块，得到分块后的矩阵
Figure FDA0002620616830000011
和
Figure FDA0002620616830000012
其中，i表示第i次积分时间，D表示有数据比特影响时本地伪随机码、ND表示没有数据比特影响时本地伪随机码，M、N表示分块矩阵参量；
S2.将S1得到矩阵进行码相位搜索处理，得到待检测变量
Figure FDA0002620616830000013
和
Figure FDA0002620616830000014
再通过相干积分检测，得到码相位搜索处理的相干积分值
Figure FDA0002620616830000015
其中i0为估算的比特数据反转位置，
Figure FDA0002620616830000016
为估计频率，dr为码相位搜索处理估算码相位参数值，i表示第i次积分时间，D表示有数据比特影响时本地伪随机码、ND表示没有数据比特影响时本地伪随机码，fd,k表示搜索频率，dc表示搜索的码片偏移；
S3.对接收信号进行码相位处理得到码相位处理待检测变量
Figure FDA0002620616830000017
根据S2得出的估算的比特数据反转位置i0，对
Figure FDA0002620616830000018
进行相干积分检测，最后得到待求参数码相位，最终获得估计的捕获参数码相位和频率，从而实现了GNSS信号的快速捕获。
2.如权利要求1所述的一种基于复信号相位的快速抗比特反转快速捕获方法，其特征在于，接收信号R(n)为：
Figure FDA0002620616830000019
其中，q为接收信号的码相位，fd为多普勒频率，C()为伪随机码，D()为调制数据，fs为采样频率，p(n)表示实部和虚部均服从均值为0、方差为σ2的高斯噪声，n表示采样点，n＝0,1,…,
接收信号R(n)进行下变频处理，得到下变频信号为：
Figure FDA0002620616830000021
其中，下变频信号f(n)为：
Figure FDA0002620616830000022
其中，fd,k为搜索频率，fs为采样频率，n表示采样点，
两种本地相关信号分别为：
LD(n)＝(-1)iC(n)；
LND(n)＝C(n)；
其中，n表示采样点数，n可以写成iT+k0，T为PN码周期采样点数，k0＝0,...,T-1，C()为伪随机码,i表示第i次积分时间，LD(n)和LND(n)分别表示有数据比特影响时本地伪随机码、没有数据比特影响时本地伪随机码；将本地相关信号和下变频信号进行分块：
Figure FDA0002620616830000023
Figure FDA0002620616830000024
Figure FDA0002620616830000031
其中，NL＝2(i-1)MN+dc，i＝1,...,Nc，Nc为相干积分时间，NR＝(i-1)MN，M、N表示分块矩阵参量，i表示第i次积分时间。
3.如权利要求1所述的一种基于复信号相位的快速抗比特反转快速捕获方法，其特征在于：
所述将S1得到矩阵进行码相位搜索处理，具体方法如下：
Figure FDA0002620616830000032
其中，M、N表示分块矩阵参量，i表示第i次积分时间，D表示有数据比特影响时本地伪随机码、ND表示没有数据比特影响时本地伪随机码，fd,k表示搜索频率，dc表示搜索的码片偏移，当dc-1<0，dc-1＝dc-1+M；如果dc+1＞M，dc+1＝dc+1-M，通过相干积分检测：
Figure FDA0002620616830000033
Figure FDA0002620616830000034
Figure FDA0002620616830000035
其中，i1表示数据比特符号可能的翻转位置，D表示有数据比特影响时本地伪随机码、ND表示没有数据比特影响时本地伪随机码，fd,k表示搜索频率，dc表示搜索的码片偏移，dr表示码相位搜索处理预测的码偏移量，i0表示预测数据比特翻转位置，fd,k0为估计的信号频率。
4.如权利要求1所述的一种基于复信号相位的快速抗比特反转快速捕获方法，其特征在于：
所述对接收信号进行码相位处理的码相位估计处理，具体方法如下：
Figure FDA0002620616830000041
Figure FDA0002620616830000042
其中，mod((c2+n),2N2)表示(c2+n)除以2N2的余数；D表示有数据比特影响时本地伪随机码、ND表示没有数据比特影响时本地伪随机码，
Figure FDA0002620616830000043
表示在码相位处理的码相位估计时没有比特翻转时第iN次积分时间码片偏移为d2时接收信号积分值，
Figure FDA0002620616830000044
表示在码相位处理的码相位估计时有比特翻转时第iN次积分时间码片偏移为d2时接收信号积分值，iN表示没有比特翻转时第iN次积分时间，iD＝i0+1表示有数据比特翻转时第iD积分时间，M、N表示分块矩阵参量，i0表示预测数据比特翻转位置，N＝M2N2，M2、N2表示码相位预测参量，
再进行相干积分检测：
Figure FDA0002620616830000045
Figure FDA0002620616830000051
Figure FDA0002620616830000052
其中，dr2表示码相位处理的码相位预测参量，
Figure FDA0002620616830000053
码相位处理的码相位预测时的相干积分后的检测量，D表示有数据比特影响时本地伪随机码、ND表示没有数据比特影响时本地伪随机码，i0表示预测数据比特翻转位置，d2表示码相位处理的码相位预测时码片偏移，根据PN码相关性质，可得当d2＝dr2时取得最大值为：
Figure FDA0002620616830000054
其中，cr2、dr2表示码相位处理的码相位预测参量，M表示分块矩阵参量，M2、N2表示码相位预测参量,令X2(dr2)相位角为
Figure FDA0002620616830000056
则
Figure FDA0002620616830000055
其中，cr2表示码相位处理的码相位预测参量，N2表示码相位预测参量。最后，接收的估计码相位可以表示为
q＝(cr2M2+dr2)M+dr；
其中，cr2、dr2表示码相位处理的码相位预测参量，dr表示码相位搜索处理的码相位预测参量，M表示码相位搜索处理的码相位预测分块矩阵参量，M2表示码相位处理的码相位预测分块矩阵参量。
CN201811221034.8A 2018-10-19 2018-10-19 一种基于复信号相位的快速抗比特反转快速捕获方法 Active CN109239744B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811221034.8A CN109239744B (zh) 2018-10-19 2018-10-19 一种基于复信号相位的快速抗比特反转快速捕获方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811221034.8A CN109239744B (zh) 2018-10-19 2018-10-19 一种基于复信号相位的快速抗比特反转快速捕获方法

Publications (2)

Publication Number Publication Date
CN109239744A CN109239744A (zh) 2019-01-18
CN109239744B true CN109239744B (zh) 2020-10-09

Family

ID=65080672

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811221034.8A Active CN109239744B (zh) 2018-10-19 2018-10-19 一种基于复信号相位的快速抗比特反转快速捕获方法

Country Status (1)

Country Link
CN (1) CN109239744B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008107160A (ja) * 2006-10-24 2008-05-08 Toyota Motor Corp 移動体用衛星電波受信機
JP2011203234A (ja) * 2010-10-19 2011-10-13 Seiko Epson Corp 信号捕捉方法、信号捕捉装置及び電子機器
CN102608626A (zh) * 2012-03-09 2012-07-25 暨南大学 一种高灵敏度卫星导航信号捕获方法及***
CN106093981A (zh) * 2016-06-06 2016-11-09 电子科技大学 基于优化并行码相位搜索的gps捕获电路
CN107247276A (zh) * 2017-07-05 2017-10-13 西安电子科技大学 基于重叠多块补零算法的弱信号捕获方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008107160A (ja) * 2006-10-24 2008-05-08 Toyota Motor Corp 移動体用衛星電波受信機
JP2011203234A (ja) * 2010-10-19 2011-10-13 Seiko Epson Corp 信号捕捉方法、信号捕捉装置及び電子機器
CN102608626A (zh) * 2012-03-09 2012-07-25 暨南大学 一种高灵敏度卫星导航信号捕获方法及***
CN106093981A (zh) * 2016-06-06 2016-11-09 电子科技大学 基于优化并行码相位搜索的gps捕获电路
CN107247276A (zh) * 2017-07-05 2017-10-13 西安电子科技大学 基于重叠多块补零算法的弱信号捕获方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Semi-bit Frequency Compensation Differential Combining Post-Correlation Processing for GNSS Weak Signal Environment;Zhang Xiang-li et al.;《2017 2nd Asia-Pacific Conference on Intelligent Robot Systems (ACIRS)》;20170720;第129-136页 *
微弱 GPS 信号的翻转位预测捕获算法;李卫斌等;《计算机应用》;20131201;第33卷(第12期);第3473-3476页 *

Also Published As

Publication number Publication date
CN109239744A (zh) 2019-01-18

Similar Documents

Publication Publication Date Title
CN103645483B (zh) 一种弱信号环境下北斗信号捕获方法
JP5015788B2 (ja) 測位信号を受信しつつコヒーレント積分長を長くするための方法及び装置
CN102435999B (zh) Gps接收机基带模块和gps信号捕获及跟踪方法
CN111090108B (zh) 欺骗信号生成方法及装置
US9270323B2 (en) Wireless communication synchronization system
CN101903794A (zh) 用于所接收卫星定位***信号的多路径效应的抑制
CN106230475B (zh) 一种基于Tong检测器伪码多普勒补偿捕获方法及装置
CN113009524B (zh) 长时间相干积分捕获的导航电文比特翻转估计方法及***
CN105553506B (zh) 一种长码扩频信号的快速捕获方法及装置
CN104898136A (zh) 一种北斗二代卫星b1信号的捕获方法及其***
CN103809192A (zh) 一种gnss接收机的动态校正算法
CN104199067A (zh) 多径环境下全球卫星导航***接收机无模糊处理方法
US20050265431A1 (en) Performing an acquisition in a receiver
CN1690726B (zh) 按比例缩减装置及方法,全球定位***同步方法及接收器
US20110050496A1 (en) Energy Domain Based Peak Reconstruction Methods And Apparatuses
CN101614804A (zh) 弱信号下gps导航电文数据位同步的方法及装置
CN109239744B (zh) 一种基于复信号相位的快速抗比特反转快速捕获方法
CN111158026B (zh) 一种boc信号的边峰消除方法
JP5519223B2 (ja) 衛星信号受信装置
JP3787592B2 (ja) 衛星測位方法及び衛星測位システム
KR102155083B1 (ko) 이진 오프셋 반송파 신호를 위한 비 모호 상관 함수 생성 방법 및 장치
US9977131B2 (en) Fast fix using perfect reference in a satellite-based positioning system
CN108627863B (zh) 基于dft的gnss导航电文数据解调方法、装置及gnss接收终端
CN104993844B (zh) 一种频域搜索方法及装置
CN103576169A (zh) 一种多径环境下的cboc调制信号边峰消除方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant