CN111786707B - 十字天线阵干扰抑制方法和系统 - Google Patents

十字天线阵干扰抑制方法和系统 Download PDF

Info

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
Application number
CN202010629259.8A
Other languages
English (en)
Other versions
CN111786707A (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.)
Shanghai Qingshen Technology Development Co ltd
Tsinghua University
Original Assignee
Shanghai Qingshen Technology Development Co ltd
Tsinghua 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 Shanghai Qingshen Technology Development Co ltd, Tsinghua University filed Critical Shanghai Qingshen Technology Development Co ltd
Priority to CN202010629259.8A priority Critical patent/CN111786707B/zh
Publication of CN111786707A publication Critical patent/CN111786707A/zh
Application granted granted Critical
Publication of CN111786707B publication Critical patent/CN111786707B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413MIMO systems
    • H04B7/0456Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details 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/02Transmitters
    • H04B1/04Circuits
    • H04B1/0475Circuits with means for limiting noise, interference or distortion
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details 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/06Receivers
    • H04B1/10Means associated with receiver for limiting or suppressing noise or interference
    • H04B1/1081Reduction of multipath noise
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413MIMO systems
    • H04B7/0456Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
    • H04B7/046Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting taking physical layer constraints into account
    • H04B7/0469Selection 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
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/08Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
    • H04B7/0837Diversity 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/0842Weighted combining
    • H04B7/0848Joint weighting
    • H04B7/0857Joint 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={θ12,…θQ},q=1,2,…,Q。
图3是根据本发明实施例提供的一种等效阵列关系示意图。如图3所示,该阵列关系示意图仅表示阵元、输入信号、输出信号的加权关系,不代表阵列空间结构,阵列中共有Ntotal个阵元,
Figure BDA0002565400440000073
分别为Ntotal个阵元接收的期望信号,
Figure BDA0002565400440000074
分别为Ntotal个阵元接收的干扰信号,
Figure BDA0002565400440000075
分别为Ntotal个阵元的权值,加权求和后得到阵列输出。
具体地,步骤S102包括如下步骤:
步骤S1021,基于十字天线阵列的入射信号的入射角,计算入射信号的方向向量。
步骤S1022,基于方向向量、入射信号的信号强度和噪声信号的信号强度,分别计算入射信号的自相关矩阵和噪声信号的自相关矩阵。
步骤S1023,基于入射信号的自相关矩阵和噪声信号的自相关矩阵,确定十字天线阵列的信干噪比关于信号加权矢量的目标数学关系。
在本发明实施例中,首先获取十字天线阵列的阵元间距和入射信号的信号波长;例如,阵元间距为d、入射信号的信号波长为λ。
然后,基于阵元间距的、信号波长λ和入射角θ0,计算十字天线阵列的阵元相对于目标参考点的相位差
Figure BDA0002565400440000071
Figure BDA0002565400440000072
k=1,2,…,max{M,N,P,R};目标参考点为十字天线阵列的两个相互垂直的直线阵列的交叉点。
最后,基于相位差和入射角,计算入射信号的方向向量:
Figure BDA0002565400440000081
q=1,2,…,Q,j为虚数单位。
之后,由期望信号的方向向量a(θ0)和干扰信号的方向向量a(θq)组成矩阵As,结合期望信号s0(t)、干扰信号sq(t)和噪声信号n(t),通过算式:
Figure BDA0002565400440000082
计算十字天线阵所接收的总信号强度为X(t);As表示期望信号的方向向量,Xs(t)表示期望信号的信号强度,Ai表示干扰信号的方向向量,Xi(t)表示干扰信号的信号强度,Xn(t)表示噪声信号的信号强度;
通过算式:
Figure BDA0002565400440000083
Figure BDA0002565400440000084
Figure BDA0002565400440000085
分别计算期望信号的自相关矩阵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为噪声信号功率;
基于第二数学关系和第三数学关系,确定十字天线阵列的信干噪比关于信号加权矢量的目标数学关系:
Figure BDA0002565400440000091
SINR表示十字天线阵列的信干噪比。
可选地,步骤S104包括如下步骤:
步骤S1041,计算期望信号的自相关矩阵相对于干扰噪声矩阵的广义特征值;干扰噪声矩阵为干扰信号的自相关矩阵与噪声信号的自相关矩阵之和;
步骤S1042,将广义特征值的最大值代入到目标数学关系,计算得到目标加权矢量。
具体地,通过Rsx=λRi+nx计算Rs相对于Ri+n的广义特征值λ12…,λn,按照大小排列λ1≤λ2≤λ3≤…≤λn,则
Figure BDA0002565400440000101
计算λn对应的广义特征向量x,最终得到最优权矢量(即目标加权矢量)wopt=x,从而通过加权抑制干扰信号。
其中,
Figure BDA0002565400440000102
表示干扰噪声矩阵。
本发明实施例提供的一种十字天线阵列干扰信号抑制方法,利用所得最优权矢量(即目标加权矢量)wopt对十字天线阵列中的阵元进行加权,可以在抑制干扰信号对十字天线阵列影响的同时使所述阵列天线的输出信干噪比最大。
实施例二:
本发明实施例结合某四阵元对称十字交叉阵列模型的参数,通过具体实施方式对本发明提供的方法作进一步详细的描述。
图4是根据本发明实施例提供的一种四阵元对称十字交叉阵列模型示意图。如图4所示,该阵列模型由十字交叉阵列模型取M=N=P=R=1所得,每个方向上只有一个阵元,阵元与中心间距为
Figure BDA0002565400440000103
整个十字阵列阵元数目为Ntotal=M+N+P+R=4。四阵元十字阵列模型输入为一个方向为θ0=0°的期望信号s0(t)、一个方向为θ1=30°的干扰信号s1(t)、噪声n(t),信号强度
Figure BDA0002565400440000104
利用实施例一提供的十字天线阵干扰抑制方法,输出四阵元对称十字交叉阵列模型的最优加权系数wopt,具体地,方法流程如下:
根据所述阵列模型间距
Figure BDA0002565400440000105
计算某一阵元相对于参考点的相位差
Figure BDA0002565400440000106
根据相位差
Figure BDA0002565400440000107
及期望信号入射角θ0和干扰信号入射角θ1,计算方向向量a(θ0)、a(θ1):
Figure BDA0002565400440000111
然后,在期望信号导向矢量a(θ0)和干扰信号导向矢量a(θ1)的基础上,计算期望信号、干扰信号和噪声信号的自相关矩阵Rs、Ri+n
Figure BDA0002565400440000112
Figure BDA0002565400440000113
结合信号强度
Figure BDA0002565400440000114
根据(Rs-λRi+n)x=0计算Rs相对于Ri+n的广义特征值λ1≤λ2≤λ3≤λ4
λ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.一种十字天线阵干扰抑制方法,应用于十字天线阵列,所述十字天线阵列包括两个相互垂直的直线阵列,每个方向上包括至少一个阵元,其特征在于,所述方法包括:
基于所述十字天线阵列的入射信号的第一信号参数和所述十字天线阵列的噪声信号的第二信号参数,确定所述十字天线阵列的信干噪比关于信号加权矢量的目标数学关系;所述入射信号包括:期望信号和干扰信号,所述第一信号参数包括:入射角,信号强度,所述第二信号参数包括:信号强度;
基于所述目标数学关系,计算在所述信干噪比为最大值时的信号加权矢量,得到目标加权矢量;
利用所述目标加权矢量对所述十字天线阵列的各个阵元所接收到的信号进行加权求和,得到目标输出信号;
基于所述十字天线阵列的入射信号的第一信号参数和所述十字天线阵列的噪声信号的第二信号参数,确定所述十字天线阵列的信干噪比关于信号加权矢量的目标数学关系,包括:
基于所述十字天线阵列的入射信号的入射角,计算所述入射信号的方向向量;
基于所述方向向量、所述入射信号的信号强度和噪声信号的信号强度,分别计算所述入射信号的自相关矩阵和所述噪声信号的自相关矩阵;
基于所述入射信号的自相关矩阵和所述噪声信号的自相关矩阵,确定所述十字天线阵列的信干噪比关于信号加权矢量的目标数学关系;
基于所述十字天线阵列的入射信号的入射角,计算所述入射信号的方向向量,包括:
获取所述十字天线阵列的阵元间距和所述入射信号的信号波长;
基于所述阵元间距、所述信号波长和所述入射角,计算所述十字天线阵列的阵元相对于目标参考点的相位差;所述目标参考点为所述十字天线阵列的两个相互垂直的直线阵列的交叉点;
基于所述相位差和所述入射角,计算所述入射信号的方向向量;
基于所述目标数学关系,计算在所述信干噪比为最大值时的信号加权矢量,得到目标加权矢量,包括:
计算所述期望信号的自相关矩阵相对于干扰噪声矩阵的广义特征值;所述干扰噪声矩阵为所述干扰信号的自相关矩阵与所述噪声信号的自相关矩阵之和;
将所述广义特征值的最大值代入到所述目标数学关系,计算得到目标加权矢量;
所述目标数学关系为:
Figure FDA0003172582870000021
SINR表示所述十字天线阵列的信干噪比,Ps为期望信号的信号功率,Pi为干扰信号的信号功率,Pn为噪声信号功率,w为信号加权矢量,RS为期望信号的自相关矩阵、Ri为干扰信号的自相关矩阵,Rn为噪声信号的自相关矩阵,H表示矩阵的共轭转置矩阵;
通过Rsx=λRi+nx计算Rs相对于Ri+n的广义特征值λ1,λ2…,λn,按照大小排列λ1≤λ2≤λ3≤…≤λn,则
Figure FDA0003172582870000022
计算λn对应的广义特征向量x,得到目标加权矢量wopt=x;
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任一项所述方法。
CN202010629259.8A 2020-07-01 2020-07-01 十字天线阵干扰抑制方法和系统 Active CN111786707B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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 展讯通信(上海)有限公司 无线通信接收机及其干扰消除方法与装置、信号解调方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
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