0 引言
在超聲波流量計中,，時延是影響測量精度的重要因素之一。目前,，超聲波流量計中常用的時延估計方法可以分為閾值法、相互關(guān)函數(shù)法,、頻域相位檢測法和時頻分析辦法,。閾值法簡略且輕易操作、實時性好,，但該方法很難正確地斷定回波峰值的呈現(xiàn)時刻,，故難以對實際的時間差進行準確估計。頻域相位法和時頻剖析方式都存在計算量大,、存儲容量請求高的毛病,，在超聲波流量丈量體系中并不適用。
相互關(guān)函數(shù)法是超聲波流量計中最常用的算法,。因為超聲波信號的振蕩特征,，其相關(guān)函數(shù)將在極值鄰近振蕩，從而給搜索其相關(guān)函數(shù)的準確極值帶來艱苦,，所以有必要將搜尋相關(guān)函數(shù)的極值轉(zhuǎn)化為搜索相關(guān)函數(shù)包絡(luò)的極值,。
在樹立超聲波流量計模型的基本上，分析了超聲回波信號解調(diào)的慣例算法,，并提出了基于fft的提取相關(guān)函數(shù)包絡(luò)的算法,；通過仿真分析了其計算龐雜度，同時也證實了本文提出的算法是一種計算量小,、可能精確搜尋相關(guān)函數(shù)峰值點的方法,。
1 時延法估計模型
1.1 超聲波流量計測量模型
時差法超聲波流量計測量模型管壁的高低兩端分辨有一個收發(fā)式的超聲換能器，在把持系統(tǒng)的作用下,，它們輪換發(fā)射接收固定頻率的超聲波信號,。因為管內(nèi)流體的作用，使得接收到的超聲信號的時間有一個時間差Δt,。
設(shè)超聲波在外界即空氣中的速度為c,，管壁中的速度為c1，靜止液體中速度為c2,；換能器的發(fā)射角為θ=45°,，第一次折射角為θ1，第二次折射角為θ2,；管的內(nèi)直徑為d，通過折射定理可以得到θ,、θ1和θ2之間的關(guān)聯(lián),。若超聲信號從上端發(fā)射下端接收所需要的時光為t1，超聲信號從下端發(fā)射上端吸收所須要的時間為t2,，則有：
式中：τ1,、τ2為超聲波在管壁傳輸?shù)臅r間以及硬件電路的延時；l為超聲波在管內(nèi)中聲道長度；v為流體的流動速度,。
這里的τ1,、τ2近似相等，則順流,、逆流的時間差為：
因為：
可以得到時光差為：
由式（6）可得,，液體速率為：
因為c2與溫度有關(guān)，并非一個常數(shù),，應(yīng)進行實時修改,，從而得到：
由式（8。流量計精度高,穩(wěn)定性好,品種齊全,價格實惠,是您的最佳選擇,。）可以得到修正速率為：
最后,，可以得到管內(nèi)液體的單位時間流量為：
1.2 超聲波時延估計模型
超聲波信號具備以下3個特色。
1相近性,，不同位置的回波信號波形相近,。
2相關(guān)性，對統(tǒng)一安裝跟統(tǒng)一對象,，超聲波的吸收信號隨發(fā)射間隔的轉(zhuǎn)變只有強弱的變化,，而波形變更不大，即發(fā)射接受信號之間是親密相關(guān)的,。
3窄帶性,，由于超聲波信號是以探頭諧振頻率為主頻率的衰減振蕩信號，所以信號的頻率主要散布在以換能器的諧振頻率為中央的一個較窄的頻域上,。
針對以上3個特色,，可以樹立一個超聲波的時延信號的模型，其表達式如下：
式中：s（t）為超聲換能器的發(fā)射信號,；α為,。氣體超聲波流量計價格實惠,超大量程,您的優(yōu)質(zhì)選擇。衰減因子,；n（t）為時間延遲,，n（t）為與s（t）不相關(guān)的零均值高斯白噪聲；t為回波信號觀測時間,。
由前文所述,，s（t）能夠建模存在慢起伏包絡(luò)的正弦調(diào)制信號：
式中：a（t）為發(fā)射信號包絡(luò)；f為換能器諧振頻率,；φ為初相,。
已知管的直徑為d，上真?zhèn)�超聲換能器接受到的信號為
下真?zhèn)�超聲換能器接收到的信號為
故時延估計問題可等效為：
將上式進行離散化：
式中：t為采樣時間距離,；
在已知觀測信號r1（nt）,、r2（nt）的情況下,，如何求出
估量值τ^就是接下來須要研討的時延估量問題。求出τ^之后,，則待測流體的流量可由以下公式盤算得到：
式中：c0為超聲波在流體中的速率修改值,。
1.3 相關(guān)函數(shù)峰值算法的提取
本文只斟酌τ2＞τ1，也就是τ＞0的情形,。r1（nt）,、r2（nt）的相干函數(shù)抒發(fā)式為：
式中：css（mt）為發(fā)射信號s（t）的自相關(guān)函數(shù)。
從超聲波信號以及超聲波信號的包絡(luò)特點,。進口超聲波流量計專業(yè)生產(chǎn)廠家,多年手持式超聲波流量計生產(chǎn)經(jīng)驗,。可以看出,，css（mt）有高頻振蕩的特性,，而發(fā)射信號包絡(luò)的自相關(guān)函數(shù)caa（mt）擁有遲緩衰減的特征。所以,，很難搜索c(mt)的精確峰值地位,。
1.4 基于fft的相關(guān)函數(shù)包絡(luò)提取的算法
針對以上這種情形，本文引入了一種基于fft疾速傅里葉變換的直接提取相關(guān)函數(shù)包絡(luò)的算法,。
首先,，通過fft得到一個方向發(fā)射信號r1(n)和接收信號r2(n)的頻譜r1(n)、r2(n),，計算互相關(guān)函數(shù)c(m)的頻譜c(k),，并尋找c(k)的最大值k0，計算方法分離為：
將數(shù)字低通濾波器搬移到以超聲波換能器諧振頻率為核心的正頻地帶中,，即：
計算c(k)的正頻率局部為：
將c'(k)搬移到零頻,，可得到c″(k)為：
求c″(k)的fft反變換c″(m)，并尋找c″(m)最大值的位置m″,，修正由于fir數(shù)字低通濾波器h(n)的群延遲帶來的的最大值m1和m″的偏差為：
同理,，能夠得出另一個方向上最大值的地位m2，從而算出待測流體流量q的估計值：
2 仿真成果的剖析
本文應(yīng)用matlab的仿真來驗證上述算法的準確性,，仿真成果如所示,。
仿真結(jié)果
超聲波信號的表白情勢為：
3 停止語
將基于fft的相干函數(shù)包絡(luò)的算法與傳統(tǒng)的方式相比擬，前者防止了回波信號的解調(diào)這一步驟,，大大下降了計算量,，同時也使提取包絡(luò)的精度得到了優(yōu)化。從計算進程來看,，其重要是fft跟逆fft的盤算,，并在仿真時程序運行速度方面優(yōu)于傳統(tǒng)的辦法，所以十分實用于目前以dsp,、arm為中心的便攜式超聲流量計體系,，這將會給實際的丈量帶來很大的方便。