CN104880720A - 一种弱信号卫星导航系统的位同步方法 - Google Patents
一种弱信号卫星导航系统的位同步方法 Download PDFInfo
- Publication number
- CN104880720A CN104880720A CN201510222446.3A CN201510222446A CN104880720A CN 104880720 A CN104880720 A CN 104880720A CN 201510222446 A CN201510222446 A CN 201510222446A CN 104880720 A CN104880720 A CN 104880720A
- Authority
- CN
- China
- Prior art keywords
- value
- bit
- bit synchronization
- coherent integration
- current
- 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.)
- Granted
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/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Signal Processing (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
本发明公开了一种应用于弱信号卫星导航系统的位同步方法,其步骤包括:读取跟踪通道输出的I支路1ms相干积分值大小,并按照1~20循环编号;对当前ms与前1ms的积分值做差,并将结果取绝对值后按照当前ms的编号进行分组累加;循环多个周期后,找出累加值最大的编号,由此可以判断出比特边缘的位置。本发明提出的位同步方法充分利用了相干积分数据提供的幅值和符号信息,相比只利用了相干积分值符号信息的传统位同步方法,能更准确和快速确定比特边缘位置,尤其适用于信号强度变化较快的场合,能实现导航接收机在复杂环境下的快速位同步,从而更快地实现定位功能。
Description
技术领域
本发明属于卫星导航定位技术领域,更具体的,是一种用在弱信号卫星导航定位接收机上的位同步技术。
技术背景
接收端进行卫星定位需要进行几个过程,首先要对卫星进行捕获,通过对可能存在的卫星信号进行频率和码相位搜索,确定真实的卫星信号;然后把对应卫星的码相位和多普勒频率传递到跟踪部分对卫星信号进行持续的跟踪处理;接收机成功跟踪卫星信号后,一方面会输出卫星的码相位和多普勒频率来进行伪距的生成,得到卫星和接收机的空间距离;另一方面则需要进行导航电文的解调,获取卫星的轨道和位置参数,然后结合空间伪距,最终完成位置和速度解算。
在卫星导航接收机实现定位的过程中,导航电文的同步和解调是至关重要的部分,而位同步则是导航电文同步和解调的基础。卫星信号除了载波和导航电文信息外,还经过了伪随机码调制,伪随机码的周期为1ms,而导航电文的码率为50bps,因此一个导航电文比特包含了20个伪随机码周期。接收机通过捕获和跟踪后,便完成了对伪随机码的同步过程,这时候1ms的伪随机码已经对准,接下来则需要找到ms之间的翻转点,也就是导航电文的比特边沿,这个过程就是位同步。只有实现了位同步,接收机才能准确掌握准确的比特位置,从而进行帧同步操作,最终准确解调出导航电文数据,进而完成定位。
位同步的效果不但会影响长时间相干积分的效果,而且会直接影响导航电文的解调的准确度;而位同步的速度也是接收机定位时间的一个关键因素,因此采用一种准确而且高效的位同步方法是卫星导航接收机必须考虑的问题。一般常用的位同步方法包括传统的直方图法,K-P法等。传统的直方图方法通过统计不同位置ms间的符号变化来确定最终比特边沿,具有实现简单,在强信号下准确度高等特点。但在弱信号情况下或者信号强度变化较快的环境下,由于该方法只考虑了ms相干积分值间的符号跳变,所以位同步的准确度会降低很多。而K-P法是通过比较不同位置作为起沿的20ms相干积分累加和,把其中最大的位置作为比特边界。该方法在弱信号环境下准确度相比直方图有较大的提高。但这种方法由于是统计20ms的能量,在强弱信号交替的环境下并不能很准确反映ms之间的具体能量差异。因此如果能采用一种实现简单,能体现ms间积分能量差异的位同步方法具有很高的实际意义。
发明内容
为克服上述现有技术的缺陷与不足,本发明提供一种快速和准确判别比特边沿位置的弱信号卫星导航系统位同步方法,
本发明的技术方案详细如下:
一种弱信号卫星导航系统的位同步方法,用于卫星导航定位接收机上的位同步,其特征在于,包括以下步骤:
(1)待导航定位接收机跟踪稳定后,以20ms为周期,连续读取跟踪通道I支路输出的每1ms的相干积分值,然后按照1~20的顺序进行循环编号;
(2)对于每个相干积分值,用当前1ms相干积分值减去前1ms相干积分值,将结果取绝对值后,按照当前1ms的编号进行分组累加;
(3)以20ms为周期,对每个输入的1ms相干积分值进行步骤(2)的操作,完成多个比特周期,使得比特位连续不变化的概率降到设定的可接受的水平后,得到对应编号1~20的20个积分累加值;
(4)从20个积分累加值中找出最大的积分累加值和相应的信号强度,并记录其对应的编号;把其余的19个积分累加值求平均,得到一个平均噪声累加值和噪声强度;
(5)根据步骤(4)所得到的信号强度和噪声强度的比值判断同步是否成功,如果所得比值满足条件,则认为位同步成功,输出当前最大积分值位置对应的编号即为比特跳变边沿;若条件不满足,则按照(2)~(3)的步骤再处理一个20ms的比特周期,然后按照步骤(4)求出一个新的最大积分累加值和平均噪声累加值,重复步骤(5),直到条件满足,位同步完成。
优选的,在对相干积分值编号为1与前1ms相干积分值的处理为:当前编号为1的时候,由于循环性,其前一个相干积分值为上一个处理周期的第20个编号的相干积分值。
优选的,假设当前1ms相干积分值为 ,前一个1ms的相干积分值为,后一个1ms的相干积分值为,将当前1ms积分值与前1ms积分值求差,然后取绝对值,按照当前1ms编号进行分组累加,然后待下一个1ms的到来,以此循环,当前编号上的积分累加值如公式1所示:
式中,N表示处理周期,K表示当前积分值编号,表示当前编号上的积分累加值。
优选的,首先连续处理多个20ms的比特周期,然后得到积分差值累加和的最大值max(),其对应编号g。然后将剩下的19个ms位置的积分累加和按照公式2取平均,平均噪声累加值计算如所示:
得到一个平均噪声累加值。
优选的,在保证通道连续输出的至少有1ms相干积分值99%概率发生跳变的条件下,首次处理时连续处理至少 6个20ms的比特周期,然后得到积分差值累加和的最大值max()。
进一步优选的,比较所述的积分差值累加和的最大值max()与平均噪声累加值的值,如果max()/,则认为位同步已经成功,记录max()对应的编号g,即为比特跳变边缘;如果max()/,则认为位同步成功条件未达到,重新再处理一个比特周期,再比较max()与的值,如果仍不满足max()/,则再处理一个比特周期,直到条件满足,位同步才完成。
本发明的有益效果:
本发明提供一种弱信号卫星导航系统位同步方法,首先统计跳变位置的幅度差值,再通过查找幅度差值最大的位置快速有效地判断出真实的比特边沿。该方法在发生符号跳变的位置对相邻两ms的积分值求差得到其幅度差值,相比于传统直方图只用符号信息进行统计,在信号强度变化较快,整体强度比较弱的情况下该方法能更有效和准确判别比特边沿位置。
附图说明
图 1 是本发明卫星导航系统位同步方法实现流程示意图。
具体实施方式
本发明的核心思想是通过累积符号跳变位置的能量差来判断真实的比特边沿位置。下面结合附图对该方法的具体实施方式进行详细说明。
图 1 是本发明提出的卫星导航系统位同步方法实现流程示意图,详细说明如下:
卫星导航接收机完成对卫星信号的捕获后,得到粗略的码相位和频偏,然后转入跟踪状态对卫星信号进行更进一步的精确跟踪。接收机对卫星信号的跟踪过程实质上是对卫星信号进行同步以及解调的一个过程,通过PLL以及DLL的相互作用,卫星信号的载波以及CA码得到比较完全的剥离,只剩下导航电文。伪随机码的周期为1ms,而导航电文的比特速率为50bps,因此在接收机中通常按照一个伪随机码周期作为相干积分时间的单位,而20个1ms的积分值则代表导航电文的1个比特长度。
当跟踪稳定后,本地码与接收码完全对齐,此时1ms相干积分值包含了导航电文以及噪声的能量。根据图 1 ,在开始进行位同步之前,首先对每个输出的1ms相干积分值进行连续编号。设跟踪通道I支路输出的1ms相干积分值为,k为当前相干积分值的编号,按照1~20循环编号;而N表示当前处于第N次循环。
第一个循环周期内的第一个1ms相干积分值不做任何处理,直到下一个1ms积分值到来。假设当前1ms积分值为。然后进行图1中步骤二的操作,将当前ms积分值与前1ms积分值求差然后取绝对值,按照当前ms编号进行分组累加。假设当前编号上的能量累加值为,可以由公式1表示:
当前积分值编号k = 1的时候,由于循环性,其前一个积分值为上一个处理周期的第20个编号积分值,也就是当k = 1时,;而当前积分值编号k = 2~19的时候,直接按照上述步骤,对每个1ms积分值进行处理。
在弱信号情况下,卫星信号幅值变化比较快,一般位同步方法需要经过多个20ms周期才能提高位同步的准确率,但比特之间的跳变并不固定,因此处理的周期数N将是一个比较难确定的值。本发明同时提出一种自动调节比特周期数N的方法,可以根据信号强度自动调节处理周期数,以达到在不同信号强度条件下都能快速、准确位同步的目的。
假设比特发生跳变的概率为0.5,在保证99%概率发生跳变的条件下,至少要连续处理6个比特周期,因此本发明的位同步方法为了提高位同步成功率,选取首先连续处理6个20ms的比特周期,然后得到积分差值累加和的最大值max(),其对应编号g。然后将剩下的19个ms位置的积分累加和按照下式取平均,得到一个平均噪声累加值:
比较所述的积分差值累加和的最大值max()与平均噪声累加值的值,则认为位同步已经成功,记录max()对应的编号g,即为比特跳变边缘;如果max()/,则认为位同步成功条件未达到,重新再处理一个比特周期,得到一个新的积分累加和最大值和新的平均噪声累加值,再比较max()与的值,如果仍不满足max()/,则再处理一个比特周期,直到条件满足,位同步才完成。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此。在本发明揭露的技术范围内,任意的变化和替换,都应该涵盖在本发明的保护范围之内。本发明未详细说明部分,属于本领域技术人员公知常识。
Claims (6)
1.一种弱信号卫星导航系统的位同步方法,用于卫星导航定位接收机上的位同步,其特征在于,包括以下步骤:
(1)待导航定位接收机跟踪稳定后,以20ms为周期,连续读取跟踪通道I支路输出的每1ms的相干积分值,然后按照1~20的顺序进行循环编号;
(2)对于每个相干积分值,用当前1ms相干积分值减去前1ms相干积分值,将结果取绝对值后,按照当前1ms的编号进行分组累加;
(3)以20ms为周期,对每个输入的1ms相干积分值进行步骤(2)的操作,完成多个比特周期,使得比特位连续不变化的概率降到设定的可接受的水平后,得到对应编号1~20的20个积分累加值;
(4)从20个积分累加值中找出最大的积分累加值和相应的信号强度,并记录其对应的编号;把其余的19个积分累加值求平均,得到一个平均噪声累加值和噪声强度;
(5)根据步骤(4)所得到的信号强度和噪声强度的比值判断同步是否成功,如果所得比值满足条件,则认为位同步成功,输出当前最大积分值位置对应的编号即为比特跳变边沿;若条件不满足,则按照(2)~(3)的步骤再处理一个20ms的比特周期,然后按照步骤(4)求出一个新的最大积分累加值和平均噪声累加值,重复步骤(5),直到条件满足,位同步完成。
2.根据权利要求1所述的一种弱信号卫星导航系统的位同步方法,其特征在于:在对相干积分值编号为1与前1ms相干积分值的处理为:当前编号为1的时候,由于循环性,其前一个相干积分值为上一个处理周期的第20个编号的相干积分值。
3.根据权利要求2所述的一种弱信号卫星导航系统的位同步方法,其特征在于:假设当前1ms相干积分值为 ,前一个1ms的相干积分值为,后一个1ms的相干积分值为,将当前1ms积分值与前1ms积分值求差,然后取绝对值,按照当前1ms编号进行分组累加,然后待下一个1ms的到来,以此循环,当前编号上的积分累加值如公式1所示:
式中,N表示处理周期,K表示当前积分值编号,表示当前编号上的积分累加值。
4.根据权利要求3所述的一种弱信号卫星导航系统的位同步方法,其特征在于:首先连续处理多个20ms的比特周期,然后得到积分差值累加和的最大值max(),其对应编号g;然后将剩下的19个ms位置的积分累加和按照公式2取平均,平均噪声累加值计算如所示:
得到一个平均噪声累加值。
5.根据权利要求4所述的一种弱信号卫星导航系统的位同步方法,其特征在于:在保证通道连续输出的至少有1ms相干积分值99%概率发生跳变的条件下,首次处理时连续处理至少 6个20ms的比特周期,然后得到积分差值累加和的最大值max()。
6.根据权利要求5所述的一种弱信号卫星导航系统的位同步方法,其特征在于:比较所述的积分差值累加和的最大值max()与平均噪声累加值的值,如果max()/,则认为位同步已经成功,记录max()对应的编号g,即为比特跳变边缘;如果max()/,则认为位同步成功条件未达到,重新再处理一个比特周期,再比较max()与的值,如果仍不满足max()/,则再处理一个比特周期,直到条件满足,位同步才完成。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510222446.3A CN104880720B (zh) | 2015-05-05 | 2015-05-05 | 一种弱信号卫星导航系统的位同步方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510222446.3A CN104880720B (zh) | 2015-05-05 | 2015-05-05 | 一种弱信号卫星导航系统的位同步方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104880720A true CN104880720A (zh) | 2015-09-02 |
CN104880720B CN104880720B (zh) | 2018-10-30 |
Family
ID=53948286
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510222446.3A Active CN104880720B (zh) | 2015-05-05 | 2015-05-05 | 一种弱信号卫星导航系统的位同步方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104880720B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106338746A (zh) * | 2015-12-31 | 2017-01-18 | 深圳艾科创新微电子有限公司 | 北斗卫星导航接收机及提高其跟踪灵敏度的方法、装置 |
CN107807369A (zh) * | 2017-10-11 | 2018-03-16 | 南京航空航天大学 | 一种基于最大似然估计的北斗弱信号位同步方法 |
CN109633711A (zh) * | 2018-12-24 | 2019-04-16 | 长沙北斗产业安全技术研究院有限公司 | 一种超大动态、高灵敏度的扩频测控基带接收方法及装置 |
CN110068842A (zh) * | 2019-05-06 | 2019-07-30 | 西安开阳微电子有限公司 | 一种卫星信号高精度捕获方法 |
CN111487655A (zh) * | 2020-05-15 | 2020-08-04 | 无锡奇芯科技有限公司 | 一种gnss信号位同步方法 |
CN111650619A (zh) * | 2020-07-17 | 2020-09-11 | 广东星舆科技有限公司 | 比特边沿位置的确定方法、装置、可读介质及位同步方法 |
CN111697970A (zh) * | 2020-06-22 | 2020-09-22 | 中国科学技术大学 | 低计算复杂度的周期性弱信号检测装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104483684A (zh) * | 2015-01-05 | 2015-04-01 | 中国科学院重庆绿色智能技术研究院 | 一种快速捕获北斗d1卫星导航系统弱信号的方法 |
-
2015
- 2015-05-05 CN CN201510222446.3A patent/CN104880720B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104483684A (zh) * | 2015-01-05 | 2015-04-01 | 中国科学院重庆绿色智能技术研究院 | 一种快速捕获北斗d1卫星导航系统弱信号的方法 |
Non-Patent Citations (5)
Title |
---|
刘泽民: "三系统卫星导航接收机关键技术研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
李新山等: "一种高动态、弱信号GPS比特同步方法", 《电子与信息学报》 * |
田炳丽等: "GPS软件接收机导航电文提取的研究", 《通信技术》 * |
赵琳等: "极弱信号环境下GPS位同步和载波跟踪技术", 《航空学报》 * |
郑睿等: "一种自适应最大似然GPS比特同步方法", 《系统仿真学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106338746A (zh) * | 2015-12-31 | 2017-01-18 | 深圳艾科创新微电子有限公司 | 北斗卫星导航接收机及提高其跟踪灵敏度的方法、装置 |
CN107807369A (zh) * | 2017-10-11 | 2018-03-16 | 南京航空航天大学 | 一种基于最大似然估计的北斗弱信号位同步方法 |
CN107807369B (zh) * | 2017-10-11 | 2019-11-12 | 南京航空航天大学 | 一种基于最大似然估计的北斗弱信号位同步方法 |
CN109633711A (zh) * | 2018-12-24 | 2019-04-16 | 长沙北斗产业安全技术研究院有限公司 | 一种超大动态、高灵敏度的扩频测控基带接收方法及装置 |
CN110068842A (zh) * | 2019-05-06 | 2019-07-30 | 西安开阳微电子有限公司 | 一种卫星信号高精度捕获方法 |
CN111487655A (zh) * | 2020-05-15 | 2020-08-04 | 无锡奇芯科技有限公司 | 一种gnss信号位同步方法 |
CN111697970A (zh) * | 2020-06-22 | 2020-09-22 | 中国科学技术大学 | 低计算复杂度的周期性弱信号检测装置 |
CN111697970B (zh) * | 2020-06-22 | 2023-04-07 | 中国科学技术大学 | 低计算复杂度的周期性弱信号检测装置 |
CN111650619A (zh) * | 2020-07-17 | 2020-09-11 | 广东星舆科技有限公司 | 比特边沿位置的确定方法、装置、可读介质及位同步方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104880720B (zh) | 2018-10-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104880720A (zh) | 一种弱信号卫星导航系统的位同步方法 | |
CN102486539B (zh) | 一种改进的导航卫星信号跟踪方法 | |
CN105187348A (zh) | 任意速率cpfsk信号定时同步方法 | |
CN104678419B (zh) | 一种卫星导航接收机失锁重捕后快速定位方法 | |
CN104316941B (zh) | 一种基于载波频率辅助相位矢量跟踪方法 | |
CN105974442B (zh) | 卫星导航系统硬件在环实时仿真测试平台及其方法 | |
CN105044737A (zh) | 一种导航接收机中检测观测量有效性的方法 | |
CN107765270A (zh) | 一种基于卡尔曼滤波的卫星导航接收机跟踪环 | |
CN104714241B (zh) | 一种快速gps位同步方法 | |
CN104143997B (zh) | 一种基于扩频体制的多址干扰抗错锁判决方法 | |
CN108055058A (zh) | 载波多普勒及其变化率的高精度测量方法 | |
CN101216549B (zh) | 中短波扩频导航系统距离差观测量提取方法 | |
CN105093242B (zh) | 一种卫星接收机失锁后快速捕获和定位方法 | |
CN105093243A (zh) | 一种基于随机共振算法的gnss载波环路跟踪方法 | |
CN104062668A (zh) | 基于双极化天线的gnss信号捕获、跟踪方法及系统 | |
CN108919311A (zh) | 用于北斗导航系统的抗干扰方法 | |
CN104931980B (zh) | 载波相位测量半周模糊解除方法 | |
CN107367741B (zh) | 用于gnss信号间歇性跟踪的开环卡尔曼方法 | |
CN101247138A (zh) | 一种多相关值定宽拟合伪码跟踪鉴相方法 | |
CN101630000B (zh) | 一种评估干扰信号对gps性能影响的系统 | |
CN103323863B (zh) | Gnss信号自适应快速牵引方法 | |
CN103078727A (zh) | 提高数传接收机帧同步性能的方法 | |
CN104101884B (zh) | 卫星信号的捕获方法及接收机 | |
CN106230477A (zh) | 一种适用于混合扩频的载波跟踪方法 | |
CN102841361B (zh) | 一种全球卫星导航定位系统信号普通测距码快速捕获方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |