CN101381004A - 基于大气阻力的微小卫星编队飞行控制方法及控制装置 - Google Patents
基于大气阻力的微小卫星编队飞行控制方法及控制装置 Download PDFInfo
- Publication number
- CN101381004A CN101381004A CNA2008100219812A CN200810021981A CN101381004A CN 101381004 A CN101381004 A CN 101381004A CN A2008100219812 A CNA2008100219812 A CN A2008100219812A CN 200810021981 A CN200810021981 A CN 200810021981A CN 101381004 A CN101381004 A CN 101381004A
- Authority
- CN
- China
- Prior art keywords
- satellite
- attitude
- formation
- control
- information
- 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
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 167
- 238000000034 method Methods 0.000 title claims abstract description 70
- 108091092878 Microsatellite Proteins 0.000 claims abstract description 32
- 238000001914 filtration Methods 0.000 claims description 41
- 239000013598 vector Substances 0.000 claims description 30
- 238000013461 design Methods 0.000 claims description 29
- 238000005259 measurement Methods 0.000 claims description 27
- 238000013016 damping Methods 0.000 claims description 17
- 230000008859 change Effects 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000012423 maintenance Methods 0.000 claims description 4
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 2
- 238000012546 transfer Methods 0.000 claims description 2
- 230000004927 fusion Effects 0.000 abstract description 4
- 238000005755 formation reaction Methods 0.000 description 135
- 230000001276 controlling effect Effects 0.000 description 16
- 239000011159 matrix material Substances 0.000 description 15
- 238000004088 simulation Methods 0.000 description 12
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 9
- 239000000446 fuel Substances 0.000 description 6
- 238000009499 grossing Methods 0.000 description 6
- 230000006641 stabilisation Effects 0.000 description 6
- 238000011105 stabilization Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000011109 contamination Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000008602 contraction Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 230000002035 prolonged effect Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 235000007926 Craterellus fallax Nutrition 0.000 description 1
- 240000007175 Datura inoxia Species 0.000 description 1
- 241000452638 Parasaissetia nigra Species 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Images
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
基于大气阻力的微小卫星编队飞行控制方法及控制装置,涉及微型航天器与航天测控技术领域。本发明的GPS接收机的信号输出端与轨道确定装置的信号输入端连接,太阳电池阵列和三轴磁强计的信号输出端与姿态融合滤波的输入端连接,三轴磁强计的输出端与纯磁测定姿的输入端连接,姿态融合滤波的信号输出端与星间链路的第二输入端连接,星间链路的第一输出端与编队队形控制律的输入端连接,星间链路的第二输出端与卫星姿态控制律的输入端连接,通过控制气动板来进行编队卫星飞行的队形和姿态控制。本发明实现了延长了航天器的寿命,功耗低、成本低、无污染的基于大气阻力的目的。
Description
技术领域
本发明涉及微型航天器与航天测控技术领域。
背景技术
微小卫星技术近十年来一直是国内外航天领域研究的热点,而与此相关的若干关键技术的研究已经吸引了包括高校在内的众多研究机构的关注,由于微小卫星具有重量轻、体积小、成本低、研制周期短以及功能密度高等一系列优点,在通讯、遥感、军事、行星探测、工程技术实验等领域发挥着重要作用。研究开发具有高指向精度和稳定度的姿态测控技术具有十分重要的意义。
编队飞行是微小卫星的一个重要发展方向,其兴起标志着航天技术进入一个崭新的时代,微小卫星编队飞行是由多颗彼此之间距离较近的卫星组成,具有特定的空间几何构形,通过星间链路实现协同工作,不仅能完成特定的航天任务,而且在体积、性能、费用等方面比单颗大卫星具有明显的优势。
为了较好的实现卫星编队功能,编队飞行的各卫星必须密切配合协同工作,其姿态和轨道控制是实现编队飞行的关键技术之一。传统的卫星姿态控制通常采用动量轮或推进系统提供控制力矩,编队飞行卫星的队形控制通常使用推进系统来提供推力。现有的技术主要存在以下三个问题:一是采用动量轮或者推进系统的控制系统不仅技术复杂而且成本;二是微小卫星携带燃料有限,姿态和队形控制需要燃料的持续消耗,一旦耗尽卫星就失去了价值;三是推进器的羽流对卫星的冲击不仅使其稳定性受到了干扰,而且对精密的星载仪器造成污染。本发明给出了一种利用大气阻力进行微小卫星编队飞行控制的装置,给出了相应的系统配置和姿态轨道控制方法,使用本发明将使得卫星的编队飞行任务将不受限于卫星所载燃料的多少,延长了航天器的寿命,尤其适用于低功耗、低成本、轻重量、无污染的卫星编队飞行任务,具有很好的应用前景。
发明内容
本发明目的是提供一种利用大气阻力进行微小卫星编队飞行轨道和姿态控制的方法,延长了航天器的寿命,功耗低、成本低、无污染的基于大气阻力的微小卫星编队飞行控制方法及控制装置。
本发明为实现上述目的,采用如下技术方案:
本发明包括GPS接收机、轨道确定装置、星间链路、太阳电池阵列、姿态融合滤波器、三轴磁强计、纯磁测定姿、编队队形控制律、卫星姿态控制律、气动板;GPS接收机的信号输出端与轨道确定装置的信号输入端连接,将GPS接收机接收的轨道信号经过轨道确定装置计算后得到准确的轨道状态;太阳电池阵列和三轴磁强计的信号输出端与姿态融合滤波的输入端连接,姿态融合滤波器对太阳电池阵列测量的太阳方位信息和磁强计测量的地球磁场信息进行滤波,得到高精度姿态信息;三轴磁强计的输出端与纯磁测定姿的输入端连接,将三轴磁强计采集的磁场信息经过纯磁测定姿的计算后得到卫星姿态信息;姿态融合滤波的信号输出端与星间链路的第二输入端连接,纯磁测定姿态的信号输出端与星间链路的第三输入端连接,经过计算得到各卫星相对和绝对的轨道与姿态信息;星间链路的第一输出端与编队队形控制律的输入端连接,将星间链路中的准确的轨道信息传输到编队队形律中;星间链路的第二输出端与卫星姿态控制律的输入端连接,将星间链路中的准确的姿态信息传输到卫星姿态控制律中;编队队形控制律、卫星姿态控制律的输出端分别与气动板连接,通过控制气动板来进行编队卫星飞行的队形和姿态控制。
比较好的是,本发明的小卫星为立方体结构,气动板设置在小卫星飞行方向的背面,气动板与小卫星的飞行方向垂直,通过控制气动板的伸缩来改变卫星的迎风面,一方面产生气动力矩来控制卫星姿态,另一方面使得编队飞行的各卫星之间产生差动大气阻力,用来对编队卫星的队形进行控制。
本发明基于大气阻力的微小卫星编队飞行的控制装置的控制方法:
第一步:卫星轨道确定。通过GPS接收机的信号确定卫星的轨道信息;若GPS信号丢失,则采用广义卡尔曼滤波方法进行轨道信息的递推;
第二步:卫星姿态确定;
第三步:编队飞行卫星相对轨道和姿态确定:经过第一步和第二步得到各颗卫星的轨道和姿态信息后,通过星间链路得到各卫星间的相对和绝对的轨道和姿态信息;
第四步:编队卫星姿态控制:上述第三步的相对和绝对姿态信息传输到卫星姿态控制律中,卫星姿态控制律控制气动板与大气的有效接触面积,应用气动力矩进行滚动、俯仰和偏航方向的姿态;
第五步:编队飞行队形控制:上述第三步的相对和绝对轨道信息传输到编队队形控制律中,并判断第四步的的姿态控制是否完成,若姿态控制完成后通过卫星编队队形控制律驱动气动板与大气的有效接触面积,应用各卫星间的大气差动阻力进行编队卫星的队形控制;
第六步:星间链路对各卫星之间的姿态信息与设定的姿态参数进行比对判断,若判断完成姿态控制任务,则继续监测;若判断为否,则执行上述第四步;
第七步:星间链路对各卫星之间的轨道信息与设定的编队飞行的队形进行判断,若判断完成队形控制任务,则继续监测;若判断为否,则执行上述第五步。
比较好的是,本发明的卫星姿态确定存在的以下三种情况:
情况1:当仅仅需要角速率信息时,利用三轴磁强计测量的磁场信息,通过纯磁测定,采用磁场差分的B-dot方法来间接的确定卫星的角速率变化;
情况2:当同时需要卫星的三轴姿态和角速率信息,并且卫星处于太阳光照区时,利用太阳电池阵列、三轴磁强计测量的太阳方位和地球磁场信息,以及经过第一步得到的卫星轨道信息,通过姿态融合滤波器来确定高精度的卫星姿态;
情况3:当同时需要卫星的三轴姿态和角速率信息,并且卫星处于地球阴影区时,利用三轴磁强计测量的磁场信息,以及经过第一步得到的卫星轨道信息,通过纯磁测定姿,采用EKF的方法来确定卫星的姿态信息。
比较好的是,本发明的第二步的卫星姿态确定方法,采用太阳电池阵列和磁强计组合,通过卫星动力学方程计算角速度来设计无陀螺滤波器,滤波器采用双矢量定姿q-method方法和非线性EKF滤波相融合,首先应用q-method方法将太阳矢量和磁场矢量的六维观测量转换为四元数,降低观测维数并将观测方程转换为线性,然后代入EKF进行滤波计算,得到卫星的姿态信息。
比较好的是,本发明第四步的编队飞行卫星的姿态控制,仅仅采用气动板作为执行机构,根据纯磁测差分计算的磁场变化率信息设计卫星速率阻尼控制律,根据姿态确定算法得到的姿态四元数信息和角速度信息设计卫星姿态捕获及稳定控制律,得到期望的气动力矩后进一步设计气动板驱动控制律,从而控制编队飞行卫星的姿态。
比较好的是,本发明第五步编队飞行卫星的队形控制仅仅采用大气阻力,根据编队飞行的各卫星的相对位置信息,通过建立编队卫星的相对运动模型,针对共面跟飞和共面椭圆编队的终端约束,设计最优控制律进行卫星编队飞行的队形建立、队形重构以及队形保持等功能控制。
本发明的基于大气阻力的微小卫星编队飞行控制系统,在卫星入轨后的主要工作模式及工作流程为:
1、消除星箭分离引起的初始扰动
对于编队飞行的各卫星,利用卫星本体所处的磁场变化情况,采用磁强计测量的磁场变化信息,间接的测量卫星角速率变化,采用设计的磁场变化率反馈控制律来控制气动板,应用气动力矩进行滚动、俯仰和偏航方向的速率阻尼。
2、姿态捕获及三轴姿态稳定控制
对于编队飞行的各卫星,在速率阻尼完成后,光照区时启动太阳电池阵列加磁强计的姿态确定组合滤波器,阴影区时采用纯磁测定姿方法,得到各颗卫星的绝对姿态后可以根据星间链路计算卫星间的相对姿态,根据任务的需要采用相应的姿态反馈控制律驱动气动板,应用气动力矩进行卫星各轴的绝对及相对的姿态控制。
3、卫星编队飞行的队形建立控制
对于编队飞行的各卫星建立姿态基准后,根据编队队形的任务要求,应用大气阻力建立编队队形。由于大气阻力仅可以提供沿航迹向的作用力,而不能提供轨道面法向的作用力,因此本发明的基于大气阻力编队队形控制主要针对共面编队飞行的情况,如共面的跟飞、绕飞、椭圆等编队控制。根据初始条件和最终目标设计最优控制律,根据卫星间的差动大气阻力形成期望的编队队形。
4、卫星编队飞行的队形保持控制
编队飞行的卫星,在轨运行期间会受到各种摄动力,使之不能保持预设的运行轨迹,从而破坏整个编队队形。根据各卫星间的相对轨道和姿态信息,设计最优控制律,通过控制气动板来保持编队卫星间的相对位置,从而保证编队队形的稳定性。
5、卫星编队飞行的队形重构控制
对于编队飞行的各卫星,根据任务的要求需要进行队形重构,例如对于共面编队飞行,根据需要控制卫星跟飞距离的增大或减小,控制绕飞椭圆的膨胀与收缩等,根据任务要求设计最优控制律,通过控制气动板来进行共面编队队形的重构。
本发明采用上述技术方案,与现有技术相比具有如下优点:
1.国内外卫星的轨道控制通常都采用推进系统,不仅系统复杂,而且需要燃料消耗。微小卫星由于体积、质量、功耗等限制,携带燃料有限,一旦耗尽卫星就失去了价值。同时推进器的羽流不仅会对星载仪器造成污染,而且对卫星造成冲击使其稳定性受到干扰。采用本发明可以使编队卫星飞行时间不受限于所载燃料的多少,尤其适用于低功耗、低成本、轻重量、无污染的卫星编队飞行任务。
2.对于编队飞行的相对和绝对姿态控制,以及单颗卫星的姿态控制,通常的姿态执行机构采用推进系统、动量轮(包括反作用轮和偏置动量轮)以及磁力矩器等,或者采用不同的执行机构组合,姿态控制系统设计较为复杂。本发明的姿态控制仅仅使用气动板,不仅是一种新的姿态控制方法,而且姿态控制的同时还可以进行编队控制,一物多用提高功能密度,显著减小控制系统的复杂性,符合低成本、低功耗的微小卫星设计思想。
3.本发明设计了一种融合姿态确定性算法q-method和非线性最优滤波EKF的双重姿态确定滤波器,根据太阳电池阵列和磁强计测量,结合磁场模型和太阳矢量模型进行姿态确定,不仅大大减小采用直接观测方法的计算量,而且可以提供姿态确定精度。本发明设计了基于大气阻力的编队卫星姿态控制和队形控制律,有效的进行微小卫星编队飞行控制。本发明中的姿态确定方法与编队控制方法,不仅应用于基于大气阻力的编队飞行控制,而且可以为采用其它控制方式的编队卫星控制也提供了重要参考。
4.本发明的大气阻力进行微小卫星编队飞行的控制系统,具有成本低、重量轻、功耗小、无污染、寿命长的优点,完全符合微小卫星的“快、好、省”的发展原则。本发明中采用的各项关键技术及创新思想,可以推广到空间飞行器的轨道和姿态控制系统的设计中,在中、低轨道卫星编队飞行、星座系统控制以及单颗卫星姿态轨道控制等领域的研究中,具有非常广阔的工程应用前景。
附图说明
图1是本发明的基本原理图。
图2a是本发明的卫星结构示意图。
图2b是图2的俯视图。
图2c是图2的左视图。
图2d是图2的侧视图。
图3是本发明的飞行姿态确定滤波器流程。
图4a是本发明的共面跟飞编队状态示意图。
图4b是本发明的共面椭圆编队状态示意图。
图5a是本发明的卫星飞行姿态确定仿真结果的参考运动轨迹示意图。
图5b是本发明的卫星飞行姿态确定仿真结果的直接观测EKF时四元数误差示意图。
图5c是本发明的卫星飞行姿态确定仿真结果的双矢量定姿时四元数误差示意图。
图5d是本发明的卫星飞行姿态确定仿真结果的双重滤波器时四元数误差示意图。
图5e是本发明的卫星飞行姿态控制仿真结果的速率阻尼控制角速度结果示意图。
图5f是本发明的卫星飞行姿态控制仿真结果的速率阻尼控制时气动板长度变化示意图。
图5g是本发明的卫星飞行姿态控制仿真结果的三轴姿态控制四元数结果示意图。
图5h是本发明的卫星飞行姿态控制仿真结果的三轴姿态控制时气动板长度变化示意图。
图5i是本发明的卫星飞行队形控制仿真结果的共面跟飞编队队形建立示意图。
图5j是本发明的卫星飞行队形控制仿真结果的共面椭圆编队队形建立示意图。
图5k是本发明的卫星飞行队形控制仿真结果的共面跟飞→椭圆编队重构示意图。
图5l是本发明的卫星飞行队形控制仿真结果的共面椭圆→跟飞编队重构示意图。
具体实施方式
下面结合附图对本发明的技术方案进行详细说明:
1、技术实施整体路线:
如图1所示,本发明包括GPS接收机11、轨道确定装置12、星间链路13、太阳电池阵列14、姿态融合滤波器15、三轴磁强计16、纯磁测定姿17、编队队形律18、卫星姿态控制律19、气动板;GPS接收机11的信号输出端与轨道确定装置12的信号输入端连接,将GPS接收机11的信号经过轨道确定装置12计算后得到准确的轨道状态,轨道确定装置12的第一信号输出端与星间链路13的第一信号输入端连接;太阳电池阵列14的信号输出端与姿态融合滤波器15的第一信号输入端连接,三轴磁强计16的信号输出端与姿态融合滤波器15的第二信号输入端连接,姿态融合滤波器15对太阳电池阵列14测量的太阳方位信息以及三轴磁强计16测量的磁场信息进行融合滤波计算后得到高精度姿态信息;三轴磁强计16的输出端与纯磁测定姿17的输入端连接,将三轴磁强计16测量的磁场信息经纯磁测定姿17的计算后得到姿态信息;姿态融合滤波器15的信号输出端与星间链路13的第二输入端连接,纯磁测定姿17的信号输出端与星间链路13的第三输入端连接;星间链路13的第一输出端与编队队形控制律18的输入端连接,将经过星间链路13处理后得到的相对和绝对的轨道信息传输到编队队形控制律18中;星间链路13的第二输出端与卫星姿态控制律19的输入端连接,将星间链路13处理后得到的相对和绝对姿态信息传输到卫星姿态控制律19中;编队队形控制律18、卫星姿态控制律19的输出端分别与气动板连接,通过控制气动板来进行编队卫星飞行的队形和姿态控制。
2、卫星设计、大气阻尼及气动力矩分析:
上述基于大气阻力的微小卫星编队飞行控制系统所设计的卫星结构示意如附图2a、2b、2c、2d所示。与目前国际上大多数的小卫星、纳卫星以及皮卫星类似,卫星采用长方体结构,所不同的是在飞行方向的背面上安装四块轻型的可以伸缩的气动板,作为姿态控制和队形控制的执行机构。记卫星长、宽和高的尺寸分别为a、b和c,气动板宽为e,通过调节长度d来控制气动板的面积。对于姿态控制坐标系定义为:轨道系Fo原点位于卫星质心,飞行方向为x轴,y轴指向轨道面负法向,z轴按右手规则定义;星体系Fb原点位于卫星质心,三轴姿态稳定时与轨道系重合;惯性系Fi原点位于地心,x轴指向春分点z轴沿地球自转轴,y轴由右手规则定义。给卫星各个外表面进行编号:1和2分别对应x方向的两个面,3和4对应y轴,5和6对应z轴,4块气动板分别编号为7、8、9和10。则卫星六个面的中心以及四块气动板的中心相对于卫星质心的矢径分别为:
当卫星轨道高度在120km以上时,大气运动可以看作自由分子流,此时相对于轨道系的大气阻力和气动力矩公式分别为
其中,a=(1.5ωe/ωo)sin i cos ωot,t为从升交点开始计算的时间。大气密度计算采用指数衰减大气模型
ρ=ρo exp(-(h-ho)/H) (6)
其中,ρo为参考高度为ho时的参考大气密度值,h为实际高度,H为标高。
对于本发明的编队卫星,对卫星每一个面所受的大气阻尼和气动力矩进行计算,然后进行求和得到整个卫星所受大气阻尼和气动力矩分别为
3、编队卫星的轨道姿态确定方法设计:
本发明的微小卫星编队飞行姿态确定滤波器设计流程如附图3所示。采用四元数 描述卫星的姿态,则轨道系到星体系的姿态旋转矩阵为
上式的矢量形式为 以轨道系为参考坐标系,则由四元数来描述得到的卫星姿态运动学方程为
在气动力矩控制下,卫星的动力学方程为
姿态确定方法主要有两个分支,其一为最小二乘估计问题,其二为估计滤波法。最小二乘法通常称为Wahba问题:在参考系中,有一组参考矢量在星体系中对应另一组观测矢量其中n≥2,根据这两组矢量来求解两坐标系之间的姿态矩阵C,使得如下误差能量函数最小
其中σi为权系数,函数还可以写为 则Wahba最小能量函数问题可以转换为求如下增益函数的最大值
g(C)=trace(CBT) (13)
姿态确定的估计滤波方法采用EKF,其递推方程组如下
其中,Kk为滤波增益,Pk/k-1和Pk分别为一步预测和估计均方差。式(14~15)为时间更新方程,式(16~18)为量测更新方程。
本发明的基于大气阻力的微小卫星编队飞行控制系统中,队形控制和姿态确定都应用到轨道信息,采用GPS进行轨道确定,当没有GPS信号时采用递推滤波进行轨道估计。当处于阴影区时姿态确定采用纯磁测加卫星动力学方程进行无陀螺估计;当处于光照区时应用太阳电池阵列和磁强计测量进行高精度融合估计。轨道确定和纯磁测姿态确定都采用EKF方法,滤波步骤如式(14~18),下面重点对本发明所设计的双重姿态融合滤波器进行说明。
本发明的双重滤波器设计首先采用q-method法进行双矢量定姿,它是采用四元数来求解Wahba问题的有效方法。将姿态矩阵(8)代入式(13)得
其中,K为4×4维矩阵,定义为
其中,Z=[B23-B32 B31-B13 B12-B21]T。则K的最大特征值λmax对应的特征向量即为所求的最优四元数qopt,即Kqopt=λmaxqopt。
当应用q-method算法得到四元数后,进一步设计EKF滤波器。由于四元数存在约束条件‖q‖=1,如果直接用四元数作为状态量,EKF协方差矩阵会出现奇异。本发明的滤波器设计中采用伪测量方法,在滤波的状态变量中除去四元数标量部分,且在滤波器状态更新时将四元数与其它状态量分别计算,四元数采用乘法而角速度采用加法运算,然后再统一进行滤波循环,该方法不仅可以保持四元数的正规化,而且降低状态维数从而降低计算量。观测量为由q-method计算得到四元数qk的矢量部分,此时观测方程化为线性方程,滤波增益Kk为6×3维矩阵,由四元数和角速度两部分组成即Kk=[Kq Kω]T。状态估值计算为
对于无陀螺EKF滤波器设计,上述的四元数正规化处理方法,状态量选为四元数矢量部分和角速度的6×1维矢量 则由卫星运动学与动力学模型,可得非线性状态方程为
此时观测方程为非线性,在滤波每一步中都需要由滤波估计值得到的进行坐标转换,输出方程非常复杂。采用本发明的q-method算法进行预处理,作为EKF滤波循环外的一部分,这样可以将6维观测量转换为3维观测量,滤波循环中不需要进行坐标转换,同时观测方程也化成线性方程,大大减少了计算量。记q-method法计算得到四元数为qc,选取qc的矢量部分为观测量,则观测量为Z=[qc,x qc,y qc,z]T,观测方程转化为线性方程
其中,Vc为量测噪声,即使用q-method的计算误差。根据以上的状态方程和预处理观测方程,带入EKF递推方程组(14~18)即可进行姿态估计。
4、基于大气阻力的卫星编队飞行的姿态控制律设计:
卫星的姿态控制通常需要姿态和角速率信息,本发明可以在没有角速率测量单机(如陀螺)的前提下实现卫星的状态控制,有利于微小卫星的系统简化和节省资源。本发明的编队卫星姿态控制主要分为消除星箭分离干扰的速率阻尼和进行姿态捕获建立三轴稳定的两个阶段。
其中,Kdet为正定控制增益,为简便期间,选为对角阵Kdet=diag{[kx ky kz]}。为磁场变化率,由磁强计的测量值进行差分计算得到的估计值代替。该方法借鉴于目前国际上使用磁控进行速率阻尼的B-dot控制方法,所不同的是控制机构为气动板而不是磁力矩器。可以近似等于卫星的角速率变化,由于磁强计所测磁场变化是由卫星飞行引起和卫星状态引起的地磁变化之和,由于飞行位置变化引起一个控制周期的磁场变化极小,因此主要反映的姿态角速度变化,根据的信息进行控制即可以阻尼角速率。
对于姿态捕获建立三轴稳定阶段,设计限制姿态反馈控制律
其中,Lpmax为姿态反馈限制系数,可以通过调整Lpmax来保证该控制律满足大角度姿态捕获的要求,Kp和Kd为正定反馈增益。该控制律是在PD控制的基础上,采用限制姿态反馈的方法,使得在大角度姿态偏差时,角速度控制起主要作用,当姿态控制到一定范围后进入角速度和姿态共同作用,即为小角度下的PD控制,适合整个姿态捕获和三轴稳定阶段。
当设计好速率阻尼和姿态捕获及三轴稳定的控制律之后,需要根据期望的力矩来控制4块气动板伸缩的长度di(i=1,2,3,4),从而产生需要的气动力矩气动力矩是卫星姿态q和气动板长度di的非线性函数,当卫星姿态和气动板长度已知时,由式(7)可以计算控制力矩,然而由期望的控制力矩求解气动板长度的逆运算将比较困难,本发明采用线性化处理的逆运算了来设计气动板驱动控制律。气动板的控制,实际产生的气动力矩与期望的控制力矩误差为 在一个控制周期内期望力矩为常数,因此可得到 其中为4块气动板的长度变化率,A为雅可比矩阵
其中,Kf和Kp为正定增益矩阵,均取为对角阵,为A的逆矩阵。
5、基于大气阻力的卫星编队飞行的队形控制律设计:
参见附图4本发明的基于大气阻力微小卫星编队飞行控制策略,当姿态控制完成后,通过姿态信息判断指向正确后进行编队的队形控制。编队队形控制是应用多颗卫星间的差分气动阻力进行控制:为了保持相对姿态,按同样的比例同时伸长某一卫星的4块气动板,来改变其迎风面,从而改变卫星间的气动阻力,导致卫星相对于参考坐标产生不同的加速度,进而改变卫星的相对位置。
采用Hill方程(亦称C-W方程)描述编队卫星的相对运动,根据国内文献资料的习惯重新定义轨道坐标系:x轴指向飞行方向,y轴背向地心方向,z轴垂直轨道面与x和y构成右手坐标,则Hill方程为
其中,ωo为轨道角速度, 为相对于卫星各轴的作用力,m为卫星的质量。当没有作用力即编队卫星自由飞行时,对Hill方程进行求解分析,可以将卫星运动分解为轨道面(xy平面)和垂直于轨道面(z方向)的两个独立运动,xy平面内相对运动为椭圆,z方向为简谐振动。
大气阻尼会造成卫星编队构形产生沿航迹方向的漂移,而不产生垂直于轨道面的作用力,因此本发明的编队队形控制仅考虑共面的编队飞行。卫星的编队构形主要有两类:其一为主从结构编队(leader-follower),其二为虚拟参考卫星编队。主从结构编队构型的主卫星位于参考坐标系,从卫星相对于主卫星运动;虚拟参考卫星编队的参考中心没有实质的卫星,各卫星相对于某个参考中心运动,是一种更为稳定的编队形式,不会因为主卫星失效而使整个编队系统实效。因此,下面主要针对共面虚拟参考卫星编队的两种编队队形进行控制律设计:共面跟飞和共面椭圆编队,参见附图4a、附图4b的基于大气阻力的微小卫星编队飞行队形示意图。
对于共面编队卫星,以双星编队为例设计控制律,选取控制状态量为
大气阻力在卫星各轴向的分量为 fy=0, 由于a=(1.5ωe/ωo)sin i cos ωot<1×10-4,则fz小于fx的10-4量级,故可不做考虑。记cdAi/mi=γi(i=1,2),选取控制量为u=[γ1 γ2]T,则系统状态方程为
本发明的卫星编队队形控制采用最优控制方法的LQ终端控制器进行控制律设计。采用大气阻力控制卫星时,主要通过伸缩气动板长度来改变迎风面面积,从而改变控制量u=[γ1 γ2]T。但是大气阻力增加将降低卫星速度,进而降低卫星轨道半径,因此编队控制的同时应尽量减少轨道降低。定义轨道能量函数为 其中r为轨道速度和半径,a为轨道半长轴,则减少轨道半径将降低轨道能量。对于卫星编队系统(33),设系统初始状态为X(t0)=X0,采用LQ软终端约束控制的目标是极小化指标泛函
(34)
其中,u=[cdA1/m1 cdA2/m2]T为控制量, 为正定实对称矩阵;Qf为半正定实对称矩阵,ef为终端误差向量,Mf和ψ分别为终端状态矩阵和向量。由系统状态方程和指标泛函,得哈密顿函数
H=1/2uT Ru+λT(AX+Bu) (35)
其中,λ为Lagrange乘子,则由最优性的必要条件得到最优控制律u(t)为
u(t)=-R-1BTλ(t) (36)
由控制律(36)及哈密顿函数(35)可以得到哈密顿正则微分方程,以及软终端约束LQ控制导出的哈密顿两端边值问题的边界条件分别为。
X(t0)=X0,
不同的卫星编队队形具有不同的终端约束,通过对本发明附图5的两种编队卫星任务及队形特征的分析,得到共面跟飞的终端状态矩阵和向量为
ψ=[0 l/2 0 0 0 l/2 0 0]T (39)
共面椭圆编队的终端矩阵Mf与共面跟飞编队的相同,所不同的是终端向量,为
ψ=[0 l ωol/2 0 0 l ωol/2 0]T (40)
对于本发明的基于大气阻力的共面编队飞行队形控制,根据编队任务的要求,如队形建立、队形重构以及队形保持等任务,设置初始条件和终端目标后,便可以根据上述的LQ终端控制律计算得到控制量u=[γ1 γ2]T,即可以得到控制量为[cdA1/m1 cdA2/m2]T。卫星的阻力系数cd和质量mi(i=1,2)是固定的,则由控制量u计算需要的迎风面A,而迎风面大小可以由气动板的伸缩来控制。
6、有益效果分析:
根据以上五个主要技术实施方式,对本发明的有益效果进行分析。首先针对步骤2所设计的编队卫星建立卫星相对运动模型,在此基础上对编队飞行姿态确定与控制、编队队形控制进行分析。根据目前国内外的微小卫星设计指标分析,设编队飞行微小卫星质量均为m=20kg,转动惯量I=[1.0 1.05 0.9]Tkg·m2,卫星尺寸为a×b×c=0.6×0.6×0.7m3,4块气动板的宽均为e=0.3m,长度的控制范围为0m≤d≤0.35m,编队飞行高度为600km。
对于编队卫星飞行姿态确定,根据磁强计及太阳电池阵列测量性能指标的实际特性分析,取磁强计测量误差为300nT(3σ),太阳电池阵列测量误差2°(3σ),滤波周期为1s。姿态捕获阶段的姿态运动变化较为丰富,因此选取该阶段进行测量,运动轨迹如附图5a,三轴初始姿态为[40 60 150]°,采用6维直接观测量的结果、单独q-method方法测量结果和本发明的双重滤波结果如图5b、5c、5d所示,可以看出采用直接EKF滤波器时,在0~2000s大姿态变化时误差达到20°,而使用本发明的双重滤波后误差在5°以内,在2000s以后的稳态阶段,不同滤波器的比较结果见表1所示,可以看出双重滤波器具有最高的精度,而且本发明的双重滤波器将6维观测量转化为四元数,将观测方程化为线性,计算量要比直接EKF滤波器大大减小。
表1 滤波稳定后不同无陀螺滤波器结果比较
对于基于大气阻力的编队卫星姿态控制,采用前面的卫星设计指标,气动板控制变化范围为0m≤d≤0.35m,给出速率阻尼、三轴姿态控制情况的有益效果。速率阻尼时卫星三轴初始角速度为[0.6 0.6 0.6]°/s,约为10倍的轨道角速度,控制结果如附图5e、5f所示。卫星三轴姿态控制时初始姿态为[35 30 -20]°,控制目标为三轴姿态均为0°,即姿态四元数为[1 0 0 0],控制结果如附图5g、5h所示。可以看出本发明的基于大气阻力的编队卫星姿态控制具有效果良好,通过多次实验结果表明在编队卫星处于1000km以下中、低轨道时均可有效的进行姿态控制,轨道越低则达到目标姿态的控制速度越快。
对于基于大气阻力的编队卫星队形控制,为了保持编队卫星的相对姿态,各卫星4块气动板应以相同的比例进行伸缩控制,若队形控制开始时各颗卫星的四块气动板长度相同,则队形控制过程中应保持d1=d2=d3=d4。对于附图4a的共面跟飞和附图4b共面椭圆两种队形控制,采用本发明所设计的队形控制律进行控制,队形建立时设两星在x和y方向的初始位置为[-10 10]m和[-10 10]m,初始速度为[-0.02 0.02]m/s和[0.02 -0.02]m/s,共面跟飞终端目标为两星位于同一轨道相距1000m,椭圆编队终端目标为半长轴500m,控制时间为5个轨道周期,结果如附图5i、5j,其中符号“”为编队相对参考中心;队形重构控制时,共面跟飞→椭圆编队重构的初始距离为1000m,终端目标为椭圆半长轴100m,共面椭圆→跟飞编队重构时初始椭圆长半轴500m,终端目标相距100m,结果如附图5k、5l所示。可以看出本发明的基于大气阻力的编队卫星队形控制具有效果良好,通过多次实验结果表明在编队卫星处于1000km以下中、低轨道时均可有效的进行队形建立、重构以及保持等控制。
Claims (7)
1、一种基于大气阻力的微小卫星编队飞行的控制装置,其特征在于:包括GPS接收机(11)、轨道确定装置(12)、星间链路(13)、太阳电池阵列(14)、姿态融合滤波器(15)、三轴磁强计(16)、纯磁测定姿(17)、编队队形控制律(18)、卫星姿态控制律(19)、气动板;GPS接收机(11)的信号输出端与轨道确定装置(12)的信号输入端连接,将GPS接收机(11)的信号经过轨道确定装置(12)计算后得到准确的轨道状态,轨道确定装置(12)的第一信号输出端与星间链路(13)的第一信号输入端连接;太阳电池阵列(14)的信号输出端与姿态融合滤波器(15)的第一信号输入端连接,三轴磁强计(16)的信号输出端与姿态融合滤波器(15)的第二信号输入端连接,姿态融合滤波器(15)对太阳电池阵列(14)测量的太阳方位信息以及三轴磁强计(16)测量的磁场信息进行融合滤波计算后得到高精度姿态信息;三轴磁强计(16)的输出端与纯磁测定姿(17)的输入端连接,将三轴磁强计(16)测量的磁场信息经纯磁测定姿(17)的计算后得到姿态信息;姿态融合滤波器(15)的信号输出端与星间链路(13)的第二输入端连接,纯磁测定姿(17)的信号输出端与星间链路(13)的第三输入端连接;星间链路(13)的第一输出端与编队队形控制律(18)的输入端连接,将经过星间链路(13)处理后得到的相对和绝对的轨道信息传输到编队队形控制律(18)中;星间链路(13)的第二输出端与卫星姿态控制律(19)的输入端连接,将星间链路(13)处理后得到的相对和绝对姿态信息传输到卫星姿态控制律(19)中;编队队形控制律(18)、卫星姿态控制律(19)的输出端分别与气动板连接,通过控制气动板来进行编队卫星飞行的队形和姿态控制。
2、根据权利要求1所述的基于大气阻力的微小卫星编队飞行的控制装置,其特征在于上述小卫星为立方体结构,气动板设置在小卫星飞行方向的背面,气动板与小卫星的飞行方向垂直,通过控制气动板的伸缩来改变卫星的迎风面,一方面产生气动力矩来控制卫星姿态,另一方面使得编队飞行的各卫星之间产生差动大气阻力,用来对编队卫星的队形进行控制。
3、一种实现权利要求1所述的基于大气阻力的微小卫星编队飞行的控制装置的控制方法,其特征在于:
第一步:卫星轨道确定:通过GPS接收机(11)的信号确定卫星的轨道信息;若GPS信号丢失,则采用广义卡尔曼滤波方法进行轨道信息的递推;
第二步:卫星姿态确定;
第三步:编队飞行卫星相对轨道和姿态确定:经过第一步和第二步得到各颗卫星的轨道和姿态信息后,通过星间链路(13)得到各卫星间的相对和绝对的轨道和姿态信息;
第四步:编队卫星姿态控制:上述第三步的相对和绝对姿态信息传输到卫星姿态控制律(19)中,卫星姿态控制律(19)控制气动板与大气的有效接触面积,应用气动力矩进行滚动、俯仰和偏航方向的姿态;
第五步:编队飞行队形控制:上述第三步的相对和绝对轨道信息传输到编队队形控制律(18)中,并判断第四步的的姿态控制是否完成,若姿态控制完成后通过卫星编队队形控制律(18)驱动气动板与大气的有效接触面积,应用各卫星间的大气差动阻力进行编队卫星的队形控制;
第六步:星间链路(13)对各卫星之间的姿态信息与设定的姿态参数进行比对判断,若判断完成姿态控制任务,则继续监测;若判断为否,则执行上述第四步;
第七步:星间链路(13)对各卫星之间的轨道信息与设定的编队飞行的队形进行判断,若判断完成队形控制任务,则继续监测;若判断为否,则执行上述第五步。
4、根据权利要求3所述的基于大气阻力的微小卫星编队飞行的控制装置的控制方法,其特征在于:上述第二步的卫星姿态确定存在的以下三种情况:
情况1:当仅仅需要角速率信息时,利用三轴磁强计(16)测量的磁场信息,通过纯磁测定(17),采用磁场差分的B-dot方法来间接的确定卫星的角速率变化;
情况2:当同时需要卫星的三轴姿态和角速率信息,并且卫星处于太阳光照区时,利用太阳电池阵列(14)、三轴磁强计(16)测量的太阳方位和地球磁场信息,以及经过第一步得到的卫星轨道信息,通过姿态融合滤波器(5)来确定高精度的卫星姿态;
情况3:当同时需要卫星的三轴姿态和角速率信息,并且卫星处于地球阴影区时,利用三轴磁强计(6)测量的磁场信息,以及经过第一步得到的卫星轨道信息,通过纯磁测定姿(17),采用广义卡尔曼滤波方法来确定卫星的姿态信息。
5、根据权利要求3所述的基于大气阻力的微小卫星编队飞行的控制装置的控制方法,其特征在于:上述第二步的卫星姿态确定方法,采用太阳电池阵列和磁强计组合,通过卫星动力学方程计算角速度来设计无陀螺滤波器,滤波器采用双矢量定姿q-method方法和非线性EKF滤波相融合,首先应用q-method方法将太阳矢量和磁场矢量的六维观测量转换为四元数,降低观测维数并将观测方程转换为线性,然后代入广义卡尔曼滤波方法进行滤波计算,得到卫星的姿态信息。
6、根据权利要求3所述的基于大气阻力的微小卫星编队飞行的控制装置的控制方法,其特征在于:上述第四步的编队飞行卫星的姿态控制,采用气动板作为执行机构,根据纯磁测差分计算的磁场变化率信息设计卫星速率阻尼控制律,根据姿态确定算法得到的姿态四元数信息和角速度信息设计卫星姿态捕获及稳定控制律,得到期望的气动力矩后进一步设计气动板驱动控制律,从而控制编队飞行卫星的姿态。
7、根据权利要求3所述的基于大气阻力的微小卫星编队飞行的控制装置的控制方法,其特征在于:上述第五步编队飞行卫星的队形控制仅仅采用大气阻力,根据编队飞行的各卫星的相对位置信息,通过建立编队卫星的相对运动模型,针对共面跟飞和共面椭圆编队的终端约束,设计最优控制律进行卫星编队飞行的队形建立、队形重构以及队形保持等功能控制。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100219812A CN101381004B (zh) | 2008-08-20 | 2008-08-20 | 基于大气阻力的微小卫星编队飞行控制方法及控制装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100219812A CN101381004B (zh) | 2008-08-20 | 2008-08-20 | 基于大气阻力的微小卫星编队飞行控制方法及控制装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101381004A true CN101381004A (zh) | 2009-03-11 |
CN101381004B CN101381004B (zh) | 2010-11-10 |
Family
ID=40461105
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008100219812A Expired - Fee Related CN101381004B (zh) | 2008-08-20 | 2008-08-20 | 基于大气阻力的微小卫星编队飞行控制方法及控制装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101381004B (zh) |
Cited By (37)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101793526A (zh) * | 2010-04-12 | 2010-08-04 | 哈尔滨工业大学 | 一种多信息融合编队航天器自主相对导航方法 |
CN101833285A (zh) * | 2010-06-12 | 2010-09-15 | 哈尔滨工业大学 | 编队卫星姿态协同控制仿真装置及其仿真验证方法 |
CN101853028A (zh) * | 2010-06-12 | 2010-10-06 | 哈尔滨工业大学 | 一种验证卫星编队的三轴气浮台导引方法 |
CN101979277A (zh) * | 2010-09-01 | 2011-02-23 | 南京航空航天大学 | 卫星磁测磁控系统的全实物验证平台与工作方法 |
CN102040008A (zh) * | 2010-12-13 | 2011-05-04 | 北京航空航天大学 | 一种用于编队卫星在轨运行安全的防碰控制方法 |
CN101694587B (zh) * | 2009-10-13 | 2011-08-17 | 清华大学 | 一种应用于微小型飞行器的集群控制系统 |
CN102541070A (zh) * | 2012-01-20 | 2012-07-04 | 哈尔滨工业大学 | 一种卫星编队飞行地面试验系统的碰撞规避方法 |
CN102682481A (zh) * | 2012-05-28 | 2012-09-19 | 国家卫星气象中心 | 区域观测模式下地球几何特征信息确定方法 |
CN102819266A (zh) * | 2012-07-20 | 2012-12-12 | 航天东方红卫星有限公司 | 一种拟周期j2不变相对轨道编队飞行控制方法 |
CN103064423A (zh) * | 2012-12-11 | 2013-04-24 | 北京空间飞行器总体设计部 | 多约束多航天器飞行间距预示及碰撞规避方法 |
CN103852082A (zh) * | 2012-11-30 | 2014-06-11 | 上海航天控制工程研究所 | 一种星间测量和陀螺的姿态轨道一体化滤波估计方法 |
CN103984238A (zh) * | 2014-06-05 | 2014-08-13 | 湖南大学 | 基于虚拟弹簧阻尼网络的卫星编队构形分布式控制方法 |
CN104142686A (zh) * | 2014-07-16 | 2014-11-12 | 北京控制工程研究所 | 一种卫星自主编队飞行控制方法 |
CN104251701A (zh) * | 2013-06-26 | 2014-12-31 | 上海新跃仪表厂 | 火星大气制动轨道确定方法 |
CN104252177A (zh) * | 2013-06-27 | 2014-12-31 | 上海新跃仪表厂 | 微小卫星对地凝视抗饱和跟踪控制方法 |
CN104330971A (zh) * | 2014-10-28 | 2015-02-04 | 蔡远文 | 微小卫星群编队消耗量优化方法 |
CN105094139A (zh) * | 2015-07-24 | 2015-11-25 | 上海微小卫星工程中心 | 一种近距离航天器共面椭圆编队的椭圆短半轴控制方法 |
CN105468011A (zh) * | 2015-12-31 | 2016-04-06 | 北京理工大学 | 一种辐射开环绳系卫星编队匀速自旋展开控制方法 |
CN105912020A (zh) * | 2016-05-26 | 2016-08-31 | 北京航空航天大学 | 一种弱稳定边界区域内航天器编队飞行方法 |
CN105930305A (zh) * | 2016-04-14 | 2016-09-07 | 清华大学深圳研究生院 | 一种三脉冲交会接近制导方法 |
CN105974939A (zh) * | 2016-07-25 | 2016-09-28 | 零度智控(北京)智能科技有限公司 | 无人机编队阵型自动生成方法及装置 |
CN106019950A (zh) * | 2016-08-09 | 2016-10-12 | 中国科学院软件研究所 | 一种手机卫星自适应姿态控制方法 |
CN107187615A (zh) * | 2017-04-25 | 2017-09-22 | 西北工业大学 | 卫星分布式载荷的编队方法 |
CN107219857A (zh) * | 2017-03-23 | 2017-09-29 | 南京航空航天大学 | 一种基于三维全局人工势函数的无人机编队路径规划算法 |
CN107246883A (zh) * | 2017-08-07 | 2017-10-13 | 上海航天控制技术研究所 | 一种高精度星敏感器安装矩阵在轨实时校准方法 |
CN107298186A (zh) * | 2017-06-18 | 2017-10-27 | 湖北航天技术研究院总体设计所 | 一种基于激光载荷的同轨面卫星自主相对位置保持方法 |
CN107422741A (zh) * | 2017-04-28 | 2017-12-01 | 西北工业大学 | 基于学习的保预设性能集群飞行分布式姿态追踪控制方法 |
CN108089434A (zh) * | 2017-12-11 | 2018-05-29 | 北京控制工程研究所 | 一种基于磁强计的皮纳卫星姿态捕获方法 |
CN108279699A (zh) * | 2018-01-02 | 2018-07-13 | 东南大学 | 一种时空可变流场下飞行器的球面轨道编队跟踪控制方法 |
CN108549412A (zh) * | 2018-04-08 | 2018-09-18 | 上海微小卫星工程中心 | 一种考虑太阳角变化率和控制因子权重可变的磁控太阳捕获方法 |
CN110471432A (zh) * | 2019-07-04 | 2019-11-19 | 中国科学院电子学研究所 | 一种卫星编队构型的方法、装置及存储介质 |
CN111114832A (zh) * | 2019-12-02 | 2020-05-08 | 北京遥测技术研究所 | 一种适用于分布式遥感编队信息处理星 |
CN111762345A (zh) * | 2020-06-22 | 2020-10-13 | 南京航空航天大学 | 一种用于集群模块拼接的多功能电磁链接机构 |
CN111854764A (zh) * | 2020-07-20 | 2020-10-30 | 中国科学院微小卫星创新研究院 | 基于星间测量信息的航天器姿态确定方法及系统 |
CN112393835A (zh) * | 2020-11-03 | 2021-02-23 | 西北工业大学深圳研究院 | 一种基于扩展卡尔曼滤波的小卫星在轨推力标定方法 |
CN112591150A (zh) * | 2021-01-05 | 2021-04-02 | 成都天巡微小卫星科技有限责任公司 | 一种控制超低轨道卫星姿态的大气阻力矩补偿方法及系统 |
CN113900442A (zh) * | 2021-10-25 | 2022-01-07 | 北京航空航天大学 | 航天器编队展开重构最优控制求解方法和系统 |
-
2008
- 2008-08-20 CN CN2008100219812A patent/CN101381004B/zh not_active Expired - Fee Related
Cited By (53)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101694587B (zh) * | 2009-10-13 | 2011-08-17 | 清华大学 | 一种应用于微小型飞行器的集群控制系统 |
CN101793526A (zh) * | 2010-04-12 | 2010-08-04 | 哈尔滨工业大学 | 一种多信息融合编队航天器自主相对导航方法 |
CN101833285A (zh) * | 2010-06-12 | 2010-09-15 | 哈尔滨工业大学 | 编队卫星姿态协同控制仿真装置及其仿真验证方法 |
CN101853028A (zh) * | 2010-06-12 | 2010-10-06 | 哈尔滨工业大学 | 一种验证卫星编队的三轴气浮台导引方法 |
CN101979277B (zh) * | 2010-09-01 | 2012-11-07 | 南京航空航天大学 | 卫星磁测磁控系统的全实物验证平台与工作方法 |
CN101979277A (zh) * | 2010-09-01 | 2011-02-23 | 南京航空航天大学 | 卫星磁测磁控系统的全实物验证平台与工作方法 |
CN102040008A (zh) * | 2010-12-13 | 2011-05-04 | 北京航空航天大学 | 一种用于编队卫星在轨运行安全的防碰控制方法 |
CN102541070A (zh) * | 2012-01-20 | 2012-07-04 | 哈尔滨工业大学 | 一种卫星编队飞行地面试验系统的碰撞规避方法 |
CN102682481A (zh) * | 2012-05-28 | 2012-09-19 | 国家卫星气象中心 | 区域观测模式下地球几何特征信息确定方法 |
CN102682481B (zh) * | 2012-05-28 | 2014-08-20 | 国家卫星气象中心 | 区域观测模式下地球几何特征信息确定方法 |
CN102819266A (zh) * | 2012-07-20 | 2012-12-12 | 航天东方红卫星有限公司 | 一种拟周期j2不变相对轨道编队飞行控制方法 |
CN102819266B (zh) * | 2012-07-20 | 2015-02-11 | 航天东方红卫星有限公司 | 一种拟周期j2不变相对轨道编队飞行控制方法 |
CN103852082A (zh) * | 2012-11-30 | 2014-06-11 | 上海航天控制工程研究所 | 一种星间测量和陀螺的姿态轨道一体化滤波估计方法 |
CN103852082B (zh) * | 2012-11-30 | 2017-04-19 | 上海航天控制工程研究所 | 一种星间测量和陀螺的姿态轨道一体化滤波估计方法 |
CN103064423A (zh) * | 2012-12-11 | 2013-04-24 | 北京空间飞行器总体设计部 | 多约束多航天器飞行间距预示及碰撞规避方法 |
CN103064423B (zh) * | 2012-12-11 | 2015-07-08 | 北京空间飞行器总体设计部 | 多约束多航天器飞行间距预示及碰撞规避方法 |
CN104251701A (zh) * | 2013-06-26 | 2014-12-31 | 上海新跃仪表厂 | 火星大气制动轨道确定方法 |
CN104252177A (zh) * | 2013-06-27 | 2014-12-31 | 上海新跃仪表厂 | 微小卫星对地凝视抗饱和跟踪控制方法 |
CN103984238A (zh) * | 2014-06-05 | 2014-08-13 | 湖南大学 | 基于虚拟弹簧阻尼网络的卫星编队构形分布式控制方法 |
CN103984238B (zh) * | 2014-06-05 | 2016-08-31 | 湖南大学 | 基于虚拟弹簧阻尼网络的卫星编队构形分布式控制方法 |
CN104142686A (zh) * | 2014-07-16 | 2014-11-12 | 北京控制工程研究所 | 一种卫星自主编队飞行控制方法 |
CN104142686B (zh) * | 2014-07-16 | 2016-06-29 | 北京控制工程研究所 | 一种卫星自主编队飞行控制方法 |
CN104330971A (zh) * | 2014-10-28 | 2015-02-04 | 蔡远文 | 微小卫星群编队消耗量优化方法 |
CN105094139A (zh) * | 2015-07-24 | 2015-11-25 | 上海微小卫星工程中心 | 一种近距离航天器共面椭圆编队的椭圆短半轴控制方法 |
CN105468011A (zh) * | 2015-12-31 | 2016-04-06 | 北京理工大学 | 一种辐射开环绳系卫星编队匀速自旋展开控制方法 |
CN105468011B (zh) * | 2015-12-31 | 2018-03-02 | 北京理工大学 | 一种辐射开环绳系卫星编队匀速自旋展开控制方法 |
CN105930305A (zh) * | 2016-04-14 | 2016-09-07 | 清华大学深圳研究生院 | 一种三脉冲交会接近制导方法 |
CN105930305B (zh) * | 2016-04-14 | 2017-07-21 | 清华大学深圳研究生院 | 一种三脉冲交会接近制导方法 |
CN105912020A (zh) * | 2016-05-26 | 2016-08-31 | 北京航空航天大学 | 一种弱稳定边界区域内航天器编队飞行方法 |
CN105974939A (zh) * | 2016-07-25 | 2016-09-28 | 零度智控(北京)智能科技有限公司 | 无人机编队阵型自动生成方法及装置 |
CN105974939B (zh) * | 2016-07-25 | 2018-12-28 | 天津远度科技有限公司 | 无人机编队阵型自动生成方法及装置 |
CN106019950A (zh) * | 2016-08-09 | 2016-10-12 | 中国科学院软件研究所 | 一种手机卫星自适应姿态控制方法 |
CN106019950B (zh) * | 2016-08-09 | 2018-11-16 | 中国科学院软件研究所 | 一种手机卫星自适应姿态控制方法 |
CN107219857A (zh) * | 2017-03-23 | 2017-09-29 | 南京航空航天大学 | 一种基于三维全局人工势函数的无人机编队路径规划算法 |
CN107187615B (zh) * | 2017-04-25 | 2019-06-21 | 西北工业大学 | 卫星分布式载荷的编队方法 |
CN107187615A (zh) * | 2017-04-25 | 2017-09-22 | 西北工业大学 | 卫星分布式载荷的编队方法 |
CN107422741A (zh) * | 2017-04-28 | 2017-12-01 | 西北工业大学 | 基于学习的保预设性能集群飞行分布式姿态追踪控制方法 |
CN107422741B (zh) * | 2017-04-28 | 2019-12-27 | 西北工业大学 | 基于学习的保预设性能集群飞行分布式姿态追踪控制方法 |
CN107298186A (zh) * | 2017-06-18 | 2017-10-27 | 湖北航天技术研究院总体设计所 | 一种基于激光载荷的同轨面卫星自主相对位置保持方法 |
CN107298186B (zh) * | 2017-06-18 | 2020-05-19 | 湖北航天技术研究院总体设计所 | 一种基于激光载荷的同轨面卫星自主相对位置保持方法 |
CN107246883A (zh) * | 2017-08-07 | 2017-10-13 | 上海航天控制技术研究所 | 一种高精度星敏感器安装矩阵在轨实时校准方法 |
CN108089434A (zh) * | 2017-12-11 | 2018-05-29 | 北京控制工程研究所 | 一种基于磁强计的皮纳卫星姿态捕获方法 |
CN108279699A (zh) * | 2018-01-02 | 2018-07-13 | 东南大学 | 一种时空可变流场下飞行器的球面轨道编队跟踪控制方法 |
CN108549412A (zh) * | 2018-04-08 | 2018-09-18 | 上海微小卫星工程中心 | 一种考虑太阳角变化率和控制因子权重可变的磁控太阳捕获方法 |
CN110471432B (zh) * | 2019-07-04 | 2020-09-08 | 中国科学院电子学研究所 | 一种卫星编队构型的方法、装置及存储介质 |
CN110471432A (zh) * | 2019-07-04 | 2019-11-19 | 中国科学院电子学研究所 | 一种卫星编队构型的方法、装置及存储介质 |
CN111114832A (zh) * | 2019-12-02 | 2020-05-08 | 北京遥测技术研究所 | 一种适用于分布式遥感编队信息处理星 |
CN111114832B (zh) * | 2019-12-02 | 2021-10-01 | 北京遥测技术研究所 | 一种适用于分布式遥感编队信息处理星 |
CN111762345A (zh) * | 2020-06-22 | 2020-10-13 | 南京航空航天大学 | 一种用于集群模块拼接的多功能电磁链接机构 |
CN111854764A (zh) * | 2020-07-20 | 2020-10-30 | 中国科学院微小卫星创新研究院 | 基于星间测量信息的航天器姿态确定方法及系统 |
CN112393835A (zh) * | 2020-11-03 | 2021-02-23 | 西北工业大学深圳研究院 | 一种基于扩展卡尔曼滤波的小卫星在轨推力标定方法 |
CN112591150A (zh) * | 2021-01-05 | 2021-04-02 | 成都天巡微小卫星科技有限责任公司 | 一种控制超低轨道卫星姿态的大气阻力矩补偿方法及系统 |
CN113900442A (zh) * | 2021-10-25 | 2022-01-07 | 北京航空航天大学 | 航天器编队展开重构最优控制求解方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN101381004B (zh) | 2010-11-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101381004A (zh) | 基于大气阻力的微小卫星编队飞行控制方法及控制装置 | |
CN104527994B (zh) | 异面交叉快变轨道固定时间稳定姿态指向跟踪控制方法 | |
CN107797130B (zh) | 低轨航天器多点多参数轨道上行数据计算方法 | |
CN100451898C (zh) | 微小卫星的姿态控制方法及系统 | |
CN101554926B (zh) | 航天器的姿态控制系统及方法 | |
CN104142686B (zh) | 一种卫星自主编队飞行控制方法 | |
CN104898642A (zh) | 一种用于航天器姿态控制算法的集成测试仿真系统 | |
CN103112603B (zh) | 欠驱动高速自旋卫星建立正常姿态的方法 | |
CN102819266B (zh) | 一种拟周期j2不变相对轨道编队飞行控制方法 | |
CN110672092B (zh) | 一种降低固定翼无人机平台磁干扰的航迹生成方法 | |
CN107609267B (zh) | 一种月球有限推力多次捕获轨道实现方法 | |
CN102354123A (zh) | 一种跨平台可扩展的卫星动态仿真测试系统 | |
CN107421550A (zh) | 一种基于星间测距的地球‑Lagrange联合星座自主定轨方法 | |
CN114967453B (zh) | 一种基于神经网络的卫星东西位保协态初值估计方法 | |
CN109063380A (zh) | 一种静止轨道电推进卫星故障检测方法及位置保持方法 | |
US11845573B2 (en) | Configuration construction and attitude control method for pyramid deorbit sail | |
CN112607065B (zh) | 一种基于电推进系统的高精度相位控制方法 | |
CN104309822A (zh) | 一种基于参数优化的航天器单脉冲水滴形绕飞轨迹悬停控制方法 | |
CN103645489A (zh) | 一种航天器gnss单天线定姿方法 | |
CN110053788B (zh) | 一种考虑复杂摄动的星座长期保持控制频次估计方法 | |
CN107861386A (zh) | 一种基于角速度观测器的抗干扰姿态控制地面验证系统及其控制方法 | |
CN104635495B (zh) | 一种含模型参数不确定性平流层卫星的经向偏移分级控制方法 | |
Ye et al. | Fuzzy sliding mode control of nonparallel-ground-track imaging satellite with high precision | |
CN101852605B (zh) | 基于简化自适应滤波的磁测微小卫星姿态确定方法 | |
CN105912020A (zh) | 一种弱稳定边界区域内航天器编队飞行方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20101110 Termination date: 20120820 |