CN101483280A - 一种稳健波束合成器权值求解方法 - Google Patents

一种稳健波束合成器权值求解方法 Download PDF

Info

Publication number
CN101483280A
CN101483280A CNA2009101032263A CN200910103226A CN101483280A CN 101483280 A CN101483280 A CN 101483280A CN A2009101032263 A CNA2009101032263 A CN A2009101032263A CN 200910103226 A CN200910103226 A CN 200910103226A CN 101483280 A CN101483280 A CN 101483280A
Authority
CN
China
Prior art keywords
array
beam synthesizer
arrival
vector
weight
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
Application number
CNA2009101032263A
Other languages
English (en)
Other versions
CN101483280B (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.)
Chongqing University
Original Assignee
Chongqing University
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 Chongqing University filed Critical Chongqing University
Priority to CN2009101032263A priority Critical patent/CN101483280B/zh
Publication of CN101483280A publication Critical patent/CN101483280A/zh
Application granted granted Critical
Publication of CN101483280B publication Critical patent/CN101483280B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

一种稳健波束合成器权值求解方法,其步骤为:第一步,通过对阵列接收信号x(t)进行AD采样,估计出协方差矩阵R;第二步,通过协方差矩阵进行波达方向估计,得到期望信号和M个干扰信号的波达方向;第三步,由波达方向得到M+1个子权矢量w0,…,wM;第四步,阵列总的权矢量按照公式
Figure 200910103226.3_AB_0
得到。这种稳健波束合成器的优势在于不需要通过增加约束条件对阵列误差作复杂校正,从而具有较少运算量的同时,又实现波束合成器的稳健性。

Description

一种稳健波束合成器权值求解方法
技术领域
本发明涉及波束合成领域,具体涉及一种稳健波束合成器权值求解方法。
技术背景
波束合成器在通信、雷达中广泛应用,通过一定的最优化方法,对阵列天线接收信号x(t)进行信号处理,计算出阵列的权矢量w(t),从而使阵列方向图主瓣对准期望信号,而零陷对准干扰信号,达到抑制干扰,提高***输出信干噪比的效果。但是,由于阵列位置误差,通道不一致误差,阵元耦合误差的存在,对信号的波达方向估计是存在误差的。误差产生的结果,会使得期望信号不在方向图主瓣范围内,而干扰也不在方向图零陷范围内,这样就使波束合成器失去抑制干扰的功能,不但不能提高输出信干噪比,甚至有可能还减少输出信干噪比。
为了克服误差影响,需要寻找稳健的波束合成器。已有的稳健波束合成器包括投影方法,对角负载方法,贝叶斯方法,凸集优化方法,其目的都是试图把主瓣和零陷宽度加宽,使期望信号不至于超出主瓣范围,干扰信号不至于超出零陷范围。这些方法可以获得很好误差校正效果,但是,由于增加大量约束条件,使得权矢量计算增加了非常大的运算量。
发明内容
本发明的目的,就是在基本波束合成器实现方法基础上,不增加新的约束条件和运算量,仅对零陷进行加宽,从而用较少运算量克服由于阵列各种误差引起的波达方向估计误差,使干扰信号始终位于零陷范围内,实现稳健的波束合成器。
一种稳健波束合成器权值求解方法,其实现步骤为:第一步,对阵列(I)接收信号x(t)进行K次采样,得到K个快拍x(1),...,x(K),从而在权值求解单元(II)中构成数据矩阵X=[x(1)...x(K)],按照公式 R = 1 K XX H 估计得到阵列协方差矩阵R。第二步,在权值求解单元(II)中,通过协方差矩阵R估计出空间期望信号s0(t)的波达方向θ0,以及M个干扰信号s1(t),...,sM(t)的波达方向θ1,...,θM。第三步,在权值求解单元(II)中,由期望信号波达方向θ0得到方向矢量v(θ0),从而得到一个子权矢量w0=v(θ0);由M个干扰波达方向θ1,...,θM得到M个方向矢量v(θ1),...,v(θM),从而得到M个子权矢量w1=v(θ1),...,wM=v(θM)。第四步,阵列总的权矢量按照公式 w ( t ) = w 0 - Σ m = 1 M w m 得到,并把该权矢量传递给加权求和单元(III);加权求和单元(III)按照公式y(t)=wH(t)×x(t)计算阵列加权输出。
该发明的有益效果是没有采用复杂的约束条件,但零陷宽度得到增加,从而使权值求解非常简单,并实现的波束合成的稳健性效果。
附图说明
图1是稳健波束合成器权值求解框图;
图2是稳健波束合成器结构框图;
图3是均匀线阵阵列接收信号示意图;
图4是稳健波束合成器方向图;
图5是加权求和单元结构图。
具体实施方式
方法实现步骤如图1所示。***结构如图2所示,波束合成器包括了天线阵列单元(I),权值求解单元(II)和加权求和单元(III)。
如图3所示,假设阵列天线为均匀线阵,阵元间距为载波波长一半,而空间具有一个期望信号s0(t),波达方向θ0,M个干扰信号分别为s1(t),...,sM(t),波达方向分别为θ1,...,θM,这些信号是远场窄带不相干信号。阵列接收信号矢量可以表示为 x ( t ) = Σ m = 0 M s m ( t ) v ( θ m ) , 其中v(θm)为波达方向为θm的信号对应的方向矢量。
在权值求解单元(II)中,对接收信号矢量进行AD变换,得到K个连续的采样信号,每个采样信号称为一个快拍,即x(k),k=1,...,K。把K个快拍数据矢量构成一个数据矩阵,即X=[x(1)...x(K)],按照公式 R = 1 K XX H 估计阵列接收信号的协方差矩阵R。对协方差矩阵R进行空间谱估计,计算出相应的期望信号和干扰信号的波达方向。这个空间谱估计可以采用MUSIC算法,步骤为:
第一步,用序列假设(SH)算法,或者赤池信息论准则(AIC)方法,或者最小描述长度(MDL)方法,估计出期望和干扰信号个数M+1。
第二步,对协方差矩阵进行特征分解,取其M+1个较大特征值对应特征矢量e0,...,eM,并把这M+1个特征矢量构成一个矩阵F=[e0...eM]。
第三步,定义变量θ,θ取值范围为00~1800,计算空间谱函数 P ( θ ) = v H ( θ ) v ( θ ) v H ( θ ) FF H v ( θ ) , 该函数必然有M+1个谱峰,即局部最大值,而这M+1个谱峰对应的θ,就是θ0,...,θM
利用MUSIC法估计出θ0,...,θM后,则可以计算阵列权矢量w(t)。基本的波束合成器权矢量按照公式 w ( t ) = R - 1 v ( θ 0 ) v H ( θ 0 ) R - 1 v ( θ 0 ) 计算,但由于阵列各种误差的存在,MUSIC法估计出的波达方向是有误差的,误差可能使得干扰超出方向图零陷范围,所以,基本波束合成器是不具备稳健性的。而本发明的权值求解则在估计出波达方向后,进行如下所示方法计算:
第一步,用的M+1个波达方向,分别计算M+1个方向矢量v(θ0),...,v(θM)。
第二步,令w0=v(θ0),...,wM=v(θM)。
第三步,阵列最终权矢量 w ( t ) = w 0 - Σ m = 1 M w m .
图4所示,为8阵元阵列,期望信号位于900,两个干扰位于600和1200时本发明稳健波束合成器方向图,图中可以看出零陷宽度较大,可以校正波达方向估计误差的影响,保证干扰不会超出零陷范围。
权值w(t)获得以后,传递给加权求和单元(III),对阵列接收信号矢量x(t)进行加权,并求和,即完成计算y(t)=wH(t)x(t)。阵列的输出为y(t)。加权求和单元(III)如图4所示。
本发明相对于传统波束合成器,具有克服误差的稳健性,同时相对于投影方法,对角负载方法,贝叶斯方法,凸集优化方法的稳健波束合成器,没有约束条件,运算量小的多。

Claims (3)

1、一种稳健波束合成器权值求解方法,其实现步骤为:
第一步,对阵列(I)接收信号x(t)进行K次采样,得到K个快拍x(1),...,x(K),从而在权值求解单元(II)中构成数据矩阵X=[x(1)...x(K)],按照公式 R = 1 K XX H 估计得到阵列协方差矩阵R;
第二步,在权值求解单元(II)中,通过协方差矩阵R估计出空间期望信号s0(t)的波达方向θ0,以及M个干扰信号s1(t),...,sM(t)的波达方向θ1,...,θM
第三步,在权值求解单元(II)中,由期望信号波达方向θ0得到方向矢量v(θ0),从而得到一个子权矢量w0=v(θ0),由M个干扰波达方向θ1,...,θM得到M个方向矢量v(θ1),...,v(θM),从而得到M个子权矢量w1=v(θ1),...,wM=v(θM);
第四步,阵列总的权矢量按照公式 w ( t ) = w 0 - Σ m = 1 M w m 得到,并把该权矢量传递给加权求和单元(III)。
2、根据权利要求1所述稳健波束合成器权值求解方法,其特征在于:天线阵(I)的N个阵元间距为接收信号波长的一半,阵列可以是均匀线阵,也可以是均匀圆阵。
3、根据权利要求1所述稳健波束合成器权值求解方法,其特征在于:加权求和单元(III)按照公式y(t)=wH(t)×x(t)计算阵列加权输出y(t)。
CN2009101032263A 2009-02-23 2009-02-23 一种稳健波束合成器权值求解方法 Expired - Fee Related CN101483280B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009101032263A CN101483280B (zh) 2009-02-23 2009-02-23 一种稳健波束合成器权值求解方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009101032263A CN101483280B (zh) 2009-02-23 2009-02-23 一种稳健波束合成器权值求解方法

Publications (2)

Publication Number Publication Date
CN101483280A true CN101483280A (zh) 2009-07-15
CN101483280B CN101483280B (zh) 2010-09-22

Family

ID=40880277

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009101032263A Expired - Fee Related CN101483280B (zh) 2009-02-23 2009-02-23 一种稳健波束合成器权值求解方法

Country Status (1)

Country Link
CN (1) CN101483280B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101908918A (zh) * 2010-07-26 2010-12-08 重庆大学 无线通信接收机中的一种波束合成方法
CN103117788A (zh) * 2013-02-20 2013-05-22 西安邮电大学 智能天线稳健干扰抑制方法
CN103245934A (zh) * 2013-06-03 2013-08-14 重庆大学 一种大动态期望信号下的波束合成方法
CN103634038A (zh) * 2013-12-09 2014-03-12 南京邮电大学 基于多天线的联合doa估计和波束形成的多径信号接收方法
CN107241131A (zh) * 2017-06-02 2017-10-10 电子科技大学 一种利用信号非圆特性的波束形成方法
CN111241470A (zh) * 2020-01-19 2020-06-05 河北科技大学 基于自适应零陷展宽算法的波束合成方法及装置
CN114563760A (zh) * 2022-02-07 2022-05-31 哈尔滨工程大学 一种基于sca阵型的二阶超波束形成方法、设备及介质

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100512047C (zh) * 2006-10-19 2009-07-08 北京交通大学 阵列天线mc-cdma***用户信号波达方向估计方法

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101908918B (zh) * 2010-07-26 2015-04-01 重庆大学 无线通信接收机中的一种波束合成方法
CN101908918A (zh) * 2010-07-26 2010-12-08 重庆大学 无线通信接收机中的一种波束合成方法
CN103117788A (zh) * 2013-02-20 2013-05-22 西安邮电大学 智能天线稳健干扰抑制方法
CN103245934A (zh) * 2013-06-03 2013-08-14 重庆大学 一种大动态期望信号下的波束合成方法
CN103245934B (zh) * 2013-06-03 2014-10-29 重庆大学 一种大动态期望信号下的波束合成方法
CN103634038B (zh) * 2013-12-09 2016-09-28 南京邮电大学 基于多天线的联合doa估计和波束形成的多径信号接收方法
CN103634038A (zh) * 2013-12-09 2014-03-12 南京邮电大学 基于多天线的联合doa估计和波束形成的多径信号接收方法
CN107241131A (zh) * 2017-06-02 2017-10-10 电子科技大学 一种利用信号非圆特性的波束形成方法
CN107241131B (zh) * 2017-06-02 2020-11-13 电子科技大学 一种利用信号非圆特性的波束形成方法
CN111241470A (zh) * 2020-01-19 2020-06-05 河北科技大学 基于自适应零陷展宽算法的波束合成方法及装置
CN111241470B (zh) * 2020-01-19 2023-08-18 河北科技大学 基于自适应零陷展宽算法的波束合成方法及装置
CN114563760A (zh) * 2022-02-07 2022-05-31 哈尔滨工程大学 一种基于sca阵型的二阶超波束形成方法、设备及介质
CN114563760B (zh) * 2022-02-07 2023-02-07 哈尔滨工程大学 一种基于sca阵型的二阶超波束形成方法、设备及介质

Also Published As

Publication number Publication date
CN101483280B (zh) 2010-09-22

Similar Documents

Publication Publication Date Title
CN101483280B (zh) 一种稳健波束合成器权值求解方法
CN107315160B (zh) 基于内插虚拟阵列信号原子范数最小化的互质阵列波达方向估计方法
CN105891771B (zh) 一种提高估计精度的基于连续分布的角度估计方法与设备
CN108896954B (zh) 互质阵中一种基于联合实值子空间的波达角估计方法
CN104698433A (zh) 基于单快拍数据的相干信号doa估计方法
CN103245956B (zh) 一种基于稳健波束形成算法的gps抗多径方法
CN108462521A (zh) 自适应阵列天线的抗干扰实现方法
CN106980106A (zh) 阵元互耦下的稀疏doa估计方法
CN106054123A (zh) 一种稀疏l阵及其二维doa估计方法
CN104977558A (zh) 一种基于贝叶斯压缩感知的分布源中心波达方向估计方法
CN104991236B (zh) 一种单基地mimo雷达非圆信号相干源波达方向估计方法
CN111707985A (zh) 基于协方差矩阵重构的off-grid DOA估计方法
CN107092007A (zh) 一种虚拟二阶阵列扩展的波达方向估计方法
CN107104720B (zh) 基于协方差矩阵虚拟域离散化重建的互质阵列自适应波束成形方法
CN110515033B (zh) 一种基于Toeplitz矩阵恢复的欠通道测向***及方法
CN107315162A (zh) 基于内插变换和波束形成的远场相干信号doa估计方法
CN101644760B (zh) 一种适用于高分辨阵列的快速鲁棒的信源个数检测方法
CN110824414A (zh) 一种到达角估计的装置及方法
CN111983556B (zh) 一种到达角估计的装置及方法
CN108594166A (zh) 一种二维波达方向估计方法及装置
CN104330766A (zh) 一种稳健的波达方向估计方法
CN101344582B (zh) 半盲的最小方差无失真响应波束形成方法
CN102394686B (zh) 高精度阵列天线接收***角度估计的装置及其方法
CN116224219A (zh) 一种阵列误差自校正原子范数最小化doa估计方法
CN103701515A (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
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20100922

Termination date: 20120223