CN105044745B - 一种圆轨道低轨卫星过顶剩余可见时长预测方法 - Google Patents
一种圆轨道低轨卫星过顶剩余可见时长预测方法 Download PDFInfo
- Publication number
- CN105044745B CN105044745B CN201510417367.8A CN201510417367A CN105044745B CN 105044745 B CN105044745 B CN 105044745B CN 201510417367 A CN201510417367 A CN 201510417367A CN 105044745 B CN105044745 B CN 105044745B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- satellite
- user
- mtd
- 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
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)
- Navigation (AREA)
Abstract
本发明公开了一种圆轨道低轨卫星过顶剩余可见时长预测方法。该方法能够充分利用卫星运行工作特点，通过不断计算和更新卫星空间位置以及对应的时刻信息，并以简易的导航信号的形式向地面用户广播这些信息，地面用户接收该导航信号后，再利用本身的GPS接收设备校准时钟、确定用户位置，实时计算卫星相对用户的仰角和迹角，进而求解出卫星过顶剩余可见时长，具有实现方法简单有效，对导航信号的内容需求少，对卫星星历依赖度低，在确保预测时长的精准性和高时效性前提下，又不增加用户硬件成本。
Description
技术领域
本发明涉及卫星通信领域，特别是涉及一种圆轨道低轨卫星过顶剩余可见时长预测方法。
背景技术
圆轨道低轨卫星是指运行轨道是以地心为圆心的圆形轨道，并且轨道高度在500千米到2000千米之间的卫星。这种圆轨道低轨卫星具有传播时延小、对终端设备天线尺寸和功率消耗要求低等优势，并且卫星发射成本较低，但卫星高度越低，其覆盖范围就越小，通过地面卫星终端用户上空的时间，即过顶时间就越短，因此往往需要多颗低轨卫星组网来满足用户通信需求。
当圆轨道低轨卫星用于通信时，主要存在两方面的不足：一是过顶时间有限，而用户的通信时长有可能大于单颗卫星的过顶时间；二是用户在切换时刻可以在多颗卫星之间进行选择。这两方面的不足导致用户通信过程中需要在多颗低轨卫星之间进行切换，以保证通信的连续性。因此，地面终端用户切换卫星的方法非常重要，如果用户随机选择卫星提供通信服务，有可能选择到当前具有最短可见时长的卫星或者较短可见时长的卫星，进而增加卫星切换的次数，增大了用户的掉话率，而中断当前通话相比阻塞新的通话会给用户带来更加糟糕的用户体验，从降低用户掉话率的角度来看，理想的卫星选择方式是用户在切换时刻能够从所有可见卫星中选择提供最长服务时长的卫星。因此，从用户观测的角度，有必要对卫星相对于用户的剩余可见时长做出精确预测，为用户在切换时刻选择卫星提供技术支撑。
现有技术中，预测过顶卫星剩余可见时长的方法主要是通过下载卫星的星历信息，计算卫星相对于用户的剩余可见时长，但该方法需要用户不断更新卫星的星历信息。另外，Irfan Ali等人在《Predicting the Visibility of LEO Satellites》(IEEETransactions on Aerospace and Electronic Systems,vol.35,no.4,1999)中提出了一种预测过顶卫星剩余可见时长的方法，但该方法需要输入卫星轨道半径、轨道倾角及卫星从南半球到北半球运行过程中同赤道平面的交点所处的经度及对应时刻，该方法计算复杂性较高，实时性不强，难以满足对卫星快速选择的需求。
发明内容
本发明主要解决的技术问题是提供一种圆轨道低轨卫星过顶剩余可见时长预测方法，解决现有技术中过度依赖卫星星历、计算复杂度高、实时性不强等问题。
为解决上述技术问题，本发明采用的一个技术方案是：提供一种圆轨道低轨卫星过顶剩余可见时长预测方法，包括如下步骤：第一步，用户接收卫星导航信号，包括：(1)计算卫星空间位置，由地面控制中心通过轨道预报方法，定期计算该卫星未来一段时间内的空间位置信息及其对应的时刻信息，并在该卫星经过该地面控制中心上空时，更新该卫星的空间位置信息，该卫星的空间位置信息包括经纬度、高程及其对应的时刻信息；(2)卫星校准时钟，该卫星利用自身GPS装置进行时钟校准，并根据校准的时钟信息以及存储的该卫星的空间位置信息，实时调整确定该卫星当前覆盖范围内的星下点经纬度信息；(3)卫星发出导航信号，该卫星定期发出导航信号，该导航信号包括该卫星在其所覆盖范围内的星下点经纬度信息及其对应的时刻信息；(4)用户接收导航信号，判断该卫星是接近用户还是远离用户，该用户接收该卫星发送的连续两次该导航信号，根据该导航信号中的星下点经纬度信息相对于用户位置的远近变化判断该卫星接近还是远离该用户；第二步，计算卫星的仰角；第三步，计算卫星的迹角；第四步，计算卫星相对用户的过顶剩余可见时长。
在本发明圆轨道低轨卫星过顶剩余可见时长预测方法另一实施例中，该轨道预报方法为SGP4方法；该定期计算该卫星未来一段时间内的空间位置信息及其对应的时刻信息的频度是每间隔1秒钟计算1次；该卫星定期发出导航信号的频度是每间隔1秒钟发送1次。
在本发明圆轨道低轨卫星过顶剩余可见时长预测方法另一实施例中，该计算卫星的迹角的方法包括：(1)确定该用户位置信息，得到该用户自身的经度λT和纬度ηT信息；(2)时钟校准，根据所接收导航信号中的时刻信息，将用户时钟同步校准与卫星时钟一致；(3)实时计算迹角，该用户通过接收的该导航信号，计算该用户相对该卫星的迹角γm，计算方法如下：
γm＝min{2sin-1(Ω)}
其中，和是该导航信号中包含的该卫星覆盖范围内的星下点的经度和纬度向量，Ω也为向量。
在本发明圆轨道低轨卫星过顶剩余可见时长预测方法另一实施例中，通过该用户内置的GPS设备测量得到该用户自身的经度λT、纬度ηT和高程κT信息。
在本发明圆轨道低轨卫星过顶剩余可见时长预测方法另一实施例中，该计算卫星的仰角的方法包括：(1)地心地固坐标系转换，将当前该卫星和用户的经纬度和高程坐标统一转换为地心地固坐标系中卫星地固坐标和用户地固坐标；(2)本地坐标系转换，将该卫星地固坐标和该用户地固坐标转换为本地坐标系中的卫星本地坐标和用户原点坐标；(3)实时计算仰角，利用该卫星本地坐标和该用户原点坐标，在该本地坐标系中计算该用户相对该卫星的仰角θ。
在本发明圆轨道低轨卫星过顶剩余可见时长预测方法另一实施例中，该卫星的空间位置信息为(λS,ηS,κS)，其中λS表示该卫星的经度，ηS表示该卫星的纬度，κS表示该卫星的高程，该用户的空间位置信息为(λT,ηT,κT)，其中λT表示该用户的经度，ηT表示该用户的纬度，κT表示该用户的高程，该卫星地固坐标为(XS,YS,ZS)，该用户地固坐标为(XU,YU,ZU)，该地心地固坐标系转换方法是：地球是一个椭球体，长轴a＝6.3781370e+03km，半轴b＝6.356752314245179e+03km，偏心率定义中间值 则该卫星的地固坐标(XS,YS,ZS)和该用户的地固坐标(XU,YU,ZU)分别为：
该本地坐标系转换方法是：该卫星在该以用户为原点坐标的该本地坐标系中的卫星本地坐标(xS,yS,zS)是：
在本发明圆轨道低轨卫星过顶剩余可见时长预测方法另一实施例中，该卫星的空间位置是：经度25.257000°、纬度-84.298000°、高程783.871776(km)，该用户的空间位置是：经度28.917655°、纬度-99.345589°、高程0.0(km)，经过该地心地固坐标系转换，该卫星在该地心地固坐标系中的卫星地固坐标是(XS:6.439025706247130e+05，YS:-6.448794085263775e+06，ZS:
3.039310992962332e+06)，该用户在该地心地固坐标系中的用户地固坐标是(XU:-9.073090474340906e+05，YU:-5.513096509229902e+06，ZU:3.065915618754285e+06)，单位为米；经过该本地坐标系转换，该用户位置为该本地坐标系的原点，该卫星在该本地坐标系中的位置是xS:1.682569089856995e+06，yS:-3.479338007695247e+05，zS:5.748030724080029e+05，单位为米；经过该实时计算仰角θ，得到θ＝18.497377364764110°。
在本发明圆轨道低轨卫星过顶剩余可见时长预测方法另一实施例中，该计算卫星相对用户的过顶剩余可见时长的方法包括：(1)计算卫星已过顶时长，计算从该用户相对该卫星的最小可见仰角θ0出现的时刻到当前仰角为θ(te)＝θ时所经过的时长te，计算方法如下：
其中，RE和RS分别为地球的半径和该卫星的轨道半径，ωE和ωS分别为地球和该卫星在地心惯性坐标系下的角速度，γ0为该用户相对该卫星的仰角为θ0时该卫星星下点到该用户的角距离，计算得到te有两个解值和
(2)计算卫星过顶剩余可见时长tr，计算方法如下：
在本发明圆轨道低轨卫星过顶剩余可见时长预测方法另一实施例中，该迹角γm＝0.229891909587292(rad)，θ＝18.497377364764110，RE＝6378.137km，ωE＝0.00007292115，κS＝780km，i＝86.4°，θ0＝8.2°，μ＝398601.58(km3/s2)，RS＝RE+κS，计算得到：γ0＝0.347752358714960，ψ0＝0.264117321639204，tc＝5.089429154649187e+02，和
如果该卫星远离该用户，则该剩余过顶时间为：
如果该卫星靠近该用户，则该剩余过顶时间为：
本发明的有益效果是：本发明提供的圆轨道低轨卫星过顶剩余可见时长预测方法，能够充分利用卫星运行工作特点，不断计算更新卫星空间位置及星下点位置以及对应的时刻信息，并以简易的导航信号的形式向地面用户广播，地面用户接收该导航信号后，再利用本身的GPS接收设备校准时钟、确定用户位置，实时计算卫星的仰角和迹角，进而求解出卫星过顶剩余可见时长，具有实现方法简单有效，对导航信号的内容需求少，确保了预测时长的精准性和高时效性，又不增加用户硬件成本，推广应用价值高。
附图说明
图1是本发明圆轨道低轨卫星过顶剩余可见时长预测方法一个实施例的流程图；
图2是根据本发明圆轨道低轨卫星过顶剩余可见时长预测方法中用户接收卫星导航信号的一个实施例的流程图；
图3是根据本发明圆轨道低轨卫星过顶剩余可见时长预测方法中导航信号报文的一个实施例的组成图；
图4是根据本发明圆轨道低轨卫星过顶剩余可见时长预测方法中计算卫星仰角的一个实施例的流程图；
图5是根据本发明圆轨道低轨卫星过顶剩余可见时长预测方法中卫星与地面用户位置关系一个实施例的示意图；
图6是根据本发明圆轨道低轨卫星过顶剩余可见时长预测方法中计算卫星迹角的一个实施例的流程图；
图7是根据本发明圆轨道低轨卫星过顶剩余可见时长预测方法中计算卫星剩余可见时长的一个实施例的流程图。
具体实施方式
为了便于理解本发明，下面结合附图和具体实施例，对本发明进行更详细的说明。附图中给出了本发明的较佳的实施例。但是，本发明可以以许多不同的形式来实现，并不限于本说明书所描述的实施例。相反地，提供这些实施例的目的是使对本发明的公开内容的理解更加透彻全面。
需要说明的是，除非另有定义，本说明书所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。在本发明的说明书中所使用的术语只是为了描述具体的实施例的目的，不是用于限制本发明。本说明书所使用的术语“和/或”包括一个或多个相关的所列项目的任意的和所有的组合。
图1显示了本发明圆轨道低轨卫星过顶剩余可见时长预测方法的一个实施例的流程图。首先，用户接收卫星导航信号S101，主要是地面用户接收过顶卫星发出的导航信号，这里的卫星不仅是作为通信卫星传输用户的通信信号，还要以一定的时间间隔发出导航信号，该导航信号中包含有该卫星在未来一段时间内的位置信息及对应的时刻信息，通过所接收的卫星导航信号可以判定该卫星是远离用户的卫星还是靠近用户的卫星。然后，计算出卫星相对用户的仰角S102和计算出卫星相对用户的迹角S103，这两个步骤是基于用户所处的位置信息以及用户所接收的导航信号包括的卫星位置信息来计算的，这两个步骤既可以并行完成也可以依次完成，得到的仰角和迹角用于进一步计算该卫星的过顶剩余可见时长。最后，通过计算卫星相对用户的过顶剩余可见时长S104，得到卫星相对用户的过顶剩余可见时长。
以下通过具体实施例对上述步骤S101、S102、S103和S104作进一步说明。
图2显示了本发明圆轨道低轨卫星过顶剩余可见时长预测方法中用户接收卫星导航信号的一个实施例的流程图。首先，计算和更新卫星空间位置信息S1011，其中，计算卫星空间位置信息是由地面控制中心完成的，主要是计算卫星在未来一段时间内的位置信息(包括经纬度、高程)及其对应的时刻信息，计算的方法采用SGP4(Simple GeneralPerturbation version 4)方法，美国国防部于1980年公布了SGP4的数学模型及其代码，该方法用于对近地轨道上运行周期小于255分钟的飞行物体的位置进行预测，特别是对于海拔1000千米的运动物体具有较好的预报精度。这里的未来一段时间主要是指卫星经过地面控制中心上空之后到再次经过地面控制中心上空的这段时间，而地面控制中心计算的频度是1次/秒，即每隔1秒计算1次卫星所在的空间位置。更新卫星空间位置信息则是在卫星到达地面控制中心上空时，由地面控制中心向该卫星注入下一个未来一段时间内的卫星空间位置信息；计算卫星空间位置，由地面控制中心通过轨道预报方法，定期计算所述卫星未来一段时间内的空间位置信息及其对应的时刻信息，并在所述卫星经过所述地面控制中心上空时，更新所述卫星的空间位置信息；然后，卫星校准时钟S1012，这里，卫星利用自身的GPS装置进行时钟校准，获得准确的时钟信息，并根据时钟信息及在步骤S1011中存储的卫星空间位置信息来实时确定该卫星的星下点位置信息，该星下点位置信息是指卫星当前覆盖范围内的星下点经纬度信息，这是卫星的空间位置信息在地球表面投影的位置信息；接着，卫星发出导航信号S1013，该导航信号中主要是指包括了该卫星在其所覆盖范围内各星下点组成的轨迹的经纬度信息及其对应的时钟信息，并且卫星以1次/秒的频度发送导航信号。图3显示了该导航信号的报文结构的一个实施例，其中(λ1、η1、κ1)表示在时刻timestamp1卫星空间位置对应的经度、纬度和高程，(λ2、η2、κ2)表示在时刻timestamp2卫星空间位置对应的经度、纬度和高程，以此类推，(λN、ηN、κN)表示在时刻timestampN卫星空间位置对应的经度、纬度和高程，并且时刻timestamp1、timestamp2…timestampN可以是依次相差1秒，表示卫星以1次/秒的频度发送导航信号；用户接受导航信号并判断卫星远离还是接近S1014，地面用户不断接收卫星发出的导航信号，其中可以利用连续相邻的两次卫星导航信号，计算判断该卫星是接近用户还是远离用户，该计算过程主要是通过计算两个星下点位置相对用户位置的距离远近进行判断，若第二个星下点位置比第一个星下点位置距离用户位置更近，则表明卫星在接近用户，否则卫星远离用户。
通过图2所示流程图及图3所示导航报文结构图可以看出，本发明圆轨道低轨卫星过顶剩余可见时长预测方法中的卫星能够发出导航信号，但这种导航信号具有简易性，只需包括卫星的经度、纬度和高程信息及对应的时刻信息，可以作为该卫星的一种广播信号进行播发，能够与这种通信卫星的通信体制相兼容，而不需要像GPS等定位导航卫星那样发出复杂的卫星星历信息，因此更为适合这种用于通信的卫星。
图4显示了本发明圆轨道低轨卫星过顶剩余可见时长预测方法中计算卫星仰角的一个实施例的流程图。首先，进行地心地固坐标系转换S1021，就是将根据当前卫星的空间位置信息和用户的空间位置信息(主要是的经纬度和高程信息)统一转换为地心地固(ECEF)坐标系中的位置信息，通常是用空间三维坐标表示。地心地固坐标系的具体转换方法是：卫星的空间位置信息为(λS,ηS,κS)，其中λS表示该卫星的经度，ηS表示该卫星的纬度，κS表示该卫星的高程，用户的空间位置信息为(λT,ηT,κT)，其中λT表示该用户的经度，ηT表示该用户的纬度，κT表示该用户的高程，该卫星地固坐标为(XS,YS,ZS)，该用户地固坐标为(XU,YU,ZU)，地球是一个椭球体，长轴a＝6.3781370e+03km，半轴b＝6.356752314245179e+03km，偏心率定义中间值则该卫星的地固坐标(XS,YS,ZS)和该用户的地固坐标(XU,YU,ZU)分别为：
例如，按照上述地心地固坐标系转换方法，某一时刻卫星的空间位置是：经度25.257000°，纬度-84.298000°，高程783.871776(km)，转换到地心地固(ECEF)坐标系中的卫星地固坐标为(XS:
6.439025706247130e+05，YS:-6.448794085263775e+06，ZS:
3.039310992962332e+06)，单位为米。又如，用户的空间位置是：经度28.917655°，纬度-99.345589°，高程0.0(km)，转换到地心地固(ECEF)坐标系中的用户地固坐标(XU:
-9.073090474340906e+05，YU:-5.513096509229902e+06，ZU:
3.065915618754285e+06)，单位为米。
然后，进行本地坐标系转换S1022，就是根据用户实际使用的本地坐标系情况，将S1021中的地心地固(ECEF)坐标系中的卫星和用户位置信息转换为本地坐标系(ENU)中的卫星和用户位置信息，这实际上是地心地固(ECEF)坐标系与本地坐标系的转换，转换时，将用户位置作为本地坐标系的原点，因此在本地坐标系中用户在该坐标系的原点，称为用户原点坐标，据此确定卫星在本地坐标系中相对于用户的位置坐标。转换方法是：若卫星的经度为λS，纬度为ηS，该卫星在地心地固坐标系的卫星地固坐标表示为(XS,YS,ZS)，用户地固坐标为(XU,YU,ZU)，则该卫星在以用户为原点的本地坐标系中的坐标表示为(xS,yS,zS)，则有：
例如，上述某一时刻卫星的空间位置是：经度λS是25.257000°，纬度ηS是-84.298000°，该卫星在地心地固(ECEF)坐标系中的卫星地固坐标是(XS:6.439025706247130e+05，YS:
-6.448794085263775e+06，ZS:3.039310992962332e+06)，用户的地心地固(ECEF)坐标系中的用户地固坐标是(XU:
-9.073090474340906e+05，YU:-5.513096509229902e+06，ZU:
3.065915618754285e+06)，转换到本地坐标系后，该卫星在以用户为原点的本地坐标系中的坐标为：xS:1.682569089856995e+06，yS:-3.479338007695247e+05，zS:5.748030724080029e+05，单位为米；实时计算仰角S1023，在本地坐标系中计算当前时刻用户观测卫星的仰角θ，该仰角是实时变化的，例如根据上述卫星和用户的本地坐标位置可以计算当前时刻θ＝18.497377364764110°。
由此可见，为了精确计算用户观测卫星的仰角，通过对卫星空间位置和用户空间位置进行坐标转换，将二者转换为同一个坐标系中计算仰角，这样有助于消除因坐标系不同带来的误差，并且当所选择的坐标系之间有固定的几何转换关系后，可以利用计算机进行自动转换，提高了转换计算效率。
图5显示了本发明圆轨道低轨卫星过顶剩余可见时长预测方法中卫星与地面用户位置关系一个实施例的示意图。其中，卫星S当前覆盖范围内的星下点的空间位置是以矢量形式表示，说明卫星在处于不断运动的状态，其星下点H位置是(λS、ηS、κ'S)，用户T的空间位置是(λT、ηT、κT)，γ(t)表示用户T到当前卫星S星下点H的最短弧距，而Q点表示当卫星S经过用户T上空时，该卫星在地球的星下点轨迹到用户T的最短弧距与该星下点轨迹的交点，γm为用户T的迹角，ψ(t)表示当前卫星S星下点H同Q点之间的最短弧距，θ即为当前用户T相对卫星S的仰角。
图6显示了本发明圆轨道低轨卫星过顶剩余可见时长预测方法中计算卫星迹角的一个实施例的流程图。首先，确定用户的位置信息S1031，可以通过用户设备内置的GPS接收机计算得到该用户的经纬度和高程信息，分别表示为用户经度λT、用户纬度ηT和用户高程κT；时钟校准S1032，从接收的卫星导航信号中获取精度较高的卫星时钟，将用户时钟同步校准与卫星时钟一致；计算迹角S1033，从接收的卫星导航信号中实时得到该卫星的星下点位置信息，即星下点经纬度信息，以动态向量形式表示为卫星星下点经度卫星星下点纬度表明该星下点经纬度是随时间在不断变化，用户相对卫星的迹角表示为γm，先计算得到：
再进一步计算得到：
γm＝min{2sin-1(Ω)}
例如，当用户的经度λT＝28.917655°、纬度ηT＝-99.345589°不变的情况下，而卫星的星下点经度和纬度随着卫星过顶运动而不断改变，可以通过列表的形式表示星下点的经纬度改变情况以及对应计算得到的Ω值，如表1部分所示。然后，再对所有的Ω值求其反正弦函数，找到其中的最小值，再乘以2，最终可以得到γm＝0.229891909587292(rad)。
表1卫星过顶时星下点经纬度及Ω值变化情况
图7显示了本发明圆轨道低轨卫星过顶剩余可见时长预测方法中计算卫星剩余可见时长的一个实施例的流程图。首先，计算已过顶时长S1034，就是计算从用户看到卫星的时刻起到当前仰角为θ(te)＝θ时所经过的时间te，这里用户看到卫星的时刻是指用户相对卫星的最小可见仰角θ0出现的时刻，计算方法如下：
其中，RE为地球半径，RS为卫星的圆形轨道半径，ωE是地球的角速度，ωS是卫星的角速度，γ0为用户在相对卫星的仰角为θ0时的卫星星下点到用户的角距离。
可以看出，te是通过求反正切三角函数得到的，而同一个正切函数对应的角度值有两个，因此对应的te存在两个解和
进一步，计算过顶剩余时长S1034，剩余时长表示为tr，计算方法如下：
例如，在迹角γm＝0.229891909587292(rad)和θ＝18.497377364764110的情况下，RE＝6378.137km，ωE＝0.00007292115，RS＝RE+κS，κS＝780km，i＝86.4°，θ0＝8.2°，ω＝(ωS-ωEcosi)/2，μ＝398601.58(km3/s2)，则可得到：
γ0＝0.347752358714960，ψ0＝0.264117321639204，
tc＝5.089429154649187e+02，
和单位为秒，如果卫星远离用户，则剩余过顶时间为：
如果卫星靠近用户，则剩余过顶时间为：
通过本发明圆轨道低轨卫星过顶剩余可见时长预测方法，能够根据卫星运行工作特点，通过不断计算和更新卫星空间位置、星下点位置以及对应的时刻信息，并以简易的导航信号的形式向地面用户广播这些信息，地面用户接收该导航信号后，再利用本身的GPS接收设备校准时钟、确定用户位置，实时计算卫星相对用户的仰角和迹角，进而求解出卫星过顶剩余可见时长，具有实现方法简单有效，对导航信号的内容需求少，对卫星星历依赖度低，在确保预测时长的精准性和高时效性前提下，又不增加用户硬件成本。
以上所述仅为本发明的实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构变换，或直接或间接运用在其他相关的技术领域，均包括在本发明的专利保护范围内。
Claims (9)
1.一种圆轨道低轨卫星过顶剩余可见时长预测方法，其特征在于，包括如下步骤：
第一步，用户接收卫星导航信号，包括：
(1)计算卫星空间位置，由地面控制中心通过轨道预报方法，定期计算所述卫星未来一段时间内的空间位置信息及其对应的时刻信息，并在所述卫星经过所述地面控制中心上空时，更新所述卫星的空间位置信息，所述卫星的空间位置信息包括经纬度、高程及其对应的时刻信息；
(2)卫星校准时钟，所述卫星利用自身GPS装置进行时钟校准，并根据校准的时钟信息以及存储的所述卫星的空间位置信息，实时调整确定所述卫星当前覆盖范围内的星下点经纬度信息；
(3)卫星发出导航信号，所述卫星定期发出导航信号，所述导航信号包括所述卫星在其所覆盖范围内的星下点经纬度信息及其对应的时刻信息；
(4)用户接收导航信号，判断所述卫星是接近用户还是远离用户，所述用户接收所述卫星发送的连续两次所述导航信号，根据所述导航信号中的星下点经纬度信息相对于用户位置的远近变化判断所述卫星接近还是远离所述用户；
第二步，计算卫星的仰角；
第三步，计算卫星的迹角；
第四步，计算卫星相对用户的过顶剩余可见时长。
2.根据权利要求1所述的圆轨道低轨卫星过顶剩余可见时长预测方法，其特征在于，所述轨道预报方法为SGP4方法；所述定期计算所述卫星未来一段时间内的空间位置信息及其对应的时刻信息的频度是每间隔1秒钟计算1次；所述卫星定期发出导航信号的频度是每间隔1秒钟发送1次。
3.根据权利要求1或2所述的圆轨道低轨卫星过顶剩余可见时长预测方法，其特征在于，所述计算卫星的迹角的方法包括：
(1)确定所述用户位置信息，得到所述用户自身的经度λT和纬度ηT信息；
(2)时钟校准，根据所接收导航信号中所述时刻信息，将用户时钟同步校准与卫星时钟一致；
(3)实时计算迹角，所述用户通过接收的所述导航信号，计算所述用户相对所述卫星的迹角γm，计算方法如下：
<mrow>
<mi>&Omega;</mi>
<mo>=</mo>
<mo>&lsqb;</mo>
<msqrt>
<mrow>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mover>
<msub>
<mi>&lambda;</mi>
<mi>S</mi>
</msub>
<mo>&RightArrow;</mo>
</mover>
<mo>-</mo>
<msub>
<mi>&lambda;</mi>
<mi>T</mi>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mover>
<msub>
<mi>&lambda;</mi>
<mi>S</mi>
</msub>
<mo>&RightArrow;</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&lambda;</mi>
<mi>T</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mover>
<msub>
<mi>&eta;</mi>
<mi>S</mi>
</msub>
<mo>&RightArrow;</mo>
</mover>
<mo>-</mo>
<msub>
<mi>&eta;</mi>
<mi>T</mi>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
<mo>&rsqb;</mo>
</mrow>
γm＝min{2sin-1(Ω)}
其中，和是所述导航信号中包含的所述卫星覆盖范围内的星下点的经度和纬度向量，Ω也为向量。
4.根据权利要求3所述的圆轨道低轨卫星过顶剩余可见时长预测方法，其特征在于，通过所述用户内置的GPS设备测量得到所述用户自身的经度λT、纬度ηT和高程κT信息。
5.根据权利要求4所述的圆轨道低轨卫星过顶剩余可见时长预测方法，其特征在于，所述计算卫星的仰角的方法包括：
(1)地心地固坐标系转换，将当前所述卫星和用户的空间位置信息统一转换为地心地固坐标系中卫星地固坐标和用户地固坐标；
(2)本地坐标系转换，将所述卫星地固坐标和所述用户地固坐标转换为本地坐标系中的卫星本地坐标和用户原点坐标；
(3)实时计算仰角，利用所述卫星本地坐标和所述用户原点坐标，在所述本地坐标系中计算所述用户相对所述卫星的仰角θ。
6.根据权利要求5所述的圆轨道低轨卫星过顶剩余可见时长预测方法，所述卫星的空间位置信息为(λS,ηS,κS)，其中λS表示所述卫星的经度，ηS表示所述卫星的纬度，κS表示所述卫星的高程，所述用户的空间位置信息为(λT,ηT,κT)，其中λT表示所述用户的经度，ηT表示所述用户的纬度，κT表示所述用户的高程，所述卫星地固坐标为(XS,YS,ZS)，所述用户地固坐标为(XU,YU,ZU)，其特征在于，所述地心地固坐标系转换方法是：
地球是一个椭球体，长轴a＝6.3781370e+03km，半轴b＝6.356752314245179e+03km，偏心率定义中间值则所述卫星的地固坐标(XS,YS,ZS)和所述用户的地固坐标(XU,YU,ZU)分别为：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>X</mi>
<mi>S</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>Normal</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&kappa;</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>cos&lambda;</mi>
<mi>S</mi>
</msub>
<msub>
<mi>cos&eta;</mi>
<mi>S</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Y</mi>
<mi>S</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>Normal</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&kappa;</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>cos&lambda;</mi>
<mi>S</mi>
</msub>
<msub>
<mi>sin&eta;</mi>
<mi>S</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mi>S</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>Normal</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>&epsiv;</mi>
<mn>2</mn>
</msup>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>&kappa;</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>sin&lambda;</mi>
<mi>S</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>X</mi>
<mi>U</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>Normal</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&kappa;</mi>
<mi>T</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>cos&lambda;</mi>
<mi>T</mi>
</msub>
<msub>
<mi>cos&eta;</mi>
<mi>T</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Y</mi>
<mi>U</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>Normal</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&kappa;</mi>
<mi>T</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>cos&lambda;</mi>
<mi>T</mi>
</msub>
<msub>
<mi>sin&eta;</mi>
<mi>T</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mi>U</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>Normal</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>&epsiv;</mi>
<mn>2</mn>
</msup>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>&kappa;</mi>
<mi>T</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>sin&lambda;</mi>
<mi>T</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
所述本地坐标系转换方法是：所述卫星在所述用户原点坐标的所述本地坐标系中的卫星本地坐标(xS,yS,zS)是：
<mrow>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mi>S</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>y</mi>
<mi>S</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>z</mi>
<mi>S</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>sin&eta;</mi>
<mi>S</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>cos&eta;</mi>
<mi>S</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>sin&lambda;</mi>
<mi>S</mi>
</msub>
<msub>
<mi>cos&eta;</mi>
<mi>S</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>sin&lambda;</mi>
<mi>S</mi>
</msub>
<msub>
<mi>cos&eta;</mi>
<mi>S</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>cos&lambda;</mi>
<mi>S</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>cos&lambda;</mi>
<mi>S</mi>
</msub>
<msub>
<mi>cos&eta;</mi>
<mi>S</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>cos&lambda;</mi>
<mi>S</mi>
</msub>
<msub>
<mi>cos&eta;</mi>
<mi>S</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>sin&lambda;</mi>
<mi>S</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>X</mi>
<mi>S</mi>
</msub>
<mo>-</mo>
<msub>
<mi>X</mi>
<mi>U</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Y</mi>
<mi>S</mi>
</msub>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mi>U</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mi>S</mi>
</msub>
<mo>-</mo>
<msub>
<mi>Z</mi>
<mi>U</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
7.根据权利要求6所述的圆轨道低轨卫星过顶剩余可见时长预测方法，其特征在于，所述卫星的空间位置是：经度25.257000°、纬度-84.298000°、高程783.871776(km)，所述用户的空间位置是：经度28.917655°、纬度-99.345589°、高程0.0(km)，经过所述地心地固坐标系转换，所述卫星地固坐标是(XS:6.439025706247130e+05，YS:-6.448794085263775e+06，ZS:3.039310992962332e+06)，所述用户地固坐标是(XU:-9.073090474340906e+05，YU:-5.513096509229902e+06，ZU:3.065915618754285e+06)，单位为米；经过所述本地坐标系转换，所述用户位置为所述本地坐标系的原点，所述卫星在所述本地坐标系中的位置是xS:1.682569089856995e+06，yS:-3.479338007695247e+05，zS:5.748030724080029e+05，单位为米；经过所述实时计算仰角θ，得到θ＝18.497377364764110°。
8.根据权利要求5或6所述的圆轨道低轨卫星过顶剩余可见时长预测方法，其特征在于，所述计算卫星相对用户的过顶剩余可见时长的方法包括：
(1)计算卫星已过顶时长，计算从所述用户相对所述卫星的最小可见仰角θ0出现的时刻到当前仰角为θ(te)＝θ时所经过的时长te，计算方法如下：
<mrow>
<msub>
<mi>&gamma;</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msup>
<mi>cos</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>R</mi>
<mi>E</mi>
</msub>
<msub>
<mi>R</mi>
<mi>S</mi>
</msub>
</mfrac>
<msub>
<mi>cos&theta;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<mi>&psi;</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>&psi;</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<msup>
<mi>cos</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>cos&gamma;</mi>
<mn>0</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>cos&gamma;</mi>
<mi>m</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&theta;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>e</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>tan</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mo>(</mo>
<msub>
<mi>&omega;</mi>
<mi>E</mi>
</msub>
<mi>cos</mi>
<mi>i</mi>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>t</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&psi;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
<mo>&times;</mo>
<msub>
<mi>cos&gamma;</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<mfrac>
<msub>
<mi>R</mi>
<mi>E</mi>
</msub>
<msub>
<mi>R</mi>
<mi>S</mi>
</msub>
</mfrac>
</mrow>
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<mo>(</mo>
<msub>
<mi>&omega;</mi>
<mi>E</mi>
</msub>
<mi>cos</mi>
<mi>i</mi>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mi>S</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>t</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&psi;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
<mo>&times;</mo>
<msub>
<mi>cos&gamma;</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
其中，RE和RS分别为地球的半径和所述卫星的轨道半径，ωE和ωS分别为地球和所述卫星在地心惯性坐标系下的角速度，γ0为所述用户相对所述卫星的仰角为θ0时所述卫星星下点到所述用户的角距离，计算得到te有两个解值和
(2)计算卫星过顶剩余可见时长tr，计算方法如下：
ω＝(ωS-ωEcosi)/2
9.根据权利要求8所述的圆轨道低轨卫星过顶剩余可见时长预测方法，其特征在于，所述迹角γm＝0.229891909587292(rad)，θ＝18.497377364764110，RE＝6378.137km，ωE＝0.00007292115，κS＝780km，i＝86.4°，θ0＝8.2°，μ＝398601.58(km3/s2)，RS＝RE+κS，计算得到：
γ0＝0.347752358714960，ψ0＝0.264117321639204，
tc＝5.089429154649187e+02，
和
如果所述卫星远离所述用户，则所述剩余过顶时间为：
<mrow>
<msub>
<mi>t</mi>
<mi>r</mi>
</msub>
<mo>=</mo>
<msub>
<mi>t</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>t</mi>
<mi>e</mi>
<mn>1</mn>
</msubsup>
<mo>=</mo>
<mn>1.800178832659187</mn>
<mi>e</mi>
<mo>+</mo>
<mn>02</mn>
<mo>;</mo>
</mrow>
如果所述卫星靠近所述用户，则所述剩余过顶时间为：
<mrow>
<msub>
<mi>t</mi>
<mi>r</mi>
</msub>
<mo>=</mo>
<msub>
<mi>t</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>t</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<mn>3.289250321989186</mn>
<mi>e</mi>
<mo>+</mo>
<mn>02.</mn>
</mrow>
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510417367.8A CN105044745B (zh) | 2015-07-15 | 2015-07-15 | 一种圆轨道低轨卫星过顶剩余可见时长预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510417367.8A CN105044745B (zh) | 2015-07-15 | 2015-07-15 | 一种圆轨道低轨卫星过顶剩余可见时长预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105044745A CN105044745A (zh) | 2015-11-11 |
| CN105044745B true CN105044745B (zh) | 2018-01-12 |
Family
ID=54451419
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510417367.8A Active CN105044745B (zh) | 2015-07-15 | 2015-07-15 | 一种圆轨道低轨卫星过顶剩余可见时长预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105044745B (zh) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109495160B (zh) * | 2018-12-04 | 2021-04-02 | 航天科工空间工程发展有限公司 | 一种低轨通信卫星星座与信关站连通规划方法 |
| CN110519695B (zh) * | 2019-05-31 | 2020-12-11 | 中国人民解放军国防科技大学 | 一种数据库辅助的卫星***与地面蜂窝网络频谱共享方法 |
| CN110907954B (zh) * | 2019-10-22 | 2022-07-15 | 深圳航天东方红海特卫星有限公司 | 一种标校数据广播装置 |
| CN111294962B (zh) * | 2020-02-26 | 2022-12-20 | 展讯通信（上海）有限公司 | 一种数据传输方法及装置、存储介质、终端 |
| CN111866119B (zh) * | 2020-07-16 | 2022-06-14 | 湖南斯北图科技有限公司 | 基于卫星轨道预报算法的物联网数据采集终端节能方法 |
| CN113687392B (zh) * | 2021-08-23 | 2024-04-12 | 深圳乐心信息技术有限公司 | 一种基于gnss信号非连续跟踪的导航方法 |
| CN116073878A (zh) * | 2021-10-29 | 2023-05-05 | 华为技术有限公司 | 卫星通信方法和卫星通信装置 |
| CN114828134B (zh) * | 2022-04-25 | 2023-03-24 | 广州爱浦路网络技术有限公司 | 通信卫星切换方法、计算机装置及存储介质 |
| CN115032671A (zh) * | 2022-08-11 | 2022-09-09 | 成都国星宇航科技股份有限公司 | 一种低轨卫星跟踪预报时间段计算方法和装置 |
| CN115292405B (zh) * | 2022-10-10 | 2022-12-13 | 中科星图测控技术(合肥)有限公司 | 一种计算测站跟踪预报过顶点的显示方法和*** |
| CN117421518B (zh) * | 2023-12-18 | 2024-03-22 | 南京融星智联信息技术有限公司 | 一种低轨卫星对地面覆盖时间窗口的计算方法及*** |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103453906B (zh) * | 2013-08-09 | 2016-04-27 | 清华大学 | 卫星轨道的预测方法 |
| CN104298647B (zh) * | 2014-09-30 | 2017-02-15 | 北京航空航天大学 | 基于低轨道地球卫星的地影时刻预报的星上确定方法 |
-
2015
- 2015-07-15 CN CN201510417367.8A patent/CN105044745B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105044745A (zh) | 2015-11-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105044745B (zh) | 一种圆轨道低轨卫星过顶剩余可见时长预测方法 | |
| CN1995917B (zh) | 一种电子地图俯视角度显示控制方法、装置及设备 | |
| CN106792946B (zh) | 卫星网络中一种面向快速终端的自适应切换方法 | |
| CN102176030B (zh) | 一种gps和glonass双***组合模拟器及方法 | |
| CN102135433B (zh) | 精确导航方法及位置服务设备 | |
| CN112558125B (zh) | 一种车辆定位的方法、相关装置、设备以及存储介质 | |
| CN100480728C (zh) | 确定接收机的位置和/或定位***的***时间 | |
| CN105973244A (zh) | 获取地图信息的方法、导航方法及设备 | |
| CN109031377A (zh) | 一种基于伪卫星的隧道内定位方法 | |
| CN101719802B (zh) | 一种预测短波通信最大可用频率的装置及计算方法 | |
| CN105580009A (zh) | 提供指示所推荐通航路段的信息的方法和*** | |
| CN107182124A (zh) | 一种基于网格的高精度定位方法及装置 | |
| CN1995918B (zh) | 一种导航路径的调整方法、装置及设备 | |
| US7215335B2 (en) | Image generating apparatus, image generating method, and computer program | |
| Ta | Global Positioning System | |
| CN201936013U (zh) | 一种gps和glonass双***组合模拟器 | |
| Ragothaman | Path planning for autonomous ground vehicles using GNSS and cellular LTE signal reliability maps and GIS 3-D maps | |
| CN104471439B (zh) | 位置计算方法以及位置计算装置 | |
| US20110282573A1 (en) | Route planning method | |
| JP7338691B2 (ja) | 信号受信機器、信号処理方法及びプログラム | |
| JP2004156999A (ja) | 測位方法および測位装置 | |
| Tsujino | Effectiveness of the Quasi-Zenith Satellite System in Ubiquitous Positioning | |
| CN103544677A (zh) | 空天地一体化态势表达引擎及其抖动消除方法 | |
| US20050216194A1 (en) | Method for automatically executing navigation software | |
| Sonowal et al. | Real time GPS software receiver with new fast signal tracking method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |