CN110174899A - 一种基于敏捷卫星的高精度成像姿态指向控制方法 - Google Patents

一种基于敏捷卫星的高精度成像姿态指向控制方法 Download PDF

Info

Publication number
CN110174899A
CN110174899A CN201910291483.8A CN201910291483A CN110174899A CN 110174899 A CN110174899 A CN 110174899A CN 201910291483 A CN201910291483 A CN 201910291483A CN 110174899 A CN110174899 A CN 110174899A
Authority
CN
China
Prior art keywords
quick
star
attitude
time
sts
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
Application number
CN201910291483.8A
Other languages
English (en)
Other versions
CN110174899B (zh
Inventor
赵江涛
王玉爽
蒋志雄
王淑一
古海冰
关宏
吕高见
王晋鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Control Engineering
Original Assignee
Beijing Institute of Control Engineering
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Control Engineering filed Critical Beijing Institute of Control Engineering
Priority to CN201910291483.8A priority Critical patent/CN110174899B/zh
Publication of CN110174899A publication Critical patent/CN110174899A/zh
Application granted granted Critical
Publication of CN110174899B publication Critical patent/CN110174899B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/02Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
    • G01C21/025Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means with the use of startrackers
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/10Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
    • G01C21/12Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
    • G01C21/16Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
    • G05D1/08Control of attitude, i.e. control of roll, pitch, or yaw
    • G05D1/0808Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
    • G05D1/10Simultaneous control of position or course in three dimensions
    • G05D1/101Simultaneous control of position or course in three dimensions specially adapted for aircraft

Abstract

一种基于敏捷卫星的成像姿态高精度指向控制方法,首先根据指定的地面成像目标和时间,规划姿态机动过程;机动过程中计算姿态预报信息及对应的时间并按照固定频率发送至星敏感器;星敏感器根据时间和姿态数据预报天区窗口,并保持以窗口跟踪模式尝试提取姿态;随后卫星机动至目标附近并逐渐减速,星敏感器快速高精度的捕获卫星姿态信息,系统利用双星敏光轴矢量完成姿态重置;最后利用重置后的姿态与目标姿态之差进行闭环修正,完成高精度姿态指向控制。充分利用仅在小角速度下可获得的星敏高精度姿态信息补偿具有非线性和漂移变化特性的陀螺测量信息对成像指向精度的影响。

Description

一种基于敏捷卫星的高精度成像姿态指向控制方法
技术领域
本发明提出的一种基于敏捷卫星的成像姿态高精度指向控制方法,适用于具备敏捷姿态机动能力的低轨道遥感卫星对地成像的指向控制算法,是一种从系统设计角度解决成像时姿态指向精度问题的方法。
背景技术
随着遥感卫星特别是详查卫星的不断发展,敏捷机动成像是目前比较关注的技术。敏捷机动指卫星平台具备大范围快速姿态机动的能力,普遍采用控制力矩陀螺群的控制方案,进而实现热点区域快速扫描和高分辨率成像。侦测一体化的系统要求和设计理念,同时要求卫星具备高姿态指向精度。
GF-11敏捷详查遥感卫星不仅要求具备快速机动能力(机动角速度达到 2.2°/s,后继型号超过4.0°/s),同时要求具备高稳定度和高指向精度,实现对热点区域高时间分辨率和高空间分辨率观测。而当前定姿敏感器普遍由高精度星敏感器和大量程陀螺构成,但物理原理决定随角速度的增大,高精度星敏感器测量精度下降严重甚至无法提取姿态信息,大量程陀螺的非线性增加且陀螺漂移值发生变化,测量误差的角速度相关性见图6。
由于星敏的姿态测量角速度相关特性,敏捷卫星姿态机动过程普遍采用大量程陀螺完成姿态确定,仅在姿态角速度降至一定程度后才引入星敏信息进行姿态修正。机动过程中大量程陀螺的测量误差累计从而影响系统指向精度;系统使用星敏信息修正姿态并按照修正量进行反馈控制,而系统在修正期间的修正角速度降低了系统姿态稳定度。
以上分析可知,有必要提出一种面向敏捷遥感卫星的成像姿态高精度指向控制方法,在满足卫星敏捷机动要求的前提下,保证卫星成像期间的高姿态稳定度和指向精度要求。
发明内容
本发明解决的技术问题是:克服现有技术和系统方案的不足,提出一种基于敏捷卫星的成像姿态高精度指向控制方法,该方法运用于星上控制系统方案和算法,利用姿态机动至目标附近时逐渐减速而星敏感器可快速高精度的捕获卫星姿态信息的系统特性,系统利用双星敏光轴矢量完成姿态重置,并利用高精度重置的姿态与目标姿态之差快速闭环修正完成高精度姿态指向控制。
本发明的技术方案:一种基于敏捷卫星的成像姿态高精度指向控制方法,步骤如下:
(1)根据指定的地面成像目标和成像时间以及卫星当前姿态,规划姿态机动过程:根据既定规划特征选项,得到姿态机动的欧拉转轴矢量e和对应转角χm;计算得到角速度曲线的特征时刻,即机动加速开始时间匀速滑行开始时刻tm1、机动减速开始时刻tm2以及机动到位时刻tm3;其中tm1、tm2和tm3均为相对的时间;
(2)在步骤(1)得出的姿态机动规划结果基础上,使用角加速度导数正弦曲线规划机动过程各时刻目标角速度和角加速度ar,并根据角速度和角加速度信息计算前馈力矩Tr
(3)根据步骤(1)计算得到的tm1、tm2、tm3和步骤(2)计算得到的Tr,按照规划机动欧拉转轴矢量e和对应转角χm进行姿态机动,并实时根据当前惯性姿态四元数qt和惯性角速度ωt,计算预报至各星敏的惯性四元数qSTS,i和角速度ωSTS,i;i代表星敏的编号,为正整数,取1~N,其中N代表系统配备的星敏感器个数,N≥2;
(4)系统控制器按照固定频率发送步骤(3)计算的星敏惯性四元数qSTS,i和惯性角速度ωSTS,i以及对应的时间至对应的星敏感器,各个星敏感器收到该信息后在无法自主识别姿态信息时,利用系统提供的惯性四元数qSTS,i和惯性角速度ωSTS,i,保持窗口跟踪模式并实时尝试提取卫星姿态;
(5)系统按照设计的规划路径和角速度,依次经历0~tm1时间段的加速过程、tm1~tm2时间段的匀速运动过程和tm2~tm3时间段的减速过程,机动至目标姿态附近且角速度为轨道运行状态,整个机动过程使用Tr进行前馈控制并进行常规的反馈闭环控制;
(6)当姿态机动过程结束进入稳定过程,即机动时间tm>tm3后,判断系统的姿态误差和姿态角速度误差,当满足既定的阈值条件后,再判断是否至少2 个星敏提取的姿态连续3次均有效,则认为双星敏有效;若双星敏有效,则使用光轴夹角最大且有效的2个星敏的输出姿态信息qI和星敏安装信息CSB,并对每个星敏的姿态信息进行时间差补偿,使用所选星敏的光轴信息ZIT,完成双矢量姿态计算qbi_STS;否则在设计时间段内仍使用陀螺积分姿态,等待星敏有效提取卫星姿态;所述机动时间tm=当前系统时间-
(7)根据步骤(6)得到的qbi_STS进行稳定性判断,若qbi_STS连续三拍中,每一拍与前一拍的值之差均小于设定阈值,则系统将姿态重置qt=qbi_STS,系统根据qt与目标姿态的差进行闭环修正,进而完成高精度成像姿态指向控制。
步骤(1)中计算得到角速度曲线的特征时刻的具体方法如下:
计算机动过程的特征时间:加速过程时间匀速过程时间当tv≤0时,则
匀速滑行开始时刻tm1=tacc
机动减速开始时刻tm2=tacc+tv
机动到位时刻tm3=2tacc+tv
其中为机动开始时间,amax为系统参数最大角加速度,为系统参数最大运动角速度。
步骤(2)中所使用的计算机动过程各时刻目标角速度角加速度ar方法如下:
设正弦曲线频率f=2π/tm1
如果机动过程尚未开始即tm≤0,则:ar=0,χr=0;
如果机动过程处于前期加速过程即0<tm<tm1,则:
ar=0.5amax(1-cos(f·tm)),
如果机动过程处于匀速滑行过程即tm1<tm<tm2,则:
如果机动过程处于机动后期的减速过程即tm2<tm<tm3,则: ar=-0.5amax(1-cos(f(tm-tm2))),
步骤(2)中计算前馈力矩Tr的方法如下:
qr=[exsin(χr/2),ey sin(χr/2),ezsin(χr/2),cos(χr/2)]
Cr=q2dcm(qr)
CRO=CpCrCBO0
Tr=mTr·Jsriri,pst)/ΔT
其中e为欧拉转轴矢量,ex、ey和ez分别为e的XYZ三轴分量;CRO为带偏流角修正量的目标姿态矩阵,CBO0是姿态机动前轨道系到星本体系姿态转换矩阵,Cp偏流角修正量对应的转移矩阵,为偏流角速度,ωOI为轨道角速度矢量,mTr为修正系数,Js为星体转动惯量矩阵,ωri,pst为上一控制周期的ωri,△T 为控制周期,q2dcm()是姿态四元数转为姿态矩阵的函数。
步骤(3)中计算惯性角速度ωt的方法为:
确定卫星惯性角速度ωt时使用区分正负角速度的大量程陀螺测量信息分段非线性补偿算法:
其中Δθ为陀螺输出的角度增量,KPj为正角速度第j段比例系数、θPj为正角速度第j段特征点、BPj为正角速度第j段常值系数,KNj为负角速度第j段比例系数、θNj为负角速度第j段特征点、BNj为负角速度第j段常值系数,k、l和 m为参与定姿的三个陀螺的编号,bs为陀螺常值漂移矢量。
步骤(3)中计算星敏的惯性四元数qSTS,i和角速度ωSTS,i的方法为:
ωSTS,i=CSB,iωt;qSTS,i=dcm2q(CSB,i·q2dcm(qt));
其中dcm2q()是由姿态矩阵转为姿态四元数的函数,CSB,i为星敏i的安装矩阵的转置。
步骤(6)中所使用的星敏i测量信息时间差补偿方法为:
设惯性角速度为星敏i姿态对应的时刻为tSTS,i,系统时间为t,由星敏i 输出QI,i计算星敏坐标系矢量Z、Y和X在惯性系下的分量分别为ZI,i、YI,i和XI,i
计算惯性系下补偿轴角矢量将其转换至星敏坐标系得到计算欧拉转角信息并对归一化处理,由转轴矢量和欧拉转角φΔi计算得到补偿四元素qΔi
计算得到补偿后的星敏i光轴矢量ZIt,i,其中上标 T为转置。
所述步骤(6)中的设计时间段为两倍星敏局部天区姿态识别时间。
本发明的有益效果:
(1)本发明方法采用的机动过程中使用区分正负角速度的陀螺分段非线性补偿算法,将陀螺的非线性测量误差进行补偿,减小重置后修正量,利于快速稳定。
(2)本发明方法中的姿态机动规划算法,采用角加速度导数正弦规划,进一步提高系统输出力矩的平滑度,缩减系统姿态稳定时间,特别试用于敏捷机动卫星。
(3)本发明方法采用的星敏测量信息时间差补偿方法,采用陀螺测量角速度信息和时间差信息补偿星敏测量信息的滞后特性,进一步提高指向精度。
(4)本发明采用的高精度星敏测量姿态重置系统姿态的方法,充分分析系统产品的特性,在无法改变产品物理特性的基础上,通过系统能力挖掘的方法,提升系统的性能,可以扩展到其他所有敏捷卫星的成像姿态指向控制,具有较高的市场竞争力。
附图说明
图1为本发明的系统性方法流程框图;
图2为本发明中应用的姿态规划流程图;
图3为星敏双矢量定姿计算流程框图;
图4为本发明应用的控制系统原理框图;
图5为姿态融合与修正算法原理框图;
图6为机动过程中测量误差的分析图;
图7为机动过程中的角速度特征曲线图;
具体实施方式
如图1所示,为本发明基于敏捷卫星的成像姿态高精度指向控制方法流程框图,具体步骤如下:
(1)根据指定的地面成像目标和成像时间以及卫星当前姿态,规划姿态机动过程,规划流程见图2:根据既定规划特征选项,根据目标姿态和当前姿态计算姿态机动的欧拉转轴矢量e和对应转角χm(采用教科书标准方法);计算角速度曲线的特征时刻:机动加速开始时间匀速滑行开始时刻tm1、机动减速开始时刻tm2及机动到位时刻tm3;(tm1、tm2和tm3均为相对的时间)具体计算方法为:
计算机动过程的特征时间:加速过程时间匀速过程时间当tv≤0时,则匀速滑行开始时刻tm1=tacc
机动减速开始时刻tm2=tacc+tv;机动到位时刻tm3=2tacc+tv
其中为机动开始时间,amax为系统参数最大角加速度,为系统参数最大运动角速度amax由系统角动量包络决定,amax由执行机构的力矩输出能力和当前的力矩输出需求决定,二者均为是常数。
(2)使用角加速度导数正弦曲线规划机动过程各时刻目标角速度和角加速度ar,并根据角速度和角加速度信息计算前馈力矩Tr;具体计算方法为:
a.计算目标角速度和角加速度ar
正弦曲线频率f=2π/tm1
如果机动过程尚未开始即tm≤0,则:ar=0,χr=0;
如果机动过程处于前期加速过程即0<tm<tm1,则:
ar=0.5amax(1-cos(f·tm)),
如果机动过程处于匀速滑行过程即tm1<tm<tm2,则:
如果机动过程处于机动后期的减速过程即tm2<tm<tm3,则:
ar=-0.5amax(1-cos(f(tm-tm2))),
其中:tm=当前系统时间-
b.计算前馈力矩Tr
qr=[exsin(χr/2),ey sin(χr/2),ezsin(χr/2),cos(χr/2)]
Cr=q2dcm(qr)
CRO=CpCrCBO0
Tr=mTr·Jsriri,pst)/ΔT
其中e为欧拉转轴矢量,ex、ey和ez分别为e的XYZ三轴分量;CRO为带偏流角修正量的目标姿态矩阵,CBO0是姿态机动前轨道系到星本体系姿态转换矩阵,Cp偏流角修正量对应的转移矩阵,为偏流角速度,ωOI为轨道角速度矢量,mTr为修正系数,Js为星体转动惯量矩阵,ωri,pst为上一控制周期的ωri,△T 为控制周期,q2dcm()是由q值转为姿态矩阵的函数。
(3)按照规划机动欧拉转轴矢量e和对应转角χm进行姿态机动,并实时根据当前惯性姿态四元数qt和惯性角速度ωt,计算预报至各星敏的惯性四元数 qSTS,i和角速度ωSTS,i;i代表星敏的编号,为正整数,取1~N,其中N代表系统配备的星敏感器个数,不小于2;具体方法为:
a.确定陀螺的角速度
利用非线性补偿算法公式计算ωg
其中Δθ为陀螺输出的角度增量,KPj为正角速度第j段比例系数、θPj为正角速度第j段特征点、BPj为正角速度第j段常值系数,KNj为负角速度第j段比例系数、θNj为负角速度第j段特征点、BNj为负角速度第j段常值系数。陀螺的分段数量由产品特性决定,各个敏感器的比例系数和常值系数由实际测量得到。
b.确定惯性角速度
其中k、l和m为参与定姿的三个陀螺的编号,bs为陀螺常值漂移矢量
c.计算星敏预报信息ωSTS,i=CSB,iωt;qSTS,i=dcm2q(CSB,i·q2dcm(qt));
其中dcm2q()是由姿态矩阵转为q值的函数,CSB,i为星敏i的安装矩阵的转置。
(4)系统控制器按照固定频率发送星敏惯性四元数qSTS,i和惯性角速度ωSTS,i以及对应的时间至对应的星敏感器,各个星敏感器收到该信息后在无法自主识别姿态信息时,利用系统提供的惯性四元数qSTS,i和惯性角速度ωSTS,i,保持窗口跟踪模式并实时尝试提取卫星姿态;
(5)系统按照设计的规划路径和角速度,依次经历0~tm1时间段的加速过程、tm1~tm2时间段的匀速运动过程和tm2~tm3时间段的减速过程,见图7。机动至目标姿态附近且角速度为轨道运行状态,整个机动过程使用Tr进行前馈控制并进行常规的反馈闭环控制;
(6)当姿态机动过程结束进入稳定过程即机动时间tm>tm3后,判断系统的姿态误差和姿态角速度误差,当满足既定的阈值条件后,判断是否至少2个星敏提取的姿态连续3次有效;tm=当前系统时间-
当至少双星敏均有效后,使用光轴夹角最大且有效的2个星敏的输出姿态信息qI和星敏安装信息CSB,并对每个星敏的姿态信息进行时间差补偿,使用所选星敏的光轴信息ZIt,完成双矢量姿态计算qbi_STS,流程见图3;否则在设计时间段内(一般取两倍星敏局部天区姿态识别时间)仍使用陀螺积分姿态,等待星敏有效提取卫星姿态,流程见图5;对星敏的姿态信息进行时间差补偿的方法为:
设惯性角速度为星敏i姿态对应的时刻为tSTS,i,系统时间为t,同时星敏输出信息已完成光行差修正。
a.由星敏i输出QI,i计算星敏坐标系矢量Z、Y和X在惯性系下的分量ZI,i、 YI,i和XI,i(标准教科书方法)。
b.计算惯性系下补偿轴角矢量将其转换至星敏坐标系得到计算欧拉转角信息随后对归一化处理,由转轴矢量和欧拉转角φΔi用通用方法计算补偿四元素qΔi
c.由计算得到补偿后的星敏i光轴矢量ZIt,i,公式中上标T为转置。
(7)根据qbi_STS进行稳定性判断,若qbi_STS连续三拍中,每一拍与前一拍的值之差均小于设定阈值,则系统将姿态重置qt=qbi_STS,系统根据qt与目标姿态的差进行闭环修正,进而完成高精度成像姿态指向控制。
说明:本发明针对的控制系统原理框图如图4所示,其给出构成姿态控制系统的各个功能单元,其中针对敏捷卫星的姿态融合与修正算法流程框图如图5所示。

Claims (8)

1.一种基于敏捷卫星的成像姿态高精度指向控制方法,其特征在于步骤如下:
(1)根据指定的地面成像目标和成像时间以及卫星当前姿态,规划姿态机动过程:根据既定规划特征选项,得到姿态机动的欧拉转轴矢量e和对应转角χm;计算得到角速度曲线的特征时刻,即机动加速开始时间匀速滑行开始时刻tm1、机动减速开始时刻tm2以及机动到位时刻tm3;其中tm1、tm2和tm3均为相对的时间;
(2)在步骤(1)得出的姿态机动规划结果基础上,使用角加速度导数正弦曲线规划机动过程各时刻目标角速度和角加速度ar,并根据角速度和角加速度信息计算前馈力矩Tr
(3)根据步骤(1)计算得到的tm1、tm2、tm3和步骤(2)计算得到的Tr,按照规划机动欧拉转轴矢量e和对应转角χm进行姿态机动,并实时根据当前惯性姿态四元数qt和惯性角速度ωt,计算预报至各星敏的惯性四元数qSTS,i和角速度ωSTS,i;i代表星敏的编号,为正整数,取1~N,其中N代表系统配备的星敏感器个数,N≥2;
(4)系统控制器按照固定频率发送步骤(3)计算的星敏惯性四元数qSTS,i和惯性角速度ωSTS,i以及对应的时间至对应的星敏感器,各个星敏感器收到该信息后在无法自主识别姿态信息时,利用系统提供的惯性四元数qSTS,i和惯性角速度ωSTS,i,保持窗口跟踪模式并实时尝试提取卫星姿态;
(5)系统按照设计的规划路径和角速度,依次经历0~tm1时间段的加速过程、tm1~tm2时间段的匀速运动过程和tm2~tm3时间段的减速过程,机动至目标姿态附近且角速度为轨道运行状态,整个机动过程使用Tr进行前馈控制并进行常规的反馈闭环控制;
(6)当姿态机动过程结束进入稳定过程,即机动时间tm>tm3后,判断系统的姿态误差和姿态角速度误差,当满足既定的阈值条件后,再判断是否至少2个星敏提取的姿态连续3次均有效,则认为双星敏有效;若双星敏有效,则使用光轴夹角最大且有效的2个星敏的输出姿态信息qI和星敏安装信息CSB,并对每个星敏的姿态信息进行时间差补偿,使用所选星敏的光轴信息ZIT,完成双矢量姿态计算qbi_STS;否则在设计时间段内仍使用陀螺积分姿态,等待星敏有效提取卫星姿态;所述机动时间tm=当前系统时间
(7)根据步骤(6)得到的qbi_STS进行稳定性判断,若qbi_STS连续三拍中,每一拍与前一拍的值之差均小于设定阈值,则系统将姿态重置qt=qbi_STS,系统根据qt与目标姿态的差进行闭环修正,进而完成高精度成像姿态指向控制。
2.根据权利要求1所述的一种基于敏捷卫星的成像姿态高精度指向控制方法,其特征在于:步骤(1)中计算得到角速度曲线的特征时刻的具体方法如下:
计算机动过程的特征时间:加速过程时间匀速过程时间当tv≤0时,则tv=0;
匀速滑行开始时刻tm1=tacc
机动减速开始时刻tm2=tacc+tv
机动到位时刻tm3=2tacc+tv
其中为机动开始时间,amax为系统参数最大角加速度,为系统参数最大运动角速度。
3.根据权利要求1所述的一种基于敏捷卫星的成像姿态高精度指向控制方法,其特征在于:步骤(2)中所使用的计算机动过程各时刻目标角速度角加速度ar方法如下:
设正弦曲线频率f=2π/tm1
如果机动过程尚未开始即tm≤0,则:ar=0,χr=0;
如果机动过程处于前期加速过程即0<tm<tm1,则:
ar=0.5amax(1-cos(f·tm)),
如果机动过程处于匀速滑行过程即tm1<tm<tm2,则:
ar=0,
如果机动过程处于机动后期的减速过程即tm2<tm<tm3,则:
ar=-0.5amax(1-cos(f(tm-tm2))),
4.根据权利要求3所述的一种基于敏捷卫星的成像姿态高精度指向控制方法,其特征在于:步骤(2)中计算前馈力矩Tr的方法如下:
qr=[exsin(χr/2),eysin(χr/2),ezsin(χr/2),cos(χr/2)]
Cr=q2dcm(qr)
CRO=CpCrCBO0
Tr=mTr·Jsriri,pst)/ΔT
其中e为欧拉转轴矢量,ex、ey和ez分别为e的XYZ三轴分量;CRO为带偏流角修正量的目标姿态矩阵,CBO0是姿态机动前轨道系到星本体系姿态转换矩阵,Cp偏流角修正量对应的转移矩阵,为偏流角速度,ωOI为轨道角速度矢量,mTr为修正系数,Js为星体转动惯量矩阵,ωri,pst为上一控制周期的ωri,△T为控制周期,q2dcm()是姿态四元数转为姿态矩阵的函数。
5.根据权利要求1所述的一种基于敏捷卫星的成像姿态高精度指向控制方法,其特征在于:步骤(3)中计算惯性角速度ωt的方法为:
确定卫星惯性角速度ωt时使用区分正负角速度的大量程陀螺测量信息分段非线性补偿算法:
其中Δθ为陀螺输出的角度增量,KPj为正角速度第j段比例系数、θPj为正角速度第j段特征点、BPj为正角速度第j段常值系数,KNj为负角速度第j段比例系数、θNj为负角速度第j段特征点、BNj为负角速度第j段常值系数,k、l和m为参与定姿的三个陀螺的编号,bs为陀螺常值漂移矢量。
6.根据权利要求5所述的一种基于敏捷卫星的成像姿态高精度指向控制方法,其特征在于:步骤(3)中计算星敏的惯性四元数qSTS,i和角速度ωSTS,i的方法为:
ωSTS,i=CSB,iωt;qSTS,i=dcm2q(CSB,i·q2dcm(qt));
其中dcm2q()是由姿态矩阵转为姿态四元数的函数,CSB,i为星敏i的安装矩阵的转置。
7.根据权利要求1所述的一种基于敏捷卫星的成像姿态高精度指向控制方法,其特征在于:步骤(6)中所使用的星敏i测量信息时间差补偿方法为:
设惯性角速度为星敏i姿态对应的时刻为tSTS,i,系统时间为t,由星敏i输出QI,i计算星敏坐标系矢量Z、Y和X在惯性系下的分量分别为ZI,i、YI,i和XI,i
计算惯性系下补偿轴角矢量将其转换至星敏坐标系得到计算欧拉转角信息并对归一化处理,由转轴矢量和欧拉转角φΔi计算得到补偿四元素qΔi
计算得到补偿后的星敏i光轴矢量ZIt,i,其中上标T为转置。
8.根据权利要求1所述的一种基于敏捷卫星的成像姿态高精度指向控制方法,其特征在于:所述步骤(6)中的设计时间段为两倍星敏局部天区姿态识别时间。
CN201910291483.8A 2019-04-12 2019-04-12 一种基于敏捷卫星的高精度成像姿态指向控制方法 Active CN110174899B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910291483.8A CN110174899B (zh) 2019-04-12 2019-04-12 一种基于敏捷卫星的高精度成像姿态指向控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910291483.8A CN110174899B (zh) 2019-04-12 2019-04-12 一种基于敏捷卫星的高精度成像姿态指向控制方法

Publications (2)

Publication Number Publication Date
CN110174899A true CN110174899A (zh) 2019-08-27
CN110174899B CN110174899B (zh) 2021-12-07

Family

ID=67689812

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910291483.8A Active CN110174899B (zh) 2019-04-12 2019-04-12 一种基于敏捷卫星的高精度成像姿态指向控制方法

Country Status (1)

Country Link
CN (1) CN110174899B (zh)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110750110A (zh) * 2019-10-18 2020-02-04 中国科学院上海光学精密机械研究所 用于空间激光通信的指向控制装置和控制方法
CN111121822A (zh) * 2019-12-25 2020-05-08 南京先进激光技术研究院 一种利用图像识别解决星敏相机自动修正指向的方法
CN112050807A (zh) * 2020-08-03 2020-12-08 河北汉光重工有限责任公司 一种基于时间同步补偿的sins_gnss组合导航方法
CN112061425A (zh) * 2020-09-08 2020-12-11 上海航天控制技术研究所 一种敏捷小卫星星敏感器受地气光干扰的规避方法
CN112208795A (zh) * 2020-08-10 2021-01-12 吉林大学 一种敏捷卫星非沿迹曲线成像的姿态轨迹规划方法
CN112357121A (zh) * 2020-10-30 2021-02-12 哈尔滨工业大学 一种基于预置姿态的卫星入轨快速成像方法
CN112498743A (zh) * 2020-10-23 2021-03-16 中国人民解放军战略支援部队航天工程大学 一种基于前馈与反馈的卫星姿态跟踪控制器
CN112591149A (zh) * 2020-12-01 2021-04-02 上海航天控制技术研究所 周期性干扰力矩辨识及补偿方法
CN112835071A (zh) * 2020-12-31 2021-05-25 上海交通大学 一种卫星姿态机动辅助实现窄视场载荷指向标定方法
CN112985458A (zh) * 2021-02-07 2021-06-18 北京控制工程研究所 在形变加载下进行成像的星敏感器指向测量仪及方法
CN113190028A (zh) * 2021-03-31 2021-07-30 北京控制工程研究所 一种敏捷卫星指向控制方法及系统
CN113252075A (zh) * 2021-05-12 2021-08-13 中国科学院微小卫星创新研究院 一种基于轨道周期特性的星敏感器精度修正方法
CN113465581A (zh) * 2021-05-21 2021-10-01 北京控制工程研究所 一种卫星所使用的面区域扫描观测方法
CN113720350A (zh) * 2021-08-03 2021-11-30 上海卫星工程研究所 一种多头星敏感器在轨测量精度评估方法和系统
CN115276761A (zh) * 2022-06-23 2022-11-01 中国空间技术研究院 基于高精高频姿态数据的卫星辅助数据生成方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4658361A (en) * 1983-07-25 1987-04-14 Hitachi, Ltd. Method and apparatus for determining satellite attitude by using star sensor
CN101078936A (zh) * 2007-06-08 2007-11-28 北京航空航天大学 一种基于遗传最优request和gupf的高精度组合定姿方法
CN103019247A (zh) * 2011-09-27 2013-04-03 上海航天控制工程研究所 一种火星探测器无陀螺自主空间姿态机动控制方法
CN108225337A (zh) * 2017-12-28 2018-06-29 西安电子科技大学 基于sr-ukf滤波的星敏感器和陀螺组合定姿方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4658361A (en) * 1983-07-25 1987-04-14 Hitachi, Ltd. Method and apparatus for determining satellite attitude by using star sensor
CN101078936A (zh) * 2007-06-08 2007-11-28 北京航空航天大学 一种基于遗传最优request和gupf的高精度组合定姿方法
CN103019247A (zh) * 2011-09-27 2013-04-03 上海航天控制工程研究所 一种火星探测器无陀螺自主空间姿态机动控制方法
CN108225337A (zh) * 2017-12-28 2018-06-29 西安电子科技大学 基于sr-ukf滤波的星敏感器和陀螺组合定姿方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
FAN ZHIRU: "A Research of Gyro/Star-Sensor Integrated Attitude Determination Based on Particle Filter", 《2013 THIRD INTERNATIONAL CONFERENCE ON INSTRUMENTATION》 *
TAMER MEKKY: "Fast converging with high accuracy estimates of satellite attitude and orbit based on magnetometer augmented with gyro, star sensor and GPS via extended Kalman filter", 《THE EGYPTIAN JOURNAL OF REMOTE SENSING AND SPACE SCIENCES》 *
刘一武: "星敏感器测量模型及其在卫星姿态确定系统中的应用", 《宇航学报》 *
杜宁: "卫星机动过程成像的姿态规划与控制研究", 《上海航天》 *
王淑一: "敏捷卫星快速姿态机动方法研究", 《空间控制技术与应用》 *
王玉爽: "敏捷卫星控制分系统动中成像姿态机动模式的集成测试方案设计", 《空间控制技术与应用》 *

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110750110A (zh) * 2019-10-18 2020-02-04 中国科学院上海光学精密机械研究所 用于空间激光通信的指向控制装置和控制方法
CN111121822B (zh) * 2019-12-25 2021-07-20 南京先进激光技术研究院 一种利用图像识别解决星敏相机自动修正指向的方法
CN111121822A (zh) * 2019-12-25 2020-05-08 南京先进激光技术研究院 一种利用图像识别解决星敏相机自动修正指向的方法
CN112050807A (zh) * 2020-08-03 2020-12-08 河北汉光重工有限责任公司 一种基于时间同步补偿的sins_gnss组合导航方法
CN112050807B (zh) * 2020-08-03 2023-08-18 河北汉光重工有限责任公司 一种基于时间同步补偿的sins_gnss组合导航方法
CN112208795A (zh) * 2020-08-10 2021-01-12 吉林大学 一种敏捷卫星非沿迹曲线成像的姿态轨迹规划方法
CN112208795B (zh) * 2020-08-10 2021-11-02 吉林大学 一种敏捷卫星非沿迹曲线成像的姿态轨迹规划方法
CN112061425A (zh) * 2020-09-08 2020-12-11 上海航天控制技术研究所 一种敏捷小卫星星敏感器受地气光干扰的规避方法
CN112498743A (zh) * 2020-10-23 2021-03-16 中国人民解放军战略支援部队航天工程大学 一种基于前馈与反馈的卫星姿态跟踪控制器
CN112357121A (zh) * 2020-10-30 2021-02-12 哈尔滨工业大学 一种基于预置姿态的卫星入轨快速成像方法
CN112591149A (zh) * 2020-12-01 2021-04-02 上海航天控制技术研究所 周期性干扰力矩辨识及补偿方法
CN112835071A (zh) * 2020-12-31 2021-05-25 上海交通大学 一种卫星姿态机动辅助实现窄视场载荷指向标定方法
CN112985458A (zh) * 2021-02-07 2021-06-18 北京控制工程研究所 在形变加载下进行成像的星敏感器指向测量仪及方法
CN113190028A (zh) * 2021-03-31 2021-07-30 北京控制工程研究所 一种敏捷卫星指向控制方法及系统
CN113252075A (zh) * 2021-05-12 2021-08-13 中国科学院微小卫星创新研究院 一种基于轨道周期特性的星敏感器精度修正方法
CN113465581A (zh) * 2021-05-21 2021-10-01 北京控制工程研究所 一种卫星所使用的面区域扫描观测方法
CN113720350A (zh) * 2021-08-03 2021-11-30 上海卫星工程研究所 一种多头星敏感器在轨测量精度评估方法和系统
CN113720350B (zh) * 2021-08-03 2023-09-26 上海卫星工程研究所 一种多头星敏感器在轨测量精度评估方法和系统
CN115276761A (zh) * 2022-06-23 2022-11-01 中国空间技术研究院 基于高精高频姿态数据的卫星辅助数据生成方法及系统
CN115276761B (zh) * 2022-06-23 2024-04-09 中国空间技术研究院 基于高精高频姿态数据的卫星辅助数据生成方法及系统

Also Published As

Publication number Publication date
CN110174899B (zh) 2021-12-07

Similar Documents

Publication Publication Date Title
CN110174899A (zh) 一种基于敏捷卫星的高精度成像姿态指向控制方法
CN106124170B (zh) 一种基于高精度姿态信息的相机光轴指向计算方法
CN111897357B (zh) 一种卫星对地扫描的姿态跟踪控制方法
CN107600464B (zh) 利用星敏感器信息的飞轮控制捕获太阳及对日定向方法
CN104848860B (zh) 一种敏捷卫星成像过程姿态机动规划方法
WO2020103290A1 (zh) 末子级留轨应用子系统姿态控制方法
CN108344396B (zh) 一种敏捷卫星斜条带成像模式姿态计算方法
CN108253967A (zh) 用于目标相对引导的方法和装置
CN105184002B (zh) 一种数传天线指向角度的仿真分析方法
CN105157705B (zh) 一种半捷联雷达导引头视线角速度提取方法
CN107380485B (zh) 一种微小卫星大面阵广域多模凝视成像控制方法
CN109782810B (zh) 基于图像引导的视频卫星运动目标跟踪成像方法及其装置
CN109343550A (zh) 一种基于滚动时域估计的航天器角速度的估计方法
CN111605737B (zh) 一种航天器三超控制多级协同规划与敏捷机动方法
Wu et al. Sliding-mode control for staring-mode spacecraft using a disturbance observer
CN104085539A (zh) 成像定标的姿态控制方法
CN112572835A (zh) 一种具有姿态切换的卫星在轨角动量管理及控制方法
CN105241319B (zh) 一种高速自旋制导炮弹空中实时对准方法
CN106842911A (zh) 一种二维指向镜成像跟踪的图像伺服控制方法
CN109823572B (zh) 敏捷卫星姿态往复快速摆动的执行机构配置及控制方法
CN108427429B (zh) 一种考虑动态指向约束的航天器视轴机动控制方法
CN110803304A (zh) 一种卫星姿态控制系统
CN111924142B (zh) 一种基于序列图像的软着陆高精度避障方法
CN112498743A (zh) 一种基于前馈与反馈的卫星姿态跟踪控制器
CN111947654A (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