CN110266363B - 一种基于张量的分布式扩散自适应抗干扰方法 - Google Patents
一种基于张量的分布式扩散自适应抗干扰方法 Download PDFInfo
- Publication number
- CN110266363B CN110266363B CN201910561315.6A CN201910561315A CN110266363B CN 110266363 B CN110266363 B CN 110266363B CN 201910561315 A CN201910561315 A CN 201910561315A CN 110266363 B CN110266363 B CN 110266363B
- Authority
- CN
- China
- Prior art keywords
- array
- tensor
- weight vector
- vector
- time
- 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/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
-
- 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/086—Weighted combining using weights depending on external parameters, e.g. direction of arrival [DOA], predetermined weights or beamforming
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Radio Transmission System (AREA)
Abstract
本发明属于分布式波束形成领域,具体为一种基于张量的分布式扩散自适应抗干扰方法,用以解决阵元数较大时传统分布式阵列自适应抗干扰方法复杂度变大、收敛速度变慢、实时性降低的问题。本发明将高维的全局多线性问题转化为若干个低维线性问题,具体表现为通过基于分集的张量模型,使本发明方法在子阵水平上并行处理,在自适应过程达到平稳后,选取任一节点的稳态权向量作为最终权向量,并利用该权张量对接收信号进行滤波。和传统的波束协调算法相比,本发明方法具有更快的收敛速度,更低的计算复杂度,从而具有更好的实时性。此外节点间共享的数据是所有子阵的回归矢量,而非原始的阵列接收信号,因此减少了共享的数据总量,提高了节点通信效率。
Description
技术领域
本发明属于分布式波束形成领域,主要涉及分布式自适应策略和多线性协同滤波,具体为一种基于张量的分布式扩散自适应抗干扰方法。
背景技术
阵列抗干扰应用领域由来已久,多年来已经形成很多成熟的理论。基于最小均方误差(MMSE)准则的LMS算法及其拓展而来的归一化LMS算法、变步长LMS算法,基于最小二乘准则(LS)的RLS算法等,此类算法适用于量测信号(期望信号)易获取的应用场景;基于最小方差无失真响应(MVDR)准则和线性约束最小方差(LCMV)准则的数字波束形成算法适用于部分信号或干扰方向已知的干扰抑制场景,相对应的还有基于最大信干噪比(MSINR)准则的波束形成算法。然而,在信号处理算法日趋完善的今天,已有的较为成熟的模型和理论在某些应用背景下需要进一步发展,科研工作者一直致力于在各个方面优化理论,完善方法,做到抗干扰算法高成效,低复杂度,低时耗,低成本等。近年来,基于张量的阵列信号处理理论日臻完善,也越来越多地应用到阵列抗干扰中,2016年,Lucas N.Ribeiro等人在文献《Tensor Beamforming for multilinear translation invariant arrays》中提出了一种基于MMSE准则的张量波束形成算法,有效地降低数据计算的复杂度,从而耗时小、效率高。张量信号处理的优势在于很好的将低维数据和高维数据进行分解或合并转化,从而降低问题的复杂程度或者提高算法的精度,使之符合人们的期望。
下面给出一些必要的张量运算法则:
相对于单阵列抗干扰研究,分布式阵列网络越来越受到科研工作者的青睐,分布式阵列抗干扰算法应运而生。自从2006以来,Ali H.Sayed等人对分布式自适应算法进行了大量的深入的研究;在2012年,他在文献《Beam coordination via diffusion adaptationover array network》中将分布式自适应策略应用到阵列抗干扰中,提出了用自适应和结合(ATC)算法求解出最优权向量,并利用该权向量对任一阵列的接收信号进行滤波以达到保留期望信号并抑制干扰的目的。
下面给出波束协调算法:
考虑一个包含N个节点的网络,各节点包含完全相同的天线阵列,其中各阵列阵元数为Ms;假设有一个期望的远场复窄带信号入射到天线阵列网络上,同时受到P-1个复窄带信号的干扰,则阵列k的离散复基带接收信号表示为:
其中,k=1,2,...,N,sp(t)是窄带信号的离散基带形式,而接收噪声zk(t)是方差为的零均值加性高斯白噪声向量,与不同阵列n的接收噪声zn(t)或不同时刻t1的接收噪声zk(t1)均相互独立,且独立于阵列接收信号uk(t)。假设抗干扰系统最优权向量为wo,那么线性量测信号模型如下:
发明内容
本发明的目的在于提供一种基于张量的分布式扩散自适应抗干扰方法,本发明将长矢量权估计问题转化为若干较小的线性估计问题,旨在解决阵元数较大时传统分布式阵列自适应抗干扰方法复杂度变大、收敛速度变慢、实时性降低的问题。
为实现上述目的,本发明采用的技术方案如下:
一种基于张量的分布式扩散自适应抗干扰方法,在任意阵列k上,k=1,2,...,N、N为网络中阵列总数,包括以下步骤:
步骤2.计算阵列k的子阵回归矢量:
对于阵列k的第l子阵,计算当前t时刻子阵回归矢量vk,l(t)、l=1,2,3:
其中,wk,l(t-1)为t-1时刻阵列k的第l子阵的权向量估计;
步骤3.计算阵列k的归一化因子βk(t):
步骤4.迭代更新阵列k的权向量的中间估计ψk,l(t):
其中,μk为迭代步长:0<μk<2,ε为预设常数:ε<10-6;
步骤5.扩散更新阵列k的权向量估计:
其中,系数an,k为结合矩阵A的第(n,k)个元素;
本发明的有益效果在于:
本发明提出的基于张量的分布式扩散自适应抗干扰方法,具有如下优点:
1.本发明将高维的全局多线性估计问题转化为若干个低维线性估计问题,为此类全局多线性问题提供了一种有效的的自适应求解方案。
2.相比已有的分布式抗干扰方法,本发明引入基于分集的张量模型,进行邻域多线性协同滤波,显著提高了算法收敛速度,使迭代过程仅需较少的样本就能收敛到稳态。
3.本发明可改善节点通信效率。由于节点间共享的数据不再是原始的阵列接收信号,而是所有子阵的回归矢量,因此减少了共享的数据总量,提高了节点通信效率。
4.本发明考虑了极化分集,可以有效地抑制与期望信号同方向的干扰。
5.本发明能够自适应地并行迭代所有子阵权向量估计,大大提高计算效率。因为子阵回归矢量和上一时刻的子阵权向量估计以及当前时刻阵列接收信号张量有关,所以本发明可并行求解所有子阵回归矢量,且能够并行迭代所有阵子权向量估计。
6.本发明的方法不局限于当前的信号模型,可拓展为更高维张量模型的波束形成算法。
附图说明
图1为本发明基于张量的分布式扩散自适应抗干扰方法各节点步骤流程图。
图2为本发明实施例中分布式阵列网络结构图。
图3为本发明实施例中节点网络拓扑。
图4为本发明实施例中各节点的噪声水平分布图。
图5、6为本发明实施例中本发明方法和波束协调方法MSE和SINR学习曲线对比图。
图7为本发明实施例中当期望信号方向突变时本发明方法和波束协调方法的自适应情况。
具体实施方式
下面结合附图和实施例对本发明做进一步详细说明。
本发明一种基于张量的分布式扩散自适应抗干扰方法,其分布式阵列网络如图2所示,节点网络拓扑如图3所示;假设极化空域滤波器wo张成的张量滤波器是秩1的,即:
本实施例提出一种基于张量的分布式扩散自适应抗干扰方法,流程如图1所示,在任意阵列k上的具体过程如下:
步骤1.获取阵列k接收信号张量和量测信号
设定包含N个节点的网络,各节点包含完全相同的天线阵列,其中各阵列阵元数为Ms,每个阵元由3个相互正交的电偶极子和3个相互正交的磁偶极子构成,各阵列以入射窄带信号的半波长间隔布阵;假设有一个期望的远场复窄带信号入射到天线阵列网络上,同时受到P-1个复窄带信号的干扰,则阵列k、k=1,2,...,N的离散复基带接收信号张量表示为:
其中,θp,γp和ηp分别表示第p个入射信号的方位角、极化相角和极化相位差,sp(t)为复窄带信号的离散基带形式,是方差为的零均值加性高斯白噪声张量,于不同阵列k的接收噪声或不同时刻t1的接收噪声相互独立、且独立于阵列接收信号导向张量为:
其中,φp=πsinθp,Ms=M1M2;
获取阵列k的量测信号dk(t),k=1,2,...,N:
给出分布式阵列抗干扰方法的目标函数如下:
步骤2.计算阵列k的子阵回归矢量
对于阵列k的第l子阵,计算当前t时刻子阵回归矢量vk,l(t)、l=1,2,3作为阵列k的第l子阵的接收信号;子阵回归矢量vk,l(t)由当前t时刻的阵列接收信号张量与上一时刻估计的子阵权向量集{wk,q(t-1)}q≠l的正序多线性积给出:
其中,·H表示共轭转置操作;
步骤3.分两步计算阵列k的归一化因子βk(t)
其中,ρk(t)代表t时刻阵列k的能量估计,表示阵列k的包括自身的邻域集;系数cn,k是交换矩阵C的第(n,k)个元素,给出邻居阵列n在归一化因子βk(t)中所占比重;归一化因子βk(t)表示的是阵列k及其邻域阵列在t时刻的能量加权平均;交换矩阵C满足的约束为:cn,k=0,1TC=1T,C1=1,即,交换矩阵C的列元素和为1、行元素和也为1;
步骤4.迭代更新阵列k、各子阵l的权向量的中间估计:
其中,ψk,l(t)代表t时刻阵列k第l子阵权向量的中间估计,wk,l(t-1)代表t-1时刻阵列k第l子阵权向量估计,l=1,2,3,各子阵并行地按上式更新;μk是迭代步长,控制迭代收敛的速度和稳态;ε是一个非常小的正数,用于保证分母不为0;系数cn,k是交换矩阵C的第(n,k)个元素,这里表示阵列k邻居节点n的数据在此迭代更新中的贡献比例;
步骤5.扩散更新阵列k、各子阵l的权向量估计:
其中,阵列k的所有子阵并行地按上式实现权向量估计,l=1,2,3;系数an,k是结合矩阵A的第(n,k)个元素,表示阵列k邻居阵列n的权向量中间估计ψn,l(t)在此扩散更新中的贡献量;结合矩阵A满足的约束为:an,k=0,1TA=1T,即,结合矩阵A列元素和为1;
因此,本实施例的具体实现步骤为:
本发明方法属于分布式算法,各阵列实现步骤相同,因此只给出任意阵列k上的具体实现步骤:
步骤1.相关参数和权向量初始化
初始化阵列k、各子阵l的权向量wk,l(0)为任意不为零的复向量;所有阵列k的步长相同,满足0<μk<2;给定小的正数ε<10-6;按照相应的计算规则给出交换矩阵C和结合矩阵A,rk代表阵列k的邻域包含的阵列个数,也叫阵列k的度,具体计算规则如下:
步骤3.计算阵列k的子阵回归矢量
步骤4.计算阵列k的归一化因子βk(t)
步骤5.迭代更新阵列k、各子阵l的权向量的中间估计ψk,l(t)
步骤6.扩散更新阵列k、各子阵l的权向量估计wk,l(t)
步骤7.计算平稳过程最优权值,得到输出信号。
下面通过仿真实验对比本发明方法和波束协调算法的抗干扰效果,说明本发明的可行性,优越性:
仿真实验
仿真1:20个节点互联的网络,每个节点是36个电磁矢量传感器构成的阵列,其中权分量长度M1=M2=M3=6,各节点噪声分布如图4所示,考虑一个期望单音信号,方向为30度,极化相角为20,极化相位差为-50,功率为0dB,频率为1e3,两个单音干扰信号,方位角,极化相角,极化相位差分别为(30,70,-50),(-60,20,-50),功率分别为15dB,15dB,频率分别为1.5kH z,2kH z,本发明方法两种步长分别为μk=0.65,k=1,2,..,N,μk=0.7,k=1,2,..,N,波束协调算法步长为μk=1.2,k=1,2,..,N,各节点阵列接收噪声功率均为0dB,采样率8kHz,快拍数400,500次独立重复实验,实验结果如图5、6所示。
如图5所示,均方误差随着快拍数增加逐渐下降,并最终达到稳态水平,观察学习曲线发现,在稳态水平几近重合的情形下,本发明方法收敛速度远远高于波束协调算法,实现小快拍数收敛;同理,图6所示信干噪比稳态值接近时,本发明方法收敛速度远远高于波束协调算法。由于本发明算法在单个节点上引入并行处理,很大程度上节省计算时间,加之本发明算法更小快拍数收敛,因此可大大提高算法的时效性,节省算法收敛所需的时间。
仿真2:20个节点互联的网络,每个节点是36个电磁矢量传感器构成的阵列,其中权分量长度M1=M2=M3=6,各节点量测噪声分布如图4所示,考虑一个期望单音信号,方向为30度,极化相角为20,极化相位差为-50,功率为0dB,频率为1e3,两个单音干扰信号,方位角,极化相角,极化相位差分别为(30,70,-50),(-60,20,-50),功率分别为15dB,20dB,频率分别为1.5kHz,2kHz,本发明方法两种步长分别为μk=0.65,k=1,2,..,N,μk=0.7,k=1,2,..,N,波束协调算法步长为μk=1.2,k=1,2,..,N,各节点阵列接收噪声功率均为0dB,采样率8kHz,快拍数1400,从第700快拍开始,期望信号方向发生1度偏移,500次独立重复实验,实验结果如图7所示。
如图7所示,在方向突变前,本发明方法收敛到稳态值速度远远高于波束协调算法,在第701快拍处,期望信号突变偏移原方向1度,此时算法需要重新适应到稳态值,可以看出本发明算法仍然可以更快的收敛,且稳态值没有损失,而波束协调算法收敛速度比突变前更缓慢,由此可以看出本发明方法在应对期望信号方向突变情况方面具有更好的适应能力,可应用于实时波束形成。
以上所述,仅为本发明的具体实施方式,本说明书中所公开的任一特征,除非特别叙述,均可被其他等效或具有类似目的的替代特征加以替换;所公开的所有特征、或所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以任何方式组合。
Claims (1)
1.一种基于张量的分布式扩散自适应抗干扰方法,在任意阵列k上,k=1,2,...,N、N为网络中阵列总数,包括以下步骤:
步骤2.计算阵列k的子阵回归矢量:
对于阵列k的第l子阵,计算当前t时刻子阵回归矢量vk,l(t)、l=1,2,3:
其中,wk,l(t-1)为t-1时刻阵列k的第l子阵的权向量估计;
步骤3.计算阵列k的归一化因子βk(t):
步骤4.迭代更新阵列k的权向量的中间估计ψk,l(t):
其中,μk为迭代步长:0<μk<2,e为预设常数:e<10-6,系数cn,k为交换矩阵C的第(n,k)个元素;
步骤5.扩散更新阵列k的权向量估计:
其中,系数an,k为结合矩阵A的第(n,k)个元素;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910561315.6A CN110266363B (zh) | 2019-06-26 | 2019-06-26 | 一种基于张量的分布式扩散自适应抗干扰方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910561315.6A CN110266363B (zh) | 2019-06-26 | 2019-06-26 | 一种基于张量的分布式扩散自适应抗干扰方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110266363A CN110266363A (zh) | 2019-09-20 |
CN110266363B true CN110266363B (zh) | 2021-09-24 |
Family
ID=67921795
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910561315.6A Active CN110266363B (zh) | 2019-06-26 | 2019-06-26 | 一种基于张量的分布式扩散自适应抗干扰方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110266363B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111224704B (zh) * | 2019-11-12 | 2022-10-11 | 电子科技大学 | 一种分布式自适应降秩波束形成方法 |
CN111262739B (zh) * | 2020-01-17 | 2021-04-06 | 西南大学 | 事件触发通信下的分布式自适应局部扩散控制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104868946A (zh) * | 2015-06-12 | 2015-08-26 | 哈尔滨工业大学 | 自适应加权的子阵级混合mimo-相控阵系统的干扰抑制方法 |
CN106779161A (zh) * | 2016-11-15 | 2017-05-31 | 南京航空航天大学 | 一种气动外形约束下分布式三维阵的阵列布局优化方法 |
CN107045131A (zh) * | 2017-05-24 | 2017-08-15 | 西北工业大学 | 张量分解卫星导航抗干扰方法 |
CN108462521A (zh) * | 2018-02-11 | 2018-08-28 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 自适应阵列天线的抗干扰实现方法 |
-
2019
- 2019-06-26 CN CN201910561315.6A patent/CN110266363B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104868946A (zh) * | 2015-06-12 | 2015-08-26 | 哈尔滨工业大学 | 自适应加权的子阵级混合mimo-相控阵系统的干扰抑制方法 |
CN106779161A (zh) * | 2016-11-15 | 2017-05-31 | 南京航空航天大学 | 一种气动外形约束下分布式三维阵的阵列布局优化方法 |
CN107045131A (zh) * | 2017-05-24 | 2017-08-15 | 西北工业大学 | 张量分解卫星导航抗干扰方法 |
CN108462521A (zh) * | 2018-02-11 | 2018-08-28 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 自适应阵列天线的抗干扰实现方法 |
Non-Patent Citations (2)
Title |
---|
"A Robust GNSS Polarized Space-Time Anti-Interference Method Based on Null Broadening";Guoqing Xia等;《 2018 10th International Conference on Communications, Circuits and Systems (ICCCAS)》;20181224;全文 * |
"基于最优滤波器的强混沌背景中谐波信号检测方法研究";胡进峰等著;《物理学报》;20150929;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110266363A (zh) | 2019-09-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JPH06204902A (ja) | 判定帰還形等化器 | |
CN110266363B (zh) | 一种基于张量的分布式扩散自适应抗干扰方法 | |
Varade et al. | Robust algorithms for DOA estimation and adaptive beamforming for smart antenna application | |
Xie et al. | Dictionary learning for channel estimation in hybrid frequency-selective mmWave MIMO systems | |
CN113162665A (zh) | 一种基于深度学习信道预测的预编码方法 | |
CN111537958A (zh) | 一种宽线性降秩最小熵无失真响应的波束形成方法 | |
Sohal et al. | Deep learning based analog beamforming design for millimetre wave massive MIMO system | |
CN111551892A (zh) | 一种稳健自适应波束形成方法及装置 | |
Alinezhad et al. | Broadband adaptive beamforming of conformal arrays for wireless communications based on generalized sidelobe canceller | |
CN115426013A (zh) | 基于无监督学习的ris辅助通信系统预编码与相移方法 | |
Lemma et al. | Analysis of ESPRIT based joint angle-frequency estimation | |
CN116148777A (zh) | 一种基于阵列虚拟扩展的阵列雷达主瓣干扰抑制方法 | |
KR100941699B1 (ko) | 진폭조정 적응 빔 형성 장치 | |
CN112165347B (zh) | 一种毫米波纯模拟波束成形用户调度方法 | |
Shubair | Improved smart antenna design using displaced sensor array configuration | |
JPH11308130A (ja) | 干渉波抑圧装置 | |
CN107064884B (zh) | 基于规则重叠子阵的自适应波束形成方法 | |
Thapa et al. | Performances of RLS algorithm for smart antennas in mobile communication system | |
Bouacha et al. | Modified Blind Beamforming Algorithm for Smart Antenna System. | |
Zhang et al. | Robust beamforming for coherent signals in the presence of steering vector and correlation matrix errors | |
CN112165346B (zh) | 一种毫米波通信共享幅度加权模拟波束赋形优化方法 | |
SHIRVANI et al. | A new switched-beam setup for adaptive antenna array beamforming | |
Shubair et al. | Improved adaptive beamforming using the least mean mixed-norm algorithm | |
CN112711837A (zh) | 一种低快拍下抗强干扰的波束形成方法 | |
Biedka | A method for reducing computations in cyclostationarity-exploiting beamforming |
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 |