CN115685266B - 基于全球网格的gps信号分量调整中心点确定方法 - Google Patents
基于全球网格的gps信号分量调整中心点确定方法 Download PDFInfo
- Publication number
- CN115685266B CN115685266B CN202211295087.0A CN202211295087A CN115685266B CN 115685266 B CN115685266 B CN 115685266B CN 202211295087 A CN202211295087 A CN 202211295087A CN 115685266 B CN115685266 B CN 115685266B
- Authority
- CN
- China
- Prior art keywords
- center point
- satellite
- signal component
- grid
- component adjustment
- 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
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种基于全球网格的GPS信号分量调整中心点确定方法,该方法主要是根据IGS或iGMAS跟踪站数据,得出GPS卫星信号载噪比变化情况,统计载噪比增高和降低的操作时刻,利用卫星星历计算卫星位置,再将全球划分为若干个网格点,计算网格点对该卫星的仰角,通过统计的方法确定GPS系统所选择的中心点。对于利用中心点规划卫星信号分量调整的模式,实施本方法可以找到相应的中心点。本发明首次涉及GPS系统信号分量调整所选中心点的确定方法,利用公开的数据,通过简单的统计计算即可确定出中心点;研究结果有助于增进对GPS系统信号调整时机、任务规划能力的了解。
Description
技术领域
本发明涉及卫星导航领域,特别是一种基于全球网格的GPS信号分量调整中心点确定方法。
背景技术
GPS系统能为各类用户提供高精度定位、授时以及测速服务,但由于GPS卫星与地面距离远,下行信号落地功率低,如果接收机附近的电磁环境复杂,则可能导致用户接收信号失锁。GPS现代化计划提出后,GPS IIR-M、GPS IIF以及GPS III卫星都具备星上信号分量的调整能力,能够根据需要调整卫星下播的信号分量,调整过程中可以将一部分信号的能量分配给指定的信号,由于指定信号分量获得了更大的能量,其落地信号的功率也将增高。通过信号分量调整的方式能够使GPS指定信号的下行功率增大,抵消一部分电磁环境的不利影响。例如,为减小地面对GPS下行信号的压制式干扰,2018年GPS系统曾通过调整不同信号能量配比的方式使P(Y)码落地信号功率得到增高。
中国发明专利CN201810922419.0一种卫星导航信号高精度时域性能评估方法,采集卫星导航信号,对卫星导航信号包含的N个信号分量进行数据预处理;根据预处理结果将导航信号转换为基带信号,获取需要恢复的第n个信号分量的基带波形;将恢复的基带波形与理想基带波形对比,获取卫星导航信号第n个信号分量的基带波形下降沿抖动的均值和标准差,与预设指标值对比,评估基带波形时域性能。
中国发明专利CN201980007696.1一种基于扩频码时移位置调制的卫星通信导航信号生成方法,包括:生成导频信号分量扩频码和数据信号分量扩频码(110);根据导频信号分量扩频码,调制生成导频分量扩频调制信号Spilot(t)(120);根据数据信号分量扩频码,采用扩频码时移为止调制的方式对二进制的电文或数据信息进行调制,生成数据分量扩频调制信号Sdata(t)(130);采用中心频率为fc、相位不相同的载波,将导频分量扩频调制信号Spilot(t)和数据分量扩频调制信号Sdata(t)调制到射频,得到两路射频分量信号,再将两路射频信号分量叠加,得到射频调制信号SRF(t)(140)。
GPS卫星调整信号分量其中一种模式是选取中心点,以中心点对卫星的仰角变化规划信号能量分配的时机。确定GPS系统规划信号分量调整所选择的中心点,可增进对GPS系统信号调整时机、任务规划能力的了解。
发明内容
本发明需要解决的技术问题是提供一种基于全球网格的GPS信号分量调整中心点确定方法。
为解决上述技术问题,本发明包括一种基于全球网格的GPS信号分量调整中心点确定方法,包括以下步骤:
S1:GPS卫星信号载噪比操作时刻确定,通过分布在全球各地的跟踪站确定GPS卫星载噪比增高和降低的操作时刻;
S2:载噪比操作时刻的卫星位置计算,根据星历计算卫星载噪比操作时刻的卫星位置;
S3:计算网格点对操作卫星的仰角,将全球划分为若干网格,每个网格的经度、纬度分别记为i°、j°,求出位于(i°,j°)的网格点对操作卫星的仰角α(i,j,k),
其中,i=-180,-179,…0,…179,j=-90,-89,…0,…89,k=1,2,…,N,N为信号分量调整操作的卫星编号;
S4:临界角设置,设置不同的临界角βm(m=1,2,…,15),临界角从0°开始,每次步进1°直至15°;
S5:操作卫星数量统计,统计全球不同网格点在不同临界角观测到的操作卫星数;
S6:找出观测到操作卫星数最多的网格点(iM°,jM°)和临界角βM;
S7:仰角极值差判断,如果对所有的k=1,2,…,N,α(iM,jM,k)极差值小于阈值,则认为只有一个中心点,否则,剔除已找到中心点的卫星,再次重复步骤S3-S7继续寻找其他中心点;
S8:确定更加精确的经纬度和临界角,找到中心点后,以中心点为原点,东西、南北拓展范围选择10-30°,划分网格精度选择0.1-0.3°,重复步骤S3~S7一次,以找到更加精确的经纬度和临界角。
优选的,步骤S2中所述星历采用卫星广播星历或者精密星历。
优选的,步骤S7中所述阈值不大于3°。
优选的,所述步骤S8中,以中心点为原点,东西、南北拓展范围选择20°,以0.2°间隔划分网格。
优选的,所述步骤S3中,以1°间隔将全球划分为360×180个网格,高程设置为0。
采用上述方法后,本发明具有以下优点:
(1)本发明首次涉及GPS系统信号分量调整所选中心点的确定方法,利用公开的数据,通过简单的统计计算即可确定出中心点;
(2)研究结果有助于增进对GPS系统信号调整时机、任务规划能力的了解。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明中心点确定方法示意图;
图2为本发明中心点确定方法流程图。
具体实施方式
图1为中心点确定方法的示意图,假设中心点对1号星和2号星仰角为3°时调整卫星信号分量,此时两颗卫星的位置如图1所示,曲线C1上任意一点对1号星仰角均为3°,曲线C2上任意一点对2号星仰角均为3°,两条曲线的交点为A点和B点,这两个点均可能是中心点。实际应用时,由于被操作的卫星数量较多,每一颗卫星都会在地球表面形成一条曲线,所有曲线将交汇在若干个交点,交点最为集中的点即可认为是中心点。
图2为基于全球网格的GPS信号分量调整中心点确定方法的流程图,图1所示的是在地球上画出一个个圆再找交点,图2所示的则是在地球上划分不同的网格、再找出全球观测操作卫星数最大的位置,二者实质是一致的。以载噪比增高所对应的中心点为例,结合附图对本方法各步骤说明如下:
步骤S1,分布在全球各地的跟踪站能够观测到不同卫星信号的载噪比,通过载噪比的瞬时变化可以确定卫星信号载噪比增高(以下记作“操作”)的时刻,由多个跟踪站数据确定出一颗卫星的操作时刻,按此方法统计出所有卫星的操作时刻。
步骤S2,在IGS网站可以下载广播星历或精密星历,利用星历能够计算出步骤1中操作时刻的卫星位置。利用广播星历计算GPS卫星位置见《GPS测量与数据处理(第三版)》75-78页,利用精密星历计算GPS卫星位置见《GPS测量与数据处理(第三版)》78页。
步骤S3,以1°间隔将全球划分为360×180个网格,高程设置为0,每个网格点的经度、纬度分别记为i°、j°,其中,i=-180,-179,…0,…179,j=-90,-89,…0,…89,依次将操作卫星编号为1,2,…,N。由于经过仿真分析,几公里的高程变化对于中心点确定的精度影响不大,为简化方法步骤,提高中心点确定效率,因此将高程设置为0。地面观测点对卫星仰角计算方法见《GPS原理与接收机设计(第二版)》48-49页。
步骤S4,设置临界角βm,临界角从0°开始,每次步进1°直至15°。临界角最大值设置为15°,主要是考虑到,从实用性考虑,转换时机一般设置为比较小的临界角,实际应用中15°可根据实际情况有所变更。信号分量调整所选择的临界角并不唯一,0~15°只是覆盖了临界角度可能的范围;例如,如果临界角度选择的是3°,那么,利用本方法时只要搜索的最大角度大于3°即可;再例如,如果临界角选择的是18°,则利用本方法0~15°范围去搜索,就得不到正确的临界角,此时就需要增大搜索范围。因此,15°只是以往应用本方法时得出的一个经验值。
步骤S5,当步骤3中的α(i,j,k)与步骤4中的βm差值在一定阈值th以内时(如th=1°),则认为(i°,j°)网格点观察到卫星k在临界角βm时进行了1次信号分量调整操作,记作观测到操作卫星数1颗,阈值th设置为1°主要是考虑到卫星执行信号分量调整指令可能出现延迟,导致不同卫星载噪比增高对应的仰角有细微差异,阈值th也并不仅限于1°。统计出(i°,j°)网格点在临界角βm时观察到的操作卫星数T(m,i,j)。利用问号表达式描述T(m,i,j)计算方法,从1到N遍历每颗卫星,如果α(i,j,k)与βm差值小于阈值th,则T(m,i,j)加1,否则T(m,i,j)不增加计数。
步骤S6,对于每个临界角βm(βm=0°,1°,…,15°),分别求出步骤5中T(m,i,j)的最大值T0(m)
使T0(m)达到最大的临界角记作βM,即βM满足
相应的经纬度分别为iM、jM,此时(iM°,jM°)网格点在临界角βM时观察到的操作卫星数为T(M,iM,jM)。如果有多个集中在同一区域的多个点,可对经纬度取中值,例如,记能使观测操作卫星数最大的邻近的网格点数量为p,记q1,q2,…qp为前述p个网格点的经度值,记w1,w2,…wp为前述p个网格点的纬度值,则中心点经度值为
中心点纬度值为
步骤S7,如果对所有的k=1,2,…,N,α(iM,jM,k)极差值小于3°,则认为只有一个中心点,中心点坐标为(iM°,jM°),临界角为βM,否则,剔除已找到中心点的卫星,再次重复步骤3~步骤7继续寻找其他中心点。极差值阈值设为3°的原因是规划信号分量调整时对同一中心点可能设置多个数值相近的临界角,同时考虑到卫星执行指令可能存在延迟,极差值在一定范围内即可认为中心点只有一个,阈值并不仅限于3°。
步骤S8,找到中心点后,以中心点为原点,向东、西、南、北各拓展20°,以0.2°间隔划分网格,重复步骤3~步骤7一次,找到更加精确的经纬度和临界角。通过更细小的网格点确定中心点,能够进一步确定出中心点精确的位置和临界角,网格范围、大小可结合实际进行调整。步骤S3~S7只需要重复1次,图2中“重新划分网格,精细寻找中心点和临界角”即是代表此步骤。重复一次步骤S3~S7的目的是使寻找到的中心点经纬度精度更高。实际上,进行到S7步骤,找到的中心点经纬度精度为1°,进行到S8(即重复一次S3~S7后),找到的中心点经纬度精度为0.2°
虽然以上描述了本发明的具体实施方式,但是本领域熟练技术人员应当理解,这些仅是举例说明,可以对本实施方式作出多种变更或修改,而不背离本发明的原理和实质,本发明的保护范围仅由所附权利要求书限定。
Claims (5)
1.一种基于全球网格的GPS信号分量调整中心点确定方法,其特征在于,包括以下步骤:
S1:GPS卫星信号载噪比操作时刻确定,通过分布在全球各地的跟踪站确定GPS卫星载噪比增高和降低的操作时刻;
S2:载噪比操作时刻的卫星位置计算,根据星历计算卫星载噪比操作时刻的卫星位置;
S3:计算网格点对操作卫星的仰角,将全球划分为若干网格,每个网格的经度、纬度分别记为i°、j°,求出位于(i°,j°)的网格点对操作卫星的仰角α(i,j,k),
其中,i=-180,-179,...0,...179,j=-90,-89,...0,...89,k=1,2,…,N,N为信号分量调整操作的卫星编号;
S4:临界角设置,设置不同的临界角βm(m=1,2,…,15),临界角从0°开始,每次步进1°直至15°;
S5:操作卫星数量统计,统计全球不同网格点在不同临界角观测到的操作卫星数;
S6:找出观测到操作卫星数最多的网格点(iM°,jM°)和临界角βM;
S7:仰角极值差判断,如果对所有的k=1,2,…,N,α(iM,jM,k)极差值小于阈值,则认为只有一个中心点,否则,剔除已找到中心点的卫星,再次重复步骤S3-S7继续寻找其他中心点;
S8:确定更加精确的经纬度和临界角,找到中心点后,以中心点为原点,东西、南北拓展范围选择10-30°,划分网格精度选择0.1-0.3°,重复步骤S3~S7一次,以找到更加精确的经纬度和临界角。
2.根据权利要求1所述的基于全球网格的GPS信号分量调整中心点确定方法,其特征在于:步骤S2中所述星历采用卫星广播星历或者精密星历。
3.根据权利要求1所述的基于全球网格的GPS信号分量调整中心点确定方法,其特征在于:步骤S7中所述阈值不大于3°。
4.根据权利要求1所述的基于全球网格的GPS信号分量调整中心点确定方法,其特征在于:所述步骤S8中,以中心点为原点,东西、南北拓展范围选择20°,以0.2°间隔划分网格。
5.根据权利要求1所述的基于全球网格的GPS信号分量调整中心点确定方法,其特征在于:所述步骤S3中,以1°间隔将全球划分为360×180个网格,高程设置为0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211295087.0A CN115685266B (zh) | 2022-10-21 | 2022-10-21 | 基于全球网格的gps信号分量调整中心点确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211295087.0A CN115685266B (zh) | 2022-10-21 | 2022-10-21 | 基于全球网格的gps信号分量调整中心点确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115685266A CN115685266A (zh) | 2023-02-03 |
CN115685266B true CN115685266B (zh) | 2023-07-11 |
Family
ID=85067479
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211295087.0A Active CN115685266B (zh) | 2022-10-21 | 2022-10-21 | 基于全球网格的gps信号分量调整中心点确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115685266B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108828626A (zh) * | 2018-07-02 | 2018-11-16 | 中国人民解放军战略支援部队信息工程大学 | 基于实时网格的网络rtk电离层延迟内插方法及系统 |
CN112014863A (zh) * | 2020-08-21 | 2020-12-01 | 中国人民解放军61081部队 | 基于多地面站数据的gnss功率变化时刻自动分析方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6411892B1 (en) * | 2000-07-13 | 2002-06-25 | Global Locate, Inc. | Method and apparatus for locating mobile receivers using a wide area reference network for propagating ephemeris |
US6670916B2 (en) * | 2002-02-19 | 2003-12-30 | Seiko Epson Corporation | Method and system for position calculation from calculated time |
CN112444834B (zh) * | 2019-08-29 | 2024-07-30 | 华为技术有限公司 | 一种定位方法及电子设备 |
-
2022
- 2022-10-21 CN CN202211295087.0A patent/CN115685266B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108828626A (zh) * | 2018-07-02 | 2018-11-16 | 中国人民解放军战略支援部队信息工程大学 | 基于实时网格的网络rtk电离层延迟内插方法及系统 |
CN112014863A (zh) * | 2020-08-21 | 2020-12-01 | 中国人民解放军61081部队 | 基于多地面站数据的gnss功率变化时刻自动分析方法 |
Non-Patent Citations (3)
Title |
---|
"Research on Secondary Analysis Method of Synchronous Satellite Monitoring Data of Power Grid Wildfire";Yuelun Di等;《2020 IEEE International Conference on Information Technology,Big Data and Artificial Intelligence 》;706-710 * |
"一种GPS定位实时选星改进算法研究";李士明等;《计算机仿真》;第26卷(第7期);65-68 * |
"近地全球重访星座轨道快速优化设计研究";马原野;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》(第8期);C031-230 * |
Also Published As
Publication number | Publication date |
---|---|
CN115685266A (zh) | 2023-02-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kassas et al. | I hear, therefore I know where I am: Compensating for GNSS limitations with cellular signals | |
Zhang et al. | A fast satellite selection algorithm: beyond four satellites | |
de Camargo et al. | Application of ionospheric corrections in the equatorial region for L1 GPS users | |
US9395443B2 (en) | Adaptive method for estimating the electron content of the ionosphere | |
CN103064092B (zh) | 一种导航卫星的选择方法 | |
CN110531315B (zh) | 一种基于信号强度变化率的卫星干扰源直接定位方法及定位装置 | |
Elango et al. | Design of complete software GPS signal simulator with low complexity and precise multipath channel model | |
Tobie et al. | Hybrid navigation filters performances between GPS, Galileo and 5G TOA measurements in multipath environment | |
CN113281702B (zh) | 协同短波多站角度与卫星时频的超视距目标直接定位方法 | |
CN115685266B (zh) | 基于全球网格的gps信号分量调整中心点确定方法 | |
Raghu et al. | Tracking of IRNSS, GPS and hybrid satellites by using IRNSS receiver in STK simulation | |
Qi et al. | Robust unambiguous ranging technique for BDS-3 B1 wideband composite signals | |
CN110208741B (zh) | 一种基于多圆阵测相的超视距单目标直接定位方法 | |
Liu et al. | Analysis on ionospheric delay variogram modeling in China | |
CN113970762B (zh) | 一种多级干扰源定位方法及系统 | |
CN112152651B (zh) | 5g系统面向gnss接收机的干扰源定位方法、存储介质及设备 | |
Brevick et al. | Evaluating Navigation Augmentation with LTE in the Urban Signal Environment | |
Namassivaya et al. | Modelling of fpga-particle swarm optimized gnss receiver for satellite applications | |
Wang et al. | Study on multipath model of BDS/GPS signal in urban canyon | |
Abdallah et al. | Multipath mitigation of 5g signals via reinforcement learning for navigation in urban environments | |
Liu et al. | A main peak extraction method for high-order BOC signals | |
Xu et al. | Satellite-induced multipath analysis on the cause of BeiDou code pseudorange bias | |
Rao et al. | Hybrid technique for GPS receiver position applications | |
PILEGGI | 5G positioning: preliminary analysis of early data sets | |
Zhang et al. | Receiver differential code bias estimation under disturbed ionosphere status using linear planar model based minimum standard deviation searching method with bias detection |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |