李秀麗，劉紀(jì)元，焦學(xué)峰，陳飛虎

（1.沈陽航空航天大學(xué) 電子信息工程學(xué)院，遼寧 沈陽 110136；2.中國科學(xué)院 聲學(xué)研究所綜合聲納實(shí)驗(yàn)室，北京 100080；3.北京航瑞博泰科技有限公司 北京 100102）

近年來，隨著電子技術(shù)的發(fā)展，也加快了水利系統(tǒng)現(xiàn)代化建設(shè)的發(fā)展進(jìn)程。在明渠監(jiān)測(cè)方法中，非接觸式測(cè)量方法由于以其測(cè)量速度快，安裝安全方便等特點(diǎn)逐步成為發(fā)展的主流趨勢(shì)，特別適用于含沙量大，淤泥較多的渠道。非接觸式測(cè)量的主流產(chǎn)品是電波流速儀。

電波流速儀是用雷達(dá)多普勒效應(yīng)來測(cè)量水流速度的，其主體思路是將雷達(dá)原理和信號(hào)處理結(jié)束相結(jié)合。而考慮到我們所要測(cè)定的是水流的速度，所以選 擇使用連續(xù)波雷達(dá)來實(shí)現(xiàn)。即利用發(fā)射信號(hào)與回波信號(hào)的多普勒頻率差來提取速度信息。

1 電波流速儀的原理

如圖1所示，雷達(dá)照射水面時(shí)，部分電磁波被水面散射，且只有向后散射的那部分電磁波能量構(gòu)成回波信號(hào)，其余電磁波則折射入水。由于波浪底下的水流集體作為波浪的載體，具有和波浪相同的運(yùn)動(dòng)速度，所以波浪上的回波產(chǎn)生的多普勒頻移必然會(huì)反應(yīng)水面流速，我們將接收到的信號(hào)頻率相對(duì)與發(fā)射的信號(hào)頻率關(guān)系定義為：

圖1 原理圖Fig.1 s diagram

其中，fD為多普勒頻率；f0為發(fā)射頻率；fr為接收到的回波頻率。

多普勒頻率與速度之間的關(guān)系表示為：

由式（2）可以推出：

其中，fD為多普勒頻率；f0為發(fā)射頻率；v為水流速度；c為光速；θ為水流方向與雷達(dá)發(fā)射器到運(yùn)動(dòng)目標(biāo)連線之間的角度。由式（3）可以看出，只要從回波信號(hào)中提取出多普勒頻偏就能求得水流的速度。

假設(shè)連續(xù)波雷達(dá)發(fā)射一個(gè)頻率為f0的0信號(hào)：

經(jīng)過時(shí)間t1后，雷達(dá)接收機(jī)接收到該信號(hào)，接收到的信號(hào)表示為：

將 x1（t）與本振信（t）=e-j（2πfDt+φ0）號(hào)混頻，來去除高頻分量，最終得到的中頻信號(hào)為：

經(jīng)過采樣，我們可以的到：

式中，fs表示采樣頻率，N表示采樣點(diǎn)數(shù)。進(jìn)行離散頻譜分析可知，用ip表示譜峰處對(duì)應(yīng)的譜線號(hào)，則測(cè)量所得的多普勒頻率可以表示為：

則，所得到的水流速度表示為：

從式（9）可以看出，水流速度與離散頻譜譜峰值的譜線號(hào)，采樣頻率，采樣點(diǎn)數(shù)有關(guān)，該式為多普勒雷達(dá)頻域測(cè)水流速度的基本公式[2]。

從上式我們可知，測(cè)速的精度主要取決于離散頻譜分析提取多普勒頻移的精度。一般測(cè)速雷達(dá)提取頻移的方法是進(jìn)行FFT，但是由FFT得到的離散頻譜譜線，目標(biāo)多普勒頻率對(duì)應(yīng)譜線通常會(huì)落到兩條譜線之間，從而產(chǎn)生一定的測(cè)速誤差。若用fD表示多普勒頻率測(cè)量值，用f′D表示目標(biāo)的真實(shí)多普勒頻率，則由于柵欄效應(yīng)所造成的目標(biāo)多普勒頻率最大誤差可以表示為：

如果改為用速度表示，則水流速度的最大誤差為：

從（12）式我們可以看出：

1）如果不采取適當(dāng)?shù)拇胧瑴y(cè)速的誤差會(huì)比較大的，不滿足系統(tǒng)精度要求。

2）由上式可以看出，要想減小FFT所引起的測(cè)速誤差，可以使采樣因子K的取值應(yīng)盡可能小，或是采樣頻率要盡可能的低，但前提是必須滿足采樣定理的要求，即K≥2。

3）只要適當(dāng)增大FFT的采樣點(diǎn)數(shù)N，也可以達(dá)到所要求的測(cè)速精度。但是在實(shí)際應(yīng)用時(shí)，F(xiàn)FT的采樣點(diǎn)數(shù)由于受到處理器件和系統(tǒng)實(shí)時(shí)性等因素的制約，不可能無條件地增大，所以這種方法受到一定的限制。

4）選擇合適的頻譜校正算法可以進(jìn)一步提高測(cè)速的精度。如果采用Chirp-z變換能夠?qū)夭l譜中的指定部分進(jìn)行局部細(xì)化，這樣就能夠?qū)崿F(xiàn)在采樣點(diǎn)數(shù)，運(yùn)算量增加不多的情況下，大大提高雷達(dá)的測(cè)量精度[2-3]。

2 Chirp-z變換

線性調(diào)頻z變換算法[4]是計(jì)算螺旋線周線上z變換的諸取樣，這些取樣在螺旋線的某一部分上按等角度分布，設(shè)x（n）表示N 點(diǎn)序列，其 z變換表示為 X（z），則給定點(diǎn) Zk上的 X（zk）為：

其中 A=A0ejθ0，Zk=AW-k，W=W0ejφ0，A0和 θ0分別為第一個(gè)取樣點(diǎn)位置的半徑和角度，其他采樣是沿螺旋周線按角度間隔φ0分布；W0控制周線盤旋的傾斜度，W0＞1，周線向內(nèi)盤旋，W0＜1，周線向外盤旋，W0=1，表示螺旋線是一段圓弧；A0=1，表示這段圓弧是單位圓的一部分[5]。

注：*表示卷積，所以最終可以表示為：

Chirp-z變換與FFT相比，不同之處在于，它是對(duì)一段z平面上的螺旋線進(jìn)行等間隔采樣，而FFT是對(duì)整個(gè)單位圓進(jìn)行采樣。為了在不增加運(yùn)算量的前提下提高測(cè)量精度，我們先對(duì)回波信號(hào)進(jìn)行FFT全景譜分析，進(jìn)行譜峰搜索找到最大的值，然后在最大值附近范圍進(jìn)行Chirp-z變化，從而得到更加精確的頻譜值[6]。

3 數(shù)據(jù)分析

3.1 仿真數(shù)據(jù)處理

根據(jù)所給的數(shù)據(jù)可以算出水流速度的真實(shí)值為：

用FFT處理后所測(cè)得的水流速度值為：

可得誤差為0.070 6m/s，而是用了Chirp-z變換處理后的水流速度為：

可得誤差為0.026 5 m/s，我們可以看出與FFT相比，經(jīng)過Chirp-z變換處理后的測(cè)量精度要高很多。圖2分別為經(jīng)過FFT與Chirp-z變換的matlab仿真圖。

同樣的方法，我們分別進(jìn)行了不同速度的測(cè)量，測(cè)得的4組數(shù)據(jù)如下表所示。

3.2 結(jié)果分析

通過對(duì)不同流速時(shí)采集到的回波信號(hào)進(jìn)行處理，對(duì)所得到的的測(cè)量結(jié)果進(jìn)行對(duì)比我們可以看出，在不改變采樣點(diǎn)數(shù)的情況下，使用Chirp-z進(jìn)行局部細(xì)化的測(cè)量方法所得到的測(cè)量精度比僅對(duì)回波信號(hào)使用FFT時(shí)提高很多。

圖2 經(jīng)過FFT與Chirp-z變換的頻譜仿真圖Fig.2 Simulation diagram of processing spectrum by using FFT and Chirp-z

表1 FFT與Chirp-z處理后的測(cè)試數(shù)據(jù)Tab.1 Test data of processing the signal by using FFT and Chirp-z

4 結(jié)束語

近年來，電波測(cè)流儀以其安裝安全方便，測(cè)速快等特點(diǎn)已經(jīng)在我國得到了迅速的發(fā)展和推廣。但對(duì)于測(cè)量的精度仍有可以提升的空間，所以在電子技術(shù)不斷發(fā)展的今天，電波流速儀必定可以進(jìn)一步提高完善，成為明渠測(cè)流的核心力量。

[1]呂宗城.電波流速儀在山區(qū)河流推算流量應(yīng)用分析[J].廣西水利水電，2012（6）:6-8.LV Zong-cheng.Application of electric wave flow velocity meter in mountainous area river current calculation[J].Guangxi Water Resources&Hydropower Engineering，2012（6）:6-8.

[2]丁康，謝明，楊志堅(jiān).離散頻譜分析校正理論與技術(shù)[M].北京:電子工業(yè)出版社，1995.

[3]馮志華，劉永斌，張健.用在頻譜細(xì)化中Chirp-Z變換的特性分析[J].信號(hào)處理，2006，22（5）:741-745.FENG Zhi-hua，LIU Yong-bin，ZHANG Jian.Character analyses of the Chirp-Z transformation used in the spectrum zooming[J].Signal Processing，2006，22（5）:741-745.

[4]徐飛，王延暴.Chirp-z變換在雷達(dá)信號(hào)處理中的應(yīng)用[J].現(xiàn)代電子技術(shù)，2011，34（9）:28-29.XU Fei，WANG Yan-bao.Application of Chirp-Z transform in radar signal processing[J].Modern Electronics Technique，2011，34（6）:28-29.

[5]丁康，潘成灝，李巍華.ZFFT與Chirp-z變換細(xì)化選帶的頻譜分析對(duì)比[J].振動(dòng)與沖擊，2006，25（6）:9-12.DING Kang，PAN Cheng-hao，LIWei-hua.Spectrum analysis comparison between ZFFT and Chirp-Z transform[J].Journal of Vibration and Shock，2006，25（6）:9-12.

[6]馮維婷.Chirp-Z變換在CW雷達(dá)高精度速度測(cè)量中的應(yīng)用[J].電子科技，2007（9）:6-7.FENGWei-ting.Application of Chirp-Z transform in velocity measurement of CW radar[J].Slectronic Science and Technology，2007（9）:6-7.