CN104931935A - 基于旁瓣对消的稳健波束形成方法 - Google Patents
基于旁瓣对消的稳健波束形成方法 Download PDFInfo
- Publication number
- CN104931935A CN104931935A CN201510329662.8A CN201510329662A CN104931935A CN 104931935 A CN104931935 A CN 104931935A CN 201510329662 A CN201510329662 A CN 201510329662A CN 104931935 A CN104931935 A CN 104931935A
- Authority
- CN
- China
- Prior art keywords
- matrix
- vector
- channel
- desired signal
- wanted signal
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/2813—Means providing a modification of the radiation pattern for cancelling noise, clutter or interfering signals, e.g. side lobe suppression, side lobe blanking, null-steering arrays
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radio Transmission System (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
该发明公开了一种基于旁瓣对消的稳健波束形成方法,属于自适应阵列信号处理领域。本发明利用已知期望信号来波方向所在区域来构造阻塞矩阵,并通过重新估计期望信号导向矢量来获得更精确的主、辅助通道加权,从而使旁瓣对消的波束形成器具备针对来波方向失配的稳健性的方法。当期望信号来波方向存在失配时,这种方法具有传统旁瓣对消算法不具备的稳健性,解决了主通道由于加权向量没有对准期望信号的来波方向使得主通道输出中的期望信号增益相较于没有失配时有所下降和利用失配的来波方向构造的阻塞矩阵不能将辅助通道中期望信号完全阻塞,以致主辅助通道对消后期望信号也会在一定程度上被对消这两个问题。
Description
技术领域
本发明属于自适应阵列信号处理领域的波束形成技术,具体是涉及一种基于旁瓣对消的稳健波束形成方法。
背景技术
旁瓣对消技术(sidelobe cancellation)很早就应用到了数字波束形成领域。通过自适应地调整阵列各通道的复加权系数,不仅保证了主瓣输出响应,而且最大限度地抑制了来自各旁瓣方向的非相关强干扰。
与自适应波束置零算法相比,旁瓣对消算法具有相当大的优势,旁瓣对消算法是在自适应置零的基础上发展起来的,并于1997年出现在IEEE标准中。旁瓣对消算法是在阵列中选定一部分阵元或者在阵列的周围添加若干个阵元作为辅助阵元,通过在数字域对辅助通道阵元接收到的信号进行幅度和相位加权,以形成对消信号。由于旁瓣对消算法只对辅助阵元进行幅度和相位加权,所以当阵列为大阵列时,它的优势更加明显。
旁瓣对消***由高增益的主天线(main antenna,波束宽度较窄,通常使用具有高方向性的天线)和低增益辅助天线(auxiliary antenna)阵列构成。通过自适应调整各辅助通道的权系数,在不影响期望信号接收的情况下,尽可能的抑制从主通道旁瓣进入的干扰信号。其***如图1所示,阵元个数为M+N,其中M和N分别为主阵列和辅助阵列包含的阵元个数;主阵列接收信号向量为xm(n),wm为主通道(main channel)的权向量,一般来说固定为期望信号导向矢量,主通道输出为ym(n);辅助阵接收信号向量为xa(n),wa为自适应权向量,自适应滤波器输出为ya(n);Ba为辅助通道(auxiliary channel)阻塞矩阵,用于阻止主瓣方向的信号进入辅助通道;e(n)是对消后的输出。
设期望信号s(n)在主通道的导向矢量为a,各干扰信号在主通道的导向矢量为al(l=1、2…L),期望信号在辅助通道中的导向矢量为af,各干扰信号在辅助通道中的导向矢量为afl(l=1、2…L),主通道接收到的噪声为v(n),辅助通道接收的噪声信号为vf(n)。
则主通道接收到的信号xm(n)是:
辅助通道接收到的信号xa(n)是:
令主辅通道的干扰和噪声分量分别为 和 则主阵列和辅助阵列的接收信号可以分别表示为
xm(n)=as(n)+z(n) (3)
xa(n)=afs(n)+za(n) (4)
在n时刻,xm(n)和xa(n)经权向量wm和wa加权后的输出为
旁瓣对消的原理就是使辅助通道输出ya(n)尽量只包含干扰和噪声,然后用主通道输出ym(n)去减辅助通道输出ya(n),从而得到尽可能大的期望信号响应输出,达到抑制干扰的效果。为使ya(n)中不包含期望信号s(n)的信息,以避免对主瓣的期望响应输出产生不必要的抵消作用,接收信号需要通过辅助通道的期望信号阻塞网络,设计阻塞矩阵Ba,令阻塞矩阵与辅助通道的期望信号导向矢量正交,即:
Baaf=0 (7)
以线阵为例,主通道由M个阵元组成,辅助通道由N个阵元组成,阵元间距为d,期望信号在辅助通道中的导向矢量为:
一般的旁瓣对消算法中,假定期望信号来波方向θ0是已知的,这里λ是信号的波长,因此,期望信号的导向矢量也是已知的。根据期望信号导向矢量的形式,一种常用的阻塞矩阵为:
此时***的输出为:
其中为主通道的期望信号输出成分,为主通道的噪声加干扰输出成分,旁瓣对消算法的原理是让辅助通道的输出尽量接近主通道的干扰加噪声输出,然后用主通道输出与辅助通道输出对消,从而得到尽可能大的期望信号响应输出,也就是让尽可能小,为了尽量抑制输出e(n)中的干扰和噪声,应选择wa使e(n)的平均功率最小,即最小化
令R1=E(xa(n)xm(n)H),R2=E(xa(n)xa(n)H),R3=E(xm(n)xm(n)H),对J(wa)关于wa求梯度有:
令求得最优加权系数wao:
从公式(13)可以看出,辅助通道的最优权向量wao与阻塞矩阵Ba和主通道权向量wm(即期望信号导向矢量a)有关。前面我们提到,在一般的旁瓣对消算法中,假定期望信号来波方向θ0是已知的,而阻塞矩阵Ba和期望信号导向矢量a都是由θ0确定的。在实际应用中,我们所得到的θ0其实是不准确的,那么在辅助通道最优权值的计算中带入由θ0得到的阻塞矩阵Ba和期望信号导向矢量a就会产生问题。即有
Baaf≠0 (14)
此时,干扰对消后的输出为
从该式可以看出,一是主通道由于加权向量没有对准期望信号的来波方向,因此期望信号相较于没有失配时有所下降;二是利用失配的来波方向构造的阻塞矩阵不能将辅助通道中期望信号完全阻塞,以致主辅通道对消后期望信号也会在一定程度上被对消。
发明内容
背景技术中提到,在传统的旁瓣对消算法中,期望信号的来波方向θ0是假定已知的,但是在实际应用中这个角度其实是不准确的,当期望信号来波方向θ0存在失配时,采用常规旁瓣对消技术进行干扰抑制会出现两个问题:一个是主通道由于加权向量没有对准期望信号的来波方向使得主通道输出中的期望信号增益相较于没有失配时有所下降,另一个是利用失配的来波方向构造的阻塞矩阵不能将辅助通道中期望信号完全阻塞,以致主辅助通道对消后期望信号也会在一定程度上被对消。这两者共同作用的结果会使得输出的期望信号增益被大大的消减了。针对这一现象,本发明提出了一种基于旁瓣对消的稳健波束形成方法,利用已知期望信号来波方向所在区域来构造阻塞矩阵,并通过重新估计期望信号导向矢量来获得更精确的主、辅助通道加权,从而使旁瓣对消的波束形成器具备针对来波方向失配的稳健性的方法。
尽管已知的期望信号来波方向θ0不准确,但是我们却可以准确地知道实际期望信号来波方向是存在于一个空域区间的,这个空域区间可以表示为Θ=[θ0-Δθ,θ0+Δθ],这里Δθ是角度的误差量,本发明正是基于这个条件提出的。通过旋转子空间投影的原理,估计期望信号导向矢量以及构造阻塞矩阵。
本发明解决上述技术问题所采用的技术方案包括以下步骤:
a.利用旋转子空间的Capon波束形成算法思路:在已知的来波方向空间范围内积分,针对主通道和辅助通道,分别构造正定矩阵Qm和Qa,并对矩阵进行特征分解,取主特征值对应的特征向量构造矩阵Um和Ua;
b.利用矩阵Ua构造与辅助通道中期望信号导向矢量af正交的阻塞矩阵Ba;
c.利用新构造的阻塞矩阵Ba和矩阵Um得到更加精确的主通道的期望信号导向矢量
d.将得到的期望信号导向矢量和新构造的阻塞矩阵Ba带入到旁瓣对消最优权值的计算式中,得到辅助通道最优的权值wao和主通道权值wm。
本发明的有益效果为,当期望信号来波方向存在失配时,这种方法具有传统旁瓣对消算法不具备的稳健性,解决了主通道由于加权向量没有对准期望信号的来波方向使得主通道输出中的期望信号增益相较于没有失配时有所下降和利用失配的来波方向构造的阻塞矩阵不能将辅助通道中期望信号完全阻塞,以致主辅助通道对消后期望信号也会在一定程度上被对消这两个问题。
附图说明
图1为自适应旁瓣对消***。
具体实施方式
下面结合实施例对本发明做进一步的详细描述:
本发明主要是利用旋转子空间的方法估计的期望信号导向矢量构造辅助通道阻塞矩阵,最终求出最优的权向量。
已知期望信号的来波方向在一个确定的空域区间中,且这个区间内不包括干扰信号。将这个空域区间表示为Θ=[θ0-Δθ,θ0+Δθ],在此区间内分别构造正定矩阵Qm和Qa:
分别对矩阵Qm和Qa进行特征分解,并取其中K个最大的特征值所对应的特征向量umk,k=1,…,K,和uak,k=1,…,K,将这两组特征向量分别构成M×K维的列正交矩阵Um=[um1 um2 … umK]∈CM×K和Ua=[ua1 ua2 … uaK]∈CN×K。这K个最大的特征值包含了矩阵Qm和Qa的所有特征值的大部分能量,任何来自于空域区间Θ的导向矢量都能分别由矩阵Um和Ua的列向量的线性表出。因此,主通道和辅助通道的期望信号的导向矢量可以表示为:
a=Umrm (18)
af=Uara (19)
其中,rm和ra是K×1维的旋转向量。
首先,我们使用Ua构造辅助通道中与期望信号导向矢量正交的阻塞矩阵。令矩阵它的秩为(M-K),在实际应用中根据干扰个数取矩阵T中的任意D行,(其中D≤(M-K)且D不小于干扰源个数),构造阻塞矩阵Ba。显然,辅助通道中估计的期望信号导向矢量与构造的阻塞矩阵Ba是正交的。在辅助通道的自适应权值为公式(13)的前提下,主辅两通道对消后的输出即为想要的期望信号。因此将公式(13)计算所得的最优权值带入公式(11),得到
该式即为对消后的输出信号的输出功率,可近似看为是期望信号的输出功率。主通道权值一般固定为期望信号导向矢量,因此这里用期望信号导向矢量a=Umrm替代主通道权值wm,得到
J(rm)=rm HUm H[R3-R1 HBa H(BaR2Ba H)-1BaR1]Umrm (21)
则可以通过最大化输出期望信号功率(即上式),求得旋转向量rm。为避免模糊,对旋转向量rm施加模约束则最优化问题就可以表示为
其中RBa=R3-R1 HBa H(BaR2Ba H)-1BaR1,利用Lagrange乘数法进行求解,有
其中,μ是Lagrange乘数。对公式(23)求关于向量rm的偏导,并使其等于零,解得rm是RBa最大的特征值所对应的特征向量,令矩阵RBa的最大的特征值所对应的特征向量为将的模进行处理,使其满足约束条件解得主通道估计信号导向矢量为
将估计出的新期望信号导向矢量和重新构造的阻塞矩阵带入到旁瓣对消最优权值的计算式中,得到辅助通道最优的权值。
在实际应用中,R1=E(xa(n)xm(n)H),R2=E(xa(n)xa(n)H),R3=E(xm(n)xm(n)H)是无法获知的,我们进行如下替代:
其中,I是快拍数,那么在实际应用中矩阵RBa被替代为:
辅助通道的最优权值为:
在实际应用中,具体步骤如下:
a.计算列正交矩阵Um和Ua;
在已知的期望信号来波方向空域区间Θ=[θ0-Δθ,θ0+Δθ]内由公式(16)(17)求得正矩阵Qm和Qa,分别对矩阵Qm和Qa进行特征分解,并取其中K个最大的特征值所对应的特征向量umk,k=1,…,K,和uak,k=1,…,K,将这两组特征向量分别构成M×K维的列正交矩阵Um=[um1 um2 … umK]∈CM×K和Ua=[ua1 ua2 … uaK]∈CN×K。
b.构造阻塞矩阵Ba;
令矩阵它的秩为(M-K),根据实际应用中干扰个数取矩阵T中的任意D行,(其中D≤(M-K)且D不小于干扰源个数),构造阻塞矩阵Ba。
c.估计主通道期望信号导向矢量
根据公式(25)(26)(27)(28)计算得到矩阵计算矩阵的最大的特征值所对应的特征向量为由公式(24)得到估计的期望信号导向矢量即主通道权值wm。
d.计算辅助通道最优权值wao。
将上面得到的变量带入公式(29)中得到辅助通道的最优权值。
Claims (1)
1.一种基于旁瓣对消的稳健波束形成方法,其特征在于,包括以下步骤:
a.在已知的来波方向空间范围内积分,针对主通道和辅助通道,分别构造正定矩阵Qm和Qa,并对矩阵进行特征分解,取主特征值对应的特征向量构造矩阵Um和Ua;
b.利用矩阵Ua构造与辅助通道中期望信号导向矢量af正交的阻塞矩阵Ba;
c.利用新构造的阻塞矩阵Ba和矩阵Um得到更加精确的主通道的期望信号导向矢量
d.将得到的期望信号导向矢量和新构造的阻塞矩阵Ba带入到旁瓣对消最优权值的计算式中,得到辅助通道最优的权值wao和主通道权值wm。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510329662.8A CN104931935B (zh) | 2015-06-15 | 2015-06-15 | 基于旁瓣对消的稳健波束形成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510329662.8A CN104931935B (zh) | 2015-06-15 | 2015-06-15 | 基于旁瓣对消的稳健波束形成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104931935A true CN104931935A (zh) | 2015-09-23 |
CN104931935B CN104931935B (zh) | 2017-07-11 |
Family
ID=54119188
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510329662.8A Expired - Fee Related CN104931935B (zh) | 2015-06-15 | 2015-06-15 | 基于旁瓣对消的稳健波束形成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104931935B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106154234A (zh) * | 2016-07-07 | 2016-11-23 | 上海航天测控通信研究所 | 一种旁瓣对消方法及*** |
CN106850026A (zh) * | 2016-12-21 | 2017-06-13 | 上海华为技术有限公司 | 一种数据处理的方法以及相关设备 |
CN107179528A (zh) * | 2017-05-18 | 2017-09-19 | 西安电子科技大学 | 基于调幅调相辅助天线的雷达发射波束置零方法 |
CN109581366A (zh) * | 2018-12-19 | 2019-04-05 | 西安空间无线电技术研究所 | 一种基于目标导向矢量失配的离散旁瓣杂波识别方法 |
EP3945636A1 (fr) | 2020-07-28 | 2022-02-02 | Thales | Procede de modification d'un diagramme de rayonnement d'un reseau antennaire, et radar mettant en uvre un tel procede |
CN114818793A (zh) * | 2022-04-12 | 2022-07-29 | 西北工业大学 | 一种基于辅助阵元的稳健波束形成方法 |
CN115061098A (zh) * | 2022-07-28 | 2022-09-16 | 南京华成微波技术有限公司 | 雷达远距离支援干扰和随队干扰的抑制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5600326A (en) * | 1991-12-16 | 1997-02-04 | Martin Marietta Corp. | Adaptive digital beamforming architecture and algorithm for nulling mainlobe and multiple sidelobe radar jammers while preserving monopulse ratio angle estimation accuracy |
US20030184473A1 (en) * | 2002-03-27 | 2003-10-02 | Yu Kai Bor | Adaptive digital sub-array beamforming and deterministic sum and difference beamforming, with jamming cancellation and monopulse ratio preservation |
CN102830387A (zh) * | 2012-08-23 | 2012-12-19 | 北京理工大学 | 一种基于数据预处理的协方差矩阵正交化波束形成方法 |
CN103346827A (zh) * | 2013-07-26 | 2013-10-09 | 四川大学 | 稳健gsc波束形成器及其波束优化方法 |
CN104270179A (zh) * | 2014-09-12 | 2015-01-07 | 北京理工大学 | 基于协方差重构和导向矢量补偿的自适应波束形成方法 |
-
2015
- 2015-06-15 CN CN201510329662.8A patent/CN104931935B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5600326A (en) * | 1991-12-16 | 1997-02-04 | Martin Marietta Corp. | Adaptive digital beamforming architecture and algorithm for nulling mainlobe and multiple sidelobe radar jammers while preserving monopulse ratio angle estimation accuracy |
US20030184473A1 (en) * | 2002-03-27 | 2003-10-02 | Yu Kai Bor | Adaptive digital sub-array beamforming and deterministic sum and difference beamforming, with jamming cancellation and monopulse ratio preservation |
CN102830387A (zh) * | 2012-08-23 | 2012-12-19 | 北京理工大学 | 一种基于数据预处理的协方差矩阵正交化波束形成方法 |
CN103346827A (zh) * | 2013-07-26 | 2013-10-09 | 四川大学 | 稳健gsc波束形成器及其波束优化方法 |
CN104270179A (zh) * | 2014-09-12 | 2015-01-07 | 北京理工大学 | 基于协方差重构和导向矢量补偿的自适应波束形成方法 |
Non-Patent Citations (1)
Title |
---|
徐伟: "自适应波束形成算法工程应用中的关键技术研究", 《火控雷达技术》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106154234A (zh) * | 2016-07-07 | 2016-11-23 | 上海航天测控通信研究所 | 一种旁瓣对消方法及*** |
CN106850026A (zh) * | 2016-12-21 | 2017-06-13 | 上海华为技术有限公司 | 一种数据处理的方法以及相关设备 |
CN106850026B (zh) * | 2016-12-21 | 2020-04-21 | 上海华为技术有限公司 | 一种数据处理的方法以及相关设备 |
CN107179528A (zh) * | 2017-05-18 | 2017-09-19 | 西安电子科技大学 | 基于调幅调相辅助天线的雷达发射波束置零方法 |
CN107179528B (zh) * | 2017-05-18 | 2020-02-21 | 西安电子科技大学 | 基于调幅调相辅助天线的雷达发射波束置零方法 |
CN109581366A (zh) * | 2018-12-19 | 2019-04-05 | 西安空间无线电技术研究所 | 一种基于目标导向矢量失配的离散旁瓣杂波识别方法 |
CN109581366B (zh) * | 2018-12-19 | 2021-02-09 | 西安空间无线电技术研究所 | 一种基于目标导向矢量失配的离散旁瓣杂波识别方法 |
EP3945636A1 (fr) | 2020-07-28 | 2022-02-02 | Thales | Procede de modification d'un diagramme de rayonnement d'un reseau antennaire, et radar mettant en uvre un tel procede |
FR3113201A1 (fr) | 2020-07-28 | 2022-02-04 | Thales | Procede de modification d'un diagramme de rayonnement d'un reseau antennaire, et radar mettant en oeuvre un tel procede |
CN114818793A (zh) * | 2022-04-12 | 2022-07-29 | 西北工业大学 | 一种基于辅助阵元的稳健波束形成方法 |
CN115061098A (zh) * | 2022-07-28 | 2022-09-16 | 南京华成微波技术有限公司 | 雷达远距离支援干扰和随队干扰的抑制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104931935B (zh) | 2017-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104931935A (zh) | 基于旁瓣对消的稳健波束形成方法 | |
US8754810B2 (en) | Hybrid adaptive antenna array | |
CN101369014B (zh) | 应用于多输入多输出雷达的双边约束自适应波束形成方法 | |
CN104111448B (zh) | Mimo雷达收发阵列误差的联合校正方法 | |
CN105137399B (zh) | 基于斜投影滤波的雷达自适应波束形成方法 | |
CN106324625B (zh) | 一种基于二范数多目标优化的卫星导航***自适应抗干扰方法 | |
CN107340499A (zh) | 基于协方差矩阵重建的稳健低副瓣波束形成方法 | |
CN108631851B (zh) | 基于均匀直线阵列零陷加深的自适应波束形成方法 | |
CN105655727B (zh) | 一种基于凸优化的增益损失恒定波束的形成方法及装置 | |
CN105403871B (zh) | 一种双基地mimo雷达阵列目标角度估计和互耦误差校准方法 | |
CN107315162A (zh) | 基于内插变换和波束形成的远场相干信号doa估计方法 | |
CN105354171B (zh) | 一种改进特征矢量的投影子空间估计自适应波束合成方法 | |
CN106443594A (zh) | 一种基于稀疏约束的雷达天线阵列稳健波束形成方法 | |
CN108880586B (zh) | 一种宽带弱信号增强方法与装置 | |
CN103037494B (zh) | 多波束天线功率控制方法 | |
CN103475394A (zh) | 一种星载相控阵天线自适应抑制干扰的方法 | |
CN107167804A (zh) | 一种稳健低副瓣自适应波束形成方法 | |
CN104931937B (zh) | 基于协方差矩阵归一化的子阵级正交投影波束形成方法 | |
CN109283496B (zh) | 一种抗运动干扰和抗导向失配的稳健波束形成方法 | |
CN101825715B (zh) | 一种卫星导航接收设备波束形成的方法 | |
CN103248412B (zh) | 一种卫星导航多波束形成的方法 | |
CN111817765B (zh) | 一种基于频率约束的广义旁瓣对消宽带波束形成方法 | |
Kikuchi et al. | Autocalibration algorithm for robust Capon beamforming | |
Suleesathira | Robust null broadening beamforming based on adaptive diagonal loading for look direction mismatch | |
CN114879173A (zh) | 基于mimo-ofdm信号的感知性能优化方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170711 Termination date: 20200615 |