CN111628790B - 一种基于干扰带宽检测的高精度抗干扰方法及装置 - Google Patents
一种基于干扰带宽检测的高精度抗干扰方法及装置 Download PDFInfo
- Publication number
- CN111628790B CN111628790B CN202010467307.8A CN202010467307A CN111628790B CN 111628790 B CN111628790 B CN 111628790B CN 202010467307 A CN202010467307 A CN 202010467307A CN 111628790 B CN111628790 B CN 111628790B
- Authority
- CN
- China
- Prior art keywords
- interference
- signal
- sampling
- bandwidth
- space
- 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
- 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/1027—Means associated with receiver for limiting or suppressing noise or interference assessing signal quality or detecting noise/interference for the received signal
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/309—Measuring or estimating channel quality parameters
- H04B17/345—Interference values
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Quality & Reliability (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
- Noise Elimination (AREA)
Abstract
本发明公开了一种基于干扰带宽检测的高精度抗干扰方法及装置,包括:接收天线阵列的时域信号,并进行采样处理;采用频域处理方法,将采样的时域信号转换到频域上来进行干扰频谱检测,通过设置噪声门限来判断信号的带宽,从而获得干扰信号频谱带宽;构建空时联合自适应抗干扰方法,利用基于修正的采样数据协方差矩阵求逆算法,基于获得的干扰信号频谱带宽,进行抗干扰处理并得到抗干扰输出信号,指导工程使用。本发明实现上来说采用基于采样数据协方差修正,其工程复杂度较不修正的情况没有明显的变化,而矩阵求逆其计算精达到理论计算精度,其抗干扰形成的零陷深度必然也能做到最优,在保证干扰抑制的同时,提升了抗干扰后接收机定位的精度。
Description
技术领域
本发明涉及抗干扰技术领域,尤其涉及一种基于干扰带宽检测的高精度抗干扰方法及装置。
背景技术
随着抗干扰技术的不断发展,采用单一的抗干扰方式已经无法满足接收机的需求,譬如时域、频域抗干扰不能对抗宽带干扰,目前工程上一般采用空时联合抗干扰技术,不仅能对抗宽带压制式信号,还能利用时域来增加抗干扰的个数,从而达到增加自由度的目的。但是目前接收机已经不再满足于单一的具有抗干扰功能,对抗干扰个数、样式、干信比的要求也在不断提高,并且还要满足抗干扰后接收机的定位精度。也就是说接收机在追求抗干扰性能的同时还追求定位精度,目前对抗干扰后信号定位精度的研究很少,几乎都是聚焦在具有较高的干信比和较多的干扰个数上,忽略了抗干扰后信号定位精度的分析。
为了提升抗干扰性能,采用基于采样协方差矩阵求逆的算法,从而使得空时联合抗干扰技术的性能达到理论抗干扰极限。不足之处在于干扰样式多样化,有单音、窄带、宽带等等,对于大带宽的宽带信号,需要通过增加时域级数来提升抗干扰性能,对于单音,窄带这样的小带宽信号,时域级数的增加会对整个频谱的完整性进行破坏,从而对有用信号进行恶化,影响其抗干扰后接收机的定位精度。
发明内容
本发明所要解决的技术问题是目前对抗干扰后信号定位精度的研究很少,几乎都是聚焦在具有较高的干信比和较多的干扰个数上,而抗干扰后信号定位精度差的问题,本发明提供了解决上述问题的一种基于干扰带宽检测的高精度抗干扰方法及装置,为了突破空、时联合抗干扰技术带来的抗干扰后接收机定位精度差的问题,本发明根据干扰信号带宽检测结果去修正采样数据的协方差矩阵,使得不同干扰样式下抗干扰后的频谱保持较高完整形,从而提升了抗干扰后接收机的定位精度,全面提升了抗干扰的整体性能,提高工程应用的可靠性。
本发明通过下述技术方案实现:
一种基于干扰带宽检测的高精度抗干扰方法,包括:
S0:接收天线阵列的时域信号,并进行采样处理;
S1:采用频域处理方法,将采样的时域信号转换到频域上来进行干扰频谱检测,通过设置噪声门限来判断信号的带宽,从而获得干扰信号频谱带宽;
S2:构建空时联合自适应抗干扰方法,利用基于修正的采样数据协方差矩阵求逆算法,基于步骤S1获得的干扰信号频谱带宽,进行抗干扰处理并得到抗干扰输出信号,提升抗干扰后接收机的定位精度,指导工程使用。
工作原理如下:
空时联合抗干扰技术通过空、时二维域的联合处理来极大的提升对抗压制式干扰的能力,为了达到该技术性能的最优,采用基于采样协方差矩阵求逆的算法来获得抗干扰后的滤波权值,在干扰方向上形成零陷,从而达到抑制干扰的目的。由于增加了时域滤波,使得在抑制干扰的同时对整个信号的频谱也有影响,严重时还可能导致干扰抑制的同时有用信号都损失了,无法定位。特别是对于像GPS这样只有2M带宽的导航信号,抗干扰后频谱的完整性对接收机的定位精度的影响就越明显。因此,为了提升各种干扰样式下抗干扰后接收机的定位精度,本发明采用基于干扰带宽检测的高精度抗干扰方法,根据干扰信号带宽检测结果去修正采样数据的协方差矩阵,从而达到减小抗干扰后信号频谱破坏的影响,提升各种干扰样式下抗干扰后接收机的定位精度。
本发明方法考虑到空时联合抗干扰技术,采用基于采样数据协方差矩阵求逆的抗干扰方法来实现,能使得空时联合抗干扰技术形成的干扰零陷达到理论最佳。但是却没有考虑时域滤波带来的抗干扰后频谱的影响,只保证了干扰零陷达到理论最深有效抑制干扰,忽略了抗干扰后接收机定位精度的影响。本发明为了达到较好的干扰抑制,并保持整个有用信号带宽内频谱的完整性,采用基于修正的采样数据协方差矩阵求逆抗干扰方法,提升各种干扰样式下抗干扰后接收机的定位精度,使得抗干扰后接收机的定位精度不再受限于干扰样式,增强工程实用性。
本发明基于干扰带宽检测的高精度抗干扰方法分为:干扰带宽检测和抗干扰两个部分,干扰带宽检测采用频域处理方案,将采样的时域信号转换到频域上来检测,通过设置噪声门限来判断信号的带宽;抗干扰部分采用空时联合自适应抗干扰方案,实现上采用基于修正的采样数据协方差矩阵求逆算法;具体实现上来说采用基于采样数据协方差修正,其工程复杂度较不修正的情况没有明显的变化,而矩阵求逆其计算精达到理论计算精度,其抗干扰形成的零陷深度必然也能做到最优,因此在保证干扰抑制的同时,提升了抗干扰后接收机定位的精度。基于干扰带宽检测的高精度抗干扰方法可以应用于同类或者其他产品上,不能规避此技术手段。
进一步地,步骤S1包括以下子步骤:
S101:当外界不施加干扰时,接收AD采样信号,将其从时域变换到频域;选择若干个点进行傅里叶变换FFT变换,变换后得到实部信号Fdatai(1)、Fdatai(2)...Fdatai(n),虚部信号 Fdataq(1)、Fdataq(2)...Fdataq(n),n≥512,转步骤S102;
S102:计算各个频域信号的幅度,频域信号的幅度计算公式如下
Famplitude=|Fdatai|+|Fdataq| (1)
其中| |表示对信号取绝对值;
由公式(1)得到频域信号的幅度Famplitude(1)、Famplitude(2)...Famplitude(n),转步骤 S103;
其中∑表示信号累加,转步骤S104;
S104:重复步骤S101到S103,得到十次计算的判决门限值Threshold(1)、Threshold(2)... Threshold(10),然后对10次判决门限结果进行平均,得到最终的门限值;
将得到的频域信号幅度值与公式(3)得到的门限值Th r进行比较,若满足Famplitude>(Thr+Thr/4)则标记为有效;只判断前个点,若连续超过20个点的幅度值大于该门限,则给出修正无效的标志即flag_revise_valid=0,否则flag_revise_valid=1,转步骤S2。
进一步地,步骤2包括以下子步骤:
S201:基于MVDR的抗干扰算法滤波器,计算最优权值系数,计算公式为:
其中wopt为最优权值系数;a(θ0)为[1,0,...0]';Rxx为采样数据协方差矩阵,具体地:
Rxx=E[X(t)XH(t)] (5)
其中X为采样数据;
S201具体包括:
步骤A:对AD采样信号进行预处理,得到IQ信号,转步骤B;
步骤B:将空域信号转化成空时联合的信号,将空域信号进行时域延迟,得到第m个天线接收信号形成的空时联合的信号Xm=[xm1,xm2,…,xmN],总共有M个天线,那么有 X1,X2,…,XM的路空时联合信号,转步骤C;
步骤C:将空时联合的信号进行采样协方差矩阵计算,形成采样协方差复数矩阵Rxx,计算公式如下:
其中K表示采样点数,M表示阵元个数,N表示时域滤波器级数,转步骤D;
步骤D:根据步骤S104得到flag_revise_valid,根据它的取值决定是否修正采样协方差矩阵Rxx,若flag_revise_valid=1,则Rxx修正为:
否则Rxx保持不变,转步骤E;
转步骤S203;
S204:利用生成的抗干扰权值wopt,对空时二维域信号进行滤波加权,得到抗干扰输出信号。
进一步地,步骤A中对AD采样信号进行预处理,包括去直流、数字IQ变换。
另一方面,本发明还提供了一种基于干扰带宽检测的高精度抗干扰装置,该装置支持所述的一种基于干扰带宽检测的高精度抗干扰方法,该装置包括AD采样模块、干扰带宽检测模块、抗干扰处理模块,其中:
所述AD采样模块,用于接收天线阵列的时域信号,并进行采样处理;
所述干扰带宽检测模块,用于采用频域处理方法,将所述AD采样模块采样的时域信号转换到频域上来进行干扰频谱检测,通过设置噪声门限来判断信号的带宽,从而获得干扰信号频谱带宽;
所述抗干扰处理模块,用于构建空时联合自适应抗干扰方法,利用基于修正的采样数据协方差矩阵求逆算法,基于所述干扰带宽检测模块获得的干扰信号频谱带宽,进行抗干扰处理并得到抗干扰输出信号,提升抗干扰后接收机的定位精度,指导工程使用。
本发明与现有技术相比,具有如下的优点和有益效果:
1、本发明为了突破空、时联合抗干扰技术带来的抗干扰后接收机定位精度差的问题,本发明根据干扰信号带宽检测结果去修正采样数据的协方差矩阵,使得不同干扰样式下抗干扰后的频谱保持较高完整形,从而提升了抗干扰后接收机的定位精度,全面提升了抗干扰的整体性能,提高工程应用的可靠性;
2、本发明基于干扰带宽检测的高精度抗干扰方法分为:干扰带宽检测和抗干扰两个部分,干扰带宽检测采用频域处理方案,将采样的时域信号转换到频域上来检测,通过设置噪声门限来判断信号的带宽;抗干扰部分采用空时联合自适应抗干扰方案,实现上采用基于修正的采样数据协方差矩阵求逆算法;具体实现上来说采用基于采样数据协方差修正,其工程复杂度较不修正的情况没有明显的变化,而矩阵求逆其计算精达到理论计算精度,其抗干扰形成的零陷深度必然也能做到最优,因此在保证干扰抑制的同时,提升了抗干扰后接收机定位的精度;基于干扰带宽检测的高精度抗干扰方法可以应用于同类或者其他产品上。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解,构成本申请的一部分,并不构成对本发明实施例的限定。在附图中:
图1为本发明一种基于干扰带宽检测的高精度抗干扰方法流程图。
图2为本发明预处理流程图。
图3为本发明空时二维自适应滤波器实现框图。
图4为本发明中步骤S2详细步骤流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
实施例1
如图1至图4所示,空时联合抗干扰技术通过空、时二维域的联合处理来极大的提升对抗压制式干扰的能力,为了达到该技术性能的最优,采用基于采样协方差矩阵求逆的算法来获得抗干扰后的滤波权值,在干扰方向上形成零陷,从而达到抑制干扰的目的。由于增加了时域滤波,使得在抑制干扰的同时对整个信号的频谱也有影响,严重时还可能导致干扰抑制的同时有用信号都损失了,无法定位。特别是对于像GPS这样只有2M带宽的导航信号,抗干扰后频谱的完整性对接收机的定位精度的影响就越明显。因此,为了提升各种干扰样式下抗干扰后接收机的定位精度,本发明采用基于干扰带宽检测的高精度抗干扰方法,根据干扰信号带宽检测结果去修正采样数据的协方差矩阵,从而达到减小抗干扰后信号频谱破坏的影响,提升各种干扰样式下抗干扰后接收机的定位精度。
本发明方法考虑到空时联合抗干扰技术,采用基于采样数据协方差矩阵求逆的抗干扰方法来实现,能使得空时联合抗干扰技术形成的干扰零陷达到理论最佳。但是却没有考虑时域滤波带来的抗干扰后频谱的影响,只保证了干扰零陷达到理论最深有效抑制干扰,忽略了抗干扰后接收机定位精度的影响。本发明为了达到较好的干扰抑制,并保持整个有用信号带宽内频谱的完整性,采用基于修正的采样数据协方差矩阵求逆抗干扰方法,提升各种干扰样式下抗干扰后接收机的定位精度,使得抗干扰后接收机的定位精度不再受限于干扰样式,增强工程实用性。
本发明基于干扰带宽检测的高精度抗干扰方法分为:干扰带宽检测和抗干扰两个部分,干扰带宽检测采用频域处理方案,将采样的时域信号转换到频域上来检测,通过设置噪声门限来判断信号的带宽;抗干扰部分采用空时联合自适应抗干扰方案,实现上采用基于修正的采样数据协方差矩阵求逆算法;本发明方案的具体实施流程图,见图1所示。方法中信号采样时钟为62m,下面对整个实现步骤进行详细描述。
本发明一种基于干扰带宽检测的高精度抗干扰方法,包括:
S0:接收天线阵列的时域信号,并进行采样处理;
S1:采用频域处理方法,将采样的时域信号转换到频域上来进行干扰频谱检测,通过设置噪声门限来判断信号的带宽,从而获得干扰信号频谱带宽;步骤S1具体包括:
S101:当外界不施加干扰时,接收AD采样信号,将其从时域变换到频域;由于不需要特别精准的带宽检测指标,本实施例选择512个点进行傅里叶变换FFT变换,变换后得到实部信号Fdatai(1)、Fdatai(2)...Fdatai(512),虚部信号Fdataq(1)、Fdataq(2)...Fdataq(512),转步骤S102;
S102:计算各个频域信号的幅度,频域信号的幅度计算公式如下
Famplitude=|Fdatai|+|Fdataq| (1)
其中| |表示对信号取绝对值;
由公式(1)得到频域信号的幅度Famplitude(1)、Famplitude(2)...Famplitude(512),转步骤S103;
S103:AD采样信号中频信号,基于FFT变换得到的信号频谱图是对称的,所以只选取前256个点的幅度值来计算判决门限值;
其中∑表示信号累加,n取512,转步骤S104;
S104:重复步骤S101到S103,得到十次计算的判决门限值Threshold(1)、Threshold(2)... Threshold(10),然后对10次判决门限结果进行平均,得到最终的门限值;
将得到的频域信号幅度值与公式(3)得到的门限值Thr进行比较,若满足 F amplitude>(Thr+Thr/4)则标记为有效;只判断前256个点,若连续超过20个点的幅度值大于该门限,则给出修正无效的标志即flag_revise_valid=0,否则flag_revise_valid=1,转步骤S2。
S2:构建空时联合自适应抗干扰方法,利用基于修正的采样数据协方差矩阵求逆算法,基于步骤S1获得的干扰信号频谱带宽,进行抗干扰处理并得到抗干扰输出信号,提升抗干扰后接收机的定位精度,指导工程使用。如图4所示,步骤2包括以下子步骤:
S201:空时联合二维自适应滤波器实现框图见图3,基于MVDR的抗干扰算法滤波器,计算最优权值系数,计算公式为:
其中wopt为最优权值系数;a(θ0)为[1,0,...0]';Rxx为采样数据协方差矩阵,具体地:
Rxx=E[X(t)XH(t)] (5)
其中X为采样数据;其流程图见图4,详细步骤如下。S201具体包括:
步骤A:对AD采样信号进行预处理,包括去直流、数字IQ变换,见图2,预处理后得到IQ信号,转步骤B;
步骤B:将空域信号转化成空时联合的信号,根据图2的方式,将空域信号进行时域延迟,得到第m个天线接收信号形成的空时联合的信号Xm=[xm1,xm2,…,xmN],总共有M个天线,那么有X1,X2,…,XM的路空时联合信号,转步骤C;
步骤C:将空时联合的信号进行采样协方差矩阵计算,形成采样协方差复数矩阵Rxx,计算公式如下:
其中K表示采样点数,本实施例中K为2048个采样点,M表示阵元个数,N表示时域滤波器级数,转步骤D;
步骤D:根据步骤S104得到flag_revise_valid,根据它的取值决定是否修正采样协方差矩阵Rxx,若flag_revise_valid=1,则Rxx修正为:
否则Rxx保持不变,转步骤E;
转步骤S203;
S204:利用生成的抗干扰权值wopt,对空时二维域信号进行滤波加权见图2,得到抗干扰输出信号。
本发明方法中采用如下六个核心关键设计点:
(1)干扰带宽门限设定,首先将AD采样信号进行时频转换,然后获得频域信号的幅度谱值,根据无干扰条件下信号的幅度谱值得到判决门限,以此门限来检测干扰带宽;
(2)干扰带宽门限设定,以连续的20个点幅度值大于门限1.25倍来给出采样协方差矩阵Rxx是否修正的标志。
(3)基于干扰带宽修正采样协方差矩阵Rxx,若给出修正有效,则只需修正参考阵元相关的时域元素,其余阵元的时域元素都保持不变。本发明中参考阵元由约束a(θ0)=[1,0,...0]'决定,因此修正Rxx中除去参考阵元所在的行,其余与参考阵元相关的时域形成的互相关元素置为零,自相关元素保持不变,得到公式(7)。
(4)Rxx的修正还可以扩展到其它参考阵元,若a(θ0)=[0,...1,0...0]'时,则修正与a(θ0)相应的非零的参考阵元的时域元素,具体的方法与本发明中一致。
(6)抗干扰权值wopt的选择还可以扩展的到其它列,只需要和a(θ0)中的非零元素相对应即可。
本发明根据以上六点进行处理后基于干扰带宽检测的高精度抗干扰算法,在抗干扰的同时还考虑了不同干扰带宽下抗干扰后的信号频谱,突破工程抗干扰算法接收机定位精度差的问题,提升抗干扰算法的性能,增强工程应用的实用性,尤其是对于GPS这样小带宽的信号使得其在满足干扰抑制的同时还能保持较好的频谱完整性。本发明方法实现上来说采用基于采样数据协方差修正,其工程复杂度较不修正的情况没有明显的变化,而矩阵求逆其计算精达到理论计算精度,其抗干扰形成的零陷深度必然也能做到最优,因此在保证干扰抑制的同时,提升了抗干扰后接收机定位的精度。基于干扰带宽检测的高精度抗干扰方法可以应用于同类或者其他产品上,不能规避此技术手段。
具体实施时,步骤如图1所示:
根据AD采样信号进行时频变换,得到频域的频谱幅度值;
在无干扰条件下,利用频谱幅度值获得判决门限;
若有连续的20个点的频谱幅度值高于判决门限的1.25倍,则给出修正有效的标志,否则修正无效。
信号预处理见图2,首先对接收到的阵元信号进行去直流,然后再进行数字IQ变换,得到IQ信号。因为算法需要利用信号的相位进行权值调整,在天线方向图上形成零陷,所以这一步预处理是必须的。
采样数据形成的协方差矩阵Rxx修正方法,首先利用干扰信号的带宽检测标志来判定是否修。其次利用约束条件a(θ0)来判断Rxx中的哪些元素需要修正,在不改变原有算法的复杂的前提下,只需要修正a(θ0)对应的非零元素相关的时域互相关元素的值,将其置零。方案中选择的是第一个阵元的第一个采样信号作为参考,因此只修正参考信号相关的时域互相关元素的值,即得到式子(7)。
抗干扰算法实现上采用的是基于采样协方差矩阵求逆来实现的,首先将复数矩阵升维转换为实数矩阵,然后利用高斯消元法进行逆矩阵求解。
利用生成的抗干扰权值wopt,对空时二维域信号进行滤波加权见图2,得到抗干扰输出信号。
图2为空时二维抗干扰方法的自适应滤波框图,从框图可以看出,将空域信号扩展到了空时二维域上,利用抗干扰算法的权值对空时二维域的信号进行滤波加权,输出抗干扰后的信号。
图3为基于干扰带宽检测的高精度抗干扰方法的实现流程图,关键部分有干扰带宽检测,采样数据的协方差矩阵生成、采样数据协方差矩阵的修正、复数矩阵转化为实数矩阵、权值选取、信号滤波等处理。
图4为步骤S2的详细流程图。
考虑到空时联合抗干扰技术,采用基于采样数据协方差矩阵求逆的抗干扰方法来实现,能使得空时联合抗干扰技术形成的干扰零陷达到理论最佳。但是却没有考虑时域滤波带来的抗干扰后频谱的影响,只保证了干扰零陷达到理论最深有效抑制干扰,忽略了抗干扰后接收机定位精度的影响。本发明为了达到较好的干扰抑制,并保持整个有用信号带宽内频谱的完整性,采用基于修正的采样数据协方差矩阵求逆抗干扰方法,提升各种干扰样式下抗干扰后接收机的定位精度,使得抗干扰后接收机的定位精度不再受限于干扰样式,增强工程实用性。
本发明基于干扰带宽检测的高精度抗干扰方法分为:干扰带宽检测和抗干扰两个部分,干扰带宽检测采用频域处理方案,将采样的时域信号转换到频域上来检测,通过设置噪声门限来判断信号的带宽;抗干扰部分采用空时联合自适应抗干扰方案,实现上采用基于修正的采样数据协方差矩阵求逆算法;具体实现上来说采用基于采样数据协方差修正,其工程复杂度较不修正的情况没有明显的变化,而矩阵求逆其计算精达到理论计算精度,其抗干扰形成的零陷深度必然也能做到最优,因此在保证干扰抑制的同时,提升了抗干扰后接收机定位的精度。基于干扰带宽检测的高精度抗干扰方法可以应用于同类或者其他产品上,不能规避此技术手段。
实施例2
如图1至图4所示,本实施例与实施例1的区别在于,本实施例提供了一种基于干扰带宽检测的高精度抗干扰装置,该装置支持实施例1中的一种基于干扰带宽检测的高精度抗干扰方法,该装置包括AD采样模块、干扰带宽检测模块、抗干扰处理模块,其中:
所述AD采样模块,用于接收天线阵列的时域信号,并进行采样处理;
所述干扰带宽检测模块,用于采用频域处理方法,将所述AD采样模块采样的时域信号转换到频域上来进行干扰频谱检测,通过设置噪声门限来判断信号的带宽,从而获得干扰信号频谱带宽;
所述抗干扰处理模块,用于构建空时联合自适应抗干扰方法,利用基于修正的采样数据协方差矩阵求逆算法,基于所述干扰带宽检测模块获得的干扰信号频谱带宽,进行抗干扰处理并得到抗干扰输出信号,提升抗干扰后接收机的定位精度,指导工程使用。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种基于干扰带宽检测的高精度抗干扰方法,其特征在于,包括:
S0:接收天线阵列的时域信号,并进行采样处理;
S1:采用频域处理方法,将采样的时域信号转换到频域上来进行干扰频谱检测,通过设置噪声门限来判断信号的带宽,从而获得干扰信号频谱带宽;
S2:构建空时联合自适应抗干扰方法,利用基于修正的采样数据协方差矩阵求逆算法,基于步骤S1获得的干扰信号频谱带宽,进行抗干扰处理并得到抗干扰输出信号,指导工程使用;
步骤S1包括以下子步骤:
S101:当外界不施加干扰时,接收AD采样信号,将其从时域变换到频域;选择若干个点进行傅里叶变换FFT变换,变换后得到实部信号Fdatai(1)、Fdatai(2)...Fdatai(n),虚部信号Fdataq(1)、Fdataq(2)...Fdataq(n),n≥512,转步骤S102;
S102:计算各个频域信号的幅度,频域信号的幅度计算公式如下
Famplitude=|Fdatai|+|Fdataq| (1)
其中||表示对信号取绝对值;
由公式(1)得到频域信号的幅度Famplitude(1)、Famplitude(2)...Famplitude(n),转步骤S103;
其中∑表示信号累加,转步骤S104;
S104:重复步骤S101到S103,得到十次计算的判决门限值Threshold(1)、Threshold(2)...Threshold(10),然后对10次判决门限结果进行平均,得到最终的门限值;
将得到的频域信号幅度值与公式(3)得到的门限值Thr进行比较,若满足Famplitude>(Thr+Thr/4)则标记为有效;只判断前个点,若连续超过20个点的幅度值大于该门限,则给出修正无效的标志即flag_revise_valid=0,否则flag_revise_valid=1,转步骤S2;
步骤2包括以下子步骤:
S201:基于MVDR的抗干扰算法滤波器,计算最优权值系数,计算公式为:
其中wopt为最优权值系数;a(θ0)为[1,0,...0]';Rxx为采样数据协方差矩阵,具体地:
Rxx=E[X(t)XH(t)] (5)
其中X为采样数据;
S201具体包括:
步骤A:对AD采样信号进行预处理,得到IQ信号,转步骤B;
步骤B:将空域信号转化成空时联合的信号,将空域信号进行时域延迟,得到第m个天线接收信号形成的空时联合的信号Xm=[xm1,xm2,…,xmN],总共有M个天线,那么有X1,X2,…,XM的路空时联合信号,转步骤C;
步骤C:将空时联合的信号进行采样协方差矩阵计算,形成采样协方差复数矩阵Rxx,计算公式如下:
其中K表示采样点数,M表示阵元个数,N表示时域滤波器级数,转步骤D;
步骤D:根据步骤S104得到flag_revise_valid,根据它的取值决定是否修正采样协方差矩阵Rxx,若flag_revise_valid=1,则Rxx修正为:
否则Rxx保持不变,转步骤E;
转步骤S203;
S204:利用生成的抗干扰权值wopt,对空时二维域信号进行滤波加权,得到抗干扰输出信号。
2.根据权利要求1所述的一种基于干扰带宽检测的高精度抗干扰方法,其特征在于,步骤A中对AD采样信号进行预处理,包括去直流、数字IQ变换。
3.一种基于干扰带宽检测的高精度抗干扰装置,其特征在于,该装置支持如权利要求1至2中任意一项所述的一种基于干扰带宽检测的高精度抗干扰方法,该装置包括AD采样模块、干扰带宽检测模块、抗干扰处理模块,其中:
所述AD采样模块,用于接收天线阵列的时域信号,并进行采样处理;
所述干扰带宽检测模块,用于采用频域处理方法,将所述AD采样模块采样的时域信号转换到频域上来进行干扰频谱检测,通过设置噪声门限来判断信号的带宽,从而获得干扰信号频谱带宽;
所述抗干扰处理模块,用于构建空时联合自适应抗干扰方法,利用基于修正的采样数据协方差矩阵求逆算法,基于所述干扰带宽检测模块获得的干扰信号频谱带宽,进行抗干扰处理并得到抗干扰输出信号,提升抗干扰后接收机的定位精度,指导工程使用。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010467307.8A CN111628790B (zh) | 2020-05-28 | 2020-05-28 | 一种基于干扰带宽检测的高精度抗干扰方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010467307.8A CN111628790B (zh) | 2020-05-28 | 2020-05-28 | 一种基于干扰带宽检测的高精度抗干扰方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111628790A CN111628790A (zh) | 2020-09-04 |
CN111628790B true CN111628790B (zh) | 2021-06-29 |
Family
ID=72272885
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010467307.8A Active CN111628790B (zh) | 2020-05-28 | 2020-05-28 | 一种基于干扰带宽检测的高精度抗干扰方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111628790B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6964833B2 (ja) * | 2019-08-23 | 2021-11-10 | 三菱電機株式会社 | 適応制御装置、適応信号処理装置及びアダプティブアレイアンテナシステム |
CN113037351A (zh) * | 2021-02-08 | 2021-06-25 | 西安电子科技大学 | 用于赋形天线的空时联合抗宽带干扰方法 |
WO2022183326A1 (zh) * | 2021-03-01 | 2022-09-09 | 深圳市大疆创新科技有限公司 | 滤波方法、装置、可移动平台和存储介质 |
WO2022205199A1 (zh) * | 2021-03-31 | 2022-10-06 | 华为技术有限公司 | 一种干扰处理方法和装置 |
CN115842568B (zh) * | 2023-02-14 | 2023-05-12 | 中国电子科技集团公司第十五研究所 | 去干扰通信方法和窄带无线通信芯片 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105099964A (zh) * | 2014-05-07 | 2015-11-25 | 电信科学技术研究院 | 一种功率强度估计和确定干扰小区的方法及设备 |
CN107356944A (zh) * | 2017-06-01 | 2017-11-17 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 提高卫星导航阵列天线抗干扰性能的方法 |
CN109959899A (zh) * | 2018-04-12 | 2019-07-02 | 哈尔滨工程大学 | 特征投影预处理及协方差矩阵稀疏重构主瓣干扰抑制算法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104914451B (zh) * | 2015-05-05 | 2017-03-01 | 中国人民解放军国防科学技术大学 | 一种块Toeplitz矩阵低复杂度求逆的空时抗干扰方法 |
US20200116851A1 (en) * | 2018-10-10 | 2020-04-16 | Government Of The United States, As Represented By The Secretary Of The Air Force | Robust Constant False Alarm Rate (CFAR) Detector for Interference-Plus-Noise Covariance Matrix Mismatch |
CN110261826A (zh) * | 2019-05-17 | 2019-09-20 | 北京理工大学 | 一种零陷展宽的相干干扰抑制方法 |
CN110850445B (zh) * | 2019-11-21 | 2024-04-05 | 中国人民解放军63961部队 | 一种基于空时采样协方差求逆的脉冲干扰抑制方法 |
CN111162878B (zh) * | 2019-12-25 | 2020-11-17 | 成都天奥信息科技有限公司 | 一种基于子阵降维带约束的多域联合抗干扰方法 |
-
2020
- 2020-05-28 CN CN202010467307.8A patent/CN111628790B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105099964A (zh) * | 2014-05-07 | 2015-11-25 | 电信科学技术研究院 | 一种功率强度估计和确定干扰小区的方法及设备 |
CN107356944A (zh) * | 2017-06-01 | 2017-11-17 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 提高卫星导航阵列天线抗干扰性能的方法 |
CN109959899A (zh) * | 2018-04-12 | 2019-07-02 | 哈尔滨工程大学 | 特征投影预处理及协方差矩阵稀疏重构主瓣干扰抑制算法 |
Also Published As
Publication number | Publication date |
---|---|
CN111628790A (zh) | 2020-09-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111628790B (zh) | 一种基于干扰带宽检测的高精度抗干扰方法及装置 | |
CN106772457B (zh) | 一种基于空时频架构的卫星导航抗干扰方法 | |
CN109407055B (zh) | 基于多径利用的波束形成方法 | |
CN110007148B (zh) | 一种基于离散频谱相位和幅值综合内插的单频信号频率估计方法 | |
CN106646350B (zh) | 一种单只矢量水听器各通道幅度增益不一致时的修正方法 | |
CN111948598B (zh) | 空域干扰信号检测方法与装置 | |
CN101482610A (zh) | 相控阵雷达自适应抑制特殊干扰技术 | |
CN104835503A (zh) | 一种改进gsc自适应语音增强方法 | |
CN111693971B (zh) | 一种用于弱目标检测的宽波束干扰抑制方法 | |
CN106443721A (zh) | 一种快速收敛的卫星导航天线阵抗干扰方法 | |
CN111580042B (zh) | 一种基于相位优化的深度学习测向方法 | |
CN111624573A (zh) | 一种海杂波背景下的时域自适应目标检测方法 | |
CN110196407B (zh) | 一种基于频率预估的单矢量水听器信号来波方向估计方法 | |
CN111175727B (zh) | 一种基于条件波数谱密度的宽带信号方位估计的方法 | |
CN110376545B (zh) | 一种单参数调优随机共振增强的单矢量水听器定向方法 | |
CN114545342A (zh) | 利用多通道侦察接收机的雷达脉冲信号参数测量方法 | |
CN104715758A (zh) | 一种分部处理式阵列式语音定位和增强方法 | |
CN107783096A (zh) | 一种用于方位历程图显示的二维背景均衡方法 | |
CN104880714B (zh) | 一种基于两级滤波结构的gnss天线阵抗干扰方法 | |
CN109669172B (zh) | 基于主瓣内强干扰抑制的弱目标方位估计方法 | |
CN109061597B (zh) | 基于盲源分离与时频脊波域滤波的电离层杂波抑制方法 | |
CN107255809B (zh) | 一种基于宽带聚焦矩阵的阻塞阵波束形成方法 | |
CN115980721A (zh) | 一种无误差协方差矩阵分离的阵列自校正方法 | |
CN113406672B (zh) | 一种基于自适应加窗的分段线性调频干扰消除方法 | |
CN111722178B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |