CN1297075C - 干扰阻塞空间谱估计方法 - Google Patents
干扰阻塞空间谱估计方法 Download PDFInfo
- Publication number
- CN1297075C CN1297075C CNB2003101114995A CN200310111499A CN1297075C CN 1297075 C CN1297075 C CN 1297075C CN B2003101114995 A CNB2003101114995 A CN B2003101114995A CN 200310111499 A CN200310111499 A CN 200310111499A CN 1297075 C CN1297075 C CN 1297075C
- Authority
- CN
- China
- Prior art keywords
- interference
- array
- matrix
- data
- signal
- 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.)
- Expired - Fee Related
Links
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种空间谱估计信号处理方法,根据利用广播、电视等信号工作的无源阵列雷达接收信号的特点,即无源阵列雷达接收信号包含固有的几个强干扰,且这些强相干干扰具有如下两个鲜明特点:方位已知、功率远大于飞行目标的反射信号。首先使用常规的空域滤波方法来抑制一些杂散干扰,增强目标反射信号的强度;其次利用强干扰方位已知信息形成多个干扰的抑制矩阵;最后将不同已知方位的干扰矩阵级联处理用于抑制这几个固定方位的强干扰,便于对飞行目标反射信号的检测。本发明提出的方法是一种稳健的、便于工程实现、同时兼顾抗干扰的空域处理方法,可广泛应用于阵列信号处理领域内方位已知强干扰的抑制,具有推广应用价值。
Description
本发明涉及阵列信号处理领域中的一种空间谱估计信号处理方法,特别适用于采用相控阵天线的多通道雷达信号处理系统,如无源阵列雷达信号处理系统、米波超分辨雷达信号处理系统。
空间谱估计被广泛应用于雷达、声纳、通信、地震勘探等领域,是这些领域内的一个重要研究方向,是提高参数估计性能的一种关键技术。随着科学技术的发展,现代飞行器具有机动、灵活、隐身、反辐射等高科技手段,这不仅是当代军事技术中的一项重大突破,同时也对雷达等探测器提出了严峻的挑战。对目前的雷达来说,不仅要提高对目标的探测能力,而且更应具备很强的生存能力。无源雷达正是在这一需求之下发展起来的,而利用空间复杂的电磁环境,如雷达(包括民用雷达)、通信(包括民用通信)、广播、电台等发射的信号来探测飞行目标的无源阵列雷达则是无源雷达发展的一个重要分支。
利用广播、电台等民用信号探测飞行目标的无源阵列雷达具有无源雷达的很多优点:不用主动发射机,系统构成简单,体积小;具有很好的低空性能、抗有源干扰能力;可探测隐身目标、系统工作可靠,生存能力强。但它也存在如下二个致命缺点:一、由于广播、电台的特定频率,空间多径反射复杂,阵列接收到的信号数目可能大于阵元数;二、阵列接收到的广播、电台信号中直达波信号和阵地周围特定反射点的信号强度会远大于目标反射信号。这两个缺点的存在会严重影响飞行目标反射信号的检测,因此常规的解决方法只能是增加接收系统的通道数,这不仅增加了系统的复杂性,而且会大大提高系统的成本。
本发明的目的正是针对上述背景技术中的二个缺点提出的,本发明提出的方法既避免了增加设备量的问题,又解决了有限强相干干扰抑制的问题,本发明具有工程实现方便、设备简单、且能同时兼顾抗特定有源干扰等特点。
本发明的目的是这样实现的,本发明的方法包括以下技术措施:
(1)利用自适应均衡技术对阵列的各种误差进行校正;
(2)利用空域滤波技术来抑制空间的杂散干扰,提高有用信号的信噪比。
(3)利用阵列测向的方法对干扰抑制后的数据进行测向处理。
本发明还包括以下技术措施步骤:
(1)使用天线扫描技术确定阵列工作环境中的强干扰信号方向;
(2)按如下公式构造干扰阻塞矩阵:
上式中
其中βj=ω0τj,ω0=2πf=2πc/λ,c为光速,λ为波长,d为阵元间距,J为干扰源数,f为无源阵列的工作频率,τj为第j个干扰到达等距均匀线阵相邻阵元时的波程差,当干扰入射方向为θ时,τj=dsinθj/c;
(3)利用信号相消的原理对阵列接收数据进行干扰抑制,并对干扰阻塞后的数据协方差矩阵进行去噪。
干扰抑制公式如下:
Y(t)=THX(t)
去噪公式如下:
R′=RY-THT
其中,RY=YYH/L为干扰阻塞后数据的协方差矩阵,L为数据快拍数。
本发明的目的还可以通过以下的技术措施达到:
(1)使用天线扫描技术确定阵列工作环境中的强干扰信号方向的技术措施步骤:采用精密仪器直接进行测量。
(2)使用天线扫描技术确定阵列工作环境中的强干扰信号方向还可以采用的技术措施步骤:将经过自适应均衡和空域滤波(减少进入系统的信号和干扰源数)后的数据形成协方差矩阵,利用空间平滑技术对数据进行解相干处理,再利用超分辨角度估计方法对干扰进行精确的定位。
(3)构造干扰阻塞矩阵的技术措施步骤:根据确定的干扰方向和源数分别构造TJ(共J个),注意当空间存在多个干扰时,构成的矩阵T1的维数为(M-1)×M,第二个干扰矩阵T2的维数为(M-2)×(M-1),第三个干扰矩阵T3的维数为(M-3)×(M-2),以此类推,当有J个干扰时阻塞矩阵T的维数为(M-J)×M。
本发明相比背景技术有如下优点:
(1)本发明相比常规无源阵列干扰抑制方法,不增加天线接收设备(这一设备在雷达系统中是最昂贵的),只是对系统中信号处理部分进行相关改造,所以设备简单,系统升级方便,成本低廉。
(2)本发明充分考虑空间谱估计的实际应用环境,在有效抑制杂散干扰的同时,可以实现有源干扰、多径干扰等强干扰的抑制,因此性能较好。
(3)本发明方法具有推广应用价值。
以下结合附图和实施例对本发明作进一步详细描述。
图1是本发明实施例原理方框图。
参照图1,本发明实施例由天线阵列接收系统1、自适应均衡器2、空域滤波器3、强干扰方位估计4、阻塞矩阵形成及存贮器5、强干扰对消单元6、协方差矩阵形成及去噪单元7、到达角(DOA)估计8及开关选择K组成。
实施本发明的原理如下:根据无源阵列雷达工作环境特性,接收的广播、电视信号对系统来说既是信号发生源又是干扰源,对于一个固定的阵地电视信号的直达波及周围环境的强反射点方位是确知的、且其功率远大于目标反射信号的强度。因此,检测强干扰背景下的目标必须抑制强干扰。假设无源阵列由M个阵元组成,实施例中为8个,共接收到N个远场信号(这些信号包括直达波信号、固定目标反射信号及飞行目标反射信号,为简单起见我们称前二类接收信号为干扰,只有飞行目标的反射信号才是我们希望得到的信号),且假设阵列中各阵元是各向同性的且不存在通道不一致、互耦等因素的影响,则阵列的数学模型为
X(t)=AS(t)+N(t) (1)
式中,X(t)为阵列的M×1维噪声数据矢量,N(t)为阵列的M×1维噪声数据矢量,S(t)为空间信号的N×1维矢量,A为空间阵列的M×N维流型矩阵(导向矢量阵)
A=[a1(θ1),a2(θ2),…,aN(θN)] (2)
其中,考虑等距均匀线阵时的导向矢量
ai(θ)=[1e-jβ…e-j(M-1)β]T,i=1,2,…,N (3)
式中β=2πdsin(θ)/λ,λ为波长,θ为信号和干扰的入射方向,d为阵元间距。
从式(1)可知,对消干扰最好的方法就是从式中减去干扰,而从式(2)可知从空域的角度干扰只和阵列流形矩阵A有关,而与信号波形和噪声无关,所以去除干扰的最好方法就是对A进行降秩,即将A中与干扰方向有关的导向矢量对消掉。下面以一个干扰为例说明干扰对消过程。假设干扰的角度为θ,则可以构造如下的一个矩阵
则显然有
假设式(5)中的θ=θN,则显然有β=βN,则式(5)中的最后一列全为0,则式(5)简化为
B(θ)=TH(θ)A(θ)=[b(θ1)b(θ2)…b(θN-1)0] (6)
上式中
i=1,2,…,N-1。
上述过程充分说明采用矩阵T可以实现阵列流型矩阵的降秩,即当变换矩阵左乘阵列接收数据(不考虑噪声)时有
Y(t)=THX(t)=THA(θ)S(t)
=b(θ1)s1(t)+b(θ2)s2(t)+…+b(θN-1)sN-1(t)+0 (7)
上式中Y(t)为干扰对消后的输出,显然上式中相比接收数据已经将第N个干扰对消掉了。
从上述的过程本发明的原理就是利用无源阵列环境的特性构造特殊的阻塞矩阵T(θ),利用阻塞矩阵T(θ)与阵列流型矩阵A(θ)的特殊关系实现对阵列接收数据X(t)的干扰对消。
下面结合实施例说明一下整个发明的详细步骤:
(1)由阵列接收单元1将接收到的数据X(t)存贮到系统中,这里需要注意的是,每个接收通道的快拍数L是有限制的,如果L取的过大,对后续的DOA估计是有利的,但这将造成采样数据的距离范围太大;如果L取的过小,则接收数据(特别是噪声)其统计特性会受影响,这将造成后续DOA估计性能严重下降。为了使由不满足条件引起的性能损失限制在3dB内,要求L取不少于2~3倍的系统自由度。
(2)各通道数据接收后需要经过单元2进行自适应均衡,这主要是为了校正各通道的幅相及频带不一致问题,这里采用的是常规自适应均衡技术——即32级的FIR滤波器。
(3)将经过自适应均衡后的数据送入空域滤波单元3来抑制杂散干扰,增强感兴趣目标的强度,具体的措施有二种:一是采用硬件的方法,即采用有向阵元(加工产生)或采用子阵合成方法形成;二是采用常规波束形成方法,即对阵列接收单元进行波束合成形成多波束。实施例中采用硬件的方法。
(4)根据开关选择K的情况,确定空域滤波后的数据是否进行单元4,如K选择“上”,则单元4对单元3送来的数据进行角度估计来精确估计初始角,并将估计的初始角送入单元5存贮。如果开关选择“下”,则直接通过精密仪器来确定初始方位角,并将估计的初始角送入单元5存贮。需要说明的是初始角度的估计只要开始工作时估计一次就行了。
(5)单元5根据送入的初始角度及干扰源数形成阻塞矩阵并存贮。当存在2个干扰时
注意当空间存在一个阵元时,构成的矩阵T1的维数为(M-1)×M,第二个干扰矩阵的维数为(M-2)×(M-1),第三个干扰矩阵的维数为(M-3)×(M-2),以此类推当有J个干扰时,干扰阻塞矩阵的维数为(M-J)×M。实施例中M=8。
(6)单元6从单元5中读取干扰阻塞矩阵T,并从单元3中读取阵列接收数据X(t)进行干扰相消处理。
Y(t)=THX(t),其中
(7)将干扰对消后的数据Y(t)进行数据协方差矩阵形成,并对之进行去噪声处理
(8)利用单元8进行DOA估计,具体的DOA估计算法有很多,如多重信号分类算法(MUSIC)、最小模算法(MNM)、最小方差法(MVM)、最小熵算法(MEM)、最大似然(ML)、加权子空间拟合(WSF)、旋转不变子空间(ESPRIT)等,实施例中应用MUSIC方法对角度进行估计,估计公式如下:
其中EN为R′的噪声子空间,注意式中导向矢量的维数为6×1。
对单元8的估计结果直接显示即可。
Claims (4)
1.一种干扰阻塞的空间谱估计处理方法,包括以下技术措施:
(1)利用自适应均衡技术对阵列的各种误差进行校正;
(2)利用空域滤波技术来抑制空间的杂散干扰,提高有用信号的信噪比;
(3)利用阵列测向的方法对干扰抑制后的数据进行测向处理;
其特征在于还包括以下技术措施步骤:
(1)使用天线扫描技术确定阵列工作环境中的强干扰信号方向;
(2)按如下公式构造干扰阻塞矩阵:
上式中
其中βj=ω0τj,ω0=2πf=2πc/λ,c为光速,λ为波长,d为阵元间距,
J为干扰源数,f为无源阵列的工作频率,τj为第j个干扰到达等距均
匀线阵相邻阵元时的波程差,当干扰入射方向为θ时,τj=dsinθj/c;
(3)利用信号相消的原理对阵列接收数据进行干扰抑制,并对干扰阻塞后的数据协方差矩阵进行去噪;
干扰抑制公式如下:
Y(t)=THX(t)
去噪公式如下:
R′=RY-THT
其中,X(t)为阵列接收数据,RY=YYH/L为干扰阻塞后数据的协方差矩阵,L为数据快拍数。
2.根据权利要求1所述的干扰阻塞空间谱技术方法,其特征在于使用天线扫描技术确定阵列工作环境中的强干扰信号方向的技术措施步骤:采用精密仪器直接进行测量。
3.根据权利要求1所述的干扰阻塞空间谱技术方法,其特征在于使用天线扫描技术确定阵列工作环境中的强干扰信号方向还可以采用的技术措施步骤:将经过自适应均衡和空域滤波后的数据形成协方差矩阵,以减少进入系统的信号和干扰源数,利用空间平滑技术对数据进行解相干处理,再利用超分辨角度估计方法对干扰进行精确的定位。
4.根据权利要求1所述的干扰阻塞空间谱技术方法,其特征在于构造干扰阻塞矩阵的技术措施步骤:根据确定的干扰方向和源数分别构造Tj,共J个,注意当空间存在多个干扰时,构成的矩阵T1的维数为(M-1)×M,第二个干扰矩阵T2的维数为(M-2)×(M-1),第三个干扰矩阵T3的维数为(M-3)×(M-2),以此类推,当有J个干扰时阻塞矩阵T的维数为(M-J)×M。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2003101114995A CN1297075C (zh) | 2003-12-02 | 2003-12-02 | 干扰阻塞空间谱估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2003101114995A CN1297075C (zh) | 2003-12-02 | 2003-12-02 | 干扰阻塞空间谱估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1547332A CN1547332A (zh) | 2004-11-17 |
CN1297075C true CN1297075C (zh) | 2007-01-24 |
Family
ID=34336148
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2003101114995A Expired - Fee Related CN1297075C (zh) | 2003-12-02 | 2003-12-02 | 干扰阻塞空间谱估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1297075C (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101482610B (zh) * | 2008-12-30 | 2012-11-21 | 中国人民解放军空军雷达学院 | 相控阵雷达自适应抑制特殊干扰方法 |
CN101533091B (zh) * | 2009-01-09 | 2012-09-05 | 中国人民解放军空军雷达学院 | 空时二维阻塞窄带干扰方法 |
CN101533092B (zh) * | 2009-01-09 | 2012-12-05 | 中国人民解放军空军雷达学院 | 基于功率区分的雷达杂波与多种干扰同时抑制方法 |
CN105530703B (zh) * | 2016-01-28 | 2018-10-26 | 中国铁建电气化局集团北方工程有限公司 | 一种基于波达方向估计值的gsm-r干扰源定位方法 |
CN105738860A (zh) * | 2016-03-23 | 2016-07-06 | 哈尔滨工业大学(威海) | 一种降低干扰阻塞算法伪峰的导向矢量构造方法 |
EP3881613A4 (en) | 2018-11-09 | 2022-06-08 | Lenovo (Beijing) Limited | METHODS AND DEVICES FOR POWER CONTROL |
CN111273237B (zh) * | 2019-05-29 | 2022-10-18 | 哈尔滨工程大学 | 基于空域矩阵滤波和干扰对消的强干扰抑制方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1316861A (zh) * | 2000-03-27 | 2001-10-10 | 株式会社Ntt杜可莫 | 空间的和时间的均衡器以及均衡方法 |
CN1321369A (zh) * | 1999-09-29 | 2001-11-07 | 松下电器产业株式会社 | 通信终端装置及无线通信方法 |
-
2003
- 2003-12-02 CN CNB2003101114995A patent/CN1297075C/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1321369A (zh) * | 1999-09-29 | 2001-11-07 | 松下电器产业株式会社 | 通信终端装置及无线通信方法 |
CN1316861A (zh) * | 2000-03-27 | 2001-10-10 | 株式会社Ntt杜可莫 | 空间的和时间的均衡器以及均衡方法 |
Also Published As
Publication number | Publication date |
---|---|
CN1547332A (zh) | 2004-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3589970B1 (en) | Method and system for obtaining an adaptive angle-doppler ambiguity function in mimo radars | |
US12032091B2 (en) | Signal detection and denoising systems | |
CN103901417B (zh) | L型阵列mimo雷达低复杂度空间目标二维角度估计方法 | |
CN102156279B (zh) | 基于mimo的双基地雷达地面动目标检测方法 | |
He et al. | Joint DOD and DOA estimation for MIMO array with velocity receive sensors | |
CN103901395B (zh) | 一种冲击噪声环境下相干信号波达方向动态跟踪方法 | |
CN102279387A (zh) | Mimo雷达的目标到达角估计方法 | |
CN102176008B (zh) | 一种三维地层成像的相控方位滤波方法 | |
JP2011158471A (ja) | 時空間適応処理システムにおいてターゲットを検出するための方法 | |
CN111736126B (zh) | 一种基于低成本安防雷达的信号处理方法 | |
CN101609150A (zh) | 一种提高基阵分辨力和增益的快速波束形成方法 | |
CN109765521A (zh) | 一种基于子阵划分的波束域成像方法 | |
CN112462363B (zh) | 非均匀稀疏极化阵列相干目标参数估计方法 | |
CN104698431A (zh) | 模糊分量空间角的估计及多通道sar方位解模糊方法 | |
CN109828252A (zh) | 一种mimo雷达参数估计方法 | |
CN110646765B (zh) | 一种基于黎曼距离的广义旁瓣对消算法 | |
CN1297075C (zh) | 干扰阻塞空间谱估计方法 | |
CN109471087B (zh) | 基于互质mimo雷达差集和集信号快速傅里叶变换的波达方向估计方法 | |
CN108490425B (zh) | 一种双基地mimo雷达的测角方法 | |
CN109597034B (zh) | 一种基于欧几里得距离的空时自适应处理方法 | |
CN113126076B (zh) | 基于双模式雷达的协作式距离和角度参数估计方法及装置 | |
CN109001690A (zh) | 基于馈电网络的时域空间域结合的雷达目标检测方法 | |
CN114814739A (zh) | 基于分布式fda-mimo雷达的主瓣干扰抑制方法 | |
CN113960584A (zh) | 一种基于可行点追踪的机载mimo雷达收发联合设计方法 | |
Rytel-Andrianik | Efficient matched filtering and beamforming for coherent MIMO radar |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20070124 Termination date: 20111202 |