基于雙本振同步Nyguist折疊接收機(jī)的多分量LFM信號檢測及參數(shù)估計(jì)

李 睿 陳殿仁

(長春理工大學(xué)電子信息工程學(xué)院 長春 130022)

基于雙本振同步Nyguist折疊接收機(jī)的多分量LFM信號檢測及參數(shù)估計(jì)

李 睿 陳殿仁*

(長春理工大學(xué)電子信息工程學(xué)院 長春 130022)

Nyquist折疊接收機(jī)(NYFR)為一種新穎的超寬帶接收機(jī)結(jié)構(gòu)，它可以利用單片ADC實(shí)現(xiàn)超寬帶范圍內(nèi)信號的采集，而同步Nyquist折疊接收機(jī)(SNYFR)為其一種改進(jìn)結(jié)構(gòu)。該文在此基礎(chǔ)上提出一種新的基于雙本振方案的SNYFR結(jié)構(gòu)，推導(dǎo)了多分量LFM輸入下的數(shù)學(xué)模型。并基于此結(jié)構(gòu)，提出一種基于正弦調(diào)頻匹配(SFMM)與分?jǐn)?shù)階域折疊補(bǔ)償相結(jié)合的多分量LFM信號奈奎斯特區(qū)域(NZ)判定與參數(shù)估計(jì)算法，克服了單本振SNYFR結(jié)構(gòu)處理多分量LFM輸入時(shí)的困難，提高了NZ判定及參數(shù)估計(jì)的抗噪聲性能。計(jì)算機(jī)仿真表明，該文算法可以有效地對多分量LFM信號進(jìn)行檢測并實(shí)現(xiàn)各分量參數(shù)的高精度估計(jì)。

超寬帶接收機(jī)；奈奎斯特折疊接收機(jī)；正弦調(diào)頻匹配；分?jǐn)?shù)階傅里葉變換

1 引言

超寬帶接收機(jī)需要在極寬的帶寬范圍內(nèi)實(shí)現(xiàn)輸入信號的高概率接收。而受限于當(dāng)前ADC器件的發(fā)展水平，無法依靠單片高性能的ADC完成超寬帶信號的采集[1]。傳統(tǒng)的解決方案一般采用時(shí)間交替多通道并行采樣[2]，基于濾波器組的信道化數(shù)據(jù)采集[3,4]以及基于信號基展開[5]等技術(shù)實(shí)現(xiàn)超寬帶信號的采樣。受壓縮感知理論啟發(fā)，文獻(xiàn)[6,7]提出了Nyquist折疊接收機(jī)(NYguist Folding Receiver,NYFR)結(jié)構(gòu)，實(shí)現(xiàn)了單片ADC對超寬帶信號的無失真采樣。針對該結(jié)構(gòu)，文獻(xiàn)[8]利用了正交匹配追蹤(Orthogonal Matching Pursuit, OMP)、迭代貪婪展開(Iterative Greedy Unfolding, IGU)等算法恢復(fù)原始信號，但運(yùn)算量較大。文獻(xiàn)[9]采用小波分析方法處理折疊信號，簡化了計(jì)算，更易工程實(shí)現(xiàn)。文獻(xiàn)[10]提出了基于周期線性調(diào)頻(LFM)本振的同步奈奎斯特折疊接收機(jī)(SNYFR)，該結(jié)構(gòu)用數(shù)字合成的多分量周期LFM本振代替全模擬結(jié)構(gòu)產(chǎn)生的本振信號，克服了NYFR在本振調(diào)制時(shí)刻易受噪聲影響及不能提取初始相位信息等缺點(diǎn)。文獻(xiàn)[11]引入編碼機(jī)制，以調(diào)制帶寬固定但編碼信息不同的調(diào)制本振分量劃分奈奎斯特區(qū)域(Nyquist Zone, NZ)，擴(kuò)展了監(jiān)測帶寬范圍，但其基于偽魏格納分布(PWVD)的參數(shù)估計(jì)算法在處理多分量信號時(shí)存在交叉項(xiàng)。文獻(xiàn)[12]采用正弦調(diào)頻(SFM)的本振信號，給出了單分量LFM輸入下的NZ檢測及參數(shù)估計(jì)算法。然而，該算法需要先對折疊信號樣本按本振調(diào)制周期進(jìn)行再抽取，未有效利用原始樣本的所有信息，抗噪聲性能不夠理想。并且該算法在進(jìn)行折疊信號的NZ檢測時(shí)需遍歷NZ函數(shù)的所有可能取值，在輸入為多分量LFM信號時(shí)，由于NZ函數(shù)在同一時(shí)刻的多值性，該算法將不再適用。

本文為解決SNYFR處理多分量LFM信號的困難，提出了一種基于正弦調(diào)頻(SFM)本振的雙路SNYFR結(jié)構(gòu)。給出了輸入為多分量LFM信號時(shí)的信號流程。并提出了基于正弦調(diào)頻匹配(Sinusoidal Frequency Modulation Matching, SFMM)與分?jǐn)?shù)階域折疊補(bǔ)償相結(jié)合的NZ判定與參數(shù)估計(jì)算法。最后用計(jì)算機(jī)仿真驗(yàn)證了算法的正確性及有效性。

2 雙本振SNYFR及多分量LFM信號流程

雙本振SNYFR接收機(jī)結(jié)構(gòu)如圖1所示。可以看出，該結(jié)構(gòu)在原型結(jié)構(gòu)[10]基礎(chǔ)上增加了一個支路。其中，LPF1為超寬帶預(yù)選低通濾波器，用以濾除帶外噪聲，其帶寬為B。為推導(dǎo)方便，假設(shè)原始輸入為I, Q兩路相互正交的信號。則當(dāng)輸入r( t)為多分量LFM信號時(shí)有其中，lA,fl,Kl及φl分別為第l個分量的幅度、初始頻率、調(diào)頻率及初始相位。本振合成模塊生成的兩路本振信號分別為

圖1 雙本振SNYFR接收機(jī)結(jié)構(gòu)

可以看出，兩路本振信號都含有M個分量，且相鄰分量頻譜的中心頻率間隔為fs。不同的是，式(2)的各個分量都是正弦調(diào)頻(SMF)的，其調(diào)制頻偏為固定值u0，調(diào)制頻率為vm(m=1,2,…,M)。而式(3)則為M分量的單頻信號。這樣，兩路本振將頻率軸劃分成M個以fs/2+mfs為中心，帶寬為fs的奈奎斯特區(qū)域(NZ)，如圖2所示。這兩路本振信號可由本振合成模塊中的多通道DDS分別合成各個正弦及SFM分量再經(jīng)過倍頻混頻陣列[10]生成，且所有DDS都采用同一個參考信號以保證各分量的初相一致(本文為推導(dǎo)方便令初相為0)。r( t)經(jīng)功分后，分別與式(2)與式(3)混頻并經(jīng)過通帶為的復(fù)低通濾波器LPF2得到兩路折疊信號：

圖2 輸入為LFM信號的折疊過程

式(4)和式(5)中，w′(t)由輸入信號的噪聲與本振信號混頻產(chǎn)生，仍可視為白噪聲，為后面推導(dǎo)方便忽略其影響，定義ml( t)與分別為輸入信號在本振p1( t)與p2( t)下的NZ函數(shù)，表示為

若u0＜＜fs，可近似認(rèn)為ml( t)≈。圖2為頻譜跨越3個NZ的LFM分量的時(shí)頻曲線經(jīng)折疊得到s′( t)的過程示意圖。從式(4)及圖2可以看出，分布于不同NZ的輸入信號分量的頻譜都被折疊到的范圍內(nèi)，且其正弦調(diào)頻部分的調(diào)制頻率與m( t)的取值一一對應(yīng)，而m( t)的取值又與輸入分量頻譜所處的NZ相對應(yīng)。因此只要先通過調(diào)制頻率檢測m( t)的取值，即可實(shí)現(xiàn)NZ判定。再結(jié)合折疊信號分量的參數(shù)估計(jì)便可間接得到對應(yīng)的輸入LFM信號分量的參數(shù)估計(jì)值。

雖然式(4)包含輸入分量的全部信息，但由于其分量為復(fù)合調(diào)制信號，很難直接處理，在多分量的情況下則更為復(fù)雜。因此考慮引入一個輔助信號簡化NZ判定。將式(4)和式(5)經(jīng)采樣率為fs的ADC采樣后送入DSP進(jìn)行共軛相乘計(jì)算，得到輔助信號：

由式(8)可以看出，q( n)共包含L2項(xiàng)，其中L項(xiàng)為中心頻率為0，調(diào)制帶寬為2u0的正弦調(diào)頻信號，由與輸入分量rl( t)對應(yīng)的兩個折疊分量的共軛乘積產(chǎn)生，定義為自項(xiàng)。其各自的調(diào)制頻率vml ( n )與該分量對應(yīng)的離散NZ函數(shù)ml( n)一一對應(yīng)。其余L2-L項(xiàng)為不同輸入分量對應(yīng)的折疊信號間的乘積，仍然為復(fù)合調(diào)制信號。因此，通過檢測自相調(diào)制頻率vml ( n )即可間接實(shí)現(xiàn)NZ判定。考慮到雙路結(jié)構(gòu)可能出現(xiàn)的不一致性問題，還加入了幅度比較及增益控制環(huán)節(jié)，以保持雙路信號的幅度基本相同，從而保證后面算法的可實(shí)現(xiàn)性。

3 多分量LFM信號NZ判定與參數(shù)估計(jì)

3.1 基于SFMM的NZ判定

定義離散信號x( n)的N點(diǎn)SFMM為式中，u0為恒定值；0＜v≤N。則式(8)分量qk( n) 的SFMM為其中，

n2]內(nèi)取值為vmk，則當(dāng)v=vmk時(shí)，

當(dāng)v≠vmk時(shí)，根據(jù)恒等式：

由于v-vmk≠0，則式(13)中，

所以有

同理，對于Wk( u0, v)，當(dāng)v≠vmk時(shí)，

式(17)中Δflk,ΔKlk不同時(shí)為0，因式(17)中下劃線部分為一個LFM信號或一個頻率非零的單頻信號對時(shí)間求和，若Δflk或ΔKlk的值足夠大(大于算法最小分辨率)，和值將遠(yuǎn)小于SFMM產(chǎn)生的峰值，可近似視為0。由此可以得到

綜上：Qk(u0, v)只有在v等于自項(xiàng)中的某個取值時(shí)才呈現(xiàn)匹配峰值，而在其他情況下幾乎為0。由于q( n)為L個形式與式qk( n)相同分量的疊加，q( n)的SFMMQ( u0, v)將會在v=vml(n)(l=1,2, …,L )處出現(xiàn)峰值，而在其它值時(shí)近似為0。因此只要在v∈vm(m=1,2,…,M)上對Q( u0, v)進(jìn)行1維峰值搜索即可實(shí)現(xiàn)輸入分量的NZ判定。

3.2 基于分?jǐn)?shù)階域折疊補(bǔ)償?shù)膮?shù)估計(jì)算法

作為Fourier變換的廣義形式，信號的分?jǐn)?shù)階傅里葉變換(FRFT)可以解釋為將坐標(biāo)軸在時(shí)頻平面上繞原點(diǎn)以角度α做逆時(shí)針旋轉(zhuǎn)，十分適合處理LFM類信號，并且在處理多分量信號時(shí)不存在交叉項(xiàng)。

假設(shè)Q( u0, v)最大峰值處對應(yīng)NZ函數(shù)的取值為mmax，且輸入分量rk( t)與該取值對應(yīng)部分的范圍為[n1, n2]。根據(jù)mmax對應(yīng)的本振調(diào)制頻率vmmax生成本地信號：

為減小計(jì)算量，本文采用了文獻(xiàn)[13]介紹的方法對變量(α,n)進(jìn)行迭代搜索，并采用Clean算法來消除強(qiáng)弱分量及多峰值間的干擾。

4 算法仿真與分析

設(shè)接收機(jī)監(jiān)測帶寬B=10 GHz ，所有仿真的輸入信號樣本長度T=1 μs , ADC采樣頻率fs= 1 GHz。所有仿真的Monte Carlo實(shí)驗(yàn)次數(shù)為500。

4.1單分量LFM信號NZ判定性能

設(shè)輸入信號為單位幅度的單分量LFM信號，參數(shù)初始頻率及初始相位別為1.25 GHz及π/6。SFM本振調(diào)制頻偏u0=30 MHz ，調(diào)制頻率vm為從1 MHz起始，間隔為Δv =0.5 MHz的10個取值。

本文采用正確檢測概率(Probability of Correct Decision, PCD)衡量NZ判定性能，其定義為：算法正確判定NZ取值，并以此得到正確參數(shù)估計(jì)的次數(shù)與總實(shí)驗(yàn)次數(shù)的比值。圖3表示當(dāng)調(diào)頻率K分別時(shí)的PCD隨SNR的變化曲線。加入文獻(xiàn)[12]中方法的PCD作為對比可知，由于本文方法利用了所有樣本的信息因此在K較小時(shí)，PCD在SNR為-7 dB即可達(dá)到100%，抗噪聲性能明顯更優(yōu)越。而當(dāng)K較大時(shí)，輸入LFM信號的頻譜跨越多個NZ，折疊信號分段增加，以至信號能量分散，判定性能下降。考慮到實(shí)際應(yīng)用中LFM信號調(diào)制帶寬大多為百兆數(shù)量級，因此本文算法實(shí)用性更好。

4.2 多分量LFM信號NZ判定性能與本振調(diào)制參數(shù)的關(guān)系

為考察多分量LFM信號輸入下的NZ判定性能，設(shè)輸入信號中同時(shí)存在3個LFM分量，它們的幅度、初始頻率、調(diào)頻率及初始相位分別為：(1)/5}，其中分量(1)頻譜跨越3個NZ，但其頻譜能量主要分布于NZ=2。設(shè)定以保證本振分量在樣本時(shí)間內(nèi)至少存在一個完整的調(diào)制周期，令Δv =0.5 MHz。則其折疊信號的SFMM如圖4所示，可以看到，各分量對應(yīng)的NZ處都產(chǎn)生了明顯峰值，峰值大小與信號分量的幅度強(qiáng)弱相對應(yīng)。由于分量(1)跨越3個NZ，所以NZ函數(shù)在3個取樣點(diǎn)處都產(chǎn)生了峰值，其中NZ=2處峰值最大。

為考察u0取不同值時(shí)對NZ判定性能的影響，定義PCD為算法按信號分量的幅度大小依次正確判定各分量的NZ取值并得到正確參數(shù)估計(jì)的次數(shù)與試驗(yàn)次數(shù)之比。圖5為PCD與信噪比的關(guān)系(輸入信噪比以輸入信號中最弱分量(2)為標(biāo)準(zhǔn))，可以看出u0取值過小(20 MHz)或過大(150 MHz)都會造成判定性能嚴(yán)重惡化。這是由于當(dāng)u0取值過小時(shí)，式(18)中u0/vm的值變小使得SFMM峰值外衰減變得緩慢，分辨力減弱，SFMM的結(jié)果在相鄰的取值上可能出現(xiàn)錯誤的峰值。同時(shí)，過小的調(diào)制頻偏使得本振調(diào)制深度不足，抗噪聲性能減弱，也會使判定錯誤的幾率上升。當(dāng)信噪比升高時(shí)，雖然噪聲對調(diào)制的影響減弱，但是若存在峰值判定錯誤造成的誤判，噪聲的減小則反而會使出現(xiàn)錯誤的幾率上升，PCD隨信噪比的增加呈下降趨勢。而當(dāng)u0取值過大時(shí)，的近似不再成立，若信號分量的頻譜分布接近NZ邊界，則其折疊信號的NZ函數(shù)m( t)在折疊帶寬邊界處的取值可能發(fā)生錯誤的跳變，經(jīng)SFMM后也會產(chǎn)生錯誤的峰值導(dǎo)致NZ判定出錯。并且u0取值的增大會使迭代消去時(shí)的遮蔽帶寬變寬，使下一次迭代的信號能量損失增加，同樣影響判定性能。因此，建議將u0的取值設(shè)定在采樣頻率fs的3%～10%范圍內(nèi)。

4.3 參數(shù)估計(jì)性能

在得到輸入分量NZ信息后，利用3.2節(jié)提出的算法對待估計(jì)信號g( n)進(jìn)行參數(shù)估計(jì)。由于調(diào)頻率的估計(jì)精度直接影響分量各參數(shù)的估計(jì)精度，因此本文僅以調(diào)頻率估計(jì)為例討論u0的選擇對估計(jì)精度的影響。圖6給出了不同u0時(shí)，各分量調(diào)頻率估計(jì)的歸一化均方誤差(NMSE)隨信噪比的變化情況(為消除柵欄效應(yīng)采用了文獻(xiàn)[14]中介紹的插值方法)，并與克拉美–羅界(CRB)進(jìn)行比較。可以看出，由于分量(3)跨越3個NZ，能量在頻域上較分散。受NZ判定性能影響，u0的變化對其影響最明顯，u0過大或過小都造成其參數(shù)估計(jì)精度的惡化。分量(1)與分量(2)在u0的各個實(shí)驗(yàn)取值下都達(dá)到了一定的精度，但在u0取值較為合適的情況下信噪比門限值更低且估計(jì)精度更高。由此可見，輸入分量頻譜的寬度對NZ檢測性能的影響將直接反映在參數(shù)估計(jì)的性能上。當(dāng)分量頻譜存在跨NZ情況時(shí)，其帶寬越大，能量在頻域越分散，在相同的條件下，其對應(yīng)的參數(shù)估計(jì)性能也越差。

5 結(jié)束語

本文給出了一種新的雙本振SNYFR結(jié)構(gòu)，利用兩路折疊信號的相似對消性質(zhì)，提出了正弦調(diào)頻匹配(SFMM)算法，實(shí)現(xiàn)了超寬帶范圍內(nèi)多分量LFM信號的NZ判定，并且只需在預(yù)先設(shè)定的調(diào)制頻率點(diǎn)上進(jìn)行1次搜索即可完成各分量的NZ判定，計(jì)算簡便且具有較好的抗噪聲能力。同時(shí)提出的基于分?jǐn)?shù)階域折疊補(bǔ)償?shù)膮?shù)估計(jì)算法利用迭代消去的方法依次得到各分量主要參數(shù)估計(jì)值。在信噪比滿足要求的情況下具有較高的估計(jì)精度，具有一定的工程及理論參考價(jià)值。

針對算法對較弱分量存在參數(shù)估計(jì)精度下降的不足。下一步的工作為改進(jìn)算法，減小信號分量間的影響，提高參數(shù)估計(jì)精度。另外，本振信號的具體實(shí)現(xiàn)方式也需要進(jìn)一步研究。

圖3 單分量LFM信號NZ判定PCD與SNR的關(guān)系

圖4 多分量LFM信號的 SFMM在各NZ處的值

圖5 多分量LFM信號NZ判 定PCD與SNR的關(guān)系

圖6 不同u0時(shí)調(diào)頻率的估計(jì)精度隨SNR變化曲線

[1] 曾德國. 低截獲概率雷達(dá)信號偵收技術(shù)研究[D]. [博士論文],電子科技大學(xué), 2012. Zeng De-guo. Study on the interception of low probability of intercept radar signals[D]. [Ph.D. dissertation], University of Electronic Science and Technology of China, 2012.

[2] 張俊杰, 武杰, 劉尉悅, 等. 高速交替/并行數(shù)據(jù)采集系統(tǒng)時(shí)鐘研究[J]. 中國科學(xué)技術(shù)大學(xué)學(xué)報(bào), 2006, 36(3): 281-284. Zhang Jun-jie, Wu Jie, Liu Wei-yue, et al.. Clock study of high speed interleaving/multiplexing data-acquisition system [J]. Journal of University of Science and Technology of China, 2006, 36(3): 281-284.

[3] Velazquez R, Nguyen Q, Broadstone R, et al.. Design of hybrid filter banks for analog/digital conversation[J]. IEEE Transactions on Signal Processing, 1998, 46(4): 956-967.

[4] Namgoong A. Channelized digital ultra wideband receiver[J]. IEEE Transactions on Wireless Communications, 2003, 2(3): 502-510.

[5] Hoyos S, Sadler M, and Arce R. Ultra-wideband analog-todigital conversion via signal expansion[J]. IEEE Transactions on Vehicular Technology, 2005, 54(5): 1609-1622.

[6] Fudge L, Bland E, Chivers A, et al.. A Nyquist folding analog-to-information receiver[C]. 42nd Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, 2008: 541-545.

[7] Pace E, Kusmanoff A, and Fudge, L. Nyquist folding analogto-information receiver: autonomous information recovery using quadrature mirror filtering[C]. 43rd Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, 2009: 1581-1585.

[8] Ray M, Gerald L. Fudge L, et al.. Analog-to-information and the Nyquist folding receiver[J]. IEEE Journal on Emerging and Selected Topics in Circuits and Systems, 2012, 2(3): 564-578.

[9] Olusegun O, Cajetan M, Gerald F, et al.. Application of analytic wavelet transform for signal detection in Nyquist folding analog-to-information receiver[C]. IEEE International Conference on Communications, Dresden, 2009: 1-5.

[10] 曾德國, 成昊, 唐斌, 等. 基于周期LFM本振的同步Nyquist折疊接收機(jī)多分量信號參數(shù)估計(jì)算法[J]. 航空學(xué)報(bào), 2012, 33(4): 688-695. Zeng De-guo, Cheng Hao, Tang Bin, et al.. Parameter estimation approach for multicomponent signals intercepted by synchronous Nyquist folding receiver using local periodic LFM signals[J]. Acta Astronautica Sinica, 2012, 33(4): 688-695.

[11] Ke Yu-long. Intercept of frequency agility signal using coding Nyquist folding receiver[J]. WSRAS Transations on Signal Processing, 2011, 7(2): 82-91.

[12] 曾德國, 祝俊, 唐斌. 左邊帶本振SNYFR接收的LFM信號參數(shù)估計(jì)算法[J]. 宇航學(xué)報(bào), 2012, 33(6): 781-787. Zeng De-guo, Zhu Jun, and Tang Bin. Parameter estimation of LFM signal intercepted by left sideband Synchronous Nyquist Folding Receiver[J]. Journal of Astronautics, 2012, 33(6): 781-787.

[13] 齊林, 陶然, 周思勇, 等. 基于分?jǐn)?shù)階Fourier變換的多分量LFM信號的檢測和參數(shù)估計(jì)[J]. 中國科學(xué)(E輯), 2003, 33(8): 749-759. Qi Lin, Tao Ran, Zhou Si-yong, et al.. Multi-component LFM signal detection and parameter estimation based on fractional Fourier transform[J]. SCIENCE IN CHINA (Series E), 2003, 33(8): 749-759.

[14] 宋軍, 劉渝, 劉云飛. LFM信號參數(shù)估計(jì)的插值FrFT修正算法[J]. 信號處理, 2012, 28(1): 112-117. Song Jun, Liu Yu, and Liu Yun-fei. Modified interpolation parameters estimation for LFM signals based on FrFT[J]. Signal Processing, 2012, 28(1): 112-117.

李 睿： 男，1985年生，博士生，研究方向?yàn)樾盘柼幚?

陳殿仁： 男，1952年生，教授，主要從事雷達(dá)系統(tǒng)及雷達(dá)信號處理方面的研究.

Multi-component LFM Signal Detection and Parameter Estimation Algorithm Based on Synchronous Nyquist Folding Receiver with Dual Local Oscillator

Li Rui Chen Dian-ren
(School of Electronics and Information Engineering, Changchun University of Science and Technology, Changchun 130022, China)

NYquist Folding Receiver (NYFR) is a novel ultra-wideband receiver structure, which can sample a large range wideband signal by using a single ADC. The Synchronous NYquist Folding Receiver (SNYFR) is an improved structure of NYFR. In this paper, a new SNYFR structure with dual local oscillator is proposed and its mathematical model with the input of multi-component LFM signal is deduced. meanwhile, under the input of multi-component LFM signal, a Nyquist Zone (NZ) judgment and parameters estimation algorithm based on combination of Sinusoidal Frequency Modulation Matching (SFMM) and fractional domain folding compensation is proposed to overcome the difficulties of single local oscillator SNYFR structure in dealing with multi-component LFM signal and to improve noise immunity. Computer simulations show that effective detection of multicomponent LFM signal and accurate parameters estimation can be achieved by using the proposed algorithm.

UWB receiver; NYquist Folding Receiver (NYFR); Sinusoidal Frequency Modulation Matching (SFMM); FRactional Fourier Transform (FRFT)

TN971.1

A

1009-5896(2015)01-0091-06

10.11999/JEIT140281

2014-03-15收到，2014-07-16改回

*通信作者：陳殿仁 dianrenchen@cust.edu.cn