CN112284377A - 一种应用于飞行器的地磁场测量系统及方法 - Google Patents

一种应用于飞行器的地磁场测量系统及方法 Download PDF

Info

Publication number
CN112284377A
CN112284377A CN202011163696.1A CN202011163696A CN112284377A CN 112284377 A CN112284377 A CN 112284377A CN 202011163696 A CN202011163696 A CN 202011163696A CN 112284377 A CN112284377 A CN 112284377A
Authority
CN
China
Prior art keywords
aircraft
geomagnetic
magnetic field
model
field
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.)
Pending
Application number
CN202011163696.1A
Other languages
English (en)
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.)
Naval University of Engineering PLA
Original Assignee
Naval University of Engineering PLA
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 Naval University of Engineering PLA filed Critical Naval University of Engineering PLA
Priority to CN202011163696.1A priority Critical patent/CN112284377A/zh
Publication of CN112284377A publication Critical patent/CN112284377A/zh
Pending legal-status Critical Current

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/04Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means
    • G01C21/08Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means involving use of the magnetic field of the earth
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C25/00Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Manufacturing & Machinery (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Automation & Control Theory (AREA)
  • Measuring Magnetic Variables (AREA)

Abstract

本发明提供了一种应用于飞行器的地磁场测量系统及方法,该方法包括如下步骤:包括如下步骤:S1,采用三轴磁场矢量传感器获取飞行器在地磁背景中的姿态信息,并建立飞行器的背景磁干扰补偿模型;S2,采用光泵磁强计测量地磁总强度;S3,利用飞行器背景磁干扰补偿模型,对飞行器背景磁场进行补偿,消除飞行器载体剩磁、感应磁场和涡流磁场对测量结果的影响,实现对地磁场的精确测量。本发明另一方面提供了一种应用于飞行器的地磁场测量系统,包括三轴磁场矢量传感器及光泵磁强计,所述应用于飞行器的地磁场测量系统根据权利要求上述的应用于飞行器的地磁场测量方法来实现飞行器的地磁场测量。

Description

一种应用于飞行器的地磁场测量系统及方法
技术领域
本发明涉及地磁场测量技术领域,具体涉及一种应用于飞行器的地磁场测量系统及方法。
背景技术
飞行器内部的磁场环境是非常复杂的,空间梯度大,严重影响到内部的地磁场测量装置的测量准确度,导致在飞行器导航定位时容易出现较大误差,造成极大的安全隐患。因此,有必要改进飞行器的地磁场测量方法,构建地磁场总量的高精度测量系统,用于在飞行器导航定位系统中实时测量地磁场数据。
发明内容
为了解决上述问题,本发明的目的在于提供一种应用于飞行器的地磁场测量方法,包括如下步骤:
S1,采用三轴磁场矢量传感器获取飞行器在地磁背景中的姿态信息,并建立飞行器的背景磁干扰补偿模型;
S2,采用光泵磁强计测量地磁总强度;
S3,利用飞行器背景磁干扰补偿模型,对飞行器背景磁场进行补偿,消除飞行器载体剩磁、感应磁场和涡流磁场对测量结果的影响,实现对地磁场的精确测量。
在一些实施例中,所述光泵磁强计采用铯光泵磁强计。
在一些实施例中,步骤S1中,建立飞行器的背景磁干扰补偿模型时,先分别建立飞行器背景磁场的剩磁模型、感应磁场模型及涡流磁场模型,并将三者结合作为飞行器总的背景磁干扰补偿模型。
在一些实施例中,步骤S1中,建立剩磁模型的步骤包括:
将剩磁表示为:Hp={p1 p2 p3}
则在测量点沿地磁矢量方向的投影为:
Figure BDA0002745052590000021
其中ui=cosθi为地磁矢量的方向余弦;
建立感应磁场模型的步骤包括:
将感应磁场在测量点处的矢量表示为:
Figure BDA0002745052590000022
沿地磁矢量方向的投影为:
Figure BDA0002745052590000023
其中:
Figure BDA0002745052590000024
i≠j时
a11=(i11-i33)B
a22=(i22-i33)B
令a33=0,此时
Figure BDA0002745052590000025
建立涡流磁场模型的步骤包括:
将感应磁场在测量点处的矢量表示为:
Figure BDA0002745052590000031
根据方向余弦的性质,有:
Figure BDA0002745052590000032
u1u′1+u2u′2+u3u′3=0
令:
Figure BDA0002745052590000033
则涡流磁场沿地磁矢量方向的投影可表达成:
Figure BDA0002745052590000034
因此,飞行器总的背景磁干扰补偿模型为:
Figure BDA0002745052590000035
本发明另一方面提供了一种应用于飞行器的地磁场测量系统,包括三轴磁场矢量传感器及光泵磁强计,所述应用于飞行器的地磁场测量系统根据上述的应用于飞行器的地磁场测量方法来实现飞行器的地磁场测量。
与现有技术相比,本发明提供的应用于飞行器的地磁场测量系统及方法,采用三轴磁场矢量传感器与光泵磁强计的结合,并利用飞行器背景磁干扰补偿模型,对飞行器背景磁场进行补偿,以消除飞行器载体剩磁、感应磁场和涡流磁场对测量结果的影响,实现了对地磁场的精确测量。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将对本发明的技术方案进行清楚、完整地描述。显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于所描述的本发明的实施例,本领域普通技术人员在无需创造性劳动的前提下所获得的所有其他实施例,都属于本发明的保护范围。
除非另外定义本发明使用的技术术语或者科学术语应当为本发明所属领域内具有一般技能的人士所理解的通常意义。
本发明提供了一种应用于飞行器的地磁场测量方法,包括如下步骤:
S1,采用三轴磁场矢量传感器获取飞行器在地磁背景中的姿态信息,并建立飞行器的背景磁干扰补偿模型;
S2,采用光泵磁强计测量地磁总强度;
S3,利用飞行器背景磁干扰补偿模型,对飞行器背景磁场进行补偿,消除飞行器载体剩磁、感应磁场和涡流磁场对测量结果的影响,实现对地磁场的精确测量。
优选地,光泵磁强计采用铯光泵磁强计。
进一步地,步骤S1中,建立飞行器的背景磁干扰补偿模型时,先分别建立飞行器背景磁场的剩磁模型、感应磁场模型及涡流磁场模型,并将三者结合作为飞行器总的背景磁干扰补偿模型。
具体地,步骤S1中,建立剩磁模型的步骤包括:
将剩磁表示为:Hp={p1 p2 p3}
则在测量点沿地磁矢量方向的投影为:
Figure BDA0002745052590000041
其中ui=cosθi为地磁矢量的方向余弦;
建立感应磁场模型的步骤包括:
将感应磁场在测量点处的矢量表示为:
Figure BDA0002745052590000042
沿地磁矢量方向的投影为:
Figure BDA0002745052590000051
其中:
Figure BDA0002745052590000052
i≠j时
a11=(i11-i33)B
a22=(i22-i33)B
令a33=0,此时
Figure BDA0002745052590000053
建立涡流磁场模型的步骤包括:
将感应磁场在测量点处的矢量表示为:
Figure BDA0002745052590000054
根据方向余弦的性质,有:
Figure BDA0002745052590000055
u1u′1+u2u′2+u3u′3=0
令:
Figure BDA0002745052590000056
则涡流磁场沿地磁矢量方向的投影可表达成:
Figure BDA0002745052590000061
因此,飞行器总的背景磁干扰补偿模型为:
Figure BDA0002745052590000062
可以理解的是,建立起上述的背景磁干扰补偿模型后,可利用仿真试验的方法,并生成关于模型参数的回归方程组,进而解出模型中的各个参数。
本发明另一方面提供了一种应用于飞行器的地磁场测量系统,包括三轴磁场矢量传感器及光泵磁强计,应用于飞行器的地磁场测量系统根据上述的应用于飞行器的地磁场测量方法来实现飞行器的地磁场测量。
光泵磁强计探头是地磁标量测量系统的核心部件,光泵磁强计探头将外部磁场信号转变成频率信号。目前,工程上常用的光泵磁强计探头主要有氦光泵磁强计探头和铯光泵磁强计探头两类,氦光泵磁强计探头内部工作室的工质采用氦气,在常温下是气体,且分子量小,容易泄露,因此其使用寿命较短,一般1-2年后需要重新补充氦气。铯光泵磁强计探头内部工作室的工质采用铯蒸汽,在常温下铯是固体,不会泄露,在使用时,才将其加热成铯蒸汽。铯光泵磁强计探头的使用寿命长,一般在10年以上。本发明在采用铯光泵磁强计探头的同时,结合高分辨率、高精度的FPGA鉴频器和信号采集系统,共同构成铯光泵磁强计。所采用的铯光泵磁强计的主要技术指标如下:
1)测量范围:15000-105000nT;
2)灵敏度:
Figure BDA0002745052590000063
3)梯度容限:4000nT/m;
4)方向误差:±0.25nT
5)全场绝对精度:±2.5nT
6)信号输出:3.49857Hz/nT
这些指标能完全满足地磁导航定位系统对测量精度的要求。该铯光泵磁强计铯的工作室的体积为1cm3,磁场梯度达每米4000nT,这一梯度容限基本能保证磁强计可在飞行器内部安装使用,对于长度小于50cm,直径不大于10cm的一般铁制物体,只要放置距离不小于30cm,其磁场梯度均小于磁探头的磁场梯度容限,光泵磁强计能正常工作。而对于一些磁性很强的永磁体和电磁体,只要其外部有铁制磁扼,则磁体磁场绝大部分被限制在磁扼内部,而磁扼外部的磁场分布与磁扼物体的磁场分布基本相同,因此可按铁磁性物体的质量和形状估计磁场梯度。
本发明中的FPGA鉴频器使用稳定度为10-8的40MHz恒温晶振提供基准频率,以高性能FPGA为控制部件,利用FPGA内部的PLL(数字锁相环)和DCM(数字时钟管理器),将时钟扩频到400M,由于同时保证了时钟频率和稳定度,使得由此构成的光泵磁强计在保证高分辨率、高精度的条件下,具有较高的信号响应速度(信号带宽)。
系统的频率分辨率及对应信号带宽估计如下:
直接鉴频系统频率分辨率的计算公式为:
Figure BDA0002745052590000071
其中:fs为光泵信号频率;
N为计数时间内的高频脉冲数,有:
N=T.fn,T为计数时间。
所以,有:
Figure BDA0002745052590000072
fT为信号采样频率。
不同采样频率下的测量系统分辨率如表1所示:
表1
Figure BDA0002745052590000073
表中取fs=170KHz,fN=400MHz,对应外磁场约48600nT;铯133的磁旋系数为3.498577Hz/nT。
光泵磁强计探头采用热曼效应原理工作,其信号载体为光波,频率高达数GHz,通过光谱的频率跃迁频率来测量外部磁场,从信号检测的机理来看,光泵磁强计探头测量机理实际上是测量外部磁场信号对光波信号的频率调制,由于光波频率远高于工业领域的无线电波频率,因此各种电磁辐射干扰不会因混叠效应而引起探头对低频磁场信号测量的干扰,电磁辐射干扰可以通过平滑滤波予于消除。
选用的Cs-3铯光泵探头,其频带宽度达到5MHz,在地磁测量领域具有足够的宽度,其对测量系统的影响可以不必考虑,即测量系统的频率响应特性只取决于FPGA高分辨率鉴频器的频率响应特性。
FPGA鉴频器对铯光泵磁强计探头的输出频率信号进行频率鉴别,采用直接频率鉴别方法,其工作原理是,在确定时间间隔tw内,记录铯光泵磁强计探头的输出信号周期数N,则光泵输出信号的频率为:
Figure BDA0002745052590000081
对应的磁场测量结果为:
Figure BDA0002745052590000082
其中,
Figure BDA0002745052590000083
为铯分子的旋磁比。
直接鉴频原理并不能测量任意时刻的磁场频率瞬时值、从而获得对外磁场信号的即时测量,而是测量tw时间内的平均频率,从而测量tw时间内的平均磁场,经过鉴频电路后,t时刻输出的tw时间内的平均磁场为:
Figure BDA0002745052590000084
其中:
Figure BDA0002745052590000091
为矩形窗函数,窗宽为tw
输出信号Bs(t)的频谱特性为:
Bs(ω)=Uw(ω)B(ω)
因此,鉴频电路的频率响应特性就是宽度为tw的矩形窗的频率响应特性。
Figure BDA0002745052590000092
当选择窗宽tw=0.04s时,其直流增益为1,当信号频率为10Hz时,增益为0.76,因此,鉴频电路的信号带宽约为10Hz。为了保证满足采样定理,信号采样间隔最大值为DT=0.04s,即采样间隔和窗宽可以相等。
同时改变DT和tw值,可以改变信号带宽,从表1可以看出,当采样率为25Hz,信号带宽约为10Hz时,磁场分辨率优于0.003nT(对应外磁场48600nT时)。在一个具体实施例中,研制成的硬件系统的设定信号带宽为10Hz,对应磁场分辨率优于0.01nT(当外磁场不大于10,0000nT时)。
综上,本发明提供的应用于飞行器的地磁场测量系统及方法,采用三轴磁场矢量传感器与光泵磁强计的结合,并利用飞行器背景磁干扰补偿模型,对飞行器背景磁场进行补偿,以消除飞行器载体剩磁、感应磁场和涡流磁场对测量结果的影响,实现了对地磁场的精确测量。
以上所描述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应在本发明的保护范围之内。

Claims (5)

1.一种应用于飞行器的地磁场测量方法,其特征在于,包括如下步骤:
S1,采用三轴磁场矢量传感器获取飞行器在地磁背景中的姿态信息,并建立飞行器的背景磁干扰补偿模型;
S2,采用光泵磁强计测量地磁总强度;
S3,利用飞行器背景磁干扰补偿模型,对飞行器背景磁场进行补偿,消除飞行器载体剩磁、感应磁场和涡流磁场对测量结果的影响,实现对地磁场的精确测量。
2.根据权利要求1所述的应用于飞行器的地磁场测量方法,其特征在于,所述光泵磁强计采用铯光泵磁强计。
3.根据权利要求1所述的应用于飞行器的地磁场测量方法,其特征在于,步骤S1中,建立飞行器的背景磁干扰补偿模型时,先分别建立飞行器背景磁场的剩磁模型、感应磁场模型及涡流磁场模型,并将三者结合作为飞行器总的背景磁干扰补偿模型。
4.根据权利要求3所述的应用于飞行器的地磁场测量方法,其特征在于,步骤S1中,建立剩磁模型的步骤包括:
将剩磁表示为:Hp={p1 p2 p3}
则在测量点沿地磁矢量方向的投影为:
Figure FDA0002745052580000011
其中ui=cosθi为地磁矢量的方向余弦;
建立感应磁场模型的步骤包括:
将感应磁场在测量点处的矢量表示为:
Figure FDA0002745052580000021
沿地磁矢量方向的投影为:
Figure FDA0002745052580000022
其中:
Figure FDA0002745052580000023
i≠j时
a11=(i11-i33)B
a22=(i22-i33)B
令a33=0,此时
Figure FDA0002745052580000024
建立涡流磁场模型的步骤包括:
将感应磁场在测量点处的矢量表示为:
Figure FDA0002745052580000025
根据方向余弦的性质,有:
Figure FDA0002745052580000026
u1u′1+u2u′2+u3u′3=0
令:
Figure FDA0002745052580000031
则涡流磁场沿地磁矢量方向的投影可表达成:
Figure FDA0002745052580000032
因此,飞行器总的背景磁干扰补偿模型为:
Figure FDA0002745052580000033
5.一种应用于飞行器的地磁场测量系统,其特征在于,包括三轴磁场矢量传感器及光泵磁强计,所述应用于飞行器的地磁场测量系统根据权利要求1-4任一项所述的应用于飞行器的地磁场测量方法来实现飞行器的地磁场测量。
CN202011163696.1A 2020-10-27 2020-10-27 一种应用于飞行器的地磁场测量系统及方法 Pending CN112284377A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011163696.1A CN112284377A (zh) 2020-10-27 2020-10-27 一种应用于飞行器的地磁场测量系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011163696.1A CN112284377A (zh) 2020-10-27 2020-10-27 一种应用于飞行器的地磁场测量系统及方法

Publications (1)

Publication Number Publication Date
CN112284377A true CN112284377A (zh) 2021-01-29

Family

ID=74372569

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011163696.1A Pending CN112284377A (zh) 2020-10-27 2020-10-27 一种应用于飞行器的地磁场测量系统及方法

Country Status (1)

Country Link
CN (1) CN112284377A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114184988A (zh) * 2021-11-11 2022-03-15 北京大学 包含补偿平台电流磁干扰的航磁补偿方法及其装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106959471A (zh) * 2017-04-21 2017-07-18 中国科学院电子学研究所 基于非线性航磁总场梯度补偿模型的航磁补偿方法
CN109444774A (zh) * 2018-11-30 2019-03-08 中国人民解放军海军潜艇学院 基于小信号法测量水下航行器磁场噪声系数的方法
CN110736484A (zh) * 2019-11-18 2020-01-31 北京理工大学 基于陀螺仪及磁传感器融合的背景磁场标定方法
CN111433634A (zh) * 2017-09-30 2020-07-17 中国科学院空天信息创新研究院 基于航磁补偿误差模型的磁补偿方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106959471A (zh) * 2017-04-21 2017-07-18 中国科学院电子学研究所 基于非线性航磁总场梯度补偿模型的航磁补偿方法
CN111433634A (zh) * 2017-09-30 2020-07-17 中国科学院空天信息创新研究院 基于航磁补偿误差模型的磁补偿方法
CN109444774A (zh) * 2018-11-30 2019-03-08 中国人民解放军海军潜艇学院 基于小信号法测量水下航行器磁场噪声系数的方法
CN110736484A (zh) * 2019-11-18 2020-01-31 北京理工大学 基于陀螺仪及磁传感器融合的背景磁场标定方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
FU MENGYIN 等: "Aeromagnetic Compensation Based on Tikhonov Regularization with Limited L-curve Parameter-choice Algorithm", 《PROCEEDINGS OF THE 37TH CHINESE CONTROL CONFERENCE》 *
S.H. BICKEL: "Small Signal Compensation of Magnetic Fields Resulting from Aircraft Maneuvers", 《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》 *
ZHANG NING等: "tudy on Model Compernsation of Aircraft magnetic interference", 《2010 IEEE INTERNATIONAL CONFERENCE ON INTELLIGENT COMPUTING AND INTELLIGENT SYSTEMS》 *
杨春玲等: "《现代可编程逻辑器件及SOPC应用设计》", 31 March 2005, 哈尔滨工业大学出版社 *
邹小金等: "《EDA应用技术》", 31 July 2011, 西苑出版社 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114184988A (zh) * 2021-11-11 2022-03-15 北京大学 包含补偿平台电流磁干扰的航磁补偿方法及其装置

Similar Documents

Publication Publication Date Title
Steiner et al. Uncertainty improvements of the NIST electronic kilogram
Griffith et al. Miniature atomic magnetometer integrated with flux concentrators
Ripka Advances in fluxgate sensors
Brown et al. New limit on lorentz-and C P T-violating neutron spin interactions
Dong et al. A high-precision frequency measurement algorithm for FID signal of proton magnetometer
Schönau et al. SQUID-based setup for the absolute measurement of the Earth’s magnetic field
CN110118948A (zh) 一种基于超导量子干涉仪的总场测量方法及装置
CN111551163B (zh) 四极核旋转边带惯性转动测量方法和三轴nmr陀螺装置
CN108717168A (zh) 一种基于光场幅度调制的标量磁场梯度测量装置及方法
CN110274586A (zh) 包含多光系原子磁力仪方向误差补偿的航空磁补偿方法
Liu et al. Construction of an Overhauser magnetic gradiometer and the applications in geomagnetic observation and ferromagnetic target localization
van Niekerk et al. A wireless radio frequency triggered acquisition device (WRAD) for self-synchronised measurements of the rate of change of the MRI gradient vector field for motion tracking
Sapunov et al. Ground overhauser DNP geophysical devices
Tan et al. Nonlinear compensation method based on data for increasing absolute measurement precision of FID signal
CN112284377A (zh) 一种应用于飞行器的地磁场测量系统及方法
Wan et al. Improved component compensation for geomagnetic field vector measurement using Lagrange multiplier method
Zikmund et al. Precise scalar calibration of a tri-axial Braunbek coil system
Kawai et al. Characterization and demonstration results of a SQUID magnetometer system developed for geomagnetic field measurements
CN108919368A (zh) 一种用于消除微小卫星剩磁干扰的系统及方法
CN112798993B (zh) 基于加速度计测量永磁材料温度系数的装置及测量方法
Dong et al. A high-accuracy and non-intermittent frequency measurement method for Larmor signal of optically pumped cesium magnetometer
Tan et al. A frequency measurement method using rising-falling edge of square wave for increasing proton magnetometer precision
Yu et al. AC vector magnetometer for space-based applications using low-resource magneto-impedance sensor
Zhi et al. Design and analysis of miniature tri-axial fluxgate magnetometer
CN109839610B (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
RJ01 Rejection of invention patent application after publication

Application publication date: 20210129

RJ01 Rejection of invention patent application after publication