一种改进的分数阶傅里叶变换机动弱目标检测方法技术

本发明专利技术公开了一种改进的分数阶傅里叶变换机动弱目标检测方法。本发明专利技术首先建立高速机动弱目标匀加速运动信号模型；接着对回波原始数据在快时间域进行傅里叶变换，变换后的数据经过Keystone变换后，校正由一次相位引起的距离单元走动问题，接着在快时间域再进行IFFT，在得到的最终回波数据中对每一个距离单元分别求p[0,2]内所有阶次的FRFT，形成信号能量在由分数阶p和分数阶变换域u组成二维参数平面（p,u）的二维分布，在此平面上进行峰值点的二维搜索，即可实现加速目标的检测。本发明专利技术可以实现低信噪比背景环境下机动弱目标的检测，同时具有较好的可靠性、可行性和实时性。

【技术实现步骤摘要】
【专利摘要】本专利技术公开了。本专利技术首先建立高速机动弱目标匀加速运动信号模型；接着对回波原始数据在快时间域进行傅里叶变换，变换后的数据经过Keystone变换后，校正由一次相位引起的距离单元走动问题，接着在快时间域再进行IFFT，在得到的最终回波数据中对每一个距离单元分别求p内所有阶次的FRFT，形成信号能量在由分数阶p和分数阶变换域u组成二维参数平面（p,u）的二维分布，在此平面上进行峰值点的二维搜索，即可实现加速目标的检测。本专利技术可以实现低信噪比背景环境下机动弱目标的检测，同时具有较好的可靠性、可行性和实时性。【专利说明】
本专利技术属于雷达信号处理领域，涉及。
技术介绍
目前机载预警雷达信号处理常采用长时间相参积累方式，以增加实际利用的目标回波信号能量提高雷达检测性能。传统的相参积累都是通过FFT来完成的，其前提是在波束驻留时间内，目标近似为匀速直线运动，多普勒频率为常数，且回波包络的移动不超过一个距离门。但是对于雷达机动弱目标，由于目标可能具有强机动性，在整个积累过程中多普勒频率不再为常数，不仅存在距离徙动，而且还存在多普勒频率徙动。在积累时间内，如果目标运动的距离超过半个距离单元，则会出现距离徙动，如果目标多普勒频率的变化超过一个多普勒频率单元，则会出现多普勒频率徙动。目标运动速度越快，加速度越大，多普勒徙动越明显。本专利技术同时采用了距离补偿和多普勒频率补偿技术。针对距离补偿问题，通常使用互相关法，最小熵法和检测前跟踪方法等，检测前跟踪方法是一种非相参积累方法，积累效率不如相参积累方法，互相关法和最小熵法对回波信噪比要求较高，本专利技术中使用Keystone变换方法，该方法与回波信噪比无关，适用于低信噪比时的距离包络对齐，其本质是在对信号回波数据采样时进行尺度变换，对慢时间轴的伸缩变化、伸缩幅度与频率相关。在高频时进行拉伸，幅度较大，低频时进行压缩，幅度较小。Keystone变换能将不同距离单元上的回波校正到同一距离单元上，实现了距离走动的校正和包络补偿。针对多普勒频率徙动补偿问题，目前有解线调法，多项式Wigner-vile分布(WVD)和高阶模糊函数法(HAF)和Wigner-Hough变换等。解线调法将多普勒补偿和MTD分成了互不相关的两步，在运动杂波背景下容易造成误补偿，多项式Wigner-vile分布(WVD)和高阶模糊函数法采用非线性变换存在交叉项干扰问题，Wigner-Hough与FRFT相比，积累增益达不到相参积累水平。本专利技术中使用分数阶傅里叶变换(fractional Fourier transform, FRFT),可以达到相参积累增益的水平且没有交叉项干扰问题。由于距离徙动和多普勒频率徙动均同时得到了补偿，使得分散在多个距离单元和多个多普勒频率单元的能量聚集在一个距离单元和多普勒频率单元中。因此，积累时间不再受徙动的制约，可以通过增加积累脉冲数，增加积累增益，提高对高速加速运动微弱目标的检测能力。
技术实现思路
本专利技术针对现有的雷达机动弱目标检测中同时存在的距离徙动、多普勒频率徙动、信噪比不足、补偿效果差的问题，提出了。本专利技术，具体包括以下步骤:步骤1、建立机动弱目标回波信号模型S(t’，tm)。(I)、雷达发射线性调频脉冲，其数学表达式为:【权利要求】1.，其特征在于该方法包括以下步骤: 步骤1、建立机动弱目标回波信号模型s(t’，tm)； (1)、雷达发射线性调频脉冲，其数学表达式为: 【文档编号】G01S13/50GK103675759SQ201310627133【公开日】2014年3月26日 申请日期:2013年11月27日 优先权日:2013年11月27日 【专利技术者】郭云飞, 周森山 申请人:杭州电子科技大学本文档来自技高网...

【技术保护点】
一种改进的分数阶傅里叶变换机动弱目标检测方法，其特征在于该方法包括以下步骤：步骤1、建立机动弱目标回波信号模型S(t',tm)；(1)、雷达发射线性调频脉冲，其数学表达式为： p ( t ) = rect ( t T p ) exp ( - jπγ t 2 ) - - - ( 1 ) 其中，u为变量，当时，rect(u)＝1，当u为其他时，rect(u)＝0，t为时间，B为脉冲带宽，Tp为脉冲时宽；调频率(2)、则雷达发射一组脉冲串信号可以表示为v(t‑mTr)＝p(t‑mTr)exp(‑j2πfc(t‑mTr))     (2)其中m为发射脉冲个数，Tr为脉冲重复周期，fc为载频；(3)、假设脉冲串照射到一个距离为R的弱目标上，接受回波信号记录在二维数组S(t',tm)；其中t'＝(t‑mTr)为距离维，tm＝mTr为脉冲维，m=0,1,…,M,令T＝MTr,则则信号表达式为 S ( t ′ , t m ) = Ap ( t ′ - 2 R ( t m ) c ) exp ( - j 4 π f c c ) R ( t m ) - - - ( 3 ) 其中A为常数，取决与弱目标雷达散射截面积的大小，c为电磁波传播速度，目标距离 R ( t m ) = R 0 + vt m + 1 2 at m 2 ; 其中R0为雷达到目标的初始距离，v为速度，a为目标的加速度；步骤2、对回波原始数据在快时间域进行傅里叶变换，进行匹配压缩得到S(f,tm)；对回波原始数据在快时间域进行傅里叶变换,进行匹配滤波，得 到信号数据S(f,tm)，其表达式为 S ( f , t m ) = AB γ sin c ( B ( t ′ - 2 R ( t m ) c ) ) exp ( ...

【技术特征摘要】

【专利技术属性】
技术研发人员：郭云飞，周森山，
申请(专利权)人：杭州电子科技大学，
类型：发明
国别省市：浙江;33