CN111786707B - 十字天线阵干扰抑制方法和系统 - Google Patents
十字天线阵干扰抑制方法和系统 Download PDFInfo
- Publication number
- CN111786707B CN111786707B CN202010629259.8A CN202010629259A CN111786707B CN 111786707 B CN111786707 B CN 111786707B CN 202010629259 A CN202010629259 A CN 202010629259A CN 111786707 B CN111786707 B CN 111786707B
- Authority
- CN
- China
- Prior art keywords
- signal
- antenna array
- interference
- cross antenna
- incident
- 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
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
- H04B7/0456—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/02—Transmitters
- H04B1/04—Circuits
- H04B1/0475—Circuits with means for limiting noise, interference or distortion
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/06—Receivers
- H04B1/10—Means associated with receiver for limiting or suppressing noise or interference
- H04B1/1081—Reduction of multipath noise
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/0413—MIMO systems
- H04B7/0456—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
- H04B7/046—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting taking physical layer constraints into account
- H04B7/0469—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting taking physical layer constraints into account taking special antenna structures, e.g. cross polarized antennas into account
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
- H04B7/0837—Diversity 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/0842—Weighted combining
- H04B7/0848—Joint weighting
- H04B7/0857—Joint weighting using maximum ratio combining techniques, e.g. signal-to- interference ratio [SIR], received signal strenght indication [RSS]
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Noise Elimination (AREA)
Abstract
本发明提供了一种十字天线阵干扰抑制方法和系统,应用于十字天线阵列,包括:基于十字天线阵列的入射信号的第一信号参数和十字天线阵列的噪声信号的第二信号参数,确定十字天线阵列的信干噪比关于信号加权矢量的目标数学关系;入射信号包括:期望信号和干扰信号;基于目标数学关系,计算在信干噪比为最大值时的信号加权矢量,得到目标加权矢量;利用目标加权矢量对十字天线阵列的各个阵元所接收到的信号进行加权求和,得到目标输出信号。本发明缓解了现有技术中存在的缺乏关于十字天线阵的抗干扰方法的技术问题。
Description
技术领域
本发明涉及卫星通信技术领域,尤其是涉及一种十字天线阵干扰抑制方法和系统。
背景技术
随着无线通信技术的不断发展,信号干扰技术也不断进步,高效地抑制干扰信号对接收有用信号的影响是获得成功无线通信服务的必要条件。针对信号接收机工作环境中可能出现的不同类型的干扰,需要通过多种抗干扰算法对干扰信号的强度进行抑制,保障整个通信过程的顺利进行。当前信号抗干扰广泛采用的原则是:从发射机、接收机、天线、信号处理等通信流程中的不同方面研究特定的针对性抗干扰措施。例如,对于发射机,常常采用改变发射频率、增大发射功率、优化波形设计、采用相干发射等方式抑制干扰信号;接收机则采用改变接收功率、优化接收机布局等手段对干扰信号进行处理;同时,将信号处理技术与天线技术相结合,也是信号干扰的一种有效途径。
阵列天线又称天线阵,是由许多相同的单天线按一定规律排列组成的天线系统,利用阵列天线进行干扰抑制是信号抗干扰领域的重要方式之一。目前,阵列天线自适应抗干扰技术并不完善,一是由于阵列天线的抗干扰性能取决于其采用的自适应算法准则,包括最小均方误差准则(MMSE)、线性约束最小方差准则(LCMV)、最大信干噪比准则(MSINR)等,需要根据不同的抗干扰需求选择不同的自适应算法准则。而现有的阵列信号抗干扰技术主要针对均匀线阵、圆阵、矩形阵等典型阵列信号模型进行设计,缺乏通用的普适阵列信号模型,而没有关于十字天线阵的抗干扰方法。
发明内容
有鉴于此,本发明的目的在于提供一种十字天线阵干扰抑制方法和系统,以缓解了现有技术中存在的缺乏关于十字天线阵的抗干扰方法的技术问题。
第一方面,本发明实施例提供了一种十字天线阵干扰抑制方法,应用于十字天线阵列,所述十字天线阵列包括两个相互垂直的直线阵列,每个方向上包括至少一个阵元,所述方法包括:基于所述十字天线阵列的入射信号的第一信号参数和所述十字天线阵列的噪声信号的第二信号参数,确定所述十字天线阵列的信干噪比关于信号加权矢量的目标数学关系;所述入射信号包括:期望信号和干扰信号,所述第一信号参数包括:入射角,信号强度,所述第二信号参数包括:信号强度;基于所述目标数学关系,计算在所述信干噪比为最大值时的信号加权矢量,得到目标加权矢量;利用所述目标加权矢量对所述十字天线阵列的各个阵元所接收到的信号进行加权求和,得到目标输出信号。
进一步地,基于所述十字天线阵列的入射信号的第一信号参数和所述十字天线阵列的噪声信号的第二信号参数,确定所述十字天线阵列的信干噪比关于信号加权矢量的目标数学关系,包括:基于所述十字天线阵列的入射信号的入射角,计算所述入射信号的方向向量;基于所述方向向量、所述入射信号的信号强度和噪声信号的信号强度,分别计算所述入射信号的自相关矩阵和所述噪声信号的自相关矩阵;基于所述入射信号的自相关矩阵和所述噪声信号的自相关矩阵,确定所述十字天线阵列的信干噪比关于信号加权矢量的目标数学关系。
进一步地,基于所述十字天线阵列的入射信号的入射角,计算所述入射信号的方向向量,包括:获取所述十字天线阵列的阵元间距和所述入射信号的信号波长;基于所述阵元间距、所述信号波长和所述入射角,计算所述十字天线阵列的阵元相对于目标参考点的相位差;所述目标参考点为所述十字天线阵列的两个相互垂直的直线阵列的交叉点;基于所述相位差和所述入射角,计算所述入射信号的方向向量。
进一步地,基于所述方向向量、所述入射信号的信号强度和噪声信号的信号强度,分别计算所述入射信号的自相关矩阵和所述噪声信号的自相关矩阵,包括:通过算式:X(t)=AsXs(t)+AiXi(t)+Xn(t),计算所述十字天线阵所接收的总信号强度为X(t);As表示所述期望信号的方向向量,Xs(t)表示所述期望信号的信号强度,Ai表示所述干扰信号的方向向量,Xi(t)表示所述干扰信号的信号强度,Xn(t)表示所述噪声信号的信号强度;通过算式:Rs=E[AsXs(t)(AsXs(t))H],Ri=E[AiXi(t)(AiXi(t))H]和Rn=E[AnXn(t)(AnXn(t))H],分别计算所述期望信号的自相关矩阵Rs、所述干扰信号的自相关矩阵Ri和所述噪声信号的自相关矩阵Rn;其中,E[]表示计算期望值,()H表示矩阵的共轭转置矩阵。
进一步地,基于所述入射信号的自相关矩阵和所述噪声信号的自相关矩阵,确定所述十字天线阵列的信干噪比关于信号加权矢量的目标数学关系,包括:获取所述十字天线阵列的输出信号的信号强度与信号加权矢量之间的第一数学关系;基于所述入射信号的自相关矩阵和所述第一数学关系,确定所述入射信号的信号功率与所述信号加权矢量之间的第二数学关系;基于所述噪声信号的自相关矩阵和所述第一数学关系,确定所述噪声信号的信号功率与所述信号加权矢量之间的第三数学关系;基于所述第二数学关系和所述第三数学关系,确定所述十字天线阵列的信干噪比关于信号加权矢量的目标数学关系。
进一步地,基于所述目标数学关系,计算在所述信干噪比为最大值时的信号加权矢量,得到目标加权矢量,包括:计算所述期望信号的自相关矩阵相对于干扰噪声矩阵的广义特征值;所述干扰噪声矩阵为所述干扰信号的自相关矩阵与所述噪声信号的自相关矩阵之和;将所述广义特征值的最大值代入到所述目标数学关系,计算得到目标加权矢量。
第二方面,本发明实施例还提供了一种十字天线阵干扰抑制系统,应用于十字天线阵列,所述十字天线阵列包括两个相互垂直的直线阵列,每个方向上包括至少一个阵元,所述系统包括:确定模块,计算模块和加权模块,其中,所述确定模块,用于基于所述十字天线阵列的入射信号的第一信号参数和所述十字天线阵列的噪声信号的第二信号参数,确定所述十字天线阵列的信干噪比关于信号加权矢量的目标数学关系;所述入射信号包括:期望信号和干扰信号,所述第一信号参数包括:入射角,信号强度,所述第二信号参数包括:信号强度;所述计算模块,用于基于所述目标数学关系,计算在所述信干噪比为最大值时的信号加权矢量,得到目标加权矢量;所述加权模块,用于利用所述目标加权矢量对所述十字天线阵列的各个阵元所接收到的信号进行加权求和,得到目标输出信号。
进一步地,所述确定模块包括:第一计算单元,第二计算单元和确定单元,其中,所述第一计算单元,用于基于所述十字天线阵列的入射信号的入射角,计算所述入射信号的方向向量;所述第二计算单元,用于基于所述方向向量、所述入射信号的信号强度和噪声信号的信号强度,分别计算所述入射信号的自相关矩阵和所述噪声信号的自相关矩阵;所述确定单元,用于基于所述入射信号的自相关矩阵和所述噪声信号的自相关矩阵,确定所述十字天线阵列的信干噪比关于信号加权矢量的目标数学关系。
第三方面,本发明实施例还提供了一种电子设备,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述第一方面所述的方法的步骤。
第四方面,本发明实施例还提供了一种具有处理器可执行的非易失的程序代码的计算机可读介质,所述程序代码使所述处理器执行上述第一方面所述方法。
本发明实施例提供的一种十字天线阵干扰抑制方法和系统,通过利用最大信干噪比准则得到目标加权矢量,然后利用目标加权矢量对天线阵元加权求和得到目标输出信号,通过这种方式得到的输出信号,可以使得十字天线阵的信干噪比最大,达到抑制干扰信号的技术效果。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种十字天线阵干扰抑制方法的流程图;
图2为本发明实施例提供的一种十字天线阵列模型示意图;
图3为本发明实施例提供的一种等效阵列关系示意图;
图4为本发明实施例提供的一种四阵元对称十字交叉阵列模型示意图;
图5为本发明实施例提供的一种十字天线阵干扰抑制系统的示意图;
图6为本发明实施例提供的一种十字天线阵干扰抑制系统的确定模块的示意图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例一:
图1是根据本发明实施例提供的一种十字天线阵干扰抑制方法的流程图,该方法应用于十字天线阵列,其中,十字天线阵列包括两个相互垂直的直线阵列,每个方向上包括至少一个阵元,且去掉了中心阵元。如图1所示,该方法具体包括如下步骤:
步骤S102,基于十字天线阵列的入射信号的第一信号参数和十字天线阵列的噪声信号的第二信号参数,确定十字天线阵列的信干噪比关于信号加权矢量的目标数学关系;入射信号包括:期望信号和干扰信号,第一信号参数包括:入射角,信号强度,第二信号参数包括:信号强度。
步骤S104,基于目标数学关系,计算在信干噪比为最大值时的信号加权矢量,得到目标加权矢量。
步骤S106,利用目标加权矢量对十字天线阵列的各个阵元所接收到的信号进行加权求和,得到目标输出信号。
本发明实施例提供的一种十字天线阵干扰抑制方法,通过利用最大信干噪比准则得到目标加权矢量,然后利用目标加权矢量对天线阵元加权求和得到目标输出信号,通过这种方式得到的输出信号,可以使得十字天线阵的信干噪比最大,达到抑制干扰信号的技术效果。
图2是根据本发明实施例提供的一种十字天线阵列模型示意图,如图2所示,十字天线阵列包括两个相互垂直的直线阵列,每个方向上包括至少一个阵元,且去掉十字天线阵列的中心阵元。
如图2所示,十字天线阵列以中心为参考点,顺时针分为Ⅰ、Ⅱ、Ⅲ、Ⅳ四部分,每个方向上分别有M、N、P、R个阵元等距离排列,阵元间距为d,整个十字阵列阵元数目为Ntotal=M+N+P+R。所述阵列模型输入为一个期望信号、Q个干扰信号,期望信号入射角为θ0,干扰信号入射角为θq={θ1,θ2,…θQ},q=1,2,…,Q。
图3是根据本发明实施例提供的一种等效阵列关系示意图。如图3所示,该阵列关系示意图仅表示阵元、输入信号、输出信号的加权关系,不代表阵列空间结构,阵列中共有Ntotal个阵元,分别为Ntotal个阵元接收的期望信号,分别为Ntotal个阵元接收的干扰信号,分别为Ntotal个阵元的权值,加权求和后得到阵列输出。
具体地,步骤S102包括如下步骤:
步骤S1021,基于十字天线阵列的入射信号的入射角,计算入射信号的方向向量。
步骤S1022,基于方向向量、入射信号的信号强度和噪声信号的信号强度,分别计算入射信号的自相关矩阵和噪声信号的自相关矩阵。
步骤S1023,基于入射信号的自相关矩阵和噪声信号的自相关矩阵,确定十字天线阵列的信干噪比关于信号加权矢量的目标数学关系。
在本发明实施例中,首先获取十字天线阵列的阵元间距和入射信号的信号波长;例如,阵元间距为d、入射信号的信号波长为λ。
最后,基于相位差和入射角,计算入射信号的方向向量:
之后,由期望信号的方向向量a(θ0)和干扰信号的方向向量a(θq)组成矩阵As,结合期望信号s0(t)、干扰信号sq(t)和噪声信号n(t),通过算式:计算十字天线阵所接收的总信号强度为X(t);As表示期望信号的方向向量,Xs(t)表示期望信号的信号强度,Ai表示干扰信号的方向向量,Xi(t)表示干扰信号的信号强度,Xn(t)表示噪声信号的信号强度;
通过算式:
分别计算期望信号的自相关矩阵Rs、干扰信号的自相关矩阵Ri和噪声信号的自相关矩阵Rn;其中,E[]表示计算期望值,()H表示矩阵的共轭转置矩阵。
最后,基于入射信号的自相关矩阵和噪声信号的自相关矩阵,确定十字天线阵列的信干噪比关于信号加权矢量的目标数学关系,具体地,包括:
获取十字天线阵列的输出信号的信号强度y(t)与信号加权矢量w=[w1,……wN-1,wN]之间的第一数学关系:
y(t)=wH X(t)=wH AsXs(t)+wH AiXi(t)+WH Xn(t)=S(t)+I(t)+N(t);
基于入射信号的自相关矩阵和第一数学关系,确定入射信号的信号功率与信号加权矢量之间的第二数学关系:
Ps=E[S(t)SH(t)]=E[wH AsXs(t)(wH AsXs(t))H]=wH Rs(wH)H
=wH Rsw;
Pi=E[I(t)IH(t)]=E[wH AiXi(t)(wH AiXi(t))H]=wH Ri(wH)H=wH Riw;
其中,Ps为期望信号的信号功率,Pi为干扰信号的信号功率;
基于噪声信号的自相关矩阵和第一数学关系,确定噪声信号的信号功率与信号加权矢量之间的第三数学关系:
Pn=E[N(t)NH(t)]=E[wH AnXn(t)(wH AnXn(t))H]=wH Rn(wH)H
=wH Rnw;
Pn为噪声信号功率;
基于第二数学关系和第三数学关系,确定十字天线阵列的信干噪比关于信号加权矢量的目标数学关系:
可选地,步骤S104包括如下步骤:
步骤S1041,计算期望信号的自相关矩阵相对于干扰噪声矩阵的广义特征值;干扰噪声矩阵为干扰信号的自相关矩阵与噪声信号的自相关矩阵之和;
步骤S1042,将广义特征值的最大值代入到目标数学关系,计算得到目标加权矢量。
具体地,通过Rsx=λRi+nx计算Rs相对于Ri+n的广义特征值λ1,λ2…,λn,按照大小排列λ1≤λ2≤λ3≤…≤λn,则计算λn对应的广义特征向量x,最终得到最优权矢量(即目标加权矢量)wopt=x,从而通过加权抑制干扰信号。
本发明实施例提供的一种十字天线阵列干扰信号抑制方法,利用所得最优权矢量(即目标加权矢量)wopt对十字天线阵列中的阵元进行加权,可以在抑制干扰信号对十字天线阵列影响的同时使所述阵列天线的输出信干噪比最大。
实施例二:
本发明实施例结合某四阵元对称十字交叉阵列模型的参数,通过具体实施方式对本发明提供的方法作进一步详细的描述。
图4是根据本发明实施例提供的一种四阵元对称十字交叉阵列模型示意图。如图4所示,该阵列模型由十字交叉阵列模型取M=N=P=R=1所得,每个方向上只有一个阵元,阵元与中心间距为整个十字阵列阵元数目为Ntotal=M+N+P+R=4。四阵元十字阵列模型输入为一个方向为θ0=0°的期望信号s0(t)、一个方向为θ1=30°的干扰信号s1(t)、噪声n(t),信号强度
利用实施例一提供的十字天线阵干扰抑制方法,输出四阵元对称十字交叉阵列模型的最优加权系数wopt,具体地,方法流程如下:
然后,在期望信号导向矢量a(θ0)和干扰信号导向矢量a(θ1)的基础上,计算期望信号、干扰信号和噪声信号的自相关矩阵Rs、Ri+n:
λ1=-5.07313644507983×10-15
λ2=-3.31501399804298×10-16
λ3=7.55673266833724×10-17
λ4=31.6901200817828
对应的特征向量:
v1=[-0.11-0.22i,0.48-0.46i,-0.07-0.24i,-0.66+0.00i]T
v2=[0.47+0.11i,0.11-0.38i,0.04-0.48i,0.40+0.00i]T
v3=[-0.26-0.07i,-0.03+0.19i,0.04+0.26i,-0.19+0.00i]T
v4=[-0.61+0.07i,0.28+0.20i,-0.46-0.41i,0.34+0.00i]T
最终得到最优权矢量:wopt=v4=[-0.61+0.07i,0.28+0.20i,-0.46-0.41i,0.34+0.00i]T。
本发明实施例通过最优权矢量wopt对所述四阵元对称十字交叉阵列进行加权,能够达到在降低干扰信号对阵列影响的同时提升整个阵列天线的输出信干噪比。
实施例三:
图5是根据本发明实施例提供的一种十字天线阵干扰抑制系统的示意图,该系统应用于十字天线阵列,十字天线阵列包括两个相互垂直的直线阵列,每个方向上包括至少一个阵元,且去掉了中心阵元。如图5所示,该系统包括:确定模块10,计算模块20和加权模块30。
具体地,确定模块10,用于基于十字天线阵列的入射信号的第一信号参数和十字天线阵列的噪声信号的第二信号参数,确定十字天线阵列的信干噪比关于信号加权矢量的目标数学关系;入射信号包括:期望信号和干扰信号,第一信号参数包括:入射角,信号强度,第二信号参数包括:信号强度。
计算模块20,用于基于目标数学关系,计算在信干噪比为最大值时的信号加权矢量,得到目标加权矢量。
加权模块30,用于利用目标加权矢量对十字天线阵列的各个阵元所接收到的信号进行加权求和,得到目标输出信号。
本发明实施例提供的一种十字天线阵干扰抑制系统,通过利用最大信干噪比准则得到目标加权矢量,然后利用目标加权矢量对天线阵元加权求和得到目标输出信号,通过这种方式得到的输出信号,可以使得十字天线阵的信干噪比最大,达到抑制干扰信号的技术效果。
可选地,图6是根据本发明实施例提供的一种十字天线阵干扰抑制系统的确定模块的示意图。如图6所示,确定模块10包括:第一计算单元11,第二计算单元12和确定单元13。
具体地,第一计算单元11,用于基于十字天线阵列的入射信号的入射角,计算入射信号的方向向量。
第二计算单元12,用于基于方向向量、入射信号的信号强度和噪声信号的信号强度,分别计算入射信号的自相关矩阵和噪声信号的自相关矩阵。
确定单元13,用于基于入射信号的自相关矩阵和噪声信号的自相关矩阵,确定十字天线阵列的信干噪比关于信号加权矢量的目标数学关系。
本发明实施例还提供了一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,处理器执行计算机程序时实现上述实施例一或实施例二中的方法的步骤。
本发明实施例还提供了一种具有处理器可执行的非易失的程序代码的计算机可读介质,程序代码使处理器执行上述实施例一或实施例二中的方法。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (5)
1.一种十字天线阵干扰抑制方法,应用于十字天线阵列,所述十字天线阵列包括两个相互垂直的直线阵列,每个方向上包括至少一个阵元,其特征在于,所述方法包括:
基于所述十字天线阵列的入射信号的第一信号参数和所述十字天线阵列的噪声信号的第二信号参数,确定所述十字天线阵列的信干噪比关于信号加权矢量的目标数学关系;所述入射信号包括:期望信号和干扰信号,所述第一信号参数包括:入射角,信号强度,所述第二信号参数包括:信号强度;
基于所述目标数学关系,计算在所述信干噪比为最大值时的信号加权矢量,得到目标加权矢量;
利用所述目标加权矢量对所述十字天线阵列的各个阵元所接收到的信号进行加权求和,得到目标输出信号;
基于所述十字天线阵列的入射信号的第一信号参数和所述十字天线阵列的噪声信号的第二信号参数,确定所述十字天线阵列的信干噪比关于信号加权矢量的目标数学关系,包括:
基于所述十字天线阵列的入射信号的入射角,计算所述入射信号的方向向量;
基于所述方向向量、所述入射信号的信号强度和噪声信号的信号强度,分别计算所述入射信号的自相关矩阵和所述噪声信号的自相关矩阵;
基于所述入射信号的自相关矩阵和所述噪声信号的自相关矩阵,确定所述十字天线阵列的信干噪比关于信号加权矢量的目标数学关系;
基于所述十字天线阵列的入射信号的入射角,计算所述入射信号的方向向量,包括:
获取所述十字天线阵列的阵元间距和所述入射信号的信号波长;
基于所述阵元间距、所述信号波长和所述入射角,计算所述十字天线阵列的阵元相对于目标参考点的相位差;所述目标参考点为所述十字天线阵列的两个相互垂直的直线阵列的交叉点;
基于所述相位差和所述入射角,计算所述入射信号的方向向量;
基于所述目标数学关系,计算在所述信干噪比为最大值时的信号加权矢量,得到目标加权矢量,包括:
计算所述期望信号的自相关矩阵相对于干扰噪声矩阵的广义特征值;所述干扰噪声矩阵为所述干扰信号的自相关矩阵与所述噪声信号的自相关矩阵之和;
将所述广义特征值的最大值代入到所述目标数学关系,计算得到目标加权矢量;
所述目标数学关系为:SINR表示所述十字天线阵列的信干噪比,Ps为期望信号的信号功率,Pi为干扰信号的信号功率,Pn为噪声信号功率,w为信号加权矢量,RS为期望信号的自相关矩阵、Ri为干扰信号的自相关矩阵,Rn为噪声信号的自相关矩阵,H表示矩阵的共轭转置矩阵;
Ri+n=Ri+Rn,表示干扰噪声矩阵。
2.根据权利要求1所述的方法,其特征在于,基于所述方向向量、所述入射信号的信号强度和噪声信号的信号强度,分别计算所述入射信号的自相关矩阵和所述噪声信号的自相关矩阵,包括:
通过算式:X(t)=AsXs(t)+AiXi(t)+Xn(t),计算所述十字天线阵所接收的总信号强度为X(t);As表示所述期望信号的方向向量,Xs(t)表示所述期望信号的信号强度,Ai表示所述干扰信号的方向向量,Xi(t)表示所述干扰信号的信号强度,Xn(t)表示所述噪声信号的信号强度;
通过算式:Rs=E[AsXs(t)(AsXs(t))H],Ri=E[AiXi(t)(AiXi(t))H]和Rn=E[AnXn(t)(AnXn(t))H],分别计算所述期望信号的自相关矩阵Rs、所述干扰信号的自相关矩阵Ri和所述噪声信号的自相关矩阵Rn;其中,E[]表示计算期望值,()H表示矩阵的共轭转置矩阵。
3.根据权利要求2所述的方法,其特征在于,基于所述入射信号的自相关矩阵和所述噪声信号的自相关矩阵,确定所述十字天线阵列的信干噪比关于信号加权矢量的目标数学关系,包括:
获取所述十字天线阵列的输出信号的信号强度与信号加权矢量之间的第一数学关系;
基于所述入射信号的自相关矩阵和所述第一数学关系,确定所述入射信号的信号功率与所述信号加权矢量之间的第二数学关系;
基于所述噪声信号的自相关矩阵和所述第一数学关系,确定所述噪声信号的信号功率与所述信号加权矢量之间的第三数学关系;
基于所述第二数学关系和所述第三数学关系,确定所述十字天线阵列的信干噪比关于信号加权矢量的目标数学关系;
所述第一数学关系为:
y(t)=wHX(t)=wHAsXs(t)+wHAiXi(t)+wHXn(t);
所述第二数学关系为:
Ps=E[wHAsXs(t)(wHAsXs(t))H]=wHRs(wH)H=wHRsw;
Pi=E[wHAiXi(t)(wHAiXi(t))H]=wHRi(wH)H=wHRiw;
所述第三数学关系为:
Pn=E[wHAnXn(t)(wHAnXn(t))H]=wHRn(wH)H=wHRnw;
y(t)为所述十字天线阵列的输出信号的信号强度。
4.一种电子设备,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现上述权利要求1至3任一项所述的方法的步骤。
5.一种具有处理器可执行的非易失的程序代码的计算机可读介质,其特征在于,所述程序代码使所述处理器执行所述权利要求1-3任一项所述方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010629259.8A CN111786707B (zh) | 2020-07-01 | 2020-07-01 | 十字天线阵干扰抑制方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010629259.8A CN111786707B (zh) | 2020-07-01 | 2020-07-01 | 十字天线阵干扰抑制方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111786707A CN111786707A (zh) | 2020-10-16 |
CN111786707B true CN111786707B (zh) | 2021-10-26 |
Family
ID=72759386
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010629259.8A Active CN111786707B (zh) | 2020-07-01 | 2020-07-01 | 十字天线阵干扰抑制方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111786707B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112904374B (zh) * | 2021-01-29 | 2024-03-19 | 湖南国科微电子股份有限公司 | 一种卫星信号强度评估方法、装置、gnss接收机及介质 |
CN114785426B (zh) * | 2022-03-30 | 2023-11-03 | 西安宇飞电子技术有限公司 | 多天线抗干扰方法、装置、设备和计算机可读存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104360323A (zh) * | 2014-11-01 | 2015-02-18 | 中国民航大学 | 基于十字阵列的ads-b欺骗式干扰抑制方法 |
CN104868946A (zh) * | 2015-06-12 | 2015-08-26 | 哈尔滨工业大学 | 自适应加权的子阵级混合mimo-相控阵系统的干扰抑制方法 |
CN105302936A (zh) * | 2015-08-31 | 2016-02-03 | 中国科学院声学研究所 | 基于相关计算和协方差矩阵重构的自适应波束形成方法 |
CN110429989A (zh) * | 2019-07-22 | 2019-11-08 | 清华大学 | 通信系统间的干扰抑制与共存方法及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
IL218047A (en) * | 2012-02-12 | 2017-09-28 | Elta Systems Ltd | Devices and methods are added, for spatial suppression of interference in wireless networks |
CN103905072B (zh) * | 2012-12-27 | 2016-04-13 | 展讯通信(上海)有限公司 | 无线通信接收机及其干扰消除方法与装置、信号解调方法 |
-
2020
- 2020-07-01 CN CN202010629259.8A patent/CN111786707B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104360323A (zh) * | 2014-11-01 | 2015-02-18 | 中国民航大学 | 基于十字阵列的ads-b欺骗式干扰抑制方法 |
CN104868946A (zh) * | 2015-06-12 | 2015-08-26 | 哈尔滨工业大学 | 自适应加权的子阵级混合mimo-相控阵系统的干扰抑制方法 |
CN105302936A (zh) * | 2015-08-31 | 2016-02-03 | 中国科学院声学研究所 | 基于相关计算和协方差矩阵重构的自适应波束形成方法 |
CN110429989A (zh) * | 2019-07-22 | 2019-11-08 | 清华大学 | 通信系统间的干扰抑制与共存方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN111786707A (zh) | 2020-10-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR100947650B1 (ko) | 간섭 억제를 실현하는 빔 성형 방법 | |
JP5244104B2 (ja) | Mimo通信のための変換領域フィードバック信号伝達 | |
CN111786707B (zh) | 十字天线阵干扰抑制方法和系统 | |
EP1344276B1 (en) | Base station, base station module and method for direction of arrival estimation | |
JP5073809B2 (ja) | チャネル情報予測システム及びチャネル情報予測方法 | |
WO2010088828A1 (zh) | 一种多用户波束成形方法及装置 | |
JP4481336B2 (ja) | チャネル情報予測システム及びチャネル情報予測方法 | |
Rasekh et al. | Design of an adaptive wideband beamforming algorithm for conformal arrays | |
CN107979404A (zh) | 基于虚拟阵列零陷加宽的自适应波束形成方法 | |
US10897299B2 (en) | Wireless apparatus, wireless communication method, and program | |
CN105306117A (zh) | 一种基于协方差矩阵扩展的半虚拟天线阵波束形成方法 | |
CN112543047B (zh) | 一种多波束卫星干扰抑制方法、存储介质及计算设备 | |
CN103248412B (zh) | 一种卫星导航多波束形成的方法 | |
CN117060954A (zh) | 一种基于mimo通信与感知技术的通感一体化波束设计方法 | |
CN102170303B (zh) | 智能天线的波束赋形方法、装置及系统 | |
CN116760439A (zh) | 一种基于雷达通信一体化系统的鲁棒波束成形方法 | |
Krishnamurthy et al. | Sidelobe level distribution computation for antenna arrays with arbitrary element distributions | |
US6999796B2 (en) | Base station, base station module and method for estimating parameters of uplink signals | |
CN111478723A (zh) | 一种自适应计算赋形因子的方法及装置 | |
Zalawadia et al. | Adaptive beam former design using RLS algorithm for smart antenna system | |
CN112350757B (zh) | 一种基于单站阵型映射到mimo阵型的设计方法 | |
Hogade et al. | Analysis of improved and traditional LMS beamforming algorithm for smart antenna | |
Senapati et al. | Performances of some combined algorithms for adaptive beamforming in smart antenna using linear array | |
Tang et al. | New robust adaptive beamforming method for multipath coherent signal reception | |
US20110098067A1 (en) | Mobile communication system, base station and interference removal method |
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 |