CN105301666B - 一种航磁干扰补偿系数的自适应调整方法 - Google Patents
一种航磁干扰补偿系数的自适应调整方法 Download PDFInfo
- Publication number
- CN105301666B CN105301666B CN201510751742.2A CN201510751742A CN105301666B CN 105301666 B CN105301666 B CN 105301666B CN 201510751742 A CN201510751742 A CN 201510751742A CN 105301666 B CN105301666 B CN 105301666B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- penalty coefficient
- magnetic disturbance
- mover
- 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
Landscapes
- Measuring Magnetic Variables (AREA)
Abstract
一种航磁干扰补偿系数的自适应调整方法，涉及航磁干扰补偿技术领域，具体涉及一种航磁干扰补偿系数的调整方法。为了解决现有的航磁干扰补偿技术中航磁干扰补偿系数求解矩阵存在病态性导致的航磁干扰补偿系数不准确的问题。本发明将飞机产生的磁干扰表示为HI；对T‑L模型进行延展，将地磁场HE与磁干扰HI进行融合；设定初始值θ0，P0，并根据公式θm＝θm‑1+Km(ym‑δmθm‑1)计算出第一次采样数据最终的补偿系数向量θ1；直至N次采样数据均处理完毕得到最终的补偿系数向量θN；令补偿系数向量的估计值以补偿系数向量的估计值作为航磁干扰补偿系数的自适应调整结果。本发明适用于磁干扰补偿系数的调整。
Description
技术领域
本发明涉及航磁干扰补偿技术领域，具体涉及一种航磁干扰补偿系数的调整方法。
背景技术
航磁干扰补偿技术是在航磁探测过程中去除航空运动平台自身磁干扰的一种技术。通过分析航空运动平台自身磁干扰的类型和性质，建立航空运动平台磁干扰数学模型，然后在校准飞行过程中按照规定的方法测得磁总场及三分量数据，并将其用来计算航空运动平台磁干扰数学模型的系数。在实际进行航磁探测时，利用估计出的模型系数及飞机姿态数据估计航空运动平台产生的磁干扰并将其从磁总场中去除，从而得到不含航空运动平台磁干扰的磁场数据。现有航磁干扰补偿系数计算方法多是基于T-L模型，该模型将航空运动平台磁干扰分为恒定场、感应场和涡流场三种类型，其中恒定场系数有3项，感应场系数和涡流场系数各有9项，且两者分别与地磁场的大小及变化率有关。由于可以把补偿系数的估计过程看作线性回归，因此系数求解矩阵的性质势必会影响估计结果的精度。其中，系数求解矩阵的病态性是影响补偿系数估计精度的关键因素。
发明内容
本发明为了解决现有的航磁干扰补偿技术中航磁干扰补偿系数求解矩阵存在病态性导致的航磁干扰补偿系数不准确的问题。进而提出了一种航磁干扰补偿系数的自适应调整方法。
一种航磁干扰补偿系数的自适应调整方法，包括以下步骤
步骤一、获得磁干扰HI与矩阵Δ：
在补偿飞行过程中，总场磁力仪测得的总场磁场数据为HT，是N×1列向量；地磁场为HE，是N×1列向量；飞机产生的磁干扰为HI，是N×1列向量；三者之间的关系是：
HT＝HE+HI (1)
航磁补偿的目的是通过测得的HT来确定HE，途径是通过磁干扰模型计算出HI；
T-L模型是当前被广泛使用的磁干扰模型，该模型将某一时刻飞机产生的磁干扰HI表示为：
其中，pi、aij、bij为需要估计出的航磁干扰补偿系数；ui为由地磁场与飞机轴向所成夹角的方向余弦，为ui的导数；i＝1,2,3，j＝1,2,3；
根据几何对称性，公式(2)可以转化为更为紧凑的形式：
其中，θ为补偿系数向量，是由需要估计出的航磁干扰补偿系数构成的列向量；δ是由ui、uiuj和构成的行向量；
利用N次采样数据，根据公式(4)，构造如下线性方程组：
HI＝Δθ (5)
其中，磁干扰HI是由所有N个采样时刻飞机产生的磁干扰HI构成的N×1列向量；为N×16维的矩阵；
步骤二、对T-L模型进行延展，将公式(1)中的地磁场HE与磁干扰HI进行融合而不影响干扰场的表达形式，故总场磁场数据HT表示为：
HT＝Δ′θ (6)
其中，Δ′为将HE与HI融合后的HT表示成与θ相乘对应的矩阵，即融合后的矩阵为Δ′；
由于存在噪声和其他频带信号的干扰，需要对HT与Δ进行滤波等预处理步骤，这里的HT与Δ都是预处理后数据；
步骤三、初始化补偿系数向量与协方差矩阵：
由于该算法为递推算法，在算法开始前需要给定相应的初值，初值的选取影响算法的收敛速度；
设定初始值为：
θ0＝Δ-1HT，P0＝105*E (7)
由于Δ在绝大多数情况下都是不可逆矩阵，因此此处Δ的逆是指伪逆，E为16×16维的单位对角阵；
步骤四、假设第m次测量得到的标量总场为ym，由ui、uiuj和构成的行向量为δm(1×16维)，由前m-1次测量递推所得的补偿系数向量为θm-1(16×1维)，根据公式(8)更新θ：
θm＝θm-1+Km(ym-δmθm-1) (8)
其中，Km为增益矩阵，大小为16×1维，用于修正每次递推得到的θm；
在N次采样数据中的第一次采样数据中，利用设定的初始值，根据公式(8)计算出第一次采样数据最终的补偿系数向量θ1；
步骤五、获得航磁干扰的补偿系数向量θN：
将第一次采样数据最终的补偿系数向量θ1作为第二次采样数据的初始值，根据步骤四计算出第二次采样数据最终的补偿系数向量θ2；
依次将前一次采样数据最终的补偿系数向量作为后一次采样数据的初始值，根据步骤四计算出补偿系数向量；
直至N次采样数据均处理完毕得到最终的补偿系数向量θN；
最终获得的补偿系数向量θN即是补偿系数向量的估计值，即：
其中，表示补偿系数向量θ的估计量；
以补偿系数向量的估计值作为航磁干扰补偿系数的自适应调整结果。
本发明具有以下有益效果：
本发明所提出的方法能够通过自适应的调整补偿系数来降低系数求解矩阵病态性造成的影响，得出稳定收敛的补偿结果；相比现有方法得到的航磁干扰补偿系数，本发明的航磁干扰补偿系数的准确度提高了52％以上，进而显著调高补偿精度；从而获得了相对准确的实际地磁场值。
附图说明
图1为本发明的流程图。
具体实施方式
具体实施方式一：结合图1说明本实施方式，
一种航磁干扰补偿系数的自适应调整方法，包括以下步骤
步骤一、获得磁干扰HI与矩阵Δ：
在补偿飞行过程中，总场磁力仪测得的总场磁场数据为HT，是N×1列向量；地磁场为HE，是N×1列向量；飞机产生的磁干扰为HI，是N×1列向量；三者之间的关系是：
HT＝HE+HI (1)
航磁补偿的目的是通过测得的HT来确定HE，途径是通过磁干扰模型计算出HI；
T-L模型是当前被广泛使用的磁干扰模型，该模型将某一时刻飞机产生的磁干扰HI表示为：
其中，pi、aij、bij为需要估计出的航磁干扰补偿系数；ui为由地磁场与飞机轴向所成夹角的方向余弦，为ui的导数；i＝1,2,3，j＝1,2,3；
根据几何对称性，公式(2)可以转化为更为紧凑的形式：
其中，θ为补偿系数向量，是由需要估计出的航磁干扰补偿系数构成的列向量；δ是由ui、uiuj和构成的行向量；
利用N次采样数据，根据公式(4)，构造如下线性方程组：
HI＝Δθ (5)
其中，磁干扰HI是由所有N个采样时刻飞机产生的磁干扰HI构成的N×1列向量；为N×16维的矩阵；
步骤二、对T-L模型进行延展，将公式(1)中的地磁场HE与磁干扰HI进行融合而不影响干扰场的表达形式，故总场磁场数据HT表示为：
HT＝Δ′θ (6)
其中，Δ′为将HE与HI融合后的HT表示成与θ相乘对应的矩阵，即融合后的矩阵为Δ′；
由于存在噪声和其他频带信号的干扰，需要对HT与Δ进行滤波等预处理步骤，这里的HT与Δ都是预处理后数据；
步骤三、初始化补偿系数向量与协方差矩阵：
由于该算法为递推算法，在算法开始前需要给定相应的初值，初值的选取影响算法的收敛速度；
设定初始值为：
θ0＝Δ-1HT，P0＝105*E (7)
由于Δ在绝大多数情况下都是不可逆矩阵，因此此处Δ的逆是指伪逆，E为16×16维的单位对角阵；
步骤四、假设第m次测量得到的标量总场为ym，由ui、uiuj和构成的行向量为δm(1×16维)，由前m-1次测量递推所得的补偿系数向量为θm-1(16×1维)，根据公式(8)更新θ：
θm＝θm-1+Km(ym-δmθm-1) (8)
其中，Km为增益矩阵，大小为16×1维，用于修正每次递推得到的θm；
在N次采样数据中的第一次采样数据中，利用设定的初始值，根据公式(8)计算出第一次采样数据最终的补偿系数向量θ1；
步骤五、获得航磁干扰的补偿系数向量θN：
将第一次采样数据最终的补偿系数向量θ1作为第二次采样数据的初始值，根据步骤四计算出第二次采样数据最终的补偿系数向量θ2；
依次将前一次采样数据最终的补偿系数向量作为后一次采样数据的初始值，根据步骤四计算出补偿系数向量；
直至N次采样数据均处理完毕得到最终的补偿系数向量θN；
最终获得的补偿系数向量θN即是补偿系数向量的估计值，即：
其中，表示补偿系数向量θ的估计量；
以补偿系数向量的估计值作为航磁干扰补偿系数的自适应调整结果。
具体实施方式二：
本实施方式步骤一中所述的ui是通过以下步骤实现的：
测量飞机三个轴向分量的磁场大小v1、v2、v3，然后根据公式(3)确定ui
其中，i＝1,2,3。
其它步骤和参数与具体实施方式一相同。
具体实施方式三：
本实施方式骤四中所述的Km的更新公式为：
其中，(·)T表示转置，Pm为协方差矩阵，大小为16×16维。
其它步骤和参数与具体实施方式一相同。
具体实施方式四：
本实施方式步骤四中所述的Pm的更新公式为：
其它步骤和参数与具体实施方式一相同。
Claims (4)
1.一种航磁干扰补偿系数的自适应调整方法，其特征在于包括以下步骤
步骤一、获得磁干扰HI与矩阵Δ：
在补偿飞行过程中，总场磁力仪测得的总场磁场数据为HT，是N×1列向量；地磁场为HE，是N×1列向量；磁干扰为HI，是N×1列向量；三者之间的关系是：
HT＝HE+HI (1)
T-L模型将某一时刻飞机产生的磁干扰HI表示为：
<mrow>
<msub>
<mi>H</mi>
<mi>I</mi>
</msub>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>3</mn>
</munderover>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>3</mn>
</munderover>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>3</mn>
</munderover>
<msub>
<mi>a</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
<msub>
<mi>u</mi>
<mi>j</mi>
</msub>
<mo>+</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>3</mn>
</munderover>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>3</mn>
</munderover>
<msub>
<mi>b</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>i</mi>
</msub>
<msub>
<mi>u</mi>
<mi>j</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，pi、aij、bij为需要估计出的航磁干扰补偿系数；ui为由地磁场与飞机轴向所成夹角的方向余弦，为ui的导数；i＝1,2,3，j＝1,2,3；
根据几何对称性，公式(2)可以转化为更为紧凑的形式：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>H</mi>
<mi>I</mi>
</msub>
<mo>=</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<msub>
<mi>u</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>p</mi>
<mn>3</mn>
</msub>
<msub>
<mi>u</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>11</mn>
</msub>
<msubsup>
<mi>u</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>12</mn>
</msub>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
<msub>
<mi>u</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>13</mn>
</msub>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
<msub>
<mi>u</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>22</mn>
</msub>
<msubsup>
<mi>u</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>23</mn>
</msub>
<msub>
<mi>u</mi>
<mn>2</mn>
</msub>
<msub>
<mi>u</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>11</mn>
</msub>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>12</mn>
</msub>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<msub>
<mi>u</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>13</mn>
</msub>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<msub>
<mi>u</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>21</mn>
</msub>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>22</mn>
</msub>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<msub>
<mi>u</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>23</mn>
</msub>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<msub>
<mi>u</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>31</mn>
</msub>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>3</mn>
</msub>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>32</mn>
</msub>
<msub>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>3</mn>
</msub>
<msub>
<mi>u</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mi>&delta;</mi>
<mi>&theta;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，θ为补偿系数向量，是由需要估计出的航磁干扰补偿系数构成的列向量；δ是由ui、uiuj和构成的行向量；
利用N次采样数据，根据公式(4)，构造如下线性方程组：
HI＝Δθ (5)
其中，磁干扰HI是由所有N个采样时刻飞机产生的磁干扰HI构成的N×1列向量；为N×16维的矩阵；
步骤二、对T-L模型进行延展，将公式(1)中的地磁场HE与磁干扰HI进行融合而不影响干扰场的表达形式，故总场磁场数据HT表示为：
HT＝Δ′θ (6)
其中，Δ′为将HE与HI融合后的HT表示成与θ相乘对应的矩阵，即融合后的矩阵为Δ′；
步骤三、初始化补偿系数向量与协方差矩阵：
设定初始值为：
θ0＝Δ-1HT，P0＝105*E (7)
Δ的逆是指伪逆，E为16×16维的单位对角阵；
步骤四、假设第m次测量得到的标量总场为ym，由ui、uiuj和构成的行向量为δm，由前m-1次测量递推所得的补偿系数向量为θm-1，根据公式(8)更新θ：
θm＝θm-1+Km(ym-δmθm-1) (8)
其中，Km为增益矩阵，用于修正每次递推得到的θm；
在N次采样数据中的第一次采样数据中，利用设定的初始值，根据公式(8)计算出第一次采样数据最终的补偿系数向量θ1；
步骤五、获得航磁干扰的补偿系数向量θN：
将第一次采样数据最终的补偿系数向量θ1作为第二次采样数据的初始值，根据步骤四计算出第二次采样数据最终的补偿系数向量θ2；
依次将前一次采样数据最终的补偿系数向量作为后一次采样数据的初始值，根据步骤四计算出补偿系数向量；
直至N次采样数据均处理完毕得到最终的补偿系数向量θN；
最终获得的补偿系数向量θN即是补偿系数向量的估计值，即：
<mrow>
<mover>
<mi>&theta;</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<msub>
<mi>&theta;</mi>
<mi>N</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，表示补偿系数向量θ的估计量；
以补偿系数向量的估计值作为航磁干扰补偿系数的自适应调整结果。
2.根据权利要求1所述的一种航磁干扰补偿系数的自适应调整方法，其特征在于步骤一中所述的ui是通过以下步骤实现的：
测量飞机三个轴向分量的磁场大小v1、v2、v3，然后根据公式(3)确定ui：
<mrow>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>v</mi>
<mi>i</mi>
</msub>
<msqrt>
<mrow>
<msubsup>
<mi>v</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>v</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>v</mi>
<mn>3</mn>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，i＝1,2,3。
3.根据权利要求1或2所述的一种航磁干扰补偿系数的自适应调整方法，其特征在于步骤四中所述的Km的更新公式为：
Km＝Pm-1δm T(1+δmPm-1δm T)-1 (9)
其中，(·)T表示转置，Pm为协方差矩阵。
4.根据权利要求3所述的一种航磁干扰补偿系数的自适应调整方法，其特征在于步骤四中所述的Pm的更新公式为：
Pm＝Pm-1-Pm-1δm T(1+δmPm-1δm T)-1δmPm-1 (10)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510751742.2A CN105301666B (zh) | 2015-11-05 | 2015-11-05 | 一种航磁干扰补偿系数的自适应调整方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510751742.2A CN105301666B (zh) | 2015-11-05 | 2015-11-05 | 一种航磁干扰补偿系数的自适应调整方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105301666A CN105301666A (zh) | 2016-02-03 |
| CN105301666B true CN105301666B (zh) | 2017-09-26 |
Family
ID=55199141
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510751742.2A Active CN105301666B (zh) | 2015-11-05 | 2015-11-05 | 一种航磁干扰补偿系数的自适应调整方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105301666B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106707352B (zh) * | 2016-11-28 | 2018-10-12 | 中国科学院电子学研究所 | 一种基于ε-SVR算法的航磁梯度干扰的去除方法 |
| CN106842344B (zh) * | 2017-04-24 | 2018-10-12 | 中国科学院电子学研究所 | 基于前馈网络的无人机航磁全轴梯度磁干扰补偿的方法 |
| CN107356978B (zh) * | 2017-07-11 | 2019-03-01 | 中国科学院电子学研究所 | 基于主成分分析的航磁补偿方法 |
| WO2019061403A1 (zh) * | 2017-09-30 | 2019-04-04 | 中国科学院电子学研究所 | 基于航磁补偿误差模型的磁补偿方法 |
| CN108919363B (zh) * | 2018-08-07 | 2019-09-13 | 哈尔滨工业大学 | 一种自适应依航向航磁数据处理方法 |
| CN109444774B (zh) * | 2018-11-30 | 2020-11-24 | 中国人民解放军海军潜艇学院 | 基于小信号法测量水下航行器磁场噪声系数的方法 |
| CN110333552B (zh) * | 2019-07-15 | 2020-09-25 | 中国地质大学（武汉） | 一种基于Liu估计的航磁补偿方法 |
| CN112858960B (zh) * | 2021-02-28 | 2022-07-22 | 哈尔滨工业大学 | 一种飞机尾杆振动引起的磁干扰的补偿方法及装置 |
| CN112858961B (zh) * | 2021-03-01 | 2023-02-07 | 哈尔滨工业大学 | 一种针对航空平台多源磁干扰的补偿方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2484422A1 (en) * | 2004-10-08 | 2006-04-08 | Furgro Airborne Surveys | Unmanned airborne vehicle for geophysical surveying |
| CN103577607B (zh) * | 2013-11-20 | 2017-06-20 | 哈尔滨工程大学 | 一种基于地磁异常数据形态特征的边界补偿方法 |
| CN103837845A (zh) * | 2014-01-22 | 2014-06-04 | 哈尔滨工程大学 | 一种飞行器磁干扰场模型参数解算方法 |
| CN104215259B (zh) * | 2014-08-22 | 2018-04-24 | 哈尔滨工程大学 | 一种基于地磁模量梯度和粒子滤波的惯导误差校正方法 |
-
2015
- 2015-11-05 CN CN201510751742.2A patent/CN105301666B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105301666A (zh) | 2016-02-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105301666B (zh) | 一种航磁干扰补偿系数的自适应调整方法 | |
| CN105425304B (zh) | 一种飞机航磁干扰的补偿方法 | |
| CN106487358B (zh) | 一种机动目标转弯跟踪方法 | |
| CN105222809B (zh) | 一种地磁梯度鲁棒的航磁干扰补偿系数估计的方法 | |
| CN104091064B (zh) | 基于优化解空间搜索法的PS‑DInSAR地表形变测量参数估计方法 | |
| CN107688554B (zh) | 基于自适应傅里叶分解的频域辨识方法 | |
| CN105510849A (zh) | 航磁干扰补偿方法 | |
| CN110174643B (zh) | 一种无需噪声功率信息的基于到达时间差的定位方法 | |
| CN105509737A (zh) | 一种不受地磁变化影响的航空运动平台磁干扰补偿方法 | |
| CN111989265B (zh) | 超低轨道卫星轨道自主维持方法 | |
| CN107991676B (zh) | 星载单航过InSAR系统对流层误差校正方法 | |
| CN110187393B (zh) | 一种基于广义回归神经网络的航磁补偿方法 | |
| CN106707352B (zh) | 一种基于ε-SVR算法的航磁梯度干扰的去除方法 | |
| Zhang et al. | A simplified aeromagnetic compensation model for low magnetism UAV platform | |
| CN108318219B (zh) | 自由飞试验相似律重力补偿设计方法 | |
| CN104076332A (zh) | 一种雷达均匀线性阵列幅度和相位的估计方法 | |
| CN105547326A (zh) | 陀螺与磁传感器联合标定方法 | |
| CN105466458A (zh) | 正六面体磁梯度张量系统误差直接校正方法 | |
| CN104048676A (zh) | 基于改进粒子滤波的mems陀螺随机误差补偿方法 | |
| CN110244367B (zh) | 一种基于地面多基站的ztem系统姿态补偿方法 | |
| CN110632555A (zh) | 一种基于矩阵特征值扰动的tdoa直接定位方法 | |
| Morelli | Dynamic modeling from flight data with unknown time skews | |
| CN112213678B (zh) | 一种矢量磁探仪三轴数据校正与补偿方法 | |
| CN110186483B (zh) | 提高惯性制导航天器落点精度的方法 | |
| CN104535082B (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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |