CN116430348A - 一种基于初相捷变脉冲串波形的空时自适应信号处理方法 - Google Patents
一种基于初相捷变脉冲串波形的空时自适应信号处理方法 Download PDFInfo
- Publication number
- CN116430348A CN116430348A CN202310699997.3A CN202310699997A CN116430348A CN 116430348 A CN116430348 A CN 116430348A CN 202310699997 A CN202310699997 A CN 202310699997A CN 116430348 A CN116430348 A CN 116430348A
- Authority
- CN
- China
- Prior art keywords
- pulse train
- initial phase
- signal processing
- following
- clutter
- 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.)
- Granted
Links
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 28
- 238000003672 processing method Methods 0.000 title claims abstract description 16
- 238000000034 method Methods 0.000 claims abstract description 41
- 239000011159 matrix material Substances 0.000 claims abstract description 22
- 238000012545 processing Methods 0.000 claims abstract description 18
- 239000013598 vector Substances 0.000 claims abstract description 16
- 238000012549 training Methods 0.000 claims abstract description 13
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 8
- 238000012360 testing method Methods 0.000 claims description 4
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 claims description 3
- 230000001629 suppression Effects 0.000 abstract description 6
- 230000000694 effects Effects 0.000 abstract description 3
- 238000005070 sampling Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 9
- 238000001514 detection method Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Classifications
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/414—Discriminating targets with respect to background clutter
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于机载雷达信号处理技术领域,具体涉及一种基于初相捷变脉冲串波形的空时自适应信号处理方法。本发明利用初相捷变脉冲串波形作为发射信号,与常规波形相比具有一定的距离选通特性,在复杂电磁环境中可获得更多的信号处理自由度;通过特征分解构建杂波子空间,相比于采样协方差矩阵求逆算法在训练样本数不足的情况下,对于通道间失真与非均匀杂波具有更好的鲁棒性;基于交替投影法更新自适应权矢量,可在信噪比损失与常规波形STAP处理相当的情况下,实现解距离模糊并进一步改善非均匀杂波的抑制效果。
Description
技术领域
本发明属于机载雷达信号处理技术领域,具体涉及一种基于初相捷变脉冲串波形的空时自适应信号处理方法。
背景技术
1973年,Brennan等人针对机载雷达杂波抑制问题,首次提出了空时自适应处理(Space-Time Adaptive Processing, STAP)的概念。经过几十年的探索和研究,STAP技术已经成为雷达领域的研究热点之一,并已发展为一门具有坚实理论基础的实用新技术。当前机载雷达采用的发射波形多为参数固定的常规波形,基于常规波形的STAP技术,在独立同分布的训练样本数充足的情况下,可以准确估计杂波协方差矩阵实现对杂波的有效抑制。但实际场景中由于存在地面起伏、地物遮挡、杂波内在运动与孤立干扰等因素的影响,造成杂波非均匀性增加,导致独立同分布的训练样本数难以满足对杂波协方差矩阵准确估计的需求,致使STAP的结果出现杂波剩余。且常规波形由于参数固定将导致接收回波出现距离模糊和杂波折叠,造成远距低速弱目标需要与折叠的剩余杂波竞争,严重影响雷达探测能力。针对杂波的非均匀性,研究人员提出了功率非均匀检测法、非均匀检测器、稀疏恢复STAP等方法。但功率非均匀检测法主要针对杂波功率非均匀的情况,无法解决运动目标干扰或孤立干扰的影响。非均匀检测器则要求估计的协方差矩阵与待检测单元的协方差矩阵一致,否则将会造成较大信杂噪比损失。稀疏恢复类的方法对信噪比要求较高,并且大多数方法都存在网格失配的问题。
综上,由于杂波非均匀性造成的剩余杂波折叠是当前影响机载雷达探测性能的主要原因。
发明内容
本发明的技术解决问题是:克服现有技术的不足即为改善机载雷达在非均匀杂波场景下的目标探测性能,提出一种基于初相捷变脉冲串波形的空时自适应信号处理方法,该方法利用初相捷变脉冲串波形的距离选通特性,通过交替投影更新自适应权矢量,相比常规波形,该方法可以解距离模糊,并进一步改善非均匀杂波的抑制效果。
本发明的技术解决方案是:
其中,为信号快时间维的复包络,/>为快时间,/>为脉冲序号,,/>为第/>个脉冲的调制初相,/>为脉冲数,/>为虚数单位,/>为脉冲重复周期,为方便叙述,定义一个脉冲重复周期/>对应的距离范围/>为一个距离段,即/>,其中/>为光速,初相捷变脉冲串波形的距离选通性能是指选通本距离段回波,抑制其他距离段回波的能力;
该方法的步骤包括:
步骤1,对初相捷变脉冲串波形下的杂波进行子空间估计,得到杂波子空间,进行杂波子空间估计的方法为:
当无目标先验信息的情况下,可以依次构建各距离段、多普勒频率、方位频率的子空间,通过遍历的方式检测场景中的目标。
步骤3,基于交替投影法计算自适应权矢量,具体方法为:
步骤4,利用步骤3得到的自适应权矢量与初相捷变脉冲串波形的回波做内积,得到基于初相捷变脉冲串波形的空时自适应信号处理结果。
有益效果
本发明利用初相捷变脉冲串波形作为发射信号,与常规波形相比具有一定的距离选通特性,在复杂电磁环境中可获得更多的信号处理自由度;通过特征分解构建杂波子空间,相比于采样协方差矩阵求逆算法在训练样本数不足的情况下,对于通道间失真与非均匀杂波具有更好的鲁棒性;基于交替投影法更新自适应权矢量,可在信噪比损失与常规波形STAP处理相当的情况下,实现解距离模糊并进一步改善非均匀杂波的抑制效果。
附图说明
图1a为线性调频脉冲串波形模糊函数的三维视图;
图1b为线性调频脉冲串波形模糊函数的距离-幅度投影图;
图2a为初相捷变脉冲串波形模糊函数的三维视图;
图2b为初相捷变脉冲串波形模糊函数的距离-幅度投影图;
图3a为基于线性调频脉冲串波形理想情况下机载雷达仿真回波数据STAP处理结果;
图3b为图3a的距离-幅度投影图;
图3c为基于初相捷变脉冲串波形理想情况下机载雷达仿真回波数据的空时自适应信号处理结果;
图3d为图3c的距离-幅度投影图;
图4a为基于线性调频脉冲串波形杂波非均匀情况下机载雷达仿真回波数据STAP处理结果;
图4b为图4a的距离-幅度投影图;
图4c为基于初相捷变脉冲串波形杂波非均匀情况下机载雷达仿真回波数据的空时自适应信号处理结果;
图4d为图4c的距离-幅度投影图。
具体实施方式
下面结合附图和实施例对本发明做进一步说明
实施例
当前机载雷达多为窄带脉冲多普勒雷达体制,发射波形大多选用参数固定的常规波形,如线性调频脉冲串波形。由图1a、图1b可知,参数固定的线性调频脉冲串波形各距离段间的回波是完全模糊的,这也是造成杂波非均匀情况下剩余杂波折叠的根本原因。
其中,为信号快时间维的复包络,/>为快时间,/>为脉冲序号,,/>为第/>个脉冲的调制初相,/>为脉冲数,/>为虚数单位,/>为脉冲重复周期,为方便叙述,定义一个脉冲重复周期/>对应的距离范围/>为一个距离段,即/>,其中/>为光速。由图2a、图2b可知,初相捷变脉冲串波形具有大约17 dB左右的距离选通性能。距离选通性能是指选通本距离段回波,抑制其他距离段回波的能力。
基于初相捷变脉冲串波形的距离选通性能,本发明的一种基于初相捷变脉冲串波形的空时自适应信号处理方法的流程,具体包括以下步骤:
步骤1,对初相捷变脉冲串波形下的杂波进行子空间估计,得到杂波子空间,进行杂波子空间估计的方法为:
当无目标先验信息的情况下,可以依次构建各距离段、多普勒频率、方位频率的子空间,通过遍历的方式检测场景中的目标。
步骤3,基于交替投影法计算自适应权矢量,具体方法为:
步骤4,利用步骤3得到的自适应权矢量与初相捷变脉冲串波形的回波做内积,得到基于初相捷变脉冲串波形的空时自适应信号处理结果。
对得到的处理结果的解距离模糊和改善非均匀杂波的抑制效果进行验证和说明:
为测试算法性能,首先对理想情况下的机载雷达探测场景进行仿真,其中理想情况是相对非均匀杂波场景提出的,即场景中不存在非均匀杂波。仿真中采用正侧视天线布局,设置雷达的不模糊测距范围为75km,两个远距低速点目标的距离和速度分别为183.75km、60m/s;270km、-40m/s,即分别位于第3距离段和第4距离段,其中速度正、负号分别代表目标远离或朝向雷达方向,目标信噪比约为15dB;仿真机载雷达杂波的速度分布范围约为-100m/s至100m/s,距离分布范围约为9km至300km,杂噪比约为60dB。图3a、图3b为基于线性调频脉冲串波形STAP处理的仿真结果,由图3a、图3b可知STAP可对杂波实现有效抑制,但目标距离存在模糊;图3c、图3d为本发明方法处理的仿真结果,由图3c、图3d可知本发明方法在对杂波有效抑制的同时可以实现对目标的无模糊测距。
其次对存在非均匀杂波的机载雷达探测场景进行仿真,仿真参数与上述仿真一致,但回波中加入场景快速变化和孤立干扰产生的非均匀杂波,图4a、图4b为基于线性调频脉冲串波形STAP处理的仿真结果,由图4a和图4b可见由于杂波非均匀导致STAP处理后存在杂波剩余,无法对目标进行有效检测;图4c、图4d为本发明方法处理的仿真结果,由图4c和图4d可见本发明方法可基于初相捷变脉冲串波形的距离选通性能,抑制杂波剩余对远距目标的影响,改善非均匀杂波的抑制效果,实现对目标的有效检测。
当然,本发明还可有其他多种实施例,在不背离本发明精神及其实质的情况下,熟悉本领域的技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (10)
1.一种基于初相捷变脉冲串波形的空时自适应信号处理方法,其特征在于该方法的步骤包括:
步骤1,对初相捷变脉冲串波形下的杂波进行子空间估计,得到杂波子空间;
步骤2,构建初相捷变脉冲串波形下的目标子空间;
步骤3,利用步骤1所得的杂波子空间构建杂波投影矩阵,利用步骤2所得的目标子空间构建目标投影矩阵;
步骤4,根据步骤3构建的杂波投影矩阵和目标投影矩阵,计算自适应权矢量;
步骤5,利用步骤4得到的自适应权矢量与初相捷变脉冲串波形的回波做内积,得到基于初相捷变脉冲串波形的空时自适应信号处理结果。
3.根据权利要求1或2所述的一种基于初相捷变脉冲串波形的空时自适应信号处理方法,其特征在于:
所述步骤1中,得到杂波子空间的方法为:
6.根据权利要求4或5所述的一种基于初相捷变脉冲串波形的空时自适应信号处理方法,其特征在于:
当无目标先验信息的情况下,依次构建各距离段、多普勒频率、方位频率的子空间,通过遍历的方式检测场景中的目标。
8.根据权利要求7所述的一种基于初相捷变脉冲串波形的空时自适应信号处理方法,其特征在于:
所述步骤4中,基于交替投影法计算自适应权矢量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310699997.3A CN116430348B (zh) | 2023-06-14 | 2023-06-14 | 一种基于初相捷变脉冲串波形的空时自适应信号处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310699997.3A CN116430348B (zh) | 2023-06-14 | 2023-06-14 | 一种基于初相捷变脉冲串波形的空时自适应信号处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116430348A true CN116430348A (zh) | 2023-07-14 |
CN116430348B CN116430348B (zh) | 2023-08-22 |
Family
ID=87089448
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310699997.3A Active CN116430348B (zh) | 2023-06-14 | 2023-06-14 | 一种基于初相捷变脉冲串波形的空时自适应信号处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116430348B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110187584A1 (en) * | 2010-01-29 | 2011-08-04 | Man-On Pun | Method for Suppressing Clutter in Space-Time Adaptive Processing Systems |
CN114442061A (zh) * | 2022-01-05 | 2022-05-06 | 北京理工大学 | 一种基于距离选通和交替反演的折叠杂波的抑制方法 |
CN116184347A (zh) * | 2022-12-29 | 2023-05-30 | 北京理工大学 | 基于距离选通与子空间投影的机载雷达折叠杂波抑制方法 |
CN116243251A (zh) * | 2022-12-13 | 2023-06-09 | 北京理工大学 | 一种基于波形分集与子空间投影的主瓣欺骗干扰抑制方法 |
-
2023
- 2023-06-14 CN CN202310699997.3A patent/CN116430348B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110187584A1 (en) * | 2010-01-29 | 2011-08-04 | Man-On Pun | Method for Suppressing Clutter in Space-Time Adaptive Processing Systems |
CN114442061A (zh) * | 2022-01-05 | 2022-05-06 | 北京理工大学 | 一种基于距离选通和交替反演的折叠杂波的抑制方法 |
CN116243251A (zh) * | 2022-12-13 | 2023-06-09 | 北京理工大学 | 一种基于波形分集与子空间投影的主瓣欺骗干扰抑制方法 |
CN116184347A (zh) * | 2022-12-29 | 2023-05-30 | 北京理工大学 | 基于距离选通与子空间投影的机载雷达折叠杂波抑制方法 |
Non-Patent Citations (2)
Title |
---|
刘子豪,任丽香,孙颖豪,范花玉: "一种捷变脉冲串的动态规划优化设计方法", 空军预警学院学报, vol. 35, no. 3, pages 157 - 161 * |
唐孝国;张剑云;: "机载相控阵雷达STAP原理及其干扰方法研究", 电子设计工程, vol. 20, no. 15, pages 71 - 77 * |
Also Published As
Publication number | Publication date |
---|---|
CN116430348B (zh) | 2023-08-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103616687B (zh) | 分段线性估计的多项式拟合isar包络对齐方法 | |
Huang et al. | Radar high speed small target detection based on keystone transform and linear canonical transform | |
CN109100718A (zh) | 基于贝叶斯学习的稀疏孔径isar自聚焦与横向定标方法 | |
CN107192991B (zh) | 一种射频辐射式仿真宽带脉冲雷达成像等效模拟方法 | |
CN110568434B (zh) | 一种多通道匀加速sar动目标二维速度估计方法 | |
CN103018727A (zh) | 一种基于样本训练的机载雷达非平稳杂波抑制方法 | |
CN110954885B (zh) | 基于sbl的捷变频雷达自适应目标重构方法 | |
CN109407070B (zh) | 一种高轨平台地面运动目标检测方法 | |
CN102621536B (zh) | 基于relax的空中多机动目标检测与参数估计方法 | |
CN109975775B (zh) | 一种雷达回波半实测数据仿真方法 | |
CN106872969A (zh) | 基于mtd脉冲积累及滑动处理的雷达目标角度估计方法 | |
CN110161472A (zh) | 一种基于信号复用的宽带车载毫米波雷达解速度模糊方法 | |
CN111551922A (zh) | 三维空间双/多基地雷达高速目标检测方法 | |
CN109031299B (zh) | 低信噪比条件下基于相位差分的isar平动补偿方法 | |
CN113655459A (zh) | 一种基于泊松盘采样的雷达无模糊多普勒扩展方法及装置 | |
CN115561728A (zh) | 一种fmcw雷达的单帧速度解模糊方法和装置 | |
CN108562901B (zh) | 基于最大信杂噪比准则的isar高分辨成像方法 | |
CN114545351A (zh) | 基于距离频率轴反转变换和二阶wvd的机动目标相参检测方法及系统 | |
CN113820678A (zh) | 一种外辐射源雷达高重频跳频信号相参积累方法 | |
CN109061626A (zh) | 一种步进频相参处理检测低信杂比动目标的方法 | |
CN115267721B (zh) | 一种基于双频sar的地面动目标径向速度估计方法 | |
CN116430348B (zh) | 一种基于初相捷变脉冲串波形的空时自适应信号处理方法 | |
CN106093927A (zh) | 基于雷达大脉宽信号的目标测速测距方法 | |
Webster et al. | Passive multistatic radar experiment using WiMAX signals of opportunity. Part 2: Multistatic velocity backprojection | |
CN112731392B (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 |