張亞偉，佟仕忠，趙希鳳

（遼寧石油化工大學(xué) 信息與控制工程學(xué)院，遼寧 撫順 113001）

Chirp-Z變換在LFMCW雷達(dá)液位測量中的應(yīng)用

張亞偉，佟仕忠，趙希鳳

（遼寧石油化工大學(xué) 信息與控制工程學(xué)院，遼寧 撫順 113001）

在LFMCW雷達(dá)液位測量原理的基礎(chǔ)上，分析了采用FFT算法來提高測量精度需以增加運(yùn)算量為代價(jià)的不足，引入了一種FFT-CZT算法，來提高系統(tǒng)的測量精度。最后在MATLAB環(huán)境下，仿真出雷達(dá)回波信號分別經(jīng)過FFT和FFT-CZT算法處理后的頻譜圖，經(jīng)計(jì)算比較和仿真結(jié)果分析表明，F(xiàn)FT-CZT算法不僅可以顯著提高計(jì)算效率而且有利于提高雷達(dá)液位測量精度。

LFMCW；Chirp-Z；FFT；液位測量

線性調(diào)頻連續(xù)波（LFMCW）雷達(dá)具有很高的距離分辨力、無盲區(qū)、功率小等優(yōu)點(diǎn)[1]。但傳統(tǒng)上采用的FFT方法會產(chǎn)生的“柵欄效應(yīng)”，使系統(tǒng)產(chǎn)生較大的相對誤差，而增加FFT采樣點(diǎn)數(shù)提高測量精度又會造成運(yùn)算量的成倍增加，降低了系統(tǒng)的實(shí)時(shí)性[2]。為此在FFT的基礎(chǔ)上引入Chirp-Z算法對回波中頻距離譜的主瓣進(jìn)行局部細(xì)化，來提高雷達(dá)液位的測量精度，同時(shí)又能大大減少計(jì)算量。

1 LFMCW雷達(dá)液位測量原理

LFMCW雷達(dá)測量系統(tǒng)采用鋸齒波調(diào)頻方式的原理如圖1 所示[3-4]。

發(fā)射信號的頻率在一個(gè)掃頻周期內(nèi)可表示為：

Δt后的回波信號頻率為：

由此可以得出差頻信號Δf:

圖1 LFMCW雷達(dá)測量原理Fig.1 The measurement principle of LFMCW radar

由（4）式可得:在調(diào)制參數(shù)B、T一定的條件下，目標(biāo)距離與差頻頻率成正比，因此測得差頻頻率Δf即可計(jì)算目標(biāo)距離。一般采用的方法是：對采樣后的差頻信號先做FFT處理，然后觀察差拍信號頻譜圖，其對應(yīng)峰值處的頻率即為所需的Δf。實(shí)際上由于干擾信號的存在和FFT處理中出現(xiàn)的“柵欄效應(yīng)[5]”，會對峰值頻率Δf的檢測帶來影響，而要提高精度需增加采樣點(diǎn)數(shù)，這會增加計(jì)算量延長計(jì)算時(shí)間，所以引入Chirp-Z算法解決這個(gè)問題。

2 Chirp-Z變換簡介

圖2 Chirp_Z變換在Z平面上的螺旋采樣Fig.2 Helical sampling of Chirp_Z in the Z-plane

Chirp-z變換也稱為CZT變換，是一種在Z平面上沿著螺旋線軌道計(jì)算有限時(shí)寬的Z變換方法[6]。已知序列x（n）（0≤n≤N-1）是一個(gè) N 點(diǎn)有限長序列，X（z）表示序列 x（n）的Z變換，X（k）表示序列x（n）的N點(diǎn)FFT采樣點(diǎn)。由信號處理理論知：X（k）是x（n）的 Z變換在Z平面單位圓上的點(diǎn)等間隔采樣，而Chirp-Z變換則是在Z平面的一段螺旋線上作M點(diǎn)等間隔采樣[7]，如圖2所示。

式中：M為采樣點(diǎn)的總數(shù)，不一定與x（n）的長度N相等。A為采樣軌跡的起始點(diǎn)位置，由它的半徑A0及相角θ0確定。當(dāng)A0=l，W0=1時(shí)，就可以實(shí)現(xiàn)在單位圓上局部細(xì)化采樣。θ0指定采樣的起始點(diǎn)，指定采樣間隔。

3 Chirp-Z變換算法

LFMCW雷達(dá)測量系統(tǒng)采用Chirp-Z變換原理是:首先對差頻信號做N點(diǎn)的FFT運(yùn)算，在其頻譜圖中找出最大峰值點(diǎn)K，計(jì)算（K-1）和（K+1）兩點(diǎn)頻率，在這兩點(diǎn)頻率段內(nèi)再做M點(diǎn)的CZT運(yùn)算細(xì)化頻譜，細(xì)化后的局部頻譜圖中最大點(diǎn)對應(yīng)的頻率值即為所求差頻頻率，算法步驟如下:

1）初始化：對采樣后的雷達(dá)回波差頻信號進(jìn)行去線性和加窗處理，以減少頻譜泄露；

2）計(jì)算N點(diǎn)FFT，找到幅度最大的峰值點(diǎn)K，計(jì)算該點(diǎn)的頻率值存入f；

3）計(jì)算（K-1）和（K+1）兩點(diǎn)處的頻率值，將兩點(diǎn)的頻率值存入 f1，f2；

4）在f1至f2區(qū)間內(nèi)做M點(diǎn)CZT變換，求出計(jì)算的M點(diǎn)CZT絕對值中的最大值點(diǎn)Kmax，計(jì)算該點(diǎn)的頻率值存入fmax；

5）將 fmax代入式（3）計(jì)算距離。

取參數(shù)N=1 024，M=64，顯然，由表1可以看出在保證同一測量精度的前提下，直接采用FFT的運(yùn)算量是FFT-CZT算

表1 2種不同算法在同一測量精度下的運(yùn)算量Tab.1 The calculation amount of two different algorithms under the same measurement accuracy

法的27倍。

4 仿真實(shí)驗(yàn)結(jié)果分析

該實(shí)驗(yàn)仿真參數(shù)如下：中心頻率f0=24.7 GHz，調(diào)頻帶寬B=1 GHz，調(diào)頻周期T=10 ms。在MATLAB的仿真環(huán)境下，圖3和圖4是理想回波信號分別進(jìn)行FFT和FFT-CZT算法處理的頻域波形圖。設(shè)定延遲時(shí)間為100 ns，則目標(biāo)距離的理論值為15 m，圖3是對差頻信號利用MATLAB進(jìn)行FFT仿真得到10 121.46 Hz，根據(jù)前面推出的式（4）可得到所測距離為15.182 m。圖4是對差頻信號利用MATLAB進(jìn)行FFT/CZT仿真得到的頻率值為10 026.32 Hz，所測得的距離為15.039 m。

圖3 回波信號直接經(jīng)過FFT變換后的頻譜圖Fig.3 The spectrum diagram of the FFT processed echo signals

圖4 回波信號經(jīng)過Chirp-Z變換后的頻譜圖Fig.4 The spectrum diagram of the Chirp-Z processed echo signals

由表2可以看出，使用 FFT方法，仿真的實(shí)驗(yàn)測量結(jié)果與理論距離值之間誤差較大。而采用CZT的方法細(xì)化后，通過仿真實(shí)驗(yàn)測量結(jié)果可以看出測量精度有了明顯的提高，誤差明顯減小，由此可看出用CZT方法可以提高測量精度。

表2 4種不同延遲在無噪聲情況下測量精度的仿真結(jié)果Tab.2 Simulation results of measurement accuracy of four different delay in the case of no noise

5 結(jié)束語

文中提出了一種采用Chirp-Z變換來提高LFMCW雷達(dá)液位測量精度的信號處理方法。該方法的基本原理是先采用FFT的方法測出回波中頻距離譜上主瓣的位置，然后采用Chirp-Z變換對主瓣進(jìn)行局部細(xì)化，來提高LFMCW雷達(dá)測量精度。以上分析和仿真結(jié)果都說明了FFT-CZT方法的有效性。

[1]孟欣喜，陳文會，劉小民，等.LFMCW雷達(dá)測距系統(tǒng)及其信號處理算法的設(shè)計(jì)[J].科學(xué)技術(shù)與工程，2011，11（33）：8191-8194.

MENG Xin-xi，CHEN Wen-hui，LIU Xiao-min，et al.Design of LFMCW radar ranging system and signal processing algorithm[J].Science Technology and Engineering，2011，11（33）：8191-8194.

[2]李鵬.基于FPGA的LFMCW雷達(dá)信號處理機(jī)研究 [J].火控雷達(dá)技術(shù)，2010，39（1）:33-36.

LI Peng.Research on LFMCW radar signal processor based on FPGA[J].Fire Control Radar Technology，2010，39（1）:33-36.

[3]韓輝，冷建偉，安飛.FMCW 雷達(dá)物位系統(tǒng)的算法設(shè)計(jì)[J].天津理工大學(xué)學(xué)報(bào)，2011，27（4）:27-31.

HAN Hui，LENG Jian-wei，AN Fei.Design of the algorithm method in FMCW radar measurement system[J].Journal of Tianjin University of Technology，2011，27（4）:27-31.

[4]楊萬海.雷達(dá)系統(tǒng)系統(tǒng)建模與仿真[M].西安：西安電子科技大學(xué)出版社，2007.

[5]朱慶彬.FMCW液位測量雷達(dá)系統(tǒng)控制及回波信號處理電路設(shè)計(jì)[D].大連：大連海事大學(xué)，2011.

[6]余海強(qiáng)，馬波，江志農(nóng)，等.Chirp-Z變換在振動信號處理中的應(yīng)用[J].計(jì)算機(jī)仿真，2010（1）:319-321.

YU Hai-qiang，MA Bo，JIANG Zhi-nong，et al.Application of Chirp-Z transform in vibration signal processing[J].Computer Simulation，2010（1）:319-321.

[7]馮維婷.Chirp-Z變換在CW雷達(dá)高精度速度測量中的應(yīng)用[J].電子科技，2007（9）:6-8.

FENG Wei-ting.Application of Chirp-Z transform in velocity measurement of CW radar[J].Electronic Science and Technology，2007（9）:6-8.

The application Chirp-Z transform in LFMCW radar liquid level measurement

ZHANG Ya-wei， TONG Shi-zhong， ZHAO Xi-feng
（School of Information and Control Engineering， Liaoning Shihua University， Fushun 113001， China）

Based on the principle of LFMCW radar liquid level measurement，this paper analyzed the deficiency of FFT algorithms， which improved the measuring accuracy at the expense of increasing the amount of computation， and introduced an FFT-CZT algorithm to improve the measurement accuracy of the system.Finally the author simulated the spectrogram of the FFT and FFT-CZT processed radar echo signal respectively in MATLAB.Analysis of the theoretical calculations and MATLAB simulation results showed that the FFT-CZT algorithm not only can significantly improved the computational efficiency，but also increased the measurement accuracy of LFMCW radar liquid level gauge.

LFMCW；Chirp-Z；FFT；liquid level measurement

TP23

A

1674－6236（2013）08-0102-03

2012-11-27稿件編號201211235

張亞偉（1985—），女，山東菏澤人，碩士研究生。研究方向：儀器儀表。