CN117997360A - 一种适用于大动态范围期望信号的阵列抗干扰方法 - Google Patents
一种适用于大动态范围期望信号的阵列抗干扰方法 Download PDFInfo
- Publication number
- CN117997360A CN117997360A CN202311805871.6A CN202311805871A CN117997360A CN 117997360 A CN117997360 A CN 117997360A CN 202311805871 A CN202311805871 A CN 202311805871A CN 117997360 A CN117997360 A CN 117997360A
- Authority
- CN
- China
- Prior art keywords
- interference
- signal
- frequency domain
- array
- sub
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 239000011159 matrix material Substances 0.000 claims abstract description 25
- 238000001228 spectrum Methods 0.000 claims abstract description 24
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 9
- 238000006243 chemical reaction Methods 0.000 claims abstract description 6
- 239000013598 vector Substances 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000001934 delay Effects 0.000 claims description 3
- 230000002452 interceptive effect Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000003044 adaptive effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 230000007480 spreading Effects 0.000 description 3
- 230000001629 suppression Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Landscapes
- Noise Elimination (AREA)
Abstract
本发明公开了一种适用于大动态范围期望信号的阵列抗干扰方法,属于阵列抗干扰技术领域,包括以下步骤:S1、对接收信号进行FFT变换,得到分解的频域子带;S2、基于期望信号的频谱特征,挑选用于干扰估计的子带,并估计干扰协方差矩阵;S3、对干扰协方差矩阵进行特征分解,估计干扰信源数目和干扰信号DOA;S4、在所有子带上构造干扰空间的正交空间,进行正交投影以消除干扰。本发明采用上述的一种适用于大动态范围期望信号的阵列抗干扰方法,在有效抑制干扰信号的同时,避免将期望信号误认为干扰,可有效降低对期望信号的损失。
Description
技术领域
本发明涉及阵列抗干扰技术领域,尤其是涉及一种适用于大动态范围期望信号的阵列抗干扰方法。
背景技术
自适应阵列处理可以有效抑制某个方向上的干扰信号,已经被广泛应用于导航接收机的干扰抑制中。然而,当期望信号的动态范围较大时,会出现期望信号功率高于噪声功率的情况。
针对干扰信号功率和期望信号功率都高于噪声功率的场景,基于功率倒置的阵列抗干扰方法会将期望信号作为干扰进行抑制,从而对期望信号造成损失。要保证期望信号增益的同时抑制空间上其他来向的干扰,自适应波束成形(Adaptive DigitalBeamforming,ADBF)方法是一个选择。但是该方法需要已知期望信号到达方向(DirectionofArrival,DOA),然而期望信号的DOA先验知识通常难以获得。
发明内容
本发明的目的是提供一种适用于大动态范围期望信号的阵列抗干扰方法,基于期望信号是频谱特征已知的扩频信号的先验信息,该方法利用期望信号的频谱特征,结合空频自适应处理的思想,将接收信号通过傅里叶变换分解为频域子带,通过挑选“零点”位置对应的子带进行干扰估计,在有效抑制干扰信号的同时,避免将期望信号误认为干扰,可有效降低对期望信号的损失。
为实现上述目的,本发明提供了一种适用于大动态范围期望信号的阵列抗干扰方法,包括以下步骤:
S1、对接收信号进行FFT变换,得到分解的频域子带;
S2、基于期望信号的频谱特征,挑选用于干扰估计的子带,并估计协方差矩阵;
期望信号在采样率带宽内存在主瓣和旁瓣,其中,主瓣幅度高,旁瓣幅度低,主瓣和旁瓣的连接处为频域零点,期望信号能量95%以上集中在信号通带内,频域零点附近的子带信号功率低,仅包含干扰;
选择频域零点位置对应的子带进行干扰估计,避免将期望信号误认为干扰,降低对期望信号的影响;
S3、对协方差矩阵进行特征分解,估计干扰信源数目和干扰信号DOA;
S4、在所有子带上构造干扰空间的正交空间,进行正交投影以消除干扰。
优选的,步骤S1中,假设第m个阵元在t时刻接收到的信号为:
其中,s0为阵列接收的期望信号,s1…sJ为J个干扰信号;τm,j为接收来向的信号时第m个阵元上信号j相对于参考阵元的时间延迟;n为高斯白噪声;m的范围为1~M,M为阵元数目;
对第m个阵元在t时刻接收到的信号做FFT变换,得到:
其中,fn表示频率;
将各阵元接收信号的FFT结果合并写成矩阵的形式,得到:
其中,
和/>分别是xm(t)、sj(t)和nm(t)的傅里叶变换后的结果;A(fn)中的元素是时延的频域表示;L为频域快拍数目;T表示转置。
优选的,步骤S2中,估计协方差矩阵包括:
阵列接收信号在频率fn处的协方差矩阵为:
其中,E表示均值,H表示共轭转置;
利用L个频域快拍计算得到协方差矩阵为:
将fn对应的索引记为NF,NF=1,2,…,NS,NS是分解的子带数;第NF个子带的协方差矩阵表示为:
其中,是信号的导向矢量,/>是/>信号的自相关结果,是频域噪声信号的功率,I是M维单位阵。
因此,本发明采用上述一种适用于大动态范围期望信号的阵列抗干扰方法,其技术效果如下:
(1)克服基于Power Inverse准则的抗干扰算法导致高信噪比期望信号损失的缺点以及自适应波束成形方法需要已知期望信号到达方向的限制,在保证期望信号增益的同时可抑制空间上其他来向的干扰。
(2)本发明所提供的方法在现有空频抗干扰技术的基础上,经过有限的改进即可实现大动态范围期望信号场景下的干扰信号抑制。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1为扩频信号的频域零点位置示意图;
图2为选择不同子带情况下的特征值分布;
图3是选择不同子带情况下期望信号特征值随输入SNR的变化;
图4是挑选零点位置数据处理相邻谱峰示意图;
图5是处理前各信号分量的频域FFT结果;
图6是本发明方法处理后各信号分量的频域FFT结果。
具体实施方式
以下通过附图和实施例对本发明的技术方案作进一步说明。
除非另外定义,本发明使用的技术术语或者科学术语应当为本发明所属领域内具有一般技能的人士所理解的通常意义。
实施例一
如图1所示,由于期望信号是直接序列扩频信号(后文简称扩频信号)且扩频信号的频谱特征是已知信息。根据扩频信号的频谱特征,可以发现扩频信号在采样率带宽内是存在主瓣和旁瓣的,主瓣幅度较高,旁瓣幅度较低,主瓣和旁瓣的连接处是频域零点。考虑到期望信号能量95%以上集中在信号通带内,信号频域零点附近的子带中信号功率低,可近似认为仅包含干扰。我们将这些位置定义为“零点”位置,选择这些频域“零点”位置对应的子带进行干扰估计,可以避免将信号误认为干扰,从而降低对期望信号的影响。
假设第m个阵元在t时刻接收到的信号为:
其中,s0为阵列接收的期望信号,s1…sJ为J个干扰信号;τm,j为接收来向的信号时第m个阵元上该信号相对于参考阵元的时间延迟;n为高斯白噪声。
对第m个阵元在t时刻接收到的信号做FFT变换,得到:
将各阵元接收信号的FFT结果合并写成矩阵的形式,得到:
其中
和/>分别是xm(t)、sj(t)和nm(t)的傅里叶变换后的结果;A(fn)中的元素是时延的频域表示。
阵列接收信号在频率fn处的协方差矩阵为:
利用L个频域快拍计算得到采样协方差矩阵为:
fn对应的索引记为NF,NF=1,2,…,NS,第NF个子带的协方差矩阵表示为:
其中,是信号的导向矢量,/>是/>信号的自相关结果,是频域噪声信号的功率,I是M维单位阵。经过FFT变换后的高斯白噪声在频域上仍然服从均匀分布的形式,频域幅度变成时域幅度的Ns倍。
对R(NF)进行特征分解,分解得到的特征值按降序排列得λ1≥λ2≥…≥λM,其对应的特征向量为v1,v2,…,vM,有:
挑选两个子带N1和N2进行特征分解并对比分析其特征值的差异,N1取零点位置的子带,N2取扩频信号通带内的子带。假设入射信号包括两个干扰信号和一个期望信号,干扰功率大于信号功率,输入SNR为14dB。图2为选择不同子带情况下的特征值分布,可以看出,前两个大特征值来自两个大功率干扰信号,第三个特征值来自期望信号。选取子带N2进行特征分解的结果中明显存在期望信号对应的特征值,那么对应的特征向量也会被划分至干扰空间,在后续的干扰消除中会对期望信号进行衰减;选取子带N1进行特征分解相比选取子带N2进行特征分解,期望信号分量减小,其期望信号特征值接近噪声特征值。在相同的期望信号条件下,显然选择子带N1进行干扰估计会得到比选择子带N2更好的结果。
针对图2进一步分析,设置输入SNR在0dB至25dB范围内变化,两个干扰信号不变。由于干扰功率远大于期望信号功率,所以SNR在该范围内变化时,λ3始终代表期望信号的特征值。选择不同子带情况下期望信号对应的特征值λ3随输入SNR的变化如图3所示,可以看出,选择N2子带时期望信号对应的特征值随输入SNR的增加而增大,选择N1子带时上期望信号对应的特征值随输入SNR的变化不大,且一直小于选择N2子带使期望信号对应的特征值。
在高SNR情况下,要使R(NF)中期望信号的分量最小,需要选择中/>幅度最小的子带进行干扰估计。通过观察扩频信号频谱得知,这些子带就是扩频零点位置子带。在数值上扩频信号在单边谱中第一零点频率等于扩频码周期的倒数1/Tc,即扩频码速率rc。假设一共存在P个零点,可以记零点频率为fZP(p),有
fZP(p)=prc p=1,2,…,P (12)
包含了频域零点的子带记作NZP(p),因为子带需要考虑双边谱,所以NZP(p)中p=±1,±2,…,±P。
在采样带宽内,只要干扰带宽覆盖到扩频零点位置子带,用该子带计算的协方差矩阵就能包含干扰分量,从而正确估计干扰信号信息。使用多个扩频零点子带估计会更好地捕捉到干扰信号,从而抑制干扰。因此,可以使用接收机采样率范围内的多个频域零点子带构造协方差矩阵,进行信源估计和DOA估计,并用每个零点位置估计的结果处理其相邻谱峰的数据。挑选零点位置数据处理相邻谱峰示意图如图4所示,坐标轴横轴表示数字中频信号频率,在每个零点处都进行干扰估计,根据得到的信息在附近两个峰之间的子带上进行干扰对消,这样就能更好地处理采样带宽内的多频点干扰。
此处只使用一个子带进行干扰信号估计。记该子带为NZP,记对应的协方差矩阵为RZP,有:
使用L个快拍计算得到的协方差矩阵的估计为:
对进行特征分解后划分特征空间:
较大的D个特征值及对应的特征向量组成ΛSI和USI,ΛSI=diag(λ1,λ2,…,λD),USI=[v1,v2,…,vD],表示信号和干扰空间;较小的M-D个特征值及对应的特征向量组成ΛN和UN,ΛN=diag(λD+1,…,λM),UN=[vD+1,…,vM],表示噪声空间。确定D的值就是划分信号子空间和噪声子空间的过程,常用的划分准则包括Akaike信息论准则或最小描述长度准则,这一步骤被称为信源数估计。
应用最小描述长度准则进行信源数估计,求出使下式值最小的D,D就是信源数:
接下来使用MUSIC算法进行二维DOA估计,MUSIC谱表示导向矢量与噪声空间的正交性,其表达式为:
其中是搜索导向矢量。当/>时与UN正交性最大,得到极大值,从而得到估计信号来向/>然而,MUSIC需要进行全维搜索,运算量过大。Root-MUSIC是一种使用多项式求根替代搜索过程的方法,相比MUSIC可以减少一半以上的运算量。Root-MUSIC的求根公式为:
为了求根计算方便,将其中的待求导向矢量根据沿X轴方向和Y轴方向重写为/>的形式,有:
把和/>中的/>项设为v,/>项设为u,在X和Y轴方向上等效的阵元数分别为MX和MY,则这两个方向的导向矢量分别表示为:
设z1=ej2πdu/λ,z2=ej2πdv/λ,简化求根公式得到:
求解上述方程,得到的估计分别为:
由于估计的和/>顺序并不符合原本的角度关系,还需要进行匹配,构造代价函数Θ为:
其中,需要将D个和/>依次代入aX和aY,共D2次代价函数的计算,取前D个最小值的组合,即获得了原本的/>和/>匹配关系。得到的DOA结果:
通过二维DOA估计得到D个信号的入射方向为了抑制干扰,在所有子带上根据DOA估计结果来构造子带频率对应的导向矢量/>然后取需要抑制的D个方向的构建/>有:
根据特征空间性质,干扰子空间的正交补空间的投影矩阵表示为:
然后对进行投影处理,有:
对Yf进行IFFT处理后还原为干扰抑制后的时域信号再进行后续信号处理。
下面对本发明所提方法的有效性进行matlab仿真验证,阵列模型采用半径为半波长的4阵元圆阵。期望信号的载波频率为1575.42MHz,码速率为10.23Mcps,入射方位角和俯仰角均为10°,信噪比为20dB。设置单干扰入射场景,干扰信号为中心频率1575.42MHz的随机调相宽带干扰,入射方向θ=25°,干噪比为100dB。图5给出了处理前各信号分量的频域FFT结果。图6给出了经本方法处理后各信号分量的频域FFT结果。
可以看出,经过抗干扰处理后,干扰信号抑制得到有效抑制同时,期望信号仍然保持了较高的信噪比。
因此,本发明采用上述一种适用于大动态范围期望信号的阵列抗干扰方法,基于期望信号是频谱特征已知的扩频信号的先验信息,该方法利用期望信号的频谱特征,结合空频自适应处理的思想,将接收信号通过傅里叶变换分解为频域子带,通过挑选“零点”位置对应的子带进行干扰估计,在有效抑制干扰信号的同时,避免将信号误认为干扰,可有效降低对期望信号的损失。
最后应说明的是:以上实施例仅用以说明本发明的技术方案而非对其进行限制,尽管参照较佳实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对本发明的技术方案进行修改或者等同替换,而这些修改或者等同替换亦不能使修改后的技术方案脱离本发明技术方案的精神和范围。
Claims (3)
1.一种适用于大动态范围期望信号的阵列抗干扰方法,其特征在于,包括以下步骤:
S1、对接收信号进行FFT变换,得到分解的频域子带;
S2、基于期望信号的频谱特征,挑选用于干扰估计的子带,并估计协方差矩阵;
期望信号在采样率带宽内存在主瓣和旁瓣,其中,主瓣幅度高,旁瓣幅度低,主瓣和旁瓣的连接处为频域零点,期望信号能量95%以上集中在信号通带内,频域零点附近的子带信号功率低,仅包含干扰;
选择频域零点位置对应的子带进行干扰估计,避免将期望信号误认为干扰,降低对期望信号的影响;
S3、对协方差矩阵进行特征分解,估计干扰信源数目和干扰信号DOA;
S4、在所有子带上构造干扰空间的正交空间,进行正交投影以消除干扰。
2.根据权利要求1所述的一种适用于大动态范围期望信号的阵列抗干扰方法,其特征在于,步骤S1中,假设第m个阵元在t时刻接收到的信号为:
其中,s0为阵列接收的期望信号,s1…sJ为J个干扰信号;τm,j为接收来向的信号时第m个阵元上信号j相对于参考阵元的时间延迟;n为高斯白噪声;m的范围为1~M,M为阵元数目;
对第m个阵元在t时刻接收到的信号做FFT变换,得到:
其中,fn表示频率;
将各阵元接收信号的FFT结果合并写成矩阵的形式,得到:
其中,
和/>分别是xm(t)、sj(t)和nm(t)的傅里叶变换后的结果;A(fn)中的元素是时延的频域表示;L为频域快拍数目;T表示转置。
3.根据权利要求1所述的一种适用于大动态范围期望信号的阵列抗干扰方法,其特征在于,步骤S2中,估计协方差矩阵包括:
阵列接收信号在频率fn处的协方差矩阵为:
其中,E表示均值,H表示共轭转置;
利用L个频域快拍计算得到协方差矩阵为:
将fn对应的索引记为NF,NF=1,2,…,NS,NS是分解的子带数;第NF个子带的协方差矩阵表示为:
其中,是信号的导向矢量,/>是/>信号的自相关结果,/>是频域噪声信号的功率,I是M维单位阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311805871.6A CN117997360A (zh) | 2023-12-26 | 2023-12-26 | 一种适用于大动态范围期望信号的阵列抗干扰方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311805871.6A CN117997360A (zh) | 2023-12-26 | 2023-12-26 | 一种适用于大动态范围期望信号的阵列抗干扰方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117997360A true CN117997360A (zh) | 2024-05-07 |
Family
ID=90891852
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311805871.6A Pending CN117997360A (zh) | 2023-12-26 | 2023-12-26 | 一种适用于大动态范围期望信号的阵列抗干扰方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117997360A (zh) |
-
2023
- 2023-12-26 CN CN202311805871.6A patent/CN117997360A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6208295B1 (en) | Method for processing radio signals that are subject to unwanted change during propagation | |
Leshem et al. | Multichannel interference mitigation techniques in radio astronomy | |
US6658234B1 (en) | Method for extending the effective dynamic range of a radio receiver system | |
US6697633B1 (en) | Method permitting increased frequency re-use in a communication network, by recovery of transmitted information from multiple cochannel signals | |
US7477190B2 (en) | Smart antenna beamforming device in communication system and method thereof | |
US6215983B1 (en) | Method and apparatus for complex phase equalization for use in a communication system | |
WO2011008146A1 (en) | Method and wideband antenna system to minimise the influence of interference sources | |
CN110361760B (zh) | 基于子空间跟踪的gnss接收机多波束指向抗干扰方法 | |
CN114755700A (zh) | 空时频多维域多波束导航抗干扰装置及抗干扰方法 | |
CN110320535A (zh) | 一种基于小波包变换与空域处理的卫星导航接收机多干扰抑制方法 | |
Myrick et al. | Anti-jam space-time preprocessor for GPS based on multistage nested Wiener filter | |
CN110557188A (zh) | 一种用于卫星通信***的抗干扰方法和装置 | |
CN111537958A (zh) | 一种宽线性降秩最小熵无失真响应的波束形成方法 | |
CN114879133A (zh) | 一种多径和高斯色噪声环境下的稀疏角度估计方法 | |
CN103701515A (zh) | 一种数字多波束形成方法 | |
CN110138413B (zh) | 一种空时域自适应宽线性降秩波束形成方法 | |
CN117997360A (zh) | 一种适用于大动态范围期望信号的阵列抗干扰方法 | |
CN113534198B (zh) | 基于协方差矩阵重构的卫星导航动态抗干扰方法及其*** | |
CN104601514B (zh) | 联合子空间投影与clean的干扰抑制方法 | |
CN107255809B (zh) | 一种基于宽带聚焦矩阵的阻塞阵波束形成方法 | |
Du et al. | A class of adaptive cyclostationary beamforming algorithms | |
CN113824488B (zh) | 基于判决反馈自适应对消的卫星通信非恶意干扰抑制方法 | |
CN113608179B (zh) | 一种雷达信号干扰抑制方法 | |
Hefnawi et al. | Performance analysis of wideband smart antenna systems using different frequency compensation techniques | |
CN115270869A (zh) | 一种基于二维空时广义旁瓣相消稳健波束形成算法 |
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 |