單個磁通門定位運動磁性目標范文

本站小編為你精心準備了單個磁通門定位運動磁性目標參考范文，愿這些范文能點燃您思維的火花，激發您的寫作靈感。歡迎深入閱讀并收藏。

《大連理工大學學報》2016年第三期

摘要：

運動磁性目標定位可歸結為一類非線性規劃問題的求解．針對其目標函數計算過程極為煩瑣，不利于磁性目標實時定位的特點，提出了一種混合算法：首先利用粒子群優化算法獲得初略解，再利用高斯－牛頓算法進行精確求解，在不知道目標初始位置的情況下，利用單個磁通門磁傳感器采集的磁場數據準確快速地實現了目標的定位．理論仿真和汽車實驗，驗證了算法的合理性和可行性，為工程中磁性目標的檢測與定位提供了新思路．

關鍵詞：

磁性目標定位；非線性規劃問題；粒子群優化算法；高斯－牛頓算法磁定位技術由于具有全天候、所需設備簡單、信號處理方便等特點，逐漸受到人們的重視．特別是隨著磁傳感器測量精度的不斷提高，探測磁性目標微弱的磁異常信號成為可能［1］．根據磁傳感器探測到的磁場信號，反演出目標的運動狀態和位置，可以實現運動磁性目標的定位．國外的磁定位技術起步較早，早在1975年，美國的海軍近海岸系統實驗室就開展了利用超導磁梯度計陣列進行目標定位的研究，他們利用磁場梯度張量數據提出磁偶極子定位跟蹤算法，根據靜態測站上的連續測量結果，實現了對運動的磁偶極子源的定位追蹤［2］．20年之后，又提出了在磁場梯度張量探測儀以已知速度運動的情況下對偶極子磁源的定位問題，對在已知磁場的三分量或梯度張量對場源變化率的條件下，結合磁場梯度張量可以唯一確定磁偶極子源參數進行了總結［3］．在國內，海軍工程大學龔沈光的研究團隊開展的相關工作較多，他們［4－9］對磁性目標建模和基于艦船等磁性物體定位方法及參數估計等定位問題進行了詳細研究，認為用單個運動標量磁強計不能完全確定目標位置參數，提出了先用改進的進化策略得到初值，然后用單純形法獲得精確值的一種通用磁定位方法，并稍加修改解決了矢量磁強計的磁定位問題．

隨后林春生等［10］使用均勻磁化旋轉橢球體模型對水中磁性運動目標進行實時建模，使用統計檢驗的方法，實現了對目標磁場信號的模型化檢測，同時還可以估計出目標運動參數和磁性狀態．李華等［11］使用單個三軸磁探頭的在線測量數據，利用遺傳算法良好的全局搜索能力對磁性目標進行定位計算，以解的穩定性判斷代替定位解正確與否的在線判斷，完整實現了一種在線磁定位算法．唐莉莉等［12］提出了一種新的運動磁目標定位算法，將求解目標方位角的問題轉化為求解目標磁場的方向角和兩者的偏差角的問題，利用單個靜止的平面四軸磁通門傳感器實現目標定位．何文輝［13］提出了利用三軸磁阻傳感器作為檢測元件，以減少傳感器陣列的數量，并通過改進的牛頓－拉斐森方程組全局求解算法，實現了對膠囊的位置坐標準確而快速的求解．本文提出一種混合算法，在不需要知道目標初始位置的情況下，利用單探頭磁通門，準確快速地求解磁性目標定位中的非線性規劃問題．

1單個磁通門定位運動磁性目標的基本原理

單個磁通門定位的基本方案：通過單個磁傳感器測量運動磁性目標在n個不同測點上的磁場強度構建磁場模型，利用定位算法反演出磁性目標在不同測點處相對磁傳感器的位置，即利用前n－1個測點的磁場數據，定位出第n個測點的位置．

1．1磁偶極子模型對于空間中某固定磁偶極子模型（磁矩M0不變），如果知道測點相對磁偶極子的位置（即F0已知），則該測點處的磁場強度H0就可以由式（6）求出；反之，如果知道了某一測點處的磁場強度，也可以根據上式反演出測點相對磁偶極子的位置，從而實現磁性目標的定位．

1．2測點間的坐標關系將目標從坐標系的第m點平移到第m＋1點處，則兩者之間的相對位置關系如圖2所示．其中φ為向量r和Z軸正方向的夾角；θ為向量r在XOY面上的投影與X軸正方向的夾角．

1．3磁性目標的定位問題磁性目標的定位問題實際上就是求解下面非線性無約束的最優化問題。

2混合算法下的磁性目標定位仿真算法

如圖3所示，磁性目標沿X方向運動范圍為－200m到195m，Y方向距離磁場測量點記為－30m，Z方向距離磁場測量點記為0m．磁通門x方向與坐標系X軸之間的夾角θ＝－45°，目標在坐標系XOY平面內運動，即φ＝90°．目標運動速度為10m／s，磁場數據采樣間隔為0．5s，整個過程磁傳感器采集的數據點為80個．高斯－牛頓算法（G－N算法）是一種常見的求解非線性最優化問題的方法，但是該方法對初值有一定的依賴性，需要利用函數的導數信息．當初值取得不合適時仍然會導致算法無法得到最優解．為了說明這一情況，分別取不同的初值使用G－N算法進行求解，隨機運行50次，選取位置參數（－100m，－30m，0m，－45°，90°）的仿真數據進行比較，結果如表1所示．從表1可以看出，初值距離真值一定范圍內，G－N算法能夠找到很好的最優解；但是當初值選取不合適時，該算法得不到最優解．大多數情況下，并不清楚目標的初始狀態．針對這個問題，考慮使用粒子群優化算法（PSO算法）．鑒于PSO算法耗時較長［14－15］，本文提出了一種混合算法：先使用PSO算法得到一個粗略解，為了節約計算量，PSO算法的進化代數可以適當減少，再以粗略解作為初值使用G－N算法進行精確求解，這樣將PSO算法和G－N算法混合使用，充分利用各自的優勢，在保證算法精度的前提下，盡量減小算法耗時．為了驗證該混合算法，同樣選取位置點（－100m，－30m，0m，－45°，90°）的數據進行比較得到的仿真數據分別為x＝－99．99m，y＝－30．00m，z＝－0．001m，θ＝－44．96°，φ＝89．5°，定位效果較好．

3實驗驗證

為了方便模擬目標的運動狀態，磁性目標選擇的是一輛汽車．為了減小周圍環境的干擾，實驗地點選擇在西安市郊區的一段很少有車輛通過的寬闊馬路．第一次汽車通過的情況如圖4（a）所示，X方向指向汽車行駛的方向，Z方向垂直地面向下，Y方向與X方向和Z方向符合右手坐標系．汽車長為4．55m、寬為1．80m、高為1．88m，保持約為30km／h的行駛速度正好通過磁傳感器時橫向距離約為3．8m，汽車中心距離地面的高度約為0．8m，磁通門x方向與坐標系X軸重合，即θ＝0°，φ＝90°．汽車第二次通過時的情況如圖4（b）所示，汽車仍以第一次的速度行駛，正好通過磁傳感器時的橫向距離約為4．2m，汽車中心距離地面的高度約為0．8m，磁通門x方向與坐標系X軸之間的夾角θ＝－30°，且φ＝90°．汽車兩次通過磁通門傳感器，每次截取剛好通過前后共4s的數據，數據采樣率為20Hz．由于存在背景磁場的干擾，使用巴特沃斯低通濾波器對數據進行了低通濾波處理，濾波器的通帶截止頻率為1Hz，通帶衰減為1dB，阻帶截止頻率為10Hz，阻帶衰減為60dB，處理后的數據如圖5和6所示．使用PSO和G－N的混合算法進行定位計算，目標磁場數據總長度為80，定位效果如圖7和8所示．由圖7可見，汽車距離磁傳感器較遠時，磁場信號弱，此時定位結果較差．隨著汽車距離磁傳感器越來越近，磁場信號逐漸增強，有效的磁場數據越來越多，定位結果越來越好．定位結果明顯改善的時間點出現在汽車正橫通過磁傳感器，此時磁場信號最強．可以看出：X方向在時間軸52～53s上大致是一條斜線，斜率約為8．385m／s，Y方向在時間軸53s的距離約為－3．774m，Z方向在時間軸53s的距離約為－0．673m，與X方向的夾角在時間軸53s時約為5．651°，與豎直方向的夾角在時間軸53s時約為89．89°．實測數據定位效果與圖4（a）所示的汽車實際行駛狀態吻合較好，表明了混合算法下的磁性目標定位方法的可行性．從圖8中同樣可看出：X方向在時間軸379～380s上約是一條斜線，斜率約為8．349m／s，Y方向在時間軸380s的距離約為－4．181m，Z方向在時間軸380s的距離約為－0．9494m，與X方向的夾角在時間軸380s時約為－27．43°，與豎直方向的夾角在時間軸380s時約為90°．實測數據定位效果與圖4（b）所示的汽車實際行駛狀態也吻合較好，驗證了混合算法下的磁性目標定位方法的實用性．

4結語

本文針對運動的磁性目標初始位置未知的情況，提出了一種先使用PSO算法進行粗略求解，再利用G－N算法精細求解的混合算法．理論分析和汽車實驗驗證了該算法定位結果準確可靠，既擺脫了G－N算法在磁性目標定位求解中對初始值的依賴，也克服了PSO算法求解進化迭代耗時過長的問題．同時單個磁通門磁傳感器價格低廉、易于安裝的特點，使其在工程實際中有著很好的應用前景.

作者：高翔 嚴勝剛 李斌 單位：西北工業大學 航海學院