CN107798187A - 一种高效卫星星座多目标优化方法 - Google Patents
一种高效卫星星座多目标优化方法 Download PDFInfo
- Publication number
- CN107798187A CN107798187A CN201710999619.1A CN201710999619A CN107798187A CN 107798187 A CN107798187 A CN 107798187A CN 201710999619 A CN201710999619 A CN 201710999619A CN 107798187 A CN107798187 A CN 107798187A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- mfrac
- 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
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
本发明公开的一种高效卫星星座多目标优化方法,属于航天器星座系统领域。本发明基于Walker‑δ星座构型,确定初始条件,建立星座轨道动力学方程,建立对地覆盖分析模型,建立对地观测分辨率模型;采用序列径向基函数多目标优化策略,以覆盖百分比和地面像元分辨率为目标函数,对轨道高度、轨道倾角、升交点赤经优化;基于l2加权和改进的Pareto适应度函数构造目标函数,采用RBF代理模型代替高耗时星座性能仿真模型优化设计,通过兴趣区间内序列采样对RBF代理模型更新与管理,从而获取满足工程需求的Pareto非劣解集作为卫星星座设计方案,实现星座对目标观测区域的覆盖百分比尽可能高以及有效载荷的像元分辨率尽可能小,降低卫星星座计算成本和设计成本、提高Pareto前沿搜索能力。
Description
技术领域
本发明涉及一种高效卫星星座多目标优化方法,属于航天器星座系统设计技术领域。
背景技术
卫星星座是指由多颗卫星按照一定规则和形状构成的可提供一定覆盖性能的卫星网络,是多颗卫星进行协同工作的基本形式。为了适应通信、导航和地球观测等领域的广泛需求,卫星星座技术得到了广泛重视。在这之中,卫星星座设计是卫星星座系统建立的前提和关键。卫星星座设计需要满足系统各项性能指标要求,如覆盖性能、观测性能、星间链路以及系统成本等,是一类典型的多目标设计优化问题。在星座系统概念设计阶段,目前的主要方法是基于多目标进化算法(如NSGA-II等)对星座性能指标进行优化,获取一组互不支配的Pareto非劣设计。然而,由于搜索过程的随机性,多目标进化算法往往需要大量调用星座系统分析模型以获取一组优化结果。此外,随着星座设计问题规模的增大,随机进化算法所需的模型调用次数迅速增加,呈现所谓的“维度灾难”。考虑到星座系统设计中需要调用各类计算耗时的仿真分析模型,直接采用多目标进化算法进行求解会显著降低设计效率,甚至导致星座多目标优化设计问题求解失败。因此,为了克服传统方法的局限性,十分有必要发展一种计算效率高、全局收敛性强的卫星星座多目标优化设计方法,从而在总体设计阶段能够快速实现星座方案的设计、优选与修改,为卫星星座系统方案论证与总体设计提供科学的依据与参考。
为了更好的说明本发明的技术方案,下面对所涉及到的相关数学工具进行简要介绍:
(1)径向基函数代理模型
径向基函数(RBF)是一种插值型代理模型,其基本形式如式(1)所示
其中ns为样本点数量,φ(||x-xi||),i=1,2..ns为基函数,β为RBF系数向量。常用的基函数类型如表1所示。
RBF需要在样本点处满足式(2)所示的插值条件
其中为样本点处的真实模型响应。
表1常用径向基函数类型
基函数类型 | 数学形式 |
线性基函数 | φ(r)=(r+c) |
高斯基函数 | φ(r)=exp(-cr2) |
样条基函数 | φ(r)=r2log(cr2) |
三次基函数 | φ(r)=(r+c)3 |
多二次基函数 | φ(r)=(r+c)1/2 |
逆多二次基函数 | φ(r)=(r+c)-1/2 |
(2)基于支持向量机的序列RBF优化策略
基于支持向量机的序列RBF优化策略(SRBF-SVM)是一种基于自适应RBF代理模型的高效全局优化算法。SRBF-SVM通过确定可能存在全局最优解的兴趣区间(ISR),在ISR内新增样本点并更新目标函数和约束函数的代理模型直至优化过程收敛。SRBF-SVM的关键技术是ISR辨识技术,其主要思想是基于已有样本点及其响应值信息,训练支持向量0-1分类机,并使用训练好的支持向量机对获取样本点进行分类,获取位于全局最优解附近的潜在优质样本点。通过特征空间下的模糊聚类方法(FS-FCM)计算潜在优质数据点的聚类中心,并结合当前最优解信息,辨识出ISR在设计空间中的位置。在构造ISR时,首先根据目标函数预测值阈值fthresh与样本点的目标函数真实响应f(xi),将已有样本点分为两类,如式(3)所示,
其中flagi为样本点xi的分类标识符。根据已有样本点分类标识符信息,训练目标函数的0-1SVM分类机。通过拉丁超方实验设计(LHD)方法在设计空间内获取海量样本点,并通过训练好的SVM对获取的样本点进行分类,得到一组目标函数响应值可能小于fthresh的潜在优质样本点,如图1(a)所示。在此基础上,采用式(4)所示的特征空间模糊聚类方法(FS-FCM)计算潜在优质数据点的聚类中心x*,
其中c=1为聚类中心个数,vi是第i个聚类中心,σ2为高斯核函数参数。以当前最优解x(k)为ISR中心,以x(k)与x*之间的欧氏距离确定ISR半径,从而构造出当前的ISR,如图1(b)所示。
发明内容
为克服传统多目标进化算法求解卫星星座设计问题中存在的计算量大、有限计算成本下Pareto前沿探索能力不足等问题,本发明公开的一种高效卫星星座多目标优化方法要解决的技术问题为:实现星座对目标观测区域的覆盖百分比尽可能高以及有效载荷的像元分辨率尽可能小,并降低卫星星座计算成本和设计成本、提高Pareto前沿搜索能力。
本发明的目的是通过下述技术方案实现的。
本发明公开的一种高效卫星星座多目标优化方法,基于Walker-δ星座构型,确定卫星星座多目标优化设计问题初始条件,建立星座轨道动力学方程,建立卫星星座对地覆盖分析模型,建立卫星星座对地观测分辨率模型;采用序列径向基函数(RBF)多目标优化策略,以覆盖百分比Cr和地面像元分辨率Rp为目标函数,对轨道高度h、轨道倾角i、升交点赤经R进行优化。此外,基于l2加权和改进的Pareto适应度函数构造多目标优化问题的目标函数,并采用RBF代理模型代替高耗时星座性能仿真模型进行优化设计,在此基础上通过兴趣区间内序列采样对RBF代理模型进行更新与管理,从而高效获取一组满足工程需求的Pareto非劣解集作为卫星星座设计方案,实现星座对目标观测区域的覆盖百分比尽可能高以及有效载荷的像元分辨率尽可能小,进而降低卫星星座计算成本和设计成本、提高Pareto前沿搜索能力。
本发明公开的一种高效卫星星座多目标优化方法,包括如下步骤:
步骤A:确定卫星星座多目标优化设计问题初始条件,具体实现方法如下:
步骤A-1:确定Walker-δ星座构型的固定参数,包括卫星总数T、轨道面数P和相位参数F,以及需要观测的区域经纬度范围。将轨道倾角i、轨道高度h、以及首颗卫星升交点赤经R作为设计变量,并确定各设计变量对应的取值范围[imin,imax],[hmin,hmax]和[Ωmin,Ωmax]。在此基础上,以对地观测覆盖百分比Cr和有效载荷地面像元分辨率RP为性能指标,建立卫星星座多目标优化设计问题数学模型如式(5)所示。
步骤A-2:确定序列RBF多目标优化策略参数,包括初始样本点数量、每次新增样本点数量,以及最大模型调用次数。
步骤B:考虑各种摄动力的影响下,建立高斯型星座轨道动力学方程。
其中a、e、i、Ω、ω、M为经典轨道根数,f为真近点角,u=ω+f为轨道幅角,p=a(1-e2)为轨道通径,A=[AR,AT,AN]为摄动力加速度在轨道坐标系RTN下的分量。在星座系统总体设计阶段,仅考虑地球非球形引力J2项的影响,如式(7)所示,
其中,J2为地球非球形引力J2项系数,μ为地球常数,aE为地球半径。
所述的各种摄动力包括地球包括:1)地球非球形引力摄动;2)太阳光压摄动;3)大气阻力摄动;4)地球磁场造成的磁力矩摄动等。
步骤C:建立卫星星座对地覆盖分析模型,具体建模方法如下:
步骤C-1:根据卫星轨道高度h和最小观测角σ,计算单颗卫星的覆盖角dσ、中心角ασ和覆盖带宽Swσ,在此基础上得到地面观测点T和卫星对应的地心角θ用以判断该点是否位于卫星的覆盖范围内。
步骤C-2:在观测区域内以预设经纬度间隔作网格图,将得到的网格点作为特征点。采用覆盖百分比Cr作为星座对地观测覆盖性能的指标,星座对地观测覆盖性能的指标等于观测区域内所有特征点被一颗或多颗卫星覆盖的总时间除以总仿真时间。
步骤C建立卫星星座对地覆盖分析模型如下:
步骤C-1:设卫星S某时刻高度为h,相应的星下点为G。按照式(8)分别计算考虑最小观测角σ下的覆盖角dσ、中心角ασ和覆盖带宽Swσ。
在此基础上,定义当前地面观测点的经度为λW、纬度为卫星星下点的经度为λS、纬度为按照式(9)计算地面观测点T和卫星对应的地心角θ。若满足θ≤dσ,则说明该地面观测点位于卫星S的覆盖范围内。
步骤C-2:在观测区域内以一定经纬度间隔作网格图,将得到的网格点作为特征点。定义两特征点之间的地球表面弧长为Δl,按照式(10)计算各特征点之间的经纬度间隔,其中λ为经度,为纬度。在此基础上,采用覆盖百分比Cr作为星座对地观测覆盖性能的指标,覆盖百分比Cr等于观测区域内所有特征点被一颗或多颗卫星覆盖的总时间除以总仿真时间。
步骤D:建立卫星星座对地观测分辨率模型。根据卫星轨道高度h和有效载荷像元尺寸d及光学系统焦距f,得到有效载荷的地面像元分辨率Rp作为星座对地观测分辨率性能的指标。
步骤E:采用多目标序列RBF优化方法,以覆盖百分比最大和地面像元分辨率最小为优化目标,对星座参数进行优化。
步骤E具体实现方法如下:
步骤E-1:在设计空间内采用拉丁超方设计(LHD)构造初始样本点,计算样本点处的各个目标函数响应值,即覆盖百分比Cr和地面像元分辨率Rp,并对Rp进行归一化处理得到无量纲的地面像元分辨率
步骤E-2:采用l2加权方法构造单目标优化函数F(x),如式(12)所示,
其中w1和w2分别为对应目标函数的权重,优选1。
步骤E-3:采用基于支持向量机的序列径向基函数优化策略(SRBF-SVM)对式(12)中的优化问题进行优化。首先基于已有样本点,构造目标函数的RBF代理模型在此基础上采用经典优化算法对进行优化,并根据当前优化结果信息构造兴趣采样空间(ISR)。通过在ISR内序列采样更新样本点数据库。
步骤E-4:根据式(13)构造改进的Pareto适应度函数MG(x),其中G(x)是样本点处的Pareto适应度,Gmax为所有样本点中的Pareto适应度最大值。
步骤E-5:采用SRBF-SVM优化策略对式(13)中的优化问题进行优化。首先基于已有样本点,构造改进Pareto适应度函数的RBF代理模型在此基础上采用经典优化算法对进行优化,并根据当前优化结果信息构造兴趣采样空间(ISR)。通过在ISR内序列采样更新样本点数据库。
步骤E-3和步骤E-5中的经典优化算法优选遗传算法。
步骤F:判断分析模型调用次数是否达到最大值。若未达到,则返回步骤E继续优化流程;反之,则优化流程终止,输出Pareto非劣解集作为一组卫星星座设计方案,实现星座对目标观测区域的高覆盖百分比以及较小的有效载荷的像元分辨率,进而降低卫星星座计算成本和设计成本、提高Pareto前沿搜索能力。
有益效果:
1、为克服传统多目标进化算法求解卫星星座设计问题中存在的计算量大等问题,本发明公开的一种高效卫星星座多目标优化方法,基于Walker-δ星座构型,以覆盖百分比和地面像元分辨率为目标函数,对轨道高度、轨道倾角、升交点赤经进行优化,保证星座对目标观测区域的覆盖百分比尽可能高并且有效载荷的像元分辨率尽可能小,采用多目标序列径向基函数优化方法,对目标函数构造代理模型,进而降低卫星星座计算成本和设计成本。
2、为克服了传统多目标进化算法求解卫星星座设计问题中存在的有限计算成本下Pareto前沿探索能力不足等问题,本发明公开的一种高效卫星星座多目标优化方法,基于l2加权和改进的Pareto适应度函数构造多目标优化问题的目标函数,并采用RBF代理模型代替高耗时星座性能仿真模型进行优化设计,在此基础上通过兴趣区间内序列采样对RBF代理模型进行更新与管理,从而引导优化过程快速收敛到多目标问题的Pareto前沿,提高有限计算成本下Pareto前沿探索能力。
3、本发明公开的一种高效卫星星座多目标优化方法,具有较好的分析精度、计算效率与工程实用性,适合应用于不同规模不同任务的Walker-δ星座构型优化设计,为实现总体设计阶段星座方案的快速设计、优化与论证提供有力的支撑。
附图说明
图1为兴趣采样空间示意图,其中图1(a)为优质样本点辨识示意图,图1(b)为兴趣采样空间构造示意图;
图2为高效卫星星座多目标优化设计方法流程图;
图3为星座对地覆盖分析模型示意图;
图4为案例星座多目标设计优化结果;
图5为优化后的卫星星座仿真示意图。
具体实施方式
为了更好的说明本发明的目的和优点,下面结合附图和实例对发明内容做进一步说明。
本发明提出并实现了一种基于代理模型的高效卫星星座多目标优化设计方法,该方法适用于Walker-δ星座构型的多目标优化设计问题,保证在总体设计阶段能够快速实现星座方案的设计优化与修改,为卫星星座系统方案论证和总体设计提供参考。
本实施例公开的一种高效卫星星座多目标优化方法,具体实现方法如下:
步骤A:采用Walker-δ构型作为卫星星座基准构型,其主要特征包括:各个轨道面的半长轴、轨道倾角和偏心率相同,且轨道面按升交点均匀分布;相邻轨道面间的卫星相位差相同,各轨道面内卫星数目相同,且各卫星的平近点角在该轨道面内均匀分布。在此基础上,根据实际任务需求,确定卫星总数T、轨道面数P和相位参数F,以及需要观测的区域经纬度范围,并将轨道倾角i、轨道高度h、以及首颗卫星升交点赤经R作为设计变量,确定各设计变量对应的取值范围[imin,imax],[hmin,hmax]和[Ωmin,Ωmax]。在此基础上,以对地观测覆盖百分比Cr和有效载荷地面像元分辨率RP为性能指标,建立卫星星座多目标优化设计问题数学模型如式(14)所示,
步骤B:考虑摄动因素的影响,建立高斯型星座轨道动力学方程如式(15)所示,
其中a、e、i、Ω、ω、M为经典轨道根数,f为真近点角,u=ω+f为轨道幅角,p=a(1-e2)为轨道通径,A=[AR,AT,AN]为摄动力加速度在轨道坐标系RTN下的分量。在星座系统总体设计阶段,仅考虑地球非球形引力J2项的影响,如式(16)所示,
其中,J2=1.082626836×10-3为地球非球形引力J2项系数,μ=398600.5km3/s2为地球常数,aE=6378km为地球半径。
步骤C:建立卫星星座对地覆盖分析模型,具体建模过程如下:
步骤C-1:设卫星S某时刻高度为h,相应的星下点为G。如图3所示,按照式(17)分别计算考虑最小观测角σ下的覆盖角dσ、中心角ασ和覆盖带宽Swσ。
在此基础上,设当前地面观测点的经度为λW、纬度为卫星星下点的经度为λS、纬度为按照式(18)计算地面观测点T和卫星对应的地心角θ。若满足θ≤dσ,则说明该地面观测点位于卫星S的覆盖范围内。
步骤C-2:在观测区域内以一定经纬度间隔作网格图,将得到的网格点作为特征点。设两特征点之间的地球表面弧长为Δl,按照式(19)计算各特征点之间的经纬度间隔,其中λ为经度,为纬度。在此基础上,采用覆盖百分比Cr作为星座对地观测覆盖性能的指标,其等于观测区域内所有特征点被一颗或多颗卫星覆盖的总时间除以总仿真时间。
步骤D:建立卫星星座对地观测分辨率模型。根据卫星轨道高度h和有效载荷像元尺寸d及光学系统焦距f,按照式(20)计算有效载荷的地面像元分辨率Rp作为星座对地观测分辨率性能的指标。
步骤E:采用多目标序列径向基函数(RBF)优化方法,以覆盖百分比最大和地面像元分辨率最小为优化目标,对星座参数进行优化,具体实施步骤如下:
步骤E-1:在设计空间内采用拉丁超方实验设计(LHD)方法构造初始样本点,计算样本点处的覆盖百分比Cr和地面像元分辨率Rp,并按照式(21)对Rp进行归一化处理得到无量纲的地面像元分辨率
步骤E-2:采用l2加权方法构造单目标优化函数F(x),如式(22)所示,
其中w1和w2分别为对应目标函数的权重,一般情况下取1。
步骤E-3:采用基于支持向量机的序列径向基函数优化策略(SRBF-SVM)对式(22)中的优化问题进行优化。首先基于已有样本点,构造目标函数的RBF代理模型在此基础上采用遗传算法对进行优化,并根据当前优化结果信息构造兴趣采样空间(ISR)。通过在ISR内序列采样更新样本点数据库。
步骤E-4:根据式(23)构造改进的Pareto适应度函数MG(x),其中G(x)是样本点处的Pareto适应度,Gmax为所有样本点中的Pareto适应度最大值。
步骤E-5:采用SRBF-SVM优化策略对式(23)中的优化问题进行优化。首先基于已有样本点,构造改进Pareto适应度函数的RBF代理模型在此基础上采用遗传算法对进行优化,并根据当前优化结果信息构造兴趣采样空间ISR。通过在ISR内序列采样更新样本点数据库。
步骤F:判断分析模型调用次数是否达到最大值。若未达到,则返回步骤E继续优化流程;反之,则优化流程终止,输出Pareto非劣解集作为一组卫星星座设计方案。
为了更好的体现本发明的有效性与工程实用性,下面通过具体的卫星星座多目标优化设计问题为例,结合附图与表格对本发明做进一步说明。
本案例中,需要覆盖的目标区域经纬度范围为73°E~135°E、3°N~53°N,仿真时间为1天,地面最小观测角为10°。卫星数量为T=18,轨道平面数P=3,相位参数F=1。卫星对地观测有效载荷像元尺寸d=10μm,焦距f=500mm。优化问题设计变量取值范围为:h∈[500km,2000km],i∈[30°,60°],Ω∈[0°,90°]。优化过程中,初始样本点数量为10,每次优化新增样本点数量为3,最大模型调用次数为100,得到卫星星座多目标设计优化结果如图4所示。从图中可以看出,本发明仅需要调用100次分析模型即可获得29组Pareto非劣解,而NSGA-II等传统进化算法在求解类似问题时往往需要调用5000次以上分析模型(戴光明,多目标优化算法及在卫星星座设计中的应用[M],中国地质大学出版社,2009.),其计算成本降低了90%以上。此外,图4中的Pareto非劣解分布情况表明,覆盖百分比和地面像元分辨率是一对矛盾的设计指标。随着星座覆盖百分比增加,卫星星座高度逐渐提升,从而影响地面像元分辨率。因此,设计者需要根据实际需求,在Pareto非劣解集中选择合适的设计点作为最终设计方案。从工程角度考虑,最终保留覆盖百分比高于90%的设计方案作为一组备选方案,如表2所示。实际工程中要求星座系统对地观测分辨率不大于30m,则选择方案2作为最终设计方案,其STK仿真示意图如图5所示。
表2覆盖百分比高于90%的星座设计方案
方案编号 | 轨道高度/km | 轨道倾角/° | 升交点赤经/° | 覆盖百分比/% | 地面像元分辨率/m |
1 | 1724.95 | 47.04 | 80.52 | 99.99 | 34.50 |
2 | 1396.82 | 41.75 | 23.06 | 96.85 | 27.94 |
3 | 1569.88 | 37.95 | 88.35 | 97.35 | 31.40 |
4 | 1724.32 | 55.93 | 6.44 | 98.42 | 34.49 |
上述优化设计结果表明,相比于传统方法,本发明能够以较小的计算成本获取一组满足实际工程需求并具有良好覆盖性能与分辨率指标的卫星星座设计方案,实现了预期的发明目的,验证了本发明的合理性、有效性和工程实用性。
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例,用于解释本发明,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种高效卫星星座多目标优化方法,其特征在于:包括如下步骤,
步骤A:确定卫星星座多目标优化设计问题初始条件,具体实现方法如下,
步骤A-1:确定Walker-δ星座构型的固定参数,包括卫星总数T、轨道面数P和相位参数F,以及需要观测的区域经纬度范围;将轨道倾角i、轨道高度h、以及首颗卫星升交点赤经R作为设计变量,并确定各设计变量对应的取值范围[imin,imax],[hmin,hmax]和[Ωmin,Ωmax];在此基础上,以对地观测覆盖百分比Cr和有效载荷地面像元分辨率RP为性能指标,建立卫星星座多目标优化设计问题数学模型如式(1)所示;
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>f</mi>
<mi>i</mi>
<mi>n</mi>
<mi>d</mi>
<mi> </mi>
<mi>x</mi>
<mo>=</mo>
<mo>&lsqb;</mo>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
<mo>,</mo>
<mi>&Omega;</mi>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mi>max</mi>
</mtd>
<mtd>
<msub>
<mi>C</mi>
<mi>r</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>min</mi>
</mtd>
<mtd>
<msub>
<mi>R</mi>
<mi>p</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>i</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>&le;</mo>
<mi>i</mi>
<mo>&le;</mo>
<msub>
<mi>i</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>h</mi>
<mi>min</mi>
</msub>
<mo>&le;</mo>
<mi>h</mi>
<mo>&le;</mo>
<msub>
<mi>h</mi>
<mi>max</mi>
</msub>
<mo>,</mo>
<msub>
<mi>&Omega;</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>&le;</mo>
<mi>&Omega;</mi>
<mo>&le;</mo>
<msub>
<mi>&Omega;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤A-2:确定序列RBF多目标优化策略参数,包括初始样本点数量、每次新增样本点数量,以及最大模型调用次数;
步骤B:考虑各种摄动力的影响下,建立高斯型星座轨道动力学方程;
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>a</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<mrow>
<mi>n</mi>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>e</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mi>R</mi>
</msub>
<mi>e</mi>
<mi> </mi>
<mi>sin</mi>
<mi> </mi>
<mi>f</mi>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>T</mi>
</msub>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>e</mi>
<mi> </mi>
<mi>cos</mi>
<mi> </mi>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>e</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>e</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mrow>
<mi>n</mi>
<mi>a</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mi>R</mi>
</msub>
<mi>sin</mi>
<mi> </mi>
<mi>f</mi>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>T</mi>
</msub>
<mo>(</mo>
<mi>cos</mi>
<mi> </mi>
<mi>f</mi>
<mo>+</mo>
<mi>cos</mi>
<mi> </mi>
<mi>E</mi>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>i</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>r</mi>
<mi> </mi>
<mi>cos</mi>
<mi> </mi>
<mi>u</mi>
</mrow>
<mrow>
<msup>
<mi>na</mi>
<mn>2</mn>
</msup>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>e</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mfrac>
<msub>
<mi>A</mi>
<mi>N</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>&Omega;</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>r</mi>
<mi> </mi>
<mi>sin</mi>
<mi> </mi>
<mi>u</mi>
</mrow>
<mrow>
<msup>
<mi>na</mi>
<mn>2</mn>
</msup>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>e</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mi>sin</mi>
<mi> </mi>
<mi>i</mi>
</mrow>
</mfrac>
<msub>
<mi>A</mi>
<mi>N</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>&omega;</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>e</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mrow>
<mi>n</mi>
<mi>a</mi>
<mi>e</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>A</mi>
<mi>R</mi>
</msub>
<mi>cos</mi>
<mi> </mi>
<mi>f</mi>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>T</mi>
</msub>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mi>r</mi>
<mi>p</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mi> </mi>
<mi>f</mi>
<mo>)</mo>
<mo>-</mo>
<mi>cos</mi>
<mi> </mi>
<mi>i</mi>
<mfrac>
<mrow>
<mi>d</mi>
<mi>&Omega;</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>M</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mi>n</mi>
<mo>-</mo>
<mfrac>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>e</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mi>n</mi>
<mi>a</mi>
<mi>e</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>A</mi>
<mi>R</mi>
</msub>
<mo>(</mo>
<mi>cos</mi>
<mi> </mi>
<mi>f</mi>
<mo>-</mo>
<mn>2</mn>
<mi>e</mi>
<mfrac>
<mi>r</mi>
<mi>p</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>T</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mi>r</mi>
<mi>p</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mi>sin</mi>
<mi> </mi>
<mi>f</mi>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中a、e、i、Ω、ω、M为经典轨道根数,f为真近点角,u=ω+f为轨道幅角,p=a(1-e2)为轨道通径,A=[AR,AT,AN]为摄动力加速度在轨道坐标系RTN下的分量;在星座系统总体设计阶段,仅考虑地球非球形引力J2项的影响,如式(3)所示,
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<msub>
<mi>RJ</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>3</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>J</mi>
<mn>2</mn>
</msub>
<mfrac>
<mi>&mu;</mi>
<msup>
<mi>r</mi>
<mn>4</mn>
</msup>
</mfrac>
<msubsup>
<mi>a</mi>
<mi>E</mi>
<mn>2</mn>
</msubsup>
<mrow>
<mo>(</mo>
<mn>3</mn>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>i</mi>
<mi> </mi>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>u</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<msub>
<mi>TJ</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mn>3</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>J</mi>
<mn>2</mn>
</msub>
<mfrac>
<mi>&mu;</mi>
<msup>
<mi>r</mi>
<mn>4</mn>
</msup>
</mfrac>
<msubsup>
<mi>a</mi>
<mi>E</mi>
<mn>2</mn>
</msubsup>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mi>i</mi>
<mi> </mi>
<mi>sin</mi>
<mn>2</mn>
<mi>u</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<msub>
<mi>NJ</mi>
<mn>2</mn>
</msub>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mn>3</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>J</mi>
<mn>2</mn>
</msub>
<mfrac>
<mi>&mu;</mi>
<msup>
<mi>r</mi>
<mn>4</mn>
</msup>
</mfrac>
<msubsup>
<mi>a</mi>
<mi>E</mi>
<mn>2</mn>
</msubsup>
<mi>sin</mi>
<mn>2</mn>
<mi>i</mi>
<mi> </mi>
<mi>sin</mi>
<mi>u</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,J2为地球非球形引力J2项系数,μ为地球常数,aE为地球半径;
所述的各种摄动力包括地球包括:1)地球非球形引力摄动;2)太阳光压摄动;3)大气阻力摄动;4)地球磁场造成的磁力矩摄动;
步骤C:建立卫星星座对地覆盖分析模型,具体建模方法如下:
步骤C-1:根据卫星轨道高度h和最小观测角σ,计算单颗卫星的覆盖角dσ、中心角ασ和覆盖带宽Swσ,在此基础上得到地面观测点T和卫星对应的地心角θ用以判断该点是否位于卫星的覆盖范围内;
步骤C-2:在观测区域内以预设经纬度间隔作网格图,将得到的网格点作为特征点;采用覆盖百分比Cr作为星座对地观测覆盖性能的指标,星座对地观测覆盖性能的指标等于观测区域内所有特征点被一颗或多颗卫星覆盖的总时间除以总仿真时间;
步骤D:建立卫星星座对地观测分辨率模型;根据卫星轨道高度h和有效载荷像元尺寸d及光学系统焦距f,得到有效载荷的地面像元分辨率Rp作为星座对地观测分辨率性能的指标;
<mrow>
<msub>
<mi>R</mi>
<mi>P</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>h</mi>
<mi>d</mi>
</mrow>
<mi>f</mi>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤E:采用多目标序列RBF优化方法,以覆盖百分比最大和地面像元分辨率最小为优化目标,对星座参数进行优化;
步骤F:判断分析模型调用次数是否达到最大值;若未达到,则返回步骤E继续优化流程;反之,则优化流程终止,输出Pareto非劣解集作为一组卫星星座设计方案,实现星座对目标观测区域的高覆盖百分比以及较小的有效载荷的像元分辨率,进而降低卫星星座计算成本和设计成本、提高Pareto前沿搜索能力。
2.如权利要求1所述的一种高效卫星星座多目标优化方法,其特征在于:步骤C建立卫星星座对地覆盖分析模型如下,
步骤C-1:设卫星S某时刻高度为h,相应的星下点为G;按照式(5)分别计算考虑最小观测角σ下的覆盖角dσ、中心角ασ和覆盖带宽Swσ;
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>d</mi>
<mi>&sigma;</mi>
</msub>
<mo>=</mo>
<mi>arccos</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>a</mi>
<mi>E</mi>
</msub>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mi>&sigma;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mi>E</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>&sigma;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&alpha;</mi>
<mi>&sigma;</mi>
</msub>
<mo>=</mo>
<mi>&pi;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>-</mo>
<msub>
<mi>d</mi>
<mi>&sigma;</mi>
</msub>
<mo>-</mo>
<mi>&sigma;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>w</mi>
<mi>&sigma;</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>2</mn>
<msub>
<mi>a</mi>
<mi>E</mi>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>d</mi>
<mi>&sigma;</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
在此基础上,定义当前地面观测点的经度为λW、纬度为卫星星下点的经度为λS、纬度为按照式(6)计算地面观测点T和卫星对应的地心角θ;若满足θ≤dσ,则说明该地面观测点位于卫星S的覆盖范围内;
步骤C-2:在观测区域内以一定经纬度间隔作网格图,将得到的网格点作为特征点;定义两特征点之间的地球表面弧长为Δl,按照式(7)计算各特征点之间的经纬度间隔,其中λ为经度,为纬度;在此基础上,采用覆盖百分比Cr作为星座对地观测覆盖性能的指标,覆盖百分比Cr等于观测区域内所有特征点被一颗或多颗卫星覆盖的总时间除以总仿真时间。
3.如权利要求1或2所述的一种高效卫星星座多目标优化方法,其特征在于:步骤E具体实现方法如下,
步骤E-1:在设计空间内采用拉丁超方设计(LHD)构造初始样本点,计算样本点处的各个目标函数响应值,即覆盖百分比Cr和地面像元分辨率Rp,并对Rp进行归一化处理得到无量纲的地面像元分辨率
步骤E-2:采用l2加权方法构造单目标优化函数F(x),如式(8)所示,
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>f</mi>
<mi>i</mi>
<mi>n</mi>
<mi>d</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>x</mi>
<mo>=</mo>
<mo>&lsqb;</mo>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
<mo>,</mo>
<mi>&Omega;</mi>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>min</mi>
</mtd>
<mtd>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<msub>
<mi>w</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>C</mi>
<mi>r</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>w</mi>
<mn>2</mn>
</msub>
<msubsup>
<mover>
<mi>R</mi>
<mo>~</mo>
</mover>
<mi>p</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>i</mi>
<mi>min</mi>
</msub>
<mo>&le;</mo>
<mi>i</mi>
<mo>&le;</mo>
<msub>
<mi>i</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>h</mi>
<mi>min</mi>
</msub>
<mo>&le;</mo>
<mi>h</mi>
<mo>&le;</mo>
<msub>
<mi>h</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>&Omega;</mi>
<mi>min</mi>
</msub>
<mo>&le;</mo>
<mi>&Omega;</mi>
<mo>&le;</mo>
<msub>
<mi>&Omega;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
其中w1和w2分别为对应目标函数的权重;
步骤E-3:采用基于支持向量机的序列径向基函数优化策略(SRBF-SVM)对式(8)中的优化问题进行优化;首先基于已有样本点,构造目标函数的RBF代理模型在此基础上采用经典优化算法对进行优化,并根据当前优化结果信息构造兴趣采样空间(ISR);通过在ISR内序列采样更新样本点数据库;
步骤E-4:根据式(9)构造改进的Pareto适应度函数MG(x),其中G(x)是样本点处的Pareto适应度,Gmax为所有样本点中的Pareto适应度最大值;
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>M</mi>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>G</mi>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
<mo>)</mo>
<mo>+</mo>
<mn>2</mn>
<mo>&times;</mo>
<msub>
<mi>G</mi>
<mi>max</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo><</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<munder>
<mi>max</mi>
<mrow>
<mi>j</mi>
<mo>&NotEqual;</mo>
<mi>k</mi>
</mrow>
</munder>
<mo>&lsqb;</mo>
<mi>min</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>R</mi>
<mo>~</mo>
</mover>
<mi>p</mi>
</msub>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>R</mi>
<mo>~</mo>
</mover>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>C</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>C</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>)</mo>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤E-5:采用SRBF-SVM优化策略对式(9)中的优化问题进行优化;首先基于已有样本点,构造改进Pareto适应度函数的RBF代理模型在此基础上采用经典优化算法对进行优化,并根据当前优化结果信息构造兴趣采样空间(ISR);通过在ISR内序列采样更新样本点数据库。
4.如权利要求3所述的一种高效卫星星座多目标优化方法,其特征在于:步骤E-3和步骤E-5中的经典优化算法选遗传算法。
5.一种高效卫星星座多目标优化方法,其特征在于:基于Walker-δ星座构型,确定卫星星座多目标优化设计问题初始条件,建立星座轨道动力学方程,建立卫星星座对地覆盖分析模型,建立卫星星座对地观测分辨率模型;采用序列径向基函数RBF多目标优化策略,以覆盖百分比Cr和地面像元分辨率Rp为目标函数,对轨道高度h、轨道倾角i、升交点赤经R进行优化;此外,基于l2加权和改进的Pareto适应度函数构造多目标优化问题的目标函数,并采用RBF代理模型代替高耗时星座性能仿真模型进行优化设计,在此基础上通过兴趣区间内序列采样对RBF代理模型进行更新与管理,从而高效获取一组满足工程需求的Pareto非劣解集作为卫星星座设计方案,实现星座对目标观测区域的覆盖百分比尽可能高以及有效载荷的像元分辨率尽可能小,进而降低卫星星座计算成本和设计成本、提高Pareto前沿搜索能力。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710999619.1A CN107798187B (zh) | 2017-10-24 | 2017-10-24 | 一种高效卫星星座多目标优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710999619.1A CN107798187B (zh) | 2017-10-24 | 2017-10-24 | 一种高效卫星星座多目标优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107798187A true CN107798187A (zh) | 2018-03-13 |
CN107798187B CN107798187B (zh) | 2019-08-02 |
Family
ID=61533624
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710999619.1A Active CN107798187B (zh) | 2017-10-24 | 2017-10-24 | 一种高效卫星星座多目标优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107798187B (zh) |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108984903A (zh) * | 2018-07-16 | 2018-12-11 | 北京航天自动控制研究所 | 一种制导参数的优选/优化设计方法 |
CN109615421A (zh) * | 2018-11-28 | 2019-04-12 | 安徽大学 | 一种基于多目标进化算法的个性化商品推荐方法 |
CN109656133A (zh) * | 2018-12-06 | 2019-04-19 | 上海航天控制技术研究所 | 一种针对空间走廊跟踪观测的分布式卫星群优化设计方法 |
CN109861742A (zh) * | 2019-03-01 | 2019-06-07 | 上海微小卫星工程中心 | 一种用于确定星座的构型的方法 |
CN109918761A (zh) * | 2019-03-01 | 2019-06-21 | 上海微小卫星工程中心 | 一种分析星座的覆盖性能的方法 |
CN109977576A (zh) * | 2019-04-03 | 2019-07-05 | 北京理工大学 | 一种基于代理模型的卫星星座系统多学科设计优化方法 |
CN109977456A (zh) * | 2019-01-31 | 2019-07-05 | 上海卫星工程研究所 | 面向观测任务的卫星组网方案辅助设计系统 |
CN110262241A (zh) * | 2019-06-28 | 2019-09-20 | 中国空间技术研究院 | 基于高斯过程预测控制的航天器轨道控制方法 |
CN110276159A (zh) * | 2019-07-01 | 2019-09-24 | 北京理工大学 | 一种基于多模型融合的卫星系统多学科优化方法 |
CN111177948A (zh) * | 2020-01-15 | 2020-05-19 | 东方红卫星移动通信有限公司 | 一种利用轨道要素进行描述的典型Walker星座数学模型 |
CN111259604A (zh) * | 2020-01-16 | 2020-06-09 | 中国科学院空间应用工程与技术中心 | 基于机器学习的高轨道卫星光压模型辨识方法及系统 |
CN111301715A (zh) * | 2020-01-17 | 2020-06-19 | 哈尔滨工业大学 | 基于霍曼变轨的同轨道特定相位分布的星座布局与轨道调整方法、装置及计算机存储介质 |
CN111881505A (zh) * | 2020-08-04 | 2020-11-03 | 河北工业大学 | 一种基于ga-rbf算法的既有建筑多目标优化改造决策方法 |
CN112085812A (zh) * | 2020-08-11 | 2020-12-15 | 中国科学院力学研究所 | 一种超大规模卫星星座组网通信动态节点确定方法 |
CN113011105A (zh) * | 2021-04-08 | 2021-06-22 | 东方红卫星移动通信有限公司 | 基于改进型布谷鸟算法的多目标星座模块优化方法 |
CN113591263A (zh) * | 2021-03-24 | 2021-11-02 | 长光卫星技术有限公司 | 一种基于几何解析法的广域目标搜索星座设计方法 |
CN115017386A (zh) * | 2022-08-08 | 2022-09-06 | 中国人民解放军63921部队 | 基于聚类的观测数据与目标库根数关联方法和装置 |
WO2022198748A1 (zh) * | 2021-03-23 | 2022-09-29 | 北京航空航天大学 | 一种面向araim应用的低轨卫星增强系统星座配置优化方法 |
CN115795927A (zh) * | 2023-02-10 | 2023-03-14 | 北京英视睿达科技股份有限公司 | 基于多目标观测的卫星轨道设计方法、装置、设备及介质 |
CN115994982A (zh) * | 2023-03-22 | 2023-04-21 | 中科星图测控技术股份有限公司 | 一种基于osg的巨型星座态势展示方法 |
CN117040607A (zh) * | 2023-10-10 | 2023-11-10 | 中国人民解放军战略支援部队航天工程大学 | 一种低轨通信卫星星座的设计方法 |
CN117521260A (zh) * | 2023-11-17 | 2024-02-06 | 北京工业大学 | 一种面向离散分布的多目标区域覆盖的卫星星座设计方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101887478A (zh) * | 2010-07-19 | 2010-11-17 | 北京理工大学 | 一种基于序列径向基代理模型的高效全局优化方法 |
CN101976290A (zh) * | 2010-11-01 | 2011-02-16 | 北京航空航天大学 | 基于分解思想和粒子群融合方法的导航星座优化设计平台及方法 |
CN103473424A (zh) * | 2013-09-23 | 2013-12-25 | 北京理工大学 | 基于序列径向基函数代理模型的飞行器系统优化设计方法 |
US20150071645A1 (en) * | 2012-07-13 | 2015-03-12 | Raytheon Company | High-bandwidth optical communications relay architecture |
CN104866692A (zh) * | 2015-06-18 | 2015-08-26 | 北京理工大学 | 一种基于自适应代理模型的飞行器多目标优化方法 |
CN106156417A (zh) * | 2016-07-04 | 2016-11-23 | 中国人民解放军国防科学技术大学 | 等时间间隔快速重访卫星星座构型优化设计方法 |
-
2017
- 2017-10-24 CN CN201710999619.1A patent/CN107798187B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101887478A (zh) * | 2010-07-19 | 2010-11-17 | 北京理工大学 | 一种基于序列径向基代理模型的高效全局优化方法 |
CN101976290A (zh) * | 2010-11-01 | 2011-02-16 | 北京航空航天大学 | 基于分解思想和粒子群融合方法的导航星座优化设计平台及方法 |
US20150071645A1 (en) * | 2012-07-13 | 2015-03-12 | Raytheon Company | High-bandwidth optical communications relay architecture |
CN103473424A (zh) * | 2013-09-23 | 2013-12-25 | 北京理工大学 | 基于序列径向基函数代理模型的飞行器系统优化设计方法 |
CN104866692A (zh) * | 2015-06-18 | 2015-08-26 | 北京理工大学 | 一种基于自适应代理模型的飞行器多目标优化方法 |
CN106156417A (zh) * | 2016-07-04 | 2016-11-23 | 中国人民解放军国防科学技术大学 | 等时间间隔快速重访卫星星座构型优化设计方法 |
Non-Patent Citations (3)
Title |
---|
彭磊: "桁架式静止轨道卫星平台多学科多任务设计优化", 《中国博士学位论文全文数据库工程科技Ⅱ辑》 * |
龙腾等: "多学科设计优化技术发展及在航空航天领域的应用", 《航空制造技术》 * |
龙腾等: "改进的Pareto多目标协同优化策略", 《系统工程与电子技术》 * |
Cited By (38)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108984903A (zh) * | 2018-07-16 | 2018-12-11 | 北京航天自动控制研究所 | 一种制导参数的优选/优化设计方法 |
CN108984903B (zh) * | 2018-07-16 | 2022-12-09 | 北京航天自动控制研究所 | 一种制导参数的优选/优化设计方法 |
CN109615421A (zh) * | 2018-11-28 | 2019-04-12 | 安徽大学 | 一种基于多目标进化算法的个性化商品推荐方法 |
CN109615421B (zh) * | 2018-11-28 | 2022-12-06 | 安徽大学 | 一种基于多目标进化算法的个性化商品推荐方法 |
CN109656133A (zh) * | 2018-12-06 | 2019-04-19 | 上海航天控制技术研究所 | 一种针对空间走廊跟踪观测的分布式卫星群优化设计方法 |
CN109656133B (zh) * | 2018-12-06 | 2022-01-14 | 上海航天控制技术研究所 | 一种针对空间走廊跟踪观测的分布式卫星群优化设计方法 |
CN109977456B (zh) * | 2019-01-31 | 2023-06-30 | 上海卫星工程研究所 | 面向观测任务的卫星组网方案辅助设计系统 |
CN109977456A (zh) * | 2019-01-31 | 2019-07-05 | 上海卫星工程研究所 | 面向观测任务的卫星组网方案辅助设计系统 |
CN109918761A (zh) * | 2019-03-01 | 2019-06-21 | 上海微小卫星工程中心 | 一种分析星座的覆盖性能的方法 |
CN109861742B (zh) * | 2019-03-01 | 2021-10-26 | 上海微小卫星工程中心 | 一种用于确定星座的构型的方法 |
CN109861742A (zh) * | 2019-03-01 | 2019-06-07 | 上海微小卫星工程中心 | 一种用于确定星座的构型的方法 |
CN109977576A (zh) * | 2019-04-03 | 2019-07-05 | 北京理工大学 | 一种基于代理模型的卫星星座系统多学科设计优化方法 |
CN110262241A (zh) * | 2019-06-28 | 2019-09-20 | 中国空间技术研究院 | 基于高斯过程预测控制的航天器轨道控制方法 |
CN110262241B (zh) * | 2019-06-28 | 2022-03-04 | 中国空间技术研究院 | 基于高斯过程预测控制的航天器轨道控制方法 |
CN110276159A (zh) * | 2019-07-01 | 2019-09-24 | 北京理工大学 | 一种基于多模型融合的卫星系统多学科优化方法 |
CN111177948A (zh) * | 2020-01-15 | 2020-05-19 | 东方红卫星移动通信有限公司 | 一种利用轨道要素进行描述的典型Walker星座数学模型 |
CN111259604A (zh) * | 2020-01-16 | 2020-06-09 | 中国科学院空间应用工程与技术中心 | 基于机器学习的高轨道卫星光压模型辨识方法及系统 |
CN111301715B (zh) * | 2020-01-17 | 2022-02-25 | 哈尔滨工业大学 | 基于霍曼变轨的同轨道特定相位分布的星座布局与轨道调整方法、装置及计算机存储介质 |
CN111301715A (zh) * | 2020-01-17 | 2020-06-19 | 哈尔滨工业大学 | 基于霍曼变轨的同轨道特定相位分布的星座布局与轨道调整方法、装置及计算机存储介质 |
CN111881505B (zh) * | 2020-08-04 | 2022-06-03 | 河北工业大学 | 一种基于ga-rbf算法的既有建筑多目标优化改造决策方法 |
CN111881505A (zh) * | 2020-08-04 | 2020-11-03 | 河北工业大学 | 一种基于ga-rbf算法的既有建筑多目标优化改造决策方法 |
CN112085812B (zh) * | 2020-08-11 | 2024-04-02 | 中国科学院力学研究所 | 一种超大规模卫星星座组网通信动态节点确定方法 |
CN112085812A (zh) * | 2020-08-11 | 2020-12-15 | 中国科学院力学研究所 | 一种超大规模卫星星座组网通信动态节点确定方法 |
WO2022198748A1 (zh) * | 2021-03-23 | 2022-09-29 | 北京航空航天大学 | 一种面向araim应用的低轨卫星增强系统星座配置优化方法 |
US11742935B2 (en) | 2021-03-23 | 2023-08-29 | Beihang University | Constellation configuration optimization method of LEO satellite augmentation system for ARAIM application |
CN113591263A (zh) * | 2021-03-24 | 2021-11-02 | 长光卫星技术有限公司 | 一种基于几何解析法的广域目标搜索星座设计方法 |
CN113591263B (zh) * | 2021-03-24 | 2024-08-09 | 长光卫星技术股份有限公司 | 一种基于几何解析法的广域目标搜索星座设计方法 |
CN113011105A (zh) * | 2021-04-08 | 2021-06-22 | 东方红卫星移动通信有限公司 | 基于改进型布谷鸟算法的多目标星座模块优化方法 |
CN115017386A (zh) * | 2022-08-08 | 2022-09-06 | 中国人民解放军63921部队 | 基于聚类的观测数据与目标库根数关联方法和装置 |
CN115017386B (zh) * | 2022-08-08 | 2022-11-08 | 中国人民解放军63921部队 | 基于聚类的观测数据与目标库根数关联方法和装置 |
CN115795927B (zh) * | 2023-02-10 | 2023-04-18 | 北京英视睿达科技股份有限公司 | 基于多目标观测的卫星轨道设计方法、装置、设备及介质 |
CN115795927A (zh) * | 2023-02-10 | 2023-03-14 | 北京英视睿达科技股份有限公司 | 基于多目标观测的卫星轨道设计方法、装置、设备及介质 |
CN115994982A (zh) * | 2023-03-22 | 2023-04-21 | 中科星图测控技术股份有限公司 | 一种基于osg的巨型星座态势展示方法 |
CN115994982B (zh) * | 2023-03-22 | 2023-07-07 | 中科星图测控技术股份有限公司 | 一种基于osg的巨型星座态势展示方法 |
CN117040607B (zh) * | 2023-10-10 | 2024-03-26 | 中国人民解放军战略支援部队航天工程大学 | 一种低轨通信卫星星座的设计方法 |
CN117040607A (zh) * | 2023-10-10 | 2023-11-10 | 中国人民解放军战略支援部队航天工程大学 | 一种低轨通信卫星星座的设计方法 |
CN117521260A (zh) * | 2023-11-17 | 2024-02-06 | 北京工业大学 | 一种面向离散分布的多目标区域覆盖的卫星星座设计方法 |
CN117521260B (zh) * | 2023-11-17 | 2024-05-03 | 北京工业大学 | 一种面向离散分布的多目标区域覆盖的卫星星座设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107798187B (zh) | 2019-08-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107798187A (zh) | 一种高效卫星星座多目标优化方法 | |
Shavarani et al. | A congested capacitated multi-level fuzzy facility location problem: An efficient drone delivery system | |
Eastman et al. | Multi-criteria and multi-objective decision making for land allocation using GIS | |
CN110288118B (zh) | 一种基于无人机的空中救援响应方法 | |
Mukhamediev et al. | Multi-criteria spatial decision making supportsystem for renewable energy development in Kazakhstan | |
CN105874479A (zh) | 用于定义和预测飞机轨迹的系统和方法 | |
CN109992923B (zh) | 一种基于变分辨率成本表面的分阶段输电线路路径规划方法 | |
CN108090624B (zh) | 一种改进元胞自动机的城市生态安全模拟与预测的方法 | |
CN105279581A (zh) | 基于差分进化的geo-uav双基sar路径规划方法 | |
Borhani et al. | A multicriteria optimization for flight route networks in large-scale airlines using intelligent spatial information | |
He et al. | Deriving urban dynamic evolution rules from self-adaptive cellular automata with multi-temporal remote sensing images | |
Hong et al. | UAV search-and-rescue planning using an adaptive memetic algorithm | |
Haklı et al. | Genetic algorithm supported by expert system to solve land redistribution problem | |
CN112766800A (zh) | 一种基于模糊层次分析法的城市绿色交通系统评价方法 | |
CN115081343B (zh) | 一种基于神经网络结合遗传算法的天基无源探测定轨方法 | |
Wu et al. | Dynamic regional splitting planning of remote sensing satellite swarm using parallel genetic PSO algorithm | |
Liu et al. | Path planning of mobile robots based on improved Genetic algorithm | |
Zeng et al. | Robotic global path-planning based modified genetic algorithm and A* algorithm | |
Rozek et al. | Multi-objective optimisation of NRHO-LLO orbit transfer via surrogate-assisted evolutionary algorithms | |
Xue et al. | High-capacity tube network design using the Hough transform | |
Xu et al. | A three-dimensional future land use simulation (FLUS-3D) model for simulating the 3D urban dynamics under the shared socio-economic pathways | |
CN113093246A (zh) | 地面多目标点成像快速判定及任务参数计算方法 | |
Gentile et al. | Structured-chromosome GA optimisation for satellite tracking | |
Tian et al. | Estimation model of global ionospheric irregularities: an artificial intelligence approach | |
Tang | Optimal scheduling method of transport path in coastal port international logistics park |
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 |