CN112231873B - 一种meo星座轨位冲突解决方法 - Google Patents
一种meo星座轨位冲突解决方法 Download PDFInfo
- Publication number
- CN112231873B CN112231873B CN202011101955.8A CN202011101955A CN112231873B CN 112231873 B CN112231873 B CN 112231873B CN 202011101955 A CN202011101955 A CN 202011101955A CN 112231873 B CN112231873 B CN 112231873B
- Authority
- CN
- China
- Prior art keywords
- track
- orbit
- intersection point
- track surface
- satellite
- 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
- 238000000034 method Methods 0.000 title claims abstract description 26
- 230000000630 rising effect Effects 0.000 claims description 12
- 101100221616 Halobacterium salinarum (strain ATCC 29341 / DSM 671 / R1) cosB gene Proteins 0.000 claims 2
- 235000002492 Rungia klossii Nutrition 0.000 claims 2
- 244000117054 Rungia klossii Species 0.000 claims 2
- 230000001174 ascending effect Effects 0.000 description 4
- 229920000535 Tan II Polymers 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Networks & Wireless Communication (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供一种MEO星座轨位冲突解决方法包括:步骤一:确定MEO星座各星标称轨道及其相位保持范围,步骤二:确定重叠轨位所在轨道面与其它两个轨道面的交点位置,步骤三:确定不同轨道面卫星同时过轨道面交点,步骤四:确定不同轨位偏移量与其它两个轨道面卫星的相位差,步骤五:确定新轨位及相位的保持范围,本发明提出基于相位的新轨位确定方法,解决MEO星座轨位重叠带来的频率干扰和碰撞风险,适用解决MEO星座轨位冲突问题,可为后续备份卫星轨位设计提供参考。
Description
技术领域
本发明属于在轨航天器测控管理领域,具体指一种MEO星座轨位冲突解决方法,适用于解决MEO星座轨位冲突问题。
背景技术
中地球轨道(Medium Earth Orbit,MEO)主要是指距离地球表面2000至 20000公里的地球轨道。北斗导航MEO星座共设计24颗卫星,卫星轨道高度相同,平均分布在三个轨道面,组成Walk星座构型中的δ星座。相邻平面的对应序号的卫星之间的相位差15度,每条轨道的升交点以间隔为120度均匀分布,每条轨道上的卫星按等间隔45度均匀分布。
在北斗二号向北斗三号平稳过渡及北斗三号全球导航星座组网建设过程中,由于北斗二号与北斗三号MEO卫星共用MEO星座轨位,造成北斗二号与北斗三号MEO星座卫星轨位重叠,需要重新设计MEO卫星轨位,确定新轨位相位保持要求,以避免与北斗三号MEO卫星出现频率干扰和碰撞风险。
从上面分析可以可知,有必要提供MEO星座轨位冲突解决方法来支持北斗导航星座建设。
发明内容
本发明的目的提供一种MEO星座轨位冲突解决方法,针对MEO星座卫星轨位重叠问题,在不改变轨道面的前提下,提出基于相位的新轨位确定方法,解决MEO星座轨位重叠带来的频率干扰和碰撞风险。
为实现上述目的,本发明提供如下技术方法:
步骤一:确定MEO星座各星标称轨道及相位保持要求
将出现轨位重叠轨道面定义为第一轨道面,该重叠轨位定义为该轨道面第一相位,且为整个星座的基准星,其余卫星与其相位差为Δui(i=0,1,…,n),在标称轨位的保持范围为±ε度,碰撞规避要求的相位差σ。
将整个星座中的卫星按照顺序进行编号,对于第i轨道面的第j颗卫星编号为MEOij,其与基准星之间标称相对相位Δμij为:
Δμij=mod((i-1)*15+(j-1)*45,360)
步骤二:计算重叠轨位所在轨道面与其它两个轨道面交点位置
第一轨道面与第二轨道面的交点,因两个轨道面交点相差幅角为180度,只需计算一个交点。
第一轨道面倾角矢量:
第二轨道面倾角矢量:
两轨道面倾角矢量差为
计算两轨道面交点中靠近第一轨道面升交点的交点的赤经:
l=a tan 2(Δi21x,Δi21y)
该交点在第一轨道面距离升交点的角距u12为:
cos B=cos l sin i
cos u=tan(90-B)tan i
u12=360-u
该交点在第二轨道面距离升交点的角距u21为:
u21=180+u
第一轨道面与第三轨道面的交点
第一轨道面倾角矢量:
第三轨道面倾角矢量:
两轨道面倾角矢量差为
计算两轨道面交点中靠近第一轨道面升交点的交点的赤经:
l=a tan 2(Δi21x,Δi21y)
该交点在第一轨道面距离升交点的角距u13为:
cos B=cos l sin i
cos u=tan(90-B)tan i
u13=u
该交点在第三轨道面距离升交点的角距u31为:
u31=180-u
步骤三:确定不同轨道面卫星同时过轨道面交点时具备的条件
第一个轨道面与第二轨道面卫星同时过交点的两星相位差满足条件
Δu12=u21-u12=2u-180
第一个轨道面与第三轨道面卫星同时过交点的两星相位差满足条件
Δu13=u31-u13=180-2u
步骤四:计算不同轨位偏移量与其它两个轨道面卫星的相位差
综合考虑与轨位重叠卫星及与其共面卫星的相位隔离问题,轨位偏离的范围为(ε+σ,45-ε-σ)。在此范围内,以ε+σ作为初始偏差,按照ΔP作为步长,分别计算不同轨位偏移量情况下其它轨道面卫星与其相位差
Δμijk=mod((i-1)*15+(j-1)*45,360)+ε+kΔP,k=1,2,…,n
其中n=int((45-2ε)/ΔP)
步骤五:确定新轨位及相位保持范围
根据上述计算的n组相位差Δμijk,分别计算统计与步骤三中相位差满足条件差值的最小值。
Δumin,k=min(Δμijk-Δu),k=1,2,…,n
其中若为第二轨道面卫星,Δu=Δu12;若为第三轨道面卫星,Δu=Δu13。
统计相位差最小值大于ε+σ的区间,以区间范围最大作为新轨位的保持范围,该区间的中间值作为新轨位。
与现有技术相比,本发明的有益效果是:本发明提出基于相位的新轨位确定方法,解决MEO星座轨位重叠带来的频率干扰和碰撞风险,适用解决MEO星座轨位冲突问题,可为后续备份卫星轨位设计提供参考。
附图说明
图1是本发明相位差差值最小值的分布图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方法进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例:一种MEO星座轨位冲突解决方法,包括如下步骤:
步骤一:确定MEO星座各星标称轨道及相位保持要求
北斗MEO轨位重叠轨道面定义为第一轨道面,该重叠轨位定义为该轨道面第一相位,且为整个星座的基准星。星座卫星标称轨道如下:
第一轨道面倾角为55度,升交点赤经为0度,其7颗卫星的标称相位分别为45,90,135,180,225,270,315
第二轨道面倾角为55度,升交点赤经为120度,标称相位分别为15,60, 105,150,195,240,285,330
第三轨道面倾角为55度,升交点赤经为240度,标称相位分别为30,75, 120,165,210,255,300,345
星座卫星相位保持范围在标称值的±5度范围内,碰撞规避要求的相位差 0.5度。
步骤二:计算重叠轨位所在轨道面与其它两个轨道面交点位置
第一轨道面与第二轨道面的交点
根据技术方案中步骤二的公式计算,u为45.18度,该交点在第一轨道面距离升交点的角距u12为314.82度,该交点在第二轨道面距离升交点的角距u21为 225.18度。
第一轨道面与第三轨道面的交点
根据技术方案中步骤二的公式计算,u为45.18度,该交点在第一轨道面距离升交点的角距u13为45.18度,该交点在第三轨道面距离升交点的角距u31为 134.82度。
步骤三:确定不同轨道面卫星同时过轨道面交点时具备的条件
第一个轨道面与第二轨道面卫星同时过交点的两星相位差Δu12等于-89.64 度。
第一个轨道面与第三轨道面卫星同时过交点的两星相位差Δu13等于89.64 度。
步骤四:计算不同轨位偏移量与其它两个轨道面卫星的相位差
轨位偏离的范围为(6.5,38.5),以6.5度作为初始偏差,按照0.5度作为步长,分别计算不同轨位偏移量情况下其它轨道面卫星与其相位差,具体见下表:
表1第二轨道面卫星的相位差
表2第三轨道面卫星的相位差
/>
/>
步骤五:确定新轨位及相位保持范围
根据上述计算的69组相位差,分别计算统计每一组相位差与步骤二的两星相位差的差值 的最小值,具体见附图1。
相位差最小值大于5.5度的区间见下表,因此,以偏离基准星22.5度作为新轨位,保持范围为±2度。
虽然在上文中已经参考实施例对本发明进行了描述,然而在不脱离本发明的范围的情况下,可以对其进行各种改进并且可以用等效物替换其中的部件。尤其是,只要不存在结构冲突,本发明所披露的实施例中的各项特征均可通过任意方式相互结合起来使用,在本说明书中未对这些组合的情况进行穷举性的描述仅仅是出于省略篇幅和节约资源的考虑。因此,本发明并不局限于文中公开的特定实施例,而是包括落入权利要求的范围内的所有技术方案。
Claims (2)
1.一种MEO星座轨位冲突解决方法,其特征在于,包括如下步骤:
步骤一:确定MEO星座各星标称轨道及其相位保持范围;
将出现轨位重叠轨道面定义为第一轨道面,该重叠轨位定义为该轨道面第一相位,且为整个星座的基准星,其余卫星与其相位差为Δui(i=0,1,…,n),在标称轨位的保持范围为±ε度,碰撞规避要求的相位差σ;
将整个星座中的卫星按照顺序进行编号,对于第i轨道面的第j颗卫星编号为MEOij,其与基准星之间标称相对相位Δμij为:
Δμij=mod((i-1)*15+(j-1)*45,360);
步骤二:确定重叠轨位所在轨道面与其它两个轨道面的交点位置;
第一轨道面倾角矢量:
第二轨道面倾角矢量:
两轨道面倾角矢量差为
两轨道面交点中靠近第一轨道面升交点的交点的赤经:
l=atan2(Δi21x,Δi21y)
该交点在第一轨道面距离升交点的角距u12为:
cosB=coslsini
cosu=tan(90-B)tani
u12=360-u
该交点在第二轨道面距离升交点的角距u21为:
u21=180+u
第一轨道面与第三轨道面的交点
第一轨道面倾角矢量:
第三轨道面倾角矢量:
两轨道面倾角矢量差为
两轨道面交点中靠近第一轨道面升交点的交点的赤经:
l=atan2(Δi21x,Δi21y)
该交点在第一轨道面距离升交点的角距u13为:
cosB=coslsini
cosu=tan(90-B)tani
u13=u
该交点在第三轨道面距离升交点的角距u31为:
u31=180-u;
步骤三:确定不同轨道面卫星同时过轨道面交点
第一个轨道面与第二轨道面卫星同时过交点的两星相位差满足条件:
Δu12=u21-u12=2u-180
第一个轨道面与第三轨道面卫星同时过交点的两星相位差满足条件:
Δu13=u31-u13=180-2u;
步骤四:确定不同轨位偏移量与其它两个轨道面卫星的相位差
轨位偏离的范围为(ε+σ,45-ε-σ),在此范围内,以ε+σ作为初始偏差,按照ΔP作为步长,分别确定不同轨位偏移量情况下其它轨道面卫星与其相位差
Δμijk=mod((i-1)*15+(j-1)*45,360)+ε+kΔP,k=1,2,...,n
其中n=int((45-2ε)/ΔP);
步骤五:确定新轨位及相位的保持范围;
根据上述步骤确定的n组相位差Δμijk,分别确定统计与步骤三中相位差满足条件差值的最小值,
Δumin,k=min(Δμijk-Δu),k=1,2,…,n
其中若为第二轨道面卫星,Δu=Δu12;若为第三轨道面卫星,Δu=Δu13;
统计相位差最小值大于ε+σ的区间,以区间范围最大作为新轨位的保持范围,该区间的中间值作为新轨位。
2.根据权利要求1所述的一种MEO星座轨位冲突解决方法,其特征在于,所述步骤二,重叠轨位所在轨道面与其它两个轨道面的交点位置,所述第一轨道面与第二轨道面的交点,因两个轨道面交点相差幅角为180度,只需确定一个交点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011101955.8A CN112231873B (zh) | 2020-10-15 | 2020-10-15 | 一种meo星座轨位冲突解决方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011101955.8A CN112231873B (zh) | 2020-10-15 | 2020-10-15 | 一种meo星座轨位冲突解决方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112231873A CN112231873A (zh) | 2021-01-15 |
CN112231873B true CN112231873B (zh) | 2024-03-26 |
Family
ID=74112985
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011101955.8A Active CN112231873B (zh) | 2020-10-15 | 2020-10-15 | 一种meo星座轨位冲突解决方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112231873B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108055069A (zh) * | 2017-12-11 | 2018-05-18 | 中国人民解放军战略支援部队航天工程大学 | 低轨通信和导航增强混合星座维持控制边界计算及控制方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7744037B2 (en) * | 2007-04-11 | 2010-06-29 | National Defense University | Calculating method for deducing possibilities of all possible scenarios of satellite members in low earth satellite constellation |
-
2020
- 2020-10-15 CN CN202011101955.8A patent/CN112231873B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108055069A (zh) * | 2017-12-11 | 2018-05-18 | 中国人民解放军战略支援部队航天工程大学 | 低轨通信和导航增强混合星座维持控制边界计算及控制方法 |
Non-Patent Citations (2)
Title |
---|
建链性能最优的中高轨卫星分布规律研究;胡敏;范丽;杨雪榕;;导航定位学报;20150320(第01期);全文 * |
面向GEO卫星频轨协调风险的定量化评估;周鑫林;巩应奎;韩朝晖;;计算机仿真;20180215(第02期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112231873A (zh) | 2021-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108562295B (zh) | 一种基于同步卫星二体模型的三站时差定轨方法 | |
CN108761507B (zh) | 基于短弧定轨和预报的导航卫星轨道快速恢复方法 | |
CN104076819B (zh) | 一种圆参考轨道下卫星有界伴飞的边界控制方法 | |
JP6271043B2 (ja) | 軌道制御装置および衛星 | |
CN105511483B (zh) | 鸟巢式星座及其设计方法 | |
CN111102982B (zh) | 高轨目标的抵近方法 | |
CN112231873B (zh) | 一种meo星座轨位冲突解决方法 | |
CN107506893B (zh) | 一种太阳同步轨道航天器安全管理方法 | |
CN110789739A (zh) | 一种j2摄动下长时间轨道交会最优速度增量快速估计方法 | |
RU2559371C2 (ru) | Способ автономной коллокации на геостационарной орбите | |
CN110765504A (zh) | 航天器环月轨道交会对接的轨道设计方法 | |
CN110068846B (zh) | 一种基于星载gnss接收机在星上自主确定轨道平根数的方法 | |
CN111806729B (zh) | 考虑拱线旋转的非冻结轨道多星定位编队设计方法 | |
CN113311854A (zh) | 一种月球采样返回任务中的定点着陆轨道设计方法 | |
Baranov et al. | Ballistic aspects of large-size space debris flyby at low Earth near-circular orbits | |
CN111541477B (zh) | 低轨星座系统内部频率干扰抑制方法及装置 | |
US6317661B1 (en) | Argument of perigee correction with longitude control for inclined, eccentric, geosynchronous satellites | |
CN110275165B (zh) | 一种编队geo sar的等效相位中心及积累时间优化方法 | |
CN1227026A (zh) | 通过附加采样的推式扫帚扫描获取图象的方法 | |
CN115933743A (zh) | 一种针对目标窗口交会的轨迹规划方法 | |
Razoumny | New research methodology for Earth periodic coverage and regularities in parametric localization of optimal low-Earth-orbit satellite constellations | |
RU2703696C1 (ru) | Способ автономной коллокации на околостационарной орбите | |
CN107479565B (zh) | 基于椭圆轨道像移补偿计算方法 | |
CN106005482B (zh) | 一种适用于导航倾斜轨道卫星的零偏持续天数确定方法 | |
CN114684389A (zh) | 考虑再入约束的月地转移窗口及精确转移轨道确定方法 |
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 |