CN108631851B - 基于均匀直线阵列零陷加深的自适应波束形成方法 - Google Patents
基于均匀直线阵列零陷加深的自适应波束形成方法 Download PDFInfo
- Publication number
- CN108631851B CN108631851B CN201711024309.4A CN201711024309A CN108631851B CN 108631851 B CN108631851 B CN 108631851B CN 201711024309 A CN201711024309 A CN 201711024309A CN 108631851 B CN108631851 B CN 108631851B
- Authority
- CN
- China
- Prior art keywords
- matrix
- interference
- signal
- array
- vector
- 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.)
- Active
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
- H04B7/0837—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
- H04B7/0842—Weighted combining
- H04B7/086—Weighted combining using weights depending on external parameters, e.g. direction of arrival [DOA], predetermined weights or beamforming
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
- H04B7/0456—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
- H04B7/0837—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
- H04B7/0842—Weighted combining
- H04B7/0848—Joint weighting
- H04B7/0854—Joint weighting using error minimizing algorithms, e.g. minimum mean squared error [MMSE], "cross-correlation" or matrix inversion
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明提出了一种基于均匀直线阵列零陷加深的自适应波束形成方法,主要解决现有技术中由于干扰功率起伏,天线方向图零陷的深度不能完全抑制干扰的问题。其实现步骤为:1.获取均匀天线阵列的相关参数,来波信号的相关参数;2.利用天线阵列接收数据,计算天线阵列接收数据的自相关矩阵;3.利用天线阵列的相关参数计算期望信号的方向矢量及导向矢量矩阵;4.利用矩阵奇异值分解,对自相关矩阵进行分解,构造新的干扰子空间;5.利用新的干扰子空间矩阵和期望信号的方向矢量,计算零陷加深后的自适应波束形成最优权系数。本发明加深了干扰信号的零陷,提升了抗干扰性能,可用于信号处理系统。
Description
技术领域
本发明属于信号处理技术领域,特别涉及一种自适应波束形成方法,可用于信号抗干扰处理系统。
背景技术
在通信和雷达领域,自适应天线是及其重要的部件之一。自适应天线能够自动将天线最大辐射方向对准所需信号的来波方向,将天线方向图的波瓣零陷对准干扰方向,具有极为灵活可靠的抗干扰检测能力,是空域上反电子对抗的一种有效措施,其由天线阵列和自适应处理系统构成,故又称自适应阵列。
自适应处理系统的核心是自适应波束形成,自适应波束形成是指根据一定的波束形成准则和算法实时自适应地调整权向量,使得各阵元的接收信号通过加权叠加后,输出信号的质量在所采用的准则下最优使期望信号得到最大的输出功率,零陷对准干扰方向,对干扰信号进行抑制,从而提高阵列输出的信干噪比。
自适应波束形成中的最佳滤波准则有最小均方误差MMSE准则、最小方差MV准则、最大信干噪比SINR准则等。其中:
最小均方误差MMSE准则使阵列输出信号与期望信号的均方误差最小,该准则不需要已知期望信号的波达方向,但会产生干扰信号。
最小方差准则MV在广义约束条件下使阵列输出方差最小,该准则必须已知期望信号的波达方向。
最大信干噪比SINR准则使有用信号功率与干扰信号及噪声功率比最大,该准则可以获得最大的信噪比,但必须已知噪声统计特性和期望信号的波达方向。
上述准则在自适应波束形成中应用广泛,但当干扰功率起伏时,自适应波束形成的零陷深度不足以完全抑制干扰,使得波束形成器的输出信号中仍含有干扰分量,抗干扰性能急剧下降。
发明内容
本发明的目的在于针对上述现有技术存在的不足,提出一种基于均匀直线阵列零陷加深的自适应波束形成方法,以提高干扰功率起伏时的抗干扰性能。
本发明的技术思路是:通过采用均匀直线阵列为信号的接收模型,对天线阵列接收信号的自相关矩阵进行奇异值分解,增大干扰子空间,提高干扰功率。通过计算改进后的干扰子空间矩阵、期望信号的方向矢量,得到零陷加深后的自适应波束形成最优权系数,其实现步骤包括如下:
1、一种基于均匀直线阵列零陷加深的自适应波束形成方法,其特征在于,包括:
(1)获取均匀天线阵列的阵元数N和阵元间距d,天线阵列接收数据x(t),获取期望信号的波达方向θ0,干扰信号个数M、干扰信号的波达方向θk,信号源频率f和发射信号波长λ;k=1,2,...,M;
(2)利用天线阵列接收数据x(t),计算天线阵列接收数据x(t)的自相关矩阵Rxx:
Rxx=E[x(t)·(x(t))H],
其中,E[·]表示·的数学期望,H表示共轭转置;
(3)利用步骤(1)获取阵列天线的阵元数N和阵元间距d,计算信号源在第i个阵元的时延τi,根据该时延τi,计算期望信号的方向矢量a(θ0)及均匀直线阵的导向矢量矩阵A,i=1,2,...,N;;
(4)利用矩阵奇异值分解定理,对自相关矩阵Rxx进行分解:
其中,σj表示自相关矩阵Rxx特征值分解后的第j个特征值,uji表示自相关矩阵Rxx特征值分解后的第j个特征向量中第i个矢量,j=1,2,...,N;
(5)利用分解后的自相关矩阵Rxx,构造新的干扰子空间:
(5a)利用分解后的自相关矩阵Rxx,得到干扰子空间的M个大特征值Λ及对应的特征向量U:
Λ=diag(σ1,σ2,...,σk,...σM)
其中,其中diag(σ1,σ2,...,σk,...σM)表示对角线元素σ1,σ2,...,σk,...,σM的对角矩阵,σk表示干扰子空间的第k个特征值,k=1,2,...,M。U1i表示Rxx特征值分解后的第1个特征向量中的第i个矢量,Uki表示Rxx特征值分解后的第k个特征向量中的第i个矢量;
(5b)利用干扰子空间的特征值Λ,特征向量U,计算干扰子空间新的特征值Λ1和特征向量U1:
U1=P·U
其中,P是可以增强干扰子空间但不改变波束到达方向的参数,P=α2·AAH。α1,α2是两个数值不同的调节因子,取值范围为1.5~5,A是均匀直线阵的导向矢量矩阵;
(5c)利用新的特征值Λ1和新的特征向量U1,计算新的干扰子空间矩阵Rx:
Rx=U1Λ1U1 H,
(6)利用新的干扰子空间矩阵Rx和期望信号的方向矢量a(θ0),计算零陷加深后的自适应波束形成最优权系数ω:
其中,θ0表示期望信号的波达方向,Rx -1表示Rx的逆矩阵。
本发明与现有技术相比,具有如下优点:
1.本发明利用天线阵列的接收数据,对天线阵列接收信号的自相关矩阵进行奇异值分解,增大干扰子空间,提高干扰功率,使得天线方向图的零陷加深,解决了现有技术中由于干扰功率起伏、自适应波束形成的零陷深度不能完全抑制干扰的问题,提高了系统的抗干扰性能。
2.本发明直接将自相关矩阵进行奇异值分解,得到其特征值和特征向量,构建新的干扰子空间,弥补了利用阵列分解自相关矩阵时天线损失的孔径,且不会增加天线阵列排布的复杂度。
附图说明
图1为本发明的实现流程图。
图2为本发明与传统自适应调零方法的对比曲线图。
具体实施方式
以下结合附图,对本发明作进一步详细说明。
参照图1,本发明的实现步骤包括如下:
步骤1,获取均匀天线阵列的相关参数和来波信号的相关参数。
设置来波信号的相关参数,包括干扰信号个数M、第k个信号的波达方向θk,k=1,2,...,M,期望信号的波达方向θ0、信号源频率f和发射信号波长λ;
根据天线排布情况,获取线性阵列天线的相关参数,包括线阵的阵元个数N和阵元间隔d,天线阵列接收数据x(t)。
步骤2,利用天线阵列接收数据x(t),计算天线阵列接收数据x(t)的自相关矩阵Rxx:
Rxx=E[x(t)·(x(t))H]
其中,E[·]表示·的数学期望,H表示共轭转置。
步骤3,根据阵列天线的相关参数,计算期望信号的方向矢量a(θ0)及均匀直线阵的导向矢量矩阵A。
3a)利用步骤1获取的阵列天线的阵元间隔d,计算信号源在第i个阵元的时延τi:
其中,c表示光速,θ0表示期望信号的波达方向,i=1,2,...,N;
3b)利用信号源在第i个阵元的时延τi,计算期望信号的方向矢量a(θ0):
其中,T表示转置,τi表示信号源在第i个阵元的时延,f表示信号源频率,j表示虚数单位;
3c)利用信号源在第i个阵元的时延τi,计算均匀直线阵的导向矢量矩阵A:
A=[a(θ1),a(θ2),...,a(θk),...,a(θM)],
步骤4,利用矩阵奇异值分解定理,对自相关矩阵Rxx进行分解如下:
其中,σj表示自相关矩阵Rxx特征值分解后的第j个特征值,uji表示自相关矩阵Rxx特征值分解后的第j个特征向量中第i个矢量,j=1,2,...,N。
步骤5,利用分解后的自相关矩阵Rxx,构造新的干扰子空间。
5a)利用分解后的自相关矩阵Rxx,得到干扰子空间的M个大特征值Λ及对应的特征向量U:
Λ=diag(σ1,σ2,...,σk,...σM)
其中,其中diag(σ1,σ2,...,σk,...σM)表示对角线元素σ1,σ2,...,σk,...,σM的对角矩阵,σk表示干扰子空间的第k个特征值,k=1,2,...,M。U1i表示Rxx特征值分解后的第1个特征向量中的第i个矢量,Uki表示Rxx特征值分解后的第k个特征向量中的第i个矢量;
5b)利用干扰子空间的特征值Λ,特征向量U,计算干扰子空间新的特征值Λ1和特征向量U1:
U1=P·U
其中,P是可以增强干扰子空间但不改变波束到达方向的参数,P=α2·AAH,
α1,α2是两个数值不同的调节因子,取值范围为1.5~5,A是均匀直线阵的导向矢量矩阵;
5c)利用新的特征值Λ1和新的特征向量U1,计算新的干扰子空间矩阵Rx:
Rx=U1Λ1U1 H。
步骤6,利用新的干扰子空间矩阵Rx和期望信号的方向矢量a(θ0),计算零陷加深后的自适应波束形成最优权系数ω:
其中,θ0表示期望信号的波达方向,Rx -1表示Rx的逆矩阵。
本发明的效果可以通过以下仿真进一步说明:
1.仿真条件:
设置天线阵是4天线均匀线阵,天线间距是80cm;
设置一路期望信号,其中心频率是1561.098MHz,带宽是4.092MHz,方向角是-60°;设置三路干扰信号,其中心频率分别是1561.098MHz,1200MHz,1600MHz,带宽分别是4.092MHz,2MHz,6MHz,三路信号方向角分别是-20°,30°,60°;天线接收到的四路信号不相干;
设置干噪比是50dB,环境中信噪比是-20dB。
2.仿真内容
在上述仿真条件下,将4路信号入射到由4根天线组成的均匀线阵,用本发明和传统自适应调零算法对该均匀线阵进行自适应波束形成,结果图2,从图2可见本发明加深了干扰信号的零陷深度,有效抑制干扰信号。
Claims (4)
1.一种基于均匀直线阵列零陷加深的自适应波束形成方法,其特征在于,包括:
(1)获取均匀天线阵列的阵元数N和阵元间距d,天线阵列接收数据x(t),获取期望信号的波达方向θ0,干扰信号个数M,干扰信号的波达方向θk,信号频率f和发射信号波长λ;k=1,2,...,M;
(2)利用天线阵列接收数据x(t),计算天线阵列接收数据x(t)的自相关矩阵Rxx:
Rxx=E[x(t)·(x(t))H],
其中,E[·]表示·的数学期望,H表示共轭转置;
(3)利用步骤(1)获取的阵列天线阵元数N和阵元间距d,计算信号源在第i个阵元的时延τi,根据该时延τi,计算期望信号的方向矢量a(θ0)及均匀直线阵的导向矢量矩阵A,i=1,2,...,N;
(4)利用矩阵奇异值分解定理,对自相关矩阵Rxx进行分解:
其中,σj表示自相关矩阵Rxx特征值分解后的第j个特征值,uji表示自相关矩阵Rxx特征值分解后的第j个特征向量中第i个矢量,j=1,2,...,N;
(5)利用分解后的自相关矩阵Rxx,构造新的干扰子空间:
(5a)利用分解后的自相关矩阵Rxx,得到干扰子空间的特征值Λ及对应的特征向量U:
Λ=diag(σ1,σ2,...,σk,...σM)
其中,其中diag(σ1,σ2,...,σk,...σM)表示对角线元素σ1,σ2,...,σk,...,σM的对角矩阵,σk表示干扰子空间的第k个特征值,k=1,2,...,M;uki表示Rxx特征值分解后的第k个特征向量中的第i个矢量;
(5b)利用干扰子空间的特征值Λ,特征向量U,计算干扰子空间新的特征值Λ1和特征向量U1:
U1=P·U
其中,P是增强干扰子空间但不改变波束到达方向的参数,P=α2·AAH;α1,α2是两个数值不同的调节因子,取值范围为1.5~5,A是均匀直线阵的导向矢量矩阵;
(5c)利用新的特征值Λ1和新的特征向量U1,计算新的干扰子空间矩阵Rx:
Rx=U1Λ1U1 H,
(6)利用新的干扰子空间矩阵Rx和期望信号的方向矢量a(θ0),计算零陷加深后的自适应波束形成最优权系数ω:
其中,θ0表示期望信号的波达方向,Rx -1表示Rx的逆矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711024309.4A CN108631851B (zh) | 2017-10-27 | 2017-10-27 | 基于均匀直线阵列零陷加深的自适应波束形成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711024309.4A CN108631851B (zh) | 2017-10-27 | 2017-10-27 | 基于均匀直线阵列零陷加深的自适应波束形成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108631851A CN108631851A (zh) | 2018-10-09 |
CN108631851B true CN108631851B (zh) | 2021-04-20 |
Family
ID=63705761
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711024309.4A Active CN108631851B (zh) | 2017-10-27 | 2017-10-27 | 基于均匀直线阵列零陷加深的自适应波束形成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108631851B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109683519A (zh) * | 2018-12-21 | 2019-04-26 | 南京航空航天大学 | 一种自适应算法的fpga实现方法 |
CN111863012A (zh) * | 2020-07-31 | 2020-10-30 | 北京小米松果电子有限公司 | 一种音频信号处理方法、装置、终端及存储介质 |
CN112816948A (zh) * | 2020-12-31 | 2021-05-18 | 中国电子科技集团公司第十四研究所 | 一种利用埃尔米特矩阵对称性实现自适应信号处理的方法 |
CN113472371B (zh) * | 2021-05-31 | 2022-09-02 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 自适应阵列天线数字波束合成抗干扰处理方法 |
CN113949424B (zh) * | 2021-10-22 | 2023-12-22 | 杭州红岭通信息科技有限公司 | 一种下行空分干扰消除的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005020445A2 (en) * | 2003-08-05 | 2005-03-03 | University Of Hawai'i | Microwave self-phasing antenna arrays for secure data transmission & satellite network crosslinks |
CN103984676A (zh) * | 2014-04-24 | 2014-08-13 | 北京理工大学 | 一种基于协方差矩阵重构的正交投影自适应波束形成方法 |
CN105022268A (zh) * | 2015-07-09 | 2015-11-04 | 哈尔滨工程大学 | 一种线性约束虚拟天线波束形成方法 |
CN105049382A (zh) * | 2015-06-18 | 2015-11-11 | 哈尔滨工程大学 | 一种抗期望信号导向矢量失配的零陷展宽自适应天线波束成形方法 |
CN107124216A (zh) * | 2017-04-07 | 2017-09-01 | 广东精点数据科技股份有限公司 | 一种针对阵列误差的Capon稳健自适应波束形成方法及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8773307B2 (en) * | 2010-09-09 | 2014-07-08 | Spatial Digital Systems, Inc. | Wide null Forming system with beamforming |
-
2017
- 2017-10-27 CN CN201711024309.4A patent/CN108631851B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005020445A2 (en) * | 2003-08-05 | 2005-03-03 | University Of Hawai'i | Microwave self-phasing antenna arrays for secure data transmission & satellite network crosslinks |
CN103984676A (zh) * | 2014-04-24 | 2014-08-13 | 北京理工大学 | 一种基于协方差矩阵重构的正交投影自适应波束形成方法 |
CN105049382A (zh) * | 2015-06-18 | 2015-11-11 | 哈尔滨工程大学 | 一种抗期望信号导向矢量失配的零陷展宽自适应天线波束成形方法 |
CN105022268A (zh) * | 2015-07-09 | 2015-11-04 | 哈尔滨工程大学 | 一种线性约束虚拟天线波束形成方法 |
CN107124216A (zh) * | 2017-04-07 | 2017-09-01 | 广东精点数据科技股份有限公司 | 一种针对阵列误差的Capon稳健自适应波束形成方法及系统 |
Non-Patent Citations (3)
Title |
---|
A novel adaptive nulling of interference method for array radar;Yang Gao等;《IET International Radar Conference 2015》;20151116;全文 * |
相干信号零陷加深自适应波束形成算法;黄超等;《探测与控制学报》;20160630;第38卷(第3期);全文 * |
零陷加深的自适应波束形成技术研究;朱德智等;《中国电子科学研究院学报》;20120831;第7卷(第4期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108631851A (zh) | 2018-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108631851B (zh) | 基于均匀直线阵列零陷加深的自适应波束形成方法 | |
CN107276658B (zh) | 色噪声下基于协方差矩阵重构的波束形成方法 | |
CN105137399B (zh) | 基于斜投影滤波的雷达自适应波束形成方法 | |
CN103837861B (zh) | 基于特征子空间的子阵级线性约束自适应波束形成方法 | |
CN107979404B (zh) | 基于虚拟阵列零陷加宽的自适应波束形成方法 | |
CN106324625B (zh) | 一种基于二范数多目标优化的卫星导航系统自适应抗干扰方法 | |
CN109254261A (zh) | 基于均匀圆阵epuma的相干信号零陷加深方法 | |
CN107340499A (zh) | 基于协方差矩阵重建的稳健低副瓣波束形成方法 | |
CN109283496B (zh) | 一种抗运动干扰和抗导向失配的稳健波束形成方法 | |
CN113884979A (zh) | 干扰加噪声协方差矩阵重构的稳健自适应波束形成方法 | |
CN113805139A (zh) | 基于聚焦变换的宽带信号稀疏表示波达方向估计方法 | |
CN113406576A (zh) | 一种基于特征斜投影协方差矩阵重构的主瓣干扰抑制方法 | |
Kikuchi et al. | Autocalibration algorithm for robust Capon beamforming | |
CN103248412B (zh) | 一种卫星导航多波束形成的方法 | |
CN110261814B (zh) | 基于空间谱重构和导向矢量直接估计的波束形成方法 | |
CN106685507A (zh) | 色噪声环境下基于约束Kalman波束形成方法 | |
CN113109768B (zh) | 零点约束的稳健自适应波束形成方法 | |
Reza et al. | Robust centered element concentric circular antenna array with low side lobe using variable loading and tapering windows in the presence of array imperfections | |
CN110208830B (zh) | 一种基于空时二维稀疏阵列的导航抗干扰方法 | |
CN107894597B (zh) | 基于线性虚拟阵元的自适应调零方法 | |
CN110018466B (zh) | 一种基于分子阵预处理的music波束形成方法 | |
CN111585631A (zh) | 一种提高分布式极化敏感阵列发射增益的波束形成方法 | |
CN113206696B (zh) | 空域抗干扰赋形方法 | |
MANAI et al. | A Hybrid Adaptive Beamforming Algorithm for SINR Enhancement in Massive MIMO Systems | |
CN111983647B (zh) | 一种获得卫星导航接收机的抗干扰信号的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |