CN115833894B - 基于子阵的数字模拟合成自适应抗干扰方法 - Google Patents
基于子阵的数字模拟合成自适应抗干扰方法 Download PDFInfo
- Publication number
- CN115833894B CN115833894B CN202310123138.XA CN202310123138A CN115833894B CN 115833894 B CN115833894 B CN 115833894B CN 202310123138 A CN202310123138 A CN 202310123138A CN 115833894 B CN115833894 B CN 115833894B
- Authority
- CN
- China
- Prior art keywords
- subarray
- level
- interference
- vector
- digital
- 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
- 238000000034 method Methods 0.000 title claims abstract description 38
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 13
- 238000003786 synthesis reaction Methods 0.000 title claims abstract description 13
- 239000013598 vector Substances 0.000 claims abstract description 74
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 36
- 238000004088 simulation Methods 0.000 claims abstract description 22
- 241000712899 Lymphocytic choriomeningitis mammarenavirus Species 0.000 claims abstract description 17
- 238000005070 sampling Methods 0.000 claims abstract description 10
- 238000001228 spectrum Methods 0.000 claims abstract description 6
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 39
- 230000003044 adaptive effect Effects 0.000 claims description 15
- 238000000354 decomposition reaction Methods 0.000 claims description 4
- 230000003595 spectral effect Effects 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 3
- 230000004721 adaptive immunity Effects 0.000 claims 2
- 230000002452 interceptive effect Effects 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 6
- 238000005259 measurement Methods 0.000 abstract description 3
- 238000012545 processing Methods 0.000 description 9
- 238000003491 array Methods 0.000 description 5
- 230000001629 suppression Effects 0.000 description 5
- 208000032365 Electromagnetic interference Diseases 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开了一种基于子阵的数字模拟合成自适应抗干扰方法,包括:1)将N元线性阵不重叠非均匀划分为L个子阵;2)对接收信号采样,进行空间谱估计,得到空间干扰的个数及其导向矢量;3)对空间干扰的个数进行检测与判断;4)若干扰个数大于子阵个数的一半时,采用LCMV算法,得到阵元级模拟权矢量,进而配置阵元级移相器和衰减器;5)若干扰个数小于或等于子阵个数一半时,利用Capon波束形成算法,得到子阵级数字权矢量,进而配置子阵级数字移相器和衰减器;6)合成波束输出。本发明提出的方法不仅测角精度高,抗干扰能力强,还可以减少硬件资源的消耗,提高计算的收敛速度。
Description
技术领域
本发明涉及阵列信号处理技术领域,特别是一种基于子阵的数字模拟合成自适应抗干扰方法。
背景技术
阵列天线常被称为“相控阵”,即相位控制阵列的简称。顾名思义,相控阵天线是由许多辐射单元排列而成,各个单元的馈电相位由计算机灵活控制的阵列。阵列天线通常是由多个按一定规则排列的全向天线排列而成,通过改变阵列中各个天线单元的幅度和相位关系,可以实现对阵列天线方向图的波束指向的控制。
在电子战环境中,雷达系统不可避免地受到各种无源和有源的电磁干扰,这些干扰信号通过接收天线方向图的旁瓣或主瓣进入接收系统,使接收机的信号检测能力大为降低。雷达的抗干扰性能已成为衡量其性能的重要指标,阵列天线技术是提高雷达抗干扰能力的重要手段。阵列天线的自适应数字波束形成技术(ADBF)使雷达根据干扰特性,自动采取相应的抗干扰方案,并能随其变化动态地调整工作参数,达到某种最佳性能。目前阵列天线的ADBF技术已成为雷达信号处理研究的一个重要内容。
近些年提出的各种波束形成方法,其应用大多是基于阵元级的,这对于只有几个天线阵元的小型阵列来说是可行的。但现在的阵列天线都朝着大型的稀疏阵列的方向发展,特别在相控阵雷达中,几百乃至上千个天线阵元已是屡见不鲜。如果在这种大型阵列中,在阵元级上应用波束形成方法,就得对每一个天线阵元的接收信号进行单独处理,即每一个阵元组成一个接收通道,每一个接收通道都必须包含几次放大、混频,最后是视频处理或模数(A/D)转换。可想而知,硬件成本将会成倍增加。因此必须减少信号处理器的维数,同时尽可能保持对阵列响应的良好控制。
减少信号处理器的维数的一般方法是采用基于子阵的自适应阵列信号处理。
对于一个N元自适应阵列,有N个自由度可以利用,其中一部分自由度用来满足特定的约束,以利用某些可以得到的先验信息,称为约束自由度;剩下的自由度用来自适应地抑制干扰和噪声,称为自适应自由度。自适应自由度越多,自适应算法的计算量越大,收敛速度越慢;反之,收敛速度加快。基于子阵的自适应阵列信号处理正是减少了系统的自适应自由度来加快收敛速度的,因此缺点就是可以自适应抑制干扰的能力下降了。
发明内容
本发明的目的在于提供一种基于子阵的数字模拟合成自适应抗干扰方法。
实现本发明目的的技术解决方案为:一种基于子阵的数字模拟合成自适应抗干扰方法,包括以下步骤:
步骤二:通过对子阵接收信号采样,在设定的快拍数下估计子阵级协方差矩阵,对子阵级协方差矩阵进行特征分解,得到分解后的干扰子空间和噪声子空间,并对该结果采用MUSIC算法进行空间谱估计,得到空间干扰的个数M和空间干扰导向矢量;
步骤四:由期望信号导向矢量和步骤二中空间干扰导向矢量构建约束矩阵C和约束响应向量f,采用LCMV波束形成算法,计算得到阵元级模拟权矢量,并根据该阵元级模拟权矢量配置各天线阵元的移相器和衰减器,返回步骤二;
步骤五:在子阵级采用Capon波束形成算法计算得到子阵级数字权矢量,并根据该子阵级数字权矢量配置各子阵通道的移相器和衰减器;
步骤六:将步骤五中各子阵通道的信号合成为一路波束作为输出。
进一步地,步骤一中使用不重叠子阵划分降维矩阵T:
进一步地,步骤二中采用如下方法计算子阵级协方差矩阵并得到干扰子空间和噪声子空间:
进一步地,步骤二中采用如下方法得到空间干扰导向矢量:
进一步地,步骤四中由期望信号导向矢量和空间干扰导向矢量构建约束矩阵C,约束响应向量f为:
进一步地,步骤五中在子阵级采用Capon波束形成算法得到子阵级数字权矢量的方法为:
进一步地,步骤六中得到的合成波束输出表达式如下:
本发明与现有方法相比,其显著优点是:(1)通过在子阵级对干扰信号个数进行循环检测与判断,并配置阵元级模拟权矢量,可以逐步减少输入到子阵级的干扰个数,最终可以抑制干扰的个数与全自适应阵列处理个数相当,从而相比常规的基于子阵的波束形成方法能够处理更多的干扰;(2)相比常规的全自适应阵列处理波束形成方法可以大大减少硬件资源的消耗,提高了计算收敛速度,适用于拥有成千上万个阵元天线的相控阵雷达。(3)该方法结合了在阵元级进行模拟波束形成和子阵级进行数字波束形成的优点,使得该系统相比普通的基于子阵的自适应波束形成方式抗干扰能力更强,测角精度更高。
附图说明
图1为本发明实施方式的信号处理流程图。
图2是在阵元级使用常规波束形成,在子阵级分别采用LCMV算法和Capon算法的天线方向图仿真结果。
图3为全阵列使用LCMV算法的天线方向图仿真结果。
图4为在阵元级分别使用LCMV算法和常规波束形成,在子阵级采用Capon算法的天线方向图仿真结果。
具体实施方式
下面结合附图和实施例对本发明的方法作进一步描述。
本发明的基本方案包括如下六个步骤,具体的信号处理流程如图1所示。
步骤二:通过对子阵接收信号采样,在设定的快拍数下估计子阵级协方差矩阵,对子阵级协方差矩阵进行特征分解,得到分解后的干扰子空间和噪声子空间,并对该结果采用MUSIC算法进行空间谱估计,得到空间干扰的个数M和空间干扰导向矢量;
步骤四:由期望信号导向矢量和步骤二中空间干扰导向矢量构建约束矩阵C和约束响应向量f,采用LCMV波束形成算法,计算得到阵元级模拟权矢量,并根据该阵元级模拟权矢量配置各天线阵元的移相器和衰减器,返回步骤二;
步骤五:在子阵级采用Capon波束形成算法计算得到子阵级数字权矢量,并根据该子阵级数字权矢量配置各子阵通道的移相器和衰减器;
步骤六:将步骤五中各子阵通道的信号合成为一路波束作为输出。
基于上述基本方案,其中步骤一划分子阵的模型为:
先建立信号模型,设有一线性阵列具有阵元数为N,阵元之间的距离均为d,将该线性阵列划分为L个子阵,形成了L个接收通道。现假设存在P个互不相关的干扰信号,干扰信号的复包络为,干扰信号的导向矢量矩阵为,背景噪声为/>,则阵元接收信号/>。
子阵划分基于降维矩阵,其中/>为加权系数对角阵,/>为/>的子阵形成矩阵,在其第1列的所有元素中,与第l个子阵的阵元序号对应的元素值为1,其余均为0,其中l为1-L之间的所有正整数,则子阵接收信号为/>。
其中:
基于上述基本方案,其中步骤二计算接收阵列的子阵级协方差矩阵和特征分解的方法为:
基于上述基本方案,其中步骤二采用MUSIC算法得到空间干扰的个数M及其导向矢量的方法为:
理论上,信号空间与噪声空间是正交的,因此信号空间方向矩阵A的列向量与噪声空间的特征向量也是正交的,而A的列向量与信号的到达方向一一对应,根据此性质,便可以利用噪声空间的特征向量来求解信号的到达方向。首先构建噪声矩阵即步骤二中得到的噪声子空间/>,然后定义空间谱函数/>:
基于上述基本方案,在进行完步骤二之后,得到了估计的空间干扰的导向矢量,在步骤四中为了得到阵元级模拟权矢量,需要构建约束矩阵和约束响应向量,使得算法能够自适应地抑制干扰信号并保持在波束指向/>的增益为常数,本实施方式构建LCMV算法流程如下:
根据该阵元级模拟权矢量配置各天线阵元的移相器和衰减器;这样虽然上述步骤是在子阵级进行采样估计的,但是通过阵元级的配置可以带来不少的优点,如既能实现基于子阵方式的设备量较少,计算量少的优点又能够提高该系统的自由度,即可以提高由于子阵划分造成的对干扰抑制的自由度不足的问题。该方法经过一次循环能够抑制子阵个数级别的干扰,经过若干次循环后计算得到的阵元级模拟权矢量,最终能够抑制阵元个数级别的干扰,克服了其它基于子阵的自适应阵列信号处理方法可用自由度变小的缺点。
这样在空间干扰数目过多时,通过若干次循环而配置阵元级模拟权矢量,使得大部分的干扰得到了抑制,最后一次循环时子阵采样估计出的干扰已经很少/>,可以在子阵级进行干扰抑制,从而结合了在阵元级进行模拟波束形成和子阵级进行数字波束形成的优点,使得该系统相比普通的基于子阵的自适应波束形成方式抗干扰能力更强,测角精度更高。
进行阵列方向图的仿真,仿真采用均匀线阵,仿真参数如表1所示。
表1仿真参数设置
图2是在阵元级使用常规波束形成,在子阵级分别采用LCMV算法和Capon算法的天线方向图仿真结果;
图3为全阵列使用LCMV算法的天线方向图仿真结果;
图4为在阵元级分别使用LCMV算法和常规波束形成,在子阵级采用Capon算法的天线方向图仿真结果。
由图2可以看到LCMV算法与Capon算法有各自的优点,由仿真的结果图可以看出LCMV算法由于约束条件更精细,所以对干扰的抑制效果更好;而本发明设计的系统在阵元级使用LCMV算法得到阵元级模拟权矢量,在模拟部分就可以精确地将大部分干扰抑制掉,而在子阵级采取Capon算法得到子阵级数字权矢量,相当于配合模拟部分抑制掉剩余的干扰,具有更强大的干扰抑制能力。
由图3可以看到,虽然全阵列的LCMV波束形成方式具有很强的抗干扰能力,但是基于全阵列的实现方式,一般会造成设备量大,成本高等缺点。
由图4可以看到相比于一般的基于子阵的自适应波束形成方法,如阵元CBF+子阵Capon,对于20°方向的干扰没有抑制能力,而本系统的基于子阵的数字模拟合成自适应抗干扰方法明显具有更强的抗干扰能力。
Claims (7)
1.一种基于子阵的数字模拟合成自适应抗干扰方法,其特征在于,包括如下步骤:
步骤二:通过对子阵接收信号采样,在设定的快拍数下估计子阵级协方差矩阵,对子阵级协方差矩阵进行特征分解,得到分解后的干扰子空间和噪声子空间,并对分解后的干扰子空间和噪声子空间采用MUSIC算法进行空间谱估计,得到空间干扰的个数M和空间干扰导向矢量;
步骤四:由期望信号导向矢量和步骤二中空间干扰导向矢量构建约束矩阵C和约束响应向量f,采用LCMV波束形成算法,计算得到阵元级模拟权矢量,并根据该阵元级模拟权矢量配置各天线阵元的移相器和衰减器,返回步骤二;
步骤五:在子阵级采用Capon波束形成算法计算得到子阵级数字权矢量,并根据该子阵级数字权矢量配置各子阵通道的移相器和衰减器;在子阵级采用Capon波束形成算法得到子阵级数字权矢量的方法为:
步骤六:将步骤五中各子阵通道的信号合成为一路波束作为输出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310123138.XA CN115833894B (zh) | 2023-02-16 | 2023-02-16 | 基于子阵的数字模拟合成自适应抗干扰方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310123138.XA CN115833894B (zh) | 2023-02-16 | 2023-02-16 | 基于子阵的数字模拟合成自适应抗干扰方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115833894A CN115833894A (zh) | 2023-03-21 |
CN115833894B true CN115833894B (zh) | 2023-06-16 |
Family
ID=85521616
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310123138.XA Active CN115833894B (zh) | 2023-02-16 | 2023-02-16 | 基于子阵的数字模拟合成自适应抗干扰方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115833894B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118625261A (zh) * | 2024-05-29 | 2024-09-10 | 南京理工大学 | 基于人工智能的空时频一体化检测方法及系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113325385A (zh) * | 2021-07-15 | 2021-08-31 | 中国人民解放军火箭军工程大学 | 一种相控阵-mimo雷达模式收发波束形成抗干扰方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103837861B (zh) * | 2014-03-19 | 2016-08-17 | 北京理工大学 | 基于特征子空间的子阵级线性约束自适应波束形成方法 |
CN103942449B (zh) * | 2014-05-05 | 2017-01-25 | 北京理工大学 | 一种基于信源个数估计的特征干扰相消波束形成方法 |
CN106569181A (zh) * | 2016-11-15 | 2017-04-19 | 大连大学 | 基于协方差矩阵重构稳健Capon波束形成的算法 |
CN114563764A (zh) * | 2021-12-03 | 2022-05-31 | 北京无线电测量研究所 | 一种自适应波束生成方法及系统 |
CN115102597A (zh) * | 2022-06-29 | 2022-09-23 | 中国科学技术大学 | 一种稳健的自适应波束形成方法 |
-
2023
- 2023-02-16 CN CN202310123138.XA patent/CN115833894B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113325385A (zh) * | 2021-07-15 | 2021-08-31 | 中国人民解放军火箭军工程大学 | 一种相控阵-mimo雷达模式收发波束形成抗干扰方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115833894A (zh) | 2023-03-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113472371B (zh) | 自适应阵列天线数字波束合成抗干扰处理方法 | |
CN109490820B (zh) | 一种基于平行嵌套阵的二维doa估计方法 | |
CN104991236B (zh) | 一种单基地mimo雷达非圆信号相干源波达方向估计方法 | |
CN109635240A (zh) | 大规模数字阵列零陷展宽自适应波束形成方法 | |
CN107703478B (zh) | 基于互相关矩阵的扩展孔径二维doa估计方法 | |
CN105354171B (zh) | 一种改进特征矢量的投影子空间估计自适应波束合成方法 | |
CN108872930B (zh) | 扩展孔径二维联合对角化doa估计方法 | |
CN107302391A (zh) | 基于互质阵列的自适应波束成形方法 | |
CN115833894B (zh) | 基于子阵的数字模拟合成自适应抗干扰方法 | |
CN110837075A (zh) | 一种低复杂度的极化参数估计跟踪装置及方法 | |
CN102142609A (zh) | 具有低旁瓣特性的子阵级自适应数字波束形成器 | |
CN104539331B (zh) | 一种基于改进的混合入侵杂草算法阵列天线波束合成方法 | |
CN115085784B (zh) | 一种分级加权稀疏子阵天线阵列稳健波束合成方法 | |
CN104868946B (zh) | 自适应加权的子阵级混合mimo‑相控阵系统的干扰抑制方法 | |
CN111817765B (zh) | 一种基于频率约束的广义旁瓣对消宽带波束形成方法 | |
CN111162878B (zh) | 一种基于子阵降维带约束的多域联合抗干扰方法 | |
CN113985358B (zh) | 一种基于相控阵mimo的接收adbf处理方法 | |
CN105024176A (zh) | 子阵级混合mimo-相控阵系统的方向图设计方法 | |
CN110208830B (zh) | 一种基于空时二维稀疏阵列的导航抗干扰方法 | |
CN114265004A (zh) | 一种基于子空间对消的干扰下的目标角度估计方法 | |
CN114371441A (zh) | 虚拟阵列波达方向估计方法、装置、产品及存储介质 | |
CN107064884B (zh) | 基于规则重叠子阵的自适应波束形成方法 | |
Yu et al. | Adaptive beamforming technique for large-scale arrays with various subarray selections | |
CN111308430A (zh) | 基于fpga及多核dsp硬件架构的测向及抗干扰方法及系统 | |
CN111859272A (zh) | 一种面向大规模天线的快速music谱分解方法、装置及设备 |
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 |