CN110162855A - 遥感卫星星上旋转载荷动态精度分析及误差分配方法 - Google Patents
遥感卫星星上旋转载荷动态精度分析及误差分配方法 Download PDFInfo
- Publication number
- CN110162855A CN110162855A CN201910391484.XA CN201910391484A CN110162855A CN 110162855 A CN110162855 A CN 110162855A CN 201910391484 A CN201910391484 A CN 201910391484A CN 110162855 A CN110162855 A CN 110162855A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- star
- remote sensing
- sensing satellite
- error
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Navigation (AREA)
Abstract
遥感卫星星上旋转载荷动态精度分析及误差分配方法,它属于航天器姿态精度建模领域。本发明解决了现有方法对磁悬浮轴承连接下的遥感卫星星上旋转载荷动态精度建模分析评估的准确率低以及无法进行误差分配的问题。根据遥感卫星***的工作模式和结构组成,分析干扰来源,明确星上旋转载荷精度误差环节及误差组成,建立星上旋转载荷精度误差传递链;再根据相应理论公式,建立适用于多误差环节传递的星上旋转载荷动态精度模型,以对磁悬浮轴承连接下的遥感卫星星上旋转载荷动态精度进行准确分析评估,实现根据旋转载荷对地经纬度精度需求设计分配各误差环节误差上限的任务需求。
Description
技术领域
本发明属于航天器姿态精度建模领域,具体涉及一种遥感卫星星上旋转载荷动态精度分析及误差分配方法。
背景技术
随着空间任务需求的不断增长,星上旋转载荷的指向精度及姿态稳定度的需求越来越高。结合遥感卫星工作模式和结构特点,在考虑多源干扰影响下,星上旋转载荷指向精度及姿态稳定度受多误差环节累积作用,如:轨道控制环节误差、姿态控制环节误差、载荷自旋环节误差等。
建立一种适合于多误差环节传递的星上旋转载荷动态精度模型对旋转载荷当前姿态精度进行有效准确评估是开展先进控制方法研究工作的前期基础,能够根据旋转载荷姿态精度与对地经纬度的换算结果为控制***设计提供各环节的误差分配和验证,能为评价、比较控制***性能提供衡量指标,具有很高的学术价值和工程意义。此外,磁悬浮轴承是一种利用电磁力实现转子漂浮的非接触型轴承,具有高转速、低功耗、无磨损、无需润滑以及中高频振动抑制特性,相比于传统机械轴承,将磁悬浮轴承应用至遥感卫星作为星上旋转载荷连接部件具有广泛的应用前景。
但由于磁悬浮轴承的低刚度特性导致星上旋转载荷的连接精度变差,在多源干扰作用下,星上旋转载荷相对于遥感卫星平台产生了章动进动现象,增添了误差环节,导致利用现有方法对磁悬浮轴承连接下的遥感卫星星上旋转载荷动态精度建模分析评估的准确率低以及无法进行误差分配,因此,对磁悬浮轴承连接下的遥感卫星星上旋转载荷动态精度建模方法的研究提出了更高的要求。
发明内容
本发明的目的是为解决现有方法对磁悬浮轴承连接下的遥感卫星星上旋转载荷动态精度建模分析评估的准确率低以及无法进行误差分配的问题。
本发明为解决上述技术问题采取的技术方案是:遥感卫星星上旋转载荷动态精度分析及误差分配方法,该方法包括以下步骤:
步骤一、确定遥感卫星***的结构和工作模式,定义地心赤道坐标系oxyz、轨道坐标系ooxoyozo、旋转载荷卫星***质心坐标系osxsyszs、卫星平台本体坐标系obxbybzb和旋转载荷本体坐标系opxpypzp;
步骤二、确定遥感卫星***包含的各误差环节,各误差环节构成精度误差传递链,各误差环节共同作用、累积传递至星上旋转载荷;
步骤三、分别计算实际星上旋转载荷指向和期望星上旋转载荷指向在期望的轨道坐标系oo1xo1yo1zo1下的投影zc和zd,以zd和zc的夹角ΔΦ为指标,来估计星上旋转载荷的指向精度;
步骤四、计算星上旋转载荷相对地心赤道坐标系oxyz的角速度矢量在测量的星上旋转载荷本体坐标系op1xp1yp1zp1下投影的期望值ωd和真实值ωc,利用期望值ωd与真实值ωc之差估计星上旋转载荷三轴姿态稳定精度;
将星上旋转载荷的指向精度和三轴姿态稳定精度的评估结果换算为旋转载荷对地经纬度指标,利用经纬度指标分配及验证遥感卫星***各误差环节的误差上限。
本发明的有益效果是:本发明的一种遥感卫星星上旋转载荷动态精度分析及误差分配方法,本发明根据遥感卫星***的工作模式和结构组成,分析干扰来源,明确星上旋转载荷精度误差环节及误差组成,建立星上旋转载荷精度误差传递链;再根据姿态运动学相应理论公式,构造适用于多误差环节传递的星上旋转载荷动态精度模型,体现各环节误差对星上旋转载荷精度的影响,以提高对磁悬浮轴承连接下的遥感卫星星上旋转载荷动态精度评估分析的准确率,从而实现针对某个确定的载荷对地经纬度精度指标设计分配并验证各误差环节误差上限的任务需求。
附图说明
图1是本发明的遥感卫星星上旋转载荷动态精度分析及误差分配方法的流程图;
图2是本发明的遥感卫星***的结构和工作模式的示意图;
图3是本发明进行坐标系定义的示意图;
图4是本发明的误差组成及误差环节的示意图;
图5是星上旋转载荷指向精度误差传递链的示意图;
图5中:OO1代表期望的轨道坐标系,OO2代表测量的轨道坐标系,Os1测量的遥感卫星***质心坐标系,Op1代表测量的星上旋转载荷本体坐标系,Zd代表光轴期望指向,Z代表光轴实际指向,θ0代表轨道控制误差,θb代表遥感卫星平台章动误差,θwz代表轴承转轴章动误差,θzj代表轴承转轴进动误差,θp代表星上旋转载荷章动误差,θz代表星上旋转载荷自旋角误差,θzx代表指向误差;
图6是星上旋转载荷指向精度与对地经纬度换算的示意图;
图6中L代表星下点,o代表地心。
具体实施方式
具体实施方式一:结合图1说明本实施方式。本实施方式所述的遥感卫星星上旋转载荷动态精度分析及误差分配方法,该方法包括以下步骤:
步骤一、确定遥感卫星***的结构和工作模式,定义地心赤道坐标系oxyz、轨道坐标系ooxoyozo、旋转载荷卫星***质心坐标系osxsyszs、卫星平台本体坐标系obxbybzb和旋转载荷本体坐标系opxpypzp;
步骤二、确定遥感卫星***包含的各误差环节,各误差环节构成精度误差传递链,各误差环节共同作用、累积传递至星上旋转载荷;
步骤三、分别计算实际星上旋转载荷指向和期望星上旋转载荷指向在期望的轨道坐标系oo1xo1yo1zo1下的投影zc和zd,以zd和zc的夹角ΔΦ为指标,来估计星上旋转载荷的指向精度;
步骤四、计算星上旋转载荷相对地心赤道坐标系oxyz的角速度矢量在测量的星上旋转载荷本体坐标系op1xp1yp1zp1下投影的期望值ωd和真实值ωc,利用期望值ωd与真实值ωc之差估计星上旋转载荷三轴姿态稳定精度;
将星上旋转载荷的指向精度和三轴姿态稳定精度的评估结果换算为旋转载荷对地经纬度指标,利用经纬度指标分配及验证遥感卫星***各误差环节的误差上限。
本实施方式针对磁悬浮轴承连接下受多源干扰影响的遥感卫星星上旋转载荷动态精度的建模问题,从遥感卫星工作模式和结构特点出发,根据姿态运动学,建立星上旋转载荷指向精度及姿态稳定度模型,能够将各误差环节影响具体体现在模型中,从而为控制***设计提供各环节指标分配和指标验证,为比较控制方法性能提供评价指标。
具体实施方式二:本实施方式与具体实施方式一不同的是:所述步骤一的具体过程为:
如图2所示,遥感卫星***的结构包括遥感卫星平台子***、星上载荷子***和旋转关节;所述遥感卫星平台子***又包含:遥感卫星平台,平台三轴正交飞轮和平台双侧太阳帆板;所述星上载荷子***又包含:星上旋转载荷和载荷内飞轮;所述旋转关节为磁悬浮轴承(并配置有机械轴承和驱动电机);
遥感卫星***的工作模式:遥感卫星平台通过平台三轴正交轮控***保持对地定向,平台双侧太阳帆板单轴驱动保持对日定向,以提供星上所需能源;旋转关节配备磁悬浮轴承提供连接,预留机械轴承为应急保护方案;星上旋转载荷通过载荷内飞轮控制保持匀速旋转,载荷内飞轮角动量饱和时需磁力矩器和驱动电机提供卸载方案;
定义坐标系如下:如图3所示;
以地心为坐标原点o,建立地心赤道坐标系(惯性坐标系)oxyz,所述地心赤道坐标系的x轴在J2000地球平赤道面内由地心指向J2000时刻的平春分点,z轴为J2000地球平赤道面的法线且指向北极方向,y轴同x轴、z轴构成右手系;
以遥感卫星***质心为坐标原点oo,建立轨道坐标系ooxoyozo,所述轨道坐标系的xo轴在轨道面内并指向遥感卫星***的前进方向,zo轴由遥感卫星***质心指向地心,yo轴同xo轴、zo轴构成右手系;
记期望的轨道坐标系为oo1xo1yo1zo1,期望的轨道坐标系的坐标原点同轨道坐标系ooxoyozo的坐标原点,在没有轨道控制误差时,期望的轨道坐标系oo1xo1yo1zo1的三轴同轨道坐标系ooxoyozo的三轴重合;记测量的轨道坐标系为oo2xo2yo2zo2,测量的轨道坐标系oo2xo2yo2zo2的坐标原点同轨道坐标系ooxoyozo的坐标原点,在没有轨道确定误差时,测量的轨道坐标系oo2xo2yo2zo2的三轴同轨道坐标系ooxoyozo的三轴重合;
以遥感卫星***质心为坐标原点os,建立遥感卫星***质心坐标系osxsyszs(与体固连),所述遥感卫星***质心坐标系的xs轴指向遥感卫星平台轴向,ys轴和zs轴位于遥感卫星平台轴向截面内,且ys轴、zs轴同xs轴构成右手系;在无姿态控制误差下,遥感卫星***质心坐标系的三轴同测量的轨道坐标系oo2xo2yo2zo2的三轴重合;
记测量的遥感卫星***质心坐标系为os1xs1ys1zs1,且测量的遥感卫星***质心坐标系os1xs1ys1zs1的坐标原点同遥感卫星***质心坐标系osxsyszs的坐标原点,在无姿态测量误差时,测量的遥感卫星***质心坐标系os1xs1ys1zs1的三轴同遥感卫星***质心坐标系osxsyszs的三轴重合;
以遥感卫星平台质心为坐标原点ob,建立遥感卫星平台本体坐标系obxbybzb(与体固连),遥感卫星平台本体坐标系的三轴指向同遥感卫星***质心坐标系osxsyszs的三轴指向;
记测量的遥感卫星平台本体坐标系为ob1xb1yb1zb1,且测量的遥感卫星平台本体坐标系ob1xb1yb1zb1的坐标原点同遥感卫星平台本体坐标系obxbybzb的坐标原点,测量的遥感卫星平台本体坐标系ob1xb1yb1zb1的三轴指向同测量的遥感卫星***质心坐标系os1xs1ys1zs1的三轴指向;
以星上旋转载荷质心为坐标原点op,建立星上旋转载荷本体坐标系opxpypzp(与体固连),在旋转关节连接误差为零、星上旋转载荷自旋转角为零时,星上旋转载荷本体坐标系opxpypzp的三轴指向同测量的遥感卫星***质心坐标系os1xs1ys1zs1的三轴指向;
记测量的星上旋转载荷本体坐标系为op1xp1yp1zp1,且测量的星上旋转载荷本体坐标系op1xp1yp1zp1的坐标原点同遥感卫星平台本体坐标系opxpypzp的坐标原点,在无测量误差时,测量的星上旋转载荷本体坐标系op1xp1yp1zp1的三轴同坐标系opxpypzp的三轴重合。
对星上旋转载荷的干扰来源分析如下:
根据遥感卫星工作模式和结构特点,干扰种类分为轨道干扰、姿态干扰和自旋干扰三类。
轨道干扰包含:在轨运行所受空间干扰力,影响轨道控制精度;轨道摄动干扰,影响轨道控制精度,具体有地球引力场摄动干扰、太阳引力场摄动干扰、太阳光压摄动干扰、月球引力场摄动干扰。
姿态干扰包含:载荷动静不平衡干扰,在旋转载荷生产、安装过程中不可避免地会有质心偏差和主轴偏转,即有残余动静不平衡量,且该不平衡量尚未被现有质量配平技术所抵消,在载荷自旋时不平衡量产生了不平衡干扰力矩,经磁悬浮轴承传递至遥感卫星平台影响姿态控制精度;
模型参数时变干扰,因载荷径向主惯量差异导致整星转动惯量随载荷旋转变化,体现模型参数时变性,产生姿态干扰力矩,经磁悬浮轴承传递至遥感卫星平台影响姿态控制精度;
载荷轴向偏置角动量干扰,在载荷飞轮未完全补偿载荷轴向角动量时,载荷轴向剩余角动量产生姿态干扰力矩,经磁悬浮轴承传递至遥感卫星平台影响姿态控制精度;
卫星平台振动干扰,遥感卫星姿态运动和太阳帆板驱动旋转激发太阳帆板多频次、宽频带的弹性振动,遥感卫星平台内三轴正交飞轮在考虑飞轮质心偏移时,也会产生中高频振动干扰,影响姿态控制精度。
空间环境干扰力矩,作用至遥感卫星***,影响姿态控制精度,具体考虑重力梯度力矩、地磁力矩、太阳辐射力矩干扰影响。
其中,为了隔离卫星平台中高频振动传递路径,磁悬浮轴承需控制在较低刚度量级,导致在载荷动静不平衡干扰和模型参数时变干扰作用下,载荷相对平台同样产生章动进动现象,影响载荷连接精度。
自旋干扰:磁轴承摩擦阻力干扰,阻碍轴承转轴自旋,影响载荷自旋精度;磁轴承磁阻干扰力矩,轴承转轴在磁悬浮轴承磁场中运动所受到的磁场阻力干扰,影响载荷自旋精度;电机齿槽力矩干扰,驱动电机不加电,由于永磁极影响,驱动电机存在齿槽力矩干扰影响载荷自旋精度;驱动电机电磁阻力矩,轴承转轴在电机磁场中运动所受到的磁场阻力干扰,影响载荷自旋精度。
具体实施方式三:结合图4说明本实施方式。本实施方式与具体实施方式二不同的是:所述步骤二的具体过程为:
确定遥感卫星***包含的各误差环节为:轨道误差环节、平台姿态误差环节、磁悬浮轴承连接误差环节和星上旋转载荷自旋误差环节;
轨道误差环节由两部分构成:一部分是轨道的控制误差,受轨道干扰影响,造成期望的轨道坐标系oo1xo1yo1zo1和轨道坐标系ooxoyozo之间有相对偏差;另一部分是轨道的确定误差,受轨道模型误差、星时误差、初值误差以及计算误差影响,造成轨道坐标系ooxoyozo和测量的轨道坐标系oo2xo2yo2zo2之间有相对偏差;
平台姿态误差环节由两部分构成:一部分是姿态的控制误差,受姿态干扰影响,造成遥感卫星平台本体系obxbybzb和测量的轨道坐标系oo2xo2yo2zo2之间有相对偏差;另一部分是姿态测量误差,包含星敏感器误差和光纤陀螺误差,造成遥感卫星平台本体系obxbybzb和测量的遥感卫星平台本体系ob1xb1yb1zb1之间有相对偏差,其中测量元件的误差由内部校准误差,空间环境干扰,测量元件老化,噪声干扰,测量元件安装误差构成;
磁悬浮轴承连接误差环节由两部分构成:一部分是轴隙跳动误差,导致星上旋转载荷本体坐标系opxpypzp相对测量的遥感卫星平台本体系ob1xb1yb1zb1存在径向章动进动姿态偏差;另一部分是测量误差,受测量元件影响,导致星上旋转载荷本体坐标系opxpypzp和测量的星上旋转载荷本体坐标系op1xp1yp1zp1之间存在径向姿态偏差;
星上旋转载荷自旋误差环节由两部分构成,一部分是转速控制误差,导致星上旋转载荷本体坐标系opxpypzp相对测量的遥感卫星平台本体系ob1xb1yb1zb1存在轴向转速误差;另一部分是测量误差,受微测陀螺和星敏感器测量误差影响,导致星上旋转载荷本体坐标系opxpypzp和测量的星上旋转载荷本体坐标系op1xp1yp1zp1之间存在轴向转速偏差;
轨道误差环节、平台姿态误差环节、磁悬浮轴承连接误差环节和星上旋转载荷自旋误差环节串联在一起构成精度误差传递链,各误差环节共同作用、累积传递至星上旋转载荷。体现为星上旋转载荷姿态指向误差及三轴姿态稳定度误差。
具体实施方式四:结合图5说明本实施方式。本实施方式与具体实施方式三不同的是:所述步骤三的具体过程为:
轨道误差环节导致轨道坐标系ooxoyozo相对期望的轨道坐标系oo1xo1yo1zo1存在误差,利用zxy顺序描述轨道坐标系ooxoyozo相对期望的轨道坐标系oo1xo1yo1zo1的姿态误差,所转过的欧拉角记为:欧拉角的定义为:期望的轨道坐标系oo1xo1yo1zo1绕zo1轴转过角度至中间系i1,中间系i1绕自身(即中间系i1)的x轴转过角度Δθ1至中间系i2,中间系i2绕轨道坐标系的yo轴转过Δψ1至轨道坐标系ooxoyozo;Δθ1和Δψ1均代表角度误差,没有实际的物理含义;
卫星姿态控制所用的测量轨道同真实轨道之间存在误差,即测量的轨道坐标系oo2xo2yo2zo2相对轨道坐标系ooxoyozo存在误差,利用zxy顺序描述测量的轨道坐标系oo2xo2yo2zo2相对轨道坐标系ooxoyozo的姿态误差,所转过的欧拉角记为:欧拉角的定义为:轨道坐标系ooxoyozo绕zo轴转过角度至中间系i1′,中间系i1′绕自身(即中间系i1′)的x轴转过角度Δθ2至中间系i2′,中间系i2′绕测量的轨道坐标系的yo2轴转过角度Δψ2至测量的轨道坐标系oo2xo2yo2zo2;Δθ2和Δψ2均代表角度误差;没有实际的物理含义;
遥感卫星平台控制结合姿态测量误差导致测量的遥感卫星平台本体系ob1xb1yb1zb1相对测量的轨道坐标系oo2xo2yo2zo2存在误差,利用zxy顺序描述遥感卫星平台本体系ob1xb1yb1zb1相对测量的轨道坐标系oo2xo2yo2zo2的姿态误差,所转过的欧拉角记为:欧拉角的定义为:测量的轨道坐标系oo2xo2yo2zo2绕zo2轴转过角度至中间系i1″,中间系i1″绕自身(即中间系i1″)的x轴转过角度Δθ3至中间系i2″,中间系i2″绕遥感卫星平台本体系的yb1轴转过角度Δψ3至遥感卫星平台本体系ob1xb1yb1zb1;Δθ3和Δψ3分别代表偏航角、滚转角和俯仰角;
磁悬浮轴承导致轴承转轴相对轴承定子之间存在姿态偏差,轴承转轴和星上旋转载荷固连,轴承定子和遥感卫星平台固连,再结合星上旋转载荷自旋精度和姿态测量误差,测量的星上旋转载荷本体坐标系op1xp1yp1zp1相对测量的遥感卫星平台本体系ob1xb1yb1zb1存在误差,利用zyx描述测量的星上旋转载荷本体坐标系op1xp1yp1zp1相对测量的遥感卫星平台本体系ob1xb1yb1zb1的姿态误差,所转过的欧拉角记为:欧拉角的定义为:测量的遥感卫星平台本体系ob1xb1yb1zb1绕zb1轴转过角度Δψ4至中间系i1″′,中间系i1″′绕自身的y轴转过角度Δθ4至中间系i2″′,中间系i2″′绕测量的星上旋转载荷本体坐标系的xp1轴转过角度至测量的星上旋转载荷本体坐标系op1xp1yp1zp1;Δψ4和Δθ4均代表角度误差,代表星上旋转载荷相对遥感卫星平台的期望转角,代表星上旋转载荷相对遥感卫星平台的转角偏差;
分别计算实际星上旋转载荷指向和期望星上旋转载荷指向在期望的轨道坐标系oo1xo1yo1zo1下的投影zc和zd;
规定实际星上旋转载荷指向在测量的星上旋转载荷本体坐标系op1xp1yp1zp1下的投影为zp1′=[0,0,1]T,则期望星上旋转载荷指向在期望的轨道坐标系oo1xo1yo1zo1下的投影zd为:
实际星上旋转载荷指向在期望的轨道坐标系oo1xo1yo1zo1下的投影zc为:
其中,表示轨道坐标系ooxoyozo相对于期望的轨道坐标系oo1xo1yo1zo1的姿态变换矩阵,表示测量的轨道坐标系oo2xo2yo2zo2相对于轨道坐标系ooxoyozo的姿态变换矩阵,表示测量的遥感卫星平台本体系ob1xb1yb1zb1相对于测量的轨道坐标系oo2xo2yo2zo2的姿态变换矩阵,表示测量的星上旋转载荷本体坐标系op1xp1yp1zp1相对于测量的遥感卫星平台本体系ob1xb1yb1zb1的姿态变换矩阵;
和的具体表达式如下:
结合公式(2)至(6)得到实际星上旋转载荷指向在期望的轨道坐标系下的投影zc,求解两个矢量zd和zc之间的夹角ΔΦ为:
ΔΦ=arccos(zd Tzc) (7)
即得到建立好的星上旋转载荷指向精度模型;其中:arccos(·)表示取反余弦计算;
利用夹角ΔΦ来估计星上旋转载荷的指向精度。
具体实施方式五:本实施方式与具体实施方式四不同的是:所述步骤四的具体过程为:
结合步骤三建立的星上旋转载荷指向精度模型,利用星上旋转载荷相对地心赤道坐标系oxyz的角速度矢量在测量的星上旋转载荷本体坐标系op1xp1yp1zp1下投影的期望值ωd和真实值ωc之差,来衡量星上旋转载荷三轴姿态稳定精度;
在不考虑姿态误差时,测量的遥感卫星平台本体系ob1xb1yb1zb1和期望的轨道坐标系oo1xo1yo1zo1重合,测量的星上旋转载荷本体坐标系op1xp1yp1zp1和星上旋转载荷本体坐标系opxpypzp重合,测量的星上旋转载荷本体坐标系op1xp1yp1zp1相对测量的遥感卫星平台本体系ob1xb1yb1zb1仅存在轴向转动,在这种情形下,测量的遥感卫星平台本体系ob1xb1yb1zb1相对地心赤道坐标系oxyz的角速度和期望的轨道坐标系oo1xo1yo1zo1相对地心赤道坐标系oxyz的角速度一致,且测量的遥感卫星平台本体系ob1xb1yb1zb1相对地心赤道坐标系oxyz的角速度在期望的轨道坐标系oo1xo1yo1zo1下的投影为[0 -ω0 0]T,其中:ω0表示轨道角速度,则星上旋转载荷相对地心赤道坐标系oxyz的角速度矢量在期望的轨道坐标系oo1xo1yo1zo1下的投影ωd1(在没有姿态误差下)表示为:
其中,表示对标量取时间导数;
在考虑姿态误差时,结合方程(2)得到星上旋转载荷相对地心赤道坐标系oxyz的角速度矢量在坐标系op1xp1yp1zp1下投影的期望值ωd为:
其中,是对的转置,是对的转置,是对的转置,是对的转置;同时在考虑姿态误差的真实情况下,测量的遥感卫星平台本体坐标系ob1xb1yb1zb1相对地心赤道坐标系oxyz的角速度在坐标系ob1xb1yb1zb1下的投影ωb为:ωb=[Δωbx -ω0+Δωby Δωbz]T,其中Δωbx、Δωby、Δωbz分别表示遥感卫星平台本体坐标系的角速度误差在ob1xb1yb1zb1系的xb1、yb1、zb1轴下投影。
式中:ωb可从遥感卫星平台上的光纤陀螺测量结果得到;
测量的星上旋转载荷本体坐标系op1xp1yp1zp1相对测量的遥感卫星平台本体坐标系ob1xb1yb1zb1存在姿态误差,需要从坐标系ob1xb1yb1zb1经过zyx顺序依次转过Δψ4、Δθ4、至坐标系op1xp1yp1zp1,记姿态角速率为(其中:是对Δψ4的导数,是对Δθ4的导数,是对的导数);则星上旋转载荷相对地心赤道坐标系oxyz的角速度矢量在坐标系op1xp1yp1zp1下投影的真实值ωc为:
公式(10)中的中间变量Cx、Cy和Cz分别表示为:
结合公式(9)和(10)得到期望值ωd与真实值ωc之差Δω为:
Δω=ωc-ωd (12)
利用期望值ωd与真实值ωc之差Δω估计星上旋转载荷三轴姿态稳定精度。即获得公式(12)的三轴姿态稳定精度模型。
根据式(7)和式(12)可对遥感卫星星上旋转载荷的指向精度和三轴姿态稳定度进行有效评估。
根据理论结果可知,精度模型受误差环节累计作用,考虑到各干扰作用特性,精度模型体现出强动态性;此外,三轴姿态稳定精度和指向精度之间还存在相互耦合特性。如图6所示,利用精度模型预估分析的载荷指向精度和三轴姿态稳定精度,可换算为载荷对地经纬度精度,从而实现针对某个确定的载荷对地经纬度精度需求设计分配各误差环节允许的误差上限,并验证各环节误差设计方案是否满足载荷的对地经纬度需求,还可以用来评价各控制***性能,具有很高的学术价值和工程意义。
本发明的上述算例仅为详细地说明本发明的计算模型和计算流程,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动,这里无法对所有的实施方式予以穷举,凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。
Claims (5)
1.遥感卫星星上旋转载荷动态精度分析及误差分配方法,其特征在于,该方法包括以下步骤:
步骤一、确定遥感卫星***的结构和工作模式,定义地心赤道坐标系oxyz、轨道坐标系ooxoyozo、旋转载荷卫星***质心坐标系osxsyszs、卫星平台本体坐标系obxbybzb和旋转载荷本体坐标系opxpypzp;
步骤二、确定遥感卫星***包含的各误差环节,各误差环节构成精度误差传递链,各误差环节共同作用、累积传递至星上旋转载荷;
步骤三、分别计算实际星上旋转载荷指向和期望星上旋转载荷指向在期望的轨道坐标系oo1xo1yo1zo1下的投影zc和zd,以zd和zc的夹角ΔΦ为指标,来估计星上旋转载荷的指向精度;
步骤四、计算星上旋转载荷相对地心赤道坐标系oxyz的角速度矢量在测量的星上旋转载荷本体坐标系op1xp1yp1zp1下投影的期望值ωd和真实值ωc,利用期望值ωd与真实值ωc之差估计星上旋转载荷三轴姿态稳定精度;
将星上旋转载荷的指向精度和三轴姿态稳定精度的评估结果换算为旋转载荷对地经纬度指标,利用经纬度指标分配及验证遥感卫星***各误差环节的误差上限。
2.根据权利要求1所述的遥感卫星星上旋转载荷动态精度分析及误差分配方法,其特征在于,所述步骤一的具体过程为:
遥感卫星***的结构包括遥感卫星平台子***、星上载荷子***和旋转关节;所述遥感卫星平台子***又包含:遥感卫星平台,平台三轴正交飞轮和平台双侧太阳帆板;所述星上载荷子***又包含:星上旋转载荷和载荷内飞轮;所述旋转关节为磁悬浮轴承;
遥感卫星***的工作模式:遥感卫星平台通过平台三轴正交轮控***保持对地定向,平台双侧太阳帆板单轴驱动保持对日定向,旋转关节配备磁悬浮轴承提供连接,星上旋转载荷通过载荷内飞轮控制保持匀速旋转,载荷内飞轮角动量饱和时需磁力矩器和驱动电机提供卸载方案;
定义坐标系如下:
以地心为坐标原点o,建立地心赤道坐标系oxyz,所述地心赤道坐标系的x轴在J2000地球平赤道面内由地心指向J2000时刻的平春分点,z轴为J2000地球平赤道面的法线且指向北极方向,y轴同x轴、z轴构成右手系;
以遥感卫星***质心为坐标原点oo,建立轨道坐标系ooxoyozo,所述轨道坐标系的xo轴在轨道面内并指向遥感卫星***的前进方向,zo轴由遥感卫星***质心指向地心,yo轴同xo轴、zo轴构成右手系;
记期望的轨道坐标系为oo1xo1yo1zo1,期望的轨道坐标系的坐标原点同轨道坐标系ooxoyozo的坐标原点,在没有轨道控制误差时,期望的轨道坐标系oo1xo1yo1zo1的三轴同轨道坐标系ooxoyozo的三轴重合;记测量的轨道坐标系为oo2xo2yo2zo2,测量的轨道坐标系oo2xo2yo2zo2的坐标原点同轨道坐标系ooxoyozo的坐标原点,在没有轨道确定误差时,测量的轨道坐标系oo2xo2yo2zo2的三轴同轨道坐标系ooxoyozo的三轴重合;
以遥感卫星***质心为坐标原点os,建立遥感卫星***质心坐标系osxsyszs,所述遥感卫星***质心坐标系的xs轴指向遥感卫星平台轴向,ys轴和zs轴位于遥感卫星平台轴向截面内,且ys轴、zs轴同xs轴构成右手系;在无姿态控制误差下,遥感卫星***质心坐标系的三轴同测量的轨道坐标系oo2xo2yo2zo2的三轴重合;
记测量的遥感卫星***质心坐标系为os1xs1ys1zs1,且测量的遥感卫星***质心坐标系os1xs1ys1zs1的坐标原点同遥感卫星***质心坐标系osxsyszs的坐标原点,在无姿态测量误差时,测量的遥感卫星***质心坐标系os1xs1ys1zs1的三轴同遥感卫星***质心坐标系osxsyszs的三轴重合;
以遥感卫星平台质心为坐标原点ob,建立遥感卫星平台本体坐标系obxbybzb,遥感卫星平台本体坐标系的三轴指向同遥感卫星***质心坐标系osxsyszs的三轴指向;
记测量的遥感卫星平台本体坐标系为ob1xb1yb1zb1,且测量的遥感卫星平台本体坐标系ob1xb1yb1zb1的坐标原点同遥感卫星平台本体坐标系obxbybzb的坐标原点,测量的遥感卫星平台本体坐标系ob1xb1yb1zb1的三轴指向同测量的遥感卫星***质心坐标系os1xs1ys1zs1的三轴指向;
以星上旋转载荷质心为坐标原点op,建立星上旋转载荷本体坐标系opxpypzp,在旋转关节连接误差为零、星上旋转载荷自旋转角为零时,星上旋转载荷本体坐标系opxpypzp的三轴指向同测量的遥感卫星***质心坐标系os1xs1ys1zs1的三轴指向;
记测量的星上旋转载荷本体坐标系为op1xp1yp1zp1,且测量的星上旋转载荷本体坐标系op1xp1yp1zp1的坐标原点同遥感卫星平台本体坐标系opxpypzp的坐标原点,在无测量误差时,测量的星上旋转载荷本体坐标系op1xp1yp1zp1的三轴同坐标系opxpypzp的三轴重合。
3.根据权利要求2所述的遥感卫星星上旋转载荷动态精度分析及误差分配方法,其特征在于,所述步骤二的具体过程为:
确定遥感卫星***包含的各误差环节为:轨道误差环节、平台姿态误差环节、磁悬浮轴承连接误差环节和星上旋转载荷自旋误差环节;
轨道误差环节由两部分构成:一部分是轨道的控制误差,造成期望的轨道坐标系oo1xo1yo1zo1和轨道坐标系ooxoyozo之间有相对偏差;另一部分是轨道的确定误差,造成轨道坐标系ooxoyozo和测量的轨道坐标系oo2xo2yo2zo2之间有相对偏差;
平台姿态误差环节由两部分构成:一部分是姿态的控制误差,受姿态干扰影响,造成遥感卫星平台本体系obxbybzb和测量的轨道坐标系oo2xo2yo2zo2之间有相对偏差;另一部分是姿态测量误差,包含星敏感器误差和光纤陀螺误差,造成遥感卫星平台本体系obxbybzb和测量的遥感卫星平台本体系ob1xb1yb1zb1之间有相对偏差;
磁悬浮轴承连接误差环节由两部分构成:一部分是轴隙跳动误差,导致星上旋转载荷本体坐标系opxpypzp相对测量的遥感卫星平台本体系ob1xb1yb1zb1存在径向章动进动姿态偏差;另一部分是测量误差,导致星上旋转载荷本体坐标系opxpypzp和测量的星上旋转载荷本体坐标系op1xp1yp1zp1之间存在径向姿态偏差;
星上旋转载荷自旋误差环节由两部分构成,一部分是转速控制误差,导致星上旋转载荷本体坐标系opxpypzp相对测量的遥感卫星平台本体系ob1xb1yb1zb1存在轴向转速误差;另一部分是测量误差,导致星上旋转载荷本体坐标系opxpypzp和测量的星上旋转载荷本体坐标系op1xp1yp1zp1之间存在轴向转速偏差;
轨道误差环节、平台姿态误差环节、磁悬浮轴承连接误差环节和星上旋转载荷自旋误差环节串联在一起构成精度误差传递链,各误差环节共同作用、累积传递至星上旋转载荷。
4.根据权利要求3所述的遥感卫星星上旋转载荷动态精度分析及误差分配方法,其特征在于,所述步骤三的具体过程为:
轨道坐标系ooxoyozo相对期望的轨道坐标系oo1xo1yo1zo1存在误差,利用zxy顺序描述轨道坐标系ooxoyozo相对期望的轨道坐标系oo1xo1yo1zo1的姿态误差,所转过的欧拉角记为:欧拉角的定义为:期望的轨道坐标系oo1xo1yo1zo1绕zo1轴转过角度至中间系i1,中间系i1绕自身的x轴转过角度Δθ1至中间系i2,中间系i2绕轨道坐标系的yo轴转过Δψ1至轨道坐标系ooxoyozo;
测量的轨道坐标系oo2xo2yo2zo2相对轨道坐标系ooxoyozo存在误差,利用zxy顺序描述测量的轨道坐标系oo2xo2yo2zo2相对轨道坐标系ooxoyozo的姿态误差,所转过的欧拉角记为:欧拉角的定义为:轨道坐标系ooxoyozo绕zo轴转过角度至中间系i1′,中间系i1′绕自身的x轴转过角度Δθ2至中间系i2′,中间系i2′绕测量的轨道坐标系的yo2轴转过角度Δψ2至测量的轨道坐标系oo2xo2yo2zo2;
测量的遥感卫星平台本体系ob1xb1yb1zb1相对测量的轨道坐标系oo2xo2yo2zo2存在误差,利用zxy顺序描述遥感卫星平台本体系ob1xb1yb1zb1相对测量的轨道坐标系oo2xo2yo2zo2的姿态误差,所转过的欧拉角记为:欧拉角的定义为:测量的轨道坐标系oo2xo2yo2zo2绕zo2轴转过角度至中间系i1″,中间系i1″绕自身的x轴转过角度Δθ3至中间系i2″,中间系i2″绕遥感卫星平台本体系的yb1轴转过角度Δψ3至遥感卫星平台本体系ob1xb1yb1zb1;
测量的星上旋转载荷本体坐标系op1xp1yp1zp1相对测量的遥感卫星平台本体系ob1xb1yb1zb1存在误差,利用zyx描述测量的星上旋转载荷本体坐标系op1xp1yp1zp1相对测量的遥感卫星平台本体系ob1xb1yb1zb1的姿态误差,所转过的欧拉角记为:欧拉角的定义为:测量的遥感卫星平台本体系ob1xb1yb1zb1绕zb1轴转过角度Δψ4至中间系i1″′,中间系i1″′绕自身的y轴转过角度Δθ4至中间系i2″′,中间系i2″′绕测量的星上旋转载荷本体坐标系的xp1轴转过角度至测量的星上旋转载荷本体坐标系op1xp1yp1zp1;代表星上旋转载荷相对遥感卫星平台的期望转角,代表星上旋转载荷相对遥感卫星平台的转角偏差;
分别计算实际星上旋转载荷指向和期望星上旋转载荷指向在期望的轨道坐标系oo1xo1yo1zo1下的投影zc和zd;
规定实际星上旋转载荷指向在测量的星上旋转载荷本体坐标系op1xp1yp1zp1下的投影为zp1′=[0,0,1]T,则期望星上旋转载荷指向在期望的轨道坐标系oo1xo1yo1zo1下的投影zd为:
实际星上旋转载荷指向在期望的轨道坐标系oo1xo1yo1zo1下的投影zc为:
其中,表示轨道坐标系ooxoyozo相对于期望的轨道坐标系oo1xo1yo1zo1的姿态变换矩阵,表示测量的轨道坐标系oo2xo2yo2zo2相对于轨道坐标系ooxoyozo的姿态变换矩阵,表示测量的遥感卫星平台本体系ob1xb1yb1zb1相对于测量的轨道坐标系oo2xo2yo2zo2的姿态变换矩阵,表示测量的星上旋转载荷本体坐标系op1xp1yp1zp1相对于测量的遥感卫星平台本体系ob1xb1yb1zb1的姿态变换矩阵;
和的具体表达式如下:
结合公式(2)至(6)得到实际星上旋转载荷指向在期望的轨道坐标系下的投影zc,求解两个矢量zd和zc之间的夹角ΔΦ为:
ΔΦ=arccos(zd Tzc) (7)
其中:arccos(·)表示取反余弦计算;
利用夹角ΔΦ来估计星上旋转载荷的指向精度。
5.根据权利要求4所述的遥感卫星星上旋转载荷动态精度分析及误差分配方法,其特征在于,所述步骤四的具体过程为:
在不考虑姿态误差时,测量的遥感卫星平台本体系ob1xb1yb1zb1和期望的轨道坐标系oo1xo1yo1zo1重合,测量的星上旋转载荷本体坐标系op1xp1yp1zp1和星上旋转载荷本体坐标系opxpypzp重合,测量的星上旋转载荷本体坐标系op1xp1yp1zp1相对测量的遥感卫星平台本体系ob1xb1yb1zb1仅存在轴向转动,测量的遥感卫星平台本体系ob1xb1yb1zb1相对地心赤道坐标系oxyz的角速度和期望的轨道坐标系oo1xo1yo1zo1相对地心赤道坐标系oxyz的角速度一致,且测量的遥感卫星平台本体系ob1xb1yb1zb1相对地心赤道坐标系oxyz的角速度在期望的轨道坐标系oo1xo1yo1zo1下的投影为[0 -ω0 0]T,其中:ω0表示轨道角速度,则星上旋转载荷相对地心赤道坐标系oxyz的角速度矢量在期望的轨道坐标系oo1xo1yo1zo1下的投影ωd1表示为:
其中,表示对标量取时间导数;
在考虑姿态误差时,结合方程(2)得到星上旋转载荷相对地心赤道坐标系oxyz的角速度矢量在坐标系op1xp1yp1zp1下投影的期望值ωd为:
其中,是对的转置,是对的转置,是对的转置,是对的转置;测量的遥感卫星平台本体坐标系ob1xb1yb1zb1相对地心赤道坐标系oxyz的角速度在坐标系ob1xb1yb1zb1下的投影ωb为:ωb=[Δωbx -ω0+Δωby Δωbz]T,其中Δωbx、Δωby、Δωbz分别表示遥感卫星平台本体坐标系的角速度误差在ob1xb1yb1zb1系的xb1、yb1、zb1轴下投影。
测量的星上旋转载荷本体坐标系op1xp1yp1zp1相对测量的遥感卫星平台本体坐标系ob1xb1yb1zb1存在姿态误差,需要从坐标系ob1xb1yb1zb1经过zyx顺序依次转过Δψ4、Δθ4、至坐标系op1xp1yp1zp1,记姿态角速率为则星上旋转载荷相对地心赤道坐标系oxyz的角速度矢量在坐标系op1xp1yp1zp1下投影的真实值ωc为:
公式(10)中的中间变量Cx、Cy和Cz分别表示为:
结合公式(9)和(10)得到期望值ωd与真实值ωc之差Δω为:
Δω=ωc-ωd (12)
利用期望值ωd与真实值ωc之差Δω估计星上旋转载荷三轴姿态稳定精度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910391484.XA CN110162855B (zh) | 2019-05-10 | 2019-05-10 | 遥感卫星星上旋转载荷动态精度分析及误差分配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910391484.XA CN110162855B (zh) | 2019-05-10 | 2019-05-10 | 遥感卫星星上旋转载荷动态精度分析及误差分配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110162855A true CN110162855A (zh) | 2019-08-23 |
CN110162855B CN110162855B (zh) | 2022-06-28 |
Family
ID=67634515
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910391484.XA Active CN110162855B (zh) | 2019-05-10 | 2019-05-10 | 遥感卫星星上旋转载荷动态精度分析及误差分配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110162855B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111062976A (zh) * | 2019-12-25 | 2020-04-24 | 中国科学院长春光学精密机械与物理研究所 | 基于fmt的低轨卫星太阳望远镜遥感图像配准方法 |
CN111381256A (zh) * | 2020-03-10 | 2020-07-07 | 上海卫星工程研究所 | 主动遥感卫星天线相位中心偏移误差计算的方法和*** |
CN111536991A (zh) * | 2020-04-26 | 2020-08-14 | 哈尔滨工业大学 | 遥感卫星多模式成像载荷模拟装置 |
CN112607065A (zh) * | 2020-12-23 | 2021-04-06 | 长春工业大学 | 一种基于电推进***的高精度相位控制方法 |
CN113541774A (zh) * | 2021-09-17 | 2021-10-22 | 椭圆时空(北京)科技有限公司 | 一种通遥一体化卫星*** |
CN115586305A (zh) * | 2022-10-09 | 2023-01-10 | 哈尔滨工业大学 | 一种基于载荷垂轨旋转扫描的超宽覆盖多模式碳监测卫星、构建方法和监测方法 |
CN116659543A (zh) * | 2023-06-21 | 2023-08-29 | 中国人民解放军61540部队 | 基于遥感卫星轨道根数的卫星位置姿态估计方法和装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6273237B1 (en) * | 1999-11-10 | 2001-08-14 | Pri Automation, Inc. | Rotating shuttle payload platform |
CN108762283A (zh) * | 2018-05-08 | 2018-11-06 | 哈尔滨工业大学 | 带有旋转载荷的卫星平台的姿态控制方法、装置及*** |
CN108959796A (zh) * | 2018-07-18 | 2018-12-07 | 哈尔滨工业大学 | 一种大惯量旋转载荷卫星的刚柔磁耦合动力学建模方法 |
CN109033604A (zh) * | 2018-07-18 | 2018-12-18 | 哈尔滨工业大学 | 含旋转载荷的卫星动力学建模及轴承处受力的确定方法 |
-
2019
- 2019-05-10 CN CN201910391484.XA patent/CN110162855B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6273237B1 (en) * | 1999-11-10 | 2001-08-14 | Pri Automation, Inc. | Rotating shuttle payload platform |
CN108762283A (zh) * | 2018-05-08 | 2018-11-06 | 哈尔滨工业大学 | 带有旋转载荷的卫星平台的姿态控制方法、装置及*** |
CN108959796A (zh) * | 2018-07-18 | 2018-12-07 | 哈尔滨工业大学 | 一种大惯量旋转载荷卫星的刚柔磁耦合动力学建模方法 |
CN109033604A (zh) * | 2018-07-18 | 2018-12-18 | 哈尔滨工业大学 | 含旋转载荷的卫星动力学建模及轴承处受力的确定方法 |
Non-Patent Citations (2)
Title |
---|
XIBINCAO: "Flexible platform based micro-satellite design method", 《AEROSPACE SCIENCE AND TECHNOLOGY》 * |
李明群 等: "具有载荷扰动抑制的卫星姿态控制方法及验证", 《中国空间科学技术》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111062976A (zh) * | 2019-12-25 | 2020-04-24 | 中国科学院长春光学精密机械与物理研究所 | 基于fmt的低轨卫星太阳望远镜遥感图像配准方法 |
CN111062976B (zh) * | 2019-12-25 | 2023-02-28 | 中国科学院长春光学精密机械与物理研究所 | 基于fmt的低轨卫星太阳望远镜遥感图像配准方法 |
CN111381256A (zh) * | 2020-03-10 | 2020-07-07 | 上海卫星工程研究所 | 主动遥感卫星天线相位中心偏移误差计算的方法和*** |
CN111536991A (zh) * | 2020-04-26 | 2020-08-14 | 哈尔滨工业大学 | 遥感卫星多模式成像载荷模拟装置 |
CN112607065A (zh) * | 2020-12-23 | 2021-04-06 | 长春工业大学 | 一种基于电推进***的高精度相位控制方法 |
CN113541774A (zh) * | 2021-09-17 | 2021-10-22 | 椭圆时空(北京)科技有限公司 | 一种通遥一体化卫星*** |
CN113541774B (zh) * | 2021-09-17 | 2021-12-24 | 椭圆时空(北京)科技有限公司 | 一种通遥一体化卫星*** |
CN115586305A (zh) * | 2022-10-09 | 2023-01-10 | 哈尔滨工业大学 | 一种基于载荷垂轨旋转扫描的超宽覆盖多模式碳监测卫星、构建方法和监测方法 |
CN116659543A (zh) * | 2023-06-21 | 2023-08-29 | 中国人民解放军61540部队 | 基于遥感卫星轨道根数的卫星位置姿态估计方法和装置 |
CN116659543B (zh) * | 2023-06-21 | 2024-05-07 | 中国人民解放军61540部队 | 基于遥感卫星轨道根数的卫星位置姿态估计方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN110162855B (zh) | 2022-06-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110162855A (zh) | 遥感卫星星上旋转载荷动态精度分析及误差分配方法 | |
CN104527994B (zh) | 异面交叉快变轨道固定时间稳定姿态指向跟踪控制方法 | |
CN109573105B (zh) | 末子级留轨应用子***姿态控制方法 | |
CN111099045B (zh) | 双超卫星动力学与控制气浮平台全物理仿真方法 | |
Chen et al. | Spacecraft angular rates and angular acceleration estimation using single-gimbal magnetically suspended control moment gyros | |
Takada et al. | Control moment gyro singularity-avoidance steering control based on singular-surface cost function | |
CN108871326B (zh) | 一种单轴旋转调制惯性-天文深组合导航方法 | |
CN109823571A (zh) | 一种遥感微纳卫星的多阶段姿态控制方法 | |
Shen et al. | Control of single-axis rotation INS by tracking differentiator based fuzzy PID | |
CN108181916B (zh) | 小卫星相对姿态的控制方法及装置 | |
CN106742071B (zh) | 一种利用正交安装飞轮控制避免过零的方法 | |
CN1076419A (zh) | 卫星姿态控制的测试设备及其应用 | |
JP2011042358A (ja) | パドル電流測定値フィードバックのみを使用したジャイロなし遷移軌道太陽捕捉 | |
CN102001453B (zh) | 一种增加前馈补偿改善磁控能力的方法 | |
CN103019247A (zh) | 一种火星探测器无陀螺自主空间姿态机动控制方法 | |
CN109033604B (zh) | 含旋转载荷的卫星动力学建模及轴承处受力的确定方法 | |
CN106052682A (zh) | 一种混合式惯性导航***及导航方法 | |
CN110304279A (zh) | 一种电推进卫星的质心在轨标定补偿方法 | |
CN108502209B (zh) | 一种基于陀螺积分计算太阳矢量的卫星自旋太阳捕获方法 | |
CN110803304A (zh) | 一种卫星姿态控制*** | |
CN101571720A (zh) | 一种磁悬浮飞轮三自由度动量交换控制方法 | |
CN113335567B (zh) | 一种微小卫星的轮磁混合姿态控制方法和*** | |
CN114802818A (zh) | 晨昏轨道卫星及其对日姿态计算方法、导引方法 | |
Zhu et al. | Attitude measurement method of geomagnetic/gyro combination based on difference-quotient | |
Carabellese et al. | Magnetorquer-only nonlinear attitude control for cubesats |
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 |