CN113720330B - 一种亚角秒级的遥感卫星高精度姿态确定设计与实现方法 - Google Patents
一种亚角秒级的遥感卫星高精度姿态确定设计与实现方法 Download PDFInfo
- Publication number
- CN113720330B CN113720330B CN202111282863.9A CN202111282863A CN113720330B CN 113720330 B CN113720330 B CN 113720330B CN 202111282863 A CN202111282863 A CN 202111282863A CN 113720330 B CN113720330 B CN 113720330B
- Authority
- CN
- China
- Prior art keywords
- precision
- attitude
- star
- gyroscope
- satellite
- 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
Images
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/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; 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/16—Navigation; 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
- G01C21/165—Navigation; 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 combined with non-inertial navigation instruments
-
- 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/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; 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/16—Navigation; 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
- G01C21/165—Navigation; 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 combined with non-inertial navigation instruments
- G01C21/1656—Navigation; 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 combined with non-inertial navigation instruments with passive imaging devices, e.g. cameras
-
- 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/20—Instruments for performing navigational calculations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/45—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Computer Networks & Wireless Communication (AREA)
- Gyroscopes (AREA)
Abstract
本发明涉及一种亚角秒级的遥感卫星高精度姿态确定设计与实现方法。在高分辨光学载荷背板上共基准固连安装高精度星相机、星敏感器、测振陀螺测量设备,在卫星平台上安装三轴光纤陀螺。上述设备联同高分辨率光学载荷一起均采用统一的高精度时间基准。先利用星敏感器和三轴光纤陀螺进行在轨实时组合定姿,定姿精度在轨实时达到角秒量级,再利用高精度星相机、测振陀螺测量数据和星敏感器进行地面的高精度组合定姿,定姿精度最终达到亚角秒级。本发明可使姿态确定的精度和稳定度同时提升一个数量级,而且细化到光学载荷每行扫描成像时每个成像像素点都可实现亚角秒定姿精度,非常适用于光学遥感卫星的无控制点的高精度几何定位需求。
Description
技术领域
本发明属于遥感卫星姿态控制技术领域,特别是涉及一种亚角秒级的遥感卫星高精度姿态确定设计与实现方法。
背景技术
高分辨率光学遥感卫星的高精度姿态确定是制约高精度几何处理与应用的关键问题。目前国内外在轨运行的高精度三轴稳定卫星中,一般均采用陀螺和星敏感器组成的高精度姿态测量系统,通过卫星的运动学规律建立状态方程,将陀螺和星敏感器的输出通过一定的信息融合算法,获取卫星的姿态信息。星敏感器工作于实时动态测量模式,其测量精度通常最高约为3角秒。然而,为了实现亚角秒级,即小于1角秒的姿态确定,传统的星敏感器已无法满足其测量精度要求。目前亚角秒级的甚高精度星敏感器技术还不成熟,并且价格昂贵、质量大、体积大、功耗大,以上各项成本的巨大代价并不适用于小卫星的应用。并且星敏感器的数据更新率最大也只有10Hz左右,不能实时测量相机指向的快速变化,需要采用新的措施进一步提高测量精度。另外,卫星平台受姿态调整、指向控制、太阳帆板调整、星上运动部件运动等因素影响引起的扰动,使星体产生的一种幅值较小的抖动响应,称为卫星平台颤振,颤振的存在会引起影像空间质量变化,其对卫星整体性能提升的影响日渐突出,对平台颤振进行直接测量显得十分必要。
发明内容
本发明针对现有技术的不足,提供一种亚角秒级的遥感卫星高精度姿态确定设计与实现方法。本发明提出了基于高精度星相机、星敏感器、测振陀螺测量设备和光纤陀螺的高精度组合定姿系统及方法,定姿精度最终可达亚角秒级。特别是采用了由三个单轴激光陀螺构成的可对高频振动进行测量的测振陀螺设备,克服了传统星敏感器只有几十Hz的数据更新率限制,不仅在几十Hz下可实现亚角秒的定姿精度,而且细化到光学载荷每行扫描成像时(一般线阵相机行频约10KHz量级)每个成像像素点都可对应实现亚角秒定姿精度。本发明相比传统方法可使最后获取的姿态确定精度和姿态确定稳定度均提升一个数量级,非常适用于光学遥感卫星的无控制点的高精度几何定位需求。
为了达到上述目的,本发明提供的技术方案是一种亚角秒级的遥感卫星高精度姿态确定设计与实现方法,包括以下步骤:
步骤1,在高分辨光学载荷背板上共基准固连安装2台高精度星相机、2台星敏感器和1套测振陀螺,在卫星平台上安装1台三轴光纤陀螺;
步骤2,将步骤1中安装的设备联同高分辨率光学载荷一起均采用统一的高精度时间基准,由GNSS秒脉冲进行高精度时间同步,精度达到微秒量级;
步骤3,利用2台星敏感器和三轴光纤陀螺的测量值进行组合定姿,得到角秒级的卫星在轨实时姿态确定精度;
步骤4,利用联合高精度星相机、测振陀螺测量数据和星敏感器进行地面的高精度组合定姿,得到亚角秒级的卫星姿态确定精度。
而且,所述步骤1中在高分辨光学载荷背板上固连共基准安装2台高精度星相机、2台星敏感器、1套测振陀螺,其中2台星敏感器正交90°安装,2台高精度星相机与星敏感器同指向安装。上述设备实现与高分辨率光学载荷的一体化设计,安装在卫星平台上三轴光纤陀螺不需要与载荷固连安装。高分辨光学相机载荷具有较高的结构刚度和稳定性,能够承载2台高精度星相机、2台星敏感器和1套测振陀螺,预留对应的安装接口,同时在背板三个侧面预留三个基准镜安装基面。在加工时,统一高精度星相机、基准棱镜及高分辨光学相机主镜安装连接面基准。
而且,所述步骤2中2台高精度星相机、2台星敏感器、1套测振陀螺和1台三轴光纤陀螺等测量设备均采用统一的高精度时间基准,由GNSS秒脉冲进行高精度时间同步。卫星计算单元从双频GNSS接收机获取时间整数秒信息,并通过双频GNSS接收机对内输出的秒脉冲进行时间对齐。计算单元内部维持粒度为1微秒的计时器,为整星计时,卫星需要进行授时单机包括2台星敏感器、2台高精度星相机、测振陀螺、三轴光纤陀螺和高分辨光学相机载荷。计算单元读取自身的时间信息对这些设备授时,对授时单机进行时间对齐的秒脉冲信号由双频GNSS接收机发出,用于时间对齐。
而且,所述步骤3中利用星敏感器、光纤陀螺等姿态敏感器件的输入,选用扩展卡尔曼滤波方法,对卫星的姿态进行实时估计,得到角秒级的卫星在轨实时姿态确定精度。
陀螺直接测量卫星星体相对于惯性系的旋转角速度,假定测量坐标系与卫星的体坐标系重合,则陀螺的测量模型为:
陀螺的误差满足:
星敏感器测量模型为:
得到系统的状态方程为:
则式(7)可改写为:
由于星敏感器直接输出三轴姿态角,则观测方程为:
式中:
式中:
式(13)和式(14)即为星敏感器和陀螺组合定姿状态方程。
而且,所述步骤4中将2台高精度星相机、2台星敏感器、三轴光纤陀螺以及测振陀螺的原始数据下传至地面,利用联合高精度星相机、测振陀螺测量数据和星敏感器进行地面的高精度组合定姿,选取扩展卡尔曼滤波方法,可达到亚角秒级定姿精度。
测振陀螺的测量模型为:
将高精度星相机的图片数据进行处理,可得到姿态四元数信息,则可以通过下式描述高精度星相机的测量模型:
组合定姿时姿态数据的选取,采用如下策略进行:当星相机数据有效时,测量输出
姿态四元数为星相机输出姿态四元数;当星相机数据无效且星敏感器数据有效时,测
量输出姿态四元数为星敏感器输出姿态四元数;则高精度星相机、测振陀螺测量数据
和星敏感器的组合定姿方程为:
高精度星相机的星图数据星点提取可以做到0.1像素的检测精度,测量精度可达0.3″,由于采用了测振陀螺测量数据,不仅可在几十Hz下实现亚角秒的定姿精度,而且细化到光学载荷每行扫描成像时每个成像像素点都可对应实现亚角秒定姿精度,至此完成亚角秒级的遥感卫星高精度姿态确定。
与现有技术相比,本发明具有如下优点:
1)本发明提出联合利用星敏感器、高精度星相机、光纤陀螺、测振陀螺等原始测量数据进行地面的高精度组合定姿方法,实现了星上高精度姿态测量设备的优化设计,可以抑制星敏感器光轴方向姿态测量精度低的不足,进一步提高姿态确定精度达到亚角秒级,相比传统方法可使姿态确定精度和姿态确定稳定度均提升一个数量级。
2)本发明采用了测振陀螺测量数据,克服了传统星敏感器只有几十Hz的数据更新率限制,不仅仅在几十Hz下可实现亚角秒的定姿精度,而且细化到光学载荷成像时每个成像像素点都可对应实现亚角秒定姿精度,这对光学遥感卫星的无控制点的高精度几何定位非常关键。
3)本发明提出一种采用2台高精度星相机的方式,实现了亚角秒级的遥感卫星高精度姿态测量,高精度星相机不做星上处理,仅保存星上原始星图数据下传地面处理,比传统方法使姿态确定精度提升一个数量级,而且比直接采用亚角秒级的甚高精度星敏感器(价格和重量等成本代价巨大,不适用于小卫星)在星上处理大大节约了卫星成本和重量资源。
4)本发明提出装配一台测振陀螺对平台颤振进行直接测量方式,并将测振陀螺测量数据与高精度星相机输出数据进行融合,提高了数据采用率,可使姿态确定稳定度比传统方法提升一个数量级,解决了传统遥感卫星扰动源不可知的缺陷,最终效果可以对含有平台颤振影像进行校正,提高影像几何精度。
5)本发明提出在高分辨光学载荷背板上共基准安装2台星敏感器正交安装,2台高精度星相机与星敏感器同指向安装,保证了星敏感器与高精度星相机观测天区一致,不仅大大降低了太阳光影响抑制的代价,同时又可以使星敏感器给高精度星相机提供一个引导数据基准,保持了很好的数据一致性,大大提高计算效率。
6)本发明提出在高分辨光学载荷背板上固连共基准安装2台高精度星相机、2台星敏感器、1套测振陀螺,实现与高分辨率光学载荷的一体化设计,这种共基座一体化安装大大抑制空间环境带来的变形误差影响,是保证遥感卫星高精度定位精度的关键。
7)本发明采用的2台高精度星相机、2台星敏感器、1套测振陀螺和1台三轴光纤陀螺等设备联同高分辨率光学载荷一起均采用统一的高精度时间基准,由GNSS秒脉冲进行高精度时间同步,精度达到微秒量级,高精度时间同步和时间基准也是保证遥感卫星高精度定位精度的关键。
附图说明
图1为本发明实施例配置示意图。
图2为本发明实施例星敏感器和高精度星相机安装示意图。
图3为本发明实施例统一的高精度时间基准示意图。
图4为本发明实施例星敏感器和三轴光纤陀螺的组合定姿流程图。
图5为本发明实施例角速度在体坐标系X轴分量估计值(星敏感器、光纤陀螺联合滤波定姿)。
图6为本发明实施例角速度在体坐标系Y轴分量估计值(星敏感器、光纤陀螺联合滤波定姿)。
图7为本发明实施例角速度在体坐标系Z轴分量估计值(星敏感器、光纤陀螺联合滤波定姿)。
图8为本发明实施例角速度在体坐标系X轴分量估计偏差(星敏感器、光纤陀螺联合滤波定姿)。
图9为本发明实施例角速度在体坐标系Y轴分量估计偏差(星敏感器、光纤陀螺联合滤波定姿)。
图10为本发明实施例角速度在体坐标系Z轴分量估计偏差(星敏感器、光纤陀螺联合滤波定姿)。
图11为本发明实施例偏航通道姿态欧拉角估计值(星敏感器、光纤陀螺联合滤波定姿)。
图12为本发明实施例俯仰通道姿态欧拉角估计值(星敏感器、光纤陀螺联合滤波定姿)。
图13为本发明实施例滚动通道姿态欧拉角估计值(星敏感器、光纤陀螺联合滤波定姿)。
图14为本发明实施例偏航通道姿态欧拉角估计偏差(星敏感器、光纤陀螺联合滤波定姿)。
图15为本发明实施例俯仰通道姿态欧拉角估计偏差(星敏感器、光纤陀螺联合滤波定姿)。
图16为本发明实施例滚动通道姿态欧拉角估计偏差(星敏感器、光纤陀螺联合滤波定姿)。
图17为本发明实施例总姿态欧拉角估计偏差(星敏感器、光纤陀螺联合滤波定姿)。
图18为本发明实施例高精度星相机、测振陀螺、星敏感器高精度组合定姿方法流程图。
图19为本发明实施例角速度在体坐标系X轴分量估计值(高精度星相机、测振陀螺、星敏感器联合滤波定姿)。
图20为本发明实施例角速度在体坐标系Y轴分量估计值(高精度星相机、测振陀螺、星敏感器联合滤波定姿)。
图21为本发明实施例角速度在体坐标系Z轴分量估计值(高精度星相机、测振陀螺、星敏感器联合滤波定姿)。
图22为本发明实施例角速度在体坐标系X轴分量估计偏差(高精度星相机、测振陀螺、星敏感器联合滤波定姿)。
图23为本发明实施例角速度在体坐标系Y轴分量估计偏差(高精度星相机、测振陀螺、星敏感器联合滤波定姿)。
图24为本发明实施例角速度在体坐标系Z轴分量估计偏差(高精度星相机、测振陀螺、星敏感器联合滤波定姿)。
图25为本发明实施例偏航通道姿态欧拉角估计值(高精度星相机、测振陀螺、星敏感器联合滤波定姿)。
图26为本发明实施例俯仰通道姿态欧拉角估计值(高精度星相机、测振陀螺、星敏感器联合滤波定姿)。
图27为本发明实施例滚动通道姿态欧拉角估计值(高精度星相机、测振陀螺、星敏感器联合滤波定姿)。
图28为本发明实施例偏航通道姿态欧拉角估计偏差(高精度星相机、测振陀螺、星敏感器联合滤波定姿)。
图29为本发明实施例俯仰通道姿态欧拉角估计偏差(高精度星相机、测振陀螺、星敏感器联合滤波定姿)。
图30为本发明实施例滚动通道姿态欧拉角估计偏差(高精度星相机、测振陀螺、星敏感器联合滤波定姿)。
图31为本发明实施例总姿态欧拉角估计偏差(高精度星相机、测振陀螺、星敏感器联合滤波定姿)。
具体实施方式
本发明提供一种亚角秒级的遥感卫星高精度姿态确定设计与实现方法,首先在高分辨光学载荷背板上共基准固连安装2台高精度星相机(精度为0.3角秒,星上仅拍摄星图不做星点提取,原始数据下传)、2台星敏感器(精度为3角秒)、1套测振陀螺(三轴,采样频率20KHz),另外配置1台3轴光纤陀螺安装在卫星平台上;然后利用2台星敏感器和光纤陀螺进行在轨实时组合定姿,用于卫星姿态控制,定姿精度在轨实时先达到角秒量级;最后利用星敏感器、高精度星相机、光纤陀螺、测振陀螺测量数据进行地面的高精度组合定姿,定姿精度最终达到亚角秒级。
下面结合附图和实施例对本发明的技术方案作进一步说明。
本发明实施例的流程包括以下步骤:
步骤1,在高分辨光学载荷背板上共基准固连安装2台高精度星相机、2台星敏感器和1套测振陀螺,在卫星平台上安装1台三轴光纤陀螺。
如图1所示,在高分辨光学载荷背板上固连共基准安装2台高精度星相机、2台星敏感器、1套测振陀螺。其中2台星敏感器正交90°安装,2台高精度星相机与星敏感器同指向安装,如图2所示。上述设备实现与高分辨率光学载荷的一体化设计,安装在卫星平台上三轴光纤陀螺不需要与载荷固连安装。高分辨光学相机载荷具有较高的结构刚度和稳定性,能够承载2台高精度星相机、2台星敏感器和1套测振陀螺,预留对应的安装接口,同时在背板三个侧面预留三个基准镜安装基面。在加工时,统一高精度星相机、基准棱镜及高分辨光学相机主镜安装连接面基准。
步骤2,将步骤1中安装的设备联同高分辨率光学载荷一起均采用统一的高精度时间基准,由GNSS秒脉冲进行高精度时间同步,精度达到微秒量级。
如图3所示,2台高精度星相机、2台星敏感器、1套测振陀螺和1台三轴光纤陀螺等测量设备均采用统一的高精度时间基准,由GNSS秒脉冲进行高精度时间同步。卫星计算单元从双频GNSS接收机获取时间整数秒信息,并通过双频GNSS接收机对内输出的秒脉冲进行时间对齐。计算单元内部维持粒度为1微秒的计时器,为整星计时。卫星需要进行授时单机包括2台星敏感器、2台高精度星相机、测振陀螺、三轴光纤陀螺和高分辨光学相机载荷。计算单元读取自身的时间信息对这些设备授时。对授时单机进行时间对齐的秒脉冲信号由双频GNSS接收机发出,用于时间对齐。
步骤3,利用2台星敏感器和三轴光纤陀螺的测量值进行组合定姿,得到角秒级的卫星在轨实时姿态确定精度。
利用星敏感器、光纤陀螺等姿态敏感器件的输入,选用扩展卡尔曼滤波方法,对卫星的姿态进行实时估计,得到角秒级的卫星在轨实时姿态确定精度。
陀螺直接测量卫星星体相对于惯性系的旋转角速度。假定测量坐标系与卫星的体坐标系重合,则陀螺的测量模型为:
陀螺的误差满足:
星敏感器测量模型为:
则式(7)可改写为:
由于星敏感器直接输出三轴姿态角,则观测方程为:
式中:
式中:
式(13)和式(14)即为星敏感器和陀螺组合定姿状态方程。
设卫星处于对地定向模式,保持三轴零姿态,设置姿态确定仿真参数如表1所示。
表1 姿态确定仿真参数
对卫星对地定向模式120s的结果进行统计可以看出,利用星敏感器和光纤陀螺组合定姿方法得到卫星在轨的姿态确定精度和稳定度分别如下表所示,姿态确定精度优于5″(3σ),达到角秒级的卫星姿态确定精度。
表2 星敏感器、陀螺联合滤波定姿的姿态确定精度
表3 星敏感器、陀螺联合滤波定姿的姿态角速度偏差
卫星的姿态变化和姿态确定误差曲线如图5~图17所示。
步骤4,利用联合高精度星相机、测振陀螺测量数据和星敏感器进行地面的高精度组合定姿,得到亚角秒级的卫星姿态确定精度。
将2台高精度星相机、2台星敏感器、三轴光纤陀螺以及测振陀螺的原始数据下传至地面,利用联合高精度星相机、测振陀螺测量数据和星敏感器进行地面的高精度组合定姿,选取扩展卡尔曼滤波方法,得到亚角秒级定姿精度。
测振陀螺的测量模型为:
将高精度星相机的图片数据进行处理,可得到姿态四元数信息,则可以通过下式描述高精度星相机的测量模型:
组合定姿时姿态数据的选取,采用如下策略进行:当星相机数据有效时,测量输出
姿态四元数为星相机输出姿态四元数;当星相机数据无效且星敏感器数据有效时,测
量输出姿态四元数为星敏感器输出姿态四元数;则高精度星相机、测振陀螺测量数据
和星敏感器的组合定姿方程为:
高精度星相机的星图数据星点提取可以做到0.1像素的检测精度,测量精度可达0.3″。由于采用了测振陀螺测量数据,不仅可在几十Hz下实现亚角秒的定姿精度,而且细化到光学载荷每行扫描成像时每个成像像素点都可对应实现亚角秒定姿精度,至此完成亚角秒级的遥感卫星高精度姿态确定的技术实现,该方法可使姿态确定精度和姿态确定稳定度均提升一个数量级。
与步骤3中仿真条件一致,仍设卫星处于对地定向模式,保持三轴零姿态。设置姿态确定仿真参数如表4所示。
表4 姿态确定仿真参数
对卫星对地定向模式120s的结果进行统计可以看出, 利用星敏感器、高精度星相机、测振陀螺组合定姿方法得到的姿态确定精度和稳定度分别如下表所示,姿态确定精度优于0.4″(3σ),达到亚角秒级的高精度姿态确定精度,可使姿态确定精度和姿态确定稳定度均提升一个数量级。
表5 高精度星相机、测振陀螺、星敏感器
联合滤波定姿的姿态确定精度
表6 高精度星相机、测振陀螺、星敏感器
联合滤波定姿的姿态角速度偏差
卫星的姿态确定误差曲线如图19~图31所示。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (6)
1.一种亚角秒级的遥感卫星高精度姿态确定设计与实现方法,其特征在于,包括如下步骤:
步骤1,在高分辨光学载荷背板上共基准固连安装高精度星相机、星敏感器和测振陀螺,在卫星平台上安装三轴光纤陀螺;
步骤2,将步骤1中安装的设备联同高分辨率光学载荷一起均采用统一的高精度时间基准,由GNSS秒脉冲进行高精度时间同步,精度达到微秒量级;
步骤3,利用星敏感器和三轴光纤陀螺的测量值进行组合定姿,得到角秒级的卫星在轨实时姿态确定精度;
步骤4,采用由三个单轴激光陀螺构成的测振陀螺设备对高频振动进行测量,并联合高精度星相机测量数据和星敏感器进行地面的高精度组合定姿,使用扩展卡尔曼滤波方法进行姿态解算,得到亚角秒级的卫星姿态确定精度。
2.如权利要求1所述的一种亚角秒级的遥感卫星高精度姿态确定设计与实现方法,其特征在于:步骤1中在高分辨光学载荷背板上固连共基准安装2台高精度星相机、2台星敏感器、1套测振陀螺,其中2台星敏感器正交90°安装,2台高精度星相机与星敏感器同指向安装,上述所有设备实现与高分辨率光学载荷的一体化设计,安装在卫星平台上三轴光纤陀螺不需要与高分辨率光学载荷固连安装;高分辨光学相机载荷具有较高的结构刚度和稳定性,能够承载2台高精度星相机、2台星敏感器和1套测振陀螺,预留对应的安装接口,同时在背板三个侧面预留三个基准镜安装基面;在加工时,统一高精度星相机、基准棱镜及高分辨光学相机主镜安装连接面基准。
3.如权利要求2所述的一种亚角秒级的遥感卫星高精度姿态确定设计与实现方法,其特征在于:步骤2中2台高精度星相机、2台星敏感器、1套测振陀螺和1台三轴光纤陀螺均采用统一的高精度时间基准,由GNSS秒脉冲进行高精度时间同步;卫星计算单元从双频GNSS接收机获取时间整数秒信息,并通过双频GNSS接收机对内输出的秒脉冲进行时间对齐;计算单元内部维持粒度为1微秒的计时器,为整星计时,卫星需要进行授时单机包括2台星敏感器、2台高精度星相机、测振陀螺、三轴光纤陀螺和高分辨光学相机载荷;计算单元读取自身的时间信息对这些设备授时,对授时单机进行时间对齐的秒脉冲信号由双频GNSS接收机发出,用于时间对齐。
4.如权利要求3所述的一种亚角秒级的遥感卫星高精度姿态确定设计与实现方法,其特征在于:步骤3中利用星敏感器、光纤陀螺的输入,选用扩展卡尔曼滤波方法,对卫星的姿态进行实时估计,得到角秒级的卫星在轨实时姿态确定精度。
5.如权利要求4所述的一种亚角秒级的遥感卫星高精度姿态确定设计与实现方法,其特征在于:步骤3中陀螺直接测量卫星星体相对于惯性系的旋转角速度,假定测量坐标系与卫星的体坐标系重合,则陀螺的测量模型为:
ω(t)=ωtrue(t)+b(t)+w(t) (1)
式中,ω为陀螺输出的测量角速度,即观测值;ωtrue为卫星相对于惯性系的旋转角速度,即真实值;b表示陀螺零偏,描述为一阶随机游走过程;w和n是不相关的零均值白噪声;t为时间;
陀螺的误差满足:
式中,σw和σn分别表示白噪声序列w和n的标准差,I3×3表示3×3单位阵,δ(t-τ)表示狄拉克函数:
星敏感器测量模型为:
qs=qtrue+ns (4)
式中,qs为星敏感器的测量输出姿态四元数;qtrue为真实姿态四元数;ns为星敏感器的等效噪声误差,噪声的方差矩阵为V=σvI3×3,σv为白噪声的标准方差,I3×3表示3×3单位矩阵;
定义姿态四元数为:q=[q0,q1,q2,q3]=[q0,q13],选取偏差姿态四元数的矢量部分和陀螺漂移量估计偏差的组合为状态量,即:
得到系统的状态方程为:
则式(7)改写为:
由于星敏感器直接输出三轴姿态角,则观测方程为:
Z=[I3×3 03×3]x+V (11)
式中,V是零均值白噪声观测误差,标准方差为σv;
对式(10)进行泰勒展开,取离散时间间隔为T,得到离散化方程:
Xk=Φk,k-1Xk-1+Ck,k-1Wk-1 (13)
式中:
对式(11)进行泰勒展开,取离散时间间隔为T,得到离散化观测方程为:
Zk=HkXk+Vk (14)
式中:
Hk=[I3×3 03×3] (15)
式(13)和式(14)即为星敏感器和陀螺组合定姿状态方程。
6.如权利要求5所述的一种亚角秒级的遥感卫星高精度姿态确定设计与实现方法,其特征在于:步骤4中测振陀螺的测量模型为:
ω′(t)=ωtrue(t)+b′(t)+w′(t) (16)
式中,ω′为测振陀螺输出的测量角速度,即观测值;ωtrue为卫星相对于惯性系的旋转角速度,即真实值;b′表示测振陀螺漂移,描述为一阶随机游走过程;w′和n′是不相关的零均值白噪声;t为时间;
将高精度星相机的图片数据进行处理,得到姿态四元数信息,则通过下式描述高精度星相机的测量模型:
qs′=qtrue+n′s (18)
式中,qs′为高精度星相机的测量输出姿态四元数;qtrue为真实姿态四元数;n′s为高精度星相机的等效噪声误差,用白噪声来描述,噪声的方差矩阵为V=σvI3×3,σv为白噪声的标准方差,I3×3表示3×3矩阵单位矩阵;
组合定姿时姿态数据的选取,采用如下策略进行:当星相机数据有效时,测量输出姿态四元数qs′为星相机输出姿态四元数;当星相机数据无效且星敏感器数据有效时,测量输出姿态四元数qs′为星敏感器输出姿态四元数;则高精度星相机、测振陀螺测量数据和星敏感器的组合定姿方程为:
X′k=Φk,k-1X′k-1+Ck,k-1W′k-1 (19)
Z′k=HkX′k+Vk (20)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111282863.9A CN113720330B (zh) | 2021-11-01 | 2021-11-01 | 一种亚角秒级的遥感卫星高精度姿态确定设计与实现方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111282863.9A CN113720330B (zh) | 2021-11-01 | 2021-11-01 | 一种亚角秒级的遥感卫星高精度姿态确定设计与实现方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113720330A CN113720330A (zh) | 2021-11-30 |
CN113720330B true CN113720330B (zh) | 2022-02-08 |
Family
ID=78686336
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111282863.9A Active CN113720330B (zh) | 2021-11-01 | 2021-11-01 | 一种亚角秒级的遥感卫星高精度姿态确定设计与实现方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113720330B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114413883B (zh) * | 2021-12-23 | 2023-09-05 | 上海航天控制技术研究所 | 卫星姿态确定精度的提升方法、存储介质和电子设备 |
CN115512239B (zh) * | 2022-10-18 | 2023-10-03 | 闽江学院 | 一种卫星平台多星敏感器姿态自适应融合方法 |
CN116015447B (zh) * | 2022-12-27 | 2024-03-15 | 长光卫星技术股份有限公司 | 一种提升激光通信指向精度安装方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20050063224A (ko) * | 2003-12-22 | 2005-06-28 | 한국항공우주연구원 | 별 감지부와 자이로스코프를 결합한 통합 항법 시스템 |
CN102865866A (zh) * | 2012-10-22 | 2013-01-09 | 哈尔滨工业大学 | 基于双星敏感器的卫星姿态确定方法及定姿误差分析方法 |
CN105486312A (zh) * | 2016-01-30 | 2016-04-13 | 武汉大学 | 一种星敏感器与高频角位移传感器组合定姿方法及系统 |
CN107101648A (zh) * | 2017-04-26 | 2017-08-29 | 武汉大学 | 基于区域网内恒星影像的星相机定标定姿方法及系统 |
CN108759818A (zh) * | 2018-04-26 | 2018-11-06 | 上海微小卫星工程中心 | 一种超高精度导星敏感器姿态确定的方法 |
CN112146663A (zh) * | 2020-09-27 | 2020-12-29 | 武汉大学 | 基于偏置星载加速度计的gnss卫星精密定轨增强方法及系统 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7739003B2 (en) * | 2006-06-20 | 2010-06-15 | Kara Whitney Johnson | Method of determining and controlling the inertial attitude of a spinning, artificial satellite and systems therefor |
TWI689705B (zh) * | 2018-10-18 | 2020-04-01 | 財團法人國家實驗研究院 | 一種衛星姿態資料融合系統及其方法 |
-
2021
- 2021-11-01 CN CN202111282863.9A patent/CN113720330B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20050063224A (ko) * | 2003-12-22 | 2005-06-28 | 한국항공우주연구원 | 별 감지부와 자이로스코프를 결합한 통합 항법 시스템 |
CN102865866A (zh) * | 2012-10-22 | 2013-01-09 | 哈尔滨工业大学 | 基于双星敏感器的卫星姿态确定方法及定姿误差分析方法 |
CN105486312A (zh) * | 2016-01-30 | 2016-04-13 | 武汉大学 | 一种星敏感器与高频角位移传感器组合定姿方法及系统 |
CN107101648A (zh) * | 2017-04-26 | 2017-08-29 | 武汉大学 | 基于区域网内恒星影像的星相机定标定姿方法及系统 |
CN108759818A (zh) * | 2018-04-26 | 2018-11-06 | 上海微小卫星工程中心 | 一种超高精度导星敏感器姿态确定的方法 |
CN112146663A (zh) * | 2020-09-27 | 2020-12-29 | 武汉大学 | 基于偏置星载加速度计的gnss卫星精密定轨增强方法及系统 |
Non-Patent Citations (3)
Title |
---|
Low-cost Stellar Sensor for Attitude Control of Small Satellites;Scibelli Sara,et.al;《2019 PHOTONICS & ELECTROMAGNETICS RESEARCH SYMPOSIUM》;20191231;第3500-3506页 * |
基于微小卫星的空间目标光学探测定位方法;高玉东等;《海军航空工程学院学报》;20161231;第31卷(第3期);第365-371页 * |
基于星相机的高精度卫星姿态计算算法研究;王飞;《电子科技》;20171231;第30卷(第8期);第13-21页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113720330A (zh) | 2021-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113720330B (zh) | 一种亚角秒级的遥感卫星高精度姿态确定设计与实现方法 | |
US8005635B2 (en) | Self-calibrated azimuth and attitude accuracy enhancing method and system (SAAAEMS) | |
CN101413800B (zh) | 导航/稳瞄一体化系统的导航、稳瞄方法 | |
CN103994763B (zh) | 一种火星车的sins/cns深组合导航系统及其实现方法 | |
US9541392B2 (en) | Surveying system and method | |
CN101788296B (zh) | 一种sins/cns深组合导航系统及其实现方法 | |
US7142983B2 (en) | Method for the processing of non-continuous atom interferometer intertial instrument measurements and continuous wide bandwidth instrument measurements with a gravity database | |
US6876926B2 (en) | Method and system for processing pulse signals within an inertial navigation system | |
CN107270893B (zh) | 面向不动产测量的杆臂、时间不同步误差估计与补偿方法 | |
WO2019071916A1 (zh) | 天线波束姿态控制方法和系统 | |
CN103674030B (zh) | 基于天文姿态基准保持的垂线偏差动态测量装置和方法 | |
CN104698486B (zh) | 一种分布式pos用数据处理计算机系统实时导航方法 | |
CN106767787A (zh) | 一种紧耦合gnss/ins组合导航装置 | |
CN110702143B (zh) | 基于李群描述的sins捷联惯性导航系统动基座快速初始对准方法 | |
US6711517B2 (en) | Hybrid inertial navigation method and device | |
CN109269471A (zh) | 一种新型gnss接收机倾斜测量系统及方法 | |
CN110285815B (zh) | 一种可在轨全程应用的微纳卫星多源信息姿态确定方法 | |
TW202041885A (zh) | 一種繞地衛星姿態資料融合系統及其方法 | |
CN103743378A (zh) | 一种管道检测器姿态检测系统 | |
CN116105730A (zh) | 基于合作目标卫星甚短弧观测的仅测角光学组合导航方法 | |
CN111207773A (zh) | 一种用于仿生偏振光导航的姿态无约束优化求解方法 | |
Cai et al. | A polar rapid transfer alignment assisted by the improved polarized-light navigation | |
CN113551669A (zh) | 基于短基线的组合导航定位方法及装置 | |
CN114509071B (zh) | 一种风洞试验模型姿态测量方法 | |
CN106404004B (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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230406 Address after: 619-2, 6th Floor, Jiajiayue Building, No. 102 Bayi Road, Wuchang District, Wuhan City, Hubei Province, 430000 Patentee after: Wuhan Yuncheng Satellite Technology Co.,Ltd. Address before: 430072 Hubei Province, Wuhan city Wuchang District of Wuhan University Luojiashan Patentee before: WUHAN University |