CN109633715B - 一种基于并行码相位捕获gps和北斗三代信号的方法 - Google Patents
一种基于并行码相位捕获gps和北斗三代信号的方法 Download PDFInfo
- Publication number
- CN109633715B CN109633715B CN201811647399.7A CN201811647399A CN109633715B CN 109633715 B CN109633715 B CN 109633715B CN 201811647399 A CN201811647399 A CN 201811647399A CN 109633715 B CN109633715 B CN 109633715B
- Authority
- CN
- China
- Prior art keywords
- data
- gps
- beidou
- generation
- code
- 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
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)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种基于并行码相位捕获GPS和北斗三代信号的方法,包括以下步骤:1)截获卫星信号,生成GPS L1C/A和北斗三代B1C主码组合;2)对主码组合进行FFT运算,得到主码FFT数据;3)读取10ms GPS L1C/A和B3 B1C的组合卫星数据,对其进行FFT运算,得到采样序列FFT数据;4)对采样序列FFT数据和主码FFT数据之一取共轭,而后与另外一个相乘得到IFFT;5)利用IFFT判断捕获类型,若捕获类型包括GPS L1C/A信号或北斗三代B1C信号,输出数据至相应的跟踪通道;否则,重复步骤1)‑3)。该方法能够减短了捕获的时间,减短捕获序列,且大大提高了捕获成功概率。
Description
技术领域
本发明属于信号处理技术,尤其是一种基于并行码相位捕获GPS和北斗三代信号的方法。
背景技术
目前,民用的导航芯片大都面向现行GPS的民用信号L1C/A进行设计,该信号也位于1575.42MHz。由于GPS的深入人心和广泛应用,各家卫星导航***为了推广各自的***,都尽量向GPS靠齐。因此,为了兼顾到导航芯片厂商更改便利和为用户提供更多更好的民用服务,在未来的发展计划中,各家都在1575.42MHz这个频点采用相似的信号。北斗***和GPS***是基于码分多址的扩频通信***,其基本原理是将导航电文与伪随机码(PRN)共同调制在载波上。由于伪随机码良好的自相关特性,接收机可以在接收的信号中分辨出这是哪一颗卫星的信号。接收机的第一步便是捕获,即确定卫星编号,大致的多普勒频移及码相位偏移。只关注第一项,即确定卫星编号。PRN码的特性可以概括为相同的码互相关为1,否则趋近于零;接收机便可以根据这个原理,产生每一颗卫星的编码分别与收到的信号做互相关运算,从第一颗试到最后一颗卫星,总可以找到合适的PRN,也就是发射该信号的卫星。
随着北斗三代的推广,现在通讯设备大多兼容GPS L1C/A和北斗三代B1C信号,极大提高了捕获效率,提高多模定位的稳定性。但现有联合捕获方案技术方案一般是将数据分开并列处理的模式。现有技术存在加长捕获序列、延长捕获的时间,并且降低了成功捕获概率。
发明内容
本发明的目的在于克服上述现有技术的缺点,提供一种基于并行码相位捕获GPS和北斗三代信号的方法。
为达到上述目的,本发明采用以下技术方案予以实现:
一种基于并行码相位捕获GPS和北斗三代信号的方法,包括以下步骤:
1)截获卫星信号,生成GPS L1C/A和北斗三代B1C主码组合;
2)对主码组合进行FFT运算,得到主码FFT数据;
3)读取10ms GPS L1C/A和B3B1C的组合卫星数据;
4)对读取组合卫星数据进行FFT运算,得到采样序列FFT数据;
5)对采样序列FFT数据和主码FFT数据之一取共轭,而后与另外一个相乘得到IFFT;
6)利用IFFT判断捕获类型,若捕获类型包括GPS L1C/A信号或北斗三代B1C信号,输出数据到相应的跟踪通道;
否则,重复步骤1)-3)。
进一步的,步骤1)包括:
101)产生GPS L1C/A主码和北斗三代B1C主码;
102)采用排列模式或累加模式组合GPS L1C/A主码和北斗三代B1C主码,生成GPSL1C/A和北斗三代B1C主码组合。
进一步的,步骤4)进行FFT运算之前还包括:
对读取的数据去中频,将采样数据频率降至原来的1/M,生成采样序列FFT数据。
进一步的,去中频运算过程如下:
I(t)=d(t)*cos(2πfIFt) (1)
Q(t)=d(t)*sin(2πfIFt) (2)
其中,其中,t为连续时间;I(t)为I支路数据;Q(t)为与I支路正交的Q支路数据;fIF为中频频率。
进一步的,步骤6)中的判断捕获类型具体为:
对IFFT取模数,得到模数数据;
将模数数据与GPS阈值对比,判断是否捕获GPSL1C/A;
将模数数据与B1C阈值对比,判断是否捕获北斗卫星B1C。
进一步的,将模数数据与L1C/A阈值对比,当模数数据中至少三个数值超过L1C/A阈值,且每个数值的时间差均为整数倍的ms时,判断捕获GPSL1C/A信号。
进一步的,将模数数据与B1C阈值对比,当模数数据中有一个或两个数值超过B1C阈值,即为捕获北斗卫星B1C的信号。
与现有技术相比,本发明具有以下有益效果:
一种基于并行码相位捕获GPS和北斗三代信号的方法,通过将通道之前的抓取,由择一抓取改为并行抓取,减短了捕获的时间,而对采样到的数据进行降频,能够减短捕获序列,在终端与卫星信号通讯时,多追踪通道同时运行,大大提高了捕获成功概率。
附图说明
图1为本发明的GPS对应多普勒结果图,此时码分辨率采用1/4模式且组合码采用累加模式;
图2为本发明的北斗三代B1C信号对应多普勒结果图,此时码分辨率采用1/4模式且组合码采用累加模式;
图3为本发明的GPS和北斗三代组合多普勒结果图,此时码分辨率采用1/4模式且组合码采用累加模式;
图4为本发明提供的实施例1的多普勒结果图;
图5为本发明提供的实施例2的多普勒结果图;
其中:B3表示北斗三代;横坐标为码片,单位为1/4码片;纵坐标为功率值。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、***、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
下面结合附图对本发明做进一步详细描述:
一种基于并行码相位捕获GPS和北斗三代信号的方法,包括:
1)截获卫星信号,生成GPS L1C/A和北斗三代B1C主码组合;
2)对主码组合进行FFT运算,得到主码FFT数据;
3)读取10ms GPS L1C/A和B3 B1C的组合卫星数据;
4)对读取组合卫星数据进行FFT运算,得到采样序列FFT数据;
5)对采样序列FFT数据和主码FFT数据之一取共轭,而后与另外一个相乘得到IFFT;
6)利用IFFT判断捕获类型,若捕获类型包括GPS L1C/A信号或北斗三代B1C信号,输出数据;
否则,重复步骤1)-3)。
其中,步骤1)包括:
101)产生GPS L1C/A主码;C/A信号主码的码速率为1.023Mcps,码长为1023。C/A码发生器包含两个移位寄存器G1和G2,两个移位寄存器都是10位,可产生长度为1023的序列,当G1和G2寄存器的多项式产生最大长度即为1023的码时,它们模二相加即可产生1023个码片的C/A码序列输出。根据不同的卫星输入,选取不同的G2寄存器的相位输出,通过查找表得到,最后得到不同的C/A码序列输出。
102)产生北斗三代B1C主码;B1C信号主码的码速率为1.023Mcps,码长为10230,由长度为10243的Weil码通过截断获得。Weil码通过两个legendre序列模二相加得到。不同的卫星,两个legendre序列之间的相位差不同,Weil码截断的初始值不同,通过查找表得到,最后得到不同的B1C主码码序列输出。
103)采用排列模式或累加模式组合GPS L1C/A主码和北斗三代B1C主码,生成GPSL1C/A和北斗三代B1C主码组合,下面为三个实施例。
1.当采用码分辩率为1/4模式时,采用累加模式,如下:
10ms B1C+GPS | 6ms 0 |
添加6ms的0,总长度65536。10ms的GPS L1C/A和北斗三代B1C主码直接累加,数据如下表:
10ms data | 6ms 0 |
或者数据
16ms data |
2.当采用码分辩率为1/4模式时,并列模式,如下
GPS L1C/A和北斗三代B1C主码平行放置,同时会添加(补偿至65536)64位的0.总长度65536,数据如下:
3.当采用码分辩率为1/3模式时,两者采用累加模式,如下表:
10ms B1C+GPS |
不需要添加0.总长度32768.
数据如下:
10ms data |
步骤4)进行FFT运算之前还包括:对读取的数据去中频,将采样数据频率降至原来的1/M,生成采样序列FFT数据;当采用四分之一码片模式时,M为16;当采用三分之一码片模式时,M为21。
其中,去中频运算过程如下:
I(t)=d(t)*cos(2πfIFt) (1)
Q(t)=d(t)*sin(2πfIFt) (2)
其中,t为连续时间;I(t)为I支路数据;Q(t)为与I支路正交的Q支路数据;fIF为中频频率;
采用定点运算过程如下:
I(k)=d(k)*cos_table(index_k) (3)
Q(k)=d(k)*sin_table(index_k) (4)
其中,k为离散时间;I(k)为I支路数据;Q(k)为与I支路正交的Q支路数据;index_k:离散时间对应cos表索引。
参见图1-3,图1、图2和图3分别为本发明的GPS对应多普勒结果图、北斗三代B1C信号对应多普勒结果图及GPS和北斗三代组合多普勒结果图,此时码分辨率采用1/4模式且组合码采用累加模式;其中,横坐标为码片,单位为1/4码片;纵坐标为功率值;三个图为三维图在横坐标和纵坐标组成平面的投影,与该平面垂直的第三个轴为多普勒频偏轴,从图中可以看出,当GPS L1C/A和北斗三代B1C的多普勒频偏相同时,就会出现在该平面重叠的情况,需要分别筛选出GPS L1C/A和B3B1C卫星的信号。
判断捕获类型具体如下:
B1C通过与B1C阈值对比,判断是否捕获了B1C卫星,在IFFT取模之后的数据中,当存在一个或者两个值超过阈值,即可判断捕获成功B1C卫星,将结果输入到B1C跟踪通道中。
GPS通过与GPS阈值对比,判断是否捕获GPS卫星,在IFFT取模之后的数据中,至少需要3个值超过阈值,并且三个值时间差整ms,当为1/4模式,索引差别是4092的整数倍,才可判断捕获成功GPS卫星,将结果输入到GPS跟踪通道中。
实施例1
捕获GPS L1C/A PRN 3和北斗三代PRN 21,采用主码累加模式,10ms数据,使用的比例阈值为25,最大值与均值的比例:max/mean。
参见图4,图4为本发明提供的实施例1的捕获图;图上GPS的多普勒是-1437.5Hz,最大比例值是276.5831,码相位值是86(1/4码片计数值,以下类同),其他衍生(以4092为间隔)出去的值为[276.5831,223.1504,188.1428,112.0199,93.9620,65.7849,40.8874,21.6674,9.3089,21.7136,57.9774,73.7619,94.6974,159.1306,165.3054,225.9995](按照索引排序),大于25的值的数目为13,大于3,判断GPS PRN 3捕获成功,以此进入跟踪模式。
图中B3B1C的多普勒是-1812.5Hz,码相位值是25943,最大比例值是81.8920,次峰值的比例是29.9011,但是索引不符合相差4092的规律,所以大于25的值的数目为1,小于3,判断北斗三代PRN 21捕获成功,进入跟踪模式。
实施例2
捕获GPS L1C/A PRN 3和北斗三代PRN 21,采用主码累加模式,取10ms数据,使用的是比例阈值为25,最大值与均值的比例:max/mean。
参见图5,图5为本发明提供的实施例2的捕获图;图中GPS的多普勒是3187.5Hz,最大比例值是243.2454,码相位值是2420,其他衍生(以4092为间隔)出去的值,按照索引排序如下,[232.7543,181.9386,169.9988,101.9293,81.9932,54.7541,29.3254,16.8035,11.7245,34.1128,54.3310,83.9683,91.7543,144.9247,189.2938,242.7930],其中大于25的值的数目为14,远大于3,判断这是GPS PRN 23捕获成功,进入跟踪模式。
图中B3B1C的多普勒是-2812.5Hz,码相位值是30345,最大比例值是26.4462,次峰不存在,大于25的值的数目为1,小于3,可以判断这是北斗三代PRN 22捕获成功,以此进入跟踪模式。
以上内容仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明权利要求书的保护范围之内。
Claims (3)
1.一种基于并行码相位捕获GPS和北斗三代信号的方法,其特征在于,包括以下步骤:
1)截获卫星信号,生成GPS L1C/A和北斗三代B1C主码组合;
101)产生GPS L1C/A主码和北斗三代B1C主码;
102)采用排列模式或累加模式组合GPS L1C/A主码和北斗三代B1C主码,生成GPS L1C/A和北斗三代B1C主码组合;
2)对主码组合进行FFT运算,得到主码FFT数据;
3)读取10ms GPS L1C/A和B3 B1C的组合卫星数据;
4)对读取组合卫星数据进行FFT运算,得到采样序列FFT数据;
5)对采样序列FFT数据和主码FFT数据之一取共轭,而后与另外一个相乘得到IFFT;
6)利用IFFT判断捕获类型,若捕获类型包括GPS L1C/A信号或北斗三代B1C信号,输出数据到相应的跟踪通道;
否则,重复步骤1)-3);
步骤6)中的判断捕获类型具体为:
对IFFT取模数,得到模数数据;
将模数数据与GPS阈值对比,判断是否捕获GPSL1C/A;
将模数数据与B1C阈值对比,判断是否捕获北斗卫星B1C;
将模数数据与L1C/A阈值对比,当模数数据中至少三个数值超过L1C/A阈值,且每个数值的时间差均为整数倍的ms时,判断捕获GPSL1C/A信号;
将模数数据与B1C阈值对比,当模数数据中有一个或两个数值超过B1C阈值,即为捕获北斗卫星B1C的信号。
2.根据权利要求1所述的基于并行码相位捕获GPS和北斗三代信号的方法,其特征在于,步骤4)进行FFT运算之前还包括:
对读取的数据去中频,将采样数据频率降至原来的1/M,生成采样序列FFT数据。
3.根据权利要求2所述的基于并行码相位捕获GPS和北斗三代信号的方法,其特征在于,去中频运算过程如下:
I(t)=d(t)*cos(2πfIFt) (1)
Q(t)=d(t)*sin(2πfIFt) (2)
其中,t为连续时间;I(t)为I支路数据;Q(t)为与I支路正交的Q支路数据;fIF为中频频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811647399.7A CN109633715B (zh) | 2018-12-29 | 2018-12-29 | 一种基于并行码相位捕获gps和北斗三代信号的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811647399.7A CN109633715B (zh) | 2018-12-29 | 2018-12-29 | 一种基于并行码相位捕获gps和北斗三代信号的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109633715A CN109633715A (zh) | 2019-04-16 |
CN109633715B true CN109633715B (zh) | 2021-04-30 |
Family
ID=66054879
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811647399.7A Active CN109633715B (zh) | 2018-12-29 | 2018-12-29 | 一种基于并行码相位捕获gps和北斗三代信号的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109633715B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110068841A (zh) * | 2019-05-06 | 2019-07-30 | 西安开阳微电子有限公司 | 一种卫星信号高灵敏度快速捕获方法 |
CN113093233B (zh) * | 2021-04-12 | 2024-04-12 | 湖南国科微电子股份有限公司 | 一种l6频段卫星信号的同步译码方法、装置、设备及介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101278207A (zh) * | 2005-07-01 | 2008-10-01 | 欧洲太空署 | 用于卫星导航***的扩频码 |
CN101634696A (zh) * | 2008-07-23 | 2010-01-27 | 中国科学院国家天文台 | Gps/caps双模组合接收机的信号处理方法 |
CN101639529A (zh) * | 2009-09-02 | 2010-02-03 | 东莞市泰斗微电子科技有限公司 | 一种兼容gps和北斗2导航***的卫星信号捕获方法 |
CN102116866A (zh) * | 2009-12-31 | 2011-07-06 | 和芯星通科技(北京)有限公司 | 整周载波的gps p和/或y码信号的跟踪方法和装置 |
CN103091687A (zh) * | 2013-01-08 | 2013-05-08 | 南京航空航天大学 | 北斗接收机多卫星信号联合捕获方法 |
CN105676236A (zh) * | 2015-12-30 | 2016-06-15 | 航天恒星科技有限公司 | 一种导航信号捕获方法及*** |
KR20170127219A (ko) * | 2016-05-11 | 2017-11-21 | 주식회사 두시텍 | Gps/갈릴레오 복합수신기 베이스밴드 트랙킹 모듈 성능 개선 장치 |
-
2018
- 2018-12-29 CN CN201811647399.7A patent/CN109633715B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101278207A (zh) * | 2005-07-01 | 2008-10-01 | 欧洲太空署 | 用于卫星导航***的扩频码 |
CN101634696A (zh) * | 2008-07-23 | 2010-01-27 | 中国科学院国家天文台 | Gps/caps双模组合接收机的信号处理方法 |
CN101639529A (zh) * | 2009-09-02 | 2010-02-03 | 东莞市泰斗微电子科技有限公司 | 一种兼容gps和北斗2导航***的卫星信号捕获方法 |
CN102116866A (zh) * | 2009-12-31 | 2011-07-06 | 和芯星通科技(北京)有限公司 | 整周载波的gps p和/或y码信号的跟踪方法和装置 |
CN103091687A (zh) * | 2013-01-08 | 2013-05-08 | 南京航空航天大学 | 北斗接收机多卫星信号联合捕获方法 |
CN105676236A (zh) * | 2015-12-30 | 2016-06-15 | 航天恒星科技有限公司 | 一种导航信号捕获方法及*** |
KR20170127219A (ko) * | 2016-05-11 | 2017-11-21 | 주식회사 두시텍 | Gps/갈릴레오 복합수신기 베이스밴드 트랙킹 모듈 성능 개선 장치 |
Non-Patent Citations (2)
Title |
---|
"A simultaneous multiple Beidou signal acquisition algorithm for a software-based GNSS receiver";Fei Xie等;《Optik》;20160228;第127卷(第4期);第1607-1614页 * |
"FFT快速捕获算法在GPS C/A码与P(Y)码中的应用";张博等;《电讯技术》;20081231;第48卷(第12期);第34-38页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109633715A (zh) | 2019-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2294216C (en) | Receiving a spread spectrum signal | |
CN101174853B (zh) | 起始胞元检索算法 | |
Alaqeeli et al. | Real-time acquisition and tracking for GPS receivers | |
WO1998059427A1 (en) | Processing a spread spectrum signal in a frequency adjustable system | |
AU8374298A (en) | Acquiring a spread spectrum signal | |
CN110071738A (zh) | 基于多伪码支路接收机的扩频信号解扩跟踪方法 | |
CN102162852A (zh) | 大规模频率偏移下微弱gnss信号捕获方法及其*** | |
CN109633715B (zh) | 一种基于并行码相位捕获gps和北斗三代信号的方法 | |
US20130039392A1 (en) | Satellite Signal Acquisition | |
CN104898136A (zh) | 一种北斗二代卫星b1信号的捕获方法及其*** | |
CN100581068C (zh) | 一种多星座卫星信号的混合串并行相关捕获方法 | |
CN102798871B (zh) | 基于伪码重构的伪码捕获方法及装置 | |
EP0992133A1 (en) | Bandpass processing of a spread spectrum signal | |
CN1576805A (zh) | 定位全球定位***相关峰值信号的方法和*** | |
CN111257913B (zh) | 北斗卫星信号捕获方法以及装置 | |
Manandhar et al. | GPS signal acquisition and tracking-An Approach towards development of Software-based GPS Receiver | |
CN105891853A (zh) | 信号捕获方法及装置 | |
JP3852533B2 (ja) | 初期捕捉回路 | |
JP2005513510A (ja) | 信号の受信と捕捉のための方法及び装置 | |
CN108957492B (zh) | 一种gps的l1c/a和l1c联合捕获方法 | |
CN104849730B (zh) | 伽利略卫星导航***e1频点信号无模糊抗多径方法 | |
CN103869343B (zh) | 捕获gps中l2c信号的方法及装置 | |
CN112600588B (zh) | 一种非整数周期扩频通信捕获方法 | |
CN101499844B (zh) | 用于wcdma***下行接收机的频率同步方法 | |
CN112764063A (zh) | 一种实现捕获处理的方法及接收机 |
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 |