CN110471432A - 一种卫星编队构型的方法、装置及存储介质 - Google Patents
一种卫星编队构型的方法、装置及存储介质 Download PDFInfo
- Publication number
- CN110471432A CN110471432A CN201910600674.8A CN201910600674A CN110471432A CN 110471432 A CN110471432 A CN 110471432A CN 201910600674 A CN201910600674 A CN 201910600674A CN 110471432 A CN110471432 A CN 110471432A
- Authority
- CN
- China
- Prior art keywords
- satellite
- coordinate system
- axis
- star
- main
- 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
Landscapes
- Navigation (AREA)
Abstract
本发明提供一种卫星编队构型的方法，包括：获取主卫星与参考卫星的第一基线范围数据、所述参考卫星与辅卫星的第二基线范围数据、以及所述主卫星的轨道根数；根据所述主卫星的轨道根数和所述第一基线范围，确定所述参考卫星的轨道根数；根据所述参考卫星的轨道根数和所述第二基线范围，确定所述辅卫星的轨道根数；通过两次轨道根数的确定实现主卫星和辅卫星的混合编队构型。本发明还公开了一种卫星编队构型的装置及计算机存储介质。
Description
技术领域
本发明涉及卫星控制技术，尤其涉及一种卫星编队构型的方法、装置及存储介质。
背景技术
分布式星载合成孔径雷达具有全天时、全天候的特点，可以实现数字波束形成、多角度成像、极化干涉测量、土地反演和动目标检测等功能。由于星间基线会对分布式干涉合成孔径雷达(Synthetic Aperture Radar，SAR)的广泛应用产生一定的影响，为了提高分布式SAR的系统性能，国内外诸多科研工作者在其轨道设计方面展开了深入研究，提出了一些经典的卫星编队构型，推动了分布式星载SAR的高速发展。
经典的分布式星载SAR编队构型主要包括：车轮(Cartwheel)构型、钟摆(Pendulum)构型、火车(Train)构型和螺旋(Helix)构型。这些现有的卫星编队构型中，存在可扩展性差、轨道控制难度高和有效基线覆盖范围有限等问题亟待解决。
发明内容
有鉴于此，本发明实施例提供了一种探测器地月转移轨道修正的方法、装置及存储介质。
本发明实施例提供一种卫星编队构型的方法，包括：
获取主卫星与参考卫星的第一基线范围数据、所述参考卫星与辅卫星的第二基线范围数据、以及所述主卫星的轨道根数；
根据所述主卫星的轨道根数，确定所述主卫星在轨道平面坐标系的位置；
根据所述主卫星在所述轨道平面坐标系的位置，确定所述主卫星的第一卫星平台坐标系；
根据所述第一基线范围数据，确定所述参考卫星在所述第一卫星平台坐标系的位置；
对所述参考卫星在所述第一卫星平台坐标系的位置进行空间坐标系转换，得到所述参考卫星在转动地心坐标系中的位置；
根据所述参考卫星在所述转动地心坐标系中的位置，确定所述参考卫星的轨道根数；
根据所述参考卫星的轨道根数，确定所述参考卫星在所述轨道平面坐标系的位置；
根据所述参考卫星在所述轨道平面坐标系的位置，确定所述参考卫星的第二卫星平台坐标系；
根据所述第二基线范围数据，确定所述辅卫星在所述第二卫星平台坐标系的位置；
对所述辅卫星在所述第二卫星平台坐标系的位置进行空间坐标系转换，得到所述参考卫星在所述转动地心坐标系中的位置；
根据所述辅卫星在所述转动地心坐标系中的位置，确定所述辅卫星的轨道根数。
上述方案中，所述轨道根数，包括：
半长轴、偏心率、轨道倾角、升交点赤经、近地点幅角、以及平近地点角。
上述方案中，所述根据所述主卫星的轨道根数，确定所述主卫星在轨道平面坐标系的位置，包括：
所述主卫星在轨道平面坐标系的位置表示为：
其中，xVb为在所述轨道平面坐标系x轴的坐标值，yVb为在所述轨道平面坐标系y轴的坐标值，zVb为在所述轨道平面坐标系z轴的坐标值，θb为所述主卫星的真近点角，rb为所述主卫星至地心的距离；
θb根据公式求得，rb根据公式rb＝ab(1-eb×cosEb)求得，Eb为所述主卫星的偏近点角，Eb根据公式Eb＝Mb+eb×sinEb求得，Mb为所述主卫星的平近点角，eb为所述主卫星的轨道偏心率。
上述方案中，所述根据所述主卫星在所述轨道平面坐标系的位置，确定所述主卫星的第一卫星平台坐标系，包括：
所述主卫星在所述主卫星的第一卫星平台坐标系的位置根据公式求得，其中，AVbr1为从到的转置矩阵，
其中，为绕z轴逆时针旋转角度，γb为计算过程量。
上述方案中，所述根据所述第一基线范围数据，确定所述参考卫星在所述第一卫星平台坐标系的位置，包括：
所述参考卫星在所述第一卫星平台坐标系的位置表示为：
其中，xr1c为在地球半径外方向x轴的坐标值，yr1c为在卫星运动的切线方向y轴的坐标值，zr1c为在垂直于轨道平面方向z轴的坐标值，三个坐标轴构成笛卡尔坐标系；Abc是所述参考卫星沿x轴的振幅，Bbc是所述参考卫星沿z轴的振幅，Δyr1c为所述参考卫星在所述第一卫星平台坐标系的位移，αc为所述参考卫星在轨道平面的初始位置所对应的相位，βc为所述参考卫星在垂直于轨道平面的初始位置所对应的相位，T0是一个轨道周期的时长。
上述方案中，所述根据所述参考卫星的轨道根数，确定所述参考卫星在所述轨道平面坐标系的位置，包括：
所述参考卫星在轨道平面坐标系的位置表示为：
其中，xVc为在所述轨道平面坐标系x轴的坐标值，yVc为在所述轨道平面坐标系y轴的坐标值，zVc为在所述轨道平面坐标系z轴的坐标值，θc为所述参考卫星的真近点角，rc为所述参考卫星至地心的距离；
θc根据公式求得，rc根据公式rc＝ac(1-ec×cosEc)求得，Ec为所述参考卫星的偏近点角，Ec根据公式Ec＝Mc+ec×sinEc求得，Mc为所述参考卫星的平近点角，ec为所述参考卫星的轨道偏心率。
上述方案中，所述根据所述参考卫星在所述轨道平面坐标系的位置，确定所述参考卫星的第二卫星平台坐标系，包括：
所述参考卫星在所述参考卫星的第二卫星平台坐标系的位置根据公式求得，其中，AVcr2为从到的转置矩阵，
其中，为绕z轴逆时针旋转角度，γc为计算过程量。
上述方案中，所述根据所述第二基线范围数据，确定所述辅卫星在所述第二卫星平台坐标系的位置，包括：
所述辅卫星在所述第二卫星平台坐标系的位置表示为：
其中，xr2d为在地球半径外方向x轴的坐标值，yr2d为在卫星运动的切线方向y轴的坐标值，zr2d为在垂直于轨道平面方向z轴的坐标值，三个坐标轴构成笛卡尔坐标系；Acd是所述辅卫星沿x轴的振幅，Bcd是所述辅卫星沿z轴的振幅，Δyr2d为所述辅卫星在所述第二卫星平台坐标系的位移，αd为所述辅卫星在轨道平面的初始位置所对应的相位，βd为所述辅卫星在垂直于轨道平面的初始位置所对应的相位。
上述方案中，所述辅卫星的数量为多个。
本发明实施例提供一种卫星编队构型的装置，包括：
主控模块，用于获取主卫星与参考卫星的第一基线范围数据、所述参考卫星与辅卫星的第二基线范围数据、以及所述主卫星的轨道根数；
第一运算模块，用于根据所述主卫星的轨道根数，确定所述主卫星在轨道平面坐标系的位置；根据所述主卫星在所述轨道平面坐标系的位置，确定所述主卫星的第一卫星平台坐标系；根据所述第一基线范围数据，确定所述参考卫星在所述第一卫星平台坐标系的位置；对所述参考卫星在所述第一卫星平台坐标系的位置进行空间坐标系转换，得到所述参考卫星在转动地心坐标系中的位置；根据所述参考卫星在所述转动地心坐标系中的位置，确定所述参考卫星的轨道根数；
第二运算模块，用于根据所述参考卫星的轨道根数，确定所述参考卫星在所述轨道平面坐标系的位置；根据所述参考卫星在所述轨道平面坐标系的位置，确定所述参考卫星的第二卫星平台坐标系；根据所述第二基线范围数据，确定所述辅卫星在所述第二卫星平台坐标系的位置；对所述辅卫星在所述第二卫星平台坐标系的位置进行空间坐标系转换，得到所述参考卫星在所述转动地心坐标系中的位置；根据所述辅卫星在所述转动地心坐标系中的位置，确定所述辅卫星的轨道根数。
上述方案中，所述轨道根数，包括：
半长轴、偏心率、轨道倾角、升交点赤经、近地点幅角、以及平近地点角。
上述方案中，所述第一运算模块，还用于：
所述主卫星在轨道平面坐标系的位置表示为：
其中，xVb为在所述轨道平面坐标系x轴的坐标值，yVb为在所述轨道平面坐标系y轴的坐标值，zVb为在所述轨道平面坐标系z轴的坐标值，θb为所述主卫星的真近点角，rb为所述主卫星至地心的距离；
θb根据公式求得，rb根据公式rb＝ab(1-eb×cosEb)求得，Eb为所述主卫星的偏近点角，Eb根据公式Eb＝Mb+eb×sinEb求得，Mb为所述主卫星的平近点角，eb为所述主卫星的轨道偏心率。
上述方案中，所述第一运算模块，还用于：
所述主卫星在所述主卫星的第一卫星平台坐标系的位置根据公式求得，其中，AVbr1为从到的转置矩阵，
其中，为绕z轴逆时针旋转角度，γb为计算过程量。
上述方案中，所述第一运算模块，还用于：
所述参考卫星在所述第一卫星平台坐标系的位置表示为：
其中，xr1c为在地球半径外方向x轴的坐标值，yr1c为在卫星运动的切线方向y轴的坐标值，zr1c为在垂直于轨道平面方向z轴的坐标值，三个坐标轴构成笛卡尔坐标系；Abc是所述参考卫星沿x轴的振幅，Bbc是所述参考卫星沿z轴的振幅，Δyr1c为所述参考卫星在所述第一卫星平台坐标系的位移，αc为所述参考卫星在轨道平面的初始位置所对应的相位，βc为所述参考卫星在垂直于轨道平面的初始位置所对应的相位，T0是一个轨道周期的时长。
上述方案中，所述第二运算模块，还用于：
所述参考卫星在轨道平面坐标系的位置表示为：
其中，xVc为在所述轨道平面坐标系x轴的坐标值，yVc为在所述轨道平面坐标系y轴的坐标值，zVc为在所述轨道平面坐标系z轴的坐标值，θc为所述参考卫星的真近点角，rc为所述参考卫星至地心的距离；
θc根据公式求得，rc根据公式rc＝ac(1-ec×cosEc)求得，Ec为所述参考卫星的偏近点角，Ec根据公式Ec＝Mc+ec×sinEc求得，Mc为所述参考卫星的平近点角，ec为所述参考卫星的轨道偏心率。
上述方案中，所述第二运算模块，还用于：
所述参考卫星在所述参考卫星的第二卫星平台坐标系的位置根据公式求得，其中，AVcr2为从到的转置矩阵，
其中，为绕z轴逆时针旋转角度，γc为计算过程量。
上述方案中，所述第二运算模块，还用于：
所述辅卫星在所述第二卫星平台坐标系的位置表示为：
其中，xr2d为在地球半径外方向x轴的坐标值，yr2d为在卫星运动的切线方向y轴的坐标值，zr2d为在垂直于轨道平面方向z轴的坐标值，三个坐标轴构成笛卡尔坐标系；Acd是所述辅卫星沿x轴的振幅，Bcd是所述辅卫星沿z轴的振幅，Δyr2d为所述辅卫星在所述第二卫星平台坐标系的位移，αd为所述辅卫星在轨道平面的初始位置所对应的相位，βd为所述辅卫星在垂直于轨道平面的初始位置所对应的相位。
上述方案中，所述辅卫星的数量为多个。
本发明实施例提供一种计算机存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述实施例任一种卫星编队构型的方法的步骤。
本发明实施例提供的卫星编队构型的方法，通过获取主卫星与参考卫星的第一基线范围数据、所述参考卫星与辅卫星的第二基线范围数据、以及所述主卫星的轨道根数，根据所述主卫星的轨道根数和所述第一基线范围，确定所述参考卫星的轨道根数；根据所述参考卫星的轨道根数和所述第二基线范围，确定所述辅卫星的轨道根数；实现了对主卫星、参考卫星和辅卫星的混合编队构型，使主卫星和辅卫星构成的编队构型同时具备了主卫星和参考卫星、以及参考卫星和辅卫星之间的编队构型的特点；在进行轨道优化时，只需先对主卫星和参考卫星的构型进行优化，进而对参考卫星和辅卫星的构型进行优化，降低了基线优化难度和轨道控制难度；通过调整辅卫星的数量，能够增加有效基线的数目，进而扩大有效基线范围，增强可扩展性。
附图说明
附图以示例而非限制的方式大体示出了本文中所讨论的各个实施例。
图1为本发明实施例一种卫星编队构型的方法的流程示意图；
图2为本发明实施例轨道根数示意图；
图3为本发明实施例所采用分布式SAR轨道仿真平台界面图；
图4为本发明实施例中包含三个辅卫星的Cartwheel构型示意图；
图5为本发明实施例中包含三个辅卫星的Pendulum构型示意图；
图6为本发明实施例ST构型示意图；
图7为本发明实施例中四种编队构型在半个轨道周期内仿真结果；
图8为本发明实施例中四种编队构型的有效基线覆盖率对比图；
图9为本发明实施例一种卫星编队构型的装置结构示意图；
图10为本发明实施例提供的电子设备的结构示意图。
具体实施方式
为了能够更加详尽地了解本发明实施例的特点与技术内容，下面结合附图对本发明实施例的实现进行详细阐述，所附附图仅供参考说明之用，并非用来限定本发明实施例。
图1为本发明实施例一种卫星编队构型的方法的流程示意图，
步骤101，获取主卫星与参考卫星的第一基线范围数据、所述参考卫星与辅卫星的第二基线范围数据、以及所述主卫星的轨道根数,。
轨道根数，也称为轨道要素或轨道参数，是用来描述天体在其轨道运行状态的一组参数。所述轨道根数包括：半长轴、偏心率、轨道倾角、升交点赤经、近地点幅角、以及平近地点角。图2为轨道根数的示意图，如图2所示，具体地：
半长轴a，即轨道半长轴a的具体数学含义可用远地点到地心的距离Rapogee和近地点到地心的距离Rperigee来表示，即：
a＝(Rapogee+Rperigee)/2
Rapogee＝R+ha
Rperigee＝R+hp
其中，R为地球半径，ha为远地点的高度，hp为近地点的高度。
轨道的周期T，它表示卫星沿轨道运行一圈的时间，与半长轴成比例，可表示为
其中，G为引力常数，Mearth为地球质量。
偏心率e，即轨道偏心率e的具体数学含义可以用远地点的高度ha，近地点的高度hp和地球半径R来表示，即：
轨道倾角i，轨道倾角i描述了轨道相对于赤道面的倾斜度。根据度数的不同，可以将卫星轨道分为赤道轨道、顺行轨道、逆行轨道等类型。
升交点赤经Ω，升交点赤经Ω是行星轨道升交点(升交点是轨道从南到北经过赤道的点)的黄道经度，自X轴(春分点)方向在赤道平面内沿逆时针方向度量到升交点的地心夹角，它的变化范围是0°到360°。近地点幅角是从升交点沿行星运动轨道逆时针度量到近地点的角度。轨道的具体方位可以由轨道倾角、升交点赤经和近地点幅角三个根数来唯一确定。
近地点幅角ω，为从升交点沿行星运动轨道逆时针量到近地点的角度，即自天体轨道升交点在天体轨道平面内沿天体运动方向度量到近地点的角度，从升交点顺轨道运行方向度量，范围为0-360°。
真近点角θ/平近地点角M，卫星的瞬时位置通常由真近点角θ来描述，真近点角θ是某个时刻卫星位置与近地点之间的地心角，以卫星的运动方向为正。在六个轨道根数中，不考虑摄动的情况下，只有真近点角θ是受时间影响的。但是在工程应用中出于问题简化、使用方便的原因，真近点角θ通常用平近地点角M代替，平近地点角M是一个与真卫星轨道周期相同的假想卫星在轨道的外切圆上作匀速运动，在某一时刻假想卫星与近地点的地心张角。
基线范围，包括顺轨基线范围和交轨基线范围；以主卫星和辅卫星为例说明，顺轨基线范围为主卫星和辅卫星之间沿卫星轨道运行方向距离的范围，交轨基线范围为主卫星轨道和辅卫星轨道之间距离的范围。
在一些实施例中，参考卫星为卫星编队过程中虚拟的参照物，在进行轨道优化时，分别对主卫星和虚拟的参考卫星、虚拟的参考卫星和辅卫星进行优化，从而实现对整个卫星编队构型的优化，降低了优化难度。
在一些实施例中，所述辅卫星的数量为多个。通过主卫星、参考卫星和多个辅卫星的编队构型，可以使主卫星与辅卫星之间的编队构型同时具备主卫星和参考卫星的编队构型、以及参考卫星和辅卫星的编队构型的特性，提高卫星编队构型的可扩展性。
步骤102、根据所述主卫星的轨道根数，确定所述主卫星在轨道平面坐标系的位置。
在一些实施例中，所述根据所述主卫星的轨道根数，确定所述主卫星在轨道平面坐标系的位置，包括：
所述主卫星在轨道平面坐标系的位置表示为：
其中，xVb为在所述轨道平面坐标系x轴的坐标值，yVb为在所述轨道平面坐标系y轴的坐标值，zVb为在所述轨道平面坐标系z轴的坐标值，θb为所述主卫星的真近点角，rb为所述主卫星至地心的距离；
θb根据公式求得，rb根据公式rb＝ab(1-eb×cosEb)求得，Eb为所述主卫星的偏近点角，Eb根据公式Eb＝Mb+eb×sinEb求得，Mb为所述主卫星的平近点角，eb为所述主卫星的轨道偏心率。
步骤103，根据所述主卫星在所述轨道平面坐标系的位置，确定所述主卫星的第一卫星平台坐标系。
在一些实施例中，所述根据所述主卫星在所述轨道平面坐标系的位置，确定所述主卫星的第一卫星平台坐标系，包括：
所述主卫星在所述主卫星的第一卫星平台坐标系的位置根据公式求得，其中，AVbr1为从到的转置矩阵，
其中，为绕z轴逆时针旋转角度，γb为计算过程量。
步骤104，根据所述第一基线范围数据，确定所述参考卫星在所述第一卫星平台坐标系的位置。
在一些实施例中，所述根据所述第一基线范围数据，确定所述参考卫星在所述第一卫星平台坐标系的位置，包括：
所述参考卫星在所述第一卫星平台坐标系的位置表示为：
其中，xr1c为在地球半径外方向x轴的坐标值，yr1c为在卫星运动的切线方向y轴的坐标值，zr1c为在垂直于轨道平面方向z轴的坐标值，三个坐标轴构成笛卡尔坐标系；Abc是所述参考卫星沿x轴的振幅，Bbc是所述参考卫星沿z轴的振幅，Δyr1c为所述参考卫星在所述第一卫星平台坐标系的位移，αc为所述参考卫星在轨道平面的初始位置所对应的相位，βc为所述参考卫星在垂直于轨道平面的初始位置所对应的相位，T0是一个轨道周期的时长；其中，所述第一基线范围数据中，顺轨基线范围为2Abc，交轨基线范围为Bbc。
在一些实施例中，Abc＝Bbc，βc根据编队卫星需要照射的地球的纬度范围和主卫星的平近点角和近地点幅角而确定，
假设主卫星的近地点幅角为在t＝0时的真近点角为θ1，目标所在的纬度值为Latitude1。如果要经常对Latitude1纬度的目标进行高精度的顺轨和交轨干涉，即要求在观测目标时的顺轨和交轨基线同时尽可能大。当等式被满足时，在观测区域的顺轨和交轨基线可同时实现最大化。
那么假设卫星第一次运动到Latitude1纬度(目标所在纬度)的上方时的时间为t1，此时的|Z(t1)|＝Bbc。那么有其中k∈Z。那么对k赋值使得βc∈(-π,π]。即可求得βc值，那么αc＝(βc+π/2)％(2π)，其中αc∈(-π,π]。
确定t1的过程如下：
假设卫星的地心角速度为常数轨道倾角为i。那么可得：
当Latitude1/sin i>(ωb-θ1)时，
当Latitude1/sin i<(ωb-θ1)时，
步骤105，对所述参考卫星在所述第一卫星平台坐标系的位置进行空间坐标系转换，得到所述参考卫星在转动地心坐标系中的位置。
对于空间坐标系转换，具体如下：
星载干涉合成孔径雷达(Synthetic Aperture Radar，SAR)平台空间变换常用的4个直角坐标系，分别为：转动地心坐标系、不转动地心坐标系、卫星轨道平面坐标系和卫星平台坐标系。
为卫星在转动地心坐标系的位置，为卫星在不转动地心坐标系的位置，为卫星在卫星轨道平面坐标系的位置，为卫星在卫星平台坐标系的位置；转换关系如下：
其中，Amn和Anm互为逆矩阵(Amn·Anm＝I)，m,n＝o,g,v,r,。
具体地：
从转动地心坐标系——→不转动地心坐标系
不转动地心坐标系，XY平面与地球赤道面重合。X轴指向天球(CelestialSphere，一种假想的无限大的球，它和地球同心。所以ECI坐标系不受地球旋转的影响)的某个位置。Z轴与XY平面垂直并指向北极。将绕Z轴逆时针旋转过一个春分点的格林威治时间HG就得到转动的地心坐标系HG＝ωe(t-t0)，其中t0时刻和重合，ωe代表地球自转角速度，那么到的转化矩阵为Ago为
从不转动地心坐标系——→轨道平面坐标系
轨道平面坐标系，与卫星轨道六根数有关。的原点为地球中心。XY平面也与地球赤道面重合。不过其X轴指向0经度方向，Y轴指向东经90度的方向。所以轨道平面坐标系实际上是随着地球一起旋转的。不转动地心坐标系需经三次旋转得到轨道平面坐标系。首先，将不转动地心坐标系绕Z轴逆时针旋转一个角Ω；其次，将得到的坐标系绕X轴逆时针旋转一个角度i，最后，再将得到的坐标系绕Z轴逆时针旋转一个角度ω得到卫星轨道平面坐标系那么到的转化矩阵为Aov为
其中的Ω，i，ω分别为卫星轨道六根数的升交点赤经、卫星轨道面倾角以及近地点幅角。
从轨道平面坐标系——→卫星平台坐标系
卫星平台坐标系的X轴与卫星速度方向平行，将绕Z轴逆时针旋转可以得到卫星平台坐标系是卫星姿态的参考坐标系，到的转化矩阵为Avr为
其中θ和e分别为卫星轨道六根数的真近心角和卫星的偏心率。
步骤106，根据所述参考卫星在所述转动地心坐标系中的位置，确定所述参考卫星的轨道根数。
在一些实施例中，根据所述参考卫星在所述转动地心坐标系中的位置，确定所述参考卫星的轨道根数，包括：
根据所述参考卫星在所述转动地心坐标系中的位置，确定所述参考卫星的位置矢量和速度矢量常系数μ＝GMearth＝3.986005e+14，其中，G为引力常数，Mearth为地球质量，由所述参考卫星的位置矢量和速度矢量计算所述参考卫星的轨道根数，具体如下：
所述参考卫星的半长轴ac根据公式求得；其中，rc为所述参考卫星至地心的距离
所述参考卫星的轨道倾角ic根据公式cosic＝hcz/hc求得，其中，hc为的模值，hcz为在z轴的投影，为过程向量且 为所述参考卫星至地心的距离向量。
所述参考卫星的升交点赤经Ωc根据公式求得，其中，ncx为在x轴的投影，ncy为在y轴的投影，为计算过程量，且其中，为沿z轴的单位矢量，
所述参考卫星的偏心率ec根据公式p＝hc 2/μ＝ac(1-ec 2)求得，其中，p为曲线半通径。
所述参考卫星的近地点幅角ωc根据公式ωc＝u-θc求得，其中，uc为计算过程量，
所述参考卫星的真近点角θc根据公式
求得；
所述参考卫星的平近地点角及偏近心角根据公式和Ec＝Mc+ec×sinEc求得，bc为所述参考卫星的轨道半短轴。
步骤107，根据所述参考卫星的轨道根数，确定所述参考卫星在所述轨道平面坐标系的位置。
在一些实施例中，所述根据所述参考卫星的轨道根数，确定所述参考卫星在所述轨道平面坐标系的位置，包括：
所述参考卫星在轨道平面坐标系的位置表示为：
其中，xVc为在所述轨道平面坐标系x轴的坐标值，yVc为在所述轨道平面坐标系y轴的坐标值，zVc为在所述轨道平面坐标系z轴的坐标值，θc为所述参考卫星的真近点角，rc为所述参考卫星至地心的距离；
θc根据公式求得，rc根据公式rc＝ac(1-ec×cos Ec)求得，Ec为所述参考卫星的偏近点角，Ec根据公式Ec＝Mc+ec×sin Ec求得，Mc为所述参考卫星的平近点角，ec为所述参考卫星的轨道偏心率。
步骤108，根据所述参考卫星在所述轨道平面坐标系的位置，确定所述参考卫星的第二卫星平台坐标系。
在一些实施例中，所述根据所述参考卫星在所述轨道平面坐标系的位置，确定所述参考卫星的第二卫星平台坐标系，包括：
所述参考卫星在所述参考卫星的第二卫星平台坐标系的位置根据公式求得，其中，AVcr2为从到的转置矩阵，
其中，为绕z轴逆时针旋转角度，γc为计算过程量。
步骤109，根据所述第二基线范围数据，确定所述辅卫星在所述第二卫星平台坐标系的位置。
在一些实施例中，所述根据所述第二基线范围数据，确定所述辅卫星在所述第二卫星平台坐标系的位置，包括：
所述辅卫星在所述第二卫星平台坐标系的位置表示为：
其中，xr2d为在地球半径外方向x轴的坐标值，yr2d为在卫星运动的切线方向y轴的坐标值，zr2d为在垂直于轨道平面方向z轴的坐标值，三个坐标轴构成笛卡尔坐标系；Acd是所述辅卫星沿x轴的振幅，Bcd是所述辅卫星沿z轴的振幅，Δyr2d为所述辅卫星在所述第二卫星平台坐标系的位移，αd为所述辅卫星在轨道平面的初始位置所对应的相位，βd为所述辅卫星在垂直于轨道平面的初始位置所对应的相位。
在一些实施例中，辅卫星的数量为N个，N个辅卫星中第i个辅卫星沿x轴的振幅为Ai，N为正整数，1≤i≤N，所述第i个辅卫星沿z轴的振幅为Bi，沿x轴的振幅的取值为：
A1＝A2＝…＝AN＝0
Bi和Δyi根据现有的卫星的控轨精度而定。
由于Ai为0，那么α值可以忽略。
其中i∈{2,3,…,N}
β1根据编队卫星需要照射的地球的纬度范围和参考卫星的真近点角和近地点幅角而确定。
当N为奇数时，可得
对k取值使得β1∈(-π,π]即可。
当N为偶数时，可得
(1)N/2仍为奇数
对k取值使得β1∈(-π,π]即可。
(2)N/2仍为偶数
对k取值使得β1∈(-π,π]即可。
步骤110，对所述辅卫星在所述第二卫星平台坐标系的位置进行空间坐标系转换，得到所述参考卫星在所述转动地心坐标系中的位置；
步骤111，根据所述辅卫星在所述转动地心坐标系中的位置，确定所述辅卫星的轨道根数。
在一些实施例中，根据所述辅卫星在所述转动地心坐标系中的位置，确定所述辅卫星的轨道根数，包括：
根据所述辅卫星在所述转动地心坐标系中的位置，确定所述辅卫星的位置矢量和速度矢量常系数μ＝GMearth＝3.986005e+14，其中，G为引力常数，Mearth为地球质量，由所述辅卫星的位置矢量和速度矢量计算所述辅卫星的轨道根数，具体如下：
所述辅卫星的半长轴ad根据公式求得；其中，rd为所述辅卫星至地心的距离
所述辅卫星的轨道倾角id根据公式cosid＝hdz/hd求得，其中，hd为的模值，hdz为在z轴的投影，为过程向量且 为所述辅卫星至地心的距离向量。
所述辅卫星的升交点赤经Ωd根据公式求得，其中，ndx为在x轴的投影，ndy为在y轴的投影，为计算过程量，且其中，为沿z轴的单位矢量，
所述辅卫星的偏心率ed根据公式p＝hd 2/μ＝ad(1-ed 2)求得，其中，p为曲线半通径。
所述辅卫星的近地点幅角ωd根据公式ωd＝ud-θd求得，其中，ud为计算过程量，
所述辅卫星的真近点角θd根据公式和求得；
所述辅卫星的平近地点角Md及偏近心角Ed根据公式
和Ed＝Md+ed×sinEd求得，其中，bd为所述辅卫星的轨道半短轴。
通过确定参考卫星的轨道根数和辅卫星的轨道根数，实现卫星编队构型。
在一些实施例中，根据上述实施例一种卫星编队构型的方法，通过SAR星座轨道仿真平台，获取和处理卫星构型的轨道数据，进而实现对编队构型的优化。所述SAR星座轨道仿真平台的界面如图3所示，由C#和C++语言混合编制而成，能够直观地对卫星编队进行构型和调整，具有良好的用户体验效果。
本发明实施例的提供了上述实施例一种卫星编队构型的方法的应用，一种超多面体(Super-Polyhedron，SP)构型，主卫星和一颗虚拟的参考辅卫星构成一个Helix构型，同时N颗伴随辅卫星相对于这颗虚拟的参考辅卫星构型一个Pendulum构型。在进行轨道优化时，只需先对Helix构型进行优化，进而对Pendulum进行优化，从而实现对SP构型的优化。由于Helix构型和Pendulum构型的基线优化较为容易，因此降低了SP构型的基线优化难度。
通过依次微调伴随辅卫星相对于参考主卫星的近地点幅角和平近点角(也可以直接调整虚拟的参考辅卫星与主卫星的近地点幅角和平近点角)，进而调整伴随辅卫星相对于主卫星的顺轨基线范围，将基线控制在顺轨干涉所要求的范围内，扩大顺轨基线的覆盖率。
通过依次微调伴随辅卫星相对于参考主卫星的升交点赤经和轨道倾角(也可以直接调整虚拟的参考辅卫星与主卫星的升交点赤经和轨道倾角)，可以调整伴随辅卫星相对于主卫星的交轨基线范围，将基线控制在交轨干涉所要求的可用基线范围内，扩大交轨基线的覆盖率。
可知N颗参考辅卫星的平近点角和升交点赤经略有差别，并且主卫星和N颗参考辅卫星可近似看成N对Helix构型。由于Helix构型的轨道控制难度低，故SP构型的轨道控制难度大大降低。
N颗伴随辅卫星和主卫星可近似视为N对Helix编队飞行。只需要调整待添加的伴随辅卫星相对于虚拟参考辅卫星的平近点角与升交点赤经即可实现对卫星数目的扩展，其可扩展性大大增加。
N颗伴随辅卫星和主卫星可近似视为N对Helix编队飞行。这N对Helix编队卫星均处于不同的圆筒壁上，彼此环绕飞行，提升了SP构型星座的安全性。
具体地，本发明实施例的提供了一种超四面体(Super-Tetrahedron，ST)构型，超四面体是超多面体构型的一个特例，所述ST构型的轨道根数如表1所示。
| 轨道根数 | 主卫星 | 辅卫星1 | 辅卫星2 | 辅卫星3 |
| 半长轴 | a | a | a | a |
| 偏心率 | e<sub>1</sub> | e<sub>2</sub> | e<sub>2</sub> | e<sub>2</sub> |
| 轨道倾角 | i | i | i | i |
| 升交点赤经 | Ω<sub>1</sub> | Ω<sub>2</sub> | Ω<sub>2</sub> | Ω<sub>2</sub> |
| 近地点幅角 | ω<sub>1</sub> | ω<sub>2</sub> | ω<sub>2</sub> | ω<sub>2</sub> |
| 平近点角 | M<sub>1</sub> | M<sub>2</sub> | M<sub>3</sub> | M<sub>4</sub> |
表1
在实验仿真部分，图4为包含三个伴随辅卫星和一个主卫星的Cartwheel构型示意图，图5为包含三个伴随辅卫星和一个主卫星的Pendulum构型示意图，图6为ST构型示意图，主要通过对由三个伴随辅卫星和一个主卫星组成的Cartwheel构型、Pendulum构型、四星(Four)-Helix构型和ST构型四种卫星构型在如图3所示的分布式SAR轨道仿真平台上进行仿真。
四种卫星构型在半个轨道周期内顺轨和交轨基线的覆盖区域的仿真结果如图7所示，其中a1为Cartwheel构型顺轨基线，a2为Cartwheel构型交轨基线，b1为Pendulum构型顺轨基线，b2为Pendulum构型交轨基线，c1为Four-Helix构型顺轨基线，c2为Four-Helix构型交轨基线，d1为ST构型顺轨基线，d2为ST构型交轨基线，实线表示有效基线的覆盖区域，虚线表示无效覆盖区域。可见在四种编队构型中，只有ST构型的顺轨和交轨基线可以同时满足80％以上有效区域的全球覆盖。其可以在保证顺轨基线最大覆盖的同时也可以实现有效交轨基线的较大覆盖，且其交轨基线明显大于其他三种编队构型。
为了突出显示ST构型的优势，四种卫星编队构型交轨基线的有效覆盖范围在世界地图中标注出来，四种编队构型的有效基线覆盖率对比如图8所示，四种不同的颜色代表四种不同编队构型的覆盖区域，ST构型的覆盖区域明显大于其他三种编队构型。
实验仿真结果表明，超多面体编队构型具有：有效基线覆盖比例大、可扩展性强、卫星轨道控制难度(能源消耗)小、星间基线优化难度低和星座构型安全性高等特点。
为实现上述实施例卫星编队构型的方法，本发明实施例提供一种卫星编队构型的装置，所述装置的结构如图9所示，包括：
主控模块91，用于获取主卫星与参考卫星的第一基线范围数据、所述参考卫星与辅卫星的第二基线范围数据、以及所述主卫星的轨道根数。
轨道根数，也称为轨道要素或轨道参数，是用来描述天气在其轨道运行状态的一组参数。所述轨道根数包括：半长轴、偏心率、轨道倾角、升交点赤经、近地点幅角、以及平近地点角。图2为轨道根数的示意图，如图2所示，具体地：
半长轴a，即轨道半长轴a的具体数学含义可用远地点到地心的距离Rapogee和近地点到地心的距离Rperigee来表示，即：
a＝(Rapogee+Rperigee)/2
Rapogee＝R+ha
Rperigee＝R+hp
其中，R为地球半径，ha为远地点的高度，hp为近地点的高度。
轨道的周期T，它表示卫星沿轨道运行一圈的时间，与半长轴成比例，可表示为
其中，G为引力常数，Mearth为地球质量。
偏心率e，即轨道偏心率e的具体数学含义可以用远地点的高度ha，近地点的高度hp和地球半径R来表示，即：
轨道倾角i，轨道倾角i描述了轨道相对于赤道面的倾斜度。根据度数的不同，可以将卫星轨道分为赤道轨道、顺行轨道、逆行轨道等类型。
升交点赤经Ω，升交点赤经Ω是行星轨道升交点(升交点是轨道从南到北经过赤道的点)的黄道经度，自X轴(春分点)方向在赤道平面内沿逆时针方向度量到升交点的地心夹角，它的变化范围是0°到360°。近地点幅角是从升交点沿行星运动轨道逆时针度量到近地点的角度。轨道的具体方位可以由轨道倾角、升交点赤经和近地点幅角三个根数来唯一确定。
近地点幅角ω，为从升交点沿行星运动轨道逆时针量到近地点的角度，即自天体轨道升交点在天体轨道平面内沿天体运动方向度量到近地点的角度，从升交点顺轨道运行方向度量，范围为0-360°。
真近点角θ/平近地点角M，卫星的瞬时位置通常由真近点角θ来描述，真近点角θ是某个时刻卫星位置与近地点之间的地心角，以卫星的运动方向为正。在六个轨道根数中，不考虑摄动的情况下，只有真近点角θ是受时间影响的。但是在工程应用中出于问题简化、使用方便的原因，真近点角θ通常用平近地点角M代替，平近地点角M是一个与真卫星轨道周期相同的假想卫星在轨道的外切圆上作匀速运动，在某一时刻假想卫星与近地点的地心张角。
基线范围，包括顺轨基线范围和交轨基线范围；以主卫星和辅卫星为例说明，顺轨基线范围为主卫星和辅卫星之间沿卫星轨道运行方向距离的范围，交轨基线范围为主卫星轨道和辅卫星轨道之间距离的范围。
在一些实施例中，参考卫星为卫星编队过程中虚拟的参照物，通过虚拟的参考卫星对卫星编队构型，在进行轨道优化时，分别对主卫星和参考卫星、参考卫星和辅卫星进行优化，从而实现对整个卫星编队构型的优化，降低了优化难度。
在一些实施例中，所述辅卫星的数量为多个。通过主卫星、参考卫星和多个辅卫星的编队构型，可以使主卫星与辅卫星之间的编队构型同时具备主卫星和参考卫星的编队构型、以及参考卫星和辅卫星的编队构型的特性，提高卫星编队构型的可扩展性。
第一运算模块92，用于根据所述主卫星的轨道根数，确定所述主卫星在轨道平面坐标系的位置；根据所述主卫星在所述轨道平面坐标系的位置，确定所述主卫星的第一卫星平台坐标系；根据所述第一基线范围数据，确定所述参考卫星在所述第一卫星平台坐标系的位置；对所述参考卫星在所述第一卫星平台坐标系的位置进行空间坐标系转换，得到所述参考卫星在转动地心坐标系中的位置；根据所述参考卫星在所述转动地心坐标系中的位置，确定所述参考卫星的轨道根数。
第二运算模块93，用于根据所述参考卫星的轨道根数，确定所述参考卫星在所述轨道平面坐标系的位置；根据所述参考卫星在所述轨道平面坐标系的位置，确定所述参考卫星的第二卫星平台坐标系；根据所述第二基线范围数据，确定所述辅卫星在所述第二卫星平台坐标系的位置；对所述辅卫星在所述第二卫星平台坐标系的位置进行空间坐标系转换，得到所述参考卫星在所述转动地心坐标系中的位置；根据所述辅卫星在所述转动地心坐标系中的位置，确定所述辅卫星的轨道根数。
在一些实施例中，所述第一运算模块，还用于：
所述主卫星在轨道平面坐标系的位置表示为：
其中，xVb为在所述轨道平面坐标系x轴的坐标值，yVb为在所述轨道平面坐标系y轴的坐标值，zVb为在所述轨道平面坐标系z轴的坐标值，θb为所述主卫星的真近点角，rb为所述主卫星至地心的距离；
θb根据公式求得，rb根据公式rb＝ab(1-eb×cos Eb)求得，Eb为所述主卫星的偏近点角，Eb根据公式Eb＝Mb+eb×sin Eb求得，Mb为所述主卫星的平近点角，eb为所述主卫星的轨道偏心率。
在一些实施例中，所述第一运算模块，还用于：
所述主卫星在所述主卫星的第一卫星平台坐标系的位置根据公式求得，其中，AVbr1为从到的转置矩阵，
其中，为绕z轴逆时针旋转角度，γb为计算过程量。
在一些实施例中，所述第一运算模块，还用于：
所述参考卫星在所述第一卫星平台坐标系的位置表示为：
其中，xr1c为在地球半径外方向x轴的坐标值，yr1c为在卫星运动的切线方向y轴的坐标值，zr1c为在垂直于轨道平面方向z轴的坐标值，三个坐标轴构成笛卡尔坐标系；Abc是所述参考卫星沿x轴的振幅，Bbc是所述参考卫星沿z轴的振幅，Δyr1c为所述参考卫星在所述第一卫星平台坐标系的位移，αc为所述参考卫星在轨道平面的初始位置所对应的相位，βc为所述参考卫星在垂直于轨道平面的初始位置所对应的相位，T0是一个轨道周期的时长；其中，所述第一基线范围数据中，顺轨基线范围为2Abc，交轨基线范围为Bbc。
在一些实施例中，Abc＝Bbc，βc根据编队卫星需要照射的地球的纬度范围和主卫星的平近点角和近地点幅角而确定，
假设主卫星的近地点幅角为在t＝0时的真近点角为θ1，目标所在的纬度值为Latitude1。如果要经常对Latitude1纬度的目标进行高精度的顺轨和交轨干涉，即要求在观测目标时的顺轨和交轨基线同时尽可能大。当等式被满足时，在观测区域的顺轨和交轨基线可同时实现最大化。
那么假设卫星第一次运动到Latitude1纬度(目标所在纬度)的上方时的时间为t1，此时的|Z(t1)|＝Bbc。那么有其中k∈Z。那么对k赋值使得βc∈(-π,π]。即可求得βc值，那么αc＝(βc+π/2)％(2π)，其中αc∈(-π,π]。
确定t1的过程如下：
假设卫星的地心角速度为常数轨道倾角为i。那么可得：
当Latitude1/sin i>(ωb-θ1)时，
当Latitude1/sin i<(ωb-θ1)时，
在一些实施例中，所述第二运算模块，还用于：
所述参考卫星在轨道平面坐标系的位置表示为：
其中，xVc为在所述轨道平面坐标系x轴的坐标值，yVc为在所述轨道平面坐标系y轴的坐标值，zVc为在所述轨道平面坐标系z轴的坐标值，θc为所述参考卫星的真近点角，rc为所述参考卫星至地心的距离；
θc根据公式求得，rc根据公式rc＝ac(1-ec×cos Ec)求得，Ec为所述参考卫星的偏近点角，Ec根据公式Ec＝Mc+ec×sin Ec求得，Mc为所述参考卫星的平近点角，ec为所述参考卫星的轨道偏心率。
在一些实施例中，所述第二运算模块，还用于：
所述参考卫星在所述参考卫星的第二卫星平台坐标系的位置根据公式求得，其中，AVcr2为从到的转置矩阵，
其中，为绕z轴逆时针旋转角度，γc为计算过程量。
在一些实施例中，所述第二运算模块，还用于：
所述辅卫星在所述第二卫星平台坐标系的位置表示为：
其中，xr2d为在地球半径外方向x轴的坐标值，yr2d为在卫星运动的切线方向y轴的坐标值，zr2d为在垂直于轨道平面方向z轴的坐标值，三个坐标轴构成笛卡尔坐标系；Acd是所述辅卫星沿x轴的振幅，Bcd是所述辅卫星沿z轴的振幅，Δyr2d为所述辅卫星在所述第二卫星平台坐标系的位移，αd为所述辅卫星在轨道平面的初始位置所对应的相位，βd为所述辅卫星在垂直于轨道平面的初始位置所对应的相位。
在一些实施例中，辅卫星的数量为N个，N个辅卫星中第i个辅卫星沿x轴的振幅为Ai，N为正整数，1≤i≤N，所述第i个辅卫星沿z轴的振幅为Bi，沿x轴的振幅的取值为：
A1＝A2＝…＝AN＝0
Bi和Δyi根据现有的卫星的控轨精度而定。
由于Ai为0，那么α值可以忽略。
其中i∈{2,3,…,N}
β1根据编队卫星需要照射的地球的纬度范围和参考卫星的真近点角和近地点幅角而确定。
当N为奇数时，可得
对k取值使得β1∈(-π,π]即可。
当N为偶数时，可得
(1)N/2仍为奇数
对k取值使得β1∈(-π,π]即可。
(2)N/2仍为偶数
对k取值使得β1∈(-π,π]即可。
为了实现本发明实施例的卫星编队构型的方法，本发明实施例提供了一种如图10所示结构示意图的的电子设备，如图10所示，本发明实施例的电子设备610包括：处理器61和用于存储能够在处理器上运行的计算机程序的存储器62，其中，
所述处理器61用于运行所述计算机程序时，执行本发明实施例任一种卫星编队构型的方法的步骤。
当然，实际应用时，如图6所示，该电子设备还可以包括至少一个通信接口63。电子设备中的各个组件通过总线系统64耦合在一起。可理解，总线系统64用于实现这些组件之间的连接通信。总线系统64除包括数据总线之外，还包括电源总线、控制总线和状态信号总线。但是为了清楚说明起见，在图6中将各种总线都标为总线系统64。
其中，通信接口63，用于与其它设备进行交互。
具体来说，所述处理器61可以通过通信接口63向对应所述被调用方应用的应用服务器发送操作结果查询请求，获取所述应用服务器发送的所述被调用方应用的操作结果。
本领域技术人员可以理解，存储器62可以是易失性存储器或非易失性存储器，也可包括易失性和非易失性存储器两者。其中，非易失性存储器可以是只读存储器(ROM，ReadOnly Memory)、可编程只读存储器(PROM，Programmable Read-Only Memory)、可擦除可编程只读存储器(EPROM，Erasable Programmable Read-Only Memory)、电可擦除可编程只读存储器(EEPROM，Electrically Erasable Programmable Read-Only Memory)、磁性随机存取存储器(FRAM，ferromagnetic random access memory)、快闪存储器(Flash Memory)、磁表面存储器、光盘、或只读光盘(CD-ROM，Compact Disc Read-Only Memory)；磁表面存储器可以是磁盘存储器或磁带存储器。易失性存储器可以是随机存取存储器(RAM，RandomAccess Memory)，其用作外部高速缓存。通过示例性但不是限制性说明，许多形式的RAM可用，例如静态随机存取存储器(SRAM，Static Random Access Memory)、同步静态随机存取存储器(SSRAM，Synchronous Static Random Access Memory)、动态随机存取存储器(DRAM，Dynamic Random Access Memory)、同步动态随机存取存储器(SDRAM，SynchronousDynamic Random Access Memory)、双倍数据速率同步动态随机存取存储器(DDRSDRAM，Double Data Rate Synchronous Dynamic Random Access Memory)、增强型同步动态随机存取存储器(ESDRAM，Enhanced Synchronous Dynamic Random Access Memory)、同步连接动态随机存取存储器(SLDRAM，SyncLink Dynamic Random Access Memory)、直接内存总线随机存取存储器(DRRAM，Direct Rambus Random Access Memory)。本发明实施例描述的存储器62旨在包括但不限于这些和任意其它适合类型的存储器。
本发明实施例中，还提供了一种计算机可读存储介质，用于存储上述实施例中提供的计算程序，以完成前述卫星编队构型的方法所述步骤。计算机可读存储介质可以是FRAM、ROM、PROM、EPROM、EEPROM、Flash Memory、磁表面存储器、光盘、或CD-ROM等存储器；也可以是包括上述存储器之一或任意组合的各种设备，如移动电话、计算机、智能家电、服务器等。
需要说明的是：本发明实施例所记载的技术方案之间，在不冲突的情况下，可以任意组合。
本申请所提供的几个方法实施例中所揭露的方法，在不冲突的情况下可以任意组合，得到新的方法实施例。
本申请所提供的几个产品实施例中所揭露的特征，在不冲突的情况下可以任意组合，得到新的产品实施例。
本申请所提供的几个方法或装置实施例中所揭露的特征，在不冲突的情况下可以任意组合，得到新的方法实施例或装置实施例。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。
Claims (19)
1.一种卫星编队构型的方法，其特征在于，所述方法包括：
获取主卫星与参考卫星的第一基线范围数据、所述参考卫星与辅卫星的第二基线范围数据、以及所述主卫星的轨道根数；
根据所述主卫星的轨道根数，确定所述主卫星在轨道平面坐标系的位置；
根据所述主卫星在所述轨道平面坐标系的位置，确定所述主卫星的第一卫星平台坐标系；
根据所述第一基线范围数据，确定所述参考卫星在所述第一卫星平台坐标系的位置；
对所述参考卫星在所述第一卫星平台坐标系的位置进行空间坐标系转换，得到所述参考卫星在转动地心坐标系中的位置；
根据所述参考卫星在所述转动地心坐标系中的位置，确定所述参考卫星的轨道根数；
根据所述参考卫星的轨道根数，确定所述参考卫星在所述轨道平面坐标系的位置；
根据所述参考卫星在所述轨道平面坐标系的位置，确定所述参考卫星的第二卫星平台坐标系；
根据所述第二基线范围数据，确定所述辅卫星在所述第二卫星平台坐标系的位置；
对所述辅卫星在所述第二卫星平台坐标系的位置进行空间坐标系转换，得到所述参考卫星在所述转动地心坐标系中的位置；
根据所述辅卫星在所述转动地心坐标系中的位置，确定所述辅卫星的轨道根数。
2.根据权利要求1所述的方法，其特征在于，所述轨道根数，包括：
半长轴、偏心率、轨道倾角、升交点赤经、近地点幅角、以及平近地点角。
3.根据权利要求2所述的方法，其特征在于，所述根据所述主卫星的轨道根数，确定所述主卫星在轨道平面坐标系的位置，包括：
所述主卫星在轨道平面坐标系的位置表示为：
其中，xVb为在所述轨道平面坐标系x轴的坐标值，yVb为在所述轨道平面坐标系y轴的坐标值，zVb为在所述轨道平面坐标系z轴的坐标值，θb为所述主卫星的真近点角，rb为所述主卫星至地心的距离；
θb根据公式求得，rb根据公式rb＝ab(1-eb×cosEb)求得，Eb为所述主卫星的偏近点角，Eb根据公式Eb＝Mb+eb×sinEb求得，Mb为所述主卫星的平近点角，eb为所述主卫星的轨道偏心率。
4.根据权利要求3所述的方法，其特征在于，所述根据所述主卫星在所述轨道平面坐标系的位置，确定所述主卫星的第一卫星平台坐标系，包括：
所述主卫星在所述主卫星的第一卫星平台坐标系的位置根据公式求得，其中，AVbr1为从到的转置矩阵，
其中，为绕z轴逆时针旋转角度，γb为计算过程量。
5.根据权利要求4所述的方法，其特征在于，所述根据所述第一基线范围数据，确定所述参考卫星在所述第一卫星平台坐标系的位置，包括：
所述参考卫星在所述第一卫星平台坐标系的位置表示为：
其中，xr1c为在地球半径外方向x轴的坐标值，yr1c为在卫星运动的切线方向y轴的坐标值，zr1c为在垂直于轨道平面方向z轴的坐标值，三个坐标轴构成笛卡尔坐标系；Abc是所述参考卫星沿x轴的振幅，Bbc是所述参考卫星沿z轴的振幅，Δyr1c为所述参考卫星在所述第一卫星平台坐标系的位移，αc为所述参考卫星在轨道平面的初始位置所对应的相位，βc为所述参考卫星在垂直于轨道平面的初始位置所对应的相位，T0是一个轨道周期的时长。
6.根据权利要求5所述的方法，其特征在于，所述根据所述参考卫星的轨道根数，确定所述参考卫星在所述轨道平面坐标系的位置，包括：
所述参考卫星在轨道平面坐标系的位置表示为：
其中，xVc为在所述轨道平面坐标系x轴的坐标值，yVc为在所述轨道平面坐标系y轴的坐标值，zVc为在所述轨道平面坐标系z轴的坐标值，θc为所述参考卫星的真近点角，rc为所述参考卫星至地心的距离；
θc根据公式求得，rc根据公式rc＝ac(1-ec×cosEc)求得，Ec为所述参考卫星的偏近点角，Ec根据公式Ec＝Mc+ec×sinEc求得，Mc为所述参考卫星的平近点角，ec为所述参考卫星的轨道偏心率。
7.根据权利要求6所述的方法，其特征在于，所述根据所述参考卫星在所述轨道平面坐标系的位置，确定所述参考卫星的第二卫星平台坐标系，包括：
所述参考卫星在所述参考卫星的第二卫星平台坐标系的位置根据公式求得，其中，AVcr2为从到的转置矩阵，
其中，为绕z轴逆时针旋转角度，γc为计算过程量。
8.根据权利要求7所述的方法，其特征在于，所述根据所述第二基线范围数据，确定所述辅卫星在所述第二卫星平台坐标系的位置，包括：
所述辅卫星在所述第二卫星平台坐标系的位置表示为：
其中，xr2d为在地球半径外方向x轴的坐标值，yr2d为在卫星运动的切线方向y轴的坐标值，zr2d为在垂直于轨道平面方向z轴的坐标值，三个坐标轴构成笛卡尔坐标系；Acd是所述辅卫星沿x轴的振幅，Bcd是所述辅卫星沿z轴的振幅，Δyr2d为所述辅卫星在所述第二卫星平台坐标系的位移，αd为所述辅卫星在轨道平面的初始位置所对应的相位，βd为所述辅卫星在垂直于轨道平面的初始位置所对应的相位。
9.根据权利要求8所述的方法，其特征在于，所述辅卫星的数量为多个。
10.一种卫星编队构型的装置，其特征在于，所述装置包括：
主控模块，用于获取主卫星与参考卫星的第一基线范围数据、所述参考卫星与辅卫星的第二基线范围数据、以及所述主卫星的轨道根数；
第一运算模块，用于根据所述主卫星的轨道根数，确定所述主卫星在轨道平面坐标系的位置；根据所述主卫星在所述轨道平面坐标系的位置，确定所述主卫星的第一卫星平台坐标系；根据所述第一基线范围数据，确定所述参考卫星在所述第一卫星平台坐标系的位置；对所述参考卫星在所述第一卫星平台坐标系的位置进行空间坐标系转换，得到所述参考卫星在转动地心坐标系中的位置；根据所述参考卫星在所述转动地心坐标系中的位置，确定所述参考卫星的轨道根数；
第二运算模块，用于根据所述参考卫星的轨道根数，确定所述参考卫星在所述轨道平面坐标系的位置；根据所述参考卫星在所述轨道平面坐标系的位置，确定所述参考卫星的第二卫星平台坐标系；根据所述第二基线范围数据，确定所述辅卫星在所述第二卫星平台坐标系的位置；对所述辅卫星在所述第二卫星平台坐标系的位置进行空间坐标系转换，得到所述参考卫星在所述转动地心坐标系中的位置；根据所述辅卫星在所述转动地心坐标系中的位置，确定所述辅卫星的轨道根数。
11.根据权利要求10所述的装置，其特征在于，所述轨道根数，包括：
半长轴、偏心率、轨道倾角、升交点赤经、近地点幅角、以及平近地点角。
12.根据权利要求11所述的装置，其特征在于，所述第一运算模块，还用于：
所述主卫星在轨道平面坐标系的位置表示为：
其中，xVb为在所述轨道平面坐标系x轴的坐标值，yVb为在所述轨道平面坐标系y轴的坐标值，zVb为在所述轨道平面坐标系z轴的坐标值，θb为所述主卫星的真近点角，rb为所述主卫星至地心的距离；
θb根据公式求得，rb根据公式rb＝ab(1-eb×cos Eb)求得，Eb为所述主卫星的偏近点角，Eb根据公式Eb＝Mb+eb×sinEb求得，Mb为所述主卫星的平近点角，eb为所述主卫星的轨道偏心率。
13.根据权利要求12所述的装置，其特征在于，所述第一运算模块，还用于：
所述主卫星在所述主卫星的第一卫星平台坐标系的位置根据公式求得，其中，AVbr1为从到的转置矩阵，
其中，为绕z轴逆时针旋转角度，γb为计算过程量。
14.根据权利要求13所述的装置，其特征在于，所述第一运算模块，还用于：
所述参考卫星在所述第一卫星平台坐标系的位置表示为：
其中，xr1c为在地球半径外方向x轴的坐标值，yr1c为在卫星运动的切线方向y轴的坐标值，zr1c为在垂直于轨道平面方向z轴的坐标值，三个坐标轴构成笛卡尔坐标系；Abc是所述参考卫星沿x轴的振幅，Bbc是所述参考卫星沿z轴的振幅，Δyr1c为所述参考卫星在所述第一卫星平台坐标系的位移，αc为所述参考卫星在轨道平面的初始位置所对应的相位，βc为所述参考卫星在垂直于轨道平面的初始位置所对应的相位，T0是一个轨道周期的时长。
15.根据权利要求14所述的装置，其特征在于，所述第二运算模块，还用于：
所述参考卫星在轨道平面坐标系的位置表示为：
其中，xVc为在所述轨道平面坐标系x轴的坐标值，yVc为在所述轨道平面坐标系y轴的坐标值，zVc为在所述轨道平面坐标系z轴的坐标值，θc为所述参考卫星的真近点角，rc为所述参考卫星至地心的距离；
θc根据公式求得，rc根据公式rc＝ac(1-ec×cosEc)求得，Ec为所述参考卫星的偏近点角，Ec根据公式Ec＝Mc+ec×sinEc求得，Mc为所述参考卫星的平近点角，ec为所述参考卫星的轨道偏心率。
16.根据权利要求15所述的装置，其特征在于，所述第二运算模块，还用于：
所述参考卫星在所述参考卫星的第二卫星平台坐标系的位置根据公式求得，其中，AVcr2为从到的转置矩阵，
其中，为绕z轴逆时针旋转角度，γc为计算过程量。
17.根据权利要求16所述的装置，其特征在于，所述第二运算模块，还用于：
所述辅卫星在所述第二卫星平台坐标系的位置表示为：
其中，xr2d为在地球半径外方向x轴的坐标值，yr2d为在卫星运动的切线方向y轴的坐标值，zr2d为在垂直于轨道平面方向z轴的坐标值，三个坐标轴构成笛卡尔坐标系；Acd是所述辅卫星沿x轴的振幅，Bcd是所述辅卫星沿z轴的振幅，Δyr2d为所述辅卫星在所述第二卫星平台坐标系的位移，αd为所述辅卫星在轨道平面的初始位置所对应的相位，βd为所述辅卫星在垂直于轨道平面的初始位置所对应的相位。
18.根据权利要求17所述的装置，其特征在于，所述辅卫星的数量为多个。
19.一种计算机存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至9任一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910600674.8A CN110471432B (zh) | 2019-07-04 | 2019-07-04 | 一种卫星编队构型的方法、装置及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910600674.8A CN110471432B (zh) | 2019-07-04 | 2019-07-04 | 一种卫星编队构型的方法、装置及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110471432A true CN110471432A (zh) | 2019-11-19 |
| CN110471432B CN110471432B (zh) | 2020-09-08 |
Family
ID=68506825
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910600674.8A Active CN110471432B (zh) | 2019-07-04 | 2019-07-04 | 一种卫星编队构型的方法、装置及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110471432B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111409864A (zh) * | 2020-03-20 | 2020-07-14 | 上海航天控制技术研究所 | 一种编队卫星系统星间气动外形差异在轨标定方法 |
| CN111505635A (zh) * | 2020-04-16 | 2020-08-07 | 北京理工大学重庆创新中心 | 用于相干层析的geo sar双星编队构型设计方法 |
| CN113184220A (zh) * | 2021-04-21 | 2021-07-30 | 中国人民解放军63923部队 | 一种地球同步轨道通信卫星的轨道控制方法及装置 |
| CN116559917A (zh) * | 2023-05-10 | 2023-08-08 | 四川大学 | 一种定位海面动目标的被动电侦卫星编队构型设计方法 |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030029968A1 (en) * | 1997-05-21 | 2003-02-13 | Toshihide Maeda | Artificial satellite with an orbit having a long staying time in a zenith direction, an orbit control method and a communication system therewith |
| CN101381004A (zh) * | 2008-08-20 | 2009-03-11 | 南京航空航天大学 | 基于大气阻力的微小卫星编队飞行控制方法及控制装置 |
| CN101414003A (zh) * | 2008-11-28 | 2009-04-22 | 北京航空航天大学 | 一种基于星地坐标转换的星载sar图像地理编码方法 |
| CN101520511A (zh) * | 2009-03-13 | 2009-09-02 | 北京航空航天大学 | 一种分布式卫星合成孔径雷达编队构形方法 |
| CN104142686A (zh) * | 2014-07-16 | 2014-11-12 | 北京控制工程研究所 | 一种卫星自主编队飞行控制方法 |
| WO2015057102A1 (ru) * | 2013-10-18 | 2015-04-23 | Юрий Николаевич РАЗУМНЫЙ | Космическая обслуживающая система и способ ее построения |
| CN107450578A (zh) * | 2017-07-26 | 2017-12-08 | 中国人民解放军空军工程大学 | 一种绕地球同步轨道卫星分布式共轨飞行的卫星群轨道设计方法 |
| CN107633142A (zh) * | 2017-09-22 | 2018-01-26 | 上海卫星工程研究所 | 相对运动轨道构型的模拟方法 |
| CN108490973A (zh) * | 2018-04-19 | 2018-09-04 | 哈尔滨工业大学 | 航天器编队相对轨道确定方法及装置 |
| CN109240322A (zh) * | 2018-09-30 | 2019-01-18 | 南京航空航天大学 | 一种面向对地超幅宽成像的卫星编队实现方法 |
| CN109459041A (zh) * | 2018-11-27 | 2019-03-12 | 上海航天控制技术研究所 | 一种微纳星群变迁规划与控制方法 |
-
2019
- 2019-07-04 CN CN201910600674.8A patent/CN110471432B/zh active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030029968A1 (en) * | 1997-05-21 | 2003-02-13 | Toshihide Maeda | Artificial satellite with an orbit having a long staying time in a zenith direction, an orbit control method and a communication system therewith |
| CN101381004A (zh) * | 2008-08-20 | 2009-03-11 | 南京航空航天大学 | 基于大气阻力的微小卫星编队飞行控制方法及控制装置 |
| CN101414003A (zh) * | 2008-11-28 | 2009-04-22 | 北京航空航天大学 | 一种基于星地坐标转换的星载sar图像地理编码方法 |
| CN101520511A (zh) * | 2009-03-13 | 2009-09-02 | 北京航空航天大学 | 一种分布式卫星合成孔径雷达编队构形方法 |
| WO2015057102A1 (ru) * | 2013-10-18 | 2015-04-23 | Юрий Николаевич РАЗУМНЫЙ | Космическая обслуживающая система и способ ее построения |
| CN104142686A (zh) * | 2014-07-16 | 2014-11-12 | 北京控制工程研究所 | 一种卫星自主编队飞行控制方法 |
| CN107450578A (zh) * | 2017-07-26 | 2017-12-08 | 中国人民解放军空军工程大学 | 一种绕地球同步轨道卫星分布式共轨飞行的卫星群轨道设计方法 |
| CN107633142A (zh) * | 2017-09-22 | 2018-01-26 | 上海卫星工程研究所 | 相对运动轨道构型的模拟方法 |
| CN108490973A (zh) * | 2018-04-19 | 2018-09-04 | 哈尔滨工业大学 | 航天器编队相对轨道确定方法及装置 |
| CN109240322A (zh) * | 2018-09-30 | 2019-01-18 | 南京航空航天大学 | 一种面向对地超幅宽成像的卫星编队实现方法 |
| CN109459041A (zh) * | 2018-11-27 | 2019-03-12 | 上海航天控制技术研究所 | 一种微纳星群变迁规划与控制方法 |
Non-Patent Citations (4)
| Title |
|---|
| DAAN C. MAESSEN: "Relative Orbital Element Estimation and Observability Analysis for Formation Flying Satellites using Inter-Satellite Range Measurements Only", 《AIAA GUIDANCE, NAVIGATION, AND CONTROL CONFERENCE》 * |
| JUNFENG LI: "Study on relative orbital configuration in satellite formation flying", 《ACTA MECH SINICA》 * |
| 邓云凯等: "分布式星载SAR系统时间同步和波束同步误差分析", 《雷达学报》 * |
| 陆宏伟: "一种小卫星绕飞编队的运动学设计新方法研究", 《空间科学学报》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111409864A (zh) * | 2020-03-20 | 2020-07-14 | 上海航天控制技术研究所 | 一种编队卫星系统星间气动外形差异在轨标定方法 |
| CN111409864B (zh) * | 2020-03-20 | 2021-10-01 | 上海航天控制技术研究所 | 一种编队卫星系统星间气动外形差异在轨标定方法 |
| CN111505635A (zh) * | 2020-04-16 | 2020-08-07 | 北京理工大学重庆创新中心 | 用于相干层析的geo sar双星编队构型设计方法 |
| CN111505635B (zh) * | 2020-04-16 | 2022-05-03 | 北京理工大学重庆创新中心 | 用于相干层析的geo sar双星编队构型设计方法 |
| CN113184220A (zh) * | 2021-04-21 | 2021-07-30 | 中国人民解放军63923部队 | 一种地球同步轨道通信卫星的轨道控制方法及装置 |
| CN113184220B (zh) * | 2021-04-21 | 2021-11-19 | 中国人民解放军63923部队 | 一种地球同步轨道通信卫星的轨道控制方法及装置 |
| CN116559917A (zh) * | 2023-05-10 | 2023-08-08 | 四川大学 | 一种定位海面动目标的被动电侦卫星编队构型设计方法 |
| CN116559917B (zh) * | 2023-05-10 | 2023-12-12 | 四川大学 | 一种定位海面动目标的被动电侦卫星编队构型设计方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110471432B (zh) | 2020-09-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110471432A (zh) | 一种卫星编队构型的方法、装置及存储介质 | |
| Canuto et al. | Spacecraft dynamics and control: the embedded model control approach | |
| Newhall et al. | DE 102-A numerically integrated ephemeris of the moon and planets spanning forty-four centuries | |
| Konopliv et al. | A global solution for the Mars static and seasonal gravity, Mars orientation, Phobos and Deimos masses, and Mars ephemeris | |
| CN107607947B (zh) | 基于卡尔曼滤波的星载雷达成像参数在线估计方法 | |
| CN109738919B (zh) | 一种用于gps接收机自主预测星历的方法 | |
| Xu | Orbits | |
| Tillerson et al. | Advanced guidance algorithms for spacecraft formation-keeping | |
| CN112629543A (zh) | 一种大椭圆轨道及小倾角圆轨道的轨道规划方法 | |
| De Simone et al. | Attitude controller design for micro-satellites | |
| Mead | Magnetic coordinates for the Pioneer 10 Jupiter encounter | |
| Davies et al. | The unified lunar control network: 1994 version | |
| D’Errico et al. | Relative trajectory design | |
| Djachkova et al. | Characterization of the Luna-25 landing sites | |
| GB2350596A (en) | Three axis position control for low-orbiting satellites | |
| Sun et al. | The active magnetic control algorithm for HITSAT‐1 | |
| Heppenheimer | Achromatic trajectories and lunar material transport for space colonization | |
| Nazirov et al. | Mission Design Problems for the Spectrum-Roentgen-Gamma Project | |
| Vassar et al. | Orbit selection for the Stanford relativity gyroscope experiment | |
| Abdel-Aziz et al. | Equilibria of a charged artificial satellite subject to gravitational and Lorentz torques | |
| Antreasian | Precision radiation force modeling for the TOPEX/POSEIDON mission | |
| Banshchikova et al. | Technique for Conducting a Space Experiment with Auroral Imagers in Satellite Orbits | |
| Akhmetshin | Multiorbit transfers to a geostationary orbit with switching low thrust off in the shadow region | |
| Jochim | Satellite equivalence orbits | |
| Abezyaev | On One Possibility of Orientation of a Near-Earth Spacecraft in the Earth’s Magnetic Field |
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 |