GPS卫星可视性研究_高建东.pdf

河南科技Henan Science and Technology电气与信息工程总第800期第6期2023年3月GPS卫星可视性研究高建东1邢尹2（1.江苏省测绘工程院，江苏南京210013；2.苏州科技大学地理科学与测绘工程学院，江苏苏州215009）摘要：【目的目的】研究相对于地面观测点的 GPS 卫星星座的可视性，从而改善观测数据质量，提高测量精度。【方法方法】本研究基于标准RINEX数据，通过假设卫星开普勒参数，拟定地面观测点及观测时间，并根据全球导航卫星系统定位原理来计算各 GPS卫星坐标及相对于地面观测点的高度角，从而筛选出符合观测条件的 GPS 卫星。【结果结果】在计算出坐标及相对于地面观测点的方位角和高度角的 24颗 GPS卫星中，选出 7颗高度角大于 15的卫星作为有效的观测卫星。【结论结论】本研究所述的方法可按照一定原则对 GPS卫星星座进行筛选，从而获取质量更佳的卫星观测数据，为相关研究提供参考。关键词：GPS；RINEX；卫星可视性；卫星高度角中图分类号：P228文献标志码：A文章编号：1003-5168（2023）06-0018-05DOI：10.19968/ki.hnkj.1003-5168.2023.06.003Research on Visibility of GPS SatelliteGAO Jiandong1XING Yin2(1.Jiangsu Institute of Surveying and Mapping Engineering,Nanjing 210013,China;2.School of GeographicScience and Surveying Engineering,Suzhou University of Science and Technology,Suzhou 215009,China)Abstract:Purposes This paper study the visibility of GPS satellite constellation relative to ground observation points,so as to improve the quality of observation data and improve the measurement accuracy.Methods In this study,based on the standard RINEX data and by assuming the Kepler parameters ofthe satellite,the ground observation points and observation time are proposed,and the GPS satellite coordinates and the elevation angle relative to the ground observation points are calculated according to thepositioning principle of the global navigation satellite system,so as to select the GPS satellite that meetsthe observation conditions.Findings In accordance with the coordinates of 24 GPS satellites and theazimuth and elevation angles relative to the ground observation points,7 satellites with an elevationangle greater than 15 are selected as effective observation satellites.Conclusions The method described in this study can sift the GPS satellite constellation according to certain principles,so as to obtain better quality satellite observation data and provide reference for related research.Keywords:GPS;RINEX;satellite visibility;satellite altitude angle0引言全球导航卫星系统（Global Navigation SatelliteSystem，GNSS）是可为地球用户或空间用户提供定位和授时信息的卫星星座1，其在很大程度上能反映一个国家的综合科技实力水平，是衡量一个国家科技发展的重要指标。全球定位系统（GPS）是美国在二十世纪七十年代开始建设的、世界上收稿日期：2022-12-27作者简介：高建东（1989），男，硕士，工程师，研究方向：大地测量、卫星导航定位。第6期19高建东，等.GPS卫星可视性研究第一个用于导航定位的全球系统，经过几十年的不断发展和更新，目前 GPS 已成为全球星座组网最完善、定位精度最高、用户数量最多的卫星导航系统2。在利用GNSS来进行控制测量时，导航卫星的分布状态、高度角等都会影响测量精度。根据观测点的位置，计算并预测卫星的高度角、方位角，从而筛选出有效的观测卫星，能有效改善观测数据的质量，极大提高GPS卫星对地面测量的精度。随着卫星导航技术的快速发展，对各类卫星的方位角和高度角的计算也越来越重要。在卫星可视性分析领域中，国内外相关学者已进行非常深入的研究。早在2000年，西安测绘研究所的贾小林教授就对卫星对地的可视性进行研究，其对卫星系统的可用性及如何计算进行详细论述。2007年10月“嫦娥一号”成功发射后，李晓杰等3就发表了 环月卫星可见时段的计算和分析，详细深入地解释了可见时段的计算原理，并根据地面站的布设方式的不同，采用不同的方法计算。除此之外，2013年魏东岩等4在第四届中国卫星导航学术年会上发表了 针对高轨平台的GNSS卫星可视性分析研究，根据高轨平台的定位需求，从多种角度出发，对卫星信号的可视性进行分析。1卫星的可视性分析原理1.1标准RINEX格式由于生产厂家的不同，导致 GPS 卫星接收到的数据在格式上存在巨大的差异，从而使数据处理较为困难。为了对GPS卫星接收数据的格式进行统一，瑞士伯尔尼大学天文学院的 Gurtner 于1989 年提出 Receiver Independent Exchange Format（RINEX）5。之后，所有的卫星数据便能以统一的形式呈现在世人的面前，方便对各种型号接收机的 数 据 进 行 综 合 处 理，极 大 提 高 工 作 效 率。RINEX 文件是一个统一的格式文件，其包含了卫星观测数据的文件、导航电文文件和气象数据文件6。其中，导航电文是 RINEX 文件中负责保存GPS 星历信息的文件，且轨道参数、星历参考时刻、钟差改正数及摄动改正数都存放在 RINEX 文件中。1.2卫星坐标计算导航电文包含卫星某时刻的具体空间位置信息7，从中读取所需的参数，再通过以下计算流程，就能得到卫星的坐标。1.2.1解算卫星运行的平均角速度n0。由开普勒第三定律可知，平均角速度n0的计算公式见式（1）。n0=GM/a3=/()a3（1）式中：G 为重力常数；在 WGS-84 坐标系中，为地球的引力常数，=3.986 0051014m3/s2。卫星的平均角速度 n通常为 n0与摄动改正数n之和，计算公式见式（2）。n=n0+n（2）1.2.2解算归化时间tk。卫星轨道参数受参考时刻的影响，其会随着选择时刻的改变而发生变化。先求出参考时刻与观测时刻的差值，从而得到归化时间tk。观测时刻卫星瞬时的平近点角Mk的计算见式（3）。Mk=M0+ntk（3）式中：M0为参考时刻的平近点角；tk为归化时间。根据公式（3）可求出平近点角Mk。偏近点角Ek的计算见式（4）。Ek=Mk+esinEk（4）可利用迭代法来求出偏近点角Ek，使Ek=Mk，求出Ek，再代入即可求出Mk。根据公式（4）中偏近点角Ek，可计算出真近点角Vk，其计算见式（5）。Vk=tan-1|()1-e2 sinEk/()cosEk-e（5）根据公式（5）算出的真近点角Vk，升交距角k的计算见式（6）。k=Vk+（6）式中：为近地点角距。对摄动改正项进行求解，计算公式见式（7）。|u=Cuc cos()2k+Cus sin()2kr=Crc cos()2k+Crs sin()2ki=Cic cos()2k+Cis sin()2k（7）式中：u为升交距角u的摄动量；r为矢径r的摄动量；i为倾角i的摄动量。对 改 正 后 的 升 交 距 角 uk、rk、ik进 行 计 算，见式（8）。|uk=k+urk=()1-ecosEk+rik=i0+r+Itk（8）对卫星在轨道平面坐标系的坐标进行计算，见式（9）。20第6期xk=rkcosukyk=rksinuk（9）根据公式（8）和公式（9）计算出的结果，可算出观测时刻升交点精度k，计算公式见式（10）。k=0+()-etk-etoe（10）式中：e为地球自转速率，e=7.292 115 6710-5rad/s；为升交点对时间的变化率。1.2.3卫星坐标计算。根据公式（10）的计算结果，可计算出卫星在WGS-84坐标系中的坐标。对轨道平面坐标系中的坐标进行变换，就能得到在WGS-84坐标系中的坐标，见式（11）。|xk=xkcosk-ykcosiksinkyk=xksink+ykcosikcoskzk=yksinik（11）经过以上步骤的计算，可得到卫星在某时刻的具体空间坐标。2实例及分析2.1参数计算在对卫星的位置坐标进行计算前，要拟定好24颗卫星的开普勒参数。首先要计算出卫星运行的平均角速度，在对坐标进行计算时，有一个非常重要的步骤，即卫星归化时间的计算，采用之前已拟定好的时刻（2021年4月29日4时0分0秒），将时间系统统一后再进行计算，在转化好时间后，可计算出每个卫星的平近点角和偏近点角，要注意的是，偏近点角要进行迭代计算。在此之后，还要计算出地心矢径、真近点角、升交点角距、倾角，最后可通过已知的公式计算得到卫星在轨道坐标系中的坐标，并通过转换得到其在WGS-84大地坐标系中的坐标。当得到卫星在 WGS-84 坐标系中的坐标后，拟定观测点的地心坐标为（3 757 014 m，980 682 m，6 378 137 m），经 过 转 换 后，得 到 的 经 度 为119.190 042 494、纬度为 33.830 618 060。然后利用已知的测站坐标和24颗卫星的地心坐标，将卫星的地心坐标通过计算公式转换成站心坐标，在求出各卫星的站心坐标后，可计算出卫星的方位角和高度角。2.2计算实例及分析为了计算出卫星的坐标，要利用到开普勒参数，拟定的开普勒参数见表1。在表1给出的开普勒参数基础上，计算得到卫星在轨道平面内的坐标，结果见表2。经过计算公式的转换，可得到卫星在WGS-84高建东，等.GPS卫星可视性研究卫星123456789101112131415161718192021222324e0.007 30.007 30.007 30.007 30.007 30.007 30.007 30.007 30.007 30.007 30.007 30.007 30.007 30.007 30.007 30.007 30.007 30.007 30.007 30.007 30.007 30.007 30.007 30.007 3/m26 522 50026 522 50026 522 50026 522 50026 522 50026 522 50026 52