CN106767846B - 三轴稳定卫星不用陀螺的姿态获取方法和系统 - Google Patents
三轴稳定卫星不用陀螺的姿态获取方法和系统 Download PDFInfo
- Publication number
- CN106767846B CN106767846B CN201710146943.9A CN201710146943A CN106767846B CN 106767846 B CN106767846 B CN 106767846B CN 201710146943 A CN201710146943 A CN 201710146943A CN 106767846 B CN106767846 B CN 106767846B
- Authority
- CN
- China
- Prior art keywords
- satellite
- attitude
- angular speed
- angle
- sensor
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
Abstract
本发明公开了一种三轴稳定卫星不用陀螺的姿态估计方法,包含:利用卫星姿态动力学方程对卫星惯性角速度进行估计,因此可以代替测量卫星角速度的惯性基准单元;利用卫星运动学方程及星敏感器、红外地球敏感器和数字太阳敏感器等对姿态确定值进行滤波估计;同时利用星敏感器、红外地球敏感器和数字太阳敏感器的姿态确定角进行卫星姿态角的估计,并将星敏感器姿态测量值作为基准。本发明取得了不用惯性基准单元,采用星敏感器、红外地球敏感器和数字太阳敏感器等进行卫星姿态角和姿态角速度的估计,提高卫星姿轨控系统可靠性的有益效果。
Description
技术领域
本发明涉及卫星姿态确定技术领域,具体涉及一种三轴稳定卫星不用陀螺的姿态获取方法。
背景技术
高精度高稳定度三轴稳定卫星姿态确定系统一般配置高精度的惯性基准单元来测量卫星惯性角速度,同时还配置高精度的光学敏感器如星敏感器等。对于具有转速部件(陀螺)的惯性基准单元来说,很难达长寿命要求,或者需要较高的代价。而光学敏感器较易达到长寿命。从国内高精度高稳定度三轴稳定卫星在轨运行情况来看,惯性基准单元的陀螺发生故障的概率较星敏感器等光学敏感器的高。为提高姿轨控系统的可靠性和寿命,不用陀螺的姿态确定技术是一种必要的备份措施。
不用陀螺获得卫星惯性角速度的方法有多种,本发明利用卫星姿态动力学模型和多种光学敏感器的姿态信息,利用卡尔曼估计方法确定卫星角速度估计算法和姿态角估计算法,并对多敏感器的姿态信息进行了融合处理。
目前,不用陀螺的三轴卫星姿态确定与控制系统一般仅用星敏感器估计卫星姿态和姿态角速度,在星敏感器十几分钟数据无效的情况下,受限于卫星姿态动力学模型的精确度,卫星姿态会逐渐漂大,可能会影响卫星正常工作。本发明融合了多种姿态敏感器,在某一种敏感器数据无效时也不可保证卫星姿态稳定控制。
发明内容
本发明提供一种三轴稳定卫星不用陀螺的姿态获取方法和系统,可在惯性基准单元的陀螺故障情况下,利用星敏感器、红外地球敏感器和数字太阳敏感器进行卫星惯性角速度的估计,以保持卫星正常的姿态控制。
为实现上述目的,本发明提供一种三轴稳定卫星不用陀螺的姿态获取方法,其特点是,该方法包含:
采用卫星姿态动力学方程作为关于卫星惯性角速度的状态方程,采用用卫星姿态角获取值和卫星姿态运动学方程得到惯性角速度的观测方程,通过卡尔曼获取方法得到角速度获取增益系数矩阵和惯性角速度获取的误差方差阵;
采用卫星姿态运动学方程作为关于卫星姿态角的状态方程,采用星敏感器、红外地球敏感器和数字太阳敏感器的姿态确定角得到姿态观测方程,通过卡尔曼获取方法得到姿态角获取增益系数矩阵和姿态角获取的误差方差阵;
同时利用星敏感器、红外地球敏感器和数字太阳敏感器的姿态确定角进行卫星姿态角的获取,并将星敏感器姿态测量值作为基准。
上述卫星姿态动力学方程如式(1):
式(1)中,h为飞轮角动量、ω为星体惯性角速度、I为星体惯量矩阵、Td为干扰力矩;
卫星在稳态运行时,卫星姿态运动学方程可以简化如式(2):
ωi(i=x,y,z)为卫星三轴惯性角速度;为卫星三轴姿态角速度;φ,θ,ψ为卫星三轴姿态角,ω0为卫星轨道角速度;
设r为角速度观测误差,则星体惯性角速度的观测方程如式(3):
yr=ω+r (3)
式(3)中yr为星体惯性角速度的观测值。
上利要求1所述的三轴稳定卫星不用陀螺的姿态获取方法,其特征在于,所述卫星在稳态运行时,卫星姿态运动学方程中惯性角速度由惯性角速度获取值代替,如式(4):
式(4)中,是惯性角速度获取值,wi(i=x,y,z)表示运动学模型误差,主要来源于角速度获取误差;
卫星姿态角可以由姿态敏感器的输出信息得到,因此角度获取器的观测方程如式(5):
xm=x+rm (5)
其中,x为卫星姿态角,rm为姿态测量误差。
一种三轴稳定卫星不用陀螺的姿态获取系统,其特点是,该系统包含:
卫星惯性角速度获取模块,其采用卫星姿态动力学方程作为关于卫星惯性角速度的状态方程,采用卫星姿态角获取值和卫星姿态运动学方程得到惯性角速度的观测方程,通过卡尔曼获取方法得到角速度获取增益系数矩阵和惯性角速度获取的误差方差阵;
卫星姿态确定值获取模块,其采用卫星姿态运动学方程作为关于卫星姿态角的状态方程,采用星敏感器、红外地球敏感器和数字太阳敏感器的姿态确定角得到姿态观测方程,通过卡尔曼获取方法得到姿态角获取增益系数矩阵和姿态角获取的误差方差阵;
卫星姿态获取模块,其连接卫星惯性角速度获取模块和卫星姿态确定值获取模块的输出,同时利用星敏感器、红外地球敏感器和数字太阳敏感器的姿态确定角进行卫星姿态角的获取,并将星敏感器姿态测量值作为基准。
本发明三轴稳定卫星不用陀螺的姿态获取方法和系统和现有技术相比,其优点在于,本发明提出了不用陀螺仅用光学敏感器进行姿态估计的方法,在惯性基准单元陀螺故障的情况下,可以利用星上的光学敏感器资源、卫星姿态动力学和运动学方程进行卡尔曼估计,得到卫星的惯性角速度和姿态角,为卫星的姿态控制提供基础,通过对多种敏感器信息的融合可以保证系统连续有效的姿态信息,有利于提高卫星角速度的估计精度,因此可以作为卫星姿轨控系统姿态确定的一种备份方案,有利于提高卫星姿轨控系统的可靠性;
本发明取得了不用惯性基准单元,星敏感器、红外地球敏感器和数字太阳敏感器的姿态确定角可同时用于卫星姿态角估计,并将星敏感器姿态确定值作为姿态测量基准值,因此可以提高卫星姿态角测量信息的连续有效性,有利于提高卫星惯性角速度估计的精度;
本发明充分利用光学敏感器资源可以得到连续有效的姿态信息,对卫星姿态动力学模型的精确性要求降低;
本发明仅用光学敏感器进行卫星惯性角速度估计,可作为三轴稳定卫星姿轨控系统姿态确定的备份方案,可提高卫星姿轨控系统的可靠性。
附图说明
图1为本发明三轴稳定卫星不用陀螺的姿态获取方法的流程图。
具体实施方式
以下结合附图,进一步说明本发明的具体实施例。
如图1所示,本发明公开了一种三轴稳定卫星不用陀螺的姿态获取方法,该方法包含:
S1、采用卫星姿态动力学方程作为关于卫星惯性角速度的状态方程,采用用卫星姿态角获取值和卫星姿态运动学方程得到惯性角速度的观测方程,通过卡尔曼获取方法得到角速度获取增益系数矩阵和惯性角速度获取的误差方差阵。
带有大型柔性附件的卫星在姿态稳定时对卫星的耦合干扰力矩可以忽略,因此卫星的姿态动力学方程描述为刚体模型,所以上述卫星姿态动力学方程如式(1):
式(1)中,h为飞轮角动量、ω为星体惯性角速度、I为星体惯量矩阵、Td为干扰力矩;
卫星在稳态运行时,卫星姿态运动学方程可以简化如式(2):
ωi(i=x,y,z)为卫星三轴惯性角速度;为卫星三轴姿态角速
度;φ,θ,ψ为卫星三轴姿态角,ω0为卫星轨道角速度;
得到角速度估计器的状态方程如:其中B′=I-1,v为干扰力矩的随机误差。
设r为角速度观测误差,则星体惯性角速度的观测方程如式(3):
yr=ω+r (3)
式(3)中yr为星体惯性角速度的观测值。
根据敏感器的姿态确定误差可以得到观测误差的方差阵QN=E{rr'}。
利用卡尔曼估计器设计函数可以求得角速度估计器的增益系数阵K和角速度估计误差方差阵P,于是得到卫星惯性角速度估计值:
其中,
其中,Tg(k)为卫星受到的环境干扰力矩,ω0为卫星轨道角速度,Ts为计算周期,姿态估计值由姿态估计器得到。
S2、采用卫星姿态运动学方程作为关于卫星姿态角的状态方程,采用星敏感器、红外地球敏感器和数字太阳敏感器的姿态确定角得到姿态观测方程,通过卡尔曼获取方法得到姿态角获取增益系数矩阵和姿态角获取的误差方差阵。
卫星在稳态运行时,卫星姿态运动学方程中惯性角速度由惯性角速度获取值代替,如式(4):
式(4)中,是惯性角速度获取值,wi(i=x,y,z)表示运动学模型误差,主要来源于角速度获取误差。
由卫星姿态运动学方程得到姿态角度估计器的状态方程:
其中,状态量输入量
卫星姿态角可以由姿态敏感器的输出信息得到,因此角度获取器的观测方程如式(5):
xm=x+rm (5)
其中,x为卫星姿态角,rm为姿态测量误差,可近似为姿态敏感器的姿态确定误差。
利用卡尔曼估计器设计函数可以求得角度估计器的增益系数阵L和角度估计误差方差阵P′,从而得到卡尔曼估计的姿态角,如下:
其中,Lx=L(1,1),Ly=L(2,2),Lz=L(3,3)。
S3、在S1和S2的基础上,同时利用星敏感器、红外地球敏感器和数字太阳敏感器的姿态确定角进行卫星姿态角的获取,并将星敏感器姿态测量值作为基准。
星敏感器、红外地球敏感器和数字太阳敏感器均可利用时,以星敏感器测量信息作为姿态基准,则姿态估计器可以写为:
其中,Lmx、Lmx、Lmx和LST均为增益系数,LBx、LBy和LBz为姿态基准修正系数。星敏感器积分限幅值根据基准偏差的最大值及修正系数LBx、LBy和LBz确定。
本发明还公开了一种三轴稳定卫星不用陀螺的姿态获取系统,该系统包含:卫星惯性角速度获取模块、卫星姿态确定值获取模块、连接卫星惯性角速度获取模块和卫星姿态确定值获取模块输出的卫星姿态获取模块。
卫星惯性角速度获取模块用于采用卫星姿态动力学方程作为关于卫星惯性角速度的状态方程,采用卫星姿态角获取值和卫星姿态运动学方程得到惯性角速度的观测方程,通过卡尔曼获取方法得到角速度获取增益系数矩阵和惯性角速度获取的误差方差阵;
卫星姿态确定值获取模块用于采用卫星姿态运动学方程作为关于卫星姿态角的状态方程,采用星敏感器、红外地球敏感器和数字太阳敏感器的姿态确定角得到姿态观测方程,通过卡尔曼获取方法得到姿态角获取增益系数矩阵和姿态角获取的误差方差阵;
卫星姿态获取模块用于连接卫星惯性角速度获取模块和卫星姿态确定值获取模块的输出,同时利用星敏感器、红外地球敏感器和数字太阳敏感器的姿态确定角进行卫星姿态角的获取,并将星敏感器姿态测量值作为基准。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (4)
1.一种三轴稳定卫星不用陀螺的姿态获取方法,其特征在于,该方法包含:
采用卫星姿态动力学方程作为关于卫星惯性角速度的状态方程,采用卫星姿态角获取值和卫星姿态运动学方程得到惯性角速度的观测方程,通过卡尔曼获取方法得到角速度获取增益系数矩阵和惯性角速度获取的误差方差阵;
卫星在稳态运行时,卫星姿态运动学方程中惯性角速度由惯性角速度获取值代替,卫星姿态角可以由姿态敏感器的输出信息得到;
采用卫星姿态运动学方程作为关于卫星姿态角的状态方程,采用星敏感器、红外地球敏感器和数字太阳敏感器的姿态确定角得到姿态观测方程,通过卡尔曼获取方法得到姿态角获取增益系数矩阵和姿态角获取的误差方差阵;
同时利用星敏感器、红外地球敏感器和数字太阳敏感器的姿态确定角进行卫星姿态角的获取,并将星敏感器姿态测量值作为基准。
2.如权利要求1所述的三轴稳定卫星不用陀螺的姿态获取方法,其特征在于,所述卫星姿态动力学方程如式(1):
式(1)中,h为飞轮角动量、ω为星体惯性角速度、I为星体惯量矩阵、Td为干扰力矩;
卫星在稳态运行时,卫星姿态运动学方程可以简化如式(2):
ωi(i=x,y,z)为卫星三轴惯性角速度;为卫星三轴姿态角速度;φ,θ,ψ为卫星三轴姿态角,ω0为卫星轨道角速度;
设r为角速度观测误差,则星体惯性角速度的观测方程如式(3):
yr=ω+r (3)
式(3)中yr为星体惯性角速度的观测值。
3.如权利要求1所述的三轴稳定卫星不用陀螺的姿态获取方法,其特征在于,所述卫星在稳态运行时,卫星姿态运动学方程中惯性角速度由惯性角速度获取值代替,如式(4):
式(4)中,是惯性角速度获取值,wi(i=x,y,z)表示运动学模型误差,主要来源于角速度获取误差,为卫星三轴姿态角速度;φ,θ,ψ为卫星三轴姿态角,ω0为卫星轨道角速度;
卫星姿态角可以由姿态敏感器的输出信息得到,因此角度获取器的观测方程如式(5):
xm=x+rm (5)
其中,x为卫星姿态角,rm为姿态测量误差。
4.一种三轴稳定卫星不用陀螺的姿态获取系统,其特征在于,该系统包含:
卫星惯性角速度获取模块,其采用卫星姿态动力学方程作为关于卫星惯性角速度的状态方程,采用卫星姿态角获取值和卫星姿态运动学方程得到惯性角速度的观测方程,通过卡尔曼获取方法得到角速度获取增益系数矩阵和惯性角速度获取的误差方差阵;
卫星姿态确定值获取模块,其采用卫星姿态运动学方程作为关于卫星姿态角的状态方程,采用星敏感器、红外地球敏感器和数字太阳敏感器的姿态确定角得到姿态观测方程,通过卡尔曼获取方法得到姿态角获取增益系数矩阵和姿态角获取的误差方差阵;
卫星姿态获取模块,其连接卫星惯性角速度获取模块和卫星姿态确定值获取模块的输出,同时利用星敏感器、红外地球敏感器和数字太阳敏感器的姿态确定角进行卫星姿态角的获取,并将星敏感器姿态测量值作为基准。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710146943.9A CN106767846B (zh) | 2017-03-13 | 2017-03-13 | 三轴稳定卫星不用陀螺的姿态获取方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710146943.9A CN106767846B (zh) | 2017-03-13 | 2017-03-13 | 三轴稳定卫星不用陀螺的姿态获取方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106767846A CN106767846A (zh) | 2017-05-31 |
CN106767846B true CN106767846B (zh) | 2019-10-25 |
Family
ID=58962528
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710146943.9A Active CN106767846B (zh) | 2017-03-13 | 2017-03-13 | 三轴稳定卫星不用陀螺的姿态获取方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106767846B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107869999B (zh) * | 2017-10-31 | 2020-02-14 | 上海航天控制技术研究所 | 利用角速率估计信息的单套陀螺故障诊断方法 |
CN108423195B (zh) * | 2018-02-28 | 2021-10-01 | 上海航天控制技术研究所 | 火星环绕器及环绕器全姿态无陀螺对地通信链路建立方法 |
CN109625335B (zh) * | 2018-12-10 | 2021-12-07 | 上海航天控制技术研究所 | 一种基于角速度估计信息和太阳敏感器的捕获太阳方法 |
CN111288987B (zh) * | 2020-01-17 | 2021-07-13 | 北京控制工程研究所 | 一种基于软同步的姿态测量系统 |
CN113686334B (zh) * | 2021-07-07 | 2023-08-04 | 上海航天控制技术研究所 | 一种提高星敏感器和陀螺在轨联合滤波精度的方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4930085A (en) * | 1986-10-16 | 1990-05-29 | Litef Gmbh | Method for determining the heading of an aircraft |
US5308024A (en) * | 1992-07-20 | 1994-05-03 | General Electric Co. | Disturbance torque compensated three axis yaw control system |
CN101706512A (zh) * | 2009-11-25 | 2010-05-12 | 哈尔滨工业大学 | 基于星敏感器的姿态测量信息和飞轮的角动量测量信息的航天器伪速率的估计方法 |
CN102252673A (zh) * | 2011-06-03 | 2011-11-23 | 哈尔滨工业大学 | 一种星敏感器在轨光行差的修正方法 |
CN102980580A (zh) * | 2012-11-16 | 2013-03-20 | 北京理工大学 | 基于张量积多胞鲁棒h2滤波的无陀螺卫星姿态确定方法 |
CN103019247A (zh) * | 2011-09-27 | 2013-04-03 | 上海航天控制工程研究所 | 一种火星探测器无陀螺自主空间姿态机动控制方法 |
CN103019252A (zh) * | 2011-09-27 | 2013-04-03 | 上海航天控制工程研究所 | 一种火星探测器自主对地定向控制方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9643740B2 (en) * | 2013-02-15 | 2017-05-09 | The Boeing Company | Gyroless three-axis sun acquisition using sun sensor and unscented kalman filter |
-
2017
- 2017-03-13 CN CN201710146943.9A patent/CN106767846B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4930085A (en) * | 1986-10-16 | 1990-05-29 | Litef Gmbh | Method for determining the heading of an aircraft |
US5308024A (en) * | 1992-07-20 | 1994-05-03 | General Electric Co. | Disturbance torque compensated three axis yaw control system |
CN101706512A (zh) * | 2009-11-25 | 2010-05-12 | 哈尔滨工业大学 | 基于星敏感器的姿态测量信息和飞轮的角动量测量信息的航天器伪速率的估计方法 |
CN102252673A (zh) * | 2011-06-03 | 2011-11-23 | 哈尔滨工业大学 | 一种星敏感器在轨光行差的修正方法 |
CN103019247A (zh) * | 2011-09-27 | 2013-04-03 | 上海航天控制工程研究所 | 一种火星探测器无陀螺自主空间姿态机动控制方法 |
CN103019252A (zh) * | 2011-09-27 | 2013-04-03 | 上海航天控制工程研究所 | 一种火星探测器自主对地定向控制方法 |
CN102980580A (zh) * | 2012-11-16 | 2013-03-20 | 北京理工大学 | 基于张量积多胞鲁棒h2滤波的无陀螺卫星姿态确定方法 |
Non-Patent Citations (3)
Title |
---|
Gyroless attitude and rate estimation of small satellites using singular value decomposition and extended Kalman filter;Chingiz Hajiyev,et al.;《IEEE:Proceedings of the 2015 16th International Carpathian Control Conference (ICCC)》;20151231;159-164 * |
三轴稳定卫星无陀螺姿态确定算法研究;钱勇等;《飞行力学》;20011231;第19卷(第4期);45-48,53 * |
无陀螺条件下的卫星姿态确定方法;王文兴等;《数字技术与应用》;20141231;206-208 * |
Also Published As
Publication number | Publication date |
---|---|
CN106767846A (zh) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106767846B (zh) | 三轴稳定卫星不用陀螺的姿态获取方法和系统 | |
Han et al. | A novel method to integrate IMU and magnetometers in attitude and heading reference systems | |
CN109682377B (zh) | 一种基于动态步长梯度下降的姿态估计方法 | |
CN110017837B (zh) | 一种姿态抗磁干扰的组合导航方法 | |
CN108225370B (zh) | 一种运动姿态传感器的数据融合与解算方法 | |
CN102809377A (zh) | 飞行器惯性/气动模型组合导航方法 | |
CN103712598B (zh) | 一种小型无人机姿态确定方法 | |
CN107247459B (zh) | 抗干扰飞行控制方法及装置 | |
CN107063254B (zh) | 一种陀螺地磁组合的姿态解算方法 | |
CN103674021A (zh) | 基于捷联惯导与星敏感器的组合导航系统及方法 | |
CN108344413B (zh) | 一种水下滑翔器导航系统及其低精度与高精度转换方法 | |
CN101393028B (zh) | 斜置imu安装角的快速估计与补偿系统 | |
CN106979781A (zh) | 基于分布式惯性网络的高精度传递对准方法 | |
CN110285815A (zh) | 一种可在轨全程应用的微纳卫星多源信息姿态确定方法 | |
CN111340851A (zh) | 基于双目视觉与imu融合的slam方法 | |
CN108592911A (zh) | 一种四旋翼飞行器动力学模型/机载传感器组合导航方法 | |
CN107702710B (zh) | 一种多陀螺表头常值漂移实时估计方法 | |
CN107084722B (zh) | 一种用于提高惯性-地磁组合静动态综合性能的方法 | |
CN106909165B (zh) | 基于导引头多传感器的旋转导弹弹体姿态信息提取方法 | |
CN108801250B (zh) | 基于水下机器人的实时姿态获取方法及装置 | |
Guan et al. | Sensor fusion of gyroscope and accelerometer for low-cost attitude determination system | |
Bao et al. | Aerodynamic model/INS/GPS failure-tolerant navigation method for multirotor UAVs based on federated Kalman Filter | |
Zhe et al. | Adaptive complementary filtering algorithm for imu based on mems | |
Hong et al. | Application of EKF for missile attitude estimation based on “SINS/CNS” integrated guidance system | |
CN107132850B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |