CN105306117A - 一种基于协方差矩阵扩展的半虚拟天线阵波束形成方法 - Google Patents
一种基于协方差矩阵扩展的半虚拟天线阵波束形成方法 Download PDFInfo
- Publication number
- CN105306117A CN105306117A CN201510664784.2A CN201510664784A CN105306117A CN 105306117 A CN105306117 A CN 105306117A CN 201510664784 A CN201510664784 A CN 201510664784A CN 105306117 A CN105306117 A CN 105306117A
- Authority
- CN
- China
- Prior art keywords
- array
- covariance matrix
- theta
- antenna array
- interference
- 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
-
- 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/0426—Power distribution
- H04B7/043—Power distribution using best eigenmode, e.g. beam forming or beam steering
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开了一种基于协方差矩阵扩展的半虚拟天线阵波束形成方法。包括以下步骤,步骤一:将阵列接收数据的协方差矩阵进行扩展,得到协方差扩展矩阵;步骤二:将天线阵列导向矢量进行扩展,得到半虚拟天线阵列的导向矢量;步骤三:计算基于协方差矩阵扩展的Capon波束形成器的最优权;步骤四:利用最优权进行波束赋形,得到阵列天线的输出数据。本发明通过构造虚拟阵元,获得了更高的输出信干噪比,通过约束条件抑制方向已知的干扰,实现抑制超自由度干扰。
Description
技术领域
本发明属于天线阵列的控制领域,尤其涉及一种基于协方差矩阵扩展的半虚拟天线阵波束形成方法。
背景技术
自适应天线波束成形技术能够在期望信号方向形成最大的增益,在干扰方向形成零陷,被广泛应用于雷达,声纳,移动通信等领域。最小方差无畸变相应(MVDR,又名Capon)波束形成器可以提高自适应天线的输出信干噪比,是一种最常见的波束形成器,但其旁瓣较高,而且无法在干扰数超过天线阵自由度的情况下工作。
虚拟天线阵就是针对自适应天线无法在干扰数超过天线阵自由度情况下工作的问题,在自适应天线基础上发展起来的一种天线形式,主要研究增加虚拟阵元,将实际天线阵变换到虚拟天线阵的方法,既用于DOA估计,也用于波束形成,主要包括高阶统计量法(参考文献:基于高阶统计量的自适应非线性波束形成新算法[J].电子学报,1997,25(3):102-104)和内插变换法(Aneffectivetechniqueforenhancinganti-interferenceperformanceofadaptivevirtualantennaarray[J].APPLIEDCOMPUTATIONALELECTROMAGNETICSSOCIETYJOURNAL,2011,26(3):234-240)。同自适应天线波束形成相比,虚拟天线阵波束形成不仅适用于干扰数超过天线阵自由度的情况,还可以获得更高的输出信干噪比,也就是说虚拟天线阵可以比自适应天线抑制更多的干扰和获得更好的天线阵性能。
虚拟天线技术通过构造虚拟阵元,获得了更好的天线阵性能,这是因为虚拟阵元的存在相当于增加了天线阵的阵元数,存在额外的处理增益。
发明内容
本发明的目的是提供一种既能获得比自适应天线更好的性能,又能抑制超过天线阵自由度干扰的,基于协方差矩阵扩展的半虚拟天线阵波束形成方法。
一种基于协方差矩阵扩展的半虚拟天线阵波束形成方法,包括以下步骤,
步骤一:将阵列接收数据的协方差矩阵R进行扩展,得到协方差扩展矩阵
步骤二:将天线阵列导向矢量a(θ)进行扩展,得到半虚拟天线阵列的导向矢量
步骤三:计算基于协方差矩阵扩展的Capon波束形成器的最优权为:
其中,b(θd)为半虚拟天线阵的期望信号导向矢量;
步骤四:利用最优权进行波束赋形,得到阵列天线的输出数据。
有益效果:
针对自适应阵列天线波束形成在干扰数超过天线阵自由度的情况下失效的问题。本发明通过对协方差矩阵的扩展和阵列导向矢量的扩展,形成虚拟阵元,能够构建更多的约束条件,可以通过约束条件抑制超过天线阵自由度的方向已知的干扰,实现在干扰数超过天线阵自由度的情况下工作。
本发明由于构造了虚拟阵元,相当于增加了阵元数,获得了额外的处理增益,输出信干噪比与虚拟天线阵波束形成方法相近,高于自适应天线波束形成方法,获得了比自适应天线更好的性能。
附图说明
图1是本发明的实现步骤示意图;
图2为干扰数小于阵元数时的波束图;
图3为快拍数改变时的输出信干噪比曲线;
图4为输入信干噪比改变时的输出信干噪比曲线。
图5为方向未知的干扰加方向已知的干扰数大于阵元数时的波束图。
具体实施方式
下面将结合附图对本发明做进一步详细说明。
本发明的目的在于提供一种半虚拟天线波束形成方法,同虚拟天线一样,既获得比自适应天线更好的性能,又能抑制超过天线阵自由度的干扰。
本发明的目的是这样实现的:
本发明包括如下步骤:
(1)采用单位矩阵I与阵列接收数据的协方差矩阵R做克罗内克积,得到协方差扩展矩阵
(2)通过克罗内克积运算扩展导向矢量
(3)计算基于协方差矩阵扩展的Capon(CME-Capon)波束形成器的最优权为:
(4)进行波束赋形,得到阵列天线的输出数据。
其中:
I:表示单位矩阵;
R:表示阵列接收数据的二阶协方差矩阵;
表示克罗内克积运算;
表示协方差扩展矩阵;
a(θ):表示天线阵列的导向矢量;
b(θ):表示半虚拟天线阵的导向矢量;
(·)-1:表示求逆运算;
(·)H:表示共轭转置;
Wopt:表示Capon波束形成器的最优权矢量;
表示CME-Capon波束形成器的最优权矢量;
a(θd):表示自适应天线期望信号导向矢量;
b(θd):表示半虚拟天线阵期望信号导向矢量;
(·)*:表示共轭。
步骤1中,协方差扩展矩阵也可表示为步骤2中,导向矢量也相应表示为
本方法不仅适用于Capon波束形成器的性能改进,还适用于任何期望信号导向矢量已知的自适应阵列天线波束成形性能改进。
本方法不但适用于单用户阵列天线波束成形,也适用于多用户阵列天线波束成形。
本发明涉及的是一种天线阵列的控制方法,提出了一种基于协方差矩阵扩展的半虚拟天线阵波束形成方法,该方法与自适应天线波束形成相比,构造了虚拟阵元,相当于增加了阵元数,能够获得额外的处理增益,极大的提高了天线阵的波束形成性能。
本发明包括如下步骤:首先用一个单位矩阵与阵列的采样协方差矩阵做克罗内克积,得到协方差扩展矩阵,同时将导向矢量通过克罗内克积运算相应扩展,相当于构建了虚拟阵元,将自适应天线变换到半虚拟天线阵,再根据变换后的半虚拟天线阵进行波束形成。与自适应天线波束形成相比,本发明通过构造虚拟阵元,获得了更高的输出信干噪比,还可以构建更多的约束条件,甚至在方向未知的干扰数不超过天线阵自由度的情况下,通过约束条件抑制方向已知的干扰,实现抑制超自由度干扰。
参照图1,本发明的具体实施步骤如下:
具体实施方式:
步骤1:采用单位矩阵I与阵列接收数据的协方差矩阵R做克罗内克积,得到协方差扩展矩阵如:
步骤2:通过克罗内克积运算扩展导向矢量如:
其中a(θ)表示天线阵列导向矢量,b(θ)表示半虚拟天线阵的导向矢量。
步骤3:计算基于协方差矩阵扩展的Capon(CME-Capon)波束形成器的最优权为
其中表示CME-Capon波束形成器的最优权,b(θd)表示半虚拟天线阵的期望信号导向矢量。
步骤4:进行波束赋形,得到阵列天线的输出数据。
本发明的方法采用下述仿真进行验证:
仿真条件:考虑远场窄带情况,自适应天线为N元均匀直线阵,阵元间距d为半波长,各干扰相互独立,信号与干扰也独立,阵元数N=9,1个期望信号入射角为0°,信噪比SNR=0dB,2个干扰入射角分别为-50°和40°,干噪比INR=30dB,快拍数K=200,比较CME-Capon波束形成器和Capon波束形成器的性能。
图2为两种波束形成器的波束图,图3为不同快拍数对应的输出信干噪比曲线,图4为不同输入信噪比对应的输出信干噪比曲线。由图2可知,CME-Capon波束形成器比Capon波束形成器的旁瓣低很多,主瓣能量更集中。由图3可知,CME-Capon波束形成器比Capon波束形成器的输出信干噪比高8dB左右,性能更好。由图4可知,在输入信干噪比较低的情况下,即输入信干噪比低于10dB时,CME-Capon波束形成器的输出信干噪比超过了自适应天线的理论最优值;在输入信干噪比改变的情况下,CME-Capon波束形成器的输出信干噪比始终比Capon波束形成器高8dB左右。
仿真条件:天线阵元数N=4;1个期望信号入射角为0°,信噪比SNR=0dB;5个干扰,其中2个方向已知的干扰入射角为-20°和35°,3个方向未知的干扰入射角为-45°,25°和65°,干噪比INR=30dB;快拍数K=200,比较CME-Capon波束形成器和Capon波束形成器在增加约束条件时的性能。
两种波束形成器均通过约束条件抑制2个方向已知的干扰,其波束图的仿真结果见图5。从图5可以看出,CME-Capon波束形成器通过约束条件抑制2个方向已知的干扰后,仍能抑制3个方向未知的干扰,而Capon波束形成器则不能。这说明CME-Capon波束形成器与Capon波束形成器相比,由于虚拟阵元的存在,能构建更多的约束条件,具有抑制特定超自由度干扰的能力。
Claims (1)
1.一种基于协方差矩阵扩展的半虚拟天线阵波束形成方法,其特征在于:包括以下步骤,
步骤一:将阵列接收数据的协方差矩阵R进行扩展,得到协方差扩展矩阵
步骤二:将天线阵列导向矢量a(θ)进行扩展,得到半虚拟天线阵列的导向矢量
步骤三:计算基于协方差矩阵扩展的Capon波束形成器的最优权为:
其中,b(θd)为半虚拟天线阵的期望信号导向矢量;
步骤四:利用最优权进行波束赋形,得到阵列天线的输出数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510664784.2A CN105306117B (zh) | 2015-10-15 | 2015-10-15 | 一种基于协方差矩阵扩展的半虚拟天线阵波束形成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510664784.2A CN105306117B (zh) | 2015-10-15 | 2015-10-15 | 一种基于协方差矩阵扩展的半虚拟天线阵波束形成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105306117A true CN105306117A (zh) | 2016-02-03 |
CN105306117B CN105306117B (zh) | 2018-08-17 |
Family
ID=55202942
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510664784.2A Active CN105306117B (zh) | 2015-10-15 | 2015-10-15 | 一种基于协方差矩阵扩展的半虚拟天线阵波束形成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105306117B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107621623A (zh) * | 2016-07-13 | 2018-01-23 | 智易科技股份有限公司 | 信号方向检测方法及应用其的波束形成器 |
CN113206696A (zh) * | 2021-02-05 | 2021-08-03 | 西安宇飞电子技术有限公司 | 空域抗干扰赋形装置 |
CN113422629A (zh) * | 2021-06-17 | 2021-09-21 | 长安大学 | 一种协方差矩阵重构自适应波束形成方法及系统 |
CN113569192A (zh) * | 2021-08-05 | 2021-10-29 | 阳光学院 | 一种多相位分级的嵌套阵列天线波束合成方法 |
CN115441918A (zh) * | 2022-08-31 | 2022-12-06 | 重庆邮电大学 | 一种抗快速运动干扰的稳健虚拟阵列变换方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101770022A (zh) * | 2009-12-30 | 2010-07-07 | 南京航空航天大学 | 基于遗传算法的mimo雷达阵列位置误差自校正方法 |
CN102279387A (zh) * | 2011-07-18 | 2011-12-14 | 西安电子科技大学 | Mimo雷达的目标到达角估计方法 |
US20120288021A1 (en) * | 2011-05-09 | 2012-11-15 | Electronics And Telecommunications Research Institute | Secondary transmitter included in cognitive radio communication system and communication method of the secondary transmitter |
CN103744061A (zh) * | 2014-01-15 | 2014-04-23 | 西安电子科技大学 | 基于迭代最小二乘方法的mimo雷达doa估计方法 |
CN103901417A (zh) * | 2014-04-02 | 2014-07-02 | 哈尔滨工程大学 | L型阵列mimo雷达低复杂度空间目标二维角度估计方法 |
CN103983958A (zh) * | 2014-05-16 | 2014-08-13 | 哈尔滨工程大学 | 基于多测量矢量稀疏表示的mimo雷达连续目标角度估计方法 |
CN104270179A (zh) * | 2014-09-12 | 2015-01-07 | 北京理工大学 | 基于协方差重构和导向矢量补偿的自适应波束形成方法 |
-
2015
- 2015-10-15 CN CN201510664784.2A patent/CN105306117B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101770022A (zh) * | 2009-12-30 | 2010-07-07 | 南京航空航天大学 | 基于遗传算法的mimo雷达阵列位置误差自校正方法 |
US20120288021A1 (en) * | 2011-05-09 | 2012-11-15 | Electronics And Telecommunications Research Institute | Secondary transmitter included in cognitive radio communication system and communication method of the secondary transmitter |
CN102279387A (zh) * | 2011-07-18 | 2011-12-14 | 西安电子科技大学 | Mimo雷达的目标到达角估计方法 |
CN103744061A (zh) * | 2014-01-15 | 2014-04-23 | 西安电子科技大学 | 基于迭代最小二乘方法的mimo雷达doa估计方法 |
CN103901417A (zh) * | 2014-04-02 | 2014-07-02 | 哈尔滨工程大学 | L型阵列mimo雷达低复杂度空间目标二维角度估计方法 |
CN103983958A (zh) * | 2014-05-16 | 2014-08-13 | 哈尔滨工程大学 | 基于多测量矢量稀疏表示的mimo雷达连续目标角度估计方法 |
CN104270179A (zh) * | 2014-09-12 | 2015-01-07 | 北京理工大学 | 基于协方差重构和导向矢量补偿的自适应波束形成方法 |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107621623A (zh) * | 2016-07-13 | 2018-01-23 | 智易科技股份有限公司 | 信号方向检测方法及应用其的波束形成器 |
CN107621623B (zh) * | 2016-07-13 | 2020-04-10 | 智易科技股份有限公司 | 信号方向检测方法及应用其的波束形成器 |
CN113206696A (zh) * | 2021-02-05 | 2021-08-03 | 西安宇飞电子技术有限公司 | 空域抗干扰赋形装置 |
CN113206696B (zh) * | 2021-02-05 | 2023-09-29 | 西安宇飞电子技术有限公司 | 空域抗干扰赋形方法 |
CN113422629A (zh) * | 2021-06-17 | 2021-09-21 | 长安大学 | 一种协方差矩阵重构自适应波束形成方法及系统 |
CN113422629B (zh) * | 2021-06-17 | 2023-02-21 | 长安大学 | 一种协方差矩阵重构自适应波束形成方法及系统 |
CN113569192A (zh) * | 2021-08-05 | 2021-10-29 | 阳光学院 | 一种多相位分级的嵌套阵列天线波束合成方法 |
CN113569192B (zh) * | 2021-08-05 | 2024-03-12 | 阳光学院 | 一种多相位分级的嵌套阵列天线波束合成方法 |
CN115441918A (zh) * | 2022-08-31 | 2022-12-06 | 重庆邮电大学 | 一种抗快速运动干扰的稳健虚拟阵列变换方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105306117B (zh) | 2018-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101369014B (zh) | 应用于多输入多输出雷达的双边约束自适应波束形成方法 | |
CN108462521B (zh) | 自适应阵列天线的抗干扰实现方法 | |
CN102830387B (zh) | 一种基于数据预处理的协方差矩阵正交化波束形成方法 | |
CN105137399B (zh) | 基于斜投影滤波的雷达自适应波束形成方法 | |
CN105306117A (zh) | 一种基于协方差矩阵扩展的半虚拟天线阵波束形成方法 | |
CN106972882B (zh) | 基于虚拟域空间功率谱估计的互质阵列自适应波束成形方法 | |
CN107104720B (zh) | 基于协方差矩阵虚拟域离散化重建的互质阵列自适应波束成形方法 | |
CN103414504B (zh) | 一种基于凸组合的自适应波束形成方法 | |
CN102983896B (zh) | 一种投影虚拟天线波束形成方法 | |
CN103312346B (zh) | 一种调零天线 | |
CN105204008A (zh) | 一种基于协方差矩阵扩展的自适应天线波束形成零陷展宽方法 | |
CN105306123A (zh) | 一种抗阵列系统误差的稳健波束形成方法 | |
CN111400919A (zh) | 阵列天线中的低旁瓣波束设计方法 | |
CN105022268A (zh) | 一种线性约束虚拟天线波束形成方法 | |
CN104931937B (zh) | 基于协方差矩阵归一化的子阵级正交投影波束形成方法 | |
CN108880586A (zh) | 一种宽带弱信号增强方法与装置 | |
CN105490763A (zh) | 一种端到端的宽带移动mimo传播信道模型及建模方法 | |
Nicolaescu et al. | Smart antennas for wireless communications systems | |
CN106842135A (zh) | 基于干扰加噪声协方差矩阵重构的自适应波束形成方法 | |
CN104346532B (zh) | 一种mimo雷达降维自适应波束形成方法 | |
US9893788B2 (en) | Node in a wireless communication system with four beam ports and corresponding method | |
CN110138413B (zh) | 一种空时域自适应宽线性降秩波束形成方法 | |
CN105187105B (zh) | 智能天线波束赋形中采用中心天线改善零陷的优化方法 | |
CN104076336A (zh) | 基于雷达折面阵规则数字子阵的和/差波束的形成方法 | |
CN113595604B (zh) | 一种部分连接架构下的多用户毫米波通信波束成形方法 |
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 |