CN115270869A - 一种基于二维空时广义旁瓣相消稳健波束形成算法 - Google Patents

一种基于二维空时广义旁瓣相消稳健波束形成算法 Download PDF

Info

Publication number
CN115270869A
CN115270869A CN202210882587.8A CN202210882587A CN115270869A CN 115270869 A CN115270869 A CN 115270869A CN 202210882587 A CN202210882587 A CN 202210882587A CN 115270869 A CN115270869 A CN 115270869A
Authority
CN
China
Prior art keywords
time
signal
space
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
Application number
CN202210882587.8A
Other languages
English (en)
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.)
Ningbo Qise Jia Metal Products Co ltd
Original Assignee
Ningbo Qise Jia Metal Products Co ltd
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 Ningbo Qise Jia Metal Products Co ltd filed Critical Ningbo Qise Jia Metal Products Co ltd
Priority to CN202210882587.8A priority Critical patent/CN115270869A/zh
Publication of CN115270869A publication Critical patent/CN115270869A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/17Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Radio Transmission System (AREA)

Abstract

本发明涉及一种基于二维空时广义旁瓣相消稳健波束形成算法,主要使用二维数组在广义旁瓣相消器(Generalized Sidelobe Canceller,GSC)架构下,解决指向误差及信号源数目过度估测问题。首先将现有技术之算法使用在GSC作指向误差校正,但此方法的问题依旧存在。为使校正能力更为稳健本发明提出新方法,基于特征分析处理之稳健空时广义旁瓣相消器(RobustSpace‑Time GSC,RST‑GSC),此方法利用CDMA技术结合二维数组天线,在有指向误差以及信号源过度估测发生时,可自适应让导向矢量有实时校正的能力。

Description

一种基于二维空时广义旁瓣相消稳健波束形成算法
技术领域
本发明涉及智能型天线指向误差以及信号源过度估测发生时,可自适应让导向矢量有实时校正的能力的技术领域,具体为一种基于二维空时广义旁瓣相消稳健波束形成算法。
背景技术
智能型天线技术的应用,对信号强度与系统容量的提升有很大的帮助,尤其在强调以CDMA技术核心的系统之中,对于消除多重接取所产生的干扰信号有很大的帮助,因此在未来6G/7G通讯的系统的规格中,建议使用智能型天线的技术来提升系统的容量与传输效率。
智能型天线是一种以多传感器所组成的天线阵列,利用天线阵列的波束形成技术,产生多个独立的波束,來追踪每一个用户,而应用智能天线系统可将天线信号导向期望方向,产生空间定向波束,使天线主波束对准用户信号到达方向(Direction of Arrival,DOA),旁波瓣或置零点对准干扰信号到达方向,可以达到充分有效利用行动用户信号并删除或抑制干扰信号的目的。因此,要实现智能型天线系统使用到的阵列天线理論,主要可分成三部分:信号源数目估测、到达角度估测与波束构成。
在文献[8]中,有两个信号源数目估测方法如Akaike信息准则(AkaikeInformation Criterion,AIC)和最小化描述长度(Minimum Description Length, MDL)方法已经被提出,由于信号源数目少估会比过度估测付出更高的代价,因此估测信号个数时常会有信号源过度估测之情况产生。而系统在实施到达角度估测时,会因为许多环境的变动因素造成角度估测不准确。当到达角度估测并不准确时,期望信号入射角度与接收机预估之入射角度并不吻合,将产生指向误差的问题,此时期望信号无法落在主波束最大响应位置,造成期望信号消除。当上述两问题同时存在且波束形成器使用特征分析算法,期望信号成分会被归类到干扰子空间,使系统效能降低。
智能型天线在时间-空间系统下线性数组信号模型[9],并将分别介绍未经过解扩频及经过解扩频之后的信号模型,为了简化后面章节对天线数组的分析,首先作了以下假设:
1.天线单元间的距离足够小,因此可以假设不同天线单元间所接收到的信号振幅变化可以忽略,只考虑路径引起的相位差。
2.任何天线单元间不会发生耦合效应且入射信号的带宽远小于载波频率,可以将入射信号视为窄频带信号。
3.假设信号源与天线数组间的距离足够远,如此便可以将此信号视为远场信号源,并且将入射波视为一平面波。
考虑一个含有K个信号源的线性系统,使用正规化扩频码c1,c2,…,cK,来传送二进制相移键控,二进制相移键控信号周期为Tb,因此接收机所接收到的用户信号可以表示为
Figure RE-GDA0003851072190000021
其中bk(z)∈±1代表第k个使用者的第z个数据位被扩频码扩频,εk为所接收到第k个使用者的功率,假设扩频码长度为L,此时第k个使用者经标准化的扩频码为
Figure RE-GDA0003851072190000031
n(z)为信号通过信道所产生的噪声,在此假设为白色高斯噪声。
假设一UPA数组天线含有M×N个全向性天线阵列,各天线阵列间距为d=0.5λ并且座落在X-Y平面上的
Figure RE-GDA0003851072190000032
此时 m=1,2,…,M且n=1,2,…,N,假设信号源的入射高低角及水平角分别为θ、φ,假设天线的振幅为单位增益,且相位增益为exp{jπ(m-1)μ+(n-1)υ}、μ=sin(θ)cos(φ)、υ=sin(θ)sin(φ)。如图1的数组天线所示,则第(m,n)个天线阵列所接收到的信号可以表示成:
Figure RE-GDA0003851072190000033
将平面型数组输出之信号表示为:
Figure RE-GDA0003851072190000034
上式中{·}T代表转置、ac(uk)=[1,exp{jπμk},…,exp{jπ(M-1)μk}]T, ark)=[1,exp{jπυk},…,exp{jπ(N-1)υk}]T,N为所有天线所收到之噪声所组成之矩阵,为了使自由度提高,将二维均匀平面型天线数组用一维表示,因此输出信号为:
Figure RE-GDA0003851072190000035
此时
Figure RE-GDA0003851072190000036
为第k个使用者的导向矢量,
Figure RE-GDA0003851072190000037
定义为Kronecker乘法[附录A],
Figure RE-GDA0003851072190000038
为二维数组展开后与时间特征合并之噪声矩阵,假设第k个使用者就是期望使用者k=d,入射角度为(μdd),此时有K-1 个干扰,并且将期望信号的导向矢量表示为ad。此时将
Figure RE-GDA0003851072190000041
以二阶统计量自相关矩阵表示:
Figure RE-GDA0003851072190000042
理论上,在稳态环境下天线数组输出自相关矩阵应由整体平均得到,但于实际操作上,可以由时间平均来获得估测的自相关矩阵[12]:
Figure RE-GDA0003851072190000043
其中Z为总共观察的位数。本发明将使用空间-时间之特征空间,因此将接收信号之扩频码与导向矢量再表示成时空数组形式[6]:
Figure RE-GDA0003851072190000044
其中
Figure RE-GDA0003851072190000045
Figure RE-GDA0003851072190000046
为将空间-时间因素合并后所之噪声矢量,而时空信号所形成的自相关矩阵表示为:
Figure RE-GDA0003851072190000047
总之,智能型天线是一种以多传感器所组成的天线阵列,利用天线阵列的波束形成技术,产生多个独立的波束,來追踪每一个用户,而应用智能天线系统可将天线信号导向期望方向,产生空间定向波束,使天线主波束对准用户信号到达方向(Direction ofArrival,DOA),旁瓣或零陷对准干扰信号到达方向,可以达到充分有效利用行动用户信号并删除或抑制干扰信号的目的。
但系统在实施到达角度估测时,会因为许多环境的变动因素造成角度估测不准确。当到达角度估测并不准确时,期望信号入射角度与接收机预估之入射角度并不吻合,将产生指向误差的问题,此时期望信号无法落在主波束最大响应位置,造成期望信号消除。
为了解决指向误差之问题,文献[1]中提出了线性限制最小变异波束构成器(Linear Constrained Minimum Variance Beamformer,LCMVB),此方法以附加约束法来减轻数组的敏感性,其中微分限制是被用于减轻由指向误差所引起的效能衰退。经过微分限制主波束将会变宽,当指向误差不大时期望信号被消除并不严重,而零陷技术用于抑制非固定的干扰。无论如何,此方法有两个缺点,当干扰入射角度在主波束范围内将无法对此干扰做有效抑制,且利用此一降低变量之自由度以求得最优权重矢量的方法,所获得的结果会使数组效能在干扰和噪声抑制上产生衰退。
后来陆续提出其他广义特征空间波束形成器(Generalized Eigenspace-BasedBeamformer,GEIB)文献[2]-[4],这个侦测器对较小的指向误差具有稳健性,侦测方式是将接收到的信号先作自相关矩阵并利用特征空间的观念将此信号自相关矩阵分解成信号子空间和噪声子空间并结合 LCMVB,但使用此方法有一个严重的缺点就是必须正确知道信号源的数目。否则,当期望信号功率小于噪声功率时,此侦测方法将造成干扰源数目过度估测的现象而使整个系统的性能降低。
为了解决此问题,文献[11]中提出了稳健的广义特征空间波束形成器,此形成器利用传送信号及接受信号的交互相关矩阵,经过处理会和正确导向矢量之间有一倍数关系,利用此关系式估出正确的信号源入射角度,并重建导向矢量以找出最优权重值。但此方法亦无法对抗靠近主波瓣宽度边缘附近的较大指向误差问题,当误差角度过大效能将会崩溃。
参见图1-图7所示,广义旁瓣相消器(Generalized Sidelobe Canceller, GSC)。GSC为基本架构,如图1所示。以下将介绍GSC基本原理、GSC 在二维UPA平面上之应用并仿真此方法在完美环境下之效能[6][7]。
基于特征分析之广义旁瓣相消器原理
首先GSC架构的权重矢量为:
wGSC=wq-Bwa (2.8)
其中
Figure RE-GDA0003851072190000061
导向矢量a1之维度为MN×1,B为阻塞矩阵也就是a1的零矩阵满足BHa1=0,维度为MN×(MN-1)并满足BHB=I。此时最优权重解wa可经由计算所有数组输出最小化所得,如下式:
Figure RE-GDA0003851072190000062
由式(2.30)求得最优适应性权重解wa
Figure RE-GDA0003851072190000063
上式中
Figure RE-GDA0003851072190000064
设U(z)=BHX(z),RB为由U(z)所形成之二阶统计量自相关矩阵表示如下:
Figure RE-GDA0003851072190000065
欲使用基于特征分析广义旁瓣相消器必须得到RB之特征矢量,因此对 RB作特征值分解如下:
Figure RE-GDA0003851072190000066
此时
Figure RE-GDA0003851072190000067
为所有的特征值,而ei对应到第i个特征值λi。由干扰特征值所形成的对角矩阵为ΛI=diag{λ12,…,λK-1},维度为(K-1)×(K-1);由噪声特征值所形成之对角矩阵为Λn=diag{λKK+1,…,λMN-1},维度为(MN-K)×(MN-K)。最后对应每个干扰特征值可以形成干扰子空间EI=[e1,e2,…,eK-1];对应每个噪声特征值可以形成噪声子空间En=[eK,eK+1,…,eMN-1]。从第(2.11)式可以知道En和阻塞矩阵B为正交
Figure RE-GDA0003851072190000073
在文献[10]中也提到此时p只包含干扰成分,所以此时
Figure RE-GDA0003851072190000074
由式(2.11)、(2.12)经过传统波束形成器权重解可以改写为:
Figure RE-GDA0003851072190000071
由式(2.13)可以得知适应性权重最优解wao矢量落在干扰子空间,因此无论期望信号成份存在与否,式(2.13)为GSC架构下适应性波束形成器之最优解。最后假设在有限的观察时间内,由所有数组获得最优权重解为:
Figure RE-GDA0003851072190000072
计算机仿真与分析:仿真当信号源入射角度及信号源个数完美估测时,使用特征GSC作为波束形成器之架构并使用wGSC作为权重矢量。假设一 6×6的均匀平面型数组,天线之间的距离为0.5λ,此时有4个信号源,分别由角度(μ,υ)=(0,0)、(-0.55,0)、(0.53,0)、(0.57,0)同时入射,而第一个信号源为期望信号,信号强度为SNR=10d,其他信号源为干扰,干扰信号为等强度INR=30dB。下图2、3、4为模拟图。首先,图2为二维均匀平面数组的波束方向图。图3为扩频信号的功率频谱的υ=0,变化μ从-1至+1来观察一维波束方向图,由此图可以看出在完美环境下,GSC可以将期望信号稳健住且将干扰抑制,即使在干扰源非常接近的情况下效果依旧很好。图4 令μ=0,变化υ从-1至+1的一维波束方向图。
信号源数目过度估测对系统效能之影响:当波束形成器使用特征空间算法时,一旦有信号源数目过度估测之情形发生,在撷取信号子空间时会将部分噪声子空间纳入范围内,使信号无法完整地投影在信号子空间,造成效能降低。但一般来说,只要在没有指向误差之环境下,信号源数目过度估测对性能不会有太大的影响。
指向误差对系统效能之影响:指向误差表示系统接收端所接收信号入射角度与预估的信号入射角度产生误差,导致天线数组的波束指向错误方向,而真正的信号被当成干扰消除掉,使系统效能崩溃。在某些算法下,当期望信号越强甚至会被消除得更彻底。此问题的发生,也有可能是因为天线受环境的影响有扰动的情形,使得期望信号的入射角度和预估的角度产生误差,以下为计算机仿真。
计算机仿真与分析:仿真当环境有指向误差及信号源数目过度估测时,使用GSC作为波束形成器之架构并使用wGSC作为权重矢量。假设一6×6的均匀平面型数组,天线之间的距离为0.5λ,此时有4个信号源,分别由角度 (μ,υ)=(0,0)、(-0.55,0)、(0.53,0)、(0.57,0)同时入射,而第一个信号源为期望信号,信号强度SNR=10dB;其他信号源为干扰,干扰信号为等强度 INR=30dB。在有指向误差(μ,υ)=(0.2,0)及信号源数目过度估测Je=J+1 的情形下,下图5、6、7为模拟图。首先,图5为二维均匀平面数组的波束方向图。图6为将图5的υ固定在0,变化μ从-1至+1来观察一维波束方向图,由此图可以看出GSC不会对指向误差做校正,因此会将期望信号当作干扰消除,造成系统效能降低。图7令μ=0,变化υ从-1至+1的一维波束方向图。此方法将在下面与其他稳健性方法作比较。
上述介绍了数组的基本观念及一些基本的波束形成器,并在系统模型结合时间处理及空间处理成为时空信号。最后仿真GSC在有信号源数目过度估测及指向误差环境下的系统,可以得知GSC对于干扰依然有很好的消除效果,但由于无法处理指向误差,造成系统效能崩溃。
上述GSC会因波束形成算法遭遇指向误差以及信号源数目过度估测时造成系统效能崩溃。为了解决期望信号强度较强以及指向误差较大时崩溃之问题。因此本发明提出稳健空时广义旁瓣相消器(Robust Space-Time GSC, RST-GSC),利用二维空时矩阵估出期望导向矢量。由仿真可得知此方法无论在信号功率较强或指向误差很大时,都有良好之性能。
参见图8-18所示,现有技术中的稳健广义旁瓣相消器:稳健广义旁瓣相消器使用一基于特征空间之导向矢量估测器结合GSC,架构图如下图8所示。此算法用在ULA时可估出信号源入射角度;但当此算法使用在UPA上时,由于无法同时估测俯仰角及水平角,只用来估测导向矢量。以下为RGSC 使用基于特征空间之导向矢量估测器介绍。
基于特征空间之导向矢量估测器:假设ad为正确导向矢量,a1为估测的导向矢量,当有指向误差存在时ad≠a1;入射角度(μdd)≠(μ11)。此时经由GSC所得到的权重并不完全落在信号子空间,会有期望信号消除的情形发生,使性能崩溃。为解决此问题可因此经由下式:
Figure RE-GDA0003851072190000091
对上式(3.1)使用拉格朗日乘法(Lagrange Multipliers)求出v之最优解:
Figure RE-GDA0003851072190000101
其中I为维度为MN×MN的单位矩阵,En是由(2.26)式之
Figure RE-GDA0003851072190000102
作特征值分解得到噪声子空间,如下:
Figure RE-GDA0003851072190000103
此时
Figure RE-GDA0003851072190000104
为所有的特征值,而ei对应到第i个特征值λi。由干扰特征值所形成的对角矩阵为ΛI=diag{λ23,…,λK},维度为(K-1)×(K-1);由噪声特征值所形成之对角矩阵为Λn=diag{λK+1K+2,…,λMN},维度为(MN-K)(MN-K)。最后对应每个干扰特征值可以形成干扰子空间EI=[e2,e3,…,eK];对应每个噪 0声特征值可以行程噪声子空间En=[eK+1,eK+2,…,eMN]。
由式(3.2)可知,真正的导向矢量第一个元素为1且每ㄧ个元素的振幅也都为1,因此对此矢量作标准化得到一矢量f=(v1)-1v,v1代表v的第一个元素。令fi为f的第i个元素,因此可以从此关系式估出正确的导向矢量:
Figure RE-GDA0003851072190000105
Figure RE-GDA0003851072190000106
代表
Figure RE-GDA0003851072190000107
第i个元素。最后将
Figure RE-GDA0003851072190000108
代入GSC架构中的阻塞矩阵
Figure RE-GDA0003851072190000109
以及权重矢量,得到新的权重矩阵:
Figure RE-GDA00038510721900001010
其中
Figure RE-GDA00038510721900001011
Figure RE-GDA00038510721900001012
为所欲求得的适应性权重
Figure RE-GDA00038510721900001013
此时的
Figure RE-GDA00038510721900001014
经由上述之替换得到之最终权重矢量解为
Figure RE-GDA0003851072190000111
式子中
Figure RE-GDA0003851072190000112
这时接收机所接受之信号通过阻塞矩阵能够被阻隔较完全,但由于 RGSC的估测能力并不好,使效能只能获得些许提升。
计算机仿真与分析:本节之模拟使用
Figure RE-GDA0003851072190000113
作为波束形成器之权重,假设一8×8的均匀平面型数组,天线之间的距离为0.5λ,此时有4个信号源分别由角度(μ,υ)=(0,0)、(-0.55,0)、(0.53,0)、(0.57,0)同时入射,而第一个信号源为期望信号,信号强度为SNR=0dB;其他信号源为干扰,干扰信号皆为等强度INR=30dB。图9为在有指向误差(μ,υ)=(0.2,0)但信号源数目正确估测的情形下之二维波束方向图。图10令图9的υ=0,变化μ从-1至 +1来观察一维波束方向图,由图可看出此方法对于强干扰在信号源正确估测时,有不错的校正能力。图13令μ=0,变化υ从-1至+1的一维波束方向图;而图12、13、14是在有指向误差(μ,υ)=(0.2,0)且信号源数目过度估测 Je=J+1的情形下,使用RGSC作为波束形成器算法之波束方向图。由图 13能发现即使信号源过度估测,指向误差仍然会被校正,但由于噪声子空间不再与信号子空间正交,造成旁波隆起且噪声放大。
这个部分为在有指向误差(μ,υ)=(0.2,0)且信号源数目过度估测 Je=J+1的情形下,利用图15、16、17、18观察RGSC与GSC的性能比较。图15将GSC及RGSC之波束方向图放在一起比较,可看出传统GSC不会对指向误差做校正,因此会将期望信号当作干扰消除;而RGSC会对指向误差做校正,但结果并不好且在信号源过度估测时会有噪声放大的情形产生。图16变化输入数据之位数由100到103来观察输出的SINR,由此图可看出两个方法的收敛速度差不多,大约在102位就能够收敛。图17变化期望信号之 SNR值由1到20来观察输出的SINR,由此图可看出由于RGSC估测不够准确,当SNR提高期望信号反而消除越厉害,因此两个方法在SNR提高的情形下,都会使系统性能崩溃。图18变化指向误差由μ从-0.2至0.2,由此图可看出在SNR=0的情况下,RGSC对抗指向误差能力较GSC好。
附录A:Kronecker乘法简介
定义:假设A是一个M×N的矩阵,B为一个P×Q矩阵此时经过 Kronecker乘法会产生一个MP×NQ的矩阵如下。
Figure RE-GDA0003851072190000121
更详细的表示方法为
Figure RE-GDA0003851072190000122
使用Kronecker乘法同时也会满足下列特性:
1.
Figure RE-GDA0003851072190000131
2.
Figure RE-GDA0003851072190000132
3.
Figure RE-GDA0003851072190000133
4.
Figure RE-GDA0003851072190000134
5.
Figure RE-GDA0003851072190000135
6.
Figure RE-GDA0003851072190000136
上式中k为常数,以上为Kronecker乘法的ㄧ些基本概念。
参考现有技术文献:
[1]B.D.Van Veen,“Eigenstructure based partially adaptive arraydesign,” IEEE Trans.Antennas Propagation,vol.36,no.1,pp.357-362,March 1988.
[2]A.M.Haimovich and Y.Bar-Ness,“An eigenanalysis interferencecanceler,”IEEE Trans.Signal Processing,vol.39,no.1,pp.76-84,Jan. 1991.
[3]D.D.Feldman and L.J.Griffiths,“A constraint projection approachfor robust adaptive beamforming,”Proceedings of IEEE International Conferenceon Acoustics,Speech,and Signal Processing,Toronto,Canada, pp.1381-1384,June1991.
[4]L.Chang and C.C.Yeh,“Performance of DMI and eigenspace-basedbeamformers,”IEEE Trans.Antennas Propagation,vol.40,no.11,pp. 1336-1347,Nov.1992.
[5]S.Haykin,Adaptive Filter Theory,4th edition,Prentice Hall,NewJersey, 1996.
[6]L.J.Griffiths and C.W.Jim,“An alternative approach to linearlyconstrained adaptive beamforming,”IEEE Trans.Acoustics,Speech,and SignalProcessing,vol.30,no.1,pp.27-34,Jan.1982.
[7]K.M.Buckley and L.J.Griffiths,“An adaptive generalized sidelobecanceller with derivative constraints,”IEEE Trans.Antennas Propagation,vol.34,no.3,pp.311-319,March 1986.
[8]M.Wax and T.Kailath,“Detection of signal by information theoreticcriteria,”IEEE Trans.Acoustics,Speech Signal Processing,vol.33,no.2, pp.387-392,April 1985.
[9]A.Chkeif,K.Abed-Meraim,G.Kawas-Kaleh,and Y.Hua, “Spatio-temporalblind adaptive multiuser detection,”IEEE Trans. Communications,vol.48,no.5,pp.729-732,May 2000.
[10]N.K.Jablon,“Steady state analysis of the generalized sidelobecanceller by adaptive noise cancelling techniques,”IEEE Trans.AntennasPropagation, vol.34,no.3,pp.330-337,March 1986.
[11]Y.H.Chen and A.C.Chang,“A generalized eigenspace-based beamformerwith robust capabilities.”Journal of Chung Cheng Institute of Technology,vol.31,no.1,Oct.2002.
[12]D.R.Brillinger,Time Series:Data Analysis and Theory,expandeded.San Francisco:Holden-Day,1981.
发明内容
本发明设计了一种基于二维空时广义旁瓣相消稳健波束形成算法,其解决的技术问题是如何使用二维数组在广义旁瓣相消器(Generalized Sidelobe Canceller,GSC)架构下,解决指向误差及信号源数目过度估测问题
为了解决上述存在的技术问题,本发明采用了以下方案:
一种基于二维空时广义旁瓣相消稳健波束形成算法,包括以下步骤:
步骤1、当智能天线系统接收器发生指向误差或是信号源过度估测的情况下,将使用空间-时间之特征空间,因此将接收信号之扩频码与导向矢量再表示成时空数组所接收之信号为:
Figure RE-GDA0003851072190000151
其中
Figure RE-GDA0003851072190000152
Figure RE-GDA0003851072190000153
为将空间-时间因素合并后所之噪声矢量;
而时空信号所形成的自相关矩阵表示为:
Figure RE-GDA0003851072190000154
其中,空间-时间信号表示为
Figure RE-GDA0003851072190000155
这边定义空时自相关矩阵二维数组
Figure RE-GDA0003851072190000156
所形成P≡VHV,V可表示为V=[v1,v2,…,vk];此时的自相关矩阵P以Schur形式表示成两个矩阵相乘,例如:P=D·Rc
D包含了空间特征之间的相关性,例如:
Figure RE-GDA0003851072190000157
而这时的Rc为时间处理的自相关矩阵,例如:
Figure RE-GDA0003851072190000158
所以此时将空-时自相关矩阵(1.2) 式之
Figure RE-GDA0003851072190000159
作特征值分解:
Figure RE-GDA00038510721900001510
其中的In为一维度为(MNL-K)×(MNL-K)的单位矩阵,此时的Us是由V所组成的信号子空间所形成的正交基底;而Un为正交于信号子空间的噪声子空间所形成的基底;对角矩阵Λs包含了由
Figure RE-GDA0003851072190000161
分解出的K个最大特征值;
步骤2、对于CDMA接收机而言,v1中的部份内容是能够被知道,扩频码c1是知道的,但导向矢量
Figure RE-GDA0003851072190000162
并不知道;因此这边使用信号子空间正交于噪声子空间的方法;详细地说,就是使用一单位特征矢量
Figure RE-GDA0003851072190000163
来估测
Figure RE-GDA0003851072190000164
此时:
Figure RE-GDA0003851072190000165
上式中IMNL为一个维度MNL×MNL的单位矩阵,由此可看出Q的单位特征矢量能够被拿来当估测正确导向矢量
Figure RE-GDA0003851072190000166
的依据;
此时,令h为由Q找出之单位特征矢量且hi定义为h的第i个元素, i=1,2,…,MNL;由于正确导向矢量ad的第一个元素为1且每个元素的振幅也都为1,且由h可以找到正确之导向矢量,结合上述两条件写成下列解最优值之问题:
Figure RE-GDA0003851072190000167
这边
Figure RE-GDA0003851072190000168
Figure RE-GDA0003851072190000169
第i个元素,经过数学式简化,上式可以重新写成下式:
Figure RE-GDA00038510721900001610
Re{x}代表取x的实数部份,解出式1.6的最优解写成一多项式为
Figure RE-GDA00038510721900001611
Figure RE-GDA00038510721900001612
为它第ith元素;
步骤3、将
Figure RE-GDA00038510721900001613
代入GSC架构中的阻塞矩阵
Figure RE-GDA00038510721900001614
以及权重矢量,得到新的权重矩阵:
Figure RE-GDA0003851072190000171
其中
Figure RE-GDA0003851072190000172
Figure RE-GDA0003851072190000173
为期望求得的适应性权重:
Figure RE-GDA0003851072190000174
上式的
Figure RE-GDA0003851072190000175
经由上述之替换得到之最终权重矢量解为:
Figure RE-GDA0003851072190000176
这时CDMA接收机所接受之信号通过阻塞矩阵能够被阻隔完全,使效能提升。
优选地,UPA数组天线含有M×N个全向性天线单元,各天线单元间距为d=0.5λ并且座落在X-Y平面上的
Figure RE-GDA0003851072190000177
此时 m=1,2,…,M且n=1,2,…,N,假设信号源的入射高低角及水平角分别为θ、φ,假设天线的振幅为单位增益,且相位增益为exp{jπ(m-1)μ+(n-1)υ}、μ=sin(θ)cos(φ)、υ=sin(θ)sin(φ);如GSC的数组天线所示,则第(m,n)个天线单元所接收到的信号表示成:
Figure RE-GDA0003851072190000178
将平面型数组输出之信号表示为:
Figure RE-GDA0003851072190000179
上式中{·}T代表转置、ac(uk)=[1,exp{jπμk},…,exp{jπ(M-1)μk}]T, ark)=[1,exp{jπυk},…,exp{jπ(N-1)υk}]T,N为所有天线所收到之噪声所组成之矩阵,为了使自由度提高,将二维均匀平面型天线数组用一维表示,因此输出信号为:
Figure RE-GDA0003851072190000181
此时
Figure RE-GDA0003851072190000182
为第k个使用者的导引向量,
Figure RE-GDA0003851072190000183
定义为Kronecker乘法,
Figure RE-GDA0003851072190000184
为二维数组展开后与时间特征合并之噪声矩阵,假设第k 个使用者就是所欲使用者k=d,入射角度为(μdd),此时有K-1个干扰, 并且将所欲信号的导引向量表示为ad。此时将
Figure RE-GDA0003851072190000185
以二阶统计量自相关矩阵表示:
Figure RE-GDA0003851072190000186
理论上,在稳态环境下天线数组输出自相关矩阵应由整体平均得到,但于实际操作上,由时间平均来获得估测的自相关矩阵:
Figure RE-GDA0003851072190000187
其中Z为总共观察的位数;
将使用空间-时间之特征空间,因此将接收信号之展频码与导引向量再表示成步骤1中公式1.1时空数组形式:
Figure RE-GDA0003851072190000188
其中,
Figure RE-GDA0003851072190000189
Figure RE-GDA00038510721900001810
为将空间-时间因素合并后所之噪声向量。
一种计算机可读存储介质,其上存储有计算机程序,其特征在于:所述计算机程序被处理器执行时实现上述的基于二维空时广义旁瓣相消稳健波束形成算法。
该基于二维空时广义旁瓣相消稳健波束形成算法具有以下有益效果:
(1)本发明主要使用二维数组在广义旁瓣相消器(Generalized SidelobeCanceller,GSC)架构下,解决指向误差及信号源数目过度估测问题。首先将文献[11]之算法使用在GSC作指向误差校正,但此方法的问题依旧存在。为使校正能力更为稳健本发明提出新方法,基于特征分析处理之稳健空时广义旁瓣相消器(Robust Space-Time GSC,RST-GSC),此方法利用CDMA技术结合二维数组天线,在有指向误差以及信号源过度估测发生时,可自适应让导向矢量有实时校正的能力。
(2)本发明算法是在期望信号在指向误差的情况下,可利用我们提出的空时导向矢量算法去调整自适应权重,使得信号保持住,干扰消除、抑制噪声。
附图说明
图1:现有技术的广义旁瓣相消器架构图;
图2:现有技术完美环境下之二维均匀平面型数组波束方向图;
图3:图2中υ=0,变化μ之一维波束方向图;
图4:图2中令μ=0,变化υ之一维波束方向图;
图5:现有技术指向误差及信号源数目过度估测环境之二维波束方向图;
图6:图5中令υ=0,变化μ之ㄧ维波束方向图;
图7:图5中令μ=0,变化υ之ㄧ维波束方向图;
图8:现有技术稳健广义旁瓣相消器架构图;
图9:现有技术稳健广义旁瓣相消器在指向误差环境之二维波束方向图;
图10:图9中令υ=0,变化μ之一维波束方向图;
图11:图9中令μ=0,变化υ之一维波束方向图;
图12:现有技术稳健广义旁瓣相消器在有指向误差及信号源数目过度估测之二维波束方向图;
图13:图12中令υ=0,变化μ之一维波束方向图;
图14:图12中令μ=0,变化υ之一维波束方向图;
图15:图12中令μ=0,变化υ比较稳健广义旁瓣相消器和广义旁瓣相消器的波束方向图;
图16:图12中变化输入数据之位数来观察输出的SINR值示意图;
图17:图12中变化期望信号之SNR来观察输出的SINR值示意图;
图18:图12中变化指向误差观察输出SINR值示意图;
图19:本发明稳健二维空时广义旁瓣相消器架构图;
图20:本发明RST-GSC令μ=0,变化υ的一维波束方向图;
图21:本发明变化输入数据之位数来观察输出的SINR值示意图;
图22:本发明变化期望信号之SNR来观察输出的SINR值示意图
图23:本发明变化指向误差观察输出SINR值示意图;
图24:本发明变化干扰角度观察输出SINR值示意图。
具体实施方式
下面结合图19至图24,对本发明做进一步说明:
以下公式中出现以下参数的含义如下:
b(t) 数据信号;
bk[i] 第k个用户的第i个数据位
Figure RE-GDA0003851072190000201
匹配滤波器的输出判断;
B 阻塞矩阵;
c(t) 扩频信号波形;
C 所有用户扩频码所组合而成的矩阵;
D 空间特征之间的相关性;
EI 干扰子空间;
En 噪声子空间;
Es 信号子空间;
f 频率;
H 共轭转置;
I 单位矩阵;
K 使用者数目;
L 扩频码长度;
M 平面行数组天线行数;
n(z) 白色高斯噪声;
N 平面数组天线列数;
r(t) 接收端之接收信号;
R 时间处理的自相关矩阵;
Figure RE-GDA0003851072190000211
时间-空间自相关矩阵;
t 时间;
T 转置;
Tc 扩频切片周期;
Tc 码片的周期时间;
T 转换矩阵;
ωo 载波频率;
W 权重向量(矢量);
x(t) 传送端之传送信号;
xk(t) 第k个用户在传送端的输出信号;
Figure RE-GDA0003851072190000221
接收信号;
X(t) 调变后信号;
Xc(t) 扩频信号;
λ 特征值;
σ 变异数;
ω 角频率;
κ 正规化系数;
Λ 特征值矩阵;
εk 第k个使用者的功率;
θ 水平方位角;
φ 高低俯仰角;
μ sin(θ)cos(φ);
υ sin(θ)sin(φ)。
如图19所示稳健二维空时广义旁瓣相消器。
由背景技术中的稳健广义旁瓣相消器可得知,RGSC在有指向误差且信号源数目过度估测的情形下仍有许多缺点,估测并不精确、噪声放大,尤其当SNR提高,性能反而逐渐崩溃。因此,本发明改进后提出一种更加稳健的二维空时估测器。
稳健二维空时广义旁瓣相消器理论基础:二维空时广义旁瓣相消器使用一基于特征空间之二维空时导向矢量估测器结合GSC,架构图如下图19所示。以下为RST-GSC使用基于特征空间之二维空时导向矢量估测器介绍。
由上述式(2.7)可知空时数组系统所接收之信号为:
Figure RE-GDA0003851072190000231
其中空间-时间信号表示为
Figure RE-GDA0003851072190000232
这边定义空时自相关矩阵二维数组
Figure RE-GDA0003851072190000233
所形成P≡VHV,V可表示为V=[v1,v2,…,vk]。此时的自相关矩阵P可以以Schur形式[5]表示成两个矩阵相乘,例如:P=D·Rc
前面的D包含了空间特征之间的相关性,例如:
Figure RE-GDA0003851072190000234
而这时的Rc为时间处理的自相关矩阵,例如:
Figure RE-GDA00038510721900002313
所以此时将空-时自相关矩阵 (1.2)式之
Figure RE-GDA0003851072190000235
作特征值分解:
Figure RE-GDA0003851072190000236
其中的In为一维度为(MNL-K)×(MNL-K)的单位矩阵。此时的Us是由V所组成的信号子空间所形成的正交基底。而Un为正交于信号子空间的噪声子空间所形成的基底。对角矩阵Λs包含了由
Figure RE-GDA0003851072190000237
分解出的K个最大特征值。
对于CDMA接收机而言,v1中的部份内容是可以被知道,扩频码c1是知道的,但导向矢量
Figure RE-GDA00038510721900002312
并不知道。因此这边使用信号子空间正交于噪声子空间的方法。详细地说,就是使用一单位特征矢量
Figure RE-GDA0003851072190000238
来估测
Figure RE-GDA0003851072190000239
此时
Figure RE-GDA00038510721900002310
上式中IMNL为一个维度MNL×MNL的单位矩阵,由此可看出Q的单位特征矢量能够被拿来当估测正确导向矢量
Figure RE-GDA00038510721900002311
的依据。此时,令h为由Q找出之单位特征矢量且hi定义为h的第i个元素,i=1,2,…,MNL。由于正确导向矢量ad的第一个元素为1且每个元素的振幅也都为1,且由h可以找到正确之导向矢量,结合上述两条件可以写成下列解最优值之问题
Figure RE-GDA0003851072190000241
这边
Figure RE-GDA0003851072190000242
Figure RE-GDA0003851072190000243
第i个元素,经过数学式简化,上式可以重新写成下式
Figure RE-GDA0003851072190000244
Re{x}代表取x的实数部份,解出式(3.5)的最优解可以写成一多项式为
Figure RE-GDA0003851072190000245
Figure RE-GDA0003851072190000246
为它第ith元素。
最后将
Figure RE-GDA0003851072190000247
代入GSC架构中的阻塞矩阵
Figure RE-GDA0003851072190000248
以及权重矢量,得到新的权重矩阵:
Figure RE-GDA0003851072190000249
其中
Figure RE-GDA00038510721900002410
Figure RE-GDA00038510721900002411
为所欲求得的适应性权重:
Figure RE-GDA00038510721900002412
上式的
Figure RE-GDA00038510721900002413
经由上述之替换得到之最终权重矢量解为:
Figure RE-GDA00038510721900002414
这时接收机所接受之信号通过阻塞矩阵能够被阻隔完全,使效能提升。
计算机仿真与分析:本发明在GSC架构下使用稳健空时估测方法来校正导向矢量并使用
Figure RE-GDA00038510721900002415
作为波束形成器之权重,仿真实现此估测器是否能准确估出正确导向矢量。假设一8×8的均匀平面型数组,天线之间的距离为0.5λ,此时有4个信号源,分别由角度(μ,υ)=(0,0)、(-0.55,0)、(0.53,0)、(0.57,0)同时入射,而第一个信号源为期望信号,信号强度为SNR=10dB;其他信号源为干扰,干扰信号为等强度INR=30dB。图20、21、22、23及24为在有指向误差(μ,υ)=(0.2,0)且信号源数目过度估测Je=J+1的情形下之模拟图。图20令υ=0,变化μ从-1到+1来观察一维波束方向图,此图并与GSC、RGSC作比较。由图可看出RST-GSC对导向矢量可以有完美估测且不影响干扰的抑制能力,而GSC以及RGSC都有期望信号消除的问题。图21变化输入数据之位数由100到103,藉由输出的SINR来观察不同方法的收敛速度并加入完美环境下的GSC作比较。由此图可看出RST-GSC在不影响收敛速度的情况下,大约在102位就能够收敛,且与完美环境下的GSC有一样稳健的效能。图22变化期望信号之SNR值由0到20来观察不同的SNR 对输出的SINR是否有影响,由于RST-GSC的准确估测,因此效能仍然逼近完美的GSC,而其他两个方法在此环境下性能将会崩溃。图23变化指向误差μ从-2到+2以观察不同方法对抗指向误差之能力,可看出RST-GSC对抗指向误差能力非常好且稳定:而在SNR提高的情况下,可看出RGSC对抗指向误差的能力逐渐衰减。图24变化第三个干扰源的角度μ从-1到+1,藉由输出SINR值评估当干扰由不同角度入射时是否对系统效能有影响。可看出只有当干扰与期望信号由同一角度入射,此时阻塞矩阵会将此干扰信号当成期望信号一同阻隔,造成系统效能崩溃,而当干扰由不同角度入射时,即使双干扰源由同一角度入射对系统也不会有太大的影响。
本发明作了稳健广义旁瓣相消器及稳健空时广义旁瓣相消器的基本理论探讨,且在假设智能型天线系统有信号源数目过度估测、指向误差的环境下作了性能比较。由结果可得知稳健空时广义旁瓣相消器对于导向矢量估测有很稳健的能力,因此在信号源数目过度估测、指向误差很大的环境下系统依然能够保持良好效能。最后由模拟图与本算法作比较可得知可解决问题并提升智能型天线系统效能。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一非易失性计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁盘、光盘、只读存储记忆体(Read-Only Memory,ROM)等。
如此处所使用的对存储器、存储、数据库或其它介质的任何引用可包括非易失性和/或易失性存储器。合适的非易失性存储器可包括只读存储器 (ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM (EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM),它用作外部高速缓冲存储器。作为说明而非局限,RAM以多种形式可得,诸如静态RAM (SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDR SDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)。
上面结合附图对本发明进行了示例性的描述,显然本发明的实现并不受上述方式的限制,只要采用了本发明的方法构思和技术方案进行的各种改进,或未经改进将本发明的构思和技术方案直接应用于其它场合的,均在本发明的保护范围内。

Claims (3)

1.一种基于二维空时广义旁瓣相消稳健波束形成算法,包括以下步骤:
步骤1、当智能天线系统接收器发生指向误差或是信号源过度估测的情况下,将使用空间-时间之特征空间,因此将接收信号之扩频码与导向矢量再表示成时空数组所接收之信号为:
Figure FDA0003764773150000011
其中
Figure FDA0003764773150000012
Figure FDA0003764773150000013
为将空间-时间因素合并后所之噪声矢量;
而时空信号所形成的自相关矩阵表示为:
Figure FDA0003764773150000014
其中,空间-时间信号表示为
Figure FDA0003764773150000015
这边定义空时自相关矩阵二维数组
Figure FDA0003764773150000016
所形成P≡VHV,V可表示为V=[v1,v2,…,vk];此时的自相关矩阵P以Schur形式表示成两个矩阵相乘,例如:P=D·Rc
D包含了空间特征之间的相关性,例如:
Figure FDA0003764773150000017
而这时的Rc为时间处理的自相关矩阵,例如:
Figure FDA0003764773150000018
所以此时将空-时自相关矩阵(1.2)式之
Figure FDA0003764773150000019
作特征值分解:
Figure FDA00037647731500000110
其中的In为一维度为(MNL-K)×(MNL-K)的单位矩阵,此时的Us是由V所组成的信号子空间所形成的正交基底;而Un为正交于信号子空间的噪声子空间所形成的基底;对角矩阵Λs包含了由
Figure FDA00037647731500000111
分解出的K个最大特征值;
步骤2、对于CDMA接收机而言,v1中的部份内容是能够被知道,扩频码c1是知道的,但导向矢量
Figure FDA00037647731500000112
并不知道;因此这边使用信号子空间正交于噪声子空间的方法;详细地说,就是使用一单位特征矢量
Figure FDA00037647731500000113
来估测
Figure FDA00037647731500000114
此时:
Figure FDA00037647731500000115
上式中IMNL为一个维度MNL×MNL的单位矩阵,由此可看出Q的单位特征矢量能够被拿来当估测正确导向矢量
Figure FDA0003764773150000021
的依据;
此时,令h为由Q找出之单位特征矢量且hi定义为h的第i个元素,i=1,2,…,MNL;由于正确导向矢量ad的第一个元素为1且每个元素的振幅也都为1,且由h可以找到正确之导向矢量,结合上述两条件写成下列解最优值之问题:
Figure FDA0003764773150000022
这边
Figure FDA0003764773150000023
Figure FDA0003764773150000024
第i个元素,经过数学式简化,上式可以重新写成下式:
Figure FDA0003764773150000025
Re{x}代表取x的实数部份,解出式1.6的最优解写成一多项式为
Figure FDA0003764773150000026
Figure FDA0003764773150000027
为它第ith元素;
步骤3、将
Figure FDA0003764773150000028
代入GSC架构中的阻塞矩阵
Figure FDA0003764773150000029
以及权重矢量,得到新的权重矩阵:
Figure FDA00037647731500000210
其中
Figure FDA00037647731500000211
Figure FDA00037647731500000212
为期望求得的适应性权重:
Figure FDA00037647731500000213
上式的
Figure FDA00037647731500000214
经由上述之替换得到之最终权重矢量解为:
Figure FDA00037647731500000215
这时CDMA接收机所接受之信号通过阻塞矩阵能够被阻隔完全,使效能提升。
2.根据权利要求1所述的基于二维空时广义旁瓣相消稳健波束形成算法,其特征在于:
UPA数组天线含有M×N个全向性天线单元,各天线单元间距为d=0.5λ并且座落在X-Y平面上的
Figure FDA00037647731500000216
此时m=1,2,…,M且n=1,2,…,N,假设信号源的入射高低角及水平角分别为θ、φ,假设天线的振幅为单位增益,且相位增益为exp{jπ(m-1)μ+(n-1)υ}、μ=sin(θ)cos(φ)、υ=sin(θ)sin(φ);如GSC的数组天线所示,则第(m,n)个天线单元所接收到的信号表示成:
Figure FDA0003764773150000031
将平面型数组输出之信号表示为:
Figure FDA0003764773150000032
上式中{·}T代表转置、ac(uk)=[1,exp{jπμk},…,exp{jπ(M-1)μk}]T,ark)=[1,exp{jπυk},…,exp{jπ(N-1)υk}]T,N为所有天线所收到之噪声所组成之矩阵,为了使自由度提高,将二维均匀平面型天线数组用一维表示,因此输出信号为:
Figure FDA0003764773150000033
此时
Figure FDA0003764773150000034
为第k个使用者的导引向量,
Figure FDA0003764773150000035
定义为Kronecker乘法,
Figure FDA0003764773150000036
为二维数组展开后与时间特征合并之噪声矩阵,假设第k个使用者就是所欲使用者k=d,入射角度为(μdd),此时有K-1个干扰,并且将所欲信号的导引向量表示为ad;此时将
Figure FDA0003764773150000037
以二阶统计量自相关矩阵表示:
Figure FDA0003764773150000038
理论上,在稳态环境下天线数组输出自相关矩阵应由整体平均得到,但于实际操作上,由时间平均来获得估测的自相关矩阵:
Figure FDA0003764773150000039
其中Z为总共观察的位数;
将使用空间-时间之特征空间,因此将接收信号之展频码与导引向量再表示成步骤1中公式1.1时空数组形式:
Figure FDA0003764773150000041
其中,
Figure FDA0003764773150000042
Figure FDA0003764773150000043
为将空间-时间因素合并后所之噪声向量。
3.一种计算机可读存储介质,其上存储有计算机程序,其特征在于:所述计算机程序被处理器执行时实现权利要求1或2所述的基于二维空时广义旁瓣相消稳健波束形成算法。
CN202210882587.8A 2022-07-26 2022-07-26 一种基于二维空时广义旁瓣相消稳健波束形成算法 Pending CN115270869A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210882587.8A CN115270869A (zh) 2022-07-26 2022-07-26 一种基于二维空时广义旁瓣相消稳健波束形成算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210882587.8A CN115270869A (zh) 2022-07-26 2022-07-26 一种基于二维空时广义旁瓣相消稳健波束形成算法

Publications (1)

Publication Number Publication Date
CN115270869A true CN115270869A (zh) 2022-11-01

Family

ID=83770545

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210882587.8A Pending CN115270869A (zh) 2022-07-26 2022-07-26 一种基于二维空时广义旁瓣相消稳健波束形成算法

Country Status (1)

Country Link
CN (1) CN115270869A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117579451A (zh) * 2024-01-15 2024-02-20 成都中科合迅科技有限公司 控制噪声分布的数字相控阵天线多方向的噪声调制方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050254347A1 (en) * 2004-05-14 2005-11-17 Mitel Networks Corporation Parallel gcs structure for adaptive beamforming under equalization constraints
US20090304200A1 (en) * 2008-06-09 2009-12-10 Samsung Electronics Co., Ltd. Adaptive mode control apparatus and method for adaptive beamforming based on detection of user direction sound
CN103399303A (zh) * 2013-07-22 2013-11-20 西安电子科技大学 机载雷达抗密集转发式欺骗干扰方法与系统
CN113504549A (zh) * 2021-07-15 2021-10-15 西安电子科技大学 基于广义旁瓣相消器的导航空时抗干扰方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050254347A1 (en) * 2004-05-14 2005-11-17 Mitel Networks Corporation Parallel gcs structure for adaptive beamforming under equalization constraints
US20090304200A1 (en) * 2008-06-09 2009-12-10 Samsung Electronics Co., Ltd. Adaptive mode control apparatus and method for adaptive beamforming based on detection of user direction sound
CN103399303A (zh) * 2013-07-22 2013-11-20 西安电子科技大学 机载雷达抗密集转发式欺骗干扰方法与系统
CN113504549A (zh) * 2021-07-15 2021-10-15 西安电子科技大学 基于广义旁瓣相消器的导航空时抗干扰方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
AMMAR CHKEIF ET AL: "Spatio-Temporal Blind Adaptive Multiuser Detection", IEEE TRANSACTIONS ON COMMUNICATIONS, vol. 48, no. 5, pages 729 - 732, XP011010836, DOI: 10.1109/26.843180 *
CHENG-CHOU LEE ET AL: "Eigenspace-Based Adaptive Array Beamforming with Robust Capabilities", IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, vol. 45, no. 12, pages 1711 - 1716, XP011003120 *
JHIH CHUNG CHANG: "A Robust Space-Time Generalized Sidelobe Canceller", WIRELESS PERS COMMUN, vol. 70, pages 129 *
JU-HONG LEE ET AL: "Two-Dimensional Adaptive Array Beamforming With Multiple Beam Constraints Using a Generalized Sidelobe Canceller", IEEE TRANSACTIONS ON SIGNAL PROCESSING, vol. 53, no. 9, 30 September 2005 (2005-09-30), pages 3517 - 3529, XP011137860, DOI: 10.1109/TSP.2005.853155 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117579451A (zh) * 2024-01-15 2024-02-20 成都中科合迅科技有限公司 控制噪声分布的数字相控阵天线多方向的噪声调制方法
CN117579451B (zh) * 2024-01-15 2024-04-05 成都中科合迅科技有限公司 控制噪声分布的数字相控阵天线多方向的噪声调制方法

Similar Documents

Publication Publication Date Title
CN105137399B (zh) 基于斜投影滤波的雷达自适应波束形成方法
CN111651719A (zh) 基于收缩估计协方差矩阵重构稳健自适应波束形成方法
CN106646388B (zh) 基于嵌套阵列的mimo雷达抗干扰方法
US10764091B2 (en) Robust adaptive method for suppressing interference in the presence of a signal of interest
CN105158741A (zh) 基于矩阵重构的自适应抗干扰多径多波束形成方法
Mu et al. Robust MVDR beamforming based on covariance matrix reconstruction
Wang et al. Optimum adaptive beamformer design with controlled quiescent pattern by antenna selection
CN106842135B (zh) 基于干扰加噪声协方差矩阵重构的自适应波束形成方法
CN115270869A (zh) 一种基于二维空时广义旁瓣相消稳健波束形成算法
Yu et al. Novel robust beamformers for coherent interference suppression with direction-of-arrival estimation errors
CN112668155B (zh) 一种基于二次重构的稳健波束形成方法及系统
CN116148777A (zh) 一种基于阵列虚拟扩展的阵列雷达主瓣干扰抑制方法
Zhang et al. Adaptive multiple-beamformers for reception of coherent signals with known directions in the presence of uncorrelated interferences
Lee et al. Adaptive array beamforming with robust capabilities under random sensor position errors
Chen et al. Finite data performance analysis of LCMV antenna array beamformers with and without signal blocking
He et al. Two-Dimensional Adaptive Beamforming Based on Atomic-Norm Minimization
CN109633563B (zh) 基于多径信息的自适应相干波束形成方法
Abohamra et al. Direction of Arrival algorithms for user identification in cellular networks
Chen et al. Robust adaptive beamforming based on matched spectrum processing with little prior information
Choi Adaptive nulling beamformer for rejection of coherent and noncoherent interferences
Suleesathira Robust LCSS Beamformer against DOA mismatch
CN117478186B (zh) 基于协方差矩阵估计和子空间构造的波束形成方法
Singh et al. Beamforming Showing Effect on BER with Change in Antenna Configuration
Liu et al. Robust adaptive beamforming against the array pointing error
Lee et al. Optimum beamformers for monopulse angle estimation using overlapping subarrays

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