2013年7月川大學(xué)學(xué)報(bào)(自然科學(xué)版)第50卷第4期Journal of Sichuan University (Natural Science EditionVol 50 No 4doi:103969iss.0490-6756.2013.04.022基于FFT的空間重采樣算法的性能分析賈巨庫,何培宇,張勇(四川大學(xué)電子信息學(xué)院,成都610065)摘要:空間重釆樣方法是一種不需要角度預(yù)估計(jì)的寬帶DOA估計(jì)方法,它利用相位聚焦的方法實(shí)現(xiàn)寬帶信號(hào)中不同信號(hào)子帶的信號(hào)子空間的一致.而基于FFT變換的空間重采樣算法利用傅里葉變換特性實(shí)現(xiàn)空間采樣頻率的變換,但此方法對(duì)陣元的選擇以及參考頻率的選擇等都有一定的要求,本文從陣元個(gè)數(shù)、參考頻率、子帶劃分、信號(hào)入射角度以及信號(hào)帶寬等多個(gè)方面對(duì)基于FFT空間重采樣算法進(jìn)行了全面細(xì)致的性能仿真分析,為該方法的應(yīng)用推廣、工程實(shí)現(xiàn)提供了較為有意義的參考關(guān)鍵詞:寬帶DOA;傅里葉變換;空間重采樣;性能分析中圖分類號(hào):TN912文獻(xiàn)標(biāo)識(shí)碼:A文章編號(hào):0490-6756(2013)4-0781-06The performance analysis of spatial resampling based on FFTJIA Ju-Ku, HE Pei-Yu, ZHANG Yong(School of Electronics and Information Sichuan U niversity Chengdu 610065, ChinaAbstract: Spatial resampling is a wideband direction of arrival (DOA )estimation method with no need topre-estimate the angle of incidence. It uses phase focusing to uniformize the subspaces of wideband signals sub-bands. The spatial resampling method based on FFT uses the Fourier transform to adjust thespatial resampling rate. Yet, the performance of the spatial resampling method based on FFT largely depends on such factors as the amount of array elements choice of the reference frequency, etc. A detailedperformance analy sis of the spatial resampling method based on FFT is given. This analysis is givenfrom multiple aspect, e.g. number of array elements, reference frequency, sub-band division, angles ofincidence, signal bandw idth, etc, w hich makes it a valuable reference for follow-up study and engineering implementationKey words wide-band DOA, fourier transform, spatial resampling performance analy sis引言的預(yù)估計(jì),且最終測向結(jié)果很大程度上受到預(yù)估計(jì)結(jié)果準(zhǔn)確性的影響目前,對(duì)寬帶信號(hào)的DOA估計(jì),主要分為有為了避免因初始方位角的估計(jì)而引入誤差,文兩類方法:基于非相干信號(hào)的處理方法(ISM)和基獻(xiàn)[1提出了廣義陣列流型內(nèi)插算法,引入了指向于相干信號(hào)處理方法(CSM).ISM算法主要問題矢量的貝塞爾函數(shù)近似表達(dá)式,從而將方向矩陣變在于計(jì)算量大,無法估計(jì)相干信號(hào)源CsM算法可換中國煤化工有關(guān),另一部分與頻率分辨相干信號(hào)源,分辨率較高,但需要進(jìn)行方向角有HCNMHG函數(shù)表示導(dǎo)向矢量時(shí)收稿日期:2012-06-11基金項(xiàng)目:國家自然科學(xué)基金(6作者簡介:賈巨庫(1987—),男,河北衡水,2010級(jí)碩士研究生,主要研究方向?yàn)樾盘?hào)處理理論通訊作者:何培宇,E-mlil:spysbsya163.com四川大學(xué)學(xué)報(bào)(自然科學(xué)版)第50卷引入了誤差,而且要獲得聚焦矩陣需要在每個(gè)子帶a(內(nèi)進(jìn)行特征值分解,因此運(yùn)算量很大文獻(xiàn)2道首先1cxp(-jgn)…cxp(j(M-1)gn)提出將一組陣元輸出數(shù)據(jù)看做是空間信號(hào)的離散k=1,2…,K;p=1,2…,P采樣值,以此作為空間頻率的概念,信號(hào)通過頻率映射濾波器達(dá)到空間頻率的變換,即陣元間距的變式中,gn=2ml/sin/e,4為第p個(gè)目標(biāo)的入射換,從而得到虛擬陣列輸出角文獻(xiàn)[3]同樣應(yīng)用此概念提岀了插值矩陣的方對(duì)于寬帶信號(hào),在信號(hào)與噪聲相互獨(dú)立的條件法文獻(xiàn)[4利用FFT變換實(shí)現(xiàn)虛擬陣元間距的變下,子帶f的協(xié)方差矩陣表示為R()滿足換,文獻(xiàn)[5對(duì)文獻(xiàn)[2]和[4中的兩種方法進(jìn)行了R(A)=E X(X"(R)H比較,得出基于FFT變換的空間重采樣算法更優(yōu)A(fuR(fu)A"(f*)+ Ro(fa的結(jié)論從目前的文獻(xiàn)來看還沒有對(duì)基于FFT的式中,R(fA)=ES(f)S"(f),R()為噪聲空間重采樣算法較全面而細(xì)致的性能分析,本文從協(xié)方差矩陣參考頻率、陣元個(gè)數(shù)、劃分子帶個(gè)數(shù)、信號(hào)入射角度以及信號(hào)帶寬等多個(gè)方面對(duì)基于FT空間重采樣3空間重采樣算法進(jìn)行了全面細(xì)致的分析,以期為該方法的推廣由式(1)和式(2)可以看到,方向矩陣不僅與信和工程應(yīng)用提供較為有意義的參考號(hào)的入射角度aaG有關(guān),而且與信號(hào)的2信號(hào)模型頻率有關(guān),那么不同子帶的陣列流型隨頻率的變化而變化,這就造成了寬帶信號(hào)各子帶的信號(hào)子空間由M個(gè)陣元組成的均勻線陣,各陣元間距為不同,對(duì)單獨(dú)一個(gè)信號(hào)而言,其陣列流型中的第md,為避免相位模糊,d為入射信號(hào)最高頻率所對(duì)行第k列元素為應(yīng)波長的一半空間中有P個(gè)遠(yuǎn)場寬帶信號(hào),占有exp((j(n-1)2dfisin0/c))相同的帶寬,其入射方向分別為a,a…,鼻,則第(4)m個(gè)陣元上接收到的信號(hào)可表示為由式(4)可以看到,要保持不同頻率處的陣列流型tn(t)s(t+(m-1) dsin a/c)+m(t)相同,需滿足dfk=C(C為常數(shù)2…,M;p=1,2…,P這樣,在不同頻率處可以得到相同的陣列流(1)型式中,s(t)為p方向入射的寬帶信號(hào),ma(t)為加通常我們選用既定的某一參考頻率f及其實(shí)性高斯噪聲,c為波的傳播速度際物理陣元間距d的乘積作為常數(shù)C,即若將寬帶信號(hào)劃分為K個(gè)子帶,則其f處的dfk= fo do= c頻域表達(dá)式則與f頻率對(duì)應(yīng)的陣元間距可以表示為對(duì)應(yīng)頻率tm(ft)的函數(shù)S,(fk )exp(j2 Tf(m-1)dsin B/c)+ nm(L)d(fa)=d fo/fs(7)k=1,2…,K1,2……,MP0=f0/f由式(7)可得:為了使得不同的頻率子帶對(duì)應(yīng)則陣列的輸出矩陣形式可表示為相同的信號(hào)子空間,當(dāng)子帶的中心頻率壓縮p倍X(f)=x(f)x(f)…xw(f)時(shí),所對(duì)應(yīng)的陣列間距應(yīng)擴(kuò)展p倍A(&S(f)+N(f)對(duì)于均勻線陣作Y中國煤化子果排可以將A(f)=[a(f)a(A)…a(n)同1不則,到的一組陣元數(shù)據(jù)看(f)=(f)s(f)P(f4)作對(duì)入射信號(hào)進(jìn)行了多次時(shí)域采樣.因此可以應(yīng)用A(/)是M×P維的方向矩陣,其中元素表示為信號(hào)的時(shí)域處理方法處理這組陣列接收數(shù)據(jù)4期賈巨庫,等:基于FFT的空間重采樣算法的性能分析b)a>1時(shí),D= Round(M/a)+M,對(duì)X:()插零個(gè)數(shù)為(D-M)×N點(diǎn)其中 Round(·)表示取最接近的整數(shù),4)獲得虛擬陣的采樣輸出a)p≤1時(shí),對(duì)插零后的X(f)進(jìn)行D點(diǎn)的逆FFT變換,直接得到虛擬輸出X(f)圖1陣列信號(hào)接收示意圖b)p>1時(shí),對(duì)插零后的X(f)進(jìn)行D點(diǎn)Fig 1 Schematic of receiving array signal的逆FFT變換,并按照( Round(p)+1)間隔取行數(shù)據(jù)得到X(f);傅里葉變換性質(zhì):然后,取X(f)的前D行作為虛擬陣的采樣x(d) X(u)(8)輸出X.(f)5)利用虛擬陣列輸出,估計(jì)聚焦后的采樣協(xié)信息擴(kuò)展p倍同樣,也可以通過變換信號(hào)的頻域方差矩陣信息來實(shí)現(xiàn)其時(shí)域信號(hào)采樣率的偽變化R A2X,(F)X,(fi)這里,可將陣元接收的一組快拍數(shù)據(jù)看作入射信號(hào)的短時(shí)序列,傅里葉變換到頻域,在頻域?qū)ζ?)利用窄帶的 MUSIC算法進(jìn)行DOA估計(jì)進(jìn)行補(bǔ)零擴(kuò)展,傅里葉反變換后,實(shí)現(xiàn)了信號(hào)空間4實(shí)驗(yàn)仿真采樣率的偽變化,即虛擬陣元間距的變化,如圖2所示實(shí)驗(yàn)一信號(hào)入射角度對(duì)算法影響實(shí)驗(yàn)條件:接收陣列是一個(gè)15陣元的均勻線陣,陣元間距為信號(hào)最高頻率對(duì)應(yīng)波長的一半,人射信號(hào)為相干信號(hào),中心頻率為2000Hz,相對(duì)帶寬50%,快拍數(shù)N=200,劃分子帶21個(gè),空間重采樣算法的參考頻率選擇為信號(hào)的最低頻率1500Hz,重復(fù)獨(dú)立試驗(yàn)次數(shù)為30次,采樣率為40kHz噪聲為高斯白噪聲,信噪比SNR=10dB圖2空間重采樣變換后陣列示意Schematic of array af ter spatial resampling將寬帶信號(hào)劃分為K個(gè)互不重疊的子帶,每個(gè)子帶滿足窄帶信號(hào)的假設(shè),子帶中心頻率為f(k=1,2…,K),參考頻率為∫m(一般設(shè)為信號(hào)的最低頻率),為了獲得各個(gè)子帶信號(hào)子空間的致需要的陣列輸出個(gè)數(shù)D=Mf1/fm(f為信號(hào)的最低頻率),對(duì)于第k個(gè)子帶,陣列的輸出為X(f)=[X1(f),X2(f)…,X(f),那么利用FFT插值的寬帶信號(hào)方位估計(jì)步驟如下:-50-40-30-20-10010203040501)對(duì)每個(gè)子帶的采樣輸出X(f)做M點(diǎn)的時(shí)角度(°)FFT變換,得到X()中國煤化工2)計(jì)算當(dāng)前頻率與參考頻率的對(duì)應(yīng)關(guān)系:aCNMHG重采樣算法的影響g 3 Incident angles impact on spatial resampling3)對(duì)X:()進(jìn)行插零由圖3可知,沿陣元法線方向入射信號(hào)的估計(jì)a)p≤1時(shí),D= Round(M/),對(duì)x()誤差最小,隨偏離法線程度的增大,計(jì)誤差逐漸插零個(gè)數(shù)為(D-M)XN點(diǎn)784四川大學(xué)學(xué)報(bào)(自然科學(xué)版)第50卷變大,當(dāng)超過40°時(shí),估計(jì)誤差急劇增大,因此此算法的估計(jì)角度范圍較窄,對(duì)其工程應(yīng)用范圍產(chǎn)生定限制實(shí)驗(yàn)二參考頻率對(duì)算法的影響實(shí)驗(yàn)條件:接收陣列是一個(gè)15陣元的均勻線陣,陣元間距為信號(hào)最高頻率對(duì)應(yīng)波長的一半,入射信號(hào)為相干信號(hào),中心頻率為2000Hz,相對(duì)帶寬50%,入射角度分別為10°、14°,劃分子帶21個(gè),快拍數(shù)N=200,重復(fù)獨(dú)立試驗(yàn)次數(shù)為30次,采樣率為40kHz,噪聲為高斯白噪聲,信噪比SNR=10dB.空間重采樣算法采用信號(hào)最低頻1500170019002100率作為參考頻率參考頻率(Hz)圖4中平均誤差值為同一參考頻率下的三個(gè)圖4參考頻率對(duì)空間重采樣算法的影響入射信號(hào)估計(jì)角度的誤差平均值,由圖4可知,參ig. 4 Reference frequency's impact on spatial考頻率為信號(hào)最低頻率1500Hz時(shí),均方誤差只有amp20dB,算法在參考頻率為1700Hz時(shí),估計(jì)誤差最小,由第二部分的操作步驟可知,虛擬陣元個(gè)數(shù)D= Round(M//f),虛擬陣元D的估計(jì)與參考頻率f有關(guān),當(dāng)f/f為整數(shù)時(shí),不存在D的近似估計(jì),所以參考頻率f的選取對(duì)于估計(jì)的準(zhǔn)確性有一定的影響.同時(shí),在不以信號(hào)最低頻率作為參考頻率fm時(shí),則p>1,如果入射信號(hào)帶寬過大,且算法不以最低頻率作為參考頻率時(shí),在進(jìn)行操作步驟中的第4步時(shí)就可能出現(xiàn)只能獲取單個(gè)采樣輸出的現(xiàn)象,即只獲取一個(gè)虛擬陣元數(shù)據(jù),就無法進(jìn)行DOA估計(jì),因此選取高頻作為算法的參考頻率時(shí),會(huì)對(duì)入射信號(hào)的帶寬有一定的要求,0.4應(yīng)用范圍變窄相對(duì)帶寬實(shí)驗(yàn)三信號(hào)帶寬對(duì)算法的影響圖5信號(hào)帶寬對(duì)空間重采樣算法的影響實(shí)驗(yàn)條件:接收陣列是一個(gè)15陣元的均勻線Fig 5 Signal bandw idths impact on spatial resampling陣,陣元間距為信號(hào)最高頻率對(duì)應(yīng)波長的一半,入射信號(hào)為相干信號(hào),中心頻率為2000Hz,角度分實(shí)驗(yàn)條件:接收陣列是一個(gè)15陣元的均勻線別為-20°、-10°、14°,快拍數(shù)N=200,劃分出的陣,陣元間距為信號(hào)最高頻率對(duì)應(yīng)波長的一半,入子帶頻率間隔為100Hz,空間重采樣算法的參考射信號(hào)為相干信號(hào),中心頻率為2000Hz,相對(duì)帶頻率選擇為1500Hz,重復(fù)獨(dú)立試驗(yàn)次數(shù)為30次,寬50%6,度分別為-20°、-10°、14°,快拍數(shù)N采樣率為40kHz,噪聲為高斯白噪聲,信噪比200,重復(fù)獨(dú)立試驗(yàn)次數(shù)為30次,采樣率為40SNR=IO dBkHz,噪聲為高斯白噪聲,信噪比SNR=10dB.空圖5中,均方誤差值為同一帶寬下3個(gè)信號(hào)估間平田位是征頻率作為參考頻率計(jì)角度均方誤差的平均值,由圖5可知,在相對(duì)帶中國煤化工帶子帶劃分個(gè)數(shù)下寬為0.1時(shí),均方誤差只有—30dB左右,相對(duì)帶3CNMH的平均值,由圖6可寬為06處,算法取得最佳估計(jì)效果,基于FFT空知,相同帶寬條件下,在子帶數(shù)為10與50處算法間重采樣算法的估計(jì)效果隨著帶寬變大而變好,當(dāng)取得最好估計(jì)效果,子帶數(shù)為20與30處,效果較相對(duì)帶寬超過一定界限時(shí),估計(jì)效果開始變差差,由第二部分的操作步驟可知,虛擬陣元個(gè)數(shù)Dk實(shí)驗(yàn)四子帶劃分對(duì)算法的影響= Round(Mf/f),從仿真條件可計(jì)算得,參考頻4期賈巨庫,等:基于FFT的空間重采樣算法的性能分析率f=1500Hz,頻帶寬度為1000Hz,當(dāng)子帶數(shù)f、fm確定后,陣元個(gè)數(shù)M變換導(dǎo)致了估計(jì)虛擬為10時(shí),Mf/f計(jì)算結(jié)果均為整數(shù),即計(jì)算虛擬陣元D時(shí)的誤差不同,所以性能存在差別,因此算陣元個(gè)數(shù)D時(shí)不存在近似計(jì)算問題,當(dāng)子帶個(gè)數(shù)法對(duì)陣元個(gè)數(shù)也有一定的要求,并非陣元個(gè)數(shù)越多分別為20、30、40、50時(shí),Mf/f計(jì)算所得結(jié)果中越好有小數(shù)存在,因此需要對(duì)其進(jìn)行取整處理,這就造成了D的估計(jì)誤差不同,在子帶數(shù)為20處取整時(shí)造成誤差較大,最終DOA估計(jì)效果較差,而子帶數(shù)為50處時(shí)取整造成誤差較小,則所得DOA估計(jì)效果也較好.因此,當(dāng)陣元個(gè)數(shù)以及參考頻率確定時(shí),虛擬陣元D估計(jì)誤差的大小與各子帶中心頻率f的選取有關(guān),當(dāng)劃分子帶個(gè)數(shù)合理時(shí),能有效減小甚至消除虛擬陣元D的估計(jì)誤差,從而提-100高DOA的估計(jì)性能140陣元個(gè)數(shù)圖7陣元個(gè)數(shù)對(duì)空間重采樣算法的影響Fig 7 Number of array eler4結(jié)論基于FFT的空間重采樣算法與CSM算法的思想相同,均為利用聚焦思想將各子帶的信號(hào)子空0.20.40.6相對(duì)帶寬間對(duì)齊后疊加,因此此算法保持了解相干的能力,同時(shí)也不需要進(jìn)行角度的預(yù)估計(jì).空間重采樣算法圖6子帶劃分個(gè)數(shù)對(duì)空間重采樣算法的影響的有效估計(jì)范圍在-40°~40°之間,估計(jì)范圍較Fig 6 Sub-band divisions impact on spatial resampli窄,其工程應(yīng)用范圍產(chǎn)生一定限制.因?yàn)樵肼暈楠?dú)實(shí)驗(yàn)五陣元個(gè)數(shù)對(duì)算法的影響立的,利用時(shí)頻域變換的方法進(jìn)行原始信號(hào)的插值實(shí)驗(yàn)條件:接收陣列為均勻線陣,陣元間距為估計(jì)不可避免的出現(xiàn)誤差,同時(shí)由D=信號(hào)最高頻率對(duì)應(yīng)波長的一半,入射信號(hào)為相干信Round(Mf/f)可以看出,插零個(gè)數(shù)與陣元個(gè)數(shù)、號(hào),中心頻率為200Hz,相對(duì)帶寬50%,角度分別參考頻率以及子帶中心頻率有關(guān),且需要近似取整為-20°、-10°快拍數(shù)N=20,劃分頻率子帶數(shù),這樣也不可避免的引入誤差因此基于FFT空21個(gè),空間重采樣算法的參考頻率選擇為150間重采樣算法對(duì)陣元個(gè)數(shù)、參考頻率、子帶劃分個(gè)數(shù)等都有一定的要求,在選擇陣元個(gè)數(shù)、參考頻率Hz,重復(fù)獨(dú)立試驗(yàn)次數(shù)為30次,采樣率為40kHz,及子帶劃分?jǐn)?shù)時(shí)盡量保證計(jì)算出的虛擬陣元數(shù)噪聲為高斯白噪聲,信噪比SNR=10dB圖7中均方誤差值為同一陣元個(gè)數(shù)下3個(gè)信M為整數(shù),進(jìn)而減小誤差,在保證以上條件情況號(hào)估計(jì)角度均方誤差的平均值,由圖7可知,在低下,此算法能取得較好的估計(jì)效果陣元個(gè)數(shù)時(shí)算法估計(jì)效果較差,隨著陣元個(gè)數(shù)增參中國煤化工加,估計(jì)效果逐漸變好,且均方誤差基本保持在CNMHG[1 Doron MA, Doron E. Weiss A J Coherent wide-80dB左右,而在陣元數(shù)為16時(shí),估計(jì)效果達(dá)到最rayJ好從計(jì)算虛擬陣元個(gè)數(shù)D= Round(Mf/f)可IEEE Trans. Signal Processing, 1993, 41: 414看出,虛擬陣元D的估計(jì)不僅與各子帶頻率∫、[2] Krolik j, Swingler d. Focused wide-band array pro-參考頻率∫m有關(guān),也與實(shí)際陣元個(gè)數(shù)M有關(guān)cessing by spatial resampling [J]. IEEE Trans AS四川大學(xué)學(xué)報(bào)(自然科學(xué)版)第50卷SP,1990,38(2):33 Friedlander B, Weiss A J. Direction finding for wide[71 Debasis Kundu Modified M USIC algorithm for estiband signals using an interpolated array [J. IEEEating dOa of signals [J]. IEEE Signal ProcessingTrans Signal Processing, 1993, 41(4): 16181996,48:85[4]朱維杰,孫進(jìn)才.基于快速富式變換插值的寬帶信號(hào)[8] Prasad K p, Satyanarayana P. Fast interpolation al-方位估計(jì)[J],聲學(xué)學(xué)報(bào),2002,27(6):513orithm using FFT [J. Electronics Letters, 1986[5]田野,趙春暉.兩種實(shí)現(xiàn)寬帶信號(hào)空間重采樣測向算法的比較[J]現(xiàn)代電子技術(shù),2008,272(9):1456]趙春暉,黃光亞,李剛.基于FFT變換的空間重采樣[責(zé)任編輯:伍少梅]寬帶Root- Music算法].應(yīng)用科技,2007,32(3):中國煤化工CNMHG