0 引言

    雷達測速的工作原理是以多普勒效應為基礎的，在科學和工程中有廣泛的應用^[1]。運動目標的雷達回波會產生多普勒頻移，而多普勒頻移與目標運動速度成正比，通過對多普勒頻移的精確測量，可實現對目標運動速度的精確測量。

    IVS-179是Innosent公司推出的一款K-波段帶VCO的雷達前端模塊，其發射頻率為24 GHz～24.25 GHz；發射功率20 dBm；采用平面微帶天線陣列結構，波束寬度水平方向為7°，垂直方向為28°；旁瓣抑制比大于15 dB。該雷達模塊性能優良，操控簡單，非常易于系統集成，應用十分廣泛。

    基于時移相位差的測頻算法是一種常用的頻率測量算法，但是由于常規的FFT不能有效地克服頻譜泄漏，且對信號頻率、相位的測量都依賴于一些近似的估計，因而測量精度不能達到令人滿意的程度^[2]。

    全相位FFT(apFFT)具有優良的抑制頻譜泄漏性能^[3]。尤其突出的是，apFFT具有“初始相位不變”的重要特性，即對單頻正弦信號作apFFT變換，最高譜線的初始相位嚴格等于該正弦信號的初始相位，而不受頻譜擴散的影響。因此，利用apFFT的“初相不變性”可以大大提高時移相位差法的頻率測量精度，從而大大提高多普勒雷達的測速精度。

1 apFFT時移相位差法的算法原理

1.1 apFFT算法原理

    apFFT算法可以分解為數據預處理和FFT兩部分，其算法原理如圖1所示，其中N是FFT點數。

    設一個長度為(2N-1)的輸入序列：

    數據預處理是指利用卷積窗W，將該輸入序列按圖1所示的運算映射為長度為N的數據向量：

    然后，對y序列做N點的FFT，從而實現apFFT變換。其中，根據數據處理的要求選擇不同的卷積窗W，卷積窗類型可分為無窗、單窗和雙窗。

1.2 基于apFFT的時移相位差法

    設單頻信號序列為：

2 雷達測速系統的組成與工作原理

    雷達測速系統組成如圖2所示。系統主要包括：多普勒雷達前端模塊、中頻信號調理電路、模數轉換電路和數模轉換電路、FPGA數字信號處理平臺、以太網接口電路、外部大容量存儲器等。

    在系統上電或復位之后，FPGA首先通過參數配置接口初始化雷達工作參數，再通過DAC電路產生雷達工作所需控制信號。

    在啟動測速后，雷達模塊輸出的多普勒中頻信號首先經過中頻信號調理電路的處理，再由16 bit的ADC完成模數轉換。其中，中頻信號調理電路由有源高通濾波器、程控衰減器和有源低通濾波器組成，FPGA可以通過程控衰減器控制輸入信號的幅度。

    采樣后的數據經FPGA緩存、加窗處理、頻譜變換、主譜識別和cordic求解相位^[4]等處理，最終得到主譜相位，并計算出多普勒頻率及相應的速度值。

    大容量存儲器由2片NAND Flash構成，每片容量8 Gb。在現場實測時，可將ADC采樣的數據存儲到NAND Flash中，之后，再帶回實驗室做詳細的數據分析。由于存儲器容量大，可存儲幾個小時的原始數據。

    以太網電路采用以太網接口控制芯片W5300。通過以太網完成控制命令的下發和測量結果的上傳工作。控制命令包括雷達工作模式、雷達啟閉、程控衰減量等，上傳的數據包括ADC采樣數據和測頻結果。

3 雷達測速系統的數據采集與處理

    ADC采樣后的數據分為兩路，一路數據通過以太網接口上傳至PC，在PC中存為TXT文本文件，并在MATLAB環境下作數字信號處理^[5]；另一路數據直接在FPGA平臺上進行計算，計算后的結果再通過以太網上傳到PC。二者處理的結果可作對比實驗。

3.1 數據的采樣與存取

    數據采樣與處理參數：

    采樣點數3N：6 144；采樣率：200 kS/s；ADC數據量化位數：16 bit；FFT點數N：2 048；FFT的頻率分辨率：97.66 Hz；窗類型：Blackman雙窗。

    經ADC采樣后的數據存儲到FPGA內部的RAM中，RAM的大小為3N字。緩存區操作流程：若寫地址大于2N-1，表明前2N-1個點寫滿，開始讀取地址為0～2N-2的數據；若寫地址大于0小于N，則讀取地址為N～3N-2的數據，不滿足讀數據的條件，則等待。在讀寫RAM操作過程中，讀數據時鐘遠遠高于寫數據時鐘，以保證讀取新數據之前，上一次數據已處理完成。

3.2 FPGA數據的處理流程

    FPGA數據處理流程主要包括：讀寫RAM、加窗處理、apFFT運算、查找主譜、計算主譜相位和計算結果的緩存等。

    為了抑制頻譜泄漏，突出頻譜主瓣，對數據采取加窗處理。

    經過乘法器的輸出數據過寬，需要截位。在模塊控制器的控制下，前N-1點截位后，壓入FIFO中緩存；后N點到來后，FIFO中數據出棧，并和截位后的數據相加，完成apFFT算法的預處理過程。

    查找主譜，即尋找頻譜中最高譜線。確定最高譜線后，需要計算其相位，計算主譜相位通過cordic算法實現。當主譜相位計算完成以后，緩存一次，緩存的結果就是每3N點的前后2N-1點的兩組apFFT主譜序號及其相位值。

    每完成一次3N點計算后，結果被存儲到RAM中，每完成一次速度測量，再將RAM中數據讀出并通過以太網傳到上位機。

    圖3為FPGA的apFFT時移相位差法測頻結構圖。

4 測速實驗與數據分析

4.1 實驗室環境下的模擬測速實驗與數據分析

    利用信號源產生一個頻率為1.2 kHz的單頻正弦波，用該信號模擬多普勒雷達測速中頻信號，則對應的速度值為27 km/h。

    對FPGA平臺計算后所得的上傳結果，利用MATLAB繪制成速度曲線如圖4所示。

    由圖4可以看到，速度曲線在27 km/h附近上下波動。其中，最大偏差速度值為27.008 km/h，相對于真實值的誤差為0.3‰。

    由此可見，基于apFFT的時移相位差法能實現高精度的速度測量。

4.2 外場實際測速實驗與數據分析

    測速器安裝在公路的人行過街天橋上，選取單一車輛高速行駛時進行測量。

    圖5為雷達多普勒回波中頻信號，圖6為雷達回波中頻前2N-1點apFFT頻譜，圖7為雷達回波中頻后2N-1點apFFT頻譜，圖8為速度測量曲線。

    從雷達回波角度看，當被測目標速度穩定且不存在其他目標干擾時，雷達回波穩定，輸出的多普勒中頻信號質量好。從頻譜角度分析，頻譜主瓣突出。

    從速度曲線來看，曲線平穩，末端有下降趨勢，這是由于被測目標靠近雷達時，被測物體運動方向和目標與雷達連線的夾角逐漸增大導致多普勒頻率減小所引起的。

    在描繪速度曲線的數據中，選擇一段連續數據制成表1，根據表1中的結果發現，最大與最小頻率差Δf=3.95 Hz，其變化率僅為1.4‰。

    同樣的，由頻率對應計算的速度值相對變化量也比較小，最大速度與最小速度之差僅為0.09 km/h，變化率仍為1.4‰。

    因此，在實際的應用中，基于apFFT的時移相位差法同樣能實現高精度的速度測量。

5 結論

    本文設計了以FPGA為處理核心的多普勒雷達測速系統，并將基于apFFT的時移相位差法這一高精度測頻算法應用于實際測速中；利用MATLAB的信號處理功能和GUI設計了上位機的可視化界面，以便于及時觀察時頻域數據和測速結果。多次實測結果表明，該雷達測速系統完全達到設計要求，能夠實現目標運動速度的高精度測量。

參考文獻

[1] 劉麗華．多普勒雷達測速系統研究[D]．武漢：華中科技大學，2007.

[2] 賀同，陳星，洪龍龍.基于FPGA的全相位FFT高精度頻率測量[J].電子測量技術，2013，36(8)：80-83.

[3] 王兆華，黃翔東.數字信號全相位譜分析與濾波技術[M].北京：電子工業出版社，2009.

[4] 王夢源，王書省，陳星．流水線雙模CORDIC算法的FPGA實現[J]．電子測量技術，2007，30(9)：184-186.

[5] 陳垚光，毛濤濤，王正林，等.精通MATLAB GUI設計[M].第2版．北京：電子工業出版社，2011.

作者信息:

張  明，陳  星，王宇飛

（北京航空航天大學 電子信息工程學院，北京100191）