CN106156417A - 等时间间隔快速重访卫星星座构型优化设计方法 - Google Patents
等时间间隔快速重访卫星星座构型优化设计方法 Download PDFInfo
- Publication number
- CN106156417A CN106156417A CN201610514705.4A CN201610514705A CN106156417A CN 106156417 A CN106156417 A CN 106156417A CN 201610514705 A CN201610514705 A CN 201610514705A CN 106156417 A CN106156417 A CN 106156417A
- Authority
- CN
- China
- Prior art keywords
- satellite
- target
- constellation
- orbit
- optimization
- 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
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 74
- 238000013461 design Methods 0.000 title claims abstract description 67
- 238000000034 method Methods 0.000 claims abstract description 30
- 230000002068 genetic effect Effects 0.000 claims abstract description 18
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 14
- 230000003287 optical effect Effects 0.000 claims abstract description 4
- 238000005265 energy consumption Methods 0.000 claims description 16
- 238000003384 imaging method Methods 0.000 claims description 16
- 238000004088 simulation Methods 0.000 claims description 10
- 238000012634 optical imaging Methods 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 3
- 230000007774 longterm Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 101100063069 Caenorhabditis elegans deg-1 gene Proteins 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 229910052741 iridium Inorganic materials 0.000 description 1
- GKOZUEZYRPOHIO-UHFFFAOYSA-N iridium atom Chemical compound [Ir] GKOZUEZYRPOHIO-UHFFFAOYSA-N 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification 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/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Radio Relay Systems (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供一种等时间间隔快速重访卫星星座构型优化设计方法,该方法以实现对任务目标的等时间间隔快速重访为设计目标,采用低倾角回归Walker星座为基本构型,将卫星数目、轨道高度、重访时间间隔、以及单颗卫星的覆盖特性、机动能力以及光学成像性能等指标作为优化约束条件,基于多岛遗传算法以较少的卫星实现了对指定目标的等时间间隔快速重访。采用该方法设计的卫星星座具有对目标区域实现等时间间隔快速重访的优点,优化得到的卫星星座能够符合任务需求和各项设计约束,可作为快速响应卫星星座部署的参考方案。
Description
技术领域
本发明涉及航天技术领域,具体的涉及一种等时间间隔快速重访卫星星座构型优化设计方法。
背景技术
卫星星座主要应用于通信、导航、对地观测空间任务中,利用卫星星座能够实现对目标的连续和多层覆盖。由于单颗卫星只能观察地球的一个球冠区域,即仅能实现对地瞬时覆盖,单颗卫星的覆盖面积与卫星轨道高度相关。对于低轨轨道的卫星而言,如果仅使用单颗卫星,此时该卫星的地球覆盖面积远远小于地球整体的表面积。在低轨道,单卫星的情况下,可借助卫星星座来实现对地球某个区域或纬度带的连续、多层或等时间间隔的覆盖。
星座构型设计的任务主要是采用优化设计方法对星座构型参数进行合理设计,待设计的构型参数通常指卫星数目、卫星轨道平面、平面内的卫星数目、卫星相位角、卫星轨道高度等参数。现有优化方法包括以下步骤:首先将星座的设计准则转换为目标函数之一或者是约束条件之一;然后建立星座构型和星座性能的优化设计模型;最后根据设计任务的优化目标,选择合适的优化算法来完成卫星星座构型的优化设计。
为实现对目标的连续覆盖,现有的技术方案通常需要使用大量的低轨卫星组成星座。如Iridium星座,为了实现对地球表面的全覆盖,需要66颗卫星分别部署在高度677km、轨道倾角为86.4°的6条轨道上。从卫星制造、发射成本以及运营管理的角度来说,这都是不可接受的。
发明内容
针对上述不足,本发明提供了一种等时间间隔快速重访卫星星座构型优化设计方法。
本发明提供一种等时间间隔快速重访卫星星座构型优化设计方法,包括以下步骤:
步骤S100:构建星座构型模型、回归轨道模型和性能优化模型;
步骤S200:设置星座构型优化设计参数、星座构型设计输入约束条件、星座构型设计输出约束条件和优化目标函数;
步骤S300:将各项参数代入公式(5)~(15),通过多岛遗传算法对卫星星座构型进行优化设计,得到能实现优化目标的等时间间隔快速重访的卫星星座构型设计结果;
其中步骤S100包括以下步骤:
步骤S110:构建星座构型模型:
参考Walker星座构型,建立对地观测卫星星座构型模型:N颗卫星均匀分布在P个轨道面上,每个轨道面上设置S颗卫星,且所形成的星座中所有卫星均具有相同的轨道高度、轨道偏心率、轨道倾角和近地点幅角;所设置的星座中各卫星之间的各轨道面之间的升交点赤经差,用△Ω表示,所设置的星座中各卫星之间的各轨道面的首颗卫星之间的平近点角差,用△M表示;
假设第1轨道面第1颗卫星的轨道参数为(a0、e0、i0、Ω0、ω0、M0),对于圆轨道,卫星轨道偏心率e0和近地点幅角ω0取值:
e0=0,ω0=0 (1)
则所构建星座中的第j个轨道面中第k颗卫星的轨道参数按公式(2)计算得到:
其中j=[1,2,…,P],k=[1,2,…,S],P为轨道面数目,S为每个轨道面的卫星数目;
所构建星座中的卫星总数N满足公式(3):
N=PS (3)
步骤S120:构建回归轨道模型:
所构建星座中卫星需满足公式(4)所示的回归特性方程:
其中N*是回归圈数,D*是回归周期,ωe是地球自转的角速度,是卫星轨道的升交点赤经的长期变化速率,TΩ为轨道交点周期;
步骤S130:性能优化模型:
低轨光学卫星需考虑的性能指标包括卫星成本、发射费用、入轨能耗和卫星对目标的覆盖性能、卫星对目标的成像性能进行优化设计;
(1)卫星成本
影响卫星星座系统成本的主要因素是卫星轨道数目N,建立卫星成本优化目标模型为:
fcost=min(N) (5)
(2)发射费用
建立卫星发射费用成本模型优化目标模型为:
flaunch=λ1min(P)+λ2min(h) (6)
其中,λ1和λ2分别是轨道面数目和轨道高度的加权系数,P为轨道面数目和h为卫星轨道高度;
(3)入轨能耗
影响入轨能耗的主要因素是卫星轨道高度;
轨道越低,则发射卫星所需的能耗就越低;
建立卫星发射入轨能耗优化目标模型为:
fpropulsion=min(△V) (7)
其中,△V为卫星入轨所需的速度增量:
其中RE为地球半径,μ为地球引力系数,al发射轨道半长轴,a0为卫星运行轨道半长轴;且al为
式中
其中βL为发射角;
(4)对目标的覆盖性能
星座对目标的观测性能指标:每天对目标的观测次数Ctimes、每次平均观测时间Tmean、对目标的最大观测间隔Tmaxgap、对目标的最小观测间隔Tmingap、对目标的平均观测间隔Tavegap、对目标的最大观测间隔的均方差△Tmaxgap、最小观测间隔的均方差△Tmingap,其中,
△Tmaxgap=|Tmaxgap-Tavegap| (12)
△Tmingap=|Tavegap-Tmingap| (13)
建立星座对目标的观测性能优化目标模型为:
fobersvation=λ1max(Ctimes)+λ2max(Tmean)+λ3min(Tmaxgap)+
λ4min(Tmingap)+λ5min(Tavegap)+λ6min(△Tmaxgap)+λ7min(△Tmingap) (14)
其中λ1、λ2、λ3、λ4、λ5、λ6和λ7分别是每天对目标的观测次数、每次平均观测时间、对目标的最大观测间隔、对目标的最小观测间隔、对目标的平均观测间隔和对目标观测间隔的均方差的加权系数值;
(5)卫星对目标成像性能
卫星目标的成像性能包括:空间分辨率L、幅宽E和视场角;
其中空间分辨率和幅宽与轨道高度、成像载荷密切相关;
建立星座卫星对目标成像性能优化目标模型:
fimage=λ1min(L)+λ2max(E) (15)
其中,λ1和λ2分别是空间分辨率和幅宽的加权系数值;
基于以上步骤所得的各个优化设计模型即公式(5)~(15),可得到针对低倾角回归轨道的,对目标位置的性能最优的卫星星座性能的优化设计模型:
fperformance=ξ1fcost+ξ2flaunch+ξ3fpropulsion+ξ4fobersvation+ξ5fimage (16)
其中ξ1、ξ2、ξ3、ξ4和ξ5分别为各性能项的加权系数值。
进一步地,星座构型优化设计参数设置为:目标点经度设置为135°,目标点纬度设置为30°,每天仿真步数86400步,仿真时长为回归轨道的回归天数D*,轨道高度下限250km,轨道高度上限1200km;低轨卫星光学成像载荷视频摄像机的参数设置:CMOS像元大小6μm,镜头焦距2000mm,CMOS分辨率设置为4096×3072,面阵对角线长度为30mm。
进一步地,多岛遗传算法中的优化参数:子代大小设置为5、岛数设置为5、遗传代数设置为200。
相对现有技术,本发明的技术效果:
1、本发明提供等时间间隔快速重访卫星星座构型优化设计方法,利用多岛遗传算法完成等时间间隔快速重访卫星星座构型的优化设计,该方法的仅需利用低轨低倾角星座发射,降低了发射成本,便于工程实现。
2、本发明提供等时间间隔快速重访卫星星座构型优化设计方法,综合考虑了卫星星座的观测性能、星座发射费用、卫星成本等各方面问题,将星座构型设计转变成了多目标、多约束优化问题,使其更符合实际情况。
3、本发明提供等时间间隔快速重访卫星星座构型优化设计方法,所设计的卫星星座构型能够实现对目标区域的等时间间隔快速重访,提高了重访效率。
具体请参考根据本发明的等时间间隔快速重访卫星星座构型优化设计方法提出的各种实施例的如下描述,将使得本发明的上述和其他方面优势显而易见。
附图说明
图1是本发明提供的等时间间隔快速重访卫星星座构型优化设计方法流程示意图;
图2是本发明提供的等时间间隔快速重访卫星星座构型优化设计方法所得卫星星座对目标点的观测时间分布示意图;
图3是本发明提供的等时间间隔快速重访卫星星座构型优化设计方法所得卫星星座目标点的覆盖百分率示意图。
具体实施方式
构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。
本发明以实现对任务目标区域的等时间间隔快速重访为设计目标,综合考虑星座构型、轨道和性能等约束条件,采用多岛遗传算法完成了可实现对特定目标的等时间间隔快速重访卫星星座构型的优化设计。对所设计的等时间间隔快速重访卫星星座的性能分析结果表明,优化得到的卫星星座能够符合任务需求和各项设计约束,且并克服了星座卫星数目大、成本高的缺点。
本发明的技术方案是:以低倾角回归Walker星座为基本构型,以较少的卫星实现对指定目标进行连续的或间隙时间最短、性能最优的观测,要满足这些观测性能时通常需要设置包括卫星数目、轨道高度、重访时间间隔等在内的多个性能指标为主要优化目标;优化约束条件除考虑卫星数目、轨道高度、重访时间间隔等参数以外,还考虑了单颗卫星的覆盖特性、机动能力以及光学成像性能等指标;从而解决对指定目标的光学成像卫星星座的多约束、多目标优化设计问题。
参见图1,本发明提供的等时间间隔快速重访卫星星座构型优化设计方法包括以下步骤:
步骤S100:构建星座构型模型、回归轨道模型和性能优化模型;
步骤S200:设置星座构型优化设计参数、星座构型设计输入约束条件、星座构型设计输出约束条件和优化目标函数;
步骤S300:将各项参数代入公式(5)~(15),通过多岛遗传算法对卫星星座构型进行优化设计,得到能实现优化目标的等时间间隔快速重访的卫星星座构型设计结果;
其中步骤S100包括以下步骤:
步骤S110:构建星座构型模型:
参考Walker星座构型,建立对地观测卫星星座构型模型:N颗卫星均匀分布在P个轨道面上,每个轨道面上设置S颗卫星,且所形成的星座中所有卫星均具有相同的轨道高度、轨道偏心率、轨道倾角和近地点幅角。所设置的星座中各卫星之间的差异主要是不同轨道面之间的升交点赤经差,用△Ω表示,以及不同轨道面的首颗卫星之间的平近点角差,用△M表示。
假设第1轨道面第1颗卫星的轨道参数为(a0、e0、i0、Ω0、ω0、M0),对于圆轨道,卫星轨道偏心率e0和近地点幅角ω0取值如下:
e0=0,ω0=0 (1)
则所构建星座中的第j个轨道面中第k颗卫星的轨道参数按公式(2)计算得到:
其中j=[1,2,…,P],k=[1,2,…,S],P为轨道面数目,S为每个轨道面的卫星数目。
所构建星座中的卫星总数N满足公式(3):
N=PS (3)
步骤S120:构建回归轨道模型:
所构建星座中卫星需满足公式(4)所示的回归特性方程:
其中N*是回归圈数,D*是回归周期,ωe是地球自转的角速度,是卫星轨道的升交点赤经的长期变化速率,TΩ为轨道交点周期。
步骤S130:性能优化模型:
低轨光学卫星需考虑的性能指标包括卫星成本、发射费用、入轨能耗和卫星对目标的覆盖性能、卫星对目标的成像性能进行优化设计。
(1)卫星成本
影响卫星星座系统成本的主要因素是卫星轨道数目N,建立卫星成本优化目标模型为:
fcost=min(N) (5)
(2)发射费用
由于不同轨道面上的卫星无法采用同一颗火箭发射入轨,因此影响发射费用的主要因素是轨道面数目P和卫星轨道高度h。建立卫星发射费用成本模型优化目标模型为:
flaunch=λ1min(P)+λ2min(h) (6)
其中,λ1和λ2分别是轨道面数目和轨道高度的加权系数。
(3)入轨能耗
影响入轨能耗的主要因素是卫星轨道高度。轨道越低,则发射卫星所需的能耗就越低。建立卫星发射入轨能耗优化目标模型为:
fpropulsion=min(△V) (7)
其中,△V为卫星入轨所需的速度增量:
其中RE为地球半径,μ为地球引力系数,al发射轨道半长轴,a0为卫星运行轨道半长轴;且al为:
式中
其中βL为发射角。
(4)对目标的覆盖性能
星座对目标的观测性能主要包括如下指标:每天对目标的观测次数Ctimes、每次平均观测时间Tmean、对目标的最大观测间隔Tmaxgap、对目标的最小观测间隔Tmingap、对目标的平均观测间隔Tavegap、对目标的最大/最小观测间隔的均方差△Tmaxgap和△Tmingap
其中,
△Tmaxgap=|Tmaxgap-Tavegap| (12)
△Tmingap=|Tavegap-Tmingap| (13)
建立星座对目标的观测性能优化目标模型为:
fobersvation=λ1max(Ctimes)+λ2max(Tmean)+λ3min(Tmaxgap)+
λ4min(Tmingap)+λ5min(Tavegap)+λ6min(△Tmaxgap)+λ7min(△Tmingap) (14)
其中λ1、λ2、λ3、λ4、λ5、λ6和λ7分别是每天对目标的观测次数、每次平均观测时间、对目标的最大观测间隔、对目标的最小观测间隔、对目标的平均观测间隔和对目标观测间隔的均方差的加权系数值。
(5)卫星对目标成像性能
卫星目标的成像性能主要包括:空间分辨率L、幅宽E和视场角。其中空间分辨率和幅宽是与轨道高度、成像载荷密切相关的。而视场角主要影响镜头的设计,不作为星座构型设计的优化目标。
建立星座卫星对目标成像性能优化目标模型:
fimage=λ1min(L)+λ2max(E) (15)
其中,λ1和λ2分别是空间分辨率和幅宽的加权系数值。
基于以上步骤所得的各个优化设计模型即公式(5)~(15),即为针对低倾角回归轨道的,对目标位置的覆盖、入轨、重访、分辨率、幅宽等性能最优的卫星星座性能的优化设计模型:
fperformance=ξ1fcost+ξ2flaunch+ξ3fpropulsion+ξ4fobersvation+ξ5fimage (16)
其中ξ1、ξ2、ξ3、ξ4和ξ5分别为各性能项的加权系数值;
该程序能够根据星座构型优化设计的设置参数和变量参数,计算得到了星座性能仿真结果。
优选的,星座构型优化设计参数设置:目标点经度设置为135°,目标点纬度设置为30°,每天仿真步数86400步,仿真时长为回归轨道的回归天数D*,是回归轨道,轨道高度下限250km,轨道高度上限1200km;低轨卫星光学成像载荷视频摄像机的参数设置:CMOS像元大小6μm,镜头焦距2000mm,CMOS分辨率设置为4096×3072,面阵对角线长度为30mm。按此设置进行星座构型优化,能提高优化效率,缩短优化时间。
此处的输入约束条件和输出约束条件按常规方法进行设置,具体可以参看具体算例。星座构型设计优化目标可以根据实际需要进行定义,例如可以为具体算例中所列。
此处的多岛遗传算法为本领域常用的多岛遗传算法。优选的,多岛遗传算法的优化参数参数:子代大小设置为5、岛数设置为5、遗传代数设置为200,则优化的总步数为5×5×200=5000;此时计算效率达到最高。所需计算步骤最少。
具体算例如下:
根据低轨成像卫星任务特点,考虑卫星星座载荷为视频摄像机,利用卫星星座实现对目标的等时间间隔快速重访观测。对地观测目标位置的经纬度为(135°,30°),对地观测最小仰角为5°。任务的主要约束为实现目标的每次观测的间隔时间不超过30分钟,平均每次的观测时间不小于5分钟,卫星入轨初始时间为2006年1月1日。
步骤一:完成卫星星座性能计算程序设计。
步骤二:星座构型优化设计参数设置如表1所示。低轨卫星光学成像载荷视频摄像机的参数设置如表2所示。
表1 优化设计参数表
序号 | 参数 | 设置 |
1 | 目标点经度/° | 135 |
2 | 目标点纬度/° | 30 |
3 | 每天仿真步数 | 86400 |
4 | 仿真时长/天 | 回归天数D |
5 | 是否回归轨道 | 是 |
6 | 轨道高度下限/km | 250 |
7 | 轨道高度上限/km | 1200 |
表2 摄像机载荷参数表
序号 | 参数 | 设置 |
1 | CMOS像元大小 | 6μm |
2 | 镜头焦距 | 2000mm |
3 | CMOS分辨率 | 4096×3072 |
4 | 面阵对角线长度 | 30mm |
输入约束条件如表3所示。输出约束条件如表4所示。
表3 星座构型设计输入约束条件表
表4 星座构型设计输出约束条件表
输出约束 | 类型 | 下界 | 当前值 | 上界 |
轨道高度/km | REAL | 250 | 0 | 1200 |
卫星数目 | INTEGER | / | 0 | 14 |
每天观测次数 | REAL | 36 | 0 | / |
每次平均观测时间/s | REAL | 300.0 | 0 | / |
平均观测间隔时间/s | REAL | / | 0 | 1800 |
对地分辨率/m | REAL | / | 0 | 2.5 |
幅宽1/km | REAL | 4.0 | 0 | / |
幅宽2/km | REAL | 3.0 | 0 | / |
观测间隔偏差/s | REAL | / | 0 | 300 |
步骤三:星座构型优化设计目标函数组合如表5所示。
表5 星座构型设计优化目标定义表
目标参数 | 权重 | 刻度因子 | 类型 | 方向 | 当前值 |
轨道倾角/deg | 1 | 1 | REAL | minimize | 0 |
轨道高度/km | 1 | 1 | REAL | minimize | 0 |
卫星数目 | 1 | 1 | INTEGER | maximize | 0 |
每天观测次数 | 20 | 1 | INTEGER | maximize | 0 |
每次平均观测时间/s | 1 | 1 | REAL | minimize | 0 |
最大观测间隔时间/s | 0.01 | 1 | REAL | minimize | 0 |
最小观测间隔时间/s | 1 | 1 | REAL | maximize | 0 |
平均观测间隔时间/s | 1.0 | 1 | REAL | minimize | 0 |
最大观测间隔时间差/s | 0.01 | 1 | REAL | minimize | 0 |
最小观测间隔时间差/s | 0.01 | 1 | REAL | minimize | 0 |
星地距离/km | 1 | 1 | REAL | minimize | 0 |
入轨能耗/km/s | 1 | 1 | REAL | minimize | 0 |
空间分辨率/m | 1 | 1 | REAL | minimize | 0 |
幅宽1/km | 1 | 1 | REAL | maximize | 0 |
幅宽2/km | 1 | 1 | REAL | maximize | 0 |
最大最小观测间隔时间差/s | 20 | 1 | REAL | minimize | 0 |
步骤四:通过多岛遗传算法实现卫星星座构型的优化设计,多岛遗传算法的优化参数包括子代大小、岛数和遗传代数,分别设置为5、5和200,则优化的总步数为5×5×200=5000。
设置完成后,启动优化设计程序开始优化设计。程序运行结束后,得到了一组能够实现对目标的等时间间隔快速重访的卫星星座设计结果,如表6所示。
表6 卫星轨道优化设计结果表
序号 | 参数 | 数值 |
1 | 轨道面数目 | 3 |
2 | 每个轨道面卫星数目 | 4 |
3 | 回归圈数 | 15 |
4 | 回归周期/d | 1 |
5 | 轨道高度/km | 491.3632 |
6 | 轨道偏心率 | 0 |
7 | 轨道倾角/deg | 44.194 |
8 | 升交点赤经/deg | 22.9451 |
9 | 近地点幅角/deg | 0 |
10 | 平近点角/deg | 163.2001 |
11 | 相邻面升交点赤经差/deg | 96.96721 |
12 | 相邻面卫星平近点角差/deg | 7.574958 |
对优化设计的卫星星座进行轨道特性和覆盖方面的性能分析。卫星的星下点轨迹每1恒星日重复,卫星星座的所有卫星均匀分布在三条星下点轨迹上。考虑最小仰角5°,得到每个回归周期内,卫星星座对目标点的观测情况如图2所示,可以看出,卫星均匀等间隔实现对目标位置的快速重访,在每个回归周期内,卫星星座对目标可见82次,平均每次观测时间为463.573s,最小观测时间为35.690s,最大观测时间为576s。卫星星座对目标从发送请求到实现对观测目标的覆盖百分率如图3所示。可以看出,从观测任务请求发送,到卫星星座实现对目标观测的最大时间间隔小于1200s,也即卫星星座对该目标每隔不超过20分钟就可以实现一次观测。
本领域技术人员将清楚本发明的范围不限制于以上讨论的示例,有可能对其进行若干改变和修改,而不脱离所附权利要求书限定的本发明的范围。尽管己经在附图和说明书中详细图示和描述了本发明,但这样的说明和描述仅是说明或示意性的,而非限制性的。本发明并不限于所公开的实施例。
通过对附图,说明书和权利要求书的研究,在实施本发明时本领域技术人员可以理解和实现所公开的实施例的变形。在权利要求书中,术语“包括”不排除其他步骤或元素,而不定冠词“一个”或“一种”不排除多个。在彼此不同的从属权利要求中引用的某些措施的事实不意味着这些措施的组合不能被有利地使用。权利要求书中的任何参考标记不构成对本发明的范围的限制。
Claims (3)
1.一种等时间间隔快速重访卫星星座构型优化设计方法,其特征在于,包括以下步骤:
步骤S100:构建星座构型模型、回归轨道模型和性能优化模型;
步骤S200:设置星座构型优化设计参数、星座构型设计输入约束条件、星座构型设计输出约束条件和优化目标函数;
步骤S300:将各项参数代入公式(5)~(15),通过多岛遗传算法对卫星星座构型进行优化设计,得到能实现优化目标的等时间间隔快速重访的卫星星座构型设计结果;
其中步骤S100包括以下步骤:
步骤S110:构建星座构型模型:
参考Walker星座构型,建立对地观测卫星星座构型模型:N颗卫星均匀分布在P个轨道面上,每个轨道面上设置S颗卫星,且所形成的星座中所有卫星均具有相同的轨道高度、轨道偏心率、轨道倾角和近地点幅角;所设置的星座中各卫星之间的各轨道面之间的升交点赤经差,用△Ω表示,所设置的星座中各卫星之间的各轨道面的首颗卫星之间的平近点角差,用△M表示;
假设第1轨道面第1颗卫星的轨道参数为a0、e0、i0、Ω0、ω0、M0,对于圆轨道,卫星轨道偏心率e0和近地点幅角ω0取值:
e0=0,ω0=0 (1)
则所构建星座中的第j个轨道面中第k颗卫星的轨道参数按公式(2)计算得到:
其中j=[1,2,…,P],k=[1,2,…,S],P为轨道面数目,S为每个轨道面的卫星数目;
所构建星座中的卫星总数N满足公式(3):
N=PS (3)
步骤S120:构建回归轨道模型:
所构建星座中卫星需满足公式(4)所示的回归特性方程:
其中N*是回归圈数,D*是回归周期,ωe是地球自转的角速度,是卫星轨道的升交点赤经的长期变化速率,TΩ为轨道交点周期;
步骤S130:性能优化模型:
低轨光学卫星需考虑的性能指标包括卫星成本、发射费用、入轨能耗和卫星对目标的覆盖性能、卫星对目标的成像性能进行优化设计;
(1)卫星成本
影响卫星星座系统成本的主要因素是卫星轨道数目N,建立卫星成本优化目标模型为:
fcost=min(N) (5)
(2)发射费用
建立卫星发射费用成本模型优化目标模型为:
flaunch=λ1min(P)+λ2min(h) (6)
其中,λ1和λ2分别是轨道面数目和轨道高度的加权系数,P为轨道面数目和h为卫星轨道高度;
(3)入轨能耗
影响入轨能耗的主要因素是卫星轨道高度;
轨道越低,则发射卫星所需的能耗就越低;
建立卫星发射入轨能耗优化目标模型为:
fpropulsion=min(△V) (7)
其中,△V为卫星入轨所需的速度增量:
其中RE为地球半径,μ为地球引力系数,al发射轨道半长轴,a0为卫星运行轨道半长轴;且al为
式中
其中βL为发射角;
(4)对目标的覆盖性能
星座对目标的观测性能指标:每天对目标的观测次数Ctimes、每次平均观测时间Tmean、对目标的最大观测间隔Tmaxgap、对目标的最小观测间隔Tmingap、对目标的平均观测间隔Tavegap、对目标的最大观测间隔的均方差△Tmaxgap、最小观测间隔的均方差△Tmingap,其中,
△Tmaxgap=|Tmaxgap-Tavegap| (12)
△Tmingap=|Tavegap-Tmingap| (13)
建立星座对目标的观测性能优化目标模型为:
fobersvation=λ1max(Ctimes)+λ2max(Tmean)+λ3min(Tmaxgap)+
λ4min(Tmingap)+λ5min(Tavegap)+λ6min(△Tmaxgap)+λ7min(△Tmingap) (14)
其中λ1、λ2、λ3、λ4、λ5、λ6和λ7分别是每天对目标的观测次数、每次平均观测时间、对目标的最大观测间隔、对目标的最小观测间隔、对目标的平均观测间隔和对目标观测间隔的均方差的加权系数值;
(5)卫星对目标成像性能
卫星目标的成像性能包括:空间分辨率L、幅宽E和视场角;
其中空间分辨率和幅宽与轨道高度、成像载荷密切相关;
建立星座卫星对目标成像性能优化目标模型:
fimage=λ1min(L)+λ2max(E) (15)
其中,λ1和λ2分别是空间分辨率和幅宽的加权系数值;
基于以上步骤所得的各个优化设计模型即公式(5)~(15),可得到针对低倾角回归轨道的,对目标位置的性能最优的卫星星座性能的优化设计模型:
fperformance=ξ1fcost+ξ2flaunch+ξ3fpropulsion+ξ4fobersvation+ξ5fimage (16)
其中ξ1、ξ2、ξ3、ξ4和ξ5分别为各性能项的加权系数值。
2.根据权利要求1所述的等时间间隔快速重访卫星星座构型优化设计方法,其特征在于,所述星座构型优化设计参数设置为:目标点经度设置为135°,目标点纬度设置为30°,每天仿真步数86400步,仿真时长为回归轨道的回归天数D*,轨道高度下限250km,轨道高度上限1200km;低轨卫星光学成像载荷视频摄像机的参数设置:CMOS像元大小6μm,镜头焦距2000mm,CMOS分辨率设置为4096×3072,面阵对角线长度为30mm。
3.根据权利要求1所述的等时间间隔快速重访卫星星座构型优化设计方法,其特征在于,所述多岛遗传算法中的优化参数:子代大小设置为5、岛数设置为5、遗传代数设置为200。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610514705.4A CN106156417B (zh) | 2016-07-04 | 2016-07-04 | 等时间间隔快速重访卫星星座构型优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610514705.4A CN106156417B (zh) | 2016-07-04 | 2016-07-04 | 等时间间隔快速重访卫星星座构型优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106156417A true CN106156417A (zh) | 2016-11-23 |
CN106156417B CN106156417B (zh) | 2019-07-09 |
Family
ID=57351056
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610514705.4A Expired - Fee Related CN106156417B (zh) | 2016-07-04 | 2016-07-04 | 等时间间隔快速重访卫星星座构型优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106156417B (zh) |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106802667A (zh) * | 2016-12-30 | 2017-06-06 | 深圳航天东方红海特卫星有限公司 | 一种基于双停泊轨道的Walker星座部署方法 |
CN107329146A (zh) * | 2017-07-05 | 2017-11-07 | 中国人民解放军装备学院 | 一种导航卫星低轨监测星座的优化设计方法 |
CN107370535A (zh) * | 2017-06-26 | 2017-11-21 | 航天东方红卫星有限公司 | 基于微小卫星系统的天地一体化及时响应信息获取方法 |
CN107609292A (zh) * | 2017-09-19 | 2018-01-19 | 中国地质大学(武汉) | 一种基于指标评价的遥感卫星星座优化方法 |
CN107798187A (zh) * | 2017-10-24 | 2018-03-13 | 北京理工大学 | 一种高效卫星星座多目标优化方法 |
CN108932384A (zh) * | 2018-07-04 | 2018-12-04 | 北京市遥感信息研究所 | 基于一次覆盖条带最大化的成像卫星区域目标覆盖方法 |
CN109492252A (zh) * | 2018-10-08 | 2019-03-19 | 湖北航天技术研究院总体设计所 | 一种面向多目标区域的分步式星座构型优化设计方法 |
CN109656133A (zh) * | 2018-12-06 | 2019-04-19 | 上海航天控制技术研究所 | 一种针对空间走廊跟踪观测的分布式卫星群优化设计方法 |
CN109861742A (zh) * | 2019-03-01 | 2019-06-07 | 上海微小卫星工程中心 | 一种用于确定星座的构型的方法 |
CN110017815A (zh) * | 2019-04-30 | 2019-07-16 | 北京空间飞行器总体设计部 | 一种天基空间目标监视光电篱笆系统的构建方法 |
CN111177948A (zh) * | 2020-01-15 | 2020-05-19 | 东方红卫星移动通信有限公司 | 一种利用轨道要素进行描述的典型Walker星座数学模型 |
CN111324851A (zh) * | 2020-02-19 | 2020-06-23 | 中国科学院遥感与数字地球研究所 | 一种赤道轨道卫星大气层边界观测面积的计算方法 |
CN111795690A (zh) * | 2020-04-27 | 2020-10-20 | 中国人民解放军战略支援部队航天工程大学 | 基于随机时间Petri网的星座备份策略评估方法及系统 |
CN112596119A (zh) * | 2020-11-20 | 2021-04-02 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种电离层探测小卫星网络规划方法 |
CN113098590A (zh) * | 2021-04-19 | 2021-07-09 | 西安电子科技大学 | 面向非均匀业务覆盖的卫星星座构型参数优化方法 |
CN113572513A (zh) * | 2021-07-23 | 2021-10-29 | 中国人民解放军战略支援部队航天工程大学 | 低轨大规模通信卫星星座构型性能分析方法及装置 |
CN113589832A (zh) * | 2021-06-30 | 2021-11-02 | 北京宇航系统工程研究所 | 对地表固定区域目标稳定观测覆盖的星座快速设计方法 |
CN114026386A (zh) * | 2019-06-04 | 2022-02-08 | 弗劳恩霍夫应用研究促进协会 | 确定包含在卫星中的传感器访问目标区域的可能性的方法及卫星访问系统 |
CN114248949A (zh) * | 2021-11-22 | 2022-03-29 | 浙江大学 | 异面质比星座的部署方法及装置、电子设备、存储介质 |
CN114819258A (zh) * | 2022-03-11 | 2022-07-29 | 北京航空航天大学 | 一种面向巨型星座的分阶段部署优化方法和系统 |
JP7507923B2 (ja) | 2020-03-10 | 2024-06-28 | 三菱電機株式会社 | 宇宙情報レコーダー、情報管理方法、ssa事業装置、およびオープンアーキテクチャデータリポジトリ |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101976290A (zh) * | 2010-11-01 | 2011-02-16 | 北京航空航天大学 | 基于分解思想和粒子群融合方法的导航星座优化设计平台及方法 |
CN104835011A (zh) * | 2015-05-13 | 2015-08-12 | 中国西安卫星测控中心 | 基于地面站布局约束的导航星座缓变星间链路规划方法 |
CN105335541A (zh) * | 2014-08-12 | 2016-02-17 | 中国人民解放军装备学院 | 导航卫星星座的工程设计方法 |
CN105573118A (zh) * | 2015-12-16 | 2016-05-11 | 中国人民解放军国防科学技术大学 | 快速重访卫星轨道设计方法 |
-
2016
- 2016-07-04 CN CN201610514705.4A patent/CN106156417B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101976290A (zh) * | 2010-11-01 | 2011-02-16 | 北京航空航天大学 | 基于分解思想和粒子群融合方法的导航星座优化设计平台及方法 |
CN105335541A (zh) * | 2014-08-12 | 2016-02-17 | 中国人民解放军装备学院 | 导航卫星星座的工程设计方法 |
CN104835011A (zh) * | 2015-05-13 | 2015-08-12 | 中国西安卫星测控中心 | 基于地面站布局约束的导航星座缓变星间链路规划方法 |
CN105573118A (zh) * | 2015-12-16 | 2016-05-11 | 中国人民解放军国防科学技术大学 | 快速重访卫星轨道设计方法 |
Non-Patent Citations (3)
Title |
---|
CUI HONGZHENG等: "Satellite Constellation Configuration Design with Rapid Performance Calculation and Ordinal Optimization", 《CHINESE JOURNAL OF AERONAUTICS》 * |
宋伟等: "空间跟踪与监视系统卫星星座设计研究", 《战术导弹技术》 * |
蒙波等: "区域覆盖特殊椭圆轨道星座优化设计", 《北京航空航天大学学报》 * |
Cited By (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106802667A (zh) * | 2016-12-30 | 2017-06-06 | 深圳航天东方红海特卫星有限公司 | 一种基于双停泊轨道的Walker星座部署方法 |
CN106802667B (zh) * | 2016-12-30 | 2020-03-10 | 深圳航天东方红海特卫星有限公司 | 一种基于双停泊轨道的Walker星座部署方法 |
CN107370535B (zh) * | 2017-06-26 | 2019-11-29 | 航天东方红卫星有限公司 | 基于微小卫星系统的天地一体化及时响应信息获取方法 |
CN107370535A (zh) * | 2017-06-26 | 2017-11-21 | 航天东方红卫星有限公司 | 基于微小卫星系统的天地一体化及时响应信息获取方法 |
CN107329146A (zh) * | 2017-07-05 | 2017-11-07 | 中国人民解放军装备学院 | 一种导航卫星低轨监测星座的优化设计方法 |
CN107329146B (zh) * | 2017-07-05 | 2021-06-15 | 中国人民解放军战略支援部队航天工程大学 | 一种导航卫星低轨监测星座的优化设计方法 |
CN107609292A (zh) * | 2017-09-19 | 2018-01-19 | 中国地质大学(武汉) | 一种基于指标评价的遥感卫星星座优化方法 |
CN107609292B (zh) * | 2017-09-19 | 2020-04-10 | 中国地质大学(武汉) | 一种基于指标评价的遥感卫星星座优化方法 |
CN107798187A (zh) * | 2017-10-24 | 2018-03-13 | 北京理工大学 | 一种高效卫星星座多目标优化方法 |
CN107798187B (zh) * | 2017-10-24 | 2019-08-02 | 北京理工大学 | 一种高效卫星星座多目标优化方法 |
CN108932384A (zh) * | 2018-07-04 | 2018-12-04 | 北京市遥感信息研究所 | 基于一次覆盖条带最大化的成像卫星区域目标覆盖方法 |
CN109492252A (zh) * | 2018-10-08 | 2019-03-19 | 湖北航天技术研究院总体设计所 | 一种面向多目标区域的分步式星座构型优化设计方法 |
CN109656133B (zh) * | 2018-12-06 | 2022-01-14 | 上海航天控制技术研究所 | 一种针对空间走廊跟踪观测的分布式卫星群优化设计方法 |
CN109656133A (zh) * | 2018-12-06 | 2019-04-19 | 上海航天控制技术研究所 | 一种针对空间走廊跟踪观测的分布式卫星群优化设计方法 |
CN109861742A (zh) * | 2019-03-01 | 2019-06-07 | 上海微小卫星工程中心 | 一种用于确定星座的构型的方法 |
CN109861742B (zh) * | 2019-03-01 | 2021-10-26 | 上海微小卫星工程中心 | 一种用于确定星座的构型的方法 |
CN110017815B (zh) * | 2019-04-30 | 2021-08-10 | 北京空间飞行器总体设计部 | 一种天基空间目标监视光电篱笆系统的构建方法 |
CN110017815A (zh) * | 2019-04-30 | 2019-07-16 | 北京空间飞行器总体设计部 | 一种天基空间目标监视光电篱笆系统的构建方法 |
CN114026386A (zh) * | 2019-06-04 | 2022-02-08 | 弗劳恩霍夫应用研究促进协会 | 确定包含在卫星中的传感器访问目标区域的可能性的方法及卫星访问系统 |
CN111177948A (zh) * | 2020-01-15 | 2020-05-19 | 东方红卫星移动通信有限公司 | 一种利用轨道要素进行描述的典型Walker星座数学模型 |
CN111324851B (zh) * | 2020-02-19 | 2023-11-24 | 中国科学院遥感与数字地球研究所 | 一种赤道轨道卫星大气层边界观测面积的计算方法 |
CN111324851A (zh) * | 2020-02-19 | 2020-06-23 | 中国科学院遥感与数字地球研究所 | 一种赤道轨道卫星大气层边界观测面积的计算方法 |
JP7507923B2 (ja) | 2020-03-10 | 2024-06-28 | 三菱電機株式会社 | 宇宙情報レコーダー、情報管理方法、ssa事業装置、およびオープンアーキテクチャデータリポジトリ |
CN111795690A (zh) * | 2020-04-27 | 2020-10-20 | 中国人民解放军战略支援部队航天工程大学 | 基于随机时间Petri网的星座备份策略评估方法及系统 |
CN112596119A (zh) * | 2020-11-20 | 2021-04-02 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种电离层探测小卫星网络规划方法 |
CN113098590A (zh) * | 2021-04-19 | 2021-07-09 | 西安电子科技大学 | 面向非均匀业务覆盖的卫星星座构型参数优化方法 |
CN113098590B (zh) * | 2021-04-19 | 2022-03-04 | 西安电子科技大学 | 面向非均匀业务覆盖的卫星星座构型参数优化方法 |
CN113589832A (zh) * | 2021-06-30 | 2021-11-02 | 北京宇航系统工程研究所 | 对地表固定区域目标稳定观测覆盖的星座快速设计方法 |
CN113589832B (zh) * | 2021-06-30 | 2024-03-19 | 北京宇航系统工程研究所 | 对地表固定区域目标稳定观测覆盖的星座快速设计方法 |
CN113572513A (zh) * | 2021-07-23 | 2021-10-29 | 中国人民解放军战略支援部队航天工程大学 | 低轨大规模通信卫星星座构型性能分析方法及装置 |
CN114248949B (zh) * | 2021-11-22 | 2022-12-27 | 浙江大学 | 异面质比星座的部署方法及装置、电子设备、存储介质 |
CN114248949A (zh) * | 2021-11-22 | 2022-03-29 | 浙江大学 | 异面质比星座的部署方法及装置、电子设备、存储介质 |
CN114819258A (zh) * | 2022-03-11 | 2022-07-29 | 北京航空航天大学 | 一种面向巨型星座的分阶段部署优化方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106156417B (zh) | 2019-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106156417B (zh) | 等时间间隔快速重访卫星星座构型优化设计方法 | |
McElwain et al. | The James Webb Space Telescope mission: optical telescope element design, development, and performance | |
Amiri et al. | The emirates Mars mission | |
Tauber et al. | Planck pre-launch status: The Planck mission | |
JP6199021B2 (ja) | サービスを提供するよう適合された少なくとも2つの衛星からなる1組の衛星を制御する方法およびシステム | |
US20110178756A1 (en) | Integrated, Predictive, Radiance Sensor Apparatus and Method | |
Paek | Reconfigurable satellite constellations for geo-spatially adaptive Earth observation missions | |
CN112649006A (zh) | 一种太阳同步圆轨道的轨道规划方法 | |
Kubik et al. | First in-flight results of Pleiades 1A innovative methods for optical calibration | |
Polanskey et al. | Dawn science planning, operations and archiving | |
Fujita et al. | Attitude maneuvering sequence design of high-precision ground target tracking control for multispectral Earth observations | |
Stahl et al. | Designing astrophysics missions for NASA’s Space Launch System | |
Salvignol et al. | The Ariel Mission: A mission of the European Space Agency for the characterization of exoplanets | |
Cervone et al. | Phase a design of the lumio spacecraft: a cubesat for observing and characterizing micro-meteoroid impacts on the lunar far side | |
Mukhayadi | Efficient and high precision momentum bias attitude control for small satellite | |
CN115062416A (zh) | 一种基于离散空间组合优化的回归轨道快速设计方法 | |
Han et al. | Dawn’s final mission at Ceres: navigation and mission design experience | |
Polanskey | Efficacy of the Dawn Vesta science plan | |
Speretta et al. | LUMIO: A CubeSat to Monitor Micro-meteroid Impacts on the Lunar Farside | |
Wertz et al. | NanoEye—Military Relevant Surveillance for Less Than $5 Million Total Recurring Mission Cost | |
Dankanich et al. | Planetary Balloon-Based Science Platform Evaluation and Program Implementation | |
Kardashev et al. | Orbit design for the Spektr-R spacecraft of the ground-space interferometer | |
Nag et al. | Subsystem support feasibility for formation flight measuring Bi-directional Reflectance | |
CN115048817B (zh) | 一种对地威慑分析的决策辅助方法及系统 | |
Mortari | Optimization of flower constellations for dual use |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190709 |