CN113514857A - 伪卫星电文生成方法、装置、计算机设备及存储介质 - Google Patents

伪卫星电文生成方法、装置、计算机设备及存储介质 Download PDF

Info

Publication number
CN113514857A
CN113514857A CN202110790889.8A CN202110790889A CN113514857A CN 113514857 A CN113514857 A CN 113514857A CN 202110790889 A CN202110790889 A CN 202110790889A CN 113514857 A CN113514857 A CN 113514857A
Authority
CN
China
Prior art keywords
pseudolite
latitude
coordinate
longitude
target
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
CN202110790889.8A
Other languages
English (en)
Other versions
CN113514857B (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.)
Shanghai Azimuth Data Technology Co ltd
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN202110790889.8A priority Critical patent/CN113514857B/zh
Publication of CN113514857A publication Critical patent/CN113514857A/zh
Application granted granted Critical
Publication of CN113514857B publication Critical patent/CN113514857B/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/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • 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

Abstract

本发明实施例公开了伪卫星电文生成方法、装置、计算机设备及存储介质。该方法包括根据地球半径及预置半长轴转换函数对伪卫星的轨道半长轴进行转换，得到目标半长轴；获取伪卫星的地固系坐标，基于目标半长轴及地固系坐标进行地心经纬度计算得到第一经纬度，对第一经纬度的纬度进行转换得到第一纬度，并根据第一纬度和第一经纬度的经度对地固系坐标进行转换，得到第一坐标表征；根据第一坐标表征及预置轨道速度函数对伪卫星的轨道速度进行计算，得到目标轨道速度；基于第一坐标表征、目标轨道速度及预置北斗三号卫星协议进行伪卫星电文转换，得到目标伪卫星电文。该方法保证了伪卫星电文被接收机识别的有效性。

Description

伪卫星电文生成方法、装置、计算机设备及存储介质
技术领域
本发明涉及通信技术领域，尤其涉及一种伪卫星电文生成方法、装置、计算机设备及存储介质。
背景技术
目前，全球卫星导航系统(GNSS)发展已经比较成熟，可以满足大部分室外场景下的定位、导航、授时需求。但是，GNSS系统卫星轨道高，以提供广域服务为基本目标，到达用户的信号强度低，导致导航信号很容易被遮挡，特别是对于城市峡谷、隧道、地下、室内、水下等区域，卫星导航信号完全被遮挡，无法提供导航服务。
伪卫星技术是解决室内、隧道、水下等多种卫星导航信号遮挡场景下定位问题的一种有效方法。伪卫星系统采用与卫星导航系统相同的技术原理，在上述场景中布设若干导航信号发射设备，发射与真实导航卫星信号体制相同的无线电信号，对无卫星导航信号的区域进行覆盖。终端通过接收伪卫星节点发射的信号，完成距离测量，进而实现高精度定位。
根据三球交汇定位原理，准确高效的表征伪卫星位置是实现定位解算的前提。GNSS系统通过广播1个参考时间参数、6个开普勒轨道根数和11个轨道摄动修正参数准确表征导航卫星轨道，用户终端根据时间信息可准确计算信号发射时刻的卫星位置。但是，伪卫星为固定于地球表面的一个固定参考点，不符合卫星的轨道动力学，其坐标在地固系为一固定参数。因此，GNSS系统伪卫星电文无法直接表征伪卫星坐标。
发明内容
本发明实施例提供了一种伪卫星电文生成方法、装置、计算机设备及存储介质，旨在解决现有技术中伪卫星电文格式无法兼容全球卫星导航系统的问题。
第一方面，本发明实施例提供了一种伪卫星电文生成方法，其包括：
根据地球半径及预置半长轴转换函数对所述伪卫星的轨道半长轴进行转换，得到目标半长轴；
获取伪卫星的地固系坐标，基于所述目标半长轴及地固系坐标进行地心经纬度计算得到第一经纬度，对所述第一经纬度的纬度进行转换得到第一纬度，并根据所述第一纬度和所述第一经纬度的经度对所述地固系坐标进行转换，得到第一坐标表征；
根据所述第一坐标表征及预置轨道速度函数对所述伪卫星的轨道速度进行计算，得到目标轨道速度；
基于所述第一坐标表征、目标轨道速度及预置北斗三号卫星协议进行伪卫星电文转换，得到目标伪卫星电文。
第二方面，本发明实施例提供了一种伪卫星电文生成装置，其包括：
半长轴转换模块，用于根据地球半径及预置半长轴转换函数对所述伪卫星的轨道半长轴进行转换，得到目标半长轴；
坐标转换模块，用于获取伪卫星的地固系坐标，基于所述目标半长轴及地固系坐标进行地心经纬度计算得到第一经纬度，对所述第一经纬度的纬度进行转换得到第一纬度，并根据所述第一纬度和所述第一经纬度的经度对所述地固系坐标进行转换，得到第一坐标表征；
速度计算模块，用于根据所述第一坐标表征及预置轨道速度函数对所述伪卫星的轨道速度进行计算，得到目标轨道速度；
基于所述第一坐标表征、目标轨道速度及预置北斗三号卫星协议进行伪卫星电文转换，得到目标伪卫星电文。
第三方面，本发明实施例又提供了一种计算机设备，其包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述第一方面所述的伪卫星电文生成方法。
第四方面，本发明实施例还提供了一种计算机可读存储介质，其中所述计算机可读存储介质存储有计算机程序，所述计算机程序当被处理器执行时使所述处理器执行上述第一方面所述的伪卫星电文生成方法。
本发明实施例提供了一种伪卫星电文生成方法、装置、计算机设备及存储介质。该方法包括根据地球半径及预置半长轴转换函数对伪卫星的轨道半长轴进行转换，得到目标半长轴；获取伪卫星的地固系坐标，基于目标半长轴及地固系坐标进行地心经纬度计算得到第一经纬度，对第一经纬度的纬度进行转换得到第一纬度，并根据第一纬度和第一经纬度的经度对地固系坐标进行转换，得到第一坐标表征；根据第一坐标表征及预置轨道速度函数对伪卫星的轨道速度进行计算，得到目标轨道速度；基于第一坐标表征、目标轨道速度及预置北斗三号卫星协议进行伪卫星电文转换，得到目标伪卫星电文。该方法通过将伪卫星地固系坐标变换为导航卫星轨道，并根据其轨道生成伪卫星电文，使得到的电文与伪卫星电文的格式完全相同，保证了伪卫星电文的有效性。
附图说明
为了更清楚地说明本发明实施例技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的伪卫星电文生成方法的流程示意图；
图2为本发明实施例提供的伪卫星电文解析方法的流程示意图；
图3为本发明实施例提供的伪卫星电文生成装置的示意性框图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
应当理解，当在本说明书和所附权利要求书中使用时，术语“包括”和“包含”指示所描述特征、整体、步骤、操作、元素和/或组件的存在，但并不排除一个或多个其它特征、整体、步骤、操作、元素、组件和/或其集合的存在或添加。
还应当理解，在此本发明说明书中所使用的术语仅仅是出于描述特定实施例的目的而并不意在限制本发明。如在本发明说明书和所附权利要求书中所使用的那样，除非上下文清楚地指明其它情况，否则单数形式的“一”、“一个”及“该”意在包括复数形式。
还应当进一步理解，在本发明说明书和所附权利要求书中使用的术语“和/或”是指相关联列出的项中的一个或多个的任何组合以及所有可能组合，并且包括这些组合。
请参阅图1所示，图1为本发明实施例提供的伪卫星电文生成方法的流程示意图，该方法包括步骤S110～S140。
步骤S110、根据地球半径及预置半长轴转换函数对所述伪卫星的轨道半长轴进行转换，得到目标半长轴；
本实施例中，考虑到伪卫星位于地球表面，向径长度
Figure BDA0003160783530000041
为地球半径，地球半径约为6378km，远小于MEO卫星轨道半长轴。因此，首先根据地球半径及预置半长轴转换函数对伪卫星的轨道半长轴进行转换：
预置半长轴转换函数具体如下：
rps-c＝rps+Δr，
其中，Δr表示修正量，用于保证目标半长轴rps-c为MEO卫星轨道半长轴，rps-c表示目标半长轴，rps表示地球半径。若如与北斗三号卫星导航系统兼容，Δr可取为21528km。
步骤S120、获取伪卫星的地固系坐标，基于所述目标半长轴及地固系坐标进行地心经纬度计算得到第一经纬度，对所述第一经纬度的纬度进行转换得到第一纬度，并根据所述第一纬度和所述第一经纬度的经度对所述地固系坐标进行转换，得到第一坐标表征；
本实施例中，考虑到伪卫星可在全球任意地区部署，其纬度在某些地区大于卫星导航系统卫星轨道倾角标称值inormal时，无法构建轨道倾角为inormal的卫星轨道通过伪卫星的地固系坐标。因此，在半长轴转换的基础上需要对伪卫星纬度进行转换。考虑到卫星轨道倾角定义为轨道面与赤道之间的夹角，因此直接对伪卫星的地心纬度进行修正。伪卫星的地心经纬度可按照如下公式计算：
Figure BDA0003160783530000042
其中，(x,y,z)表示伪卫星的地固系坐标，λps，
Figure BDA0003160783530000043
分别表示伪卫星的第一经纬度的经度和纬度；
按以下公式对纬度进行转换：
Figure BDA0003160783530000044
其中，inormal表示卫星导航系统MEO卫星轨道标称轨道倾角，对于北斗三号卫星导航系统inormal＝55°；Δi表示轨道倾角保护余量，Δi＝5°；
Figure BDA0003160783530000051
表示第一纬度；
所述第一坐标表征如下：
Figure BDA0003160783530000052
步骤S130、根据所述第一坐标表征及预置轨道速度函数对所述伪卫星的轨道速度进行计算，得到目标轨道速度；
本实施例中，只有第一坐标表征无法确定一个轨道面，因此需要计算第一坐标表征对应位置处伪卫星的速度矢量，构建一个过第一坐标表征且轨道倾角为inormal、轨道偏心率为0的轨道，计算过程如下：
取伪卫星当前位置矢量的单位矢量：
Figure BDA0003160783530000053
构建轨道速度如下：
Figure BDA0003160783530000054
其中，为保证偏心率为0，取速度为
Figure BDA0003160783530000055
式中μ为地心引力常数，μ＝3.986004418×1014m3/s2；为保证轨道倾角inormal、轨道偏心率为0，且过第一坐标表征对应的位置时从南向北飞，取速度矢量的单位矢量为：
Figure BDA0003160783530000056
步骤S140、基于所述第一坐标表征、目标轨道速度及预置北斗三号卫星协议进行伪卫星电文转换，得到目标伪卫星电文。
本实施例中，基于所述第一坐标表征、目标轨道速度及预置北斗三号卫星协议进行伪卫星电文转换，得到目标导航电。其中，预置北斗三号卫星协议为北斗三号卫星导航系统伪卫星电文参数定义及北斗三号卫星导航系统伪卫星电文参数格式说明。伪卫星电文由伪卫星的轨道根数及轨道摄动参数组成。具体如下：
计算参考时刻第二赤道坐标系下的轨道根数，首先计算动力学常量：
Figure BDA0003160783530000057
Figure BDA0003160783530000058
式中，
Figure BDA0003160783530000059
表示动量矩常量，E表示卫星机械能常量；
计算轨道半长轴：
Figure BDA0003160783530000061
半通径和动量矩大小满足以下关系：
Figure BDA0003160783530000062
计算偏心率：
Figure BDA0003160783530000063
轨道倾角为动量矩
Figure BDA0003160783530000064
与坐标系Z轴的夹角：
Figure BDA0003160783530000065
升交点经度为轨道节线矢量与坐标系X轴的夹角，则：
Ω0＝arctan2(Ny,Nx)，
其中，轨道单位节线矢量为
Figure BDA0003160783530000066
计算伪卫星的平均角速度：
Figure BDA0003160783530000067
计算偏近点角：
Figure BDA0003160783530000068
由偏近点角Ek计算当前时刻平近点角：
M0＝E0-esinE0；
计算纬度幅角值：
Figure BDA0003160783530000069
式中，
Figure BDA00031607835300000610
为
Figure BDA00031607835300000611
的三维分量；
计算真近点角：
Figure BDA00031607835300000612
计算近地点幅角：
ω＝u-f；
根据预置北斗三号卫星协议对第二赤道坐标系中的轨道根数进行转换，得到伪卫星电文中的轨道根数为：
Figure BDA00031607835300000613
取伪卫星电文参考时刻为当前时刻tcur对应的整小时时刻，将伪卫星电文中的轨道摄动参数全部取零，得到：
Figure BDA0003160783530000071
将式(1)与式(2)组合，并添加伪卫星轨道类型信息，得到目标伪卫星电文。
实施例一，将位于北京附近的某个伪卫星坐标表征为与GNSS系统导航电文兼容的伪卫星电文，该伪卫星在地固系坐标为[-2158577.461,4381498.809,4087970.331]，按照本申请中的伪卫星电文生成方法，每个步骤的计算结果如下：
(1)伪卫星半长轴修正
修正之前伪卫星坐标向径长度为：rps＝6369339.832m；
取修正量，修正之后卫星半长轴为：rps-c＝27897339.832m。
(2)伪卫星纬度修正
修正之前伪卫星的经纬度为：[116.227500°，39.927674°]；
修正之后伪卫星的经纬度为：[116.227500°，22.182041°]；
修正之后坐标为
Figure BDA0003160783530000072
(3)轨道面确定
当前时刻为保证过
Figure BDA0003160783530000073
点的轨道倾角为55°、偏心率为零，计算轨道当前时刻速度为：
Figure BDA0003160783530000074
(4)导航电文生成
进而按照轨道动力学，将
Figure BDA0003160783530000075
转化为轨道根数：
[a,e,i0,Ω0,ω,M0]＝[27897339.832，0,0.959931,1.739031,0.667992,6.094207]；
(5)格式调整
将轨道根数按照三号卫星协议进行范围变换、单位统一，即可在伪卫星上播放给用户，提供GNSS导航信号缺失场景下的导航服务。
请参阅图2所示，图2为本发明实施例提供的伪卫星电文解析方法的流程示意图，该方法包括步骤H110～H140。
步骤H110、根据预置提取规则从接收的电文中提取伪卫星电文，并根据预置北斗三号卫星协议对所述伪卫星电文进行轨道根数解析，得到第二坐标表征；
本实施例中，电文接收端按照全球导航卫星系统(GNSS)的导航电文物理帧格式解析导航电文之后，先判断接收的导航电文是否为伪卫星电文。由于常规导航电文的轨道摄动参数不可能同时为零，因此根据轨道摄动参数判断导航电文是否为伪卫星电文。若接收的导航电文为伪卫星电文，则根据北斗三号卫星协议对伪卫星电文进行轨道根数解析，得到第二坐标表征，具体如下：
计算纬度幅角值：
uk＝M0+ω；
计算轨道面半径：
rk＝a；
轨道面中伪卫星的x坐标：
xp＝rkcosuk；
轨道面中伪卫星的y坐标：
yp＝rksinuk；
地固坐标系中伪卫星的x坐标：
xps-c＝xpcosΩ0-ypcosi0sinΩ0；
地固坐标系中伪卫星的y坐标：
yps-c＝xpsinΩ0-yp|cosi0cosΩ0；
地固坐标系中伪卫星的z坐标：
zps-c＝ypsini0；
则伪卫星的第二坐标表征为：
Figure BDA0003160783530000081
步骤H120、根据所述第二坐标表征及预置经纬度转换函数进行经纬度转换得到第二经纬度，并基于所述第二经纬度进行纬度修正，得到第二纬度；
步骤H130、根据所述第二经纬度的经度和第二纬度进行坐标表征转换，得到所述第二坐标表征的目标单位矢量；
本实施例中，基于第二坐标表征及预置经纬度转换函数进行经纬度转换得到第二经纬度。其中，基于第二坐标表征，按如下公式进行经纬度转换：
Figure BDA0003160783530000082
式中，
Figure BDA0003160783530000083
(x,y,z)表示伪卫星的地固系坐标，λps-c，
Figure BDA0003160783530000084
分别表示伪卫星的第二经纬度的经度和纬度。
接着基于第二经纬度进行纬度修正，得到第二纬度：
Figure BDA0003160783530000091
式中，inormal表示卫星导航系统MEO卫星轨道标称轨道倾角，对于北斗三号卫星导航系统inormal＝55°；Δi表示轨道倾角保护余量，Δi＝5°，
Figure BDA0003160783530000092
表示第一纬度；
转换纬度后得到第二坐标表征的单位矢量为：
Figure BDA0003160783530000093
步骤H140、对所述伪卫星的轨道半长轴进行还原，基于还原的轨道半长轴、所述目标单位矢量和第二坐标表征进行坐标转换，得到所述伪卫星电文对应的地固系坐标。
本实施例中，为了计算伪卫星电文对应的地固系坐标，首先对伪卫星的轨道半长轴进行还原，基于还原的轨道半长轴、目标单位矢量和第二坐标表征进行坐标变换，得到伪卫星电文对应的地固系坐标。其中，按以下技术函数将伪卫星的轨道半长轴进行还原：
rps＝rps-c-Δr，
Figure BDA0003160783530000094
式中，Δr为半长轴修正量。
实施例二，将实施例一种的伪卫星电文图2的伪卫星电文解析方法进行解析，得到伪卫星的地固系坐标。具体如下：
(1)卫星地固系坐标计算
将导航电文进行固定转换生成卫星坐标，转换之后的坐标为：
Figure BDA0003160783530000095
(2)伪卫星纬度修正
修正之前伪卫星的经纬度为：[116.227500°，22.182041°]；
修正之后伪卫星的经纬度为：[116.227500°，39.927674°]；
修正之后的卫星向径单位矢量为
Figure BDA0003160783530000096
(3)伪卫星半长轴修正
修正之前卫星半长轴为：rps-c＝27897339.832m；
修正之后伪卫星坐标向径长度为：rps＝6369339.832m；
(4)伪卫星坐标解析
融合伪卫星单位矢量
Figure BDA0003160783530000101
和修正之后伪卫星坐标向径长度rps，可计算伪卫星位置为：
Figure BDA0003160783530000102
与待表征的坐标参数几乎相等。
该方法通过将伪卫星地固系坐标变换为导航卫星轨道，并根据其轨道生成伪卫星电文，使得到的电文与伪卫星电文的格式完全相同，保证了伪卫星电文的有效性。
本发明实施例还提供一种伪卫星电文生成装置，该伪卫星电文生成装置用于执行前述伪卫星电文生成方法的任一实施例。具体地，请参阅图3，图3是本发明实施例提供的伪卫星电文生成装置的示意性框图。该伪卫星电文生成装置100可以配置于服务器中。
如图3所示，伪卫星电文生成装置100包括半长轴转换模块110、坐标转换模块120、速度计算模块130、电文转换模块140。
半长轴转换模块110，用于根据地球半径及预置半长轴转换函数对所述伪卫星的轨道半长轴进行转换，得到目标半长轴；
坐标转换模块120，用于获取伪卫星的地固系坐标，基于所述目标半长轴及地固系坐标进行地心经纬度计算得到第一经纬度，对所述第一经纬度的纬度进行转换得到第一纬度130，并根据所述第一纬度和所述第一经纬度的经度对所述地固系坐标进行转换，得到第一坐标表征；
速度计算模块140，用于根据所述第一坐标表征及预置轨道速度函数对所述伪卫星的轨道速度进行计算，得到目标轨道速度；
电文转换模块，用于基于所述第一坐标表征、目标轨道速度及预置北斗三号卫星协议进行伪卫星电文转换，得到目标伪卫星电文。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，其中，所述处理器执行所述计算机程序时实现如上所述的伪卫星电文生成方法。
在本发明的另一实施例中提供计算机可读存储介质。该计算机可读存储介质可以为非易失性的计算机可读存储介质。该计算机可读存储介质存储有计算机程序，所述计算机程序当被处理器执行时使所述处理器执行如上所述的基于伪卫星电文生成方法。
所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，上述描述的设备、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
在本发明所提供的几个实施例中，应该理解到，所揭露的设备、装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为逻辑功能划分，实际实现时可以有另外的划分方式，也可以将具有相同功能的单元集合成一个单元，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另外，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口、装置或单元的间接耦合或通信连接，也可以是电的，机械的或其它的形式连接。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本发明实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以是两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分，或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到各种等效的修改或替换，这些修改或替换都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。

Claims (10)

1.一种伪卫星电文生成方法，其特征在于，包括：
根据地球半径及预置半长轴转换函数对所述伪卫星的轨道半长轴进行转换，得到目标半长轴；
获取伪卫星的地固系坐标，基于所述目标半长轴及地固系坐标进行地心经纬度计算得到第一经纬度，对所述第一经纬度的纬度进行转换得到第一纬度，并根据所述第一纬度和所述第一经纬度的经度对所述地固系坐标进行转换，得到第一坐标表征；
根据所述第一坐标表征及预置轨道速度函数对所述伪卫星的轨道速度进行计算，得到目标轨道速度；
基于所述第一坐标表征、目标轨道速度及预置北斗三号卫星协议进行伪卫星电文转换，得到目标伪卫星电文。
2.根据权利要求1所述的伪卫星电文生成方法，其特征在于，预置半长轴转换函数具体如下：
rps-c＝rps+Δr，
其中，Δr表示修正量，rps-c表示目标半长轴，rps表示地球半径。
3.根据权利要求2所述的伪卫星电文生成方法，其特征在于，所述获取伪卫星的地固系坐标，基于所述目标半长轴及地固系坐标进行地心经纬度计算得到第一经纬度，对所述第一经纬度的纬度进行转换得到第一纬度，并根据所述第一纬度和所述第一经纬度的经度对所述地固系坐标进行转换，得到第一坐标表征，包括：
按以下公式计算所述伪卫星的地心经纬度：
Figure FDA0003160783520000011
其中，(x,y,z)表示伪卫星的地固系坐标，λps，
Figure FDA0003160783520000012
分别表示伪卫星的第一经纬度的经度和纬度；
按以下公式对纬度进行转换：
Figure FDA0003160783520000013
其中，inormal表示轨道倾角，Δi表示轨道倾角保护余量，
Figure FDA0003160783520000014
表示第一纬度；
所述第一坐标表征如下：
Figure FDA0003160783520000021
4.根据权利要求3所述的伪卫星电文生成方法，其特征在于，所述根据所述第一坐标表征及预置轨道速度函数对所述伪卫星的轨道速度进行计算，得到目标轨道速度，包括：
取伪卫星当前坐标矢量的单位矢量：
Figure FDA0003160783520000022
构建轨道速度如下：
Figure FDA0003160783520000023
其中，速度为
Figure FDA0003160783520000024
μ为地心引力常数，取速度矢量的单位矢量为：
Figure FDA0003160783520000025
5.根据权利要求4所述的伪卫星电文生成方法，其特征在于，所述基于所述第一坐标表征、目标轨道速度及预置北斗三号卫星协议进行伪卫星电文转换，得到目标伪卫星电文，包括：
计算参考时刻第二赤道坐标系下的轨道根数，首先计算动力学常量：
Figure FDA0003160783520000026
Figure FDA0003160783520000027
式中，
Figure FDA0003160783520000028
表示动量矩常量，E表示卫星机械能常量；
计算轨道半长轴：
Figure FDA0003160783520000029
半通径和动量矩大小满足以下关系：
Figure FDA00031607835200000210
计算偏心率：
Figure FDA00031607835200000211
轨道倾角为动量矩
Figure FDA00031607835200000212
与坐标系Z轴的夹角：
Figure FDA00031607835200000213
升交点经度为轨道节线矢量与坐标系X轴的夹角，则：
Ω0＝arctan2(Ny,Nx)，
其中，轨道单位节线矢量为
Figure FDA0003160783520000031
计算伪卫星的平均角速度：
Figure FDA0003160783520000032
计算偏近点角：
Figure FDA0003160783520000033
由偏近点角Ek计算当前时刻平近点角：
M0＝E0-esinE0；
计算纬度幅角值：
Figure FDA0003160783520000034
式中，
Figure FDA0003160783520000035
为
Figure FDA0003160783520000036
的三维分量；
计算真近点角：
Figure FDA0003160783520000037
计算近地点幅角：
ω＝u-f；
根据预置北斗三号卫星协议对第二赤道坐标系中的轨道根数进行转换，得到伪卫星电文中的轨道根数为：
Figure FDA0003160783520000038
取伪卫星电文参考时刻为当前时刻tcur对应的整小时时刻，将伪卫星电文中的轨道摄动参数全部取零，得到：
Figure FDA0003160783520000039
将式(1)与式(2)组合，并添加伪卫星轨道类型信息，得到目标伪卫星电文。
6.一种伪卫星伪卫星电文解析方法，其特征在于，包括：
根据预置提取规则从接收的电文中提取伪卫星电文，并根据预置北斗三号卫星协议对所述伪卫星电文进行轨道根数解析，得到第二坐标表征；
根据所述第二坐标表征及预置经纬度转换函数进行经纬度转换得到第二经纬度，并基于所述第二经纬度进行纬度修正，得到第二纬度；
根据所述第二经纬度的经度和第二纬度进行坐标表征转换，得到所述第二坐标表征的目标单位矢量；
对所述伪卫星的轨道半长轴进行还原，基于还原的轨道半长轴、所述目标单位矢量和第二坐标表征进行坐标转换，得到所述伪卫星电文对应的地固系坐标。
7.根据权利要求6所述的伪卫星电文生成方法，其特征在于，所述根据预置提取规则从接收的电文中提取伪卫星电文，并根据预置北斗三号卫星协议对所述伪卫星电文进行轨道根数解析，得到第二坐标表征，包括：
计算纬度幅角值：
uk＝M0+ω；
计算轨道面半径：
rk＝a；
轨道面中伪卫星的x坐标：
xp＝rkcosuk；
轨道面中伪卫星的y坐标：
yp＝rksinuk；
地固坐标系中伪卫星的x坐标：
xps-c＝xp cosΩ0-yp cosi0sinΩ0；
地固坐标系中伪卫星的y坐标：
yps-c＝xpsinΩ0-yp cosi0cosΩ0；
地固坐标系中伪卫星的z坐标：
zps-c＝ypsini0；
则伪卫星的第二坐标表征为：
Figure FDA0003160783520000041
8.一种伪卫星电文生成装置，其特征在于，包括：
半长轴转换模块，用于根据地球半径及预置半长轴转换函数对所述伪卫星的轨道半长轴进行转换，得到目标半长轴；
坐标转换模块，用于获取伪卫星的地固系坐标，基于所述目标半长轴及地固系坐标进行地心经纬度计算得到第一经纬度，对所述第一经纬度的纬度进行转换得到第一纬度，并根据所述第一纬度和所述第一经纬度的经度对所述地固系坐标进行转换，得到第一坐标表征；
速度计算模块，用于根据所述第一坐标表征及预置轨道速度函数对所述伪卫星的轨道速度进行计算，得到目标轨道速度；
电文转换模块，用于基于所述第一坐标表征、目标轨道速度及预置北斗三号卫星协议进行伪卫星电文转换，得到目标伪卫星电文。
9.一种计算机设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1至5中任一项所述的伪卫星电文生成方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机程序，所述计算机程序当被处理器执行时使所述处理器执行如权利要求1至5任一项所述的伪卫星电文生成方法。
CN202110790889.8A 2021-07-13 2021-07-13 伪卫星电文生成方法、装置、计算机设备及存储介质 Active CN113514857B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110790889.8A CN113514857B (zh) 2021-07-13 2021-07-13 伪卫星电文生成方法、装置、计算机设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110790889.8A CN113514857B (zh) 2021-07-13 2021-07-13 伪卫星电文生成方法、装置、计算机设备及存储介质

Publications (2)

Publication Number Publication Date
CN113514857A true CN113514857A (zh) 2021-10-19
CN113514857B CN113514857B (zh) 2024-04-26

Family

ID=78066824

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110790889.8A Active CN113514857B (zh) 2021-07-13 2021-07-13 伪卫星电文生成方法、装置、计算机设备及存储介质

Country Status (1)

Country Link
CN (1) CN113514857B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130002484A1 (en) * 2011-07-03 2013-01-03 Daniel A. Katz Indoor navigation with gnss receivers
CN104749588A (zh) * 2015-03-30 2015-07-01 北京华云智联科技有限公司 一种与北斗系统实时同步生成伪卫星信号的方法
CN105093247A (zh) * 2015-07-09 2015-11-25 交通信息通信技术研究发展中心 一种基于北斗的地基导航信号网络系统
CN106291643A (zh) * 2016-07-21 2017-01-04 芜湖航飞科技股份有限公司 一种基于北斗导航定位的伪卫星系统
CN106772487A (zh) * 2016-12-22 2017-05-31 格星微电子科技成都有限公司 一种基于伪点的伪卫星室内导航方法
CN112698361A (zh) * 2021-03-24 2021-04-23 航天宏图信息技术股份有限公司 一种基于伪卫星的定位方法及装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130002484A1 (en) * 2011-07-03 2013-01-03 Daniel A. Katz Indoor navigation with gnss receivers
CN104749588A (zh) * 2015-03-30 2015-07-01 北京华云智联科技有限公司 一种与北斗系统实时同步生成伪卫星信号的方法
CN105093247A (zh) * 2015-07-09 2015-11-25 交通信息通信技术研究发展中心 一种基于北斗的地基导航信号网络系统
CN106291643A (zh) * 2016-07-21 2017-01-04 芜湖航飞科技股份有限公司 一种基于北斗导航定位的伪卫星系统
CN106772487A (zh) * 2016-12-22 2017-05-31 格星微电子科技成都有限公司 一种基于伪点的伪卫星室内导航方法
CN112698361A (zh) * 2021-03-24 2021-04-23 航天宏图信息技术股份有限公司 一种基于伪卫星的定位方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王莹;等;: "全球定位系统信号模拟中导航电文生成与测试", 太赫兹科学与电子信息学报, no. 03 *

Also Published As

Publication number Publication date
CN113514857B (zh) 2024-04-26

Similar Documents

Publication Publication Date Title
Morales et al. Tightly coupled inertial navigation system with signals of opportunity aiding
US8624779B2 (en) Global navigation satellite system (GNSS) reference station integrity monitoring and assurance
US10598792B2 (en) Information processing device and positioning device
CN100437143C (zh) 区域卫星导航系统及方法
JP2007524089A (ja) 高度なナビゲーション性能のための方法およびシステム
Mahmoud et al. VANETs positioning in urban environments: A novel cooperative approach
JP4723932B2 (ja) 測位システム
CN108344415A (zh) 一种组合导航信息融合方法
KR101419339B1 (ko) 위성항법 보정시스템의 보정정보 처리 장치 및 그 방법
CN111123295B (zh) 基于ssr的定位方法及装置、定位系统
Tabti et al. Improvement of single-frequency GPS positioning performance based on EGNOS corrections in Algeria
Shen Improving ambiguity convergence in carrier phase-based precise point positioning
JP2006090912A (ja) 測位装置，情報配信装置，測位方法及び情報配信方法
RietDorf et al. Precise positioning in real-time using navigation satellites and telecommunication
CN113514857A (zh) 伪卫星电文生成方法、装置、计算机设备及存储介质
KR101344426B1 (ko) 위성항법 보정시스템의 보정정보 처리 장치 및 그 방법
JP2010145178A (ja) 移動体位置特定装置
JP3042456B2 (ja) 測位衛星を用いた測位システム
Kitamura et al. Development of Prototype Dual-frequency Multi-constellation SBAS for MSAS
CN117782080B (zh) 基于PPP-B2b/INS的实时空基导航系统及方法
Layh et al. Gps-denied navigator for small uavs
WO2023216766A1 (zh) 一种定位方法以及相关装置
Nabhan et al. Client-server based LBS architecture: A novel positioning module for improved positioning performance
KR20170115786A (ko) 선박기준국 추가에 의한 sbas 시스템을 위한 정지 궤도 위성의 정밀 궤도 결정 방법
Crews et al. Long-term future of GPS

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
TA01 Transfer of patent application right

Effective date of registration: 20220726

Address after: 200000 Building 1, Lane 1333, Jiangnan Avenue, Changxing Town, Chongming District, Shanghai (Shanghai Changxing marine equipment industrial base)

Applicant after: Shanghai azimuth data Technology Co.,Ltd.

Address before: 518000 1103, building B, langtinghaoyuan, No. 7, Guihua Road, Futian District, Shenzhen, Guangdong Province

Applicant before: Yan Jingyu

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant