基于單快拍數(shù)據(jù)的相干信號doa估計方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及能同時指示不同信號方向的陣列天線系統(tǒng)或方法技術(shù)領(lǐng)域，尤其涉及 一種基于單快拍數(shù)據(jù)的相干信號D0A估計方法。
【背景技術(shù)】
[0002] 以多重信號分類算法（MUSIC,multiplesignalclassification)和基于旋轉(zhuǎn)不 變技術(shù)的信號參數(shù)估計算法（ESPRIT,estimationofsignalparameterviarotational inviancetechniques)等經(jīng)典算法為代表的超分辨波達(dá)方向（DOA)技術(shù)是現(xiàn)代空間譜估 計的一個重要研宄點(diǎn)。其對信號空間到達(dá)角估計的精確性、實時性以及其廣泛的應(yīng)用前景 引起人們的極大關(guān)注。在經(jīng)典的D0A估計算法中，MUSIC算法和ESPRIT算法的實現(xiàn)都依 賴于接收數(shù)據(jù)的協(xié)方差矩陣，接收數(shù)據(jù)協(xié)方差矩陣只能通過有限次快拍數(shù)據(jù)計算來近似得 到，且要求該協(xié)方差矩陣的秩等于信號源的數(shù)目。在實際應(yīng)用中，在處理短時突發(fā)數(shù)據(jù)或接 收信號經(jīng)過相干積累后只有單快拍數(shù)據(jù)可用時，該協(xié)方差矩陣的秩降為1，則經(jīng)典的MUSIC 算法和ESPRIT算法失效。單快拍條件下的DOA估計是空間譜估計亟需解決的一個實際問 題。
[0003]目前，針對單快拍條件下的D0A估計的方法包括：直接數(shù)據(jù)域類方法、加權(quán)求和方 法和基于數(shù)據(jù)互相關(guān)預(yù)處理的相關(guān)算法等。在直接數(shù)據(jù)類算法中，大多數(shù)此類算法都只使 用奇數(shù)個陣元的接收數(shù)據(jù)構(gòu)造偽協(xié)方差矩陣，若陣元個數(shù)為偶數(shù)，則會造成接收數(shù)據(jù)信息 的浪費(fèi)；還有一些此類算法則在構(gòu)造偽協(xié)方差矩陣時增加的對信號源形式的約束，當(dāng)信號 源形式不滿足約束條件時，算法失效。加權(quán)求和方法利用部分接收數(shù)據(jù)求和后的數(shù)據(jù)構(gòu)造 偽協(xié)方差矩陣，該方法通過增加求和的信號數(shù)目來獲得信噪比的提升，同時減小了偽協(xié)方 差矩陣的自由度?；跀?shù)據(jù)互相關(guān)預(yù)處理的算法能獲得較好的D0A估計性能，但在數(shù)據(jù)預(yù) 處理方面增加了額外的計算量，且D0A估計性能受參考點(diǎn)接收數(shù)據(jù)的影響。目前單快拍條 件下的D0A估計算法中構(gòu)造的偽協(xié)方差矩陣大多都是方陣，其他形式的偽協(xié)方差矩陣尚未 得到充分應(yīng)用。

【發(fā)明內(nèi)容】

[0004] 本發(fā)明所要解決的技術(shù)問題是提供一種基于單快拍數(shù)據(jù)的相干信號D0A估計方 法，所述方法利用了單快拍下的所有接收數(shù)據(jù)，提高了D0A估計的準(zhǔn)確性。
[0005] 為解決上述技術(shù)問題，本發(fā)明所采取的技術(shù)方案是：一種基于單快拍數(shù)據(jù)的相干 信號D0A估計方法，其特征在于所述方法包括以下步驟：
[0006] 步驟一：天線的均勻直線陣列的陣元數(shù)為M，空間有N個相關(guān)性未知的遠(yuǎn)場窄帶信 號入射到所述均勻直線陣列上，則t時刻各陣元的輸出數(shù)據(jù)矩陣X(t) =A( 0 )S(t)+N(t) 為NX1的復(fù)矩陣，其中A(0) = [aplapj…a(0N)]，為MXN的陣列導(dǎo)向矢量矩陣， S(t)表示源信號矢量矩陣；N(t)表示陣列輸出的噪聲均值為零方差為〇2的加性高斯白噪 聲，且與源信號不相關(guān)；
[0007]步驟二：t時刻各陣元的輸出數(shù)據(jù)矩陣X(t) = [xjt),x2(t),???，xM(t)]T，利用每 個陣元的輸出數(shù)據(jù)構(gòu)造出偽協(xié)方差矩陣札和R2，其中札和R2的定義如下：
[0008] 當(dāng)均勻直線陣列的陣元數(shù)M為奇數(shù)時，Ri表示為
【主權(quán)項】
1. 一種基于單快拍數(shù)據(jù)的相干信號DOA估計方法，其特征在于所述方法包括以下步 驟： 步驟一：天線的均勻直線陣列的陣元數(shù)為M，空間有N個相關(guān)性未知的遠(yuǎn)場窄帶信號 入射到所述均勻直線陣列上，則t時刻各陣元的輸出數(shù)據(jù)矩陣X(t) =A( 0 )S(t)+N(t)為 NX1的復(fù)矩陣，其中A(0) = [a(01)，a(02)...a(0N)]，為MXN的陣列導(dǎo)向矢量矩陣， S(t)表示源信號矢量矩陣；N(t)表示陣列輸出的噪聲均值為零方差為〇2的加性高斯白噪 聲，且與源信號不相關(guān)； 步驟二：t時刻各陣元的輸出數(shù)據(jù)矩陣X(t) = [11(1:)，12(1:)，"％111(1:)]1，利用每個陣 元的輸出數(shù)據(jù)構(gòu)造出偽協(xié)方差矩陣札和r2，其中札和r2的定義如下： 當(dāng)均勻直線陣列的陣元數(shù)M為奇數(shù)時，Ri表示為
式中，叉為副對角線上元素全為1的方陣，維數(shù)為[(M+l)/2]X[(M+l)/2]; 當(dāng)均勻直線陣列的陣元數(shù)M為偶數(shù)時，&表示為 R2表不為
步驟三：針對所使用的天線陣列陣元的奇偶狀況，計算出步驟二所構(gòu)造出的偽協(xié)方差 矩陣札的轉(zhuǎn)置，即R/; 步驟四：構(gòu)造出新的擴(kuò)展的偽協(xié)方差矩陣R=[IVR2]; 步驟五：對R進(jìn)行奇異值分解，在奇數(shù)個天線陣元的條件下，分解出（M+1)/2個特征值 為入i彡人2彡…彡人N彡人N+1 =…=人tt+1)/2= 〇 2,在偶數(shù)個天線陣元的條件下，分解出 (M+2)/2個特征值為A#A2彡…彡AAN+1 =…=A_/2= 〇 2,通過判斷大特征值 的個數(shù)來估計信號源個數(shù)，并根據(jù)所對應(yīng)的特征向量分別得到信號子空間Us和噪聲子空 間矩陣UN;
步驟六：利用MUSIC算法構(gòu)建空間譜函； 1為源信號 的空間到達(dá)角，當(dāng)M為奇數(shù)時，^(0)表示a(0)的后（M+l)/2行，當(dāng)M為偶數(shù)時，表 示a(0)的后（M+2)/2行，使空間到達(dá)角0在（-90°，90° )范圍內(nèi)變化，找出空間譜PiuSrc( 9 )極大值點(diǎn)所對應(yīng)的角度即為源信號的D0A。
【專利摘要】本發(fā)明公開了一種基于單快拍數(shù)據(jù)的相干信號DOA估計方法，涉及能同時指示不同信號方向的陣列天線系統(tǒng)或方法技術(shù)領(lǐng)域。所述方法通過對陣列的單次快拍接收數(shù)據(jù)進(jìn)行重排，得到兩個偽協(xié)方差矩陣，然后使用這兩個偽協(xié)方差矩陣來擴(kuò)展重構(gòu)子空間算法的偽協(xié)方差矩陣。對新的偽協(xié)方差矩陣進(jìn)行奇異值分解得到信號子空間和噪聲子空間，再利用MUSIC譜估計方法對相干來波信號進(jìn)行DOA估計。本方法在單快拍條件下能夠解除源信號間相干性，同時進(jìn)一步提高了DOA的估計精度，主要應(yīng)用于單快拍條件下對相干信號的來波方向的快速估計。
【IPC分類】G01S3-74
【公開號】CN104698433
【申請?zhí)枴緾N201510114319
【發(fā)明人】鄧超升, 朱立東, 陳天宇
【申請人】電子科技大學(xué)
【公開日】2015年6月10日
【申請日】2015年3月16日