本發(fā)明涉及通信領(lǐng)域，尤指一種多星約束穩(wěn)健空頻抗干擾處理方法和裝置。

背景技術(shù)：

為了在復(fù)雜電磁環(huán)境中抑制掉強(qiáng)壓制式干擾，陣列天線波束形成技術(shù)廣泛應(yīng)用于擴(kuò)頻通信系統(tǒng)。然而自適應(yīng)波束形成技術(shù)對(duì)于期望信號(hào)導(dǎo)向矢量的失配非常敏感，由于存在這種失配會(huì)使陣列波束形成器把期望信號(hào)進(jìn)行一定程度抑制，因此被稱為信號(hào)的自相消。

傳統(tǒng)的波束形成算法不能改善未知傳感器互耦、近場(chǎng)波前失真、源擴(kuò)展以及相干和非相干的局部散射等影響，因此穩(wěn)健的波束形成技術(shù)越來(lái)越成為波束形成研究中的重點(diǎn)和熱點(diǎn)。

技術(shù)實(shí)現(xiàn)要素：

為了解決上述技術(shù)問(wèn)題，本發(fā)明提供了一種多星約束穩(wěn)健空頻抗干擾處理方法和裝置，可以解決由于期望信號(hào)導(dǎo)向矢量的失配等原因引起的輸出信干噪比惡化問(wèn)題，增強(qiáng)了算法穩(wěn)健性。

為了達(dá)到本發(fā)明目的，本發(fā)明提供了一種多星約束穩(wěn)健空頻抗干擾處理方法，包括：

將接收到的時(shí)域信號(hào)通過(guò)傅立葉變換，轉(zhuǎn)換為頻域數(shù)據(jù)；

根據(jù)所述頻域數(shù)據(jù)計(jì)算信號(hào)帶寬內(nèi)頻點(diǎn)的協(xié)方差矩陣，確定協(xié)方差矩陣修正因子，根據(jù)所述協(xié)方差矩陣修正因子計(jì)算權(quán)向量；

使用所述權(quán)向量對(duì)所述頻域數(shù)據(jù)進(jìn)行濾波處理；

對(duì)濾波后的頻域數(shù)據(jù)進(jìn)行傅里葉逆變換，轉(zhuǎn)換為時(shí)域數(shù)據(jù)并輸出。

可選地，所述根據(jù)所述頻域數(shù)據(jù)計(jì)算信號(hào)帶寬內(nèi)頻點(diǎn)的協(xié)方差矩陣，確定協(xié)方差矩陣修正因子，根據(jù)所述協(xié)方差矩陣修正因子計(jì)算權(quán)向量的步驟中，

對(duì)頻域數(shù)據(jù)采用分段處理的方式，通過(guò)設(shè)置所述權(quán)向量中部分矢量為零，以將期望信號(hào)帶外部分進(jìn)行清零處理。

可選地，所述確定協(xié)方差矩陣修正因子的步驟中，

所述協(xié)方差矩陣修正因子α/β通過(guò)下式計(jì)算得到：

α/β＝γtr(R_k)/M，γ＝1-e^-0.1M

其中，R_k為所述協(xié)方差矩陣，M為所述協(xié)方差矩陣的維數(shù)，tr()表示矩陣的跡。

可選地，所述方法還包括：

利用先驗(yàn)信息估計(jì)衛(wèi)星波達(dá)方向，確定空頻導(dǎo)向矢量。

可選地，所述根據(jù)所述協(xié)方差矩陣修正因子計(jì)算權(quán)向量的步驟中，分別對(duì)每個(gè)頻點(diǎn)的權(quán)向量進(jìn)行計(jì)算，所述權(quán)向量w(f_k)通過(guò)下式計(jì)算得到：

其中，f_k表示第k個(gè)頻點(diǎn)，為樣本協(xié)方差矩陣，α/β為協(xié)方差矩陣修正因子，I為單位陣，A為空頻導(dǎo)向矢量，b為N×1維的全1列向量，N為所述時(shí)域信號(hào)的個(gè)數(shù)。

本發(fā)明還提供一種多星約束穩(wěn)健空頻抗干擾處理裝置，應(yīng)用于空頻自適應(yīng)處理器，包括：

第一信號(hào)轉(zhuǎn)換模塊，用于將接收到的時(shí)域信號(hào)通過(guò)傅立葉變換，轉(zhuǎn)換為頻域數(shù)據(jù)；

計(jì)算模塊，用于根據(jù)所述頻域數(shù)據(jù)計(jì)算信號(hào)帶寬內(nèi)頻點(diǎn)的協(xié)方差矩陣，確定協(xié)方差矩陣修正因子，根據(jù)所述協(xié)方差矩陣修正因子計(jì)算權(quán)向量；

濾波模塊，用于使用所述權(quán)向量對(duì)所述頻域數(shù)據(jù)進(jìn)行濾波處理；

第二信號(hào)轉(zhuǎn)換模塊，用于對(duì)濾波后的頻域數(shù)據(jù)進(jìn)行傅里葉逆變換，轉(zhuǎn)換為時(shí)域數(shù)據(jù)并輸出。

可選地，所述計(jì)算模塊，進(jìn)一步用于對(duì)頻域數(shù)據(jù)采用分段處理的方式，通過(guò)設(shè)置所述權(quán)向量中部分矢量為零，以將期望信號(hào)帶外部分進(jìn)行清零處理。

可選地，所述計(jì)算模塊，進(jìn)一步用于通過(guò)下式計(jì)算得到所述協(xié)方差矩陣修正因子α/β：

α/β＝γtr(R_k)/M，γ＝1-e^-0.1M

其中，R_k為所述協(xié)方差矩陣，M為所述協(xié)方差矩陣的維數(shù)，tr()表示矩陣的跡。

可選地，所述計(jì)算模塊，還用于利用先驗(yàn)信息估計(jì)衛(wèi)星波達(dá)方向，確定空頻導(dǎo)向矢量。

可選地，所述計(jì)算模塊，進(jìn)一步用于分別對(duì)每個(gè)頻點(diǎn)的權(quán)向量進(jìn)行計(jì)算，通過(guò)下式計(jì)算得到權(quán)向量w(f_k)：

其中，f_k表示第k個(gè)頻點(diǎn)，為樣本協(xié)方差矩陣，α/β為協(xié)方差矩陣修正因子，I為單位陣，A為空頻導(dǎo)向矢量，b為N×1維的全1列向量，N為所述時(shí)域信號(hào)的個(gè)數(shù)。

與現(xiàn)有技術(shù)相比，本發(fā)明實(shí)施例通過(guò)采用協(xié)方差矩陣修正因子計(jì)算權(quán)向量，可以保證多顆衛(wèi)星信號(hào)數(shù)據(jù)失配時(shí)仍能有效地抑制寬帶干擾，增強(qiáng)了空域?yàn)V波算法的穩(wěn)健性。進(jìn)一步地，通過(guò)對(duì)帶外信號(hào)進(jìn)行清零處理，在保證濾波效果的前提下，簡(jiǎn)化了工程運(yùn)算量。

附圖說(shuō)明

附圖用來(lái)提供對(duì)本發(fā)明技術(shù)方案的進(jìn)一步理解，并且構(gòu)成說(shuō)明書的一部分，與本申請(qǐng)的實(shí)施例一起用于解釋本發(fā)明的技術(shù)方案，并不構(gòu)成對(duì)本發(fā)明技術(shù)方案的限制。

圖1為本發(fā)明實(shí)施例的多星約束穩(wěn)健空頻抗干擾處理方法流程圖；

圖2為本發(fā)明應(yīng)用示例的示意圖；

圖3為本發(fā)明實(shí)施例的多星約束穩(wěn)健空頻抗干擾處理裝置示意圖；

圖4為本發(fā)明應(yīng)用示例的空頻自適應(yīng)處理器內(nèi)部處理過(guò)程示意圖。

具體實(shí)施方式

為使本發(fā)明的目的、技術(shù)方案和優(yōu)點(diǎn)更加清楚明白，下文中將結(jié)合附圖對(duì)本發(fā)明的實(shí)施例進(jìn)行詳細(xì)說(shuō)明。需要說(shuō)明的是，在不沖突的情況下，本申請(qǐng)中的實(shí)施例及實(shí)施例中的特征可以相互任意組合。

在附圖的流程圖示出的步驟可以在諸如一組計(jì)算機(jī)可執(zhí)行指令的計(jì)算機(jī)系統(tǒng)中執(zhí)行。并且，雖然在流程圖中示出了邏輯順序，但是在某些情況下，可以以不同于此處的順序執(zhí)行所示出或描述的步驟。

如圖1所示，本發(fā)明實(shí)施例的多星約束穩(wěn)健空頻抗干擾處理方法包括如下步驟：

步驟101，將接收到的時(shí)域信號(hào)通過(guò)傅立葉變換，轉(zhuǎn)換為頻域數(shù)據(jù)；

所述傅立葉變換可以是快速傅立葉變換(FFT)或者離散傅立葉變換(DFT)。

對(duì)于寬帶接收機(jī)，考慮N個(gè)遠(yuǎn)場(chǎng)的寬帶信號(hào)s_i,i＝1,2,...N(包括期望信號(hào)和干擾)以θ_i,i＝1,2,...N方向入射到空間M元均勻線陣上，通常假設(shè)第1項(xiàng)s₁對(duì)應(yīng)于感興趣的期望信號(hào)，而s_i,i＝2,...N對(duì)應(yīng)于N-1個(gè)干擾信號(hào)。陣列單元間距為接收信號(hào)最高頻率對(duì)應(yīng)的波長(zhǎng)的一半。假設(shè)信號(hào)的帶寬為B，則第m個(gè)陣元的接收數(shù)據(jù)

其中τ_mi表示第m個(gè)陣元上接收的第i個(gè)信號(hào)的延時(shí)，n_m(t)為第m個(gè)陣元通道的熱噪聲。

在步驟101之前，本發(fā)明實(shí)施例還可包括：利用先驗(yàn)信息估計(jì)衛(wèi)星波達(dá)方向，確定空頻導(dǎo)向矢量。

其中，利用先驗(yàn)信息估計(jì)當(dāng)前N顆可見(jiàn)衛(wèi)星的波達(dá)方向θ_i,i＝1,2,...N，求解出對(duì)應(yīng)空頻導(dǎo)向矢量A(f_k)＝[a(θ₁,f_k)a(θ₂,f_k)...a(θ_N,f_k)]，以均勻線陣為例，波達(dá)方向?yàn)棣?sub>1的衛(wèi)星信號(hào)空頻導(dǎo)向矢量為

其中d為陣元間距，M為陣元數(shù)，c為光速常數(shù)。

如圖2所示，采用M陣元的空頻二維處理器結(jié)構(gòu)，對(duì)各陣元中頻AD(模擬信號(hào)轉(zhuǎn)換為數(shù)字信號(hào))采樣信號(hào)進(jìn)行K點(diǎn)傅立葉變換變換得到頻域數(shù)組x₁₁,x₁₂,…x_1K；x₂₁,x₂₂,…x_2K；…x_M1,x_M2,…x_MK，數(shù)據(jù)中包含衛(wèi)星信號(hào)、干擾以及噪聲。

步驟102，根據(jù)所述頻域數(shù)據(jù)計(jì)算信號(hào)帶寬內(nèi)頻點(diǎn)的協(xié)方差矩陣，確定協(xié)方差矩陣修正因子，根據(jù)所述協(xié)方差矩陣修正因子計(jì)算權(quán)向量；

采樣數(shù)據(jù)經(jīng)過(guò)K點(diǎn)的傅立葉變換后，有1＜B₁＜B₂＜K，B₁、B₂、K均為正整數(shù)，其中區(qū)間[B₁,B₂]為期望信號(hào)帶寬。為了對(duì)數(shù)據(jù)進(jìn)行有效的處理，同時(shí)降低計(jì)算復(fù)雜度，可以對(duì)頻域數(shù)據(jù)采用分段處理的方式，對(duì)于期望信號(hào)帶外部分進(jìn)行清零處理，減少不必要的計(jì)算工作。

即：對(duì)于頻點(diǎn)f_k，k∈[1,B₁)or(B₂,K]的頻域數(shù)據(jù)，有w(f_k)＝O_M×1，O_M×1為M×1維的零矢量。

對(duì)于頻點(diǎn)f_k，k∈[B₁,B₂]的頻域數(shù)據(jù)，進(jìn)行后續(xù)處理。

根據(jù)頻域數(shù)組，計(jì)算信號(hào)帶寬內(nèi)頻點(diǎn)的協(xié)方差矩陣R_k＝E[X(f_k)X^H(f_k)](M×M維)，其中X(f_k)＝[x_1k,x_2k,...,x_Mk]^T k∈[B₁,B₂]；

其中，R_k為期望的協(xié)方差矩陣，在實(shí)際應(yīng)用中，R通常由樣本協(xié)方差矩陣代替，在空頻二維濾波器中，對(duì)應(yīng)于頻點(diǎn)f_k的采樣數(shù)據(jù)的頻域協(xié)方差矩陣可表示為如下M×M維的矩陣：

其中L表示頻域數(shù)據(jù)的快拍數(shù)，而x_mk(l)表示第m個(gè)陣元采樣數(shù)據(jù)經(jīng)過(guò)傅立葉變換之后對(duì)應(yīng)頻點(diǎn)f_k上的第l個(gè)快拍數(shù)。

修正協(xié)方差矩陣為α、β為大于0的組合參數(shù)，使修正協(xié)方差矩陣逼近期望的協(xié)方差矩陣，根據(jù)最小方差準(zhǔn)則有，

tr()表示矩陣的跡，求解得

求解協(xié)方差矩陣修正因子α/β，由于實(shí)時(shí)性要求，信號(hào)帶內(nèi)頻域數(shù)據(jù)協(xié)方差積累時(shí)間有限，當(dāng)陣元數(shù)M較大，空間自由度多于衛(wèi)星數(shù)，導(dǎo)向矢量估計(jì)值偏差大，權(quán)值訓(xùn)練的數(shù)據(jù)與權(quán)值應(yīng)用的數(shù)據(jù)之間失配嚴(yán)重，α/β取值越大，反之越小。

在本實(shí)施例中，α/β＝γtr(R_k)/M，γ＝1-e^-0.1M。

多星約束方程如下：

Min E{|Y(f_k)|²}＝w(f_k)^HR(f_k)w(f_k)

s.t.A^H(f_k)w(f_k)＝b

則求解權(quán)向量表達(dá)式為，

根據(jù)協(xié)方差矩陣修正因子計(jì)算權(quán)向量其中I是單位陣，b為N×1維的全1列向量。

在本步驟中，通過(guò)采用協(xié)方差矩陣修正因子，對(duì)樣本協(xié)方差矩陣進(jìn)行修正，可以保證多顆衛(wèi)星信號(hào)數(shù)據(jù)失配時(shí)仍能有效地抑制寬帶干擾，增強(qiáng)了空域?yàn)V波算法的穩(wěn)健性。

步驟103，使用所述權(quán)向量對(duì)所述頻域數(shù)據(jù)進(jìn)行濾波處理；

利用權(quán)值w(f_k)對(duì)頻域數(shù)組進(jìn)行濾波處理，設(shè)Y(f_k)為第k個(gè)頻點(diǎn)空域?yàn)V波的輸出，則

在本步驟中，由于對(duì)于頻點(diǎn)f_k，k∈[1,B₁)or(B₂,K]的頻域數(shù)據(jù)，有w(f_k)＝O_M×1，O_M×1為M×1維的零矢量，所以可以將期望信號(hào)帶外部分進(jìn)行清零處理，降低了運(yùn)算量和運(yùn)算復(fù)雜度。

步驟104，對(duì)濾波后的頻域數(shù)據(jù)進(jìn)行傅里葉逆變換，轉(zhuǎn)換為時(shí)域數(shù)據(jù)并輸出。

其中，可以對(duì)空域?yàn)V波后的頻域數(shù)據(jù)Y(f_k)進(jìn)行K點(diǎn)IFFT(或IDFT)變換，得到抗干擾處理后的時(shí)域中頻數(shù)據(jù)。

如圖2所示，所述時(shí)域中頻數(shù)據(jù)發(fā)送至導(dǎo)航處理單元進(jìn)行下一步處理。

如圖3所示，為本發(fā)明實(shí)施例的多星約束穩(wěn)健空頻抗干擾處理裝置，應(yīng)用于空頻自適應(yīng)處理器，包括：

第一信號(hào)轉(zhuǎn)換模塊21，用于將接收到的時(shí)域信號(hào)通過(guò)傅立葉變換，轉(zhuǎn)換為頻域數(shù)據(jù)；

所述傅立葉變換可以是快速傅立葉變換(FFT)或者離散傅立葉變換(DFT)。

計(jì)算模塊22，用于根據(jù)所述頻域數(shù)據(jù)計(jì)算信號(hào)帶寬內(nèi)頻點(diǎn)的協(xié)方差矩陣，確定協(xié)方差矩陣修正因子，根據(jù)所述協(xié)方差矩陣修正因子計(jì)算權(quán)向量；

可選地，所述計(jì)算模塊22，還用于利用先驗(yàn)信息估計(jì)衛(wèi)星波達(dá)方向，確定空頻導(dǎo)向矢量。

其中，利用先驗(yàn)信息估計(jì)當(dāng)前N顆可見(jiàn)衛(wèi)星的波達(dá)方向θ_i,i＝1,2,...N，求解出對(duì)應(yīng)空頻導(dǎo)向矢量A(f_k)＝[a(θ₁,f_k)a(θ₂,f_k)...a(θ_N,f_k)]，以均勻線陣為例，波達(dá)方向?yàn)棣?sub>1的衛(wèi)星信號(hào)空頻導(dǎo)向矢量為

其中d為陣元間距，M為陣元數(shù)，c為光速常數(shù)。

可采用M陣元的空頻二維處理器結(jié)構(gòu)，對(duì)各陣元中頻AD采樣信號(hào)進(jìn)行K點(diǎn)傅立葉變換變換得到頻域數(shù)組x₁₁,x₁₂,…x_1K；x₂₁,x₂₂,…x_2K；…x_M1,x_M2,…x_MK，數(shù)據(jù)中包含衛(wèi)星信號(hào)、干擾以及噪聲。

采樣數(shù)據(jù)經(jīng)過(guò)K點(diǎn)的傅立葉變換后，有1＜B₁＜B₂＜K，B₁、B₂、K均為正整數(shù)，其中區(qū)間[B₁,B₂]為期望信號(hào)帶寬。為了對(duì)數(shù)據(jù)進(jìn)行有效的處理，同時(shí)降低計(jì)算復(fù)雜度，可以對(duì)頻域數(shù)據(jù)采用分段處理的方式，對(duì)于期望信號(hào)帶外部分進(jìn)行清零處理，減少不必要的計(jì)算工作。

即：對(duì)于頻點(diǎn)f_k，k∈[1,B₁)or(B₂,K]的頻域數(shù)據(jù)，有w(f_k)＝O_M×1，O_M×1為M×1維的零矢量。

對(duì)于頻點(diǎn)f_k，k∈[B₁,B₂]的頻域數(shù)據(jù)，進(jìn)行后續(xù)處理。

根據(jù)頻域數(shù)組，計(jì)算信號(hào)帶寬內(nèi)頻點(diǎn)的協(xié)方差矩陣R_k＝E[X(f_k)X^H(f_k)](M×M維)，其中X(f_k)＝[x_1k,x_2k,...,x_Mk]^T k∈[B₁,B₂]；

其中，R_k為期望的協(xié)方差矩陣，在實(shí)際應(yīng)用中，R通常由樣本協(xié)方差矩陣代替，在空頻二維濾波器中，對(duì)應(yīng)于頻點(diǎn)f_k的采樣數(shù)據(jù)的頻域協(xié)方差矩陣可表示為如下M×M維的矩陣：

其中L表示頻域數(shù)據(jù)的快拍數(shù)，而x_mk(l)表示第m個(gè)陣元采樣數(shù)據(jù)經(jīng)過(guò)FFT之后對(duì)應(yīng)頻點(diǎn)f_k上的第l個(gè)快拍數(shù)。

修正協(xié)方差矩陣為α、β為大于0的組合參數(shù)，使修正協(xié)方差矩陣逼近期望的協(xié)方差矩陣，根據(jù)最小方差準(zhǔn)則有，

tr()表示矩陣的跡，求解得

求解協(xié)方差矩陣修正因子α/β，由于實(shí)時(shí)性要求，信號(hào)帶內(nèi)頻域數(shù)據(jù)協(xié)方差積累時(shí)間有限，當(dāng)陣元數(shù)M較大，空間自由度多于衛(wèi)星數(shù)，導(dǎo)向矢量估計(jì)值偏差大，權(quán)值訓(xùn)練的數(shù)據(jù)與權(quán)值應(yīng)用的數(shù)據(jù)之間失配嚴(yán)重，α/β取值越大，反之越小。

在本實(shí)施例中，α/β＝γtr(R_k)/M，γ＝1-e^-0.1M。

多星約束方程如下：

Min E{|Y(f_k)|²}＝w(f_k)^HR(f_k)w(f_k)

s.t.A^H(f_k)w(f_k)＝b

則求解權(quán)向量表達(dá)式為，

根據(jù)協(xié)方差矩陣修正因子計(jì)算權(quán)向量其中I是單位陣，b為N×1維的全1列向量。

所述計(jì)算模塊22，通過(guò)采用協(xié)方差矩陣修正因子，對(duì)樣本協(xié)方差矩陣進(jìn)行修正，可以保證多顆衛(wèi)星信號(hào)數(shù)據(jù)失配時(shí)仍能有效地抑制寬帶干擾，增強(qiáng)了空域?yàn)V波算法的穩(wěn)健性。

濾波模塊23，用于使用所述權(quán)向量對(duì)所述頻域數(shù)據(jù)進(jìn)行濾波處理；

利用權(quán)值w(f_k)對(duì)頻域數(shù)組進(jìn)行濾波處理，設(shè)Y(f_k)為第k個(gè)頻點(diǎn)空域?yàn)V波的輸出，則

由于對(duì)于頻點(diǎn)f_k，k∈[1,B₁)or(B₂,K]的頻域數(shù)據(jù)，有w(f_k)＝O_M×1，O_M×1為M×1維的零矢量，所以可以將期望信號(hào)帶外部分進(jìn)行清零處理，降低了運(yùn)算量和運(yùn)算復(fù)雜度。

第二信號(hào)轉(zhuǎn)換模塊24，用于對(duì)濾波后的頻域數(shù)據(jù)進(jìn)行傅里葉逆變換，轉(zhuǎn)換為時(shí)域數(shù)據(jù)并輸出。

其中，可以對(duì)空域?yàn)V波后的頻域數(shù)據(jù)Y(f_k)進(jìn)行K點(diǎn)IFFT(或IDFT)變換，得到抗干擾處理后的時(shí)域中頻數(shù)據(jù)。

如圖4所示，為本發(fā)明應(yīng)用示例的空頻自適應(yīng)處理器內(nèi)部處理過(guò)程示意圖，該空頻自適應(yīng)處理器也可稱為頻域?qū)拵Рㄊ纬善鳌?/p>

本應(yīng)用示例中傅里葉變換采用離散傅里葉變換(DFT)。

寬帶信號(hào)可以通過(guò)傅里葉變換分為若干個(gè)頻點(diǎn)，對(duì)于每個(gè)頻點(diǎn)按窄帶波束形成的方法進(jìn)行處理。由K點(diǎn)離散傅里葉變換把寬帶信號(hào)在頻域上劃分為K個(gè)窄子帶，并對(duì)每一個(gè)窄子帶進(jìn)行窄帶波束形成，然后把波束輸出轉(zhuǎn)化為時(shí)域輸出，所以頻域處理是一種劃分子帶的處理方法。如果將一個(gè)觀察時(shí)間子段的觀察數(shù)據(jù)進(jìn)行K點(diǎn)的離散傅立葉變換，可以得到如下的寬帶模型：

X(f_k)＝A(f_k)S(f_k)+N(f_k)k＝1,2,...,K

A(f_k)＝[a(θ₁,f_k)a(θ₂,f_k)...a(θ_N,f_k)]

式中X(f_k)，S(f_k)，N(f_k)分別為對(duì)應(yīng)某頻率的接收數(shù)據(jù)矢量、信號(hào)矢量及噪聲矢量的離散傅立葉變換，K是指將帶寬為B的信號(hào)劃分為K個(gè)子帶。其中，

X(f_k)＝[x_1k,x_2k,...,x_Mk]^T k＝1,2,...,K

用M×1維向量w(f_k)表示對(duì)應(yīng)于接收數(shù)據(jù)f_k子帶的處理器權(quán)矢量，則

w(f_k)＝[w_1k,w_2k,...,w_Mk]^T k＝1,2,...,K

使用前述所述方法計(jì)算w(f_k)，即：

其中，f_k表示第k個(gè)頻點(diǎn)，為樣本協(xié)方差矩陣，α/β為協(xié)方差矩陣修正因子α/β＝γtr(R_k)/M，γ＝1-e^-0.1M，I為單位陣，A為空頻導(dǎo)向矢量，b為N×1維的全1列向量。

則空頻二維處理器的濾波輸出信號(hào)為

綜上所述，傳統(tǒng)空域自適應(yīng)濾波方法很難避免期望信號(hào)自相消現(xiàn)象，惡化系統(tǒng)輸出信干噪比，嚴(yán)重情況下，常規(guī)方法可能完全失效。而本發(fā)明實(shí)施例利用空頻自適應(yīng)處理結(jié)構(gòu)，保證多星空頻導(dǎo)向矢量約束，對(duì)帶外信號(hào)進(jìn)行清零處理，對(duì)帶內(nèi)信號(hào)提出了協(xié)方差矩陣修正因子的具體表達(dá)式。本發(fā)明實(shí)施例既簡(jiǎn)化工程運(yùn)算量，又保證多顆衛(wèi)星信號(hào)數(shù)據(jù)失配時(shí)仍能有效地抑制寬帶干擾，增強(qiáng)了空域?yàn)V波算法的穩(wěn)健性。

雖然本發(fā)明所揭露的實(shí)施方式如上，但所述的內(nèi)容僅為便于理解本發(fā)明而采用的實(shí)施方式，并非用以限定本發(fā)明。任何本發(fā)明所屬領(lǐng)域內(nèi)的技術(shù)人員，在不脫離本發(fā)明所揭露的精神和范圍的前提下，可以在實(shí)施的形式及細(xì)節(jié)上進(jìn)行任何的修改與變化，但本發(fā)明的專利保護(hù)范圍，仍須以所附的權(quán)利要求書所界定的范圍為準(zhǔn)。