CN105158741A - 基于矩阵重构的自适应抗干扰多径多波束形成方法 - Google Patents
基于矩阵重构的自适应抗干扰多径多波束形成方法 Download PDFInfo
- Publication number
- CN105158741A CN105158741A CN201510339083.1A CN201510339083A CN105158741A CN 105158741 A CN105158741 A CN 105158741A CN 201510339083 A CN201510339083 A CN 201510339083A CN 105158741 A CN105158741 A CN 105158741A
- Authority
- CN
- China
- Prior art keywords
- signal
- interference
- coherent
- signals
- matrix
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Noise Elimination (AREA)
- Radio Transmission System (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于矩阵重构的自适应抗干扰多径多波束形成方法,属于数字信号处理领域。本发明首先通过矩阵重构的方法得到干扰信号的自相关矩阵,进而获得干扰子空间的正交补空间,然后将阵列接收数据在干扰子空间的正交补空间上进行投影,达到抑制不相干干扰信号和获得只含相干信号的数据的目的,最后通过特征值分解的方法得到该自适应抗干扰多径多波束形成器的阵列加权矢量。本发明能够在已知干扰信号的个数和入射信号中各个相干信号的入射角度的范围的情况下,在相干信号的入射方向形成多波束把相干信号有效的进行相参积累,充分利用相干信号的信息;并在干扰信号的入射方向形成零陷,达到抑制干扰和提高阵列输出信干噪比的目的。
Description
技术领域
本发明属于数字信号处理领域,涉及阵列信号处理领域中的波束形成技术,具体是指当期望信号存在多径入射,干扰信号与期望信号相互独立时,该波束形成器能在各个相干信号的入射方向形成多个波束把相干信号有效的进行相参积累,并且自适应的在干扰信号的入射方向形成零陷,以达到抑制干扰和提高阵列输出信干噪比的目的。
背景技术
将多个传感器接收到的信号,在空间上增强期望信号,抑制干扰和噪声的处理过程就是波束形成,该技术已成为阵列信号处理中最重要的内容之一。波束形成分为模拟波束形成和数字波束形成(DBF)两部分。数字波束形成(DBF)因为具有自适应波束置零、超分辨定向、超低副瓣、波束切换快速灵活、空时自适应处理等优点,得到了国内外众多研究者的重视。
自适应数字波束形成技术是根据一定的最佳准则,对空间每个天线接收数据进行加权处理后合并,自适应地将天线阵主瓣指向期望方向,而在干扰到达方向形成零点。相对普通相控阵波束形成技术,它有阵列波束图控制更加灵活方便、干扰抑制能力更好、输出信干噪比更高、空间分辨能力更强等优点。尤其是计算机技术和电子技术的快速发展,为自适应波束形成技术的实现和应用提供了条件,使得自适应数字波束形成技术在近年来得到了巨大的发展。然而以往的自适应波束形成主要都是在假设期望信号和干扰信号相互独立的条件下进行的。但是在阵列信号处理中,往往存在有相干干扰信号,如多径反射、智能干扰等。特别是在当代通信系统中,由于多径反射而形成的相干干扰更是普遍存在,在这种情况下,常规自适应波束形成器会引起期望信号对消,使得波束形成性能急剧下降,因此存在相干干扰时的自适应波束形成技术引起了人们越来越多的注意。为了避免期望信号对消,有人提出了不仅能对不相关干扰信号抑制还能对相干干扰信号抑制的自适应波束形成方法。但是,在有些情况下,例如期望信号较弱,甚至是没有期望信号的直达波,只有期望信号的多径入射信号时,为了能充分利用相干信号的信息,波束形成器需要把期望信号和相干信号有效地组合起来,增大阵列输出信干噪比,而不是把相干干扰信号抑制掉。
如图1所示,一M阵元的均匀线阵,阵元间距为λ/2,假设有K个互不相关的干扰信号和一组包含P+1个信号源的相干信号入射,即相干信号包括一个期望信号和P个相干信号。设K个互不相关的干扰信号的入射方向为θuk,k=1,2,…,K,P+1个相干信号的入射方向为θdi,i=0,1,…,P,其中θd0和θdi,i=1,2,…,P分别表示期望信号的入射方向和第i个相干信号的入射方向,且所有信号均为远场窄带信号。那么阵列的接收信号可表示为:
其中为相干信号的合成导向矢量,sd(t)和sk(t)分别为期望信号和第k个不相关干扰信号的复包络;ρi(i=1,2,…,P)表示第i个相干信号与期望信号复包络的相对值,为阵列对于方向θ的导向矢量,n(t)为阵元接收噪声。假设sd(t),{sk(t),k=1,2,…,K}和n(t)互不相关,那么阵列接收信号相关矩阵为:
其中为期望信号功率,E[·]表示求数学期望,A=[au1,au2,…,auK], 为第k个不相关信号的导向矢量;其中s=[s1(t),s2(t),…,sK(t)]T,T表示矩阵转置,为噪声功率,I为单位阵,表示干扰和噪声的自相关矩阵。按如下线性约束最小方差准则来求权矢量:
s.twHbd=1
由拉格朗日法解得权向量为:
w1=μR-1bd(4)
μ为一常数,利用矩阵求逆引理可得:
w1=βRu -1bd(5)
其中β也为一常数。可以证明在最大输出信干噪比(SINR)准则下,权向量w1为最优解。以w1为权矢量的波束形成器能最优地把期望信号和相干干扰信号组合起来,并通过在干扰信号方向形成零陷有效地抑制不相关干扰信号,使阵列输出信干噪比达到最大。因此,该波束形成器被称为最优波束形成器。然而要得到最优波束形成器w1需要已知Ru和bd,要获得bd不仅要知道各相干信号的准确来波方向还要知道各个来波方向的相干信号相对于期望信号复包络的相对值ρi(i=1,2,…,P)。而对于Ru,实际一般是只知道接收数据的协方差矩阵R,要从R中分离出Ru是非常困难的。
当入射信号中没有不相关干扰即K=0时,接收信号的自相关矩阵变为
对R进行特征值分解可得:
是相应的M个特征值,um(m=1,2,…,M)是对应的特征向量。由信号理论的知识可知,可以用R的最大的特征值对应的特征向量u1作为bd的估计值
发明内容
在实际中,由于不相干干扰的存在,会使得相关矩阵有不止一个大的特征值,因而不能用相关矩阵R的最大特征值作为合成导向矢量bd的估计值,另外由于快拍数是有限的,因此式(4)和式(5)实际上是不等价的,所以要想使用式(5)来计算w1获得最优波束形成器是很困难的。在这里本发明提出了一种通过矩阵重构的方法来进行自适应抗干扰多径多波束形成。在实际中,尽管已知的相关信号来波方向θdi,i=0,1,…,P不准确,但是我们却可以准确地知道相关信号来波方向是存在于一个空域区间的,这个空域区间可以表示为Θdi=[θdi-Δθdi,θdi+Δθdi]表示第i个信号的入射角度的范围,这里Δθdi是角度的误差量。首先通过接收信号的相关矩阵R得到入射信号的谱函数P(θ),θ表示空间角度;然后在Θu上通过重构出只含不相关干扰信号的相关矩阵RJ,其中Θu=Θ-Θd,Θ表示整个扫描角度范围。通过不相关干扰信号的相关矩阵RJ得到干扰信号的正交补空间UJ ⊥,再通过将接收数据x(t)投影到正交补空间UJ ⊥达到抑制干扰和得到不含不相关干扰信号的数据,最后通过特征分解求得该波束形成器的阵列加权矢量。
本发明解决上述技术问题所采用的技术方案是:
一种基于矩阵重构的自适应抗干扰多径多波束形成方法,包括以下步骤:
a,通过阵列接收的N次快拍数据x(n)(这里的x(n)与背景技术中的x(t)都是表示阵列的接收信号),估计信号的空间相关矩阵并由得到信号的谱函数P(θ)。
b,根据来重构不相关信号的自相关矩阵,其中Θu=Θ-Θd,Θ表示整个扫描角度范围,Θdi表示第i相干信号的入射角度的范围。
c,对由b得到的RJ进行特征值分解得到K个大的特征值对应的特征向量uk,k=1,2,…,K,由UJ=[u1,u2,…,uK]得到不相关干扰信号的信号子空间UJ。
d,根据UJ ⊥=I-UJ(UJ HUJ)-1UJ H得到干扰子空间的正交补空间,由P=UJ ⊥(UJ ⊥)H得到投影矩阵。
e,将阵列接收数据通过y(n)=Px(n)投影后,得到只含相关信号的数据。
f,根据估计相关信号的空间相关矩阵
g,对由f得到的相关矩阵进行特征值分解,得到最大的特征值对应的特征向量u1就是该波束形成器的阵列加权矢量。
本发明的有益效果是:提出了一种基于矩阵重构的自适应抗干扰多径多波束形成方法,利用本发明获得的阵列加权矢量形成的波束形成器在期望信号存在相干信号的环境里,可以自适应的在各相干信号的入射方向形成多波束,将各相干信号有效的结合起来,充分利用相干信号的信息;并且自适应的在不相干干扰信号的入射方向形成零陷,抑制不相干干扰,提高阵列输出信干噪比;该波束形成器由于不需要已知各信号的精确入射角度信息,只需知道各相干信号的入射角度范围,因而相较于其他需要精确已知信号入射角度的波束形成器有更好的稳健性。另外由于该波束形成器是通过将阵列接收数据投影到干扰子空间的正交补空间上来抑制干扰的,因而它不仅适用于强期望环境,在弱期望环境也能很好的抑制干扰。由于不相关干扰信号的自相关矩阵是通过矩阵重构获得的,因而该波束形成器有较快的收敛速度。
附图说明
图1为阵元数为M的阵列的结构示意图。
具体实施方式
下面结合实施例对本发明的方案作进一步的描述。
实施例:
该方法的前提是已经获得了不相干干扰信号的个数K和相干信号的入射角度的范围Θdi表示第i个相干信号的入射角度的范围。
本例中的基于矩阵重构的自适应抗干扰多径多波束形成方法包括以下具体实现步骤:
a,通过阵列接收的N次快拍数据,根据估计信号的空间相关矩阵
快拍次数N的大小影响着信号的空间相关矩阵估计的精确性,N越大,空间相关矩阵的估计就越精确。因而在实际的应用中,在允许的范围内,应尽可能的增加快拍数,相关矩阵的获得不仅限于上述方法,也可通过其他方法来获得接收信号的相关矩阵。
b,根据来重构不相关干扰信号的自相关矩阵,其中Θu=Θ-Θd,Θ表示整个扫描角度范围,Θdi表示第i个相干信号的入射角度的范围。
在这里我们通过MVDR谱来重构不相关信号的自相关矩阵,但是不相关信号的自相关矩阵的重构不应只限于通过MVDR谱,还可以通过其他的谱函数如MUSIC谱等来重构不相关信号的自相关矩阵。
c,对由b得到的RJ进行特征值分解得到K个大的特征值对应的特征向量uk,k=1,2,…,K,由UJ=[u1,u2,…,uK]得到不相关干扰信号的信号子空间UJ。
d,根据UJ ⊥=I-UJ(UJ HUJ)-1UJ H得到干扰子空间的正交补空间,由P=UJ ⊥(UJ ⊥)H得到投影矩阵。
e,将阵列接收数据通过y(n)=Px(n)进行投影,得到只含相关信号的数据。
f,根据估计相关信号的空间相关矩阵
g,对由f得到的相关矩阵进行特征值分解,得到最大的特征值对应的特征向量u1就是该多径多波束形成器的阵列加权矢量。
Claims (1)
1.一种基于矩阵重构的自适应抗干扰多径多波束形成方法,其特征在于,包括以下步骤:
a,通过阵列接收的N次快拍数据x(n),估计信号的空间相关矩阵并由得到信号的谱函数P(θ);
b,根据来重构不相关信号的自相关矩阵,其中Θu=Θ-Θd,Θ表示整个扫描角度范围,Θdi表示第i相干信号的入射角度的范围;
c,对由b得到的RJ进行特征值分解得到K个大的特征值对应的特征向量uk,k=1,2,…,K,由UJ=[u1,u2,…,uK]得到不相关干扰信号的信号子空间UJ;
d,根据UJ ⊥=I-UJ(UJ HUJ)-1UJ H得到干扰子空间的正交补空间,由P=UJ ⊥(UJ ⊥)H得到投影矩阵;
e,将阵列接收数据通过y(n)=Px(n)投影后,得到只含相关信号的数据;
f,根据估计相关信号的空间相关矩阵
g,对由f得到的相关矩阵进行特征值分解,得到最大的特征值对应的特征向量u1就是波束形成器的阵列加权矢量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510339083.1A CN105158741B (zh) | 2015-06-17 | 2015-06-17 | 基于矩阵重构的自适应抗干扰多径多波束形成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510339083.1A CN105158741B (zh) | 2015-06-17 | 2015-06-17 | 基于矩阵重构的自适应抗干扰多径多波束形成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105158741A true CN105158741A (zh) | 2015-12-16 |
CN105158741B CN105158741B (zh) | 2018-07-27 |
Family
ID=54799659
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510339083.1A Expired - Fee Related CN105158741B (zh) | 2015-06-17 | 2015-06-17 | 基于矩阵重构的自适应抗干扰多径多波束形成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105158741B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107894597A (zh) * | 2017-10-27 | 2018-04-10 | 西安电子科技大学 | 基于线性虚拟阵元的自适应调零方法 |
CN107918112A (zh) * | 2017-11-10 | 2018-04-17 | 西安电子科技大学 | 基于有源‑无源雷达协同的联合收发波束形成方法 |
WO2018119891A1 (zh) * | 2016-12-29 | 2018-07-05 | 深圳天珑无线科技有限公司 | 适用于多天线系统的阻塞信号配置方法及装置 |
CN108572347A (zh) * | 2017-03-09 | 2018-09-25 | 上海交通大学 | 基于通信信号信道状态响应信息的面阵二维测角方法及系统 |
CN109633692A (zh) * | 2018-11-26 | 2019-04-16 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | Gnss导航卫星信号抗干扰处理方法 |
CN109633563A (zh) * | 2019-01-22 | 2019-04-16 | 西安电子科技大学 | 基于多径信息的自适应相干波束形成方法 |
CN109901131A (zh) * | 2019-03-24 | 2019-06-18 | 西安电子科技大学 | 基于斜投影的多径利用相干波束形成方法 |
CN109901132A (zh) * | 2019-03-24 | 2019-06-18 | 西安电子科技大学 | 一种多径利用的自适应相干波束形成方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5260707A (en) * | 1988-12-22 | 1993-11-09 | Hughes Aircraft Company | Phase coherent interference signal suppression system and method |
GB2318914B (en) * | 1996-10-30 | 2001-05-09 | Motorola Inc | Method and intelligent digital beam forming system with improved signal quality communications |
EP1225455A2 (en) * | 1994-08-10 | 2002-07-24 | Mcdonnell Douglas Corporation | Radar apparatus with RFI suppression |
JP2003014836A (ja) * | 2001-07-02 | 2003-01-15 | Mitsubishi Electric Corp | レーダ装置 |
CN103984676A (zh) * | 2014-04-24 | 2014-08-13 | 北京理工大学 | 一种基于协方差矩阵重构的正交投影自适应波束形成方法 |
-
2015
- 2015-06-17 CN CN201510339083.1A patent/CN105158741B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5260707A (en) * | 1988-12-22 | 1993-11-09 | Hughes Aircraft Company | Phase coherent interference signal suppression system and method |
EP1225455A2 (en) * | 1994-08-10 | 2002-07-24 | Mcdonnell Douglas Corporation | Radar apparatus with RFI suppression |
GB2318914B (en) * | 1996-10-30 | 2001-05-09 | Motorola Inc | Method and intelligent digital beam forming system with improved signal quality communications |
JP2003014836A (ja) * | 2001-07-02 | 2003-01-15 | Mitsubishi Electric Corp | レーダ装置 |
CN103984676A (zh) * | 2014-04-24 | 2014-08-13 | 北京理工大学 | 一种基于协方差矩阵重构的正交投影自适应波束形成方法 |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018119891A1 (zh) * | 2016-12-29 | 2018-07-05 | 深圳天珑无线科技有限公司 | 适用于多天线系统的阻塞信号配置方法及装置 |
CN108572347A (zh) * | 2017-03-09 | 2018-09-25 | 上海交通大学 | 基于通信信号信道状态响应信息的面阵二维测角方法及系统 |
CN107894597B (zh) * | 2017-10-27 | 2021-01-01 | 西安电子科技大学 | 基于线性虚拟阵元的自适应调零方法 |
CN107894597A (zh) * | 2017-10-27 | 2018-04-10 | 西安电子科技大学 | 基于线性虚拟阵元的自适应调零方法 |
CN107918112A (zh) * | 2017-11-10 | 2018-04-17 | 西安电子科技大学 | 基于有源‑无源雷达协同的联合收发波束形成方法 |
CN107918112B (zh) * | 2017-11-10 | 2019-11-15 | 西安电子科技大学 | 基于有源-无源雷达协同的联合收发波束形成方法 |
CN109633692A (zh) * | 2018-11-26 | 2019-04-16 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | Gnss导航卫星信号抗干扰处理方法 |
CN109633692B (zh) * | 2018-11-26 | 2022-07-08 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | Gnss导航卫星信号抗干扰处理方法 |
CN109633563A (zh) * | 2019-01-22 | 2019-04-16 | 西安电子科技大学 | 基于多径信息的自适应相干波束形成方法 |
CN109633563B (zh) * | 2019-01-22 | 2023-03-10 | 西安电子科技大学 | 基于多径信息的自适应相干波束形成方法 |
CN109901131A (zh) * | 2019-03-24 | 2019-06-18 | 西安电子科技大学 | 基于斜投影的多径利用相干波束形成方法 |
CN109901132A (zh) * | 2019-03-24 | 2019-06-18 | 西安电子科技大学 | 一种多径利用的自适应相干波束形成方法 |
CN109901131B (zh) * | 2019-03-24 | 2022-12-02 | 西安电子科技大学 | 基于斜投影的多径利用相干波束形成方法 |
CN109901132B (zh) * | 2019-03-24 | 2022-12-02 | 西安电子科技大学 | 一种多径利用的自适应相干波束形成方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105158741B (zh) | 2018-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105158741A (zh) | 基于矩阵重构的自适应抗干扰多径多波束形成方法 | |
CN107315162B (zh) | 基于内插变换和波束形成的远场相干信号doa估计方法 | |
CN101369014B (zh) | 应用于多输入多输出雷达的双边约束自适应波束形成方法 | |
CN109407055B (zh) | 基于多径利用的波束形成方法 | |
CN109031231B (zh) | 雷达低空目标时间反演相干角度估计方法 | |
CN109254261B (zh) | 基于均匀圆阵epuma的相干信号零陷加深方法 | |
CN104730491A (zh) | 一种基于l型阵的虚拟阵列doa估计方法 | |
Chen et al. | Adaptive beamforming with sensor position errors using covariance matrix construction based on subspace bases transition | |
CN105676168A (zh) | 一种声矢量阵方位估计方法 | |
CN103837861A (zh) | 基于特征子空间的子阵级线性约束自适应波束形成方法 | |
Bingbing | DOA estimation of the coherent signals using beamspace matrix reconstruction | |
CN104156553B (zh) | 无需信源数估计的相干信号波达方向估计方法及系统 | |
Luo et al. | Mainlobe anti-jamming via eigen-projection processing and covariance matrix reconstruction | |
CN107390188B (zh) | 一种用于mimo雷达主瓣抗干扰的发射初相优化方法 | |
CN110196417A (zh) | 基于发射能量集中的双基地mimo雷达角度估计方法 | |
Yang et al. | Estimation of the DOAs of coherent signals in beam space processing for phased arrays | |
Wang et al. | Eigenspace-based beamforming technique for multipath coherent signals reception | |
Du et al. | Coherent signals adaptive beamforming algorithm based on eigenvalue | |
CN115575941A (zh) | 面向稀疏阵列的频控阵mimo雷达目标参数估计方法 | |
CN114184999B (zh) | 一种互耦小孔径阵列的生成式模型处理方法 | |
CN109633563B (zh) | 基于多径信息的自适应相干波束形成方法 | |
Han et al. | Direction finding for both uncorrelated and coherent signals in improved uniform circular array | |
Tang et al. | New robust adaptive beamforming method for multipath coherent signal reception | |
Cheng et al. | A reduced-complex method based on Toeplitz reconstruction for direction of arrival estimation in multiple-input multiple-output sonar | |
Tian et al. | Direction estimation of coherent signals based on the symmetry of uniform linear array |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180727 Termination date: 20210617 |