[發(fā)明專利]基于慣導(dǎo)輔助的空時(shí)域聯(lián)合抗干擾方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201310017473.8
• 申請(qǐng)日: 2013-01-17
• 公開（公告）號(hào): CN103105614A
• 公開（公告）日: 2013-05-15
• 發(fā)明（設(shè)計(jì)）人: 李彥濤
• 申請(qǐng)（專利權(quán)）人: 陜西北斗恒通信息科技有限公司
• 主分類號(hào): G01S19/21
• 分類號(hào): G01S19/21
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.基于慣導(dǎo)輔助的空時(shí)域聯(lián)合抗干擾方法，其特征在于，包括以下步驟：

步驟1、對(duì)經(jīng)過天線、射頻通道的A/D數(shù)據(jù)x(n)的固有協(xié)方差矩陣R采用分頻技術(shù)，進(jìn)行N次的累加運(yùn)算，得到Rx，令X(n)=[x₁，x₂，…，x_N]，則第m采樣點(diǎn)R(m)計(jì)算公式如下：

R(m)=x(m)^Hx(m)????(1-1)

Rx的計(jì)算可以由用公式表示為：

Rx=Σl=mm+NR(m+K(l-m))---(1-2)]]>

其中分頻的方式取決于K的取值，比如四分頻的話K＝4；

步驟2、將RX進(jìn)行放縮，用公式表示如下

其中2^N表示與Rx(1,1)最接近的整數(shù)；

步驟3、對(duì)步驟2得到的Rx'進(jìn)行矩陣求逆，該求逆的過程表示如下：

a)令Rx'可以分解為如下公式：

Rx'：LDL^H????(1-4)

其中將L設(shè)為對(duì)角線元素為1的下三角矩陣，L^H表示對(duì)L取共軛轉(zhuǎn)置，D是對(duì)角矩陣，根據(jù)(1-4)求L、D，方法采用遞推方法，遞推的法則如下表示：

令A(yù)=Rx’，m代表矩陣的每行元素個(gè)數(shù)

b)對(duì)L、D進(jìn)行矩陣求逆；

c)逆矩陣有如下公式得到：

步驟4、結(jié)合慣導(dǎo)與星歷經(jīng)過DSP方式解算得到衛(wèi)星信號(hào)的波達(dá)方向(Direction?0fArrival，DOA)，然后在用DOA根據(jù)不同天線陣列計(jì)算加權(quán)矢量W₁，W₁的計(jì)算如下：

已知得到的DOA為θ表示方位角，表示與陣列天線所在平面法線(線陣指與線陣垂面)的夾角(俯仰角)；

設(shè)d表示陣元間距，λ表示入射信號(hào)波長(zhǎng)，M表示陣元數(shù)，則對(duì)于不同的陣列，得到的加權(quán)矢量W₁不同；

對(duì)于任意陣列：其中(x_k，y_k，z_k)k＝1，2，…，M表示每個(gè)陣元在空間坐標(biāo)系的位置

再根據(jù)公式

W=(Rx′)^-1W₁/(W₁^HR^-1W₁)????(1-6)

得到最優(yōu)權(quán)值矢量W；

步驟5、根據(jù)最優(yōu)權(quán)值矢量W對(duì)每個(gè)采樣點(diǎn)的輸入信號(hào)x(n)進(jìn)行干擾消除，公式表示如下：

s(n)=W^HX(n)????(1-7)

步驟6、空域處理得到的輸出信號(hào)s(n)進(jìn)行如下操作；

設(shè)初始相關(guān)矩陣P為對(duì)角元素接近零的方陣，大小為M＝11階，對(duì)s進(jìn)行時(shí)域延拖后得到

經(jīng)過多個(gè)采樣時(shí)刻的訓(xùn)練，最優(yōu)權(quán)值w_opt訓(xùn)練成功，此時(shí)的輸出信號(hào)s_out＝(w_opt)^TS。