CN107933967A - 一种卫星转动惯量的在轨辨识方法 - Google Patents
一种卫星转动惯量的在轨辨识方法 Download PDFInfo
- Publication number
- CN107933967A CN107933967A CN201711182002.7A CN201711182002A CN107933967A CN 107933967 A CN107933967 A CN 107933967A CN 201711182002 A CN201711182002 A CN 201711182002A CN 107933967 A CN107933967 A CN 107933967A
- Authority
- CN
- China
- Prior art keywords
- inertia
- attitude maneuver
- satellite
- axis
- angular speed
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- 230000001133 acceleration Effects 0.000 claims abstract description 25
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012850 discrimination method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
- B64G1/245—Attitude control algorithms for spacecraft attitude control
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明提供一种卫星转动惯量的在轨辨识方法,涉及航天器姿态快速机动及高稳定度姿态控制领域,步骤为:(1)采集卫星当前姿态机动三轴惯性角速度和三轴控制力矩;(2)获取所述姿态机动三轴惯性角速度对应的多项式拟合公式;(3)根据所述多项式拟合公式,确定所述姿态机动三轴惯性角速度对应的姿态机动角加速度;(4)根据所述姿态机动角加速度和所述三轴控制力矩获取卫星转动惯量;(5)根据所述卫星转动惯量获取姿态机动前馈力矩并对当前姿态机动前馈力矩进行更新。该方法根据卫星当前姿态机动三轴惯性角速度和三轴控制力的在轨数据,实时更新姿态机动前馈力矩,从而可以解决现有无法实现卫星姿态快速机动和快速稳定的问题。
Description
技术领域
本发明属于航天器姿态快速机动及高稳定度姿态控制领域,涉及一种卫星转动惯量的在轨辨识方法。
背景技术
随着卫星技术的不断发展,遥感卫星的应用也越来越普遍。其中,以高分辨率对地观测和立体测绘为代表的遥感卫星,往往需要通过高精度指向控制和快速姿态机动,实现对观测地区的单景成像、区域成像、大面积成像和机动巡查。根据该卫星的任务和工作模式需求,姿轨控分***必须具备频繁快速机动和快速稳定能力。为满足卫星姿态快速机动和快速稳定的要求,采用前馈与反馈的复合控制策略是较为合理的。其中前馈控制的基本思想为按照设计的运动规律,直接驱动卫星姿态机动,以满足快速性要求。由于前馈控制回路不需要等到输出量发生变化并形成偏差才产生作用,避免了***响应的滞后效应,因此对于提高快速性是非常有利的。
然而,前馈补偿使用的前提是前馈力矩与星体转动惯量相匹配,否则将出现过补偿或者欠补偿,反而不利于机动后的快速稳定。由于目前星体转动惯量难以准确计算且随时间会发生缓慢变化,因此现有卫星转动惯量通常是根据地面收集的卫星历史数据,通过人工经验判断得到的,因此星体转动惯量的获取精度较低,进而导致与其匹配的前馈力矩的误差较大,进而导致无法实现卫星姿态快速机动和快速稳定。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供了一种卫星转动惯量的在轨辨识方法。该方法根据卫星当前姿态机动三轴惯性角速度和三轴控制力的在轨数据,实时更新姿态机动前馈力矩,从而可以解决现有无法实现卫星姿态快速机动和快速稳定的问题。
本发明的技术解决方案是:一种卫星转动惯量的在轨辨识方法,包括如下步骤:
(1)采集卫星当前姿态机动三轴惯性角速度和三轴控制力矩;
(2)获取所述姿态机动三轴惯性角速度对应的多项式拟合公式;
(3)根据所述多项式拟合公式,确定所述姿态机动三轴惯性角速度对应的姿态机动角加速度;
(4)根据所述姿态机动角加速度和所述三轴控制力矩获取卫星转动惯量;
(5)根据所述卫星转动惯量获取姿态机动前馈力矩并对当前姿态机动前馈力矩进行更新。
进一步地,根据所述姿态机动角加速度和所述三轴控制力矩获取卫星转动惯量的方法为:
根据公式进行计算,其中,为姿态机动角加速度,Isat为卫星转动惯量,为三轴控制力矩,J11-J33为卫星转动惯量对应的矩阵元素,J21与J12、J31与J13、以及J32与J23分别相等。
进一步地,根据所述卫星转动惯量获取姿态机动前馈力矩的方法为:
根据公式进行计算,其中,为姿态机动前馈力矩,为期望角加速度,为实时采集的姿态机动三轴惯性角速度。
进一步地,所述姿态机动三轴惯性角速度对应的多项式拟合公式为 为通过所述多项式拟合公式计算的姿态机动三轴惯性角速度,t为姿态机动三轴惯性角速度对应的采集时间,n为多项式拟合公式的阶数,a0-an为与阶数n对应的常数。
进一步地,根据确定的所述多项式拟合公式,确定所述姿态机动三轴惯性角速度对应的姿态机动角加速度的方法为:
对所述多项式拟合公式进行微分运算,得到所述姿态机动角加速度对应的计算公式 为姿态机动角加速度。
本发明与现有技术相比的优点在于:本发明通过实时采集卫星当前在轨数据,即采集姿态机动三轴惯性角速度和三轴控制力矩,并根据实时采集的姿态机动三轴惯性角速度和三轴控制力矩进行多项式拟合运算得到当前卫星转动惯量,从而可以提高卫星转动惯量的获取精度,此时根据卫星转动惯量获取的姿态机动前馈力矩的精度也较高,计算误差较小,进而实现了卫星姿态快速机动和快速稳定。
附图说明
图1为本发明方法的流程框图。
具体实施方式
本发明方法的流程框图如图1所示,下面对图1中本发明方法的具体步骤进行详细阐述:
(1)采集卫星当前姿态机动三轴惯性角速度和三轴控制力矩。
具体地,可以通过在卫星上配置的陀螺仪,采集不同时刻的姿态机动三轴惯性角速度。其中,为了减少后续计算误差,对采集的姿态机动三轴惯性角速度进行陀螺常值漂移补偿和滤波处理。三轴控制力矩的数值可以用三轴控制计算输出的理论值,并补偿动量轮摩擦力矩。
(2)获取所述姿态机动三轴惯性角速度对应的多项式拟合公式。
其中,所述姿态机动三轴惯性角速度对应的多项式拟合公式为 为通过所述多项式拟合公式计算的姿态机动三轴惯性角速度,t为姿态机动三轴惯性角速度对应的采集时间,n为多项式拟合公式的阶数,a0-an为与阶数n对应的常数。
需要说明的是,所述多项式拟合公式为从不同阶数多项式拟合公式中,选择的计算得到的姿态机动三轴惯性角速度与所述采集的姿态机动三轴惯性角速度之间的差值最小的多项式拟合公式。通过选择差值最小的多项式拟合公式,作为所述姿态机动三轴惯性角速度对应的多项式拟合公式,可以保证姿态机动三轴惯性角速度的计算误差最小,进而保证姿态机动三轴惯性角速度对应的姿态机动角加速度的计算精度。
(3)根据所述多项式拟合公式,确定所述姿态机动三轴惯性角速度对应的姿态机动角加速度。
具体地,对确定的所述多项式拟合公式进行微分运算,得到所述姿态机动角加速度对应的计算公式其中,姿态机动角加速度,n为从所述不同阶数多项式拟合公式中选择的姿态机动三轴惯性角速度与所述采集的姿态机动三轴惯性角速度之间的差值最小的多项式拟合公式的阶数。本发明中通过对差值最小的多项式拟合公式进行微分运算,得到所述姿态机动角加速度对应的计算公式,可以保证姿态机动角加速度的计算精度。
(4)根据所述姿态机动角加速度和所述三轴控制力矩获取卫星转动惯量。
具体地,首先建立公式其中,Isat为卫星转动惯量,T为三轴控制力矩,J11-J33为卫星转动惯量对应的矩阵元素,卫星转动惯量的矩阵元素中对角线上的元素为主惯量,对角线两侧对称的矩阵元素为惯量基,由于卫星转动惯量中的惯量基是对称的,因此J21与J12、J31与J13、以及J32与J23分别为对称矩阵元素并且分别相等,为三轴控制力矩,Tx、Ty、Tz分别为x、y、z轴的控制力矩。进一步将公式转换为的形式,此时,为卫星转动惯量,
(5)根据所述卫星转动惯量获取姿态机动前馈力矩并对当前姿态机动前馈力矩进行更新。
具体地,根据公式进行计算,其中,为姿态机动前馈力矩,为期望角加速度,期望角加速度是根据实际需要由外部输入的值,为实时采集的姿态机动三轴惯性角速度。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (5)
1.一种卫星转动惯量的在轨辨识方法,其特征在于,包括如下步骤:
(1)采集卫星当前姿态机动三轴惯性角速度和三轴控制力矩;
(2)获取所述姿态机动三轴惯性角速度对应的多项式拟合公式;
(3)根据所述多项式拟合公式,确定所述姿态机动三轴惯性角速度对应的姿态机动角加速度;
(4)根据所述姿态机动角加速度和所述三轴控制力矩获取卫星转动惯量;
(5)根据所述卫星转动惯量获取姿态机动前馈力矩并对当前姿态机动前馈力矩进行更新。
2.根据权利要求1所述的一种卫星转动惯量的在轨辨识方法,其特征在于,根据所述姿态机动角加速度和所述三轴控制力矩获取卫星转动惯量的方法为:
根据公式进行计算,其中,为姿态机动角加速度,Isat为卫星转动惯量,为三轴控制力矩,J11-J33为卫星转动惯量对应的矩阵元素,J21与J12、J31与J13、以及J32与J23分别相等。
3.根据权利要求2所述的一种卫星转动惯量的在轨辨识方法,其特征在于,根据所述卫星转动惯量获取姿态机动前馈力矩的方法为:
根据公式进行计算,其中,为姿态机动前馈力矩,为期望角加速度,为实时采集的姿态机动三轴惯性角速度。
4.根据权利要求1所述的一种卫星转动惯量的在轨辨识方法,其特征在于,所述姿态机动三轴惯性角速度对应的多项式拟合公式为 为通过所述多项式拟合公式计算的姿态机动三轴惯性角速度,t为姿态机动三轴惯性角速度对应的采集时间,n为多项式拟合公式的阶数,a0-an为与阶数n对应的常数。
5.根据权利要求4所述的一种卫星转动惯量的在轨辨识方法,其特征在于,根据确定的所述多项式拟合公式,确定所述姿态机动三轴惯性角速度对应的姿态机动角加速度的方法为:
对所述多项式拟合公式进行微分运算,得到所述姿态机动角加速度对应的计算公式 为姿态机动角加速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711182002.7A CN107933967A (zh) | 2017-11-23 | 2017-11-23 | 一种卫星转动惯量的在轨辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711182002.7A CN107933967A (zh) | 2017-11-23 | 2017-11-23 | 一种卫星转动惯量的在轨辨识方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107933967A true CN107933967A (zh) | 2018-04-20 |
Family
ID=61931007
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711182002.7A Pending CN107933967A (zh) | 2017-11-23 | 2017-11-23 | 一种卫星转动惯量的在轨辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107933967A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109632185A (zh) * | 2018-12-13 | 2019-04-16 | 上海航天控制技术研究所 | 一种大惯量转动部件动不平衡特性在轨估计方法 |
CN112326117A (zh) * | 2020-10-23 | 2021-02-05 | 北京控制工程研究所 | 一种基于cmg的航天器惯量在轨辨识方法 |
CN112684697A (zh) * | 2020-12-15 | 2021-04-20 | 上海卫星工程研究所 | 分体式卫星在轨两舱转动惯量辨识方法及*** |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3905223A (en) * | 1971-04-07 | 1975-09-16 | Karl Ludloff | Method of measuring moment of inertia |
CN102620886A (zh) * | 2012-03-27 | 2012-08-01 | 南京航空航天大学 | 两步在轨辨识组合航天器转动惯量估计方法 |
CN103457523A (zh) * | 2013-08-31 | 2013-12-18 | 西北工业大学 | 一种无刷直流电机的参数辨识方法 |
CN104290925A (zh) * | 2014-09-29 | 2015-01-21 | 北京控制工程研究所 | 一种航天器在惯性系内的角动量控制方法 |
CN104777842A (zh) * | 2015-01-06 | 2015-07-15 | 中国人民解放军装备学院 | 一种基于磁悬浮控制敏感陀螺的卫星单轴测控一体化方法 |
CN106248300A (zh) * | 2016-07-27 | 2016-12-21 | 上海航天控制技术研究所 | 基于成对推力器连续工作的卫星质心位置测量方法 |
CN106379558A (zh) * | 2016-09-09 | 2017-02-08 | 上海航天控制技术研究所 | 一种基于角加速度前馈的滑模变结构复合控制方法 |
CN107339987A (zh) * | 2017-04-21 | 2017-11-10 | 上海交通大学 | 一种基于函数迭代积分的刚体姿态解算方法 |
-
2017
- 2017-11-23 CN CN201711182002.7A patent/CN107933967A/zh active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3905223A (en) * | 1971-04-07 | 1975-09-16 | Karl Ludloff | Method of measuring moment of inertia |
CN102620886A (zh) * | 2012-03-27 | 2012-08-01 | 南京航空航天大学 | 两步在轨辨识组合航天器转动惯量估计方法 |
CN103457523A (zh) * | 2013-08-31 | 2013-12-18 | 西北工业大学 | 一种无刷直流电机的参数辨识方法 |
CN104290925A (zh) * | 2014-09-29 | 2015-01-21 | 北京控制工程研究所 | 一种航天器在惯性系内的角动量控制方法 |
CN104777842A (zh) * | 2015-01-06 | 2015-07-15 | 中国人民解放军装备学院 | 一种基于磁悬浮控制敏感陀螺的卫星单轴测控一体化方法 |
CN106248300A (zh) * | 2016-07-27 | 2016-12-21 | 上海航天控制技术研究所 | 基于成对推力器连续工作的卫星质心位置测量方法 |
CN106379558A (zh) * | 2016-09-09 | 2017-02-08 | 上海航天控制技术研究所 | 一种基于角加速度前馈的滑模变结构复合控制方法 |
CN107339987A (zh) * | 2017-04-21 | 2017-11-10 | 上海交通大学 | 一种基于函数迭代积分的刚体姿态解算方法 |
Non-Patent Citations (1)
Title |
---|
包敏凤等: "基于总体最小二乘算法的质心在轨辨识", 《中国会议》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109632185A (zh) * | 2018-12-13 | 2019-04-16 | 上海航天控制技术研究所 | 一种大惯量转动部件动不平衡特性在轨估计方法 |
CN112326117A (zh) * | 2020-10-23 | 2021-02-05 | 北京控制工程研究所 | 一种基于cmg的航天器惯量在轨辨识方法 |
CN112326117B (zh) * | 2020-10-23 | 2022-07-05 | 北京控制工程研究所 | 一种基于cmg的航天器惯量在轨辨识方法 |
CN112684697A (zh) * | 2020-12-15 | 2021-04-20 | 上海卫星工程研究所 | 分体式卫星在轨两舱转动惯量辨识方法及*** |
CN112684697B (zh) * | 2020-12-15 | 2022-07-29 | 上海卫星工程研究所 | 分体式卫星在轨两舱转动惯量辨识方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106873611B (zh) | 一种多通道线性自抗扰控制器的设计方法 | |
CN104898429B (zh) | 一种基于自抗扰控制的三旋翼姿态控制方法 | |
CN107933967A (zh) | 一种卫星转动惯量的在轨辨识方法 | |
CN110174899A (zh) | 一种基于敏捷卫星的高精度成像姿态指向控制方法 | |
CN105786036B (zh) | 一种抑制转子动不平衡扰动的控制力矩陀螺框架控制***及方法 | |
CN109823571A (zh) | 一种遥感微纳卫星的多阶段姿态控制方法 | |
CN108845588B (zh) | 一种基于非线性制导的四旋翼飞行器轨迹跟踪控制方法 | |
CN103712598B (zh) | 一种小型无人机姿态确定方法 | |
CN103955218A (zh) | 一种基于非线性控制理论的无人艇轨迹跟踪控制装置及方法 | |
CN106774374B (zh) | 一种无人机自动巡检方法及*** | |
CN104199460A (zh) | 基于图像识别跟踪的四旋翼飞行器智能侦察*** | |
CN112572835B (zh) | 一种具有姿态切换的卫星在轨角动量管理及控制方法 | |
CN104699108A (zh) | 一种多旋翼飞行器的控制分配方法 | |
CN109752954A (zh) | 一种航空相机扫描机构位置步进自抗扰控制方法 | |
WO2018166068A1 (zh) | 无人机的飞行控制方法、装置和无人机 | |
CN110941285A (zh) | 一种基于双ip核的无人机飞行控制*** | |
CN104085539A (zh) | 成像定标的姿态控制方法 | |
CN110806754A (zh) | 四旋翼飞行器姿态修正控制***及方法 | |
CN106292297B (zh) | 基于pid控制器和l1自适应控制器的姿态控制方法 | |
CN102880049B (zh) | 一种基于帆板挠性形变测量的自适应振动控制方法 | |
CN108037523A (zh) | 一种应用于无人机的电子辅助波束对准方法 | |
CN103323031B (zh) | 一种基于星敏的地平仪***误差在线补偿的方法 | |
CN102997919A (zh) | 采用隔离载体运动提升旋转式捷联惯导误差抑制效果的方法 | |
CN109445448B (zh) | 一种轮控小卫星自适应积分滑模姿态控制器 | |
CN108168545A (zh) | 一种拟牛顿法优化补偿系数的圆锥误差补偿算法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180420 |