CN106526568B - 基于短时稀疏分数阶傅里叶变换的雷达动目标检测方法 - Google Patents
基于短时稀疏分数阶傅里叶变换的雷达动目标检测方法 Download PDFInfo
- Publication number
- CN106526568B CN106526568B CN201611249041.XA CN201611249041A CN106526568B CN 106526568 B CN106526568 B CN 106526568B CN 201611249041 A CN201611249041 A CN 201611249041A CN 106526568 B CN106526568 B CN 106526568B
- Authority
- CN
- China
- Prior art keywords
- sparse
- signal
- term
- short
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 46
- 230000009466 transformation Effects 0.000 title claims abstract description 23
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 9
- 230000001052 transient effect Effects 0.000 claims abstract description 9
- 238000004458 analytical method Methods 0.000 claims abstract description 8
- 238000013461 design Methods 0.000 claims abstract description 4
- 238000005457 optimization Methods 0.000 claims abstract description 4
- 210000001367 artery Anatomy 0.000 claims description 10
- 238000009826 distribution Methods 0.000 claims description 10
- 210000003462 vein Anatomy 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 2
- 230000006698 induction Effects 0.000 claims 1
- 238000012545 processing Methods 0.000 abstract description 9
- 238000007906 compression Methods 0.000 abstract description 2
- 230000006835 compression Effects 0.000 abstract 1
- 238000009825 accumulation Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- JEIPFZHSYJVQDO-UHFFFAOYSA-N ferric oxide Chemical compound O=[Fe]O[Fe]=O JEIPFZHSYJVQDO-UHFFFAOYSA-N 0.000 description 3
- 150000001875 compounds Chemical class 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000009304 pastoral farming Methods 0.000 description 1
- 238000004080 punching Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
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/415—Identification of targets based on measurements of movement associated with the target
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及基于短时稀疏分数阶傅里叶变换(ST‑SFRFT)的雷达动目标检测方法,属于雷达信号处理和检测技术领域。其步骤包括:雷达回波解调和脉冲压缩,实现距离高分辨,并选取待检测距离单元;稀疏时频分析参数初始化,包括短时窗函数、窗长度选取和稀疏分解字典设计;ST‑SFRFT运算,进行稀疏优化求解,完成动目标回波的高分辨稀疏时频表示;遍历所有距离搜索单元,进行动目标信号稀疏域检测并估计目标运动参数。本发明采用高分辨稀疏表示方法,在时频‑稀疏域实现时变信号的高分辨、低复杂度时频表示,适于分析信号的时变非平稳特性,提高动目标信号的能量聚集性的同时抑制背景杂波,显著改善雷达动目标检测和参数估计能力。
Description
技术领域
本发明属于雷达信号处理和检测技术领域,更具体地,本发明涉及一种基于短时稀疏分数阶傅里叶变换的雷达动目标检测方法,可用于复杂环境下微弱动目标的雷达检测处理。
背景技术
雷达作为目标探测和监视的主要手段,在公共安全以及国防安全领域应用广泛。然而受复杂环境以及目标复杂运动特性的影响,目标回波极其微弱,具有低可观测性,使得雷达对杂波背景下动目标的探测性能难以满足实际需求。杂波中低可观测动目标检测技术成为关键制约因素,也是世界性难题。雷达低可观测动目标主要包括“低(低掠射角)、慢(慢速运动目标)、小(小尺寸目标)、快(机动目标)、隐(隐身目标)”等,无论是在时域还是在频域,回波信杂比(Signal-to-clutter Ratio,SCR)都很低,难以实现稳健、可靠和快速的检测。
传统的动目标检测(Moving Target Detection,MTD)方法仅适用于匀速运动目标,针对强杂波和干扰条件下的机动目标检测,雷达回波将不满足传统信号处理中平稳性要求,导致MTD方法失效。基于时频变换的雷达动目标检测方法实质是时间和多普勒维对传统MTD方法的扩展,很多方法已用于实际雷达装备中,但问题在于该类方法多为参数搜索类方法,运算效率难以满足实际要求,且参数估计精度受时频分辨率和搜索步长的限制。若目标运动特性与变换方法不相匹配,难以达到显著提升SCR的效果。目前,雷达动目标检测多从如何有效积累目标回波信号,改善SCR的角度考虑,而达到期望的增益需要大量的回波脉冲进行相参积累,对雷达实时信号处理提出严峻挑战。受目标运动以及雷达资源的限制,需要权衡积累脉冲数和积累增益,如何利用有限的回波脉冲,有效提取信号特征并用于检测成为亟需解决的问题。亟需发展和研究高时频分辨率、大数据量高效以及适用于多分量信号分析的方法和手段。
时频分析方法是将信号在一组完备的时频基上展开,如果将能够很好地刻画信号局部时频结构的时频原子构成的过完备字典替代完备基函数,则参数估计问题转化为信号的稀疏表示问题。目标雷达回波可视为少数强散射中心回波的叠加,回波具有稀疏特性。利用动目标回波信号具有稀疏性的特点,将稀疏分解的局部优化思想引入时频分析,即采用稀疏时频分布(Sparse Time-frequency Distribution,STFD)的方法对目标特性进行研究,能够有效提高算法运算效率、时频分辨率和参数估计性能,从而更有利于获得目标精细特征。因此,高分辨STFD理论为非平稳信号处理提供了新的思路和方向,也将极大提高雷达动目标的探测能力。目前,美国麻省理工学院(MIT)成立了专门的实验室对该领域进行技术攻关,但仅针对傅里叶变换的稀疏计算方法(Sparse Fourier Transform,SFT)进行了初步研究。然而,该方法不能反映信号频率的时变特性,也不能非匀速运动目标的信号处理。分数阶傅里叶变换(Fractional Fourier Transform,FRFT)对匀加速运动目标有良好的能量聚集性和检测性能,但其需要旋转角匹配搜索,运算量较大,并且缺少时域定位的功能,不能对任意时刻信号局部频率特性分析。通过构建稀疏FRFT,并加入滑动的短时窗函数,得到短时稀疏FRFT(Short-time Sparse FRFT,ST-SFRFT),能够实现高分辨时变信号时频表示的同时,改善SCR,提高复杂环境下雷达动目标检测的性能。
发明内容
本发明的目的在于突破传统时频变换采样定理的限制,提高时变信号时频分辨率,提高雷达动目标检测性能,提出一种基于短时稀疏分数阶傅里叶变换的雷达动目标检测方法。其中要解决的技术问题包括:
(1)传统MTD方法仅适用于匀速运动目标,难以有效积累杂波背景下的非匀速运动目标能量,回波信杂比低;
(2)基于时频变换的雷达动目标检测方法多为参数搜索类方法,运算效率难以满足实际要求,且参数估计精度受时频分辨率和搜索步长的限制;
(3)受目标运动以及雷达资源的限制,需要权衡积累脉冲数和积累增益,在回波脉冲有限条件下,信号积累增益低,难以有效提取目标运动特征。
本发明所述的基于短时稀疏分数阶傅里叶变换的雷达动目标检测方法,其特征在于包括以下技术措施:
步骤一、雷达回波解调和脉冲压缩,实现距离高分辨,并选取待检测距离单元;
步骤二、稀疏时频分析参数初始化,包括短时窗函数、窗长度选取和稀疏分解字典设计;
步骤三、短时稀疏分数阶傅里叶变换(ST-SFRFT)运算,进行稀疏优化求解,完成动目标回波的高分辨稀疏时频表示;
步骤四、遍历所有距离搜索单元,进行动目标信号稀疏域检测;
步骤五、目标运动参数估计。
对比现有技术,本技术方案所述的基于短时稀疏分数阶傅里叶变换的雷达动目标检测方法,有益效果在于:
(1)该方法采用高分辨稀疏表示方法,突破采样定理的限制,在时间-稀疏域实现时变信号的高分辨低复杂度时频表示;
(2)该方法在动目标的时频-稀疏域进行检测,仅保留最为稀疏的目标信号成分,因此,既保证了信号的能量聚集,又实现了时变信号的时频表示,抑制杂波,改善SCR;
(3)该方法时间-稀疏域无交叉项,适用于同一距离单元内多个运动目标的检测;
(4)该方法能够在信号时间-稀疏域,根据稀疏分解稀疏估计目标运动参数,参数估计精度高。
附图说明
附图1是基于短时稀疏分数阶傅里叶变换的雷达动目标检测方法的实施流程图。
具体实施方式
为解释本发明中涉及的动目标回波信号的稀疏时频表示,首先给出时变非平稳信号的时频分布(TFD)模型。
雷达动目标信号经过脉冲压缩后,可建模为幅度起伏的调频(FrequencyModulated,FM)信号,即
式中,a(tm)为信号的包络,tm为脉间慢时间,为信号相位。
在实际工程应用中,可根据雷达观测和信号积累时间的长短,将机动目标回波信号近似为一次调频(Linear FM,LFM)信号,描述目标的速度和加速度运动状态。尽管高阶相位信号(如,二次调频信号或多项式信号)能够更精确地描述时变信号,但后续的相参积累和参数搜索过程极其复杂,运算量大,在实际中并不常用。
对于上式的动目标信号模型,其时频分布(TFD)可表示为
式中,为信号频率的估计。
上式表示,信号可被描述为任意时刻的不同瞬时频率分量(一次或二次调频、高次调频)的叠加,在时间-频率平面内,通过稀疏分解,动目标回波信号可在处呈现冲击函数,表明能量得到积累。因此,由ρx(tm,f)构成的时频域即为信号的稀疏表示域,即稀疏时频分布。例如,对于单一分量的LFM信号,a1(tm)=1,μ为信号调频率。
以下结合说明书附图1对本发明作进一步详细描述。参照说明书附图1,本发明的处理流程分以下步骤:
(1)雷达回波解调和脉冲压缩
在相参雷达接收端,将接收并经过放大和限幅处理后得到的雷达回波数据进行距离向和方位向采样,通常距离向采样间隔等于雷达距离分辨单元,方位向采样频率等于脉冲重复频率,以保证在距离向和方位向的信号处理时间中,运动目标的回波能够被完整采集,对距离向的雷达回波数据进行解调处理,获得零中频信号sIF(t,tm),可采用雷达发射信号作为解调的参考信号
式中,t为脉内快时间,tm为脉间慢时间,tm=mTl,m=1,2,…,Tl为脉冲重复周期,sr(t,tm)为回波信号,st(t)为雷达发射信号,‘*’表示复共轭运算。将解调后的雷达回波数据进行脉冲压缩处理,得到脉内积累后的雷达回波数据sPC(t,tm),
式中,Rs(tm)为雷达与目标的视线距离,Ar(tm)是回波幅度,2Rs(tm)/c为时间延迟,B为发射信号带宽,c代表光速,tm表示脉间慢时间,λ为波长,存储距离-脉间慢时间二维数据矩阵SN×M=sPC(i,j),i=1,2,…,N,j=1,2,…,M,N为脉冲数,M为距离单元数,选取某一距离单元作为待检测单元,进行后续的处理。
(2)稀疏时频分析参数初始化
●短时窗函数、窗长度选取
窗函数可选取矩形窗,
g(τ)=1,|τ|≤Tn
Tn为窗长度,也可采用高斯窗函数,
σ为标准差,可通过调整该参数改善信号的频率分辨率。
合适的窗长度也是决定参数估计精度和分辨率的重要因素,一方面,在窗长度内需保证采样点足够多,以提高相参积累增益并保证正确估计运动参数;另一方面,窗长度越窄信号的逼近性能好,使得在有限的信号观测时间范围内,很好地描述多普勒的时变特性,因此,最佳时宽应最大程度地满足两方面的要求,对于信号
式中,为连续多项式函数,并且(n+2)阶可导,则最佳时宽为
式中,因此,对于建模为LFM信号的动目标回波,即当n=0时,窗长度应不大于调频率平方根的倒数
●稀疏分解字典设计
在雷达发射单频信号或线性调频(LFM)信号的前提下,动目标多普勒频率与速度近似成正比,复杂运动目标在一段短的时间里,可用LFM信号作为其一阶近似,因此,对于以非匀速平动为主要运动方式的动目标,如导弹、飞行器、快艇等,可采用chirp基构造稀疏分解字典,设定搜索精度和范围,假设中心频率fl的搜索范围为fl∈[0,F],搜索个数为L,中心频率分辨率为Δf=F/L,调频率μm的搜索范围为μm∈[0,K],搜索个数为M,调频率分辨率为Δμ=K/M,则构造的过完备chirp字典为L×M的矩阵
式中,gs(fl,μm)=exp(j2πflt+jπμmt2),l=1,2,...,L;m=1,2,...,M,fl=2vl/λ,vl为目标初速度分量,μm=2am/λ,am为目标加速度分量。
(3)动目标回波ST-SFRFT高分辨稀疏时频表示
对于信号的稀疏表示,集合g={gi;i=1,2,…,K},其元素是张成整个Hilbert空间H=RN的单位矢量,并且K≥N,称集合g为原子库(字典),集合中的元素为原子。对于任意信号f∈H可以展开为一组原子的线性组合,即对信号x做逼近
式中,i为原子个数,系数βi的大小表示信号与原子的相似程度。
信号x(tm)的TFD是稀疏表示的特例,其原子可由信号在TFD域中的频率估计组成,
式中,ρx(tm,f)为信号x(tm)的稀疏时频分布,h(tm)为窗函数,gi(tm,f)为稀疏表示的原子。
求解式(1)的信号稀疏表示问题可转化为最优化问题,可近似采用l1范数的最小化求解
式中,ρx∈RN,b∈RK,b为实数,o为K×N的稀疏算子。式(2)可松弛为不等约束,即
当ε=0时,式(2)和式(3)具有相同的形式。当o为分数阶傅里叶变换(FRFT)时,b为FRFT域幅值,则短时稀疏FRFT(ST-SFRFT)可表示为
式中,为ST-SFRFT时频分布,α为旋转角,u为FRFT域,(α,u)与信号中心频率f和调频率μ的关系为
式中,尺度因子fs为采样频率,则可采用步骤(2)所述的过完备chirp字典作为ST-SFRFT的稀疏分解字典。
(4)遍历所有距离搜索单元,进行动目标信号稀疏域检测
将步骤(3)的输出的二维ST-SFRFT域信号幅值作为检测统计量,与给定虚警概率下的检测门限进行比较,如果检测统计量低于检测门限,判决为没有动目标信号,继续处理后续的检测单元,若检测统计量高于检测门限,则判决为存在动目标信号,并记录ST-SFRFT域的最大峰值坐标,
式中,η为检测门限,峰值位置(α0,u0)对应动目标回波信号的中心频率和调频率(f0,μs)=(u0cscα0/S,-cotα0/S2),继续滑动时间窗h(tm),计算该时间窗内,最佳变换角α0或匹配调频率μs条件下的ST-SFRFT,从而得到不同时刻信号的瞬时频率。
(5)目标运动参数估计
根据动目标所在的ST-SFRFT域的峰值坐标估计目标的运动参数,即初速度估计值加速度估计值
Claims (5)
1.基于短时稀疏分数阶傅里叶变换的雷达动目标检测方法,其特征在于包括以下步骤:
步骤一、雷达回波解调和脉冲压缩,实现距离高分辨,并选取待检测距离单元;
步骤二、稀疏时频分析参数初始化,包括短时窗函数、窗长度选取和稀疏分解字典设计;
步骤三、短时稀疏分数阶傅里叶变换运算,进行稀疏优化求解,完成动目标回波的高分辨稀疏时频表示,设信号x(tm)的短时稀疏时频分布表示为
式中,i为原子个数,f为信号的瞬时频率,tm为脉间慢时间,系数βi(tm)的大小表示信号与原子的相似程度,h(tm)为窗函数,gi(tm,f)为稀疏表示的原子,采用l1范数的最小化求解信号的稀疏时频分布,
式中,b为实数,o为稀疏算子,上式松弛为不等约束,即
当o为分数阶傅里叶变换时,b为分数阶傅里叶变换域幅值,则得到信号短时稀疏分数阶傅里叶变换的表示,
式中,为短时稀疏分数阶傅里叶变换的时频分布,α为旋转角,u为分数阶傅里叶变换域,动目标信号将在短时稀疏分数阶傅里叶变换域形成峰值,峰值位置(α0,u0),继续滑动时间窗h(tm),计算该时间窗内,最佳变换角α0条件下的短时稀疏分数阶傅里叶变换,从而得到不同时刻信号的瞬时频率;
步骤四、遍历所有距离搜索单元,进行动目标信号稀疏域检测;
步骤五、目标运动参数估计。
2.根据权利要求1所述的基于短时稀疏分数阶傅里叶变换的雷达动目标检测方法,其特征在于步骤二所述的窗长度的选取方法为,对于信号tm为脉间慢时间,为连续多项式函数,并且(n+2)阶可导,则最佳时宽为
式中,则对于建模为线性调频信号的动目标回波,即当n=0时,窗长度Tn应不大于调频率μ平方根的倒数,
3.根据权利要求2所述的基于短时稀疏分数阶傅里叶变换的雷达动目标检测方法,其特征在于步骤二所述的稀疏分解字典设计方法为,设中心频率fl的搜索范围为fl∈[0,F],搜索个数为L,中心频率分辨率为Δf=F/L,调频率μm的搜索范围为μm∈[0,K],搜索个数为M,调频率分辨率为Δμ=K/M,则构造的过完备chirp字典为L×M的矩阵
式中,gs(fl,μm)=exp(j2πflt+jπμmt2),l=1,2,...,L;m=1,2,...,M,fl=2vl/λ,vl为目标初速度分量,μm=2am/λ,am为目标加速度分量。
4.根据权利要求1所述的基于短时稀疏分数阶傅里叶变换的雷达动目标检测方法,其特征在于步骤四所述的动目标信号稀疏域检测方法为,将二维短时稀疏分数阶傅里叶变换域信号幅值作为检测统计量,与给定虚警概率下的检测门限进行比较,如果检测统计量低于检测门限,判决为没有动目标信号,继续处理后续的检测单元,若检测统计量高于检测门限,则判决为存在动目标信号,并记录短时稀疏分数阶傅里叶变换域的最大峰值坐标,
式中,η为检测门限。
5.根据权利要求1所述的基于短时稀疏分数阶傅里叶变换的雷达动目标检测方法,其特征在于步骤五所述的目标运动参数估计方法为,检测到的动目标在短时稀疏分数阶傅里叶变换域对应的峰值位置为(α0,u0),其对应信号的中心频率和调频率(f0,μs)=(u0cscα0/S,-cotα0/S2),尺度因子Tn为窗长度,fs为采样频率,则目标的运动参数估计值为
式中,速度估计值,为加速度估计值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611249041.XA CN106526568B (zh) | 2016-12-29 | 2016-12-29 | 基于短时稀疏分数阶傅里叶变换的雷达动目标检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611249041.XA CN106526568B (zh) | 2016-12-29 | 2016-12-29 | 基于短时稀疏分数阶傅里叶变换的雷达动目标检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106526568A CN106526568A (zh) | 2017-03-22 |
CN106526568B true CN106526568B (zh) | 2019-04-12 |
Family
ID=58338530
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611249041.XA Active CN106526568B (zh) | 2016-12-29 | 2016-12-29 | 基于短时稀疏分数阶傅里叶变换的雷达动目标检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106526568B (zh) |
Families Citing this family (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106918811B (zh) * | 2017-04-05 | 2019-04-19 | 中国石油大学(华东) | 一种逆合成孔径雷达船只成像时窗选取方法 |
CN107884761A (zh) * | 2017-11-22 | 2018-04-06 | 中南大学 | 基于能量算子的多普勒穿墙雷达定位方法 |
CN108169739B (zh) * | 2017-12-27 | 2019-12-27 | 中国人民解放军战略支援部队信息工程大学 | 基于分数阶傅立叶变换和最小脉宽检测的线性调频连续波时宽比估计方法 |
CN108693509B (zh) * | 2018-04-08 | 2021-01-26 | 中国人民解放军海军航空大学 | 频控阵雷达空距频聚焦动目标积累检测方法 |
CN109283506B (zh) * | 2018-10-23 | 2022-09-30 | 中北大学 | 基于稀疏傅立叶变换的弱目标信号检测装置及方法 |
CN109919999B (zh) * | 2019-01-31 | 2021-06-11 | 深兰科技(上海)有限公司 | 一种目标位置检测的方法及装置 |
CN109856611B (zh) * | 2019-02-22 | 2020-10-27 | 清华大学 | 捷变频雷达速度-距离参数联合估计方法及装置 |
CN111983607A (zh) * | 2019-05-24 | 2020-11-24 | 西安电子科技大学 | 一种雷达三维成像方法及系统 |
CN110208785B (zh) * | 2019-07-03 | 2022-08-05 | 中国人民解放军海军航空大学 | 基于稳健稀疏分数阶傅立叶变换的雷达机动目标快速检测方法 |
CN111190157B (zh) * | 2020-01-10 | 2021-10-29 | 中国地质大学(武汉) | 一种ipix雷达回波数据时频分析方法及系统 |
CN111600821B (zh) * | 2020-04-30 | 2023-02-03 | 哈尔滨工业大学 | 一种基于分数傅里叶变换域的线性调频信号稀疏采样与重建方法 |
CN111781591B (zh) * | 2020-05-26 | 2024-02-09 | 中国空间技术研究院 | 基于分数阶傅里叶变换的雷达目标速度快速估计方法 |
CN111796242B (zh) * | 2020-06-12 | 2022-08-19 | 中国船舶重工集团公司第七二四研究所 | 一种基于功率特征值提取的改进块稀疏贝叶斯抗干扰方法 |
CN112115918A (zh) * | 2020-09-29 | 2020-12-22 | 西北工业大学 | 一种信号稀疏表示及重构的时频原子字典及信号处理方法 |
CN112462355B (zh) * | 2020-11-11 | 2023-07-14 | 西北工业大学 | 一种基于时频三特征提取的对海目标智能检测方法 |
CN112925290A (zh) * | 2021-01-22 | 2021-06-08 | 浙江大学 | 一种基于多变量本征chirp模态分解的厂级振荡检测方法 |
CN112904303B (zh) * | 2021-04-22 | 2022-06-28 | 中国人民解放军海军航空大学 | 结合杂波抑制与网格化frft处理的雷达多目标检测方法 |
CN113325407B (zh) * | 2021-06-01 | 2022-11-25 | 中国人民解放军海军航空大学 | 一种基于自适应脉冲压缩的雷达目标距离超分辨方法 |
CN113406575B (zh) * | 2021-06-17 | 2022-11-01 | 电子科技大学 | 一种基于稀疏贝叶斯学习算法的雷达距离超分辨计算方法 |
CN113534065B (zh) * | 2021-07-20 | 2024-02-02 | 中国人民解放军海军航空大学 | 一种雷达目标微动特征提取与智能分类方法及系统 |
CN113567983B (zh) * | 2021-07-23 | 2024-01-16 | 中国人民解放军海军航空大学 | 一种利用滑轨的雷达合成孔径稀疏成像方法及系统 |
CN113985381A (zh) * | 2021-10-28 | 2022-01-28 | 中国科学院空天信息创新研究院 | 回波信号的处理方法、旋翼无人机的检测方法 |
CN115840225B (zh) * | 2023-02-24 | 2023-05-05 | 中国人民解放军国防科技大学 | 基于稀疏时频分布的雷达目标高分辨成像方法和装置 |
CN117031464B (zh) * | 2023-07-25 | 2024-08-30 | 南京航空航天大学 | 一种舱内运动活体与动目标干扰区分方法和装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102866391A (zh) * | 2012-09-05 | 2013-01-09 | 中北大学 | 基于短时傅里叶变换和分数阶傅里叶变换的多目标检测方法 |
CN102914768A (zh) * | 2012-09-25 | 2013-02-06 | 中国人民解放军海军航空工程学院 | 基于形态成分分析的海面微动目标检测和特征提取方法 |
CN103176178A (zh) * | 2013-02-04 | 2013-06-26 | 中国人民解放军海军航空工程学院 | 雷达动目标Radon-分数阶傅里叶变换长时间相参积累检测方法 |
CN105827249A (zh) * | 2015-01-06 | 2016-08-03 | 南京理工大学 | 结合分数阶Fourier变换的稀疏光学采样方法 |
-
2016
- 2016-12-29 CN CN201611249041.XA patent/CN106526568B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102866391A (zh) * | 2012-09-05 | 2013-01-09 | 中北大学 | 基于短时傅里叶变换和分数阶傅里叶变换的多目标检测方法 |
CN102914768A (zh) * | 2012-09-25 | 2013-02-06 | 中国人民解放军海军航空工程学院 | 基于形态成分分析的海面微动目标检测和特征提取方法 |
CN103176178A (zh) * | 2013-02-04 | 2013-06-26 | 中国人民解放军海军航空工程学院 | 雷达动目标Radon-分数阶傅里叶变换长时间相参积累检测方法 |
CN105827249A (zh) * | 2015-01-06 | 2016-08-03 | 南京理工大学 | 结合分数阶Fourier变换的稀疏光学采样方法 |
Non-Patent Citations (2)
Title |
---|
Detection and Extraction of Target With Micromotion in Spiky Sea Clutter Via Short-Time Fractional Fourier Transform;Xiaolong Chen et al.;《IEEE Transactions on Geoscience and Remote Sensing》;20130521;1002-1018 |
稀疏域海杂波抑制与微动目标检测方法;陈小龙 等;《电子学报》;20160415;860-867 |
Also Published As
Publication number | Publication date |
---|---|
CN106526568A (zh) | 2017-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106526568B (zh) | 基于短时稀疏分数阶傅里叶变换的雷达动目标检测方法 | |
Stankovic et al. | Inverse radon transform–based micro-Doppler analysis from a reduced set of observations | |
Chen et al. | Detection and extraction of target with micromotion in spiky sea clutter via short-time fractional Fourier transform | |
Chen et al. | Time-frequency transforms for radar imaging and signal analysis | |
CN101975939A (zh) | 一种基于压缩感知理论的多普勒解模糊处理方法 | |
Pang et al. | Coherent detection algorithm for radar maneuvering targets based on discrete polynomial-phase transform | |
CN108226928B (zh) | 基于期望传播算法的逆合成孔径雷达成像方法 | |
CN102914768A (zh) | 基于形态成分分析的海面微动目标检测和特征提取方法 | |
CN106772295B (zh) | 一种用于动目标雷达截面积测量的距离与多普勒匹配方法 | |
Chen et al. | LFM signal detection and estimation based on deep convolutional neural network | |
Zhou et al. | Hough transform‐based large micro‐motion target detection and estimation in synthetic aperture radar | |
Meiyan et al. | Improved TQWT for marine moving target detection | |
CN107402382A (zh) | 基于压缩感知理论的认知被动雷达及其实现方法 | |
CN106443629A (zh) | 一种基于傅里叶基函数的雷达目标识别方法 | |
CN106093927A (zh) | 基于雷达大脉宽信号的目标测速测距方法 | |
Yu et al. | Sparse fractional Fourier transform and its applications in radar moving target detection | |
Zhang et al. | Micro-Doppler period estimation based on concentration statistics of ambiguity function | |
Chen et al. | Moving target detection at sea based on fractal characters in FRFT domain | |
Wang et al. | High range resolution profile construction exploiting modified fractional Fourier transformation | |
CN209542842U (zh) | 一种基于fpga的测风激光雷达信号处理器 | |
Xia et al. | Aircraft Target Classification Method for Conventional Narrowband Radar Based on Micro‐Doppler Effect | |
Chen et al. | Radar micro-Doppler signal detection and extraction via short-time sparse fractional Fourier transform | |
Chen et al. | Fast and Reftned Radar Processing for Maneuvering Target via Two-Stage Integration Detection | |
Cunsuo et al. | A high speed target detection approach based on STFrFT | |
Fan et al. | Target detection in non‐stationary clutter background and Riemannian geometry |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
TA01 | Transfer of patent application right | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20190314 Address after: 264001 scientific research office, two Road 188, Zhifu District, Yantai, Shandong Applicant after: Naval Aviation University of PLA Address before: 264001 scientific research department, two Road 188, Yantai, Shandong Applicant before: Naval Aeronautical Engineering Institute PLA |
|
GR01 | Patent grant | ||
GR01 | Patent grant |