基于相關性算法的超聲波回波定位

胡 暉 許浩峰 包偉華

(上海自動化儀表股份有限公司，上海 200072)

基于相關性算法的超聲波回波定位

胡 暉 許浩峰 包偉華

(上海自動化儀表股份有限公司，上海 200072)

隨著人們對產品的功耗要求不斷提升，復雜的硬件比較器方式越來越不能適應工業領域發展，所以設計了一種基于A/D采樣軟件算法的超聲波物位測量系統。軟件采用相關性算法來描述兩組信號的相似程度，并加以顯著性差異檢驗，來檢驗兩組數據是否顯著相關，從而準確地進行回波定位。現場超聲波物位儀表性能測試表明，該軟件算法能夠快速、穩定、可靠地找到回波，有一定的發展價值。

A/D采樣 低功耗 相關性 物位儀表 超聲波回波 回波定位

0 引言

在采用A/D采樣方式的測距系統中，軟件選用何種合適的算法，直接影響了測距系統的精度及實時性。根據超聲波的特性，在傳播途徑中波形會衰減，導致回波波形的幅值小于發射波形的幅值，但是從根本上講，回波是由發射波產生的。分別將發射波與回波或者2個周期的回波進行比較，兩者必然有一定的聯系。相關概念是為研究隨機信號的統計特性而引入的，相關性可以用來描述兩組信號的相似程度。如果2個周期內的回波均為真實回波，則兩者的相似程度肯定很高，具體相似程度由相關系數表示，之后再進行顯著性分析，對由相關系數得出的結論作出肯定或否定的決策[1]。

1 相關系數及顯著性算法

相關分析發端于1900年代統計學開創時期，目前仍是統計信號處理領域的研究熱點。所謂相關，是表征兩個隨機變量或兩路信號直接的統計關系強弱的指標。如果一個隨機變量隨著另外一個隨機變量的增大(減小)而增大(減小)，則這兩個隨機變量滿足正相關關系；反之，如果一個隨機變量隨著另一個隨機變量的增大(減小)而減小(增大)，則這兩個隨機變量滿足負相關關系[2,6]。

相關系數是用以反映變量之間相關關系密切程度的統計指標。相關系數是按積差方法計算，同樣以兩變量與各自平均值的離差為基礎，通過兩個離差相乘來反映兩變量之間相關程度。公式如下：

(1)

式中：x、y為兩組序列，并且兩組序列長度相等且為n，則兩者的相關系數為r。

由式(1)計算得出r值。一般來說，r取絕對值后，0～0.3為低度相關，0.3～0.8為中度相關，0.8～1.0為高度相關。

雙側：Z0.05/2=1.96,Z0.01/2=2.58

單側：Z0.05=1.645,Z0.01=2.33

Z檢驗一般是用于大樣本(即樣本容量大于30)平均值差異性檢驗的方法。它是用標準正態分布的理論來推斷差異發生的概率，從而比較兩個平均數的差異是否顯著。

統計量Z值的計算公式為：

就智能化而言，高斯（中國）為自己插上了一雙幫助自己可以飛得更高的翅膀，一是公司推出的匯印e家印機管理平臺，其以印刷設備為載體，從各類設備采集印刷狀態實時數據，應用物聯網、云計算、大數據分析等技術，為印刷企業構建全方面、多層次、全覆蓋的印刷管理體系，并賦予印刷設備感知、思考及服務的能力；二是智能化的聯線產品系統，其可以將印前、印刷和印后進行無縫連接。劉忠榮總工程師告訴我們，前不久，高斯（中國）就與青海日報社印刷廠就“匯印e家”智能化數據采集及控制進行了良好合作，并將把這樣的樣板效能推介到更多的印刷廠中。

(2)

式中：X1、X2為兩組序列，并且兩組序列長度分別為n1和n2，則統計量Z值為Z。

由式(2)可計算得出Z值。當|Z|值小于1.96時，認為這兩組數據無顯著差異;|Z|值大于1.96時，認為這兩組數據有著顯著的差異。

2 超聲波信號

超聲波在介質中傳播，遇到障礙物或物料等則會反射形成回波。本算法由于同時需要發射波的波形特征及回波的波形特征，所以CPU必須在發射波發出之前就開始采樣其電壓值。如圖1所示，這是一個測量周期內采集的數據， 包括了發射信號、回波信號及一些干擾信號，軟件需要實現的工作就是找到真實回波，排除干擾信號[3-4]。

圖1 A/D采樣數據

2.1 發射信號與回波信號

采用相關性算法來找到真實回波，發射信號特征非常重要，因為它是一個參考標準，所以待測回波或者虛假回波(干擾信號)都需要與之進行比較分析，計算出相關系數，并成為判斷是否為真實回波的一個條件。

與回波的位置不可控相比而言，可以非常準確地判斷出發射信號的位置及持續時間，這和超聲波探頭的種類及其盲區有一定的關系。從圖1可知，在某段時間內采樣數據，通過對其圖形分析，可以簡單地找到其波峰頂點，也就是超聲波發射聲強[5]。

相隔一個周期的兩個回波信號，因為是同一個發射信號產生的，兩者的相似程度必然是極高的，因此利用相關性算法，可以得出兩者的相關系數，并進行顯著性分析。而干擾信號的來源廣泛，可能是電路上的干擾，也可能是收到了其他雜波，而不知道它是干擾信號，可能將一些大的干擾信號當做回波與上一周期回波進行處理，并計算得出其相關系數。但是干擾信號的特征相似程度必然遠遠不如真實回波，所以取相關性系數最大的作為判斷是否為真實回波的一個條件。

2.3 真實回波

通過對相關系數及顯著性分析的計算，可以得出判斷為真實回波的條件：待測回波信號與上一周期真實回波信號的顯著性Z值小于1.96，認為這兩個信號無顯著差異，即是由同一個發射源產生。

3 回波快速定位

基于相關系數及顯著性算法，可以較為準確地得出回波位置。但是在數據量較大的情況下，如果采用逐點計算的方式，一組數據的處理計算可能需要套用公式幾百次，這是非常耗時的。所以需要大致地估摸出波峰的位置，并采用窗口方式及自學習曲線來快速地進行定位[7]。

3.1 初始方式

在系統開機后的第一次搜索回波或者是在上一周期丟失回波的情況下，系統往往缺少上一周期的數據來幫助搜索，因此找到第一個正確的回波很關鍵。

3.2 窗口方式

只要確定上一周期的回波位置，就可以在其周圍一定窗口范圍內搜索回波，從而大量地減少搜索范圍。如圖2所示，在預定長度的窗口范圍內找到了合適的回波，則可以直接得出結果；如果搜索失敗，則直接進行全局方式搜索。

圖2 窗口方式

3.3 全局方式

如果在窗口方式內沒有搜索到回波，那就必須擴大搜索范圍，將范圍擴大至全局。雖然全局方式耗時要比窗口方式慢，但是首先必須保證搜索到真實回波，然后才考慮搜索速度。

3.4 自學習曲線

經過多次的測量周期之后，可以很清楚地了解到整個采樣數據組的情況。例如，有些干擾可能是電路上的固定或者周期性干擾。有些干擾可能是外來的不確定干擾，通過自學習方式，就可以排除一些固定或者周期性干擾，如圖3所示，所有采樣值低于自學習曲線的信號都將視作無效或干擾信號。

圖3 自學習曲線

3.5 實驗數據

在8 MHz晶振的國產CMC芯片開發平臺下，進行了超聲波物位計產品性能測試，實驗數據如表1所示。在距離不變的情況下連續測距，計算耗時保證在300 ms以內，并且重復性及穩定性都相當不錯。

表1 實驗數據

4 軟件流程

主控芯片上電進行初始化工作。由于是低功耗產品，因此需要對激勵電路進行充電，充電完畢后開啟A/D采樣，并將儲存的電能用作激勵信號激勵超聲波探頭，隨后一直采樣到預定時間，開始進行數據分析。

為節省計算時間，主控芯片根據上一周期的計算結果來處理這一周期的計算數據。如果上一周期已經找到回波，則這一周期先以開窗的方式尋找回波；如果未找到，則以全局方式搜索回波。如果上一周期未找到回波(例如開機情況下)，則以初始方式尋找回波。隨后將計算結果通過顯示屏或者上位機進行數據顯示。軟件流程圖如圖4所示。

圖4 軟件流程圖

5 結束語

本文介紹了一種基于相關性的超聲波回波定位的方法，并采用8 MHz晶振的國產CMC芯片實現了一種超聲波物位計測量的方案。根據超聲波物位計的特點，研究了相關性算法、顯著性算法等關鍵算法，還應用窗口方式和自學習曲線來完善軟件邏輯，給出了軟件方案。經過現有產品實驗證明，算法性能已完全達到了設計要求，為未來研究打下了基礎。

[1] 劉成瑞,傅惠民.相關系數平穩序列自適應算法[J].機械強度，2003,25(6)：661-665.

[2] 趙小川,羅慶生,韓寶玲.基于偽隨機序列自相關性的新型超聲波測距系統[J].壓電與聲光，2009,31(6)：856-861.

[3] 陶建平,尹文慶,柳軍.基于DSP和單片機的超聲波測距系統[J].科學技術與工程，2010,10(6)：762-765.

[4] 蘇煒,龔壁建,潘笑.超聲波測距誤差分析[J].傳感器技術，2004,23(6)：8-11.

[5] 童峰,許水源,許天增.一種高精度超聲波測距處理方法[J].廈門大學學報:自然科學版，1998,37(4)：507-512.

[6] 程曉暢,王躍科，蘇紹，等.基于相關函數包絡峰細化的高精度超聲測距法[J].測試技術學報,2006，20(4)：320-324.

[7] 付光懷,楊惠中.基于雙正交樣條小波的超聲回波信號檢測[J].儀器儀表學報,2008,29(4)：309-310.

[8] 李戈,孟祥杰,王曉華,等.國內超聲波測距研究應用現狀[J].測繪科學，2011,36(4)：60-62.

[9] 鄒曉紅.小波分析在超聲波回波測量中的應用[J].傳感器技術,2005,24(7):78-82.

[10]陳先中，侯慶文，葉琳，等.超聲測距系統的高精度中心橢圓算法[J].北京科技大學學報，2007，29(11):1154-1157.

Ultrasonic Echo Location Based on Correlation Algorithm

As people continuously to enhance the requirement for power consumption of the products, the method of complex comparator becomes increasingly unsuited for development of industrial areas. Thus the ultrasonic material level measurement system based on A/D sampling software algorithm is designed. The system adopts the correlation algorithm to describe the degree of similarity of two groups of signals, and test the significant difference, for judging if these two groups of signals are significantly interrelated, thus precise echo location is conducted. The on-site performance test of ultrasonic material level measuring instrument show that the software algorithm is able to find the echo rapidly and stably; it possesses certain development value.

A/D sampling Low power consumption Correlation Material level instrument Ultrasonic echo Echo location

國家高技術研究發展計劃(863計劃)基金資助項目(編號：2012AA041704)。

胡暉(1989-)，男，2012年畢業于華東理工大學測控技術與儀器儀表專業，獲學士學位，助理工程師；主要從事嵌入式開發方面的研究。

TH82;TB559

A

10.16086/j.cnki.issn1000-0380.201510024

修改稿收到日期：2015-02-11。