趙繼廣,宋一鑠,張智詮

(1.裝甲兵工程學(xué)院控制工程系,北京 100072;2.裝備指揮技術(shù)學(xué)院,北京 101416)

0 引言

線性調(diào)頻連續(xù)波(Linear Frequency Modulation/Continuous Wave,LFM/CW)體制測距系統(tǒng),通過對回波信號與本振信號混頻后得到的差頻信號頻率的測量實現(xiàn)對目標(biāo)距離的測量[1]。該測距系統(tǒng)在一些距離范圍內(nèi)獲得的差頻信號的初始相位近似為0,目前一般利用傅里葉變換進(jìn)行時頻分析[2]從而得到差頻信號的頻率,但由于其測距分辨力不夠高且運算量大,因此難于實現(xiàn)實時多距離分辨。為此,本文提出利用余弦變換進(jìn)行時頻分析的方法。

1 傅里葉變換與余弦變換

傅里葉變換是進(jìn)行頻譜分析的重要方法,余弦變換是傅里葉變換的一部分。傅里葉變換具有普適性,對于一般信號都能給出正確的頻譜。具體到一些擁有特殊性質(zhì)的信號,可以考慮利用其他時頻分析方法處理,從而突出某一方面的性能。例如對于初始相位為0(或近似為0)的余弦信號,余弦變換不僅能夠給出近似正確的頻譜,還具有較小的運算量[3]和更高的頻率分辨力。

1.1 傅里葉變換

在實際應(yīng)用中受器件和環(huán)境等因素的影響,所得到的差頻信號的信噪比較低。由于傅里葉變換對噪聲抑制能力強(qiáng),所以一般用其獲得信號頻譜信息。根據(jù)傅里葉變換的定義,連續(xù)非周期信號s(t)的傅里葉變換可以表示為[4]:

將式(1)展開,得到

對上式取模,即得到信號的頻譜。由(2)式可知:將信號s(t)分別與各種頻率下的余弦和正弦信號進(jìn)行混頻,再將每一個頻率下的兩個混頻結(jié)果取模就得到了該頻率下的頻譜幅值。由于連續(xù)非周期信號等同于周期無限大的周期信號,所以得到的頻譜應(yīng)該被理解為頻譜密度。信號s(t)中包含的任一頻率分量的初始相位都是未知的,頻譜分析不應(yīng)受到未知初始相位的影響,所以傅里葉變換采用正交混頻并取模的方式,使得只要信號s(t)中存在該頻率分量,計算結(jié)果就會在該頻率點上得到峰值且不會受到其初始相位的影響。

對頻率為 fIF,時長為 T的余弦信號,無論其初始相位如何,對其進(jìn)行傅里葉變換后,其頻域?qū)⒊尸F(xiàn)一個峰值在 f IF處的SINC函數(shù),如圖1所示。

圖1 1 ms時長矩形窗函數(shù)的傅立葉變換頻譜圖Fig.1 The frequency spectrum of 1 ms-lasting rectangular window in Fourier transform

1.2 余弦變換

連續(xù)非周期信號S IF的余弦變換為[5]:

在余弦變換下,由于沒有涉及正交混頻,所以信號的初始相位不同會引起頻譜誤差,即通過余弦變換得到的頻譜會因為信號的初始相位而出現(xiàn)頻譜偏移和畸變。例如對一時長為T、頻率為 f 1且初始相位φ0不為0的余弦信號

作余弦變換,得到信號頻譜

由式(7)可知用余弦變換得到的頻譜函數(shù)在差頻信號頻率點上的幅度會受到初始相位的影響。只有當(dāng)信號初始相位為0(或近似為0)時,才能得到和傅里葉變換一致(或近似一致)的結(jié)果,即在頻域呈現(xiàn)一個峰值在 f IF處的SINC函數(shù),如圖2所示。

圖2 1 ms時長矩形窗函數(shù)的余弦變換頻譜圖Fig.2 The frequency spectrum of 1 ms-lasting rectangular window in cosine transform

2 線性調(diào)頻連續(xù)波時頻分析的余弦變換原理

LFM/CW體制測距系統(tǒng)將回波信號與本振信號混頻得到差頻信號,該差頻信號的頻率與目標(biāo)距離間存在線性關(guān)系[1]。在實際應(yīng)用中,綜合考慮信號信噪比和抗干擾能力等因素,一般利用傅里葉變換提取差頻信號的頻率。如果在一次測量中存在兩個或以上不同距離的目標(biāo),那么得到的差頻信號中也含有與這些目標(biāo)距離對應(yīng)的頻率。分辨開這些頻率就分辨開了這些不同的距離。由第一節(jié)中傅里葉變換和余弦變換的性質(zhì)可知,對于時長為 T的余弦信號,兩種變換下的頻率分辨力分別為 T-1和0.5T-1,但使用余弦變換時要求信號的初始相位近似為零以避免頻譜偏移和失真。LFM/CW體制測距系統(tǒng)得到的差頻信號由一個或多個余弦信號組成,通過改變系統(tǒng)參數(shù),這些余弦信號的初始相位在一些距離范圍內(nèi)近似為0,此時利用余弦變換可得到近似正確的信號頻譜。這使余弦變換在 LFM/CW測距中的應(yīng)用成為可能,即可利用余弦變換處理差頻信號,使得在差頻信號時長不變的情況下,將頻率分辨力提高到使用傅里葉變換時的兩倍。

3 余弦變換在LFM/CW體制測距系統(tǒng)中的應(yīng)用

3.1 余弦變換下LFM/CW體制測距系統(tǒng)的測距分辨力

具體到 LFM/CW 體制測距系統(tǒng)[6],由于一個周期內(nèi)觀察到的差頻信號的頻率為fIF,時長約等于掃頻周期T,故對其進(jìn)行傅里葉變換后,其頻域?qū)⒊尸F(xiàn)一個峰值在 fIF處的SINC函數(shù),此時SINC函數(shù)的主瓣寬度為 T-1,頻率分辨力為T-1,即如果兩個差頻信號相差小于T-1則會造成頻譜疊加而無法分辨。所以如果發(fā)射信號掃頻帶寬為B,則利用傅里葉變換處理被矩形窗截短的差頻信號時,LFM/CW體制測距系統(tǒng)的測距分辨力為:

這說明如果系統(tǒng)探測到的兩個目標(biāo)之間的距離差小于0.5cB-1,則不能正確分辨兩個目標(biāo)。相比在傅里葉變換下得到的結(jié)果,利用余弦變換得到的SINC函數(shù)主瓣寬度不再是T-1,而是0.5T-1。即對于初始相位為0(或近似為0)的余弦信號,在信號時長為T的情況下,其頻率分辨力為0.5T-1,比傅里葉變換提高了一倍,相應(yīng)的測距分辨力也提高到

3.2 差頻信號分析

LFM/CW體制測距系統(tǒng)得到的差頻信號由回波信號與本振信號經(jīng)混頻濾波后得到,在τ≤t≤T時間內(nèi),其形式由式(10)給出[6]。

式中,τ為回波時延,f 0為信號掃頻起始頻率,k為掃頻斜率。由式(10)可知,當(dāng)系統(tǒng)參數(shù)和目標(biāo)距離R一定時,差頻信號是一個頻率和初始相位都確定的余弦信號,其頻率為kτ,初始相位為2πf 0τ。由于

則初始相位為

分析式(12)可知,當(dāng)

差頻信號的初始相位為kπ,對于傅里葉變換和余弦變換而言,相當(dāng)于一系列初始相位為0的余弦信號,當(dāng)

差頻信號的初始相位為kπ-0.5π,相當(dāng)于一系列初始相位為0.5π的余弦信號。

4 傅里葉變換與余弦變換下LFM/CW體制測距系統(tǒng)距離分辨仿真

4.1 余弦變換下系統(tǒng)參數(shù)的設(shè)定

由前面的分析可知,傅里葉變換采用正交混頻并取模的方式避免了信號初始相位對頻譜分析的影響,而當(dāng)信號初始相位為0(或近似為0)時,使用余弦變換的方法同樣可以得到近似正確的頻譜。具體到LFM/CW體制測距系統(tǒng)中得到的差頻信號,由前面的分析可知,若目標(biāo)距離

則此時初始相位φ0將趨近于0,此時該距離范圍內(nèi)目標(biāo)對應(yīng)的一系列差頻信號的初始相位可以近似看作為0,這使得利用余弦變換處理差頻信號成為可能。若差頻信號頻率f IF和初始相位φ0滿足下面的條件

則目標(biāo)距離R的范圍可由式(17)得到

限制參數(shù)α、β可以根據(jù)具體要求的頻譜誤差和濾波器性能設(shè)定。根據(jù)這些限制參數(shù),令

便得到利用余弦變換可測量最大目標(biāo)距離

例如,若令α=0.05,β=0.7則保證了由目標(biāo)距離引入的初始相位最大不超過0.05π以及差頻信號為掃頻起始頻率 f 0的0.7倍,當(dāng)掃頻帶寬B=100 MHz,掃頻周期T=10 ms時,得到可測量最大目標(biāo)距離約為198 m。

4.2 測距分辨力仿真

為了驗證利用余弦變換提高LFM/CW體制測距系統(tǒng)測距分辨力的可行性,現(xiàn)對20 m處的多目標(biāo)進(jìn)行距離分辨,要求測距分辨力為0.25 m,掃頻周期T=1 ms。根據(jù)分辨力要求及參數(shù)限制,掃頻起始頻率應(yīng)為f 0=100 k Hz,掃頻帶寬B=300 MHz?，F(xiàn)假設(shè)三個目標(biāo)分別距離系統(tǒng)為19.75 m,20 m,20.25 m利用余弦變換得到的仿真結(jié)果如圖3所示。

圖3 利用余弦變換進(jìn)行20 m附近三目標(biāo)距離譜仿真Fig.3 The distance spectrum simulation of the 3 targets around 20 m using cosine transform

從圖中可以看出,在式(15)的條件下差頻信號的初始相位對測距造成的影響很小,余弦變換不僅給出了準(zhǔn)確的距離信息而且分辨出了相距0.25 m的三個距離。此時若使用傅里葉變換,則根據(jù)式(8)得到其測距分辨力為0.5 m,不能對三個距離進(jìn)行分辨,如圖4所示。

現(xiàn)將LFM/CW體制下利用余弦變換進(jìn)行距離分辨的方法總結(jié)如下:

1)根據(jù)測距分辨力要求確定掃頻帶寬B;

2)根據(jù)目標(biāo)距離及限制參數(shù),由式(17)確定起始掃頻頻率f0和掃頻周期T,由式(19)確定最大距離Rmax;

3)對差頻信號進(jìn)行余弦變換,完成分辨。

圖4 利用傅里葉變換進(jìn)行20 m附近三目標(biāo)距離譜仿真Fig.4 The distance spectrum simulation of the 3 targets around 20 m using Fourier transform

5 結(jié)論

本文提出了調(diào)頻連續(xù)波時頻分析中的余弦變換方法,該方法利用余弦變換在處理零初始相位單頻信號時具有高頻率分辨力這一優(yōu)點,將其應(yīng)用于處理LFM/CW體制測距系統(tǒng)得到的差頻信號從而提高測距分辨力。理論分析和仿真證實了在一定距離范圍內(nèi),使用余弦變換處理差頻信號可以得到近似正確的頻譜,相比使用傅里葉變換時其測距分辨力提高了一倍。

[1]丁鷺飛,耿富錄.雷達(dá)原理[M].西安:西安電子科技大學(xué)出版社,2006:177-182.

[2]高慶地,李世光,高正中,等.傅里葉變換的數(shù)學(xué)再認(rèn)識[J].數(shù)據(jù)采集與處理,2008,23(增刊):23-26.GAO Qingdi,LI Shiguang,GAO Zhengzhong,et al.Recognition mathematics based on fourier transform[J].Jourrnal of Data Acquisition&Processing,2008,23(增刊):23-26.

[3]Wang Z,Jillian G A,Miller W C.Recursive algorithms for the forward and inverse discretecosine transform with arbitrary length[J].IEEE Signal Processing Letter,1994,1(7):101-102.

[4]鄭君里,應(yīng)啟珩,楊為理.信號與系統(tǒng)(第二版)[M].北京:高等教育出版社,2004:109-113.

[5]郭大鈞.大學(xué)數(shù)學(xué)手冊[M].濟(jì)南:山東科學(xué)技術(shù)出版社,1985:541-542.

[6]ZHAO Jiguang,ZHANG Zhiquan,DENG Chen.Research on design and simulation of FM/cw lidar system[J].Laser Sensing and Imaging,2009,7 382:169-171.