CN103954286B - 微小卫星磁传感器多误差模型在轨迭代标定方法 - Google Patents
微小卫星磁传感器多误差模型在轨迭代标定方法 Download PDFInfo
- Publication number
- CN103954286B CN103954286B CN201410169216.0A CN201410169216A CN103954286B CN 103954286 B CN103954286 B CN 103954286B CN 201410169216 A CN201410169216 A CN 201410169216A CN 103954286 B CN103954286 B CN 103954286B
- Authority
- CN
- China
- Prior art keywords
- theta
- satellite
- error
- gaussmeter
- axis
- 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
- 238000000034 method Methods 0.000 title claims abstract description 18
- 108091092878 Microsatellite Proteins 0.000 title claims description 10
- 238000005286 illumination Methods 0.000 claims abstract description 10
- 238000005259 measurement Methods 0.000 claims description 27
- 239000011159 matrix material Substances 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 239000002131 composite material Substances 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 2
- 230000005389 magnetism Effects 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measuring Magnetic Variables (AREA)
- Navigation (AREA)
Abstract
本发明公开了一种基于多参数的二步微小卫星磁传感器多误差模型在轨迭代标定方法。首先在光照区利用太阳敏感器和陀螺计算得到卫星俯仰角,然后以地面标定或上一次测量的标度因数为标准,只标定零偏误差,最后重新建立多误差标定模型,采用最佳拟合的方法重新标定出零偏和标度因数误差。本发明通过标定磁强计的零偏和标度因数误差,并在姿态解算中予以补偿,从而有效地提高卫星导航精度;其不需要专门针对标定算法设计新的滤波器,也不需要实时计算,有效地降低了卫星运行成本和负载。
Description
技术领域
本发明涉及一种微小卫星磁传感器多误差模型在轨迭代标定方法,尤其是一种基于多参数的二步微小卫星磁传感器多误差模型在轨迭代标定方法。
背景技术
磁强计是微小卫星姿态测量的主要部件,磁强计的测量精度是影响微小卫星定姿性能的重要指标。
卫星的剩磁干扰是星上磁强计的主要误差来源,控制去除剩磁干扰的方法主要包括:(1)、磁补法。即将发射前的剩磁测量与补偿的方法相结合对剩磁量加以控制,但该方法作用有限,经过磁补以后的剩磁量一般还有900nT左右,对姿态测量将产生较大的误差;(2)、标定法。就是确定磁强计的误差参数,以便在姿态解算中用于计算补偿。传统的在轨标定算法中多数只标定零偏误差,但标度因数误差也应考虑,还有学者提出用多阶滤波器进行全误差模型标定,需要实时计算,增加了***的计算负担。
2011年12月,中国惯性技术学报,第19卷第6期,614-617页,作者:华冰,郁丰,公开了一种基于光照区太阳敏感器和陀螺辅助修正的微小卫星磁测技术,其采用太阳敏感器和陀螺对磁强计误差进行辅助测量与修正,推导了磁强计误差估计方法,在光照区以太阳敏感器与陀螺输出作为俯仰滤波器观测量,估计出卫星俯仰角度和角速度,再采用最小二乘方法,利用滤波输出量对磁强计误差进行估计,估计的结果进入滤波器对磁场输出进行测量修正。该方法,姿态角精度提高了1°左右,角速度精度最高提高了0.003(°)/s左右,并增强了卫星稳定性。但是该方法只能标定零偏误差,无法标定标度因数误差,导致只能标定简化的磁强计模型,实用精度并不高。
发明内容
本发明所要解决的技术的问题在于克服现有技术的缺陷,提供一种通过同时标定零偏误差和标度因数误差来提高导航定位精度的微小卫星磁传感器多误差模型在轨迭代标定方法。
为了解决上述技术问题,本发明提供的微小卫星磁传感器多误差模型在轨迭代标定方法,包括以下步骤:
1)、在光照区利用太阳敏感器和陀螺计算得到卫星俯仰角θ;
2)、以地面标定或者上一次测量的标度因数为标准,建立零偏误差粗标定模型:
HiE=Gi (14)
式(14)中,Hi为量测量和已知量计算得到的用于标定的矢量,Hi=[kz(bxosinθ+bzocosθ)kx(bzosinθ-bxocosθ)];E为磁强计零偏误差向量,Gi为量测量和已知量计算得到的用于标定的标量, 其中,θ为卫星俯仰角,为X轴磁强计测量值,为X轴磁强计零偏误差,kx为X轴磁强计标度因数,为Z轴磁强计测量值,为Z轴磁强计零偏误差,kz为Z轴磁强计标度因数,bxo为卫星轨道坐标系X轴磁场强度,bzo为卫星轨道坐标系Z轴磁场强度,kx、kz均为常数;
3)、采集N组测量数据组成测量模型:
4)、采用最小二乘估计法估计得出磁强计零偏误差向量E:
E=(HTH)-1HTG (17)
式(17)中,H为量测矩阵,G为量测矢量,HT为H的转置;
5)、定义函数f:
6)、以步骤2)中的标度因数和步骤4)得到的零偏误差为起始点,kx,kz,小量变化,采用N组测量数据对函数f进行计算,将每组等式记为fi,即得到:
F=f1 2+f2 2+…+fN 2 (19)
式(19)中,F为N组量测数据构成的代价函数;
7)、以F作为标定的目标函数,并将优化算法设置为最小值搜索,在设定标度因数和零偏参数范围内,令F为最小值的误差参数需要求解的量。
所述步骤1)中卫星俯仰角θ采用k时刻俯仰角θ(k),其具体步骤为:
11)、利用太阳敏感器计算卫星的俯仰角:
式(10)中,θs为利用太阳敏感器输出计算得到的卫星俯仰角,rx为测量的卫星本体系X轴向太阳分量,rz为测量的卫星本体系Z轴向太阳分量,rxo为轨道参考系下的X轴向太阳分量,rzo为轨道参考系下的Z轴向太阳分量;
12)、得到Y方向陀螺输出与卫星的俯仰角速度的关系:
式(11)中,为陀螺测得的卫星俯仰角速度,ωy为陀螺输出,ε为陀螺输出误差;
13)、综合式(10)与(11),得出卫星的k时刻俯仰角θ(k):
式(12)中,θs(k)为k时刻太阳敏感器计算出的俯仰角信息,为k时刻俯仰角速度信息,Δ为k-1到k时刻采样步长。
本发明的有益效果在于:(1)、本发明通过同时计算磁强计的零偏和标度因数误差,并在姿态解算中予以补偿,从而有效地提高卫星导航精度;(2)、本发明不需要专门针对标定算法设计新的滤波器,也不需要实时计算,有效地降低了卫星运行成本和负载。
附图说明
图1为本发明标定流程示意图;
图2为俯仰角速率曲线图;
图3为俯仰角曲线图;
图4为俯仰角放大曲线图。
具体实施方式
下面结合附图对本发明作进一步详细说明。
如图1所示,本发明首先在光照区利用太阳敏感器和陀螺计算得到卫星的俯仰角,然后以地面标定或上一次测量的标度因数误差为标准,只标定零偏误差,最后重新建立多误差标定模型,采用最佳拟合的方法重新标定出零偏和标度因数误差。在光照区以太阳敏感器与陀螺输出作为俯仰滤波器观测量,估计出卫星俯仰角度和角速度,利用滤波输出量对磁强计误差进行估计,估计的结果进入滤波器对磁场输出进行测量修正。
一、卫星姿态估计模型分析
1、磁强计误差模型
卫星磁强计的误差主要考虑为零偏和标度因数误差。剩磁是星上磁强计的主要误差来源,剩磁造成的磁强计误差虽然是时变的,但是短期内可以近似认为是常值,反映在误差项上即为星上磁强计零偏误差,同时考虑磁强计标度因数,则误差模型如下:
式(1)中,b为磁强计真值,k为磁强计标度因数,为磁强计测量值,bε为磁强计零偏误差,σ为磁强计随机量测噪声;
2、采用磁强计输出信息为量测量的俯仰滤波器
卫星的测量信息估计方法采用俯仰滤波器。卫星在俯仰Y轴方向上的姿态动力学与X、Z轴解耦,则有:
式(2)中,θ是卫星相对于轨道的Y向姿态角,Y向一般被称为俯仰方向,因此该姿态角简称卫星俯仰角,Ny为作用在Y轴上的控制力矩,Iyy为Y轴方向的转动惯量;
状态方程:
测量方程:
式(3)、(4)中,Δt为星载计算机滤波步长,k表示时刻,bx为X轴磁强计真值,bz为Z轴磁强计真值,bxo为卫星轨道坐标系X轴磁场强度,bzo为卫星轨道坐标系Z轴磁场强度。
二、基于多参数的磁强计迭代误差估计
1、磁强计多参数误差标定分析
由式(4)可得:
bx(bxosinθ+bzocosθ)+bz(bzosinθ-bxocosθ)=0 (5)
考虑磁强计测量误差主要为零偏误差和标度因数误差,:
式(5)、(6)中,bx为X轴磁强计真值,为X轴磁强计测量值,由磁强计测量输出获得,为X轴磁强计零偏误差,kx为X轴磁强计标度因数,由地面标定或者上一次测量的标度因数误差暂定为常数,bz为Z轴磁强计真值,为Z轴磁强计测量值,由磁强计测量输出获得,为Z轴磁强计零偏误差,kz为Z轴磁强计标度因数,由地面标定或者上一次测量的标度因数误差暂定为常数,bxo为卫星轨道坐标系X轴磁场强度,由卫星预存的磁场表查表获得,bzo为卫星轨道坐标系Z轴磁场强度,由卫星预存的磁场表查表获得,所以有:
对式(7)进一步推导:
写成矩阵形式:
该方程为多元齐次方程,存在无穷多个解,因此不能直接用最小二乘的算法来采用多组测量值求解。
2、磁强计多参数人迭代误差标定
由于磁强计的输出耦合了误差,俯仰滤波器中的姿态估算存在较大误差。在光照区间太阳敏感器可用,陀螺也可短期开机,因此可以选择精度较高的太阳敏感器和陀螺估计出姿态信息。
利用太阳敏感器信息计算卫星的俯仰角为:
式(10)中,θs为利用太阳敏感器输出计算得到的卫星俯仰角,rx为测量的卫星本体系X轴向太阳分量,rz为测量的卫星本体系Z轴向太阳分量,rxo为轨道参考系下的X轴向太阳分量,rzo为轨道参考系下的Z轴向太阳分量;
Y方向陀螺输出与卫星的俯仰角速度的关系为:
式(11)中,为陀螺测得的卫星俯仰角速度,ωy为陀螺输出,ε为陀螺输出误差(含轨道角速率)。
综合式(10)与(11),忽略误差项的影响,得到k时刻卫星俯仰角估计值为:
式(12)中,θ(k)为k时刻俯仰角估计值,θs(k)为k时刻太阳敏感器计算出的俯仰角信息,为k时刻俯仰角速度信息,Δ为k-1到k时刻采样步长。
分析磁强计的零偏和标定因数的特点,标度因数是磁强计测量与真实输出的比例系数,变化量级相对较小,而零偏受星上剩磁影响严重,变化量级较大。当考虑标度因数kx,kz为常数时,将方程(7)改写成以下零偏误差粗标定模型:
即:
HiE=Gi (14)
式(14)中,Hi为量测量和已知量计算得到的用于标定的矢量,Hi=[kz(bxosinθ+bzocosθ)kx(bzosinθ-bxocosθ)];E为磁强计零偏误差向量,Gi为量测量和已知量计算得到的用于标定的标量,
采集N组测量数据,测量值组成测量模型:
则:
HE=G (16)
式(15)、(16)中,H为量测矩阵,G为量测矢量;
利用最小二乘估计法,可以估计出磁强计零偏误差向量E:
E=(HTH)-1HTG (17)
式(17)中,HT为H的转置;
将式(7)定义为函数f:
以地面标定或上一次测量的标度因数和标定后的零偏误差为起始点,kx,kz,小量变化,采用多组测量数据对函数f进行计算;在多组测量数据下,可得到不同的函数计算值,将每组等式记为fi。
因此,记
F=f1 2+f2 2+…+fN 2 (19)
式(19)中,F为N组量测数据构成的代价函数;
将F作为标定参数问题的目标函数,并将优化算法设置为最小值搜索。则在设定标度因数和零偏参数范围内,其中标度因数范围可以取[标度因数起始点-0.1,标度因数起始点+0.1],零偏范围可以取[零偏起始点-0.1,零偏起始点+0.1],令F为最小值的误差参数为需要求解的量,即kx,kz,
在光照区,利用式(17),其中俯仰角计算采用式(12)估计值,轨道系磁场强度通过轨道信息查表求得,则可以估计出磁强计误差。由于此时估计的角度和角速度信息与剩磁误差不相关,因此计算出的俯仰角和俯仰角速度不仅可以用于姿态控制,也可以作为辅助信息进入滤波器对磁场输出进行测量修正。
3、仿真验证
卫星轨道为10∶30AM太阳同步轨道,高度496km。磁强计量测噪声均方差30nT,输出频率10Hz。设置卫星的剩余磁矩为[0.30.30.3]A·m2,为更进一步模拟磁场误差,假设磁强计标度因数为0.988,零偏误差为800nT,地面标度因数试验结果为0.991。卫星转动惯量为diag[0.80.80.8]kg·m2,姿态控制***使卫星保持对地稳定,并加入了剩磁扰动力矩,其余干扰力矩假设为白噪声,均方差为5×10-7N·m。太阳敏感器视场范围±30°×±30°,测量精度优于0.5°(3σ),陀螺的测量精度约为0.003(°)/s。
每次进入光照区开始数据采集,同时进行实时估计,估计出的误差在滤波器中用于量测修正。假设标度因数为0.991的常量,经过两个轨道周期后的数据采集,X轴磁强计零偏误差结果为988nT,Z轴磁强计零偏误差估计结果为637T。然后开始误差的迭代标定,标度因数以0.001为步长,迭代范围为[0.991-0.1,0.991+0.1],零偏以30nT为步长,迭代范围为分别为[988-600,988+600],[637-600,637+600]。重新经过两个轨道周期的数据采集计算后,kx,kz,标定结果分别为0.986、0.990、904、723。
如图2、3所示,将只标定磁场零偏与同时标定标度因数和零偏两种情况下的参数用于滤波补偿,补偿后的俯仰角和俯仰角速度曲线,可见在两种补偿方式下***都是稳定的。
如图4所示,将俯仰角曲线局部放大,b曲线为只标定磁场零偏的补偿结果,a曲线为同时标定标度因数和零偏的补偿结果,可以看到a曲线相比b曲线有明显改善,整体波动更为平缓。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下还可以做出若干改进,这些改进也应视为本发明的保护范围。
Claims (2)
1.一种微小卫星磁传感器多误差模型在轨迭代标定方法,其特征在于包括以下步骤:
1)、在光照区利用太阳敏感器和陀螺计算得到卫星俯仰角θ;
2)、以地面标定或者上一次测量的标度因数为标准,建立零偏误差粗标定模型:
HiE=Gi (14)
式(14)中,Hi为量测量和已知量计算得到的用于标定的矢量,Hi=[kz(bxo sinθ+bzo cosθ)kx(bzo sinθ-bxo cosθ)];E为磁强计零偏误差向量,Gi为量测量和已知量计算得到的用于标定的标量,其中,θ为卫星俯仰角,为X轴磁强计测量值,为X轴磁强计零偏误差,kx为X轴磁强计标度因数,为Z轴磁强计测量值,为Z轴磁强计零偏误差,kz为Z轴磁强计标度因数,bxo为卫星轨道坐标系X轴磁场强度,bzo为卫星轨道坐标系Z轴磁场强度,kx、kz均为常数;
3)、采集N组测量数据组成测量模型:
4)、采用最小二乘估计法估计得出磁强计零偏误差向量E:
E=(HTH)-1HTG (17)
式(17)中,H为量测矩阵,G为量测矢量,HT为H的转置;
5)、定义函数f:
6)、以步骤2)中的标度因数和步骤4)得到的零偏误差为起始点,kx,kz,小量变化,采用N组测量数据对函数f进行计算,将每组等式记为fi,即得到:
F=f1 2+f2 2+…+fN 2 (19)
式(19)中,F为N组量测数据构成的代价函数;
7)、以F作为标定的目标函数,并将优化算法设置为最小值搜索,在设定标度因数和零偏参数范围内,令F为最小值的误差参数为需要求解的量,即kx,kz,
2.根据权利要求1所述的微小卫星磁传感器多误差模型在轨迭代标定方法,其特征在于所述步骤1)中卫星俯仰角θ采用k时刻俯仰角θ(k),其步骤为:
11)、利用太阳敏感器计算卫星的俯仰角:
式(10)中,θs为利用太阳敏感器输出计算得到的卫星俯仰角,rx为测量的卫星本体系X轴向太阳分量,rz为测量的卫星本体系Z轴向太阳分量,rxo为轨道参考系下的X轴向太阳分量,rzo为轨道参考系下的Z轴向太阳分量;
12)、得到Y方向陀螺输出与卫星的俯仰角速度的关系:
式(11)中,为陀螺测得的卫星俯仰角速度,ωy为陀螺输出,ε为陀螺输出误差;
13)、综合式(10)与(11),得出卫星的k时刻俯仰角θ(k):
式(12)中,θs(k)为k时刻太阳敏感器计算出的俯仰角信息,为k时刻俯仰角速度信息,Δ为k-1到k时刻采样步长。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410169216.0A CN103954286B (zh) | 2014-04-24 | 2014-04-24 | 微小卫星磁传感器多误差模型在轨迭代标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410169216.0A CN103954286B (zh) | 2014-04-24 | 2014-04-24 | 微小卫星磁传感器多误差模型在轨迭代标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103954286A CN103954286A (zh) | 2014-07-30 |
CN103954286B true CN103954286B (zh) | 2016-08-24 |
Family
ID=51331596
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410169216.0A Active CN103954286B (zh) | 2014-04-24 | 2014-04-24 | 微小卫星磁传感器多误差模型在轨迭代标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103954286B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106033131B (zh) * | 2015-03-20 | 2018-09-11 | 阿里巴巴集团控股有限公司 | 一种地磁传感器校准方法、装置及智能设备 |
CN105937911A (zh) * | 2016-07-01 | 2016-09-14 | 南京理工大学 | 一种磁传感器姿态解算方法 |
CN106403932B (zh) * | 2016-08-24 | 2019-07-23 | 易文俊 | 一种弹载地磁姿态测量处理算法的验证方法 |
CN107270940B (zh) * | 2017-07-11 | 2019-11-15 | 深圳清华大学研究院 | 三轴磁强计的输出偏差估计方法 |
CN108631859B (zh) * | 2018-07-20 | 2020-11-20 | 中国空间技术研究院 | 一种通信卫星星表设备场强预估方法及*** |
CN109932672B (zh) * | 2019-03-27 | 2021-03-26 | 上海微小卫星工程中心 | 一种修正三轴磁强计的误差的方法 |
CN111256744B (zh) * | 2020-02-27 | 2021-06-29 | 苏州海之博电子科技有限公司 | 一种线性输出位置传感器标定方法 |
CN111272168B (zh) * | 2020-03-24 | 2022-03-29 | 哈尔滨工业大学 | 一种基于磁场特征矢量的定位方法、装置及*** |
CN111289933B (zh) * | 2020-05-07 | 2020-09-01 | 中航金城无人***有限公司 | 一种多旋翼飞行器的磁力计自动校准***和方法 |
CN111504301B (zh) * | 2020-05-08 | 2022-02-18 | 哈尔滨工业大学 | 一种基于磁场特征矢量的定位方法、装置及*** |
CN112034402B (zh) * | 2020-07-16 | 2021-12-28 | 中国人民解放军军事科学院国防科技创新研究院 | 微纳卫星的剩磁和剩磁矩联合标定方法 |
CN111891394B (zh) * | 2020-08-11 | 2022-01-04 | 北京控制工程研究所 | 一种卫星冷气推进***流量传感器在轨标定方法 |
CN111982128B (zh) * | 2020-08-21 | 2022-03-18 | 上海航天控制技术研究所 | 一种陀螺安装矩阵的在轨标定方法及星载软件 |
CN113109751B (zh) * | 2021-04-15 | 2021-11-30 | 中国科学院地质与地球物理研究所 | 一种矢量磁强计在轨实时标校***及方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1096581A (zh) * | 1993-06-16 | 1994-12-21 | 张学孚 | 分量磁强计校正方法和可校正装置 |
CN102607564A (zh) * | 2012-03-09 | 2012-07-25 | 北京航空航天大学 | 一种基于星光/地磁组合信息的小卫星自主导航***及其导航方法 |
CN102735268A (zh) * | 2012-07-10 | 2012-10-17 | 中国人民解放军国防科学技术大学 | 一种基于姿态优化激励的捷联三轴磁力仪标定方法 |
CN103278165A (zh) * | 2013-05-22 | 2013-09-04 | 上海新跃仪表厂 | 基于剩磁标定的磁测及星光备份的自主导航方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080243417A1 (en) * | 2007-03-26 | 2008-10-02 | Destinator Technologies, Inc. | Magnetometer normalization |
US7890262B2 (en) * | 2008-03-31 | 2011-02-15 | Honeywell International Inc. | Position estimation for navigation devices |
-
2014
- 2014-04-24 CN CN201410169216.0A patent/CN103954286B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1096581A (zh) * | 1993-06-16 | 1994-12-21 | 张学孚 | 分量磁强计校正方法和可校正装置 |
CN102607564A (zh) * | 2012-03-09 | 2012-07-25 | 北京航空航天大学 | 一种基于星光/地磁组合信息的小卫星自主导航***及其导航方法 |
CN102735268A (zh) * | 2012-07-10 | 2012-10-17 | 中国人民解放军国防科学技术大学 | 一种基于姿态优化激励的捷联三轴磁力仪标定方法 |
CN103278165A (zh) * | 2013-05-22 | 2013-09-04 | 上海新跃仪表厂 | 基于剩磁标定的磁测及星光备份的自主导航方法 |
Non-Patent Citations (4)
Title |
---|
关于剩磁对卫星姿态确定与控制影响的研究;黄琳等;《航空学报》;20060531;第27卷(第3期);全文 * |
基于光照区太阳敏感器和陀螺辅助修正的微小卫星磁测技术;华冰等;《中国惯性技术学报》;20111231;第19卷(第6期);全文 * |
微小卫星剩磁在轨标定技术研究;郁丰等;《航空学报》;20070930;第28卷(第5期);全文 * |
空间环境探测卫星用磁强计误差分析及在线标定;杨照华等;《宇航学报》;20120831;第33卷(第8期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN103954286A (zh) | 2014-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103954286B (zh) | 微小卫星磁传感器多误差模型在轨迭代标定方法 | |
CN103630137B (zh) | 一种用于导航***的姿态及航向角的校正方法 | |
CN102486377B (zh) | 一种光纤陀螺捷联惯导***初始航向的姿态获取方法 | |
CN105509737B (zh) | 一种不受地磁变化影响的航空运动平台磁干扰补偿方法 | |
CN104698485B (zh) | 基于bd、gps及mems的组合导航***及导航方法 | |
CN104165638B (zh) | 一种双轴旋转惯导***多位置自主标定方法 | |
CN103954303B (zh) | 一种用于磁力计导航***航向角动态计算及校正方法 | |
CN102252689A (zh) | 一种基于磁传感器的电子罗盘校准方法 | |
CN112833917B (zh) | 基于磁航向角和最小二乘法的三轴磁传感器标定方法 | |
CN103852085B (zh) | 一种基于最小二乘拟合的光纤捷联惯导***现场标定方法 | |
CN105510849A (zh) | 航磁干扰补偿方法 | |
CN107063254A (zh) | 一种陀螺地磁组合的姿态解算方法 | |
CN107024673A (zh) | 基于陀螺仪辅助的三轴磁强计全误差标定方法 | |
CN110274586A (zh) | 包含多光系原子磁力仪方向误差补偿的航空磁补偿方法 | |
CN103217699A (zh) | 一种基于偏振信息的组合导航***递推优化初始对准方法 | |
CN107270940A (zh) | 三轴磁强计的输出偏差估计方法 | |
Feng et al. | A new method of initial alignment and self-calibration based on dual-axis rotating strapdown inertial navigation system | |
CN106767925A (zh) | 带双轴转位机构的惯导***三位置参数辨识对准方法 | |
CN106595669B (zh) | 一种旋转体姿态解算方法 | |
EP3214404B1 (en) | Systems and methods for calibrating and adjusting a heading reference system | |
CN105865444B (zh) | 基于仿射变换的惯性/地磁匹配迭代定位方法 | |
CN107860382A (zh) | 一种在地磁异常情况下应用ahrs测量姿态的方法 | |
CN103884356A (zh) | 一种标定捷联惯性组合陀螺仪组合的方法 | |
CN105865492A (zh) | 两轴磁强计在线误差补偿方法及*** | |
CN103954299B (zh) | 一种标定捷联惯性组合陀螺仪组合的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |