晉良念， 王 燃

(1. 桂林電子科技大學(xué)信息與通信學(xué)院， 廣西桂林 541004；2. 廣西無線寬帶通信與信號處理重點(diǎn)實(shí)驗(yàn)室, 廣西桂林 541004)

0 引言

近年來，自動駕駛、智慧交通、區(qū)域監(jiān)控以及安防等領(lǐng)域的需求飆升，這些領(lǐng)域的相關(guān)應(yīng)用場景要求傳感器能夠探測環(huán)境中目標(biāo)的空間三維信息。相比于其他感知環(huán)境的傳感器如相機(jī)、激光雷達(dá)來說，商業(yè)的毫米波雷達(dá)有著獨(dú)特的優(yōu)勢。不僅因其大小適合安裝在汽車的各個位置，能很好地適應(yīng)汽車和工業(yè)領(lǐng)域?qū)走_(dá)設(shè)備體積的要求，而且探測距離遠(yuǎn)，對各種惡劣天氣環(huán)境的魯棒性強(qiáng)，在復(fù)雜工作環(huán)境下能正常運(yùn)行。而毫米波FMCW MIMO雷達(dá)高分辨的3D點(diǎn)云成像技術(shù)將有助于提高上述應(yīng)用中目標(biāo)檢測、分類、預(yù)警等功能的效率與準(zhǔn)確性。

現(xiàn)有毫米波FMCW MIMO雷達(dá)3D點(diǎn)云成像方法包括三維快速傅里葉變換(3D FFT)、快速傅里葉變換-二維多重信號分類(FFT-2D MUSIC)和三維多重信號分類(3D MUSIC)等算法。其中，3D FFT是一種基于快速傅里葉變換算法的三維參數(shù)估計方法。FFT-2D MUSIC方法利用FFT估計距離、2D MUSIC估計方位角和俯仰角，與3D FFT相比，提高了空間角分辨能力。3D MUSIC算法在此基礎(chǔ)上將2D MUSIC擴(kuò)展到3D MUSIC，提出一種距離-多普勒-方位角聯(lián)合估計方法，它通過三維張量構(gòu)造新的信號模型，并利用前向后向空間平滑技術(shù)對數(shù)據(jù)進(jìn)行疊加和平滑，獲得了三維參數(shù)的高分辨。但是，汽車和工業(yè)毫米波雷達(dá)系統(tǒng)往往方位角的視場(FOV)很大，距離分辨率高，導(dǎo)致劃分的方位角網(wǎng)格數(shù)和距離單元數(shù)多。諸如上述結(jié)合MUSIC的點(diǎn)云成像方法，既需要已知源數(shù)量和大量的快照，又要對協(xié)方差矩陣的奇異值求解，以及聯(lián)合估計的多層循環(huán)搜索，都極大增加了算法復(fù)雜度，算法運(yùn)行的時間往往成倍增加。為了滿足實(shí)時性的要求，TI公司提出了一種 FFT測距加Capon測角的三維點(diǎn)云成像方法，一方面利用Capon測角沒有特征值分解帶來的復(fù)雜運(yùn)算，另一方面通過加入CFAR檢測算法，檢測存在目標(biāo)的單元，來避免多維循環(huán)搜索，這種信號處理框架大大降低了方法運(yùn)行的時間。

盡管如此，Capon算法的角分辨不足，因而為了改善到達(dá)角(DOA)估計的性能，文獻(xiàn)[10]提出一種基于漸近最小方差的稀疏迭代方法(SAMV)。該方法在1D角度上劃分網(wǎng)格，利用采樣數(shù)據(jù)協(xié)方差矩陣與真實(shí)協(xié)方差矩陣的差異，根據(jù)AMV準(zhǔn)則，推導(dǎo)出網(wǎng)格功率的迭代公式，從而得到更好的角度譜，有著比Capon和MUSIC更佳的分辨能力。本文在SAMV的基礎(chǔ)上，將其擴(kuò)展至2D，提出一種適用于毫米波FMCW MIMO雷達(dá)的三維點(diǎn)云成像方法。首先對差拍信號做FFT，得到1D距離像；然后在每個距離單元的水平方向上劃分方位角網(wǎng)格，構(gòu)造過完備的方向矢量矩陣，從AMV準(zhǔn)則的角度出發(fā)，得出網(wǎng)格功率的迭代求解公式，由此得到距離-方位角熱圖；接著對熱圖的兩個維度分別做1D CASO-CFAR，得到目標(biāo)的距離和方位角單元；最后根據(jù)檢測出來的單元位置信息，在俯仰角上劃分網(wǎng)格，采用SAMV算法估計網(wǎng)格的信號功率，經(jīng)過檢測得到俯仰角。相較于現(xiàn)有方法，該方法在天線數(shù)量受限和有限快拍的情況下，能明顯提高方位角和俯仰角上的分辨率。

1 信號模型

FMCW MIMO雷達(dá)系統(tǒng)模型如圖1所示，有3個發(fā)射陣元～和4個接收陣元～，其中和的間距為2(為雷達(dá)信號的波長)，與和的間距均為，并向軸正方向平移2；到的陣元間距為2。為了方便后續(xù)的建模，按照文獻(xiàn)[11]，將圖1的MIMO陣列結(jié)構(gòu)等效為圖2的12個虛擬陣元收發(fā)共置來進(jìn)行處理，對應(yīng)的坐標(biāo)位置矩陣=[,…,,…,]∈，其中=?+?，∈和∈分別是發(fā)射陣元與接收陣元的3D坐標(biāo)矩陣。

圖1 MIMO陣列模型

圖2 等效虛擬陣列模型

MIMO陣列采用時分發(fā)射模式。以為例，假設(shè)在一個調(diào)頻周期內(nèi)發(fā)射的信號為

()=exp[j(2π+π)],0≤≤

(1)

式中，為信號時間，為每個Chirp的調(diào)頻周期，為調(diào)頻斜率，為中心頻率。根據(jù)圖2，接收的回波信號可以表示為

(2)

式中，為目標(biāo)的反射系數(shù)，為發(fā)射天線到接收天線的信號傳播時延，包括距離和角度時延，即

(3)

式中，為目標(biāo)和原點(diǎn)的徑向距離，=[sincos,coscos,sin]為目標(biāo)的單位方向矢量。中的和為目標(biāo)的方位角和俯仰角，取值范圍為方位角和俯仰角的有效視場角和內(nèi)，即∈(-2,2),∈(-2,2)。

將回波信號與發(fā)射信號進(jìn)行混頻，得到差拍信號為

exp[j2π(+++

(4)

式中，平方項(xiàng)和二次項(xiàng)相對于其他項(xiàng)可以忽略，又因的值太小，exp[j2π()]項(xiàng)也可以忽略，由此改寫為

()=exp[j2π(·)]×

exp[j2π(2+2)]+

()

(5)

對()進(jìn)行采樣得到離散采樣序列：

()

(6)

式中，為采樣速率，=1,2,…,為采樣點(diǎn)數(shù)。若將式(6)推廣到個目標(biāo)，設(shè)目標(biāo)距離分別為,,…,，所在的方位角和俯仰角分別為(,),(,),…,(,)，則式(6)的信號模型可以進(jìn)一步寫為

()

(7)

()

(8)

2 算法描述

為了從式(8)得到3D高分辨點(diǎn)云像，需經(jīng)過一維距離像估計、距離-方位角熱圖生成、CFAR目標(biāo)檢測和俯仰角估計等一系列的信號處理，具體流程如圖3所示。

圖3 三維點(diǎn)云成像方法流程框圖

2.1 一維距離像估計

如圖2所示，就第個虛擬陣元，先對第

=1,2,…,

(9)

2.2 距離-方位角熱圖生成

(10)

(11)

=1,…,,=1,…,

(12)

()=()+()

(13)

上述式(13)為用2D SAMV聯(lián)合方位角和俯仰角估計時的信號模型，將其用于方位角網(wǎng)格功率估計時，需令=0，而且等于1和2，等于1到4。將()重新組合，得到

(14)

()=()+()

(15)

(16)

式中，()=diag((),…,(),…,())，對應(yīng)的()為方位角上第個網(wǎng)格上的信號功率，()為噪聲功率，()和()為待求的參數(shù)， 對協(xié)方差矩陣()矢量化，即

(())=(())=()

(17)

式中，()=[(),…,(),…,(),()]，為系數(shù)矩陣，即=[(,0)?(,0),…,(,0)?(,0),…,(,0)?(,0),()]，這里的?代表Kronecker積。對矩陣()的矢量表示可以看作未知的參數(shù)矢量()中元素的線性組合。由于Chirp數(shù)有限，可以得到采樣協(xié)方差矩陣為

(18)

(19)

(20)

(21)

(22)

式中，為迭代次數(shù)。根據(jù)Capon波束形成算法對網(wǎng)格功率的估計，有=1/，則式(21)、式(22)可改寫為

(23)

(24)

在上述迭代過程中，初始值可以通過周期圖法來求解，即

(25)

(26)

=[(1),…,(),…,()]∈×

(27)

2.3 CASO-CFAR檢測

將式(27)以距離維度和方位角維度分別送入CASO-CFAR檢測器，其中=[(),…,(),…,()]，=[(1),…,(),…,()]，對每個方位角網(wǎng)格上的所有距離單元做1D CASO-CFAR，檢測出存在目標(biāo)的距離單元，再對距離單元處的所有方位角網(wǎng)格做1D CASO-CFAR，檢測出存在目標(biāo)的距離-方位角網(wǎng)格單元(,d)，得到對應(yīng)的(,)。由于在應(yīng)用中我們往往需要將距離較近的值能夠檢測出來，文獻(xiàn)[15]中表明，單元平均選小恒虛警檢測(CASO-CFAR)相對于單元平均恒虛警檢測(CA-CFAR)，使得距離較近的目標(biāo)檢測出來而不出現(xiàn)目標(biāo)遮蔽的現(xiàn)象。具體處理結(jié)構(gòu)如圖4所示。

圖4 CASO-CFAR處理結(jié)構(gòu)圖

2.4 俯仰角估計

接下來，對檢測的距離-方位角網(wǎng)格單元(,d)進(jìn)行俯仰角估計時，根據(jù)式(13)，需令=，=d，且等于1到3，等于1到4，將()重新組合，得到

(28)

()=()+()

(29)

(30)

(31)

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

3.1 雷達(dá)參數(shù)設(shè)計

為了對近距離區(qū)域內(nèi)人員等目標(biāo)的3D位置參數(shù)進(jìn)行探測，需要進(jìn)行合理的毫米波雷達(dá)參數(shù)設(shè)計。雷達(dá)參數(shù)主要包括MIMO雷達(dá)天線數(shù)量及陣列設(shè)計、調(diào)頻起始頻率、調(diào)頻斜率、ADC采樣點(diǎn)數(shù)、ADC采樣率、工作波長和每幀的調(diào)頻周期數(shù)，這些參數(shù)由應(yīng)用場景的最大測量距離、距離分辨率等指標(biāo)來確定，以下指標(biāo)分別滿足公式：

=×

(32)

(33)

(34)

(35)

雷達(dá)硬件支持的參數(shù)有上限，將其作為雷達(dá)參數(shù)設(shè)計的先決條件，包括發(fā)射功率決定的最大測量距離、最大的ADC采樣率與最小的ADC采樣率、最大的掃頻帶寬、最大的調(diào)頻斜率與最小的調(diào)頻斜率、水平及垂直方向最大視場和等。具體的參數(shù)設(shè)置如表1所示。

表1 雷達(dá)系統(tǒng)和目標(biāo)參數(shù)

3.2 仿真結(jié)果分析

設(shè)置4個點(diǎn)目標(biāo)，參數(shù)設(shè)置如表2所示。

表2 仿真目標(biāo)三維參數(shù)

圖5和圖6分別為文獻(xiàn)[9]方法和本文方法生成的距離-方位角熱圖以及經(jīng)過CASO-CFAR檢測的結(jié)果。從圖中可以看出，文獻(xiàn)[9]的方法已經(jīng)無法在分辨出方位角上相差5°的目標(biāo)1、2和目標(biāo)3，只在目標(biāo)1、2的方位角處有亮點(diǎn)，而方位角相差更大的目標(biāo)4處也有明顯的亮點(diǎn)。同時該方法的旁瓣較高，導(dǎo)致經(jīng)過CFAR檢測之后，真實(shí)距離-方位角單元的周圍許多單元也被檢測出來。相比之下，本文方法不僅在4個目標(biāo)的真實(shí)方位角位置處均有明顯的亮點(diǎn)，而且旁瓣較低，真實(shí)位置周圍只有極少的單元被檢測出來。圖7和圖8分別為兩種方法在檢測單元處估計的俯仰角網(wǎng)格功率譜圖。從圖中可以看出，文獻(xiàn)[9]方法在所有檢測單元處都只有一個譜峰，無法分辨出不同俯仰角的目標(biāo)1和2，而本文方法的俯仰角功率譜在不同的檢測單元上有兩個不同位置的亮點(diǎn)。圖9和圖10分別為兩種方法在空間直角坐標(biāo)系中生成的3D點(diǎn)云像。從圖中可以看出，文獻(xiàn)[9]方法生成的點(diǎn)云像只能分辨出兩個目標(biāo)，而本文方法4個目標(biāo)的點(diǎn)云都很清晰分辨出來，成像效果很好，只有目標(biāo)4的俯仰角估計有5°的偏移。由此可見，相比較文獻(xiàn)[9]的方法，本文提出的方法不僅在方位角上獲得了更好的分辨率，而且在只有兩排天線的條件下，俯仰角上依然可以得到一定的分辨率。

(a) 距離-方位角熱圖

(a) 距離-方位角熱圖

圖7 文獻(xiàn)[9]方法生成的俯仰網(wǎng)格功率譜圖

圖8 本文方法生成的俯仰網(wǎng)格功率譜圖

圖9 文獻(xiàn)[9]方法生成的3D點(diǎn)云像

圖10 本文方法生成的3D點(diǎn)云像

3.3 實(shí)測數(shù)據(jù)處理結(jié)果

實(shí)驗(yàn)所采用雷達(dá)系統(tǒng)由TI(德州儀器)公司的毫米波雷達(dá)IWR6843ISK評估板卡、毫米波傳感器承載卡平臺MMWAVEICBOOST和數(shù)據(jù)采集卡DCA1000組成。評估板卡可以設(shè)置雷達(dá)系統(tǒng)參數(shù)，天線系統(tǒng)如圖1所示，將差頻信號在設(shè)置的采樣率下進(jìn)行I/Q采樣，之后將采樣好的信號通過工業(yè)載板傳輸給DCA1000，DCA1000將處理好的差頻采樣復(fù)信號打包成二進(jìn)制文件通過網(wǎng)口傳輸給電腦。

實(shí)驗(yàn)時的場景如圖11所示，在距離雷達(dá)5 m處的兩個人體目標(biāo)，人體目標(biāo)的高度大約為1.7 m，人體的寬度大約為0.5 m，厚度大約為0.2 m，兩人相距為1 m，雷達(dá)的高度為1 m，圖12和圖13分別為兩種方法生成的距離方位熱圖及CFAR檢測結(jié)果。從圖中可以看出，文獻(xiàn)[9]中的方法，已無法分辨在方位角上相差5°的兩個人體目標(biāo)，而本文方法依然能夠清晰地分辨出兩個目標(biāo)。圖14和圖15分別為兩種方法在檢測單元處估計出來的俯仰角網(wǎng)格功率譜圖。從圖中可以看出，文獻(xiàn)[9]方法由于Capon測角算法的分辨率有限，且雷達(dá)在俯仰維上只有兩排天線，無法有效地估計目標(biāo)所在的俯仰角范圍，而本文方法依然在目標(biāo)的俯仰角范圍有明顯的亮斑。圖16和圖17分別為兩種方法的三維點(diǎn)云成像結(jié)果。從圖中可以看出，文獻(xiàn)[9]方法只能分辨出一個目標(biāo)的點(diǎn)云像，且因俯仰角上分辨率很低，只能取峰值位置，導(dǎo)致點(diǎn)云成像結(jié)果和目標(biāo)輪廓相差很大。而本文方法能清晰地生成兩個目標(biāo)的點(diǎn)云，同時由于SAMV方法估計出來的是網(wǎng)格的真實(shí)功率，經(jīng)過檢測算法后，三維點(diǎn)云結(jié)果和目標(biāo)的輪廓十分相似。

圖11 實(shí)驗(yàn)場景

(a) 距離-方位角熱圖

(a) 距離-方位角熱圖

圖14 文獻(xiàn)[9]方法生成的俯仰網(wǎng)格功率譜圖

圖15 本文方法生成的俯仰網(wǎng)格功率譜圖

圖16 文獻(xiàn)[9]方法生成的3D點(diǎn)云像

圖17 本文方法生成的3D點(diǎn)云像

4 結(jié)束語

本文提出一種毫米波FMCW MIMO雷達(dá)高分辨三維點(diǎn)云成像方法，通過毫米波雷達(dá)的大帶寬，得到高分辨的距離信息，引入SAMV方法，提升雷達(dá)的測角能力，結(jié)合CASO-CFAR算法，減少點(diǎn)云成像的時間。通過仿真和實(shí)驗(yàn)結(jié)果表明，相對于現(xiàn)有方法，該方法即使在天線數(shù)目有限的商業(yè)毫米波雷達(dá)上，依然有著更高的分辨率，兼顧了實(shí)時性的要求。本方法生成的近距離人體目標(biāo)三維點(diǎn)云像，點(diǎn)云輪廓更加清晰、準(zhǔn)確，便于后續(xù)基于點(diǎn)云的目標(biāo)分類、姿勢識別和計數(shù)等雷達(dá)應(yīng)用的開發(fā)。