基于聲陣列定位系統的時(shí)差信息提取方法的研究

作者  邵云峰 韓焱 中北大學(xué) 信息探測與處理技術(shù)研究所(山西 太原 030051)

*基金項目 : 國家自然科學(xué)基金(編號：60772102);國防重點(diǎn)實(shí)驗室基金(編號：9140c1204040908);國防重點(diǎn)實(shí)驗室基金(編號：9140c12040051010)

　　邵云峰(1990-)，男，碩士，研究方向：信號處理與定位算法，無(wú)線(xiàn)傳感器網(wǎng)絡(luò )微震信息探測;韓焱，男，教授，研究方向：信號與信息處理、精密儀器及機械、信息處理與重建等。

摘要：基于地震動(dòng)信號的可穿戴式被動(dòng)聲探測系統具備體積小巧、方便攜帶等優(yōu)點(diǎn)。該研究采用地面走動(dòng)信號的目標定位技術(shù)設計了鞋底聲傳感器探測陣列及其數據采集系統。利用數據采集系統完成了陣列信號的獲取、分析和處理，介紹了目標定位的方法及時(shí)延估計方法，并設計了數據的預處理與時(shí)延估計的Matlab程序實(shí)現，實(shí)現了陣列延時(shí)的檢測，針對時(shí)延估計不準確的問(wèn)題，從系統、算法精度以及頻域分析等方面做了分析，提出了改進(jìn)的方向。

引言

　　在科技飛速進(jìn)步，信息即是財富的時(shí)代，信息的地位日趨重要。軍事專(zhuān)家們預言：21世紀的戰爭將是一場(chǎng)別開(kāi)生面的信息戰?；?a class="contentlabel" href="http://dyxdggzs.com/news/listbylabel/label/地震動(dòng)信號">地震動(dòng)信號的地面被動(dòng)偵察^[1]在復雜的地形條件下，甚至在嚴密偽裝的情況下仍能充分發(fā)揮其作用，這正是光學(xué)、無(wú)線(xiàn)電等現代監視技術(shù)的盲區，尤其是被動(dòng)聲探測技術(shù)^[2]。

　　但在特殊情況下，如大風(fēng)等惡劣天氣或在植被非常茂盛的地方，聲的傳播特性就會(huì )受到很大影響，致使聲探測系統無(wú)法正常工作。因此，研究基于地面震動(dòng)信號^[3]的偵察技術(shù)具有非常重要的價(jià)值。

　　地震動(dòng)信號較其他的被動(dòng)偵察信號不同，其信噪比較差，受環(huán)境影響較大，但由于地震動(dòng)探測具有可全天候值守、抗電磁干擾能力強、功耗低、隱蔽性好、成本低等優(yōu)點(diǎn)，作為一種重要的被動(dòng)偵察手段，其地位的不可小覷。而基于地震動(dòng)信號的可穿戴式被動(dòng)聲探測系統因其體積小巧、方便攜帶，并且可使穿戴者根據情況迅速做出反應等優(yōu)點(diǎn)，正日趨受到廣泛的關(guān)注。

1 定位原理與傳感器布設

1.1 定位原理

　　聲源位置與傳感器的平面分布圖如圖1，聲源定位的基本原理是：傳感器陣列中每個(gè)傳感器距聲源的距離是不一樣的，因此，聲源發(fā)出聲音后，傳感器陣列中每個(gè)傳感器采集到的信號會(huì )存在時(shí)間差，此時(shí)間差乘以聲音在介質(zhì)中的傳播速度可求出距離差，5個(gè)傳感器中每3個(gè)傳感器可以確定出2個(gè)距離差，據此我們可以列方程算出聲源的坐標，實(shí)現聲源的定位。

1.2 傳感器布設

　　為了實(shí)現全方位信號探測，在鞋底上布設5個(gè)傳感器，前腳掌布設4個(gè)傳感器，呈菱形分布(前腳掌每個(gè)傳感器之間相距5cm)，后腳掌布設一個(gè)傳感器，距前腳掌最后面的那個(gè)傳感器距離為12cm，并對傳感器編號，如圖2所示。

2 目標定位與時(shí)延估計方法

　　目前，目標定位方法主要分為有源定位和無(wú)源定位^[4]。有源定位就是利用雷達、激光等有源設備來(lái)對目標進(jìn)行定位，但其存在容易暴露、能耗較高的弊端。無(wú)源定位是通過(guò)對目標自身產(chǎn)生的輻射源信號的接收和處理，提取目標的距離、方位、航跡等信息。無(wú)源定位克服了有源定位的缺點(diǎn)。本設計的方向是采用通過(guò)多地震動(dòng)傳感器實(shí)現目標定位的方法，屬于無(wú)源定位方法。

　　目前，無(wú)源定位算法中目標定位[5-8]方法主要有基于到達時(shí)間(TOA)、到達時(shí)間差(TDOA)和接收信號強度或能量(RSSI)法。TOA技術(shù)算法簡(jiǎn)單，但是要求移動(dòng)目標和節點(diǎn)始終精確同步，易受多徑傳播和噪聲的干擾。RSSI技術(shù)較簡(jiǎn)單，近距離精度高，受環(huán)境影響較大，室外遠距離定位精度較差。

　　時(shí)延估計^[9-10]是數字信號處理中一個(gè)非?；钴S的研究領(lǐng)域。在雷達、語(yǔ)音信號處理、地球勘探等領(lǐng)域都有較廣泛的應用。

　　直接互相關(guān)算法計算簡(jiǎn)單，要求信號、噪聲互不相關(guān);廣義互相關(guān)相當于對信號進(jìn)行了白化處理，抑制噪聲當信噪比較低時(shí)誤差較大;互雙譜算法對背景噪聲不敏感，可以抑制空間噪聲，不適合非高斯造聲的情形，且信號序列要足夠長(cháng)，計算量和估計方差較大。

　　直接互相關(guān)時(shí)延估計是通過(guò)對兩信號直接做互相關(guān)處理，然后進(jìn)行峰值檢測來(lái)估計時(shí)延。廣義互相關(guān)時(shí)延估計與基本互相關(guān)類(lèi)似，但其在小信噪比時(shí)，按平方加權;在大信噪比時(shí)，按信噪比加權。相當于對信號進(jìn)行了白化處理，一定程度上抑制了噪聲。

3 信號處理

3.1 信號預處理

　　采集獲取的傳感器信號混有大量的噪聲，因此需對其進(jìn)行預處理。本設計基于中值濾波器優(yōu)良的抑制噪聲的特點(diǎn)，首先對采集到的信號進(jìn)行中值濾波處理，以便濾除毛刺。

　　中值濾波之后，為了便于數據處理，需要截取出采集到信號的有用部分，通過(guò)觀(guān)察采集到的信號可知，有用信號的幅值遠大于噪聲的幅值。由正態(tài)分布的性質(zhì)可知，如果X服從正態(tài)分布，則：

　　其中，指的是正態(tài)分布的均值，指的是正態(tài)分布的標準差，因此我們取數據的均值加3倍的標準差的值作為閾值，據此截取出信號的有用部分。

3.2 直接相關(guān)法時(shí)延估計

　　表1是用直接相關(guān)法算出的時(shí)延。接下來(lái)我們根據坐標位置，算出理論上每?jì)蓚€(gè)傳感器到聲源的距離差，用此距離差除以估計出的時(shí)延，求出速度，圖4是用估計出的時(shí)延算出的理論速度，將其與聲速進(jìn)行比較，發(fā)現我們用估計出的時(shí)延算出的速度和聲速相差太大，造成這種結果的原因可能是因為信號之間相關(guān)性較強，致使估計出的時(shí)延不準確。因此，我們采用廣義相關(guān)算法估計時(shí)延。

3.3 用廣義互相關(guān)算法估計時(shí)延

　　在這里使用PHAT加權進(jìn)行時(shí)延估計。表2是我們用PHAT加權估計出的時(shí)延值，圖5是用根據估計出的時(shí)延算出的理論速度，將其與聲速進(jìn)行比較，發(fā)現算出的10個(gè)速度值中有5個(gè)和聲速在一個(gè)數量級上。觀(guān)察圖形可判斷出，使用算法估計出的時(shí)延值要比直接互相關(guān)估計出的時(shí)延值準確，10個(gè)速度值中有5個(gè)與聲速不在一個(gè)數量級上，造成這種結果的原因可能是在測量坐標點(diǎn)位置時(shí)人為造成的誤差，也有可能是由于多徑效應造成的。

4 結果分析

　　對截取后的數據進(jìn)行頻譜分析，其幅頻特性曲線(xiàn)如圖6所示，從圖中可以看出信號頻率成分較為復雜，造成這種結果的原因可能是多徑效應^[11]的影響。從圖中可以看出信號頻率主要集中在兩個(gè)頻段：970Hz~1600Hz，3330Hz~3880Hz，假設需要的信號是一個(gè)單頻成分，那么就需要對信號在整個(gè)頻段上進(jìn)行窄帶濾波，返回時(shí)域看其波形并進(jìn)行分頻時(shí)延估計，從中找出我們需要的信號，并進(jìn)行后續分析。

　　用信號發(fā)生器對音箱輸入1kHz，20V峰峰值的正弦信號，將音箱放到點(diǎn)(-5m,0)坐標點(diǎn)上，采樣頻率為225kHz，圖7是采集到其中一個(gè)通道信號的時(shí)域波形，圖8是其幅頻特性曲線(xiàn)圖。結合上述三張圖可以看出采集到的信號并不是一個(gè)頻率單一的標準正弦波，圖9和圖10是對其幅頻特性曲線(xiàn)中強度(50Hz左右，1950Hz左右)最大的兩個(gè)干擾的放大圖，這兩個(gè)干擾和1kHz信號的強度相差一個(gè)數量級，不能被忽略，造成這種結果的原因可能是由于多徑效應引起的，也有可能是環(huán)境中的交流干擾造成的，因此需要進(jìn)一步對數據進(jìn)行處理。

　　通過(guò)對比圖6和圖8可發(fā)現，我們采集到的信號受到的干擾要比采集一個(gè)單頻聲音信號受到的干擾多，從圖6我們可以看出強度最大的兩個(gè)峰值分別是45dB和50dB，圖9和圖10分別是對兩個(gè)峰值取下降15dB部分的放大圖，從圖中可以看出，第一個(gè)峰值處衰減15dB后，頻段集中在1250Hz~1450Hz，第二個(gè)峰值處衰減15dB后頻段集中在3250Hz~3650Hz，編一個(gè)FIR模擬帶通濾波器對以上兩個(gè)頻段進(jìn)行濾波，圖11和圖12分別是對兩個(gè)頻段作帶通濾波后的時(shí)域波形。和圖3進(jìn)行對比后，發(fā)現第二個(gè)峰值衰減15dB的部分比第一個(gè)峰值衰減15dB的部分與需要的信號的相似度高，說(shuō)明第二個(gè)峰值衰減15dB的部分中所含的信號的有用部分比第一個(gè)峰值衰減15dB的部分包含的多。

5 結論

　　本研究是在被動(dòng)聲定位技術(shù)被廣泛應用到軍事及民用中，成為定位技術(shù)領(lǐng)域的研究熱點(diǎn)背景下展開(kāi)的，主要是針對被動(dòng)聲定位技術(shù)中的目標定位方法和時(shí)延估計算法進(jìn)行了研究。完成了信號的預處理的相關(guān)程序設計，對信號進(jìn)行了中值濾波，截取出了信號的有用部分，方便信號的后續處理;完成了時(shí)延估計算法的程序實(shí)現，對結果進(jìn)行了分析，提出了改進(jìn)的方向。還有以下幾方面需要改進(jìn)：對時(shí)延估計算法的進(jìn)一步優(yōu)化進(jìn)行程序實(shí)現;針對算法計算不準確的問(wèn)題，應對信號分頻段估計時(shí)延，從中找出需要的信號，從而實(shí)現較準確地時(shí)延估計。

　　參考文獻：

　　[1]陳立學(xué),魚(yú)云岐.地面偵察與精確打擊[J]. 兵工學(xué)報,2010,S2:79-82.

　　[2]馬春庭,高萍.地面戰場(chǎng)偵察傳感器系統的發(fā)展狀況及技術(shù)分析[J].軍械工程學(xué)院學(xué)報,2003(4).

　　[3]聶偉榮.多傳感器探測與控制網(wǎng)絡(luò )技術(shù)——地面運動(dòng)目標震動(dòng)信號探測與識別[D].南京理工大學(xué),2002.

　　[4]范志平,鄧平,劉林.蜂窩網(wǎng)無(wú)線(xiàn)定位[M].北京:電子工業(yè)出版社,2002

　　[5]雷文英,陳伯孝,楊明磊,等.基于TOA和TDOA的三維無(wú)源目標定位方法[J].系統工程與電子技術(shù),2014,05:816-823.

　　[6]張馳.基于無(wú)線(xiàn)傳感器網(wǎng)絡(luò )的無(wú)源目標定位跟蹤仿真研究[D].北京郵電大學(xué),2013.

　　[7]王小雪.基于無(wú)線(xiàn)傳感器網(wǎng)絡(luò )的無(wú)源被動(dòng)式目標定位研究[D].浙江工業(yè)大學(xué),2013.

　　[8]單振興.目標無(wú)源定位與跟蹤算法研究[D].國防科學(xué)技術(shù)大學(xué),2005.

　　[9]李雪梅,陶然,王越.時(shí)延估計技術(shù)研究[J]. 雷達科學(xué)與技術(shù),2010,(04):362-371.

　　[10]行鴻彥,唐娟.時(shí)延估計方法的分析[J]. 聲學(xué)技術(shù),2008,01:110-114.

　　[11]崔瑋瑋,曹志剛,魏建強.聲源定位中的時(shí)延估計技術(shù)[J].數據采集與處理,2007,22(1):90-95.

　　本文來(lái)源于《電子產(chǎn)品世界》2017年第8期第29頁(yè)，歡迎您寫(xiě)論文時(shí)引用，并注明出處。