CN106599551A - 一种用于阵列天线足球机器人的快速自适应波束形成算法 - Google Patents
一种用于阵列天线足球机器人的快速自适应波束形成算法 Download PDFInfo
- Publication number
- CN106599551A CN106599551A CN201611061197.5A CN201611061197A CN106599551A CN 106599551 A CN106599551 A CN 106599551A CN 201611061197 A CN201611061197 A CN 201611061197A CN 106599551 A CN106599551 A CN 106599551A
- Authority
- CN
- China
- Prior art keywords
- prime
- beta
- sigma
- covariance matrix
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供了一种用于阵列天线足球机器人的快速自适应波束形成算法,其特征在于,包括以下步骤:对训练样本集合进行预处理,剔除期望信号;利用预处理后的数据估计协方差矩阵,并对协方差矩阵的列向量进行RGS正交化构造干扰子空间;将对应的静态权矢量向干扰子空间作正交投影得到自适应权矢量。相较于其它自适应波束形成算法,本发明提出的算法具有波束指向精度高、输出SINR高和对波束指向误差稳健性强等优势,是对期望信号相消现象稳健的一种实用且快速的自适应波束形成算法。
Description
技术领域
本发明涉及一种利用阵列天线协助机器人进行目标识别的算法。
背景技术
当前大多数足球机器人不能在大雾、下雨、黑夜等环境中比赛,为了能让机器人足球比赛24小时无障碍开展,需要利用阵列天线协助机器人进行目标识别。
传统的阵列天线当训练快拍数据不含期望信号时,常规SMI算法对波束指向误差、通道幅相误差等误差不敏感,但对低快拍误差不稳健。低快拍时,常规SMI算法会因噪声子空间对应的小特征值抖动,导致方向图畸变,副瓣电平抬高。而且对于实际雷达***,往往要求自适应波束形成算法具有较低的运算量和较快的收敛速度。为了提高常规自适应波束形成算法对低快拍误差的稳健性和实现的快速性,国内外学者相继提出了许多性能良好的ADBF算法。对角加载(LSMI)算法通过在常规SMI算法的目标函数上添加了对角加载项,减小了噪声扰动的影响,从而克服了低快拍误差。但LSMI方法最大的不足是对角加载量的最优值不易确定。带恶化函数的SMI(PFM-SMI)算法用自适应权值与静态权值之间的误差值构造恶化函数,在保证阵列输出SINR最大的同时,满足某些二次约束,进而求得最优权矢量。PFM-SMI算法在一定程度上避免了低快拍误差导致的旁瓣抬高,但该方法运算量大,且存在参数选择问题。正交投影算法通过干扰噪声协方差矩阵进行特征值分解构造干扰信号空间,然后将约束导向矢量向干扰子空间作正交投影,获得自适应权矢量。由于通过干扰子空间求取自适应权,因此该算法避免了噪声扰动对自适应权矢量的影响,不会使自适应波束图副瓣抬高。但该算法由于需要特征值分解,因此运算量较大,同时在构造干扰子空间时存在干扰个数的确定问题。
上述这些稳健的ADBF算法或多或少都存在一些不足,且都存在一个快速实现的问题。近年来,Hung等提出的Gram-Schmidt(GS)正交化算法是一种快速子空间投影算法,该算法在较高的干噪比下能够较好地重构干扰子空间,具有子空间投影算法的快速收敛特性,而且其运算复杂度小,便于工程实现,尤其适用于阵元数较多、干扰源数较少及强干扰源的场合,从而受到广泛的关注。
发明内容
本发明的目的是提出一种基于GS正交化的快速自适应波束形成算法。
为了达到上述目的,本发明的技术方案是提供了一种用于阵列天线足球机器人的快速自适应波束形成算法,其特征在于,包括以下步骤:
步骤1、对训练样本集合进行预处理,剔除期望信号,假设k时刻阵列接收的快拍数据矢量为x(k),则x(k)中第l个数据分量xl可以表示为,l=1,2,...,N:
式(1)中,si表示第i个导向矢量,βl,i表示第l个快拍的第i个相位调制,里面包含了DOA信息,nl表示第l个快拍的噪声矢量;
对xl作数据预处理得到xl′,
步骤2、利用预处理后的数据估计协方差矩阵,并对协方差矩阵的列向量进行RGS正交化构造干扰子空间;
步骤3、将对应的静态权矢量向干扰子空间作正交投影得到自适应权矢量。
优选地,所述步骤2包括以下步骤:
步骤2.1、将式(2)写成矩阵形式,有:
x′(k)=Bx(k) (3),式(3)中:
x′(k)=[x1′,x2′,…,xN-1′]T表示预处理数据矢量;
B表示阻塞矩阵
x′(k)的协方差矩阵表示为
步骤2.2、计算得到协方差矩阵GS正交算法的正交化自适应门限Δ′(k)
预处理数据矢量x′(k)的噪声方差近似表示为:
式(6)中,σn表示噪声方差;
协方差矩阵的噪声方差近似为:
则有正交化自适应门限Δ′(k)为:
式中,Ui及Ui′表示对的列矢量进行RGS正交化得到的正交向量,i表示第i次迭代;
步骤2.3、协方差矩阵GS正交算法的正交化过程表示为:
优选地,在所述步骤3中:
协方差矩阵GS正交算法的自适应权矢量为wMRGS,则有:
式(10)中,表示协方差矩阵GS正交算法判定的干扰源个数,wq表示静态权矢量。
相较于其它自适应波束形成算法,本发明提出的算法具有波束指向精度高、输出SINR高和对波束指向误差稳健性强等优势,是对期望信号相消现象稳健的一种实用且快速的自适应波束形成算法。
附图说明
图1为本发明的算法处理流程示意图;
图2(a)为本发明算法和RGS算法在强干扰环境下自适应方向图比较;
图2(b)为本发明算法和RGS算法在弱干扰环境下自适应方向图比较。
具体实施方式
为使本发明更明显易懂,兹以优选实施例,并配合附图作详细说明如下。
针对协方差矩阵中存在期望信号将导致常规RGS算法性能严重下降的问题,提出了基于数据预处理的协方差矩阵GS正交化算法(MRGS算法)。该算法首先对训练样本集合进行预处理,剔除期望信号;然后利用预处理后的数据估计协方差矩阵,并对协方差矩阵的列向量进行RGS正交化构造干扰子空间;最后将对应的静态权矢量向干扰子空间作正交投影得到自适应权矢量。另外,本发明还将针对MRGS算法对预处理数据的正交化自适应门限进行修正,以准确地估计干扰子空间。
假设k时刻阵列接收的快拍数据矢量为x(k),则x(k)中第l个数据分量xl可以表示为,l=1,2,...,N:
式(1)中,si表示第i个导向矢量,βl,i表示第l个快拍的第i个相位调制,里面包含了DOA信息,nl表示第l个快拍的噪声矢量;
对xl作数据预处理得到xl′,
由式(2)可以看出,对xl作数据预处理得到的xl′中只含有干扰和噪声信号成分,预处理过程起到了剔除期望信号的作用。因此,将式(2)写成矩阵形式,有:
x′(k)=Bx(k) (3),式(3)中:
x′(k)=[x1′,x2′,…,xN-1′]T表示预处理数据矢量;
B表示阻塞矩阵
x′(k)的协方差矩阵表示为
只含有干扰和噪声信号成分。对进行RGS正交化时,正交化自适应门限需要加以改进,以准确估计干扰子空间。
预处理数据矢量x′(k)的噪声方差近似表示为:
式(6)中,σn表示噪声方差;
协方差矩阵的噪声方差近似为:
则计算得到MRGS算法的正交化自适应门限Δ′(k)为:
式中,Ui及Ui′表示对的列矢量进行RGS正交化得到的正交向量,i表示第i次迭代;
因此,MRGS算法的正交化过程可以表示如下:
那么,MRGS算法的自适应权矢量wMRGS可计算如下:
式(10)中,表示协方差矩阵GS正交算法判定的干扰源个数,wq表示静态权矢量,一般为期望导向矢量。
图1给出了MRGS算法的处理流程示意图。本发明提出的MRGS算法对训练快拍数据矢量进行数据预处理,解决了常规RGS算法出现的信号相消问题;对估计协方差矩阵的列矢量进行RGS正交化,能够利用更多的样本信息减小噪声扰动影响。因此,在训练样本中含有期望信号的应用场合,本发明提出的MRGS算法确实是一种快速且稳健的自适应波束形成方法。
下面将通过仿真实验验证本发明提出的MRGS算法的性能。仿真中,假设阵元数为16的全向阵元组成等间距均匀线阵,阵元间距d=λ/2,噪声的均方差σn=1。此时,阻塞矩阵B中输出SINR为100次蒙特卡洛实验仿真平均得到,波束图为1次蒙特卡洛实验仿真得到。
实验:MRGS算法和RGS算法自适应方向图比较。假定有1个期望信号和3个干扰信号,来向分别为θ0=0°,θi=-28°、17°、41°,输入信噪比SNR=0dB,快拍数K=20。自适应波束图的仿真结果如图2所示,图2(a)为强干扰情况下的仿真结果,干噪比INR均为30dB;而图2(b)为弱干扰情况下的仿真结果,INR均为10dB。
从图2(a)可以看出,由于期望信号功率相对较弱,因此未被RGS算法当作干扰进行抑制。此时,RGS算法和MRGS算法都能在干扰方向形成零陷,在期望信号方向形成最大波束指向,但MRGS算法的零陷深度比RGS算法至少深15dB,性能更好。从图2(b)可以看出,由于此时期望信号功率相对较强,RGS算法将期望信号当作干扰进行抑制,不仅没有在期望信号方向形成最大波束指向,反而形成了一定深度的零陷,波束图严重恶化;而MRGS算法的波束图保持良好,不会在期望信号方向形成零陷,这主要是数据预处理的作用。可见,在一定的干扰功率范围内,MRGS算法能保持良好的波束性能,而RGS算法在弱干扰环境下性能就会严重恶化。
Claims (3)
1.一种用于阵列天线足球机器人的快速自适应波束形成算法,其特征在于,包括以下步骤:
步骤1、对训练样本集合进行预处理,剔除期望信号,假设k时刻阵列接收的快拍数据矢量为x(k),则x(k)中第l个数据分量xl可以表示为,l=1,2,...,N:
式(1)中,si表示第i个导向矢量,βl,i表示第l个快拍的第i个相位调制,里面包含了DOA信息,nl表示第l个快拍的噪声矢量;
对xl作数据预处理得到xl′,
步骤2、利用预处理后的数据估计协方差矩阵,并对协方差矩阵的列向量进行RGS正交化构造干扰子空间;
步骤3、将对应的静态权矢量向干扰子空间作正交投影得到自适应权矢量。
2.如权利要求1所述的一种用于阵列天线足球机器人的快速自适应波束形成算法,其特征在于,所述步骤2包括以下步骤:
步骤2.1、将式(2)写成矩阵形式,有:
x′(k)=Bx(k) (3),式(3)中:
x′(k)=[x1′,x2′,…,xN-1′]T表示预处理数据矢量;
B表示阻塞矩阵
x′(k)的协方差矩阵表示为
步骤2.2、计算得到协方差矩阵GS正交算法的正交化自适应门限Δ′(k)预处理数据矢量x′(k)的噪声方差近似表示为:
式(6)中,σn表示噪声方差;
协方差矩阵的噪声方差近似为:
则有正交化自适应门限Δ′(k)为:
式中,Ui及Ui′表示对的列矢量进行RGS正交化得到的正交向量,i表示第i次迭代;
步骤2.3、协方差矩阵GS正交算法的正交化过程表示为:
3.如权利要求2所述的一种用于阵列天线足球机器人的快速自适应波束形成算法,其特征在于,在所述步骤3中:
协方差矩阵GS正交算法的自适应权矢量为wMRGS,则有:
式(10)中,表示协方差矩阵GS正交算法判定的干扰源个数,wq表示静态权矢量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611061197.5A CN106599551A (zh) | 2016-11-25 | 2016-11-25 | 一种用于阵列天线足球机器人的快速自适应波束形成算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611061197.5A CN106599551A (zh) | 2016-11-25 | 2016-11-25 | 一种用于阵列天线足球机器人的快速自适应波束形成算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106599551A true CN106599551A (zh) | 2017-04-26 |
Family
ID=58594743
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611061197.5A Pending CN106599551A (zh) | 2016-11-25 | 2016-11-25 | 一种用于阵列天线足球机器人的快速自适应波束形成算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106599551A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107194177A (zh) * | 2017-05-23 | 2017-09-22 | 佛山职业技术学院 | 一种移动分布式波束形成算法 |
CN109462445A (zh) * | 2018-12-12 | 2019-03-12 | 河南宏泰控飞信息技术有限公司 | 一种城市环境中无人机同频段多目标空间分离方法及装置 |
CN112881973A (zh) * | 2021-01-20 | 2021-06-01 | 西北工业大学 | 一种基于rbf神经网络的自修正波束设计方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050157776A1 (en) * | 2004-01-15 | 2005-07-21 | Kil-Hyen Ryu | Method and apparatus for forming a beam |
KR100762219B1 (ko) * | 2006-06-30 | 2007-10-01 | 삼성탈레스 주식회사 | 공간분할 다중접속 방법 |
KR20100052279A (ko) * | 2008-11-10 | 2010-05-19 | 인하대학교 산학협력단 | 코그니티브 라디오 시스템에서의 직교 빔을 이용한 주파수 공유 방법 |
CN102830387A (zh) * | 2012-08-23 | 2012-12-19 | 北京理工大学 | 一种基于数据预处理的协方差矩阵正交化波束形成方法 |
CN103942449A (zh) * | 2014-05-05 | 2014-07-23 | 北京理工大学 | 一种基于信源个数估计的特征干扰相消波束形成方法 |
CN104898093A (zh) * | 2015-05-29 | 2015-09-09 | 西安电子科技大学 | 基于gs算法的mimo雷达正交相位编码信号设计方法 |
-
2016
- 2016-11-25 CN CN201611061197.5A patent/CN106599551A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050157776A1 (en) * | 2004-01-15 | 2005-07-21 | Kil-Hyen Ryu | Method and apparatus for forming a beam |
KR100762219B1 (ko) * | 2006-06-30 | 2007-10-01 | 삼성탈레스 주식회사 | 공간분할 다중접속 방법 |
KR20100052279A (ko) * | 2008-11-10 | 2010-05-19 | 인하대학교 산학협력단 | 코그니티브 라디오 시스템에서의 직교 빔을 이용한 주파수 공유 방법 |
CN102830387A (zh) * | 2012-08-23 | 2012-12-19 | 北京理工大学 | 一种基于数据预处理的协方差矩阵正交化波束形成方法 |
CN103942449A (zh) * | 2014-05-05 | 2014-07-23 | 北京理工大学 | 一种基于信源个数估计的特征干扰相消波束形成方法 |
CN104898093A (zh) * | 2015-05-29 | 2015-09-09 | 西安电子科技大学 | 基于gs算法的mimo雷达正交相位编码信号设计方法 |
Non-Patent Citations (1)
Title |
---|
杨小鹏 等: "基于数据预处理的改进GS正交化波束形成", 《北京理工大学学报》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107194177A (zh) * | 2017-05-23 | 2017-09-22 | 佛山职业技术学院 | 一种移动分布式波束形成算法 |
CN109462445A (zh) * | 2018-12-12 | 2019-03-12 | 河南宏泰控飞信息技术有限公司 | 一种城市环境中无人机同频段多目标空间分离方法及装置 |
CN112881973A (zh) * | 2021-01-20 | 2021-06-01 | 西北工业大学 | 一种基于rbf神经网络的自修正波束设计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103837861B (zh) | 基于特征子空间的子阵级线性约束自适应波束形成方法 | |
CN106569181A (zh) | 基于协方差矩阵重构稳健Capon波束形成的算法 | |
CN105302936A (zh) | 基于相关计算和协方差矩阵重构的自适应波束形成方法 | |
CN103245956B (zh) | 一种基于稳健波束形成算法的gps抗多径方法 | |
CN103984676A (zh) | 一种基于协方差矩阵重构的正交投影自适应波束形成方法 | |
CN102944870A (zh) | 协方差矩阵对角加载的稳健自适应波束形成方法 | |
CN107462872A (zh) | 一种抗主瓣干扰算法 | |
CN109407055A (zh) | 基于多径利用的波束形成方法 | |
CN105137399A (zh) | 基于斜投影滤波的雷达自适应波束形成方法 | |
CN105306123A (zh) | 一种抗阵列***误差的稳健波束形成方法 | |
CN103885045B (zh) | 基于子阵划分的循环联合自适应波束形成方法 | |
CN105182302A (zh) | 一种抗快速运动干扰稳健零陷展宽波束形成方法 | |
CN106443594A (zh) | 一种基于稀疏约束的雷达天线阵列稳健波束形成方法 | |
CN103293517B (zh) | 基于脊参数估计的对角加载稳健自适应雷达波束形成方法 | |
CN107290732B (zh) | 一种量子大***的单基地mimo雷达测向方法 | |
CN106324625A (zh) | 一种基于二范数多目标优化的卫星导航***自适应抗干扰方法 | |
CN105049382A (zh) | 一种抗期望信号导向矢量失配的零陷展宽自适应天线波束成形方法 | |
CN103728601A (zh) | 雷达信号运动干扰空域-极化域联合稳健滤波方法 | |
CN106599551A (zh) | 一种用于阵列天线足球机器人的快速自适应波束形成算法 | |
CN104330766A (zh) | 一种稳健的波达方向估计方法 | |
CN104360337B (zh) | 基于1范数约束的自适应波束形成方法 | |
CN115372925A (zh) | 一种基于深度学习的阵列稳健自适应波束形成方法 | |
CN104346532B (zh) | 一种mimo雷达降维自适应波束形成方法 | |
CN105699988A (zh) | 用于导航接收机的去噪mvdr欺骗干扰抑制方法 | |
CN105068090A (zh) | 一种利用单快拍数据进行gnss天线阵干扰抑制的方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170426 |
|
RJ01 | Rejection of invention patent application after publication |