CN111884689B - 一种基于阶梯阵列的波达方向-距离自适应波束形成方法 - Google Patents

一种基于阶梯阵列的波达方向-距离自适应波束形成方法 Download PDF

Info

Publication number
CN111884689B
CN111884689B CN202010662264.9A CN202010662264A CN111884689B CN 111884689 B CN111884689 B CN 111884689B CN 202010662264 A CN202010662264 A CN 202010662264A CN 111884689 B CN111884689 B CN 111884689B
Authority
CN
China
Prior art keywords
array
equivalent
sub
arrival
appointed
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
Application number
CN202010662264.9A
Other languages
English (en)
Other versions
CN111884689A (zh
Inventor
叶中付
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Science and Technology of China USTC
Original Assignee
University of Science and Technology of China USTC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Science and Technology of China USTC filed Critical University of Science and Technology of China USTC
Priority to CN202010662264.9A priority Critical patent/CN111884689B/zh
Publication of CN111884689A publication Critical patent/CN111884689A/zh
Application granted granted Critical
Publication of CN111884689B publication Critical patent/CN111884689B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413MIMO systems
    • H04B7/0426Power distribution
    • H04B7/043Power distribution using best eigenmode, e.g. beam forming or beam steering
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/08Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
    • H04B7/0837Diversity 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/0842Weighted combining
    • H04B7/086Weighted combining using weights depending on external parameters, e.g. direction of arrival [DOA], predetermined weights or beamforming

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于阶梯阵列的波达方向‑距离自适应波束形成方法,该方法首先设计线性子阵列,把若干相同子阵列排列成阶梯阵型,形成阶梯阵列;接着对指定的信号波达方向,每一子阵列进行方向自适应波束形成,形成对指定波达方向上所有信号的接收;然后将每一子阵看成一个等效阵元,构成等效阵列,每一子阵的输出看成是该等效阵元的接收数据,所有等效阵元的接收数据看成是等效阵列的接收数据;最后将指定波达方向上指定接收距离换算为等效阵列的波达方向,利用等效阵列进行该方向的自适应波束形成,等价于完成指定距离上的自适应波束形成。通过阶梯阵列阵型和二次级联自适应波束形成,实现了对指定波达方向和指定距离上空间信号的最佳接收。

Description

一种基于阶梯阵列的波达方向-距离自适应波束形成方法
技术领域
本发明涉及阵列信号处理领域中波束形成研究领域,尤其涉及一种基于阶梯阵列的波达方向-距离自适应波束形成方法。
背景技术
波束形成是阵列信号处理领域中的一个研究方向,波束形成器是一种空间滤波器。它是根据一定准则对每一个阵元接收到的数据进行加权处理,目的是只保留接收到的数据中的期望信号成分,而对干扰和噪声进行有效抑制。也就是要保证阵列的主瓣要对准期望信号的波达方向,并将零陷对准干扰的波达方向。自适应波束形成器是根据阵列接收数据的信号环境自适应地调节每个阵元的加权系数,保证阵列的主瓣始终对准期望信号的波达方向、零陷始终对准干扰的波达方向,在最大程度上接收期望信号的同时最有效地抑制干扰和噪声。Capon自适应波束形成器从理论上被证明在理想情况下是最优的自适应波束形成器,它保证在无失真接收期望信号的前提下、最小化阵列的输出功率来达到对干扰和噪声进行抑制的目的。最小畸变无失真响应波束形成器是另一种有效的自适应波束形成器。Capon自适应波束形成器和最小畸变无失真响应波束形成器都是通过自适应算法实现的。
然而上述的适应波束形成器均假设在指定的波达方向上只存在一个期望信号,然而在实际情况下,可能有多个信号对自于同一个波达方向的多个不同距离,上述的自适应波束形成方法只能同时接收该波达方向上的所有信号,无法区分出每一个信号。
鉴于以上分析,有必要研究新的波达方向-距离自适应波束形成方法,以实现指定波达方向指定距离上的期望信号的自适应接收。
发明内容
本发明公开了一种基于阶梯阵列的波达方向-距离自适应波束形成方法,该方法首先设计线性子阵列,把若干相同子阵列排列成阶梯阵型,形成阶梯阵列;接着对指定的信号波达方向,每一子阵列进行方向自适应波束形成,形成对指定波达方向上所有信号的接收;然后将每一子阵看成一个等效阵元,构成等效阵列,每一子阵的自适应波束形成输出看成是该等效阵元的接收数据,所有等效阵元的接收数据看成是等效阵列的接收数据;最后将指定接收距离换算为等效阵列的波达方向,利用等效阵列进行该方向的自适应波束形成,等价于完成指定距离上的自适应波束形成。通过阶梯阵列阵型和二次级联自适应波束形成,实现了对指定波达方向和指定距离上空间信号的最佳接收。
本发明的目的是通过以下技术方案实现的:一种基于阶梯阵列的波达方向-距离自适应波束形成方法,包括如下步骤:
步骤1、设计线性子阵列,把若干相同结构子阵列排列成阶梯结构,形成阶梯阵列;
步骤2、对指定的信号波达方向,每一子阵列进行方向自适应波束形成,形成对指定波达方向上所有信号的接收;
步骤3、将每一子阵看成是一个等效阵元,构成等效阵列,每一子阵的自适应波束形成输出看成是该等效阵元的接收数据,所有等效阵元的接收数据看成是等效阵列的接收数据;
步骤4、将指定波达方向上指定接收距离换算为等效阵列的波达方向,利用等效阵列进行该方向的自适应波束形成,等价于完成指定距离上的自适应波束形成。
进一步地,所述步骤1包括以下步骤:
步骤11、设计线性子阵列,子阵列的阵元数为M,相邻阵元的最小间距为d;可以采用均匀线阵、中心对称线阵和非均匀线阵,子阵列孔径为Dsub
步骤12、将L个相同结构线性子阵列,按阶梯结构排列,形成阶梯阵列,各个子阵列中心连线与线性子阵列的夹角为ψ。
进一步地,所述步骤2包括以下步骤:
步骤21、各个子阵列接收数据分别为xl(k),l=1,2,...,L,计算各个阵列接收数据的协方差矩阵
Figure BDA0002579040610000021
l=1,2,...,L;并按下式计算它们的平均值:
Figure BDA0002579040610000022
步骤22、给定需要接收的期望信号的波达方向θ0,按子阵列结构生成期望信号的导向矢量a0,计算出子阵列最佳权矢量
Figure BDA0002579040610000023
步骤23、对各个子阵列接收数据xl(k),l=1,2,...,L,用子阵列最佳权矢量wsub加权处理,各个子阵列波束形成器的输出信号分别为:
Figure BDA0002579040610000024
不是来自于波达方向θ0的其他信号被当成干扰均被零陷滤除或抑制。
进一步地,所述步骤3包括以下步骤:
步骤31、将每一子阵看成一个等效阵元,所有等效阵元组成一个等效的均匀线阵,等效阵列相邻阵元间距为两个相邻子阵列中心距离为dsub,阶梯阵列的等效孔径为D=(L-1)dsub
步骤32、等效阵元的接收数据分别为:
yl(k),l=1,2,...,L,
其中,y(k)=[y1(k),y2(k),...,yL(k)]T为等效阵列的接收数据,等效阵列接收数据的协方差矩阵为
Figure BDA0002579040610000031
进一步地,所述步骤4包括以下步骤:
步骤41、将波达方向θ0上指定接收距离r0换算为等效阵列的波达方向
Figure BDA0002579040610000034
步骤42、对于等效阵列的波达方向
Figure BDA0002579040610000035
按等效阵列生成期望信号的导向矢量b0,计算出等效阵列最佳权矢量
Figure BDA0002579040610000032
步骤43、将最佳权矢量w对等效阵列接收数据y(k)加权,等效阵列波束形成器的输出信号为:
z(k)=wHy(k),
不是来自于波达方向
Figure BDA0002579040610000033
的其他信号被当成干扰均被零陷滤除或抑制。经过二次级联自适应波束形成,等价于完成了指定波达方向θ0、距离r0上的期望信号的自适应接收。
本发明与现有技术相比的优点在于:
(1)本发明针对的是一个波达方向上存在多个信号的空间信号接收问题。
(2)本发明根据上述问题,设计了阶梯阵列阵型。
(3)本发明基于阶梯阵列阵型,采用二次级联的自适应波束形成,首先通过线性子阵列实现对自于波达方向θ0多有信号的接收,然后将波达方向θ0上指定接收距离r0换算为等效阵列的波达方向
Figure BDA0002579040610000036
通过等效阵列实现对自于波达方向
Figure BDA0002579040610000037
信号的接收,从而达到了对指定波达方向和指定距离上空间信号的最佳接收。
附图说明
图1为本发明一种基于阶梯阵列的波达方向-距离自适应波束形成方法的流程图;
图2为线性子阵列的示意图;
图3为阶梯阵列的示意图。
具体实施方式
下面结合附图以及具体实施方式进一步说明本发明。
如图1所示,一种基于阶梯阵列的波达方向-距离自适应波束形成方法,包括如下步骤:
步骤1、设计线性子阵列,把若干相同结构子阵列排列成阶梯结构,形成阶梯阵列;所述步骤1具体包括以下步骤:
步骤11、设计线性子阵列,子阵列的阵元数为M,相邻阵元的最小间距为d;可以采用均匀线阵、中心对称线阵和非均匀线阵,子阵列孔径为Dsub
步骤12、将L个相同结构线性子阵列,按阶梯结构排列,形成阶梯阵列,各个子阵列中心连线与线性子阵列的夹角为ψ。
步骤2、对指定的信号波达方向,每一子阵列进行方向自适应波束形成,形成对指定波达方向上所有信号的接收;所述步骤2具体包括以下步骤:
步骤21、各个子阵列接收数据分别为xl(k),l=1,2,...,L,计算各个阵列接收数据的协方差矩阵
Figure BDA0002579040610000041
l=1,2,...,L;并按下式计算它们的平均值:
Figure BDA0002579040610000042
步骤22、给定需要接收的期望信号的波达方向θ0,按子阵列结构生成期望信号的导向矢量a0,计算出子阵列最佳权矢量
Figure BDA0002579040610000043
步骤23、对各个子阵列接收数据xl(k),l=1,2,...,L,用子阵列最佳权矢量wsub加权处理,各个子阵列波束形成器的输出信号分别为:
Figure BDA0002579040610000044
不是来自于波达方向θ0的其他信号被当成干扰均被零陷滤除或抑制。
步骤3、将每一子阵看成是一个等效阵元,构成等效阵列,每一子阵的自适应波束形成输出看成是该等效阵元的接收数据,所有等效阵元的接收数据看成是等效阵列的接收数据;所述步骤3具体包括以下步骤:
步骤31、将每一子阵看成一个等效阵元,所有等效阵元组成一个等效的均匀线阵,等效阵列相邻阵元间距为两个相邻子阵列中心距离为dsub,阶梯阵列的等效孔径为D=(L-1)dsub
步骤32、等效阵元的接收数据分别为:
yl(k),l=1,2,...,L,
其中,y(k)=[y1(k),y2(k),...,yL(k)]T为等效阵列的接收数据,等效阵列接收数据的协方差矩阵为
Figure BDA0002579040610000051
步骤4、将指定接收距离换算为等效阵列的波达方向,利用等效阵列进行该方向的自适应波束形成,等价于完成指定距离上的自适应波束形成。所述步骤4具体包括以下步骤:
步骤41、将波达方向θ0上指定接收距离r0换算为等效阵列的波达方向
Figure BDA0002579040610000053
步骤42、对于等效阵列的波达方向
Figure BDA0002579040610000054
按等效阵列生成期望信号的导向矢量b0,计算出等效阵列最佳权矢量
Figure BDA0002579040610000052
步骤43、将最佳权矢量w对等效阵列接收数据y(k)加权,等效阵列波束形成器的输出信号为:
z(k)=wHy(k),
不是来自于波达方向
Figure BDA0002579040610000055
的其他信号被当成干扰均被零陷滤除或抑制。经过二次级联自适应波束形成,等价于完成了指定波达方向θ0、距离r0上的期望信号的自适应接收。
对于来自于同一个波达方向的多个信号,以往的自适应波束形成方法是同时接收的。由上述本发明提供的技术方案可以看出,通过设计线性子阵列,把若干相同结构子阵列排列成阶梯阵列,通过子阵列的方向自适应波束形成可以接收到同一个波达方向的多个信号,然后由子阵构成的等效阵列的方向自适应波束形成完成了指定距离上的自适应波束形成。这种基于阶梯阵列阵型和二次自适应波束形成的算法流程,解决了指定波达方向-距离上的期望信号的自适应接收问题。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到上述实施例可以通过软件实现,也可以借助软件加必要的通用硬件平台的方式来实现。基于这样的理解,上述实施例的技术方案可以以软件产品的形式体现出来,该软件产品可以存储在一个非易失性存储介质(可以是CD-ROM,U盘,移动硬盘等)中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述的方法。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书的保护范围为准。

Claims (1)

1.一种基于阶梯阵列的波达方向-距离自适应波束形成方法,其特征在于:包括如下步骤:
步骤1、设计线性子阵列,把若干相同结构子阵列排列成阶梯结构,形成阶梯阵列;
步骤2、对指定的信号波达方向,每一子阵列进行方向自适应波束形成,形成对指定波达方向上所有信号的接收;
步骤3、将每一子阵看成是一个等效阵元,构成等效阵列,每一子阵的自适应波束形成输出看成是该等效阵元的接收数据,所有等效阵元的接收数据看成是等效阵列的接收数据;
步骤4、将指定波达方向上指定接收距离换算为等效阵列的波达方向,利用等效阵列进行该方向的自适应波束形成,等价于完成指定距离上的自适应波束形成;
其中,所述步骤1包括以下步骤:
步骤11、设计线性子阵列,子阵列的阵元数为M,相邻阵元的最小间距为d;可以采用均匀线阵、中心对称线阵和非均匀线阵,子阵列孔径为Dsub
步骤12、将L个相同结构线性子阵列,按阶梯结构排列,形成阶梯阵列,各个子阵列中心连线与线性子阵列的夹角为ψ;
所述步骤2包括以下步骤:
步骤21、各个子阵列接收数据分别为xl(k),l=1,2,...,L,计算各个阵列接收数据的协方差矩阵
Figure FDA0003069115450000011
并按下式计算它们的平均值:
Figure FDA0003069115450000012
步骤22、给定需要接收的期望信号的波达方向θ0,按子阵列结构生成期望信号的导向矢量a0,计算出子阵列最佳权矢量
Figure FDA0003069115450000013
步骤23、对各个子阵列接收数据xl(k),l=1,2,...,L,用子阵列最佳权矢量wsub加权处理,各个子阵列波束形成器的输出信号分别为:
Figure FDA0003069115450000014
不是来自于波达方向θ0的其他信号被当成干扰均被零陷滤除或抑制;
所述步骤3包括以下步骤:
步骤31、将每一子阵看成一个等效阵元,所有等效阵元组成一个等效的均匀线阵,等效阵列相邻阵元间距为两个相邻子阵列中心距离为dsub,阶梯阵列的等效孔径为D=(L-1)dsub
步骤32、等效阵元的接收数据分别为:
yl(k),l=1,2,...,L,
其中,y(k)=[y1(k),y2(k),...,yL(k)]T为等效阵列的接收数据,等效阵列接收数据的协方差矩阵为
Figure FDA0003069115450000021
所述步骤4包括以下步骤:
步骤41、将波达方向θ0上指定接收距离r0换算为等效阵列的波达方向
Figure FDA0003069115450000022
步骤42、对于等效阵列的波达方向
Figure FDA0003069115450000023
按等效阵列生成期望信号的导向矢量b0,计算出等效阵列最佳权矢量
Figure FDA0003069115450000024
步骤43、将最佳权矢量w对等效阵列接收数据y(k)加权,等效阵列波束形成器的输出信号为:
z(k)=wHy(k),
不是来自于波达方向
Figure FDA0003069115450000025
的其他信号被当成干扰均被零陷滤除或抑制,经过二次级联自适应波束形成,等价于完成了指定波达方向θ0、距离r0上的期望信号的自适应接收。
CN202010662264.9A 2020-07-10 2020-07-10 一种基于阶梯阵列的波达方向-距离自适应波束形成方法 Active CN111884689B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010662264.9A CN111884689B (zh) 2020-07-10 2020-07-10 一种基于阶梯阵列的波达方向-距离自适应波束形成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010662264.9A CN111884689B (zh) 2020-07-10 2020-07-10 一种基于阶梯阵列的波达方向-距离自适应波束形成方法

Publications (2)

Publication Number Publication Date
CN111884689A CN111884689A (zh) 2020-11-03
CN111884689B true CN111884689B (zh) 2021-10-22

Family

ID=73150083

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010662264.9A Active CN111884689B (zh) 2020-07-10 2020-07-10 一种基于阶梯阵列的波达方向-距离自适应波束形成方法

Country Status (1)

Country Link
CN (1) CN111884689B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106295142A (zh) * 2016-07-29 2017-01-04 电子科技大学 一种基于概率约束的鲁棒Capon波束形成方法
CN107092007A (zh) * 2017-05-25 2017-08-25 电子科技大学 一种虚拟二阶阵列扩展的波达方向估计方法
CN108181507A (zh) * 2017-12-25 2018-06-19 中国科学技术大学 一种稳健自适应波束形成方法
CN110749855A (zh) * 2019-09-10 2020-02-04 杭州电子科技大学 一种基于协方差域零化的均匀线阵波达方向估计方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10386447B2 (en) * 2015-09-16 2019-08-20 Qatar University Method and apparatus for simple angle of arrival estimation
CN111366892A (zh) * 2020-03-24 2020-07-03 西北工业大学 一种基于神经网络Massive MIMO DOA系统及实现方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106295142A (zh) * 2016-07-29 2017-01-04 电子科技大学 一种基于概率约束的鲁棒Capon波束形成方法
CN107092007A (zh) * 2017-05-25 2017-08-25 电子科技大学 一种虚拟二阶阵列扩展的波达方向估计方法
CN108181507A (zh) * 2017-12-25 2018-06-19 中国科学技术大学 一种稳健自适应波束形成方法
CN110749855A (zh) * 2019-09-10 2020-02-04 杭州电子科技大学 一种基于协方差域零化的均匀线阵波达方向估计方法

Also Published As

Publication number Publication date
CN111884689A (zh) 2020-11-03

Similar Documents

Publication Publication Date Title
EP1429581B1 (en) Method of broadband constant directivity beamforming for non linear and non axi-symmetric sensor arrays embedded in a obstacle
CN113032989B (zh) 一种约束方向性系数的最小化阵元数目的稀布阵优化方法
CN106650104B (zh) 考虑互耦效应的宽带非频变稀疏阵列综合方法
CN106972882B (zh) 基于虚拟域空间功率谱估计的互质阵列自适应波束成形方法
CN107302391A (zh) 基于互质阵列的自适应波束成形方法
CN107104720A (zh) 基于协方差矩阵虚拟域离散化重建的互质阵列自适应波束成形方法
Le Son On the design of sparse arrays with frequency-invariant beam pattern
CN115085784B (zh) 一种分级加权稀疏子阵天线阵列稳健波束合成方法
CN111988078B (zh) 一种基于三维台阶阵列的方向-距离自适应波束形成方法
CN111884689B (zh) 一种基于阶梯阵列的波达方向-距离自适应波束形成方法
CN104539331A (zh) 一种基于改进的混合入侵杂草算法阵列天线波束合成方法
CN111487594A (zh) 一种基于粒子群优化的圆阵列波束形成方法
Elayaperumal et al. Optimal irregular subarray design for adaptive jammer suppression in phased array radar
Sun et al. Robust spherical microphone array beamforming with multi-beam-multi-null steering, and sidelobe control
CN115833894B (zh) 基于子阵的数字模拟合成自适应抗干扰方法
CN116718999A (zh) 一种数字子阵划分方法、系统、设备及介质
CN114510879B (zh) 一种任意间隔直线阵列端射增强的相位配置方法
Sun et al. On beampattern design for beamspace music
CN113255119B (zh) 基于非规则阵元重构的网络化雷达低栅副瓣波束优化方法
WO2007063298A1 (en) Spatial array
CN115642944A (zh) 一种基于序列凸近似法的鲁棒自适应波束形成方法和系统
Bailey et al. Design of digital beamforming subarrays for a multifunction radar
CN107167809A (zh) 一种基于信号子空间聚焦的宽带阻塞阵波束形成方法
CN110208830B (zh) 一种基于空时二维稀疏阵列的导航抗干扰方法
CN115133291A (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