運動狀態下磁場模量修正技術研究*

楊福杏,曹大平

(武漢大學物理科學與技術學院,武漢 430072)

運動狀態下磁場模量修正技術研究*

楊福杏,曹大平*

(武漢大學物理科學與技術學院,武漢 430072)

考慮到現有三軸磁通門傳感器的修正存在方法復雜、對設備要求較高、效果不佳等問題。根據三軸磁通門傳感器是由于三軸不正交、靈敏度不一致和零點偏移所引起的誤差進行分析和理論計算,提出了三軸磁通門傳感器修正模型并運用奇異值分解和單純型法進行實驗。實驗結果表明,該方法能夠簡單、有效的對磁場模量進行修正,能夠廣泛的應用于三軸磁通門傳感器靜態和動態的總模量修正。

三軸磁通門傳感器;誤差修正;奇異值分解;單純形

三軸磁通門傳感器是一門比較成熟的技術,用于地磁導航、礦產探測衛星姿態、艦船磁測和武器彈丸等多個領域[1-7]。但是實際上,由于加工工藝和安裝水平的限制,3個磁軸不可能嚴格正交,三軸靈敏度以及其他電氣性能也不可能完全對稱,因此還會存在靈敏度誤差和零偏誤差。由于這些誤差的引入,往往導致測量的精度在應用時是難以接受的。因此,三軸磁通門傳感器在使用之前必須進行修正。目前修正的方法有變換系數法[8],高斯牛頓迭代算法[9]、FLANN[10]、自適應遺傳算法[11]和橢球擬合[12]等。然而,這些方法有些不足之處:①算法本身具有誤差,導致結果不是很理想;②有些計算過于復雜,不利于實際應用;③靜態采樣數據,沒有對動態過程進行分析。本文基于正交變換原理,建立誤差模型,運用非線性最小二乘法、奇異值分解法、單純型法解得最優解,進而求得理想輸出。由于三軸磁傳感器在長期在空中或水下自主運行時,傳感器的參數隨時間會有緩慢漂移,需要進行動態修正,且傳感器傾斜較小,可視為Z軸繞垂直方向在較小的范圍變化,原來用于三維求解的方法不在適用,本文圍繞這一問題展開研究工作。

1 三軸磁通門傳感器靜態修正模型建立

根據理想磁通門傳感器測量環境磁場強度與姿態無關的特點,當三軸磁傳感器繞同一個軸轉動時,可以認為測的的磁場總的模量不變?；谶@個特點提出了修正模型:

B2=a1(Bx-Bx0)2+a2(By-By0)2+a3(Bz-Bz0)2+

a4(Bx-Bx0)(By-By0)+a5(By-By0)(Bz-Bz0)+

a6(Bz-Bz0)(Bx-Bx0)

(1)

式中:a1、a2、a3為X、Y、Z的輸出靈敏度系數;a4、a5、a6為XY、YZ、ZX之間的相互影響系數;Bx0、By0、Bz0為X、Y、Z的輸出零點;總磁場模量為B。利用實驗所測得的Bx、By、Bz,運用奇異值分解和單純形法對式(1)進行非線性擬合,可得到a1～a6及3個輸出零點Bx0、By0、Bz0的值。

在理想狀態下,磁傳感器3個分量應為彼此垂直,但實際安裝中不可能達到。設i、j、k為理想參考正交單位矢量,x、y、z為磁傳感器3個分量的軸向,如圖1所示。

圖1 參考坐標系與磁探頭軸向示意圖

取x方向與i方向一致,j在xy平面內與i方向垂直,k垂直于xy平面,則有:

i=d1x

(2)

j=k×x=b1x+b2y

(3)

(4)

在參考坐標系中所測得的磁場總矢量為:

B=Bxri+Byrj+Bzrk

(5)

式中:Bxr、Byr、Bzr為理想狀態下三軸輸出,B為磁場總量。

將式(2)～式(4)代入式(5)既有:

B=Bxrd1x+Byr(b1x+b2y)+Bzr(c1x+c2y+c3z)=

(Bxrd1+Byrb1+Bzrc1)x+(Byrb2+Bzrc2)y+Bzrc3z

(6)

考慮到三軸磁通門傳感器3個分量的輸出零點Bx0、By0、Bz0,又有:

B=(Bx-Bx0)x+(By-By0)y+(Bz-Bz0)z

(7)

對比式(6)與式(7),可得:

(8)

U即為三軸不正交的修正矩陣。

由式(4)可得,磁場模量為:

(9)

將式(8)代入式(1),可得:

BxrByr+(2a1b1c1+2b2c2a2+a4b1c2+a4c1b2+a5b2c3+

a6c3b1)ByrBzr+(2a1d1c1+a4d1c2+a6d1c3)BxrBzr

(10)

(11)

解式(11)可得:

(12)

綜上可得,三軸正交的修正矩陣為:

(13)

結合運用非線性擬合取得的9個參數,即可得到靜態修正后的結果。

2 運動狀態下三軸磁通門傳感器的修正

當三軸磁通門傳感器應用于水面艦船磁場的測量時,通常狀態下z軸磁場變化很小,屬于二維情況,這時只能求出與x、y相關的5個參數,即a1、a2、a4及Bx0、By0。若z軸變化很大時屬于三維情況,其求解方法與二維情況類似。這里我們只討論二維情況下的動態修正,為此將靜態的磁場模量表達式(1)修改為:

B′2=a1(Bx-Bx0)2+a2(By-By0)2+

a4(Bx-Bx0)(By-By0)

(14)

式中:

B′2=B2-a3(Bz-Bz0)2-a5(By-By0)(Bz-Bz0)-

a6(Bz-Bz0)(Bx-Bx0)

(15)

且式(15)中的參數全部用上一部分求得的靜態參數代替,這里的B是根據三軸磁通門傳感器實際采集的磁場模量的平均值作為修正的目標。然后,根據采取得樣點數,組成如下方程組:

(16)

與靜態方法一樣通過單純形與奇異值分解相結合的方法求出a1、a2、a4及Bx0、By0,將相應的原有參數命名為a1O、a2O、a4O及Bx0O、By0O,選擇參數更新因子λ=0.05,則,更新后的參數為:

(17)

將更新后的參數代入修正矩陣,得到修正后的結果。隨著三軸磁通門傳感器在某一小范圍,輕微的擾

動,修正參數不斷的進行更新以達到動態修正的效果。

3 實驗過程及結果

3.1 靜態修正實驗及結果

本實驗在武漢大學珞珈山弱磁實驗室進行,磁場穩定,受外界干擾較小。由于不可能遍歷所有方向,在試驗中我們選取26各具有代表性的點如圖2所示,讓他代表所有的球空間分布。

運用三軸無磁轉臺,結合三軸磁傳感器個分量輸出確定點的大致位置,并將此點作為一個樣本點,共26個。

圖2 樣本點位置

實驗結果如表1、表2和圖3所示。

表1 靜態修正

表2 平均值

圖3 靜態修正前后數據對比

3.2 動態修正實驗及結果

實驗場所位于珞珈山上,周圍磁場干擾較小,將三軸磁通門傳感器懸吊在半空中,z軸豎直,傳感器下面加上重物以減小傳感器運動的幅度,來模擬二維狀態下的微擾。通過數據采集模塊將數據采集。實驗結果如圖4所示。

4 結論

針對三軸磁通門傳感器三軸不正交,零點偏置和靈敏度不同建立了修正模型。由于奇異值分[13-14]解精確穩定,本身并不引入誤差而且在病態條件下也有很好的穩定性。單純型法不需要函數的一階與二階導數,僅需要函數本身的值,計算簡潔,數值穩定性高,迭代過程清晰,便于分析。運用奇異值分解和單純形法的結合進行非線性擬合,取得最優化參數。通過靜態修正磁場模量最大偏差1 561 nT校準到28 nT,經過動態修正后磁場模量最大偏差值由282 nT減小到4 nT。本文不僅在靜態條件下給出了校正過程,而且完成了動態的實驗分析。實驗不需要高精度、精密復雜的設備,試驗方法簡單有效且具有較高的實用價值。

[1]吳永亮,王田苗,梁建宏. 微小型無人機三軸磁強計現場誤差校正方法[J]. 航空學報,2011,32(2):330-336.

[2]張朝陽,肖昌漢. 海底布放磁傳感器的磁定位方法的模擬實驗研究[J]. 上海交通大學學報,2011,45(6):827-830.

[3]曹鵬,于紀言,王曉鳴,等. 基于地磁與衛星組合的高旋彈丸滾轉角高頻測量及系統誤差計算研究[J]. 兵工學報,2014,35(6):795-800.

[4]董海龍,楊冬梅,張素琴,等. 嘉峪關地震臺磁通門磁力儀記錄準確分析[J]. 地震地磁觀測與研究,2016,37(1):82-85.

[5]周鶴. 磁通門磁力儀采集系統研制[J]. 煤炭技術,2016,35(8):293-295.

[6]郞雪,杜愛民,李瓊,等. 基于PIC單片機的小型化低功耗磁通門磁強計的硬件設計和實現[J]. 地球物理學進展,2016,31(2):908-912.

[7]楊文,裴東興,張瑜,等. 彈丸動態參數測量的方法研究[J]. 傳感技術學報,2015,28(10):1467-1470.

[8]林春生,向前,龔沈光. 三軸磁強計正交誤差分析與校正[J]. 探測與控制學報,2005,27(2):10-12.

[9]龐鴻鋒,潘孟春,王偉,等. 基于高斯牛頓迭代算法的三軸磁強計校正[J]. 儀器儀表學報,2013,34(7):1506-1510.

[10]黃玉,郝燕玲. 基于FLANN的磁通門總場梯度計誤差校正[J]. 華中科技大學學報,2011,39(5):83-87.

[11]盧兆興,張金生,王仕成,等. 基于自適應遺傳算法的三軸磁強計誤差校正方法[J]. 傳感技術學報,2014,27(3):331-335.

[12]李勇,劉文怡,李杰,等. 基于橢球擬合的三軸磁傳感器誤差補償方法[J]. 傳感技術學報,2012,25(7):917-920.

[13]尹芳黎,楊雁瑩,王傳棟,等. 矩陣奇異值分解及其在高維數據處理中的應用[J]. 數學的實踐與認識,2011,41(15):171-177.

[14]趙學智,葉邦彥,陳統堅. 多分辨奇異值分解理論及其在信號處理和故障診斷中的應用[J]. 機械工程學報,2010,46(20):64-74.

楊福杏(1990-),男,黑龍江省同江市人,現為武漢大學在讀碩士研究生,目前主要從事數據采集處理,磁傳感器等方面的研究,953968487@qq.com;

曹大平(1961-),男,湖北省荊州市人,現為武漢大學物理科學與技術學院教授,博士生導師,科研活動圍繞磁傳感器與磁信號處理展開,主要包括:磁傳感器探頭磁芯磁場計算與分析;磁傳感器模擬電路設計與調試、電路噪聲與直流零點漂移抑制;矩陣計算、自適應濾波、數字濾波器設計與調試;探頭磁芯材料探索研究,dapingcao@whu.edu.cn。

Research on Calibration of Magnetic Modulus in the Motion State*

YANGFuxing,CAODaping*

(School of Physics and Technology,Wuhan University,Wuhan 430072,China)

Considering existing correction method for the triaxial flux-gate sensor it has a series of problems,such as complicated method,high precise equipment is needed,poor effect and so on. We make a careful analysis and calculation of measure error in triaxial flux-gate sensor caused by non-orthogonality,different sensitivity between three axes and zero-shift. And we put forward a three axis flux-gate sensor correction model,as well as proposed a calibration method based on singular value decomposition and the simplex method. The experimental results show that the method is simple and effective,and it can be widely applied to the correction of three-axis flux-gate sensor total static and dynamic modulus correction.

three-axis fluxgate magnetometers;calibration;SVD;simplex method

項目來源:國家自然科學基金面上項目(41274188)

2016-12-09 修改日期:2017-03-17

TH762.3

A

1004-1699(2017)08-1178-04

C:7220;7310L

10.3969/j.issn.1004-1699.2017.08.008