CN104749588A - 一种与北斗系统实时同步生成伪卫星信号的方法 - Google Patents

一种与北斗系统实时同步生成伪卫星信号的方法 Download PDF

Info

Publication number
CN104749588A
CN104749588A CN201510145678.3A CN201510145678A CN104749588A CN 104749588 A CN104749588 A CN 104749588A CN 201510145678 A CN201510145678 A CN 201510145678A CN 104749588 A CN104749588 A CN 104749588A
Authority
CN
China
Prior art keywords
satellite
pseudolite
signal
real
time
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
Application number
CN201510145678.3A
Other languages
English (en)
Other versions
CN104749588B (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.)
Xi Beidou technology development (Beijing) Co., Ltd.
Original Assignee
Beijing Hua Yunzhi Joins Science And Technology Ltd
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 Beijing Hua Yunzhi Joins Science And Technology Ltd filed Critical Beijing Hua Yunzhi Joins Science And Technology Ltd
Priority to CN201510145678.3A priority Critical patent/CN104749588B/zh
Publication of CN104749588A publication Critical patent/CN104749588A/zh
Application granted granted Critical
Publication of CN104749588B publication Critical patent/CN104749588B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/03Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/10Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing dedicated supplementary positioning signals
    • G01S19/11Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing dedicated supplementary positioning signals wherein the cooperating elements are pseudolites or satellite radio beacon positioning system signal repeaters

Abstract

本发明涉及一种与北斗系统实时同步生成伪卫星信号的方法,其步骤:获取当前的实际卫星系统星历参数,同时实现将本地信号与实际北斗时的时间同步,并判断是否实现时间同步;根据当前的可见星和用户位置关系利用选星方法进行伪卫星星号选择;根据存储的星号,将该星号对应的伪卫星信号进行仿真计算,生成与实际星历同步的电文参数和观测量参数;输出伪卫星信号的伪距参数、多普勒参数和误差参数,判断当前伪卫星是否已经不再满足选星要求;计算选出的备用卫星的电文参数和观测量参数,并判断当前时刻是否应该进行换星;满足则进行换星;反之,循环进行伪卫星的生成,实现与北斗系统的同步。本发明能够广泛应用于通用接收机,能够和在轨卫星联合定位。

Description

一种与北斗系统实时同步生成伪卫星信号的方法
技术领域
本发明涉及一种位置服务、卫星导航和物联网应用领域,特别是关于一种与北斗系统实时同步生成伪卫星信号的方法。
背景技术
导航卫星信号生成是为基于位置服务、卫星导航和物联网的应用终端提供高可靠性、高精度、可重复应用的测试信号和测试手段。基于通用接收机的大范围的卫星信号遮挡区域的高精度定位是当前导航定位领域研究的一个热点和难点,由于范围较大(达到公里量级),WiFi定位、移动通信基站定位等技术难以应用,最有效可行的方式伪卫星定位方式。
现在大部分基于伪卫星的定位技术都需要对接收机进行改进,主要的原因是卫星导航定位系统里卫星位置的星历都是采用的开普勒六参数的方式发布,接收机根据开普勒六参数及其修正参数实时计算出卫星的位置,而伪卫星的位置一般不能用开普勒六参数来表示,一般都采用地球协仪坐标的方式来表示,所以,普通的接收机无法解析这部分参数,也就无法获得卫星的位置。对于接收机来说,需要开辟专门的通道来解析伪卫星的电文信息,这也就限制了通用接收机利用伪卫星来进行定位。特别是对于现在广泛使用的移动终端的接收机来说,无法利用这种伪卫星的信号进行定位,也就限制了伪卫星技术的应用。
发明内容
针对上述问题,本发明的目的是提供一种与北斗系统实时同步生成伪卫星信号的方法,该方法能够应用于通用接收机,并且能够实现和在轨卫星联合定位。
为实现上述目的,本发明采取以下技术方案:一种与北斗系统实时同步生成伪卫星信号的方法,其包括以下步骤:1)由导航接收机接收实际卫星信号,并由导航系统控制器向导航接收机发送多个实际卫星系统星历接收指令,导航接收机获取当前的实际卫星系统星历参数,同时通过1PPS实现将本地信号与实际北斗时的时间同步,并判断是否实现时间同步,若实现时间同步后,将同步标志设置为“时间已同步”,进入步骤2);反之继续通过1PPS实现时间同步;2)根据当前的可见星和用户位置关系利用选星方法进行伪卫星星号选择,并将选择的星号进行存储;3)根据存储的星号,将该星号对应的伪卫星信号进行仿真计算,生成与实际星历同步的电文参数和观测量参数;4)通过已有的信号生成机制将伪卫星信号的伪距参数、多普勒参数和误差参数转换后并输出,信号播发后每隔1秒进行一次判断,判断当前伪卫星是否已经不再满足选星要求,如果不满足要求则重新输出伪卫星信号;满足则根据选星方法选出1颗备用卫星;5)计算选出的备用卫星的电文参数和观测量参数,并判断当前时刻是否应该进行换星;如果满足换星时刻,则进行换星处理,关闭当前卫星信号,启动新的卫星,换星过程完成;反之,返回步骤4),循环进行伪卫星的生成,实现了与北斗系统的同步。
所述步骤1)中,所述通过1PPS实现将本地信号与实际北斗时的时间同步包括以下步骤:(1)设置时钟同步结构,该时钟同步结构包括授时型导航接收机、接收天线、带有本地原子钟的时钟驯服模块和伪卫星信号发生器;(2)授时型导航接收机通过接收天线接收实际卫星信号,并读取实际卫星系统UTC时;授时型导航接收机将实际卫星系统UTC时进行读秒处理后,利用1PPS信号将时间信息传输至时钟驯服模块,同时通过串口将该时间信息传输至伪卫星信号发生器;通过时钟驯服模块同步生成10MHz时钟信号和稳定的1PPS信号;(3)授时型导航接收机通过时钟驯服模块输出的10MHz时钟信号作为伪卫星信号发生器的时钟基准信号,授时型导航接收机通过时钟驯服模块输出的下一个1PPS信号作为伪卫星信号发生器的外触发信号;(4)伪卫星信号发生器将读取的实际卫星UTC时间转换后产生伪卫星系统时间,伪卫星信号发生器工作在外触发模式下,读取实际卫星系统UTC时间并通过时钟驯服模块获取到稳定的1PPS信号后,判决是否具备从下一秒开始和实际时间对齐的条件,判决条件成立后开始计算输出下一秒的伪卫星信号,并同步利用1PPS信号和10MHz将伪卫星系统本地时间与UTC时间进行对齐,实现同步。
所述步骤4)中,所述伪卫星的选星原则如下:(1)用户范围内难以正常接收的实际卫星,即不在用户视野内的实际卫星;(2)实际卫星—伪卫星—用户所构成的转发线路相对于直线线路的误差要尽可能的小;(3)伪卫星与现有卫星组成的联合星座的GDOP值小。
所述GDOP值计算方法如下:(1)针对伪卫星的不同位置,计算所选伪卫星的星历:以用户设备接收天线位置为原点,建立测量坐标系,在该坐标系内,计算实际卫星的高度角和方位角;(a)计算实际卫星在测量坐标系的位置 式中,为实际卫星在ECEF坐标系的坐标;为用户在ECEF坐标系的坐标;(L,B)为用户的经度、纬度;RY、RX和RZ分别是绕Y轴、X轴和Z轴的旋转矩阵;(b)根据实际卫星的位置计算实际卫星的高度角E、方位角Az:
E = arctg ( y sc / D ) Az = arccos ( x sc / D ) + 0 , z sc &GreaterEqual; 0 , &pi; , z sc < 0
式中,分别是用户坐标系下卫星的x、y和z轴坐标;(2)根据实际卫星的高度角E和方位角Az计算GDOP值:
GDOP = D 11 + D 12 + D 13 + D 14 .
所述步骤2)、步骤4)中,所述伪卫星选星方法为:收集实际卫星参数,根据用户和障碍物的位置关系区分可见区域,在不可见区域内选择仰角最优的实际卫星,提取相应实际卫星星号,完成选星。
所述步骤5)中,所述换星过程如下:(1)伪卫星信号发生器选好符合要求的不可见区域的实际卫星并播发这颗卫星的信号;(2)计算不可见区域的仰角范围,实时进行伪卫星仰角值计算,根据卫星运动轨迹推算该卫星走出不可见区域的时间;(3)如果计算伪卫星离开的时间小于10分钟,重新进行选星过程,并选出当前条件下更优的伪卫星;(4)收集新的伪卫星星历,同时计算出该卫星的实时电文参数和观测量参数,包括伪距参数、多普勒参数和误差参数;(5)新的伪卫星的卫星数据准备好之后,通知伪卫星信号发生器进行换星操作,关闭当前卫星信号,并在下一个1PPS发送新的伪卫星的卫星信号;(6)换星过程完成。
本发明由于采取以上技术方案,其具有以下优点:1、本发明实现了本地时间与实际北斗时的时间同步,时间同步指伪卫星信号输出的系统时间与实际导航系统时之间的时间同步精度,本发明将系统时间同步到GPS时或者北斗时,然后再实现与UTC时间的转换;根据系统的应用要求,将伪卫星系统时间与UTC时间之间的时间差小于20ns,同时为实现伪卫星信号发生器内部时钟精度,保证伪卫星信号发生器本身的系统时间同步精度小于10ns。2、本发明实现了随时间变化的卫星选星策略。伪卫星选星是指为了提高系统的可用性,将伪卫星所模拟的实际卫星进行合理设计,使服务区域内的导航接收机能够在最大范围内实现折线误差小和良好的DOP值。根据用户的位置和所接收的实际卫星的位置,选择用于转发的卫星编号。3、本发明实现了与实际星历同步的电文生成。星历同步是指伪卫星系统的星历参数与实际卫星导航系统的星历参数保持一致。需要解决的问题是在伪卫星运行过程中长期稳定输出可靠的星历,并且需要实现实时星历的更新,保持每1个小时更新一次系统星历。保证伪卫星播发的电文信息中相关的位置信息、钟差信息和其他相关信息与现有的导航系统机制保持一致,是保证伪卫星播发的电文信息能够在现有通用接收机中正常解算,使接收机实现增强定位。4、本发明对于一些有遮挡,也能够收到1到3颗卫星,能够不干扰正常接收的卫星信号,同时能够通过卫星信号生成后发射其不能接收的卫星信号,使得导航接收机能够接收到大于4颗卫星信号,进行定位解算;同时,产生的伪卫星信号增强了终端的定位成功率和定位精度,能够应用于通用接收机,并且能够实现和在轨卫星联合定位。本发明能够广泛应用于通用接收机,并且能够和在轨卫星联合定位。
附图说明
图1是本发明的整体流程示意图;
图2是本发明的时钟同步示意图;
图3是本发明的选星示意图;
图4是本发明的选星流程示意图;
图5是本发明的实际卫星高度角和方位角示意图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示,本发明提供一种与北斗系统实时同步生成伪卫星信号的方法,与只能根据电文和星历离线产生卫星导航模拟信号的方法不同,本发明包括以下步骤:
1)由导航接收机接收实际卫星信号,并由导航系统控制器向导航接收机发送多个实际卫星系统星历接收指令,导航接收机获取当前的实际卫星系统星历参数,同时通过1PPS实现将本地信号与实际北斗时的时间同步,并判断是否实现时间同步,若实现时间同步后,将同步标志设置为“时间已同步”,进入步骤2);反之继续通过1PPS实现时间同步。
2)根据当前的可见星和用户位置关系利用选星方法进行伪卫星星号选择,并将选择的星号进行存储,以供后续使用。
3)根据存储的星号,将该星号对应的伪卫星信号通过已有仿真算法进行仿真计算,生成与实际星历同步的电文参数和观测量参数,以提供产生伪卫星信号的伪距参数、多普勒参数和误差参数等。
4)通过已有的信号生成机制将伪卫星信号的伪距参数、多普勒参数和误差参数转换后并输出,信号播发后每隔1秒进行一次判断,判断当前伪卫星是否已经不再满足选星要求,如果不满足要求则重新输出伪卫星信号;满足则根据选星方法选出1颗备用卫星。
5)计算选出的备用卫星的电文参数和观测量参数,并判断当前时刻是否应该进行换星;如果满足换星时刻,则进行换星处理,关闭当前卫星信号,启动新的卫星,换星过程完成;反之,返回步骤4),循环进行伪卫星的生成,实现了与北斗系统的同步。
上述步骤1)中,如图2所示,由于真实的GNSS系统(全球导航卫星系统)时标信号是由导航系统的原子钟组进行维持的,普通用户无法获取导航系统的原子时钟,因此本发明进行时间同步时采用通过GNSS驯服的原子钟,在地面取得一个GNSS系统时标信号,实现与UTC时间的同步,该信号用作伪卫星系统的基准时钟,该基准时钟与GNSS系统时的同步精度受限于原子钟的实现精度,本发明采用高性能的授时型导航接收机与本地原子钟联合实现高性能的本地时钟。即通过1PPS实现将本地信号与实际北斗时的时间同步具体包括以下步骤:
(1)设置时钟同步结构,该时钟同步结构包括授时型导航接收机、接收天线、带有本地原子钟的时钟驯服模块和伪卫星信号发生器。
(2)授时型导航接收机通过接收天线接收实际卫星信号,并读取实际卫星系统UTC时。授时型导航接收机将实际卫星系统UTC时进行读秒处理后,利用1PPS信号将时间信息传输至时钟驯服模块,同时通过串口将该时间信息传输至伪卫星信号发生器;通过时钟驯服模块同步生成10MHz时钟信号和稳定的1PPS信号。
(3)授时型导航接收机通过时钟驯服模块输出的10MHz时钟信号作为伪卫星信号发生器的时钟基准信号,授时型导航接收机通过时钟驯服模块输出的下一个1PPS信号作为伪卫星信号发生器的外触发信号。
(4)伪卫星信号发生器将读取的实际卫星UTC时间转换后产生伪卫星系统时间,伪卫星信号发生器工作在外触发模式下,读取实际卫星系统UTC时间并通过时钟驯服模块获取到稳定的1PPS信号后,判决是否具备从下一秒开始和实际时间对齐的条件,判决条件成立后开始计算输出下一秒的伪卫星信号,并同步利用1PPS信号和10MHz将伪卫星系统本地时间与UTC时间进行对齐,实现同步。
上述步骤(4)中,伪卫星本地时间与UTC时间同步精度取决于伪卫星信号发生器内部采样时钟频率和导航观测量和电文计算时产生的时间偏差,导航观测量和电文计算时间偏差采用的是数学计算方式,内部字长为64位浮点数计算,产生的误差可以忽略,目前内部采样时间频率要满足时间小于10ns的要求,需要将系统时钟提高到100MHz以上,为了与导航系统实现基准频率相同,本发明伪卫星信号发生器内部采样时钟为327.36MHz。
上述步骤4)中,由于伪卫星选择的基本原则是尽可能少干扰可以正常接收的实际卫星;同时要求伪卫星产生的伪距误差尽可能小;所以,伪卫星需要对所模拟的实际卫星进行选择。如图3所示,对于应用区域来说,可根据周边的环境,划定“可视区域”、“非可视区域”等;为了不干扰可视区域实际卫星的正常接收,所以伪卫星只能模拟“非可视区域”的实际卫星;同时,对于“非可视区域”的实际卫星,如果从实际卫星经伪卫星转发到用户的折线角度过大,会造成伪距误差较大,虽然可以通过修正的方法来修正,但是由于实际卫星的运动速度较快,修正的实时性难度较大;所以,尽可能选择折线角度较小的区域的实际卫星。具体地,伪卫星的选星原则如下:
(1)用户范围内难以正常接收的实际卫星,也就是不在用户视野内的实际卫星;例如:
根据应用区域的大小、形状以及伪卫星布设高度的不同,分别设置伪卫星所模拟的实际卫星高度角临界值,所模拟的实际卫星应该是整个应用区域内都不可见的实际卫星。
(2)实际卫星—伪卫星—用户所构成的转发线路相对于直线线路的误差要尽可能的小;
(3)伪卫星与现有卫星组成的联合星座的GDOP值小:在满足一定折线误差的情况下,选择能够产生最优或次优GDOP值(几何精度因子)的伪卫星,这样可以提高水平定位精度。方法如下:
根据已经收到的卫星,在不可见卫星中选择相应数量卫星,使得总卫星数达到4颗。根据用户接收机的概略坐标和卫星的概略星历,计算GDOP值,在所有不可见卫星中进行遍历,计算GDOP值最小的一组,伪卫星则模拟该组卫星中的不可见卫星。4个伪卫星分别在90度方位角的不可见区域内进行遍历搜索。
基于上述选星原则,在选星过程中还要注意以下两点:
(1)选星区域:以用户设备接收天线位置为原点,建立测量坐标系,在用户坐标系内判断实际卫星相对用户的高度角(即仰角),同时根据障碍物的位置计算不可见范围的最低高度角,如果实际卫星高度角大于最低高度角,则该实际卫星视为可见卫星。伪卫星在排除可见卫星的不可见区域内选择,实现伪卫星信号模拟。
(2)伪卫星与用户的张角θ:为了减小信号折线传输的误差,需要选择张角θ小的实际卫星,具体实现时可选择与伪卫星相对于用户的高度角、方位角接近的实际卫星。
上述步骤2)、步骤4)中涉及到的伪卫星选星方法,具体如下(如图4所示):收集实际卫星参数,根据用户和障碍物的位置关系区分可见区域,在不可见区域内选择仰角最优的实际卫星,提取相应实际卫星星号,完成选星。
上述实施例中,GDOP值计算方法如下:
(1)针对伪卫星的不同位置,计算所选伪卫星的星历:以用户设备接收天线位置为原点,建立测量坐标系,在该坐标系内,计算实际卫星的高度角和方位角(如图5所示)。
(a)计算实际卫星在测量坐标系的位置
式中,为实际卫星在ECEF坐标系的坐标;为用户在ECEF坐标系的坐标;(L,B)为用户的经度、纬度;RY、RX和RZ分别是绕Y轴、X轴和Z轴的旋转矩阵。其中:
R X ( &epsiv; X ) = 1 0 0 0 cos &epsiv; X - sin &epsiv; X 0 sin &epsiv; X cos &epsiv; X ,
R Y ( &epsiv; Y ) = cos &epsiv; Y 0 - sin &epsiv; Y 0 1 0 sin &epsiv; Y 0 cos &epsiv; Y ,
R Z ( &epsiv; Z ) = cos &epsiv; Z sin &epsiv; Z 0 - sin &epsiv; Z cos &epsiv; Z 0 0 0 1
其中,εX为B,εY为-90°,εZ为(-90°+L),单位均为度。
(b)根据实际卫星的位置计算实际卫星的高度角E、方位角Az:
E = arctg ( y sc / D ) Az = arccos ( x sc / D ) + 0 , z sc &GreaterEqual; 0 , &pi; , z sc < 0
式中,分别是用户坐标系下卫星的x、y和z轴坐标。
(2)根据实际卫星的高度角E和方位角Az计算GDOP值:
根据实际卫星的位置可以计算出H矩阵为:
H = cos E 1 cos Az 1 cos E 1 sin Az 1 sin E 1 cos E 2 cos Az 2 cos E 2 sin Az 2 si E 2 cos E 3 cos Az 3 cos E 3 sin Az 3 sin E 3 cos E 4 cos Az 4 cos E 4 sin Az 4 sin E 4 ,
式中Ei表示第i颗星的高度角,Azi表示第i颗星的方位角,i=1,2,3,4。
然后根据H矩阵计算(HTH)-1矩阵(HT矩阵是H矩阵的转置),其分量的形式表示为:
( H T H ) - 1 = D 11 D 12 D 13 D 14 D 21 D 22 D 23 D 24 D 31 D 32 D 33 D 34 D 41 D 42 D 43 D 44 ;
GDOP值能够作为(HTH)-1矩阵迹的平方根计算出:
GDOP = D 11 + D 12 + D 13 + D 14 .
上述步骤5)中,由于实际卫星的运动速度较快,实际卫星的位置是实时变化的;所以选星也要实时地进行计算,并根据卫星的位置进行所要模拟的卫星的切换。具体的换星过程如下:
(1)伪卫星信号发生器选好符合要求的不可见区域的实际卫星(星号n1)并播发这颗卫星的信号;
(2)计算不可见区域的仰角范围,实时进行伪卫星仰角值计算,根据卫星运动轨迹推算该卫星(星号n1)走出不可见区域的时间;
(3)如果计算伪卫星离开的时间小于10分钟,重新进行选星过程,并选出当前条件下更优的伪卫星(星号n2);
(4)收集新的伪卫星星历(星号n2),同时计算出该卫星(星号n2)的实时电文参数和观测量参数,包括伪距参数、多普勒参数和误差参数;
(5)星号n2的卫星数据准备好之后,通知伪卫星信号发生器进行换星操作,关闭当前卫星信号,并在下一个1PPS发送星号n2的卫星信号;
(6)换星过程完成。
步骤(1)至步骤(6)的换星过程是一个动态过程,每台伪卫星信号发生器都会独立的运行相关的换星过程,保证系统中出现的伪卫星随着时间变化能够顺利切换,同时为服务区域提供连续的伪卫星信号服务。
上述各实施例仅用于说明本发明,各步骤都是可以有所变化的,在本发明技术方案的基础上,凡根据本发明原理对个别步骤进行的改进和等同变换,均不应排除在本发明的保护范围之外。

Claims (6)

1.一种与北斗系统实时同步生成伪卫星信号的方法,其包括以下步骤:
1)由导航接收机接收实际卫星信号,并由导航系统控制器向导航接收机发送多个实际卫星系统星历接收指令,导航接收机获取当前的实际卫星系统星历参数,同时通过1PPS实现将本地信号与实际北斗时的时间同步,并判断是否实现时间同步,若实现时间同步后,将同步标志设置为“时间已同步”,进入步骤2);反之继续通过1PPS实现时间同步;
2)根据当前的可见星和用户位置关系利用选星方法进行伪卫星星号选择,并将选择的星号进行存储;
3)根据存储的星号,将该星号对应的伪卫星信号进行仿真计算,生成与实际星历同步的电文参数和观测量参数;
4)通过已有的信号生成机制将伪卫星信号的伪距参数、多普勒参数和误差参数转换后并输出,信号播发后每隔1秒进行一次判断,判断当前伪卫星是否已经不再满足选星要求,如果不满足要求则重新输出伪卫星信号;满足则根据选星方法选出1颗备用卫星;
5)计算选出的备用卫星的电文参数和观测量参数,并判断当前时刻是否应该进行换星;如果满足换星时刻,则进行换星处理,关闭当前卫星信号,启动新的卫星,换星过程完成;反之,返回步骤4),循环进行伪卫星的生成,实现了与北斗系统的同步。
2.如权利要求1所述的一种与北斗系统实时同步生成伪卫星信号的方法,其特征在于:所述步骤1)中,所述通过1PPS实现将本地信号与实际北斗时的时间同步包括以下步骤:
(1)设置时钟同步结构,该时钟同步结构包括授时型导航接收机、接收天线、带有本地原子钟的时钟驯服模块和伪卫星信号发生器;
(2)授时型导航接收机通过接收天线接收实际卫星信号,并读取实际卫星系统UTC时;授时型导航接收机将实际卫星系统UTC时进行读秒处理后,利用1PPS信号将时间信息传输至时钟驯服模块,同时通过串口将该时间信息传输至伪卫星信号发生器;通过时钟驯服模块同步生成10MHz时钟信号和稳定的1PPS信号;
(3)授时型导航接收机通过时钟驯服模块输出的10MHz时钟信号作为伪卫星信号发生器的时钟基准信号,授时型导航接收机通过时钟驯服模块输出的下一个1PPS信号作为伪卫星信号发生器的外触发信号;
(4)伪卫星信号发生器将读取的实际卫星UTC时间转换后产生伪卫星系统时间,伪卫星信号发生器工作在外触发模式下,读取实际卫星系统UTC时间并通过时钟驯服模块获取到稳定的1PPS信号后,判决是否具备从下一秒开始和实际时间对齐的条件,判决条件成立后开始计算输出下一秒的伪卫星信号,并同步利用1PPS信号和10MHz将伪卫星系统本地时间与UTC时间进行对齐,实现同步。
3.如权利要求1或2所述的一种与北斗系统实时同步生成伪卫星信号的方法,其特征在于:所述步骤4)中,所述伪卫星的选星原则如下:
(1)用户范围内难以正常接收的实际卫星,即不在用户视野内的实际卫星;
(2)实际卫星—伪卫星—用户所构成的转发线路相对于直线线路的误差要尽可能的小;
(3)伪卫星与现有卫星组成的联合星座的GDOP值小。
4.如权利要求3所述的一种与北斗系统实时同步生成伪卫星信号的方法,其特征在于:所述GDOP值计算方法如下:
(1)针对伪卫星的不同位置,计算所选伪卫星的星历:以用户设备接收天线位置为原点,建立测量坐标系,在该坐标系内,计算实际卫星的高度角和方位角;
(a)计算实际卫星在测量坐标系的位置
式中,为实际卫星在ECEF坐标系的坐标;为用户在ECEF坐标系的坐标;(L,B)为用户的经度、纬度;RY、RX和RZ分别是绕Y轴、X轴和Z轴的旋转矩阵;
(b)根据实际卫星的位置计算实际卫星的高度角E、方位角Az:
E = arctg ( y sc / D ) Az = arccos ( x sc / D ) + 0 , z sc &GreaterEqual; 0 &pi; , z sc < 0 ,
式中,xsc、ysc、zsc分别是用户坐标系下卫星的x、y和z轴坐标;
(2)根据实际卫星的高度角E和方位角Az计算GDOP值:
GDOP = D 11 + D 12 + D 13 + D 14 .
5.如权利要求1或2或4所述的一种与北斗系统实时同步生成伪卫星信号的方法,其特征在于:所述步骤2)、步骤4)中,所述伪卫星选星方法为:收集实际卫星参数,根据用户和障碍物的位置关系区分可见区域,在不可见区域内选择仰角最优的实际卫星,提取相应实际卫星星号,完成选星。
6.如权利要求1或2或4所述的一种与北斗系统实时同步生成伪卫星信号的方法,其特征在于:所述步骤5)中,所述换星过程如下:
(1)伪卫星信号发生器选好符合要求的不可见区域的实际卫星并播发这颗卫星的信号;
(2)计算不可见区域的仰角范围,实时进行伪卫星仰角值计算,根据卫星运动轨迹推算该卫星走出不可见区域的时间;
(3)如果计算伪卫星离开的时间小于10分钟,重新进行选星过程,并选出当前条件下更优的伪卫星;
(4)收集新的伪卫星星历,同时计算出该卫星的实时电文参数和观测量参数,包括伪距参数、多普勒参数和误差参数;
(5)新的伪卫星的卫星数据准备好之后,通知伪卫星信号发生器进行换星操作,关闭当前卫星信号,并在下一个1PPS发送新的伪卫星的卫星信号;
(6)换星过程完成。
CN201510145678.3A 2015-03-30 2015-03-30 一种与北斗系统实时同步生成伪卫星信号的方法 Expired - Fee Related CN104749588B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510145678.3A CN104749588B (zh) 2015-03-30 2015-03-30 一种与北斗系统实时同步生成伪卫星信号的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510145678.3A CN104749588B (zh) 2015-03-30 2015-03-30 一种与北斗系统实时同步生成伪卫星信号的方法

Publications (2)

Publication Number Publication Date
CN104749588A true CN104749588A (zh) 2015-07-01
CN104749588B CN104749588B (zh) 2017-02-01

Family

ID=53589562

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510145678.3A Expired - Fee Related CN104749588B (zh) 2015-03-30 2015-03-30 一种与北斗系统实时同步生成伪卫星信号的方法

Country Status (1)

Country Link
CN (1) CN104749588B (zh)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105974455A (zh) * 2016-07-06 2016-09-28 四川泰富地面北斗科技股份有限公司 一种基于光纤授时的陆基定位方法及装置
CN106526616A (zh) * 2016-12-22 2017-03-22 格星微电子科技成都有限公司 一种基于虚拟网格的伪卫星室内导航方法
CN106959452A (zh) * 2017-04-07 2017-07-18 湖南国科防务电子科技有限公司 一种卫星导航信号再生器、系统及方法
CN106980122A (zh) * 2017-04-07 2017-07-25 湖南国科防务电子科技有限公司 室内环境下卫星导航定位增强系统及方法
CN108508454A (zh) * 2017-02-27 2018-09-07 精工爱普生株式会社 定时信号输出装置、电子设备以及移动体
CN108768505A (zh) * 2018-06-19 2018-11-06 湖南斯北图科技有限公司 一种物联网卫星广播数据获取方法和装置
CN108897017A (zh) * 2018-07-05 2018-11-27 格星微电子科技成都有限公司 一种基于cs构架的室内伪卫星星历更新方法
CN109100746A (zh) * 2018-09-30 2018-12-28 武汉大学 一种基于转发节点的隧道定位系统及方法
CN109541650A (zh) * 2018-12-13 2019-03-29 北京无线电计量测试研究所 一种长周期伪随机序列本地同步码生成方法和装置
CN110471090A (zh) * 2019-07-03 2019-11-19 中国人民解放军63686部队 一种基于自解算的船载卫导设备导航电文数据正确性检验系统
CN110708669A (zh) * 2019-10-09 2020-01-17 广州小鹏汽车科技有限公司 一种车辆定位方法、装置、车辆和计算机可读存储介质
JP2021032561A (ja) * 2019-08-13 2021-03-01 東芝エレクトロニックシステムズ株式会社 測位環境提供システム
CN112526550A (zh) * 2020-11-26 2021-03-19 中国电子科技集团公司第五十四研究所 一种同源阵列伪卫星系统
CN112698563A (zh) * 2020-12-02 2021-04-23 四川九洲北斗导航与位置服务有限公司 卫星授时方法、装置、电子设备和存储介质
CN113514857A (zh) * 2021-07-13 2021-10-19 阎镜予 伪卫星电文生成方法、装置、计算机设备及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110037648A1 (en) * 2009-08-17 2011-02-17 Electronics And Telecommunications Research Institute Apparatus for transmitting pseudollite signal based on single clock and positioning system using the same
CN102540214A (zh) * 2012-01-12 2012-07-04 电子科技大学 一种用于导航卫星系统信号源的平滑选星方法
CN103064092A (zh) * 2012-12-28 2013-04-24 中国科学院光电研究院 一种导航卫星的选择方法
CN104122567A (zh) * 2014-07-29 2014-10-29 中国电子科技集团公司第五十四研究所 伪卫星、gps和北斗导航系统组合的定位方法
CN104133221A (zh) * 2014-07-24 2014-11-05 中国人民解放军信息工程大学 一种基于通用接收机的伪卫星定位系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110037648A1 (en) * 2009-08-17 2011-02-17 Electronics And Telecommunications Research Institute Apparatus for transmitting pseudollite signal based on single clock and positioning system using the same
CN102540214A (zh) * 2012-01-12 2012-07-04 电子科技大学 一种用于导航卫星系统信号源的平滑选星方法
CN103064092A (zh) * 2012-12-28 2013-04-24 中国科学院光电研究院 一种导航卫星的选择方法
CN104133221A (zh) * 2014-07-24 2014-11-05 中国人民解放军信息工程大学 一种基于通用接收机的伪卫星定位系统
CN104122567A (zh) * 2014-07-29 2014-10-29 中国电子科技集团公司第五十四研究所 伪卫星、gps和北斗导航系统组合的定位方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
高鹏: ""GPS模拟信号源基带算法设计与软件实现"", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105974455A (zh) * 2016-07-06 2016-09-28 四川泰富地面北斗科技股份有限公司 一种基于光纤授时的陆基定位方法及装置
CN106526616A (zh) * 2016-12-22 2017-03-22 格星微电子科技成都有限公司 一种基于虚拟网格的伪卫星室内导航方法
CN108508454A (zh) * 2017-02-27 2018-09-07 精工爱普生株式会社 定时信号输出装置、电子设备以及移动体
CN106959452A (zh) * 2017-04-07 2017-07-18 湖南国科防务电子科技有限公司 一种卫星导航信号再生器、系统及方法
CN106980122A (zh) * 2017-04-07 2017-07-25 湖南国科防务电子科技有限公司 室内环境下卫星导航定位增强系统及方法
CN108768505B (zh) * 2018-06-19 2020-09-11 湖南斯北图科技有限公司 一种物联网卫星广播数据获取方法和装置
CN108768505A (zh) * 2018-06-19 2018-11-06 湖南斯北图科技有限公司 一种物联网卫星广播数据获取方法和装置
CN108897017A (zh) * 2018-07-05 2018-11-27 格星微电子科技成都有限公司 一种基于cs构架的室内伪卫星星历更新方法
CN109100746A (zh) * 2018-09-30 2018-12-28 武汉大学 一种基于转发节点的隧道定位系统及方法
CN109541650A (zh) * 2018-12-13 2019-03-29 北京无线电计量测试研究所 一种长周期伪随机序列本地同步码生成方法和装置
CN110471090A (zh) * 2019-07-03 2019-11-19 中国人民解放军63686部队 一种基于自解算的船载卫导设备导航电文数据正确性检验系统
JP2021032561A (ja) * 2019-08-13 2021-03-01 東芝エレクトロニックシステムズ株式会社 測位環境提供システム
CN110708669A (zh) * 2019-10-09 2020-01-17 广州小鹏汽车科技有限公司 一种车辆定位方法、装置、车辆和计算机可读存储介质
CN110708669B (zh) * 2019-10-09 2021-02-19 广州小鹏汽车科技有限公司 一种车辆定位方法、装置、车辆和计算机可读存储介质
CN112526550A (zh) * 2020-11-26 2021-03-19 中国电子科技集团公司第五十四研究所 一种同源阵列伪卫星系统
CN112698563A (zh) * 2020-12-02 2021-04-23 四川九洲北斗导航与位置服务有限公司 卫星授时方法、装置、电子设备和存储介质
CN113514857A (zh) * 2021-07-13 2021-10-19 阎镜予 伪卫星电文生成方法、装置、计算机设备及存储介质

Also Published As

Publication number Publication date
CN104749588B (zh) 2017-02-01

Similar Documents

Publication Publication Date Title
CN104749588B (zh) 一种与北斗系统实时同步生成伪卫星信号的方法
CN106569239B (zh) 一种广播式网络rtk定位技术
CN102326092B (zh) 用于组合使用标准rtk系统与全球载波相位差分定位系统的导航接收器和方法
CN105044747B (zh) 一种基于多星共视和滤波的时间同步装置及其方法
CN101038337B (zh) 储存全球导航卫星系统的辅助校正数据的方法、装置及全球导航卫星系统接收器
CN105353391A (zh) 一种支持多类型定位终端的多网融合定位增强系统及方法
CN102288977B (zh) 一种基于北斗gps双模定位模块的差分定位方法
CN111856534B (zh) 智能终端的双模gnss载波精密单点定位方法及系统
CN105158780A (zh) 一种基于多种导航卫星可互换的导航定位方法
CN106405579A (zh) 一种实时卫星导航模拟方法及设备
CN100471122C (zh) 基于p2p技术的导航定位增强信息互联网实时播发方法
Jia et al. Ground vehicle navigation integrity monitoring for multi-constellation GNSS fused with cellular signals of opportunity
CN106383355A (zh) 基于高精度时间的gnss时差监测方法及系统
Li et al. Review of PPP–RTK: Achievements, challenges, and opportunities
CN107807373A (zh) 基于移动智能终端的gnss高精度定位方法
Li et al. Improving availability and accuracy of GPS/BDS positioning using QZSS for single receiver
Wolf Defining GPS and its Capabilities
CN104392108A (zh) 一种采用迭代差分算法的远程定位系统及方法
Liu Positioning performance of single-frequency GNSS receiver using Australian regional ionospheric corrections
CN112731268A (zh) 一种差分数据的处理方法和定位跟踪系统
Rocken et al. A new real-time global GPS and GLONASS precise positioning correction service: Apex
KR102511969B1 (ko) 선박기준국 추가에 의한 sbas 시스템을 위한 정지 궤도 위성의 정밀 궤도 결정 방법
Samper et al. Advantages and drawbacks of the Precise Point Positioning (PPP) technique for earthquake, tsunami prediction and monitoring
Wang et al. Design and implementation of programmable multi-mode GNSS signal simulator
Caibing et al. Research on BeiDou positioning performance in urban environments

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20160330

Address after: 100190, room 816, room 8, floor 9, North Fourth Ring Road, Haidian District, Beijing,

Applicant after: Xi Beidou technology development (Beijing) Co., Ltd.

Address before: 100096, room 2, unit 1, building 87, building No. 1303, Beijing, Xisanqi, Changping District, building materials West Road, China

Applicant before: Beijing Hua Yunzhi joins Science and Technology Ltd.

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: 20170201

Termination date: 20200330

CF01 Termination of patent right due to non-payment of annual fee