CN113630355B - 基于空时功率倒置阵的宽带干扰抑制装置及抑制方法 - Google Patents
基于空时功率倒置阵的宽带干扰抑制装置及抑制方法 Download PDFInfo
- Publication number
- CN113630355B CN113630355B CN202111184393.2A CN202111184393A CN113630355B CN 113630355 B CN113630355 B CN 113630355B CN 202111184393 A CN202111184393 A CN 202111184393A CN 113630355 B CN113630355 B CN 113630355B
- Authority
- CN
- China
- Prior art keywords
- unit
- signal
- delay
- weight
- adder
- 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
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L25/03012—Arrangements for removing intersymbol interference operating in the time domain
- H04L25/03019—Arrangements for removing intersymbol interference operating in the time domain adaptive, i.e. capable of adjustment during data reception
- H04L25/03057—Arrangements for removing intersymbol interference operating in the time domain adaptive, i.e. capable of adjustment during data reception with a recursive structure
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L2025/03592—Adaptation methods
- H04L2025/03598—Algorithms
- H04L2025/03611—Iterative algorithms
- H04L2025/03636—Algorithms using least mean square [LMS]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Noise Elimination (AREA)
Abstract
本发明提供一种基于空时功率倒置阵的宽带干扰抑制装置及抑制方法,包括空时二维功率倒置阵结构及重定时延时最小均方滤波器;空时二维功率倒置阵结构由M个阵元组成,设第一个阵元为参考阵元,其他M‑1个阵元为权值阵元;重定时延时最小均方滤波器包括M‑1组权值处理模块、两个单位延时器、第一加法器、一个单位延时器及第二加法器;在每个天线阵元后加上若干个延迟单元,使得每个阵元后构成一个FIR滤波器,从而使得阵列具有一定的频率分辨能力,通过增加FIR滤波器阶数,大大提高了功率倒置阵的自由度,从而有效提高了干扰抑制的带宽。将时域横向滤波器的总延时时间大于最大的多径传播延时,可实现多径环境下的干扰抑制。
Description
技术领域
本发明涉及信号处理领域技术领域,具体涉及一种基于空时功率倒置阵的宽带干扰抑制装置及抑制方法。
背景技术
无线通信电台广泛应用于野外科考、应急救灾、舰艇飞机通信等军民用领域,具有信号稳定、安全性高等优点,也是保障长距离视距通信的重要设备。然而,在复杂电磁环境下,通信电台易受共平台发射机或敌方干扰机产生的多样式干扰影响,通信鲁棒性较差。特别是在扩频通信体制下,需要对扩频带宽内的宽带干扰信号进行抑制。
为解决电台在扩频通信体制下的同频宽带干扰问题,多接收天线技术得到了广泛的应用。由于通信信道通常是衰落的且各向异性,无法采用阵列信号处理中的波束形成技术,而功率倒置算法可利用自适应调零天线接收信号,在强干扰来向形成波束零陷,从而实现空域滤波抑制干扰。功率倒置阵直接将阵列的输出作为误差信号,在保证参考路权系数恒定的前提下,通过自适应算法调整权值路的加权系数,使总的输出功率最小,具有工程实现价值。
然而,在实际环境中,信号传播通常具有很强的多径效应,造成信道较强的频率选择性,宽带干扰信号的不同频点对应信道的响应也不同。而目前使用的功率倒置天线阵列在时域上通常只有一个抽头,只能针对单个频点的信道进行干扰抑制,无法覆盖整个宽带信号的带宽,造成传统对消算法性能下降。从方向图角度考虑,基于传统空域功率倒置阵的自适应抗干扰算法只可以在干扰方向形成较深的零陷,进而在空域维度上对干扰进行有效的抑制,但是只能分辨来自不同方向的信号,而无法分辨相同方向不同频率的信号,因此,对多径场景下的宽带干扰,其抑制能力十分有限。
中国专利基于波形重构的宽带干扰抑制方法(申请号CN201710277014)、一种盲波束宽带干扰抑制方法及装置(申请号CN201510194436.3)、最低误码率准则的超宽带干扰抑制技术(申请号CN201010510882.8)、一种空频抗干扰方法及装置(申请号CN201611008389.X)等或采用基于时域积累的空时处理手段,或建立非线性模型,利用时域重构完成宽带干扰抑制。均存在算法实时性较差、无法高速流水运行、不利于工程实现等缺点,始终无法克服多径信道下大带宽通信干扰的抑制难题。
发明内容
本发明的目的是针对上述技术问题,提供一种具有高速、低延时性能且能够有效抑制宽带干扰的基于空时功率倒置阵的宽带干扰抑制装置及抑制方法。
为实现上述目的,本发明所设计的基于空时功率倒置阵的宽带干扰抑制装置,包括空时二维功率倒置阵结构及重定时延时最小均方滤波器;所述空时二维功率倒置阵结构由M个阵元组成,设第一个阵元为参考阵元,其他M-1个阵元为权值阵元,每个权值阵元各有N个单位延时器,每个单位延时器都有1个复数权值,使得每个权值阵元构成一个阵元FIR滤波器;
重定时延时最小均方滤波器包括M-1组权值处理模块、两个单位延时器、第一加法器、一个单位延时器及第二加法器;两个单位延时器的输入端连接参考阵元接收信号,每个权值处理模块一路输入端一一对应连接权值阵元接收信号、每个权值处理模块另一路输入端一一对应连接一个初始信号,每组权值处理模块一路输出端连接三个单位延时器的输入端、三个单位延时器的输出端连接权值处理模块;所有权值处理模块另一路输出端均连接第二加法器的输入端,第二加法器的输出端和两个单位延时器的输出端均连接第一加法器的输入端,第一加法器的输出端连接一个单位延时器的输入端,一个单位延时器的输入端分别与每个权值处理模块的输入端相连。
进一步地,所述参考阵元的权值始终为1,M个阵元的接收信号分别记为x0(n)、x1(n)、……、xM-1(n),其中,x0 (n)表示参考阵元接收信号,x1(n) 、……、xM-1(n)为权值阵元接收信号,n表示当前时刻;在每个权值阵元后加上N个单位延时器Z-1,记为xi (n)、xi (n-1)、……、xi (n-N+1),其中,i=0、……、M-1,每个单位延时器的复数权值记为wij,i=0、……、M-1,j=1、……、N,其中i表示阵元序号,j表示延时单元序号。
进一步地,每个所述权值阵元接收信号xi(n)输送至对应的权值处理模块后经过三个单位延时器形成延时信号xi(n-3D)给对应的权值处理模块输出;每个初始信号yi(n-1D)输送至对应的权值处理模块后形成延时信号yi(n-2D),所有的延时信号yi(n-2D)经过第二加法器产生权值合并输出信号y(n-2D),权值合并输出信号y(n-2D)和参考阵元接收信号x0(n)经过第一加法器形成误差信号e(n-2D),误差信号e(n-2D)经过一个单位延时器形成延时误差信号e(n-3D)输送至每个权值处理模块输出。
进一步地,每组所述权值处理模块包括四个乘加单元、两个单位延时器、第一加法器及第二加法器,四个乘加单元分别为第一乘加单元、第二乘加单元、第三乘加单元、第四乘加单元;第一乘加单元和第三乘加单元均包括上支路乘法器、下支路乘法器和加法器,上支路乘法器和下支路乘法器的输出端均与加法器的输入端相连;第二乘加单元和第四乘加单元均包括左支路乘法器、右支路乘法器和加法器,左支路乘法器和右支路乘法器的输出端均与加法器的输入端相连。
进一步地,每组所述权值处理模块的输入信号为IQ信号xiI(n)、xiQ(n),i=1、……、M-1,IQ信号xiI(n)、xiQ(n)一路经过三个单位延时器形成三个单位延时信号xiI(n-3D)、xiQ(n-3D),三个单位延时信号xiI(n-3D)、xiQ(n-3D)经过两个单位延时器形成五个单位延时信号xiI(n-5D)、xiQ(n-5D)输出;延时信号xiI(n-3D)分别连接第一乘加单元的上支路乘法器、第三乘加单元的上支路乘法器,延时信号xiQ(n-3D)分别连接第一乘加单元的下支路乘法器、第三乘加单元的下支路乘法器,第一乘加单元的加法器输出端依次经过一个单位延时器、上支路移位器、加法器、一个单位延时器与第二乘加单元的左支路乘法器和右支路乘法器输入端均相连,同理,第三乘加单元的加法器输出端依次经过一个单位延时器、下支路移位器、加法器、一个单位延时器与第四乘加单元的左支路乘法器和右支路乘法器输入端均相连;信号xiI(n)分别连接第三乘加单元的左支路乘法器和第四乘加单元的左支路乘法器,信号xiQ(n) 分别连接第三乘加单元的右支路乘法器和第四乘加单元的右支路乘法器;初始信号yiI(n-1D)输送至第一加法器的左支路,同时第二乘加单元加法器的输出端经过一个单位延时器输送至第一加法器的右支路,第一加法器左右支路相加后再经过一个单位延时器形成延时信号yiI (n-2D);初始信号yiQ(n-1D)输送至第二加法器左支路,同时第四乘加单元加法器的输出端经过一个单位延时器输送至第二加法器右支路,第二加法器左右支路相加后再经过一个单位延时器形成延时信号yiQ(n-2D)。
进一步地,每组所述权值处理模块的延时信号yiI(n-2D)经过第一加法器形成权值合并输出信号yI(n-2D),每个权值处理模块的延时信号yiQ(n-2D)经过第一加法器形成权值合并输出信号yQ(n-2D)。权值合并输出信号yI(n-2D)和参考阵元接收信号x0I (n)经过第一加法器形成误差信号eI (n-2D),误差信号eI(n-2D)经过一个单位延时器形成延时误差信号eI (n-3D)输送至每个权值处理模块,延时误差信号eI(n-3D)一路分别连接第一乘加单元的上支路乘法器和第三乘加单元的上支持乘法器,延时误差信号eI(n-3D)另一路经过权值处理模块中的一个单位延时器形成延时误差信号eI(n-4D),延时误差信号eI(n-4D)输出;同理,权值合并输出信号yQ (n-2D)和参考阵元接收信号x0Q (n)经过第一加法器形成误差信号eQ(n-2D),误差信号eQ(n-2D)经过一个单位延时器形成延时误差信号eQ(n-3D)输送至每个权值处理模块,延时误差信号eQ(n-3D) 一路分别连接第一乘加单元的下支路乘法器和第三乘加单元的下支持乘法器,延时误差信号eQ(n-3D)另一路经过权值处理模块中的一个单位延时器形成延时误差信号eQ(n-4D)输出。
还提供一种上述所述基于空时功率倒置阵的宽带干扰抑制装置的抑制方法,其特征在于:所述抑制方法如下:
S1)初始化空时二维功率倒置阵模型,包括以下子步骤:
S11)空时二维功率倒置阵结构由M个阵元组成,设第一个阵元为参考阵元1,其他M-1个阵元为权值阵元,M个阵元的接收信号分别记为x0(n)、x1 (n)、……、xM-1(n),其中,x0(n)表示参考阵元接收信号,x1(n) 、……、xM-1(n)为权值阵元接收信号,n表示当前时刻;
S12)每个权值阵元各有N个单位延时器,记为xi (n)、xi (n-1) 、……、xi (n-N+1),其中,i=0、……、M-1,每个单位延时器的复数权值记为wij,i=0、……、M-1,j=1、……、N,其中i表示阵元序号,j表示延时单元序号;
S13)初始化参考阵元x0路权值为w01=1,w0j=0,其中j=2、……、N,初始化权值路x1、……、xM-1的权值为
wij=0,其中i=1、……、M-1,j=1、……、N
记输入信号向量为
x(n)=[x1 (n)⋯x1 (n-N+1)⋯xM-1 (n)⋯xM-1(n-N+1)]T∈C(M-1)×N
记权向量为
w=[w1,1 w1,2⋯w1,N⋯wM-1,1 wM-1,2⋯wM-1,N ]T∈C(M-1)×N
初始化当前时刻为n=1;
S2)根据多径信道环境计算处理延时量,将参考阵元接收信号进行延时,将x0路信号x0 (n)延时2D,得到x0 (n-2D),其中2D要大于最大多径传播时延;
S3)利用权值处理模块计算输出信号、误差信号以及权向量的更新,具体为:
S31)对于n=1,2,3,……,利用当前权值w,计算输出信号为
y(n)=w1,1 * x1 (n)+w1,2 * x1 (n-1)+⋯+w1,N * x1 (n-N+1)+⋯+wM-1,1 * xM-1 (n)+wM-1,2 *xM-1 (n-1)+⋯+wM-1,N * xM-1 (n-N+1)=wH x(n)
计算期望延时信号为
d(n-2D) = -x0 (n-2D)
计算延时误差信号为
e(n-2D) = d(n-2D)-y(n-2D)
S32)设置步长μ,对于空时二维功率倒置阵,μ由移位器实现,μ为2-8~2-12,并进行权值的更新,权向量更新为:
w(n+1) = w(n)+μx(n-2D)e* (n-2D)
S4)令n=n+1,转子步骤S31),直到满足‖w(n+1)-w(n)‖2 2≤ε,其中ε为迭代终止阈值;
S5)为实现自适应宽带干扰抑制,对于当前时刻n,输出信号为y ̃(n)=x0(n)+wH x(n)。
与现有技术相比,本发明具有以下优点:
1)有效提高干扰抑制的带宽。在每个天线阵元后加上若干个延迟单元,使得每个阵元后构成一个FIR滤波器,从而使得阵列具有一定的频率分辨能力,通过增加FIR滤波器阶数,大大提高了功率倒置阵的自由度,从而有效提高了干扰抑制的带宽。
2)能够实现多径等复杂信道下的宽带干扰抑制。将时域多抽头的思想引入到空域天线阵列中,综合考虑接收机的工作环境,将时域横向滤波器的总延时时间大于最大的多径传播延时,可实现多径环境下的干扰抑制。
3)有效提高干扰抑制速度。相比于传统LMS算法,当权系数越多时,系统吞吐能力越低,而本发明提出的TF-RDLMS结构,权值是局部更新的,在硬件实现中可采用高度流水的结构,可以通过简单增加权值处理模块来提高空时滤波器的阶数,因此,本发明适合基于FPGA的高速低延时硬件实现。
附图说明
图1为本发明基于空时功率倒置阵的宽带干扰抑制装置原理图;
图2为图1中重定时延时最小均方滤波器原理图;
图3为图2中权值处理模块原理图;
图4为载频150MHz带宽10MHz宽带干扰信号频谱;
图5为传统功率倒置阵干扰抑制效果频谱图;
图6为本发明宽带干扰抑制效果频谱图。
具体实施方式
下面结合附图和具体实例对本发明作进一步的详细说明。
本发明基于空时功率倒置阵的宽带干扰抑制装置,包括空时二维功率倒置阵结构及重定时延时最小均方滤波器(TF-RDLMS)。
其中,如图1所示空时二维功率倒置阵结构由M个阵元组成,设第一个阵元为参考阵元,其他M-1个阵元为权值阵元,参考阵元的权值始终为1,每个权值阵元各有N个单位延时器,每个单位延时器都有1个复数权值,使得每个权值阵元构成一个阵元FIR滤波器。M个阵元的接收信号分别记为x0(n)、x1 (n)、……、xM-1(n),其中,x0 (n)表示参考阵元接收信号,x1(n) 、……、xM-1(n)为权值阵元接收信号,n表示当前时刻;在每个权值阵元后加上N个单位延时器Z-1,记为xi (n)、xi (n-1) 、……、xi (n-N+1),其中,i=0、……、M-1,每个单位延时器的复数权值记为wij,i=0、……、M-1,j=1、……、N,其中i表示阵元序号,j表示延时单元序号。
结合图2所示,重定时延时最小均方滤波器(TF-RDLMS)包括M-1组权值处理模块04、两个单位延时器01、第一加法器02、一个单位延时器03及第二加法器06。两个单位延时器01的输入端连接参考阵元接收信号x0(n),每个权值处理模块04一路输入端一一对应连接权值阵元接收信号xi(n)、每个权值处理模块04另一路输入端一一对应连接一个初始信号yi(n-1D),其中i=1、……、M-1,每组权值处理模块04一路输出端连接三个单位延时器05的输入端、三个单位延时器05的输出端连接权值处理模块04;所有权值处理模块04另一路输出端均连接第二加法器06的输入端,第二加法器06的输出端和两个单位延时器01的输出端均连接第一加法器02的输入端,第一加法器02的输出端连接一个单位延时器03的输入端,一个单位延时器03的输入端分别与每个权值处理模块04的输入端相连。每个权值阵元接收信号xi(n)输送至对应的权值处理模块04后经过三个单位延时器05形成延时信号xi(n-3D)给对应的权值处理模块04输出;每个初始信号yi(n-1D)输送至对应的权值处理模块04后形成延时信号yi(n-2D),所有的延时信号yi(n-2D)经过第二加法器06产生权值合并输出信号y(n-2D),权值合并输出信号y(n-2D)和参考阵元接收信号x0(n)经过第一加法器02形成误差信号e(n-2D),误差信号e(n-2D)经过一个单位延时器03形成延时误差信号e(n-3D)输送至每个权值处理模块04输出。
结合图3所示,每组权值处理模块04包括四个乘加单元、两个单位延时器041、第一加法器048及第二加法器049,四个乘加单元分别为第一乘加单元042、第二乘加单元044、第三乘加单元045、第四乘加单元047。第一乘加单元042和第三乘加单元045均包括上支路乘法器、下支路乘法器和加法器,上支路乘法器和下支路乘法器的输出端均与加法器的输入端相连;第二乘加单元044和第四乘加单元047均包括左支路乘法器、右支路乘法器和加法器,左支路乘法器和右支路乘法器的输出端均与加法器的输入端相连。
每组权值处理模块04的输入信号为IQ信号xiI(n)、xiQ(n)(i=1、……、M-1),IQ信号xiI(n)、xiQ(n)一路经过三个单位延时器05形成三个单位延时信号xiI(n-3D)、xiQ(n-3D),三个单位延时信号xiI(n-3D)、xiQ(n-3D)经过两个单位延时器041形成五个单位延时信号xiI(n-5D)、xiQ(n-5D)输出;延时信号xiI(n-3D)分别连接第一乘加单元042的上支路乘法器、第三乘加单元045的上支路乘法器,延时信号xiQ(n-3D)分别连接第一乘加单元042的下支路乘法器、第三乘加单元045的下支路乘法器,第一乘加单元042的加法器输出端依次经过一个单位延时器、上支路移位器043、加法器、一个单位延时器与第二乘加单元044的左支路乘法器和右支路乘法器输入端均相连,同理,第三乘加单元045的加法器输出端依次经过一个单位延时器、下支路移位器046、加法器、一个单位延时器与第四乘加单元047的左支路乘法器和右支路乘法器输入端均相连;信号xiI(n)分别连接第三乘加单元044的左支路乘法器和第四乘加单元047的左支路乘法器,信号xiQ(n) 分别连接第三乘加单元044的右支路乘法器和第四乘加单元047的右支路乘法器;初始信号yiI(n-1D)输送至第一加法器048的左支路,同时第二乘加单元044加法器的输出端经过一个单位延时器输送至第一加法器048的右支路,第一加法器048左右支路相加后再经过一个单位延时器形成延时信号yiI (n-2D);初始信号yiQ(n-1D)输送至第二加法器049左支路,同时第四乘加单元047加法器的输出端经过一个单位延时器输送至第二加法器049右支路,第二加法器049左右支路相加后再经过一个单位延时器形成延时信号yiQ (n-2D)。
每个权值处理模块04的延时信号yiI(n-2D)经过第一加法器02形成权值合并输出信号yI(n-2D),每个权值处理模块04的延时信号yiQ(n-2D)经过第一加法器02形成权值合并输出信号yQ(n-2D)。权值合并输出信号yI(n-2D)和参考阵元接收信号x0I (n)经过第一加法器02形成误差信号eI (n-2D),误差信号eI(n-2D)经过一个单位延时器03形成延时误差信号eI (n-3D)输送至每个权值处理模块04,延时误差信号eI(n-3D)一路分别连接第一乘加单元042的上支路乘法器和第三乘加单元045的上支持乘法器,延时误差信号eI(n-3D)另一路经过权值处理模块04中的一个单位延时器形成延时误差信号eI(n-4D),延时误差信号eI(n-4D)输出;同理,权值合并输出信号yQ (n-2D)和参考阵元接收信号x0Q (n)经过第一加法器02形成误差信号eQ(n-2D),误差信号eQ(n-2D)经过一个单位延时器03形成延时误差信号eQ(n-3D)输送至每个权值处理模块04,延时误差信号eQ(n-3D) 一路分别连接第一乘加单元042的下支路乘法器和第三乘加单元045的下支持乘法器,延时误差信号eQ(n-3D)另一路经过权值处理模块04中的一个单位延时器形成延时误差信号eQ(n-4D)输出。
基于空时功率倒置阵的宽带干扰抑制方法具体如下:
S1)初始化空时二维功率倒置阵模型,包括以下子步骤:
S11)空时二维功率倒置阵结构由M个阵元组成,设第一个阵元为参考阵元1,其他M-1个阵元为权值阵元,M个阵元的接收信号分别记为x0(n)、x1 (n)、……、xM-1(n),其中,x0(n)表示参考阵元接收信号,x1(n) 、……、xM-1(n)为权值阵元接收信号,n表示当前时刻;
S12)每个权值阵元各有N个单位延时器,记为xi (n)、xi (n-1) 、……、xi (n-N+1),其中,i=0、……、M-1,每个单位延时器的复数权值记为wij,i=0、……、M-1,j=1、……、N,其中i表示阵元序号,j表示延时单元序号;
S13)初始化参考阵元x0路权值为w01=1,w0j=0,其中j=2、……、N,初始化权值路x1、……、xM-1的权值为
wij=0,其中i=1、……、M-1,j=1、……、N
记输入信号向量为
x(n)=[x1 (n)⋯x1 (n-N+1)⋯xM-1 (n)⋯xM-1(n-N+1)]T∈C(M-1)×N
记权向量为
w=[w1,1 w1,2⋯w1,N⋯wM-1,1 wM-1,2⋯wM-1,N ]T∈C(M-1)×N
初始化当前时刻为n=1;
S2)根据多径信道环境计算处理延时量,将参考阵元接收信号进行延时,将x0路信号x0 (n)延时2D,得到x0 (n-2D),其中2D要大于最大多径传播时延;
S3)利用权值处理模块(04)计算输出信号、误差信号以及权向量的更新,具体为:
S31)对于n=1,2,3,……,利用当前权值w,计算输出信号为
y(n)=w1,1 * x1 (n)+w1,2 * x1 (n-1)+⋯+w1,N * x1 (n-N+1)+⋯+wM-1,1 * xM-1 (n)+wM-1,2 *xM-1 (n-1)+⋯+wM-1,N * xM-1 (n-N+1)=wH x(n)
计算期望延时信号为
d(n-2D) = -x0 (n-2D)
计算延时误差信号为
e(n-2D) = d(n-2D)-y(n-2D)
S32)设置步长μ,对于空时二维功率倒置阵,μ由移位器实现,μ为2-8~2-12,并进行权值的更新,权向量更新为:
w(n+1) = w(n)+μx(n-2D)e* (n-2D)
S4)令n=n+1,转子步骤S31),直到满足‖w(n+1)-w(n)‖2 2≤ε,其中ε为迭代终止阈值;
S5)为实现自适应宽带干扰抑制,对于当前时刻n,输出信号为y ̃(n)=x0(n)+wH x(n)。
上述中,H表示共轭转置、*表示共轭。
图4~6展示了本发明提出的宽带干扰抑制方法在AD9361开发板FPGA上实现的效果。为了更好的观察干扰抑制效果,没有加入有用信号分量。其中宽带干扰信号是10MHz带宽的BPSK调制信号,中心频率为150MHz。图4给出了原始10MHz宽带干扰信号频谱,干扰信号在整个带内具有-70dBm的功率谱。当采用传统功率倒置算法时,在干扰中心频点获得了20dB的干扰对消效果,而在频谱两侧仅有10dB干扰抑制效果,如图5所示。而当采用本发明提出的宽带干扰抑制方法时,整个带内干扰均获得了40dB以上的干扰抑制效果,如图6所示。需要注意的是频点1(150.1MHz)为开发板AD9361的本振泄露,实际应用中可以通过设置上下变频频点避开有用信号,例如当有用信号工作在150MHz时将本振设置为151MHz,前端滤波器可以有效滤除本振信号。
最后所应该说明的是:以上具体的实施步骤主要是用来说明本发明的具体实现过程,而非限制本发明的技术方案,尽管该部分针对具体的实现实例进行了详细的说明,该领域的技术人员都应该理解:依然可以对本发明进行修改或者等同替换,而不能脱离本发明的精神和范围的任何修改或者局部替换,其均应该涵盖在本发明的权利要求范围当中。
Claims (6)
1.一种基于空时功率倒置阵的宽带干扰抑制装置,其特征在于:包括空时二维功率倒置阵结构及重定时延时最小均方滤波器;所述空时二维功率倒置阵结构由M个阵元组成,设第一个阵元为参考阵元,其他M-1个阵元为权值阵元,每个权值阵元各有N个单位延时器,每个单位延时器都有1个复数权值,使得每个权值阵元构成一个阵元FIR滤波器;
重定时延时最小均方滤波器包括M-1组权值处理模块、两个单位延时器、第一加法器、一个单位延时器及第二加法器;两个单位延时器的输入端连接参考阵元接收信号,每个权值处理模块一路输入端一一对应连接权值阵元接收信号、每个权值处理模块另一路输入端一一对应连接一个初始信号,每组权值处理模块一路输出端连接三个单位延时器的输入端、三个单位延时器的输出端连接权值处理模块;所有权值处理模块另一路输出端均连接第二加法器的输入端,第二加法器的输出端和两个单位延时器的输出端均连接第一加法器的输入端,第一加法器的输出端连接一个单位延时器的输入端,一个单位延时器的输入端分别与每个权值处理模块的输入端相连;
抑制方法如下:
S1)初始化空时二维功率倒置阵模型,包括以下子步骤:
S11)空时二维功率倒置阵结构由M个阵元组成,设第一个阵元为参考阵元1,其他M-1个阵元为权值阵元,M个阵元的接收信号分别记为x0(n)、x1 (n)、……、xM-1(n),其中,x0 (n)表示参考阵元接收信号,x1(n) 、……、xM-1(n)为权值阵元接收信号,n表示当前时刻;
S12)每个权值阵元各有N个单位延时器,记为xi (n)、xi (n-1) 、……、xi (n-N+1),其中,i=0、……、M-1,每个单位延时器的复数权值记为wij,i=0、……、M-1,j=1、……、N,其中i表示阵元序号,j表示延时单元序号;
S13)初始化参考阵元x0路权值为w01=1,w0j=0,其中j=2、……、N,初始化权值路x1、……、xM-1的权值为
wij=0,其中i=1、……、M-1,j=1、……、N
记输入信号向量为
x(n)=[x1 (n)⋯x1 (n-N+1)⋯xM-1 (n)⋯xM-1(n-N+1)]T∈C(M-1)×N
记权向量为
w=[w1,1 w1,2⋯w1,N⋯wM-1,1 wM-1,2⋯wM-1,N ]T∈C(M-1)×N
初始化当前时刻为n=1;
S2)根据多径信道环境计算处理延时量,将参考阵元接收信号进行延时,将x0路信号x0(n)延时2D,得到x0 (n-2D),其中2D要大于最大多径传播时延;
S3)利用权值处理模块计算输出信号、误差信号以及权向量的更新,具体为:
S31)对于n=1,2,3,……,利用当前权值w,计算输出信号为
y(n)=w1,1 * x1 (n)+w1,2 * x1 (n-1)+⋯+w1,N * x1 (n-N+1)+⋯+wM-1,1 * xM-1 (n)+wM-1,2 * xM-1(n-1)+⋯+wM-1,N * xM-1 (n-N+1)=wH x(n)
计算期望延时信号为
d(n-2D) = -x0 (n-2D)
计算延时误差信号为
e(n-2D) = d(n-2D)-y(n-2D)
S32)设置步长μ,对于空时二维功率倒置阵,μ由移位器实现,μ为2-8~2-12,并进行权值的更新,权向量更新为:
w(n+1) = w(n)+μx(n-2D)e* (n-2D)
S4)令n=n+1,转子步骤S31),直到满足‖w(n+1)-w(n)‖2 2≤ε,其中ε为迭代终止阈值;
S5)为实现自适应宽带干扰抑制,对于当前时刻n,输出信号为y ̃(n)=x0(n)+wH x(n)。
2.根据权利要求1所述基于空时功率倒置阵的宽带干扰抑制装置,其特征在于:所述参考阵元的权值始终为1,M个阵元的接收信号分别记为x0(n)、x1 (n)、……、xM-1(n),其中,x0(n)表示参考阵元接收信号,x1(n) 、……、xM-1(n)为权值阵元接收信号,n表示当前时刻;在每个权值阵元后加上N个单位延时器Z-1,记为xi (n)、xi (n-1) 、……、xi (n-N+1),其中,i=0、……、M-1,每个单位延时器的复数权值记为wij,i=0、……、M-1,j=1、……、N,其中i表示阵元序号,j表示延时单元序号。
3.根据权利要求2所述基于空时功率倒置阵的宽带干扰抑制装置,其特征在于:每个所述权值阵元接收信号xi(n)输送至对应的权值处理模块后经过三个单位延时器形成延时信号xi(n-3D)给对应的权值处理模块输出;每个初始信号yi(n-1D)输送至对应的权值处理模块后形成延时信号yi(n-2D),所有的延时信号yi(n-2D)经过第二加法器产生权值合并输出信号y(n-2D),权值合并输出信号y(n-2D)和参考阵元接收信号x0(n)经过第一加法器形成误差信号e(n-2D),误差信号e(n-2D)经过一个单位延时器形成延时误差信号e(n-3D)输送至每个权值处理模块输出。
4.根据权利要求3所述基于空时功率倒置阵的宽带干扰抑制装置,其特征在于:每组所述权值处理模块包括四个乘加单元、两个单位延时器、第一加法器及第二加法器,四个乘加单元分别为第一乘加单元、第二乘加单元、第三乘加单元、第四乘加单元;第一乘加单元和第三乘加单元均包括上支路乘法器、下支路乘法器和加法器,上支路乘法器和下支路乘法器的输出端均与加法器的输入端相连;第二乘加单元和第四乘加单元均包括左支路乘法器、右支路乘法器和加法器,左支路乘法器和右支路乘法器的输出端均与加法器的输入端相连。
5.根据权利要求4所述基于空时功率倒置阵的宽带干扰抑制装置,其特征在于:每组所述权值处理模块的输入信号为IQ信号xiI(n)、xiQ(n),i=1、……、M-1,IQ信号xiI(n)、xiQ(n)一路经过三个单位延时器形成三个单位延时信号xiI(n-3D)、xiQ(n-3D),三个单位延时信号xiI(n-3D)、xiQ(n-3D)经过两个单位延时器形成五个单位延时信号xiI(n-5D)、xiQ(n-5D)输出;延时信号xiI(n-3D)分别连接第一乘加单元的上支路乘法器、第三乘加单元的上支路乘法器,延时信号xiQ(n-3D)分别连接第一乘加单元的下支路乘法器、第三乘加单元的下支路乘法器,第一乘加单元的加法器输出端依次经过一个单位延时器、上支路移位器、加法器、一个单位延时器与第二乘加单元的左支路乘法器和右支路乘法器输入端均相连,同理,第三乘加单元的加法器输出端依次经过一个单位延时器、下支路移位器、加法器、一个单位延时器与第四乘加单元的左支路乘法器和右支路乘法器输入端均相连;信号xiI(n)分别连接第三乘加单元的左支路乘法器和第四乘加单元的左支路乘法器,信号xiQ(n) 分别连接第三乘加单元的右支路乘法器和第四乘加单元的右支路乘法器;初始信号yiI(n-1D)输送至第一加法器的左支路,同时第二乘加单元加法器的输出端经过一个单位延时器输送至第一加法器的右支路,第一加法器左右支路相加后再经过一个单位延时器形成延时信号yiI (n-2D);初始信号yiQ(n-1D)输送至第二加法器左支路,同时第四乘加单元加法器的输出端经过一个单位延时器输送至第二加法器右支路,第二加法器左右支路相加后再经过一个单位延时器形成延时信号yiQ(n-2D)。
6.根据权利要求4所述基于空时功率倒置阵的宽带干扰抑制装置,其特征在于:每组所述权值处理模块的延时信号yiI(n-2D)经过第一加法器形成权值合并输出信号yI(n-2D),每个权值处理模块的延时信号yiQ(n-2D)经过第一加法器形成权值合并输出信号yQ(n-2D);权值合并输出信号yI(n-2D)和参考阵元接收信号x0I (n)经过第一加法器形成误差信号eI(n-2D),误差信号eI(n-2D)经过一个单位延时器形成延时误差信号eI (n-3D)输送至每个权值处理模块,延时误差信号eI(n-3D)一路分别连接第一乘加单元的上支路乘法器和第三乘加单元的上支持乘法器,延时误差信号eI(n-3D)另一路经过权值处理模块中的一个单位延时器形成延时误差信号eI(n-4D),延时误差信号eI(n-4D)输出;同理,权值合并输出信号yQ(n-2D)和参考阵元接收信号x0Q (n)经过第一加法器形成误差信号eQ(n-2D),误差信号eQ(n-2D)经过一个单位延时器形成延时误差信号eQ(n-3D)输送至每个权值处理模块,延时误差信号eQ(n-3D) 一路分别连接第一乘加单元的下支路乘法器和第三乘加单元的下支持乘法器,延时误差信号eQ(n-3D)另一路经过权值处理模块中的一个单位延时器形成延时误差信号eQ(n-4D)输出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111184393.2A CN113630355B (zh) | 2021-10-12 | 2021-10-12 | 基于空时功率倒置阵的宽带干扰抑制装置及抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111184393.2A CN113630355B (zh) | 2021-10-12 | 2021-10-12 | 基于空时功率倒置阵的宽带干扰抑制装置及抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113630355A CN113630355A (zh) | 2021-11-09 |
CN113630355B true CN113630355B (zh) | 2022-02-08 |
Family
ID=78391063
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111184393.2A Active CN113630355B (zh) | 2021-10-12 | 2021-10-12 | 基于空时功率倒置阵的宽带干扰抑制装置及抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113630355B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108462521A (zh) * | 2018-02-11 | 2018-08-28 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 自适应阵列天线的抗干扰实现方法 |
CN113472371A (zh) * | 2021-05-31 | 2021-10-01 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 自适应阵列天线数字波束合成抗干扰处理方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1989007864A1 (en) * | 1988-02-11 | 1989-08-24 | Rilling Kenneth F | Multipath reduction system |
JP2663820B2 (ja) * | 1992-12-28 | 1997-10-15 | 日本電気株式会社 | 判定帰還形等化器 |
US7714782B2 (en) * | 2004-01-13 | 2010-05-11 | Dennis Willard Davis | Phase arrays exploiting geometry phase and methods of creating such arrays |
CN103792550B (zh) * | 2014-02-11 | 2015-12-02 | 哈尔滨工程大学 | 一种基于阵列天线和gps/sins的联合抗干扰方法 |
CN110708103B (zh) * | 2019-10-30 | 2022-05-03 | 电子科技大学 | 一种无需预延迟的宽带波束形成方法 |
-
2021
- 2021-10-12 CN CN202111184393.2A patent/CN113630355B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108462521A (zh) * | 2018-02-11 | 2018-08-28 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 自适应阵列天线的抗干扰实现方法 |
CN113472371A (zh) * | 2021-05-31 | 2021-10-01 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 自适应阵列天线数字波束合成抗干扰处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113630355A (zh) | 2021-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2982504B2 (ja) | 適応受信機 | |
EP0604956B1 (en) | Decision feedback equalizer with adaptive filter array operating as feedforward filter of the equalizer | |
US8599906B2 (en) | Closed form calculation of temporal equalizer weights used in a repeater transmitter leakage cancellation system | |
EP0961416B1 (en) | Adaptive array transceiver | |
Erdem et al. | Monostatic antenna in-band full duplex radio: Performance limits and characterization | |
US10243593B2 (en) | Method of combatting interference by spatial filtering or spatio-temporal filtering in a multi-channel receiver | |
US11005540B2 (en) | Method and system for multiple input, multiple output communications in millimeter wave networks | |
EP2254270B1 (en) | Receiver and receiving method | |
CN111537958B (zh) | 一种宽线性降秩最小熵无失真响应的波束形成方法 | |
Ohgane | Characteristics of CMA adaptive array for selective fading compensation in digital land mobile radio communications | |
CN106941383B (zh) | 一种基于空时滤波的突发干扰抑制方法 | |
CN113630355B (zh) | 基于空时功率倒置阵的宽带干扰抑制装置及抑制方法 | |
CN110365421B (zh) | 一种单载波干扰抑制的多输入多输出水声通信方法 | |
CN113824488B (zh) | 基于判决反馈自适应对消的卫星通信非恶意干扰抑制方法 | |
JPH0786972A (ja) | 適応等化器 | |
Ratnam et al. | Multi-antenna FSR receivers: Low complexity, non-coherent, massive antenna receivers | |
Walkenhorst et al. | Polarization-based interference mitigation for OFDM signals in channels with polarization mode dispersion | |
WO2005002067A2 (en) | Multi-carrier spread spectrum using non-linear modification of sub-carrier bands | |
Djigan | Adaptive antenna array, shared with adaptive equalizer | |
Martone | On-board regeneration of uplink signals using a blind multichannel estimator | |
Liang et al. | A time-domain calibration scheme of channel reciprocity for TDD MIMO-OFDM system with IQ imbalance | |
Ikeda et al. | Interference Suppression Performance of CMA Using Time-Domain Inter-Symbol Spreading in Time-Varying Channel Environment | |
Kim et al. | Space-time-frequency processing for OFDM wireless mobile communications | |
Shayovitz et al. | Robust Low Complexity Digital Self Interference Cancellation for Multi Channel Full Duplex Systems | |
CN117978605A (zh) | 一种lte系统复合域自适应零陷抗干扰方法 |
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 |