CN104849708B - 基于频域多项式相位变换的高速机动目标参数估计方法 - Google Patents
基于频域多项式相位变换的高速机动目标参数估计方法 Download PDFInfo
- Publication number
- CN104849708B CN104849708B CN201510251019.8A CN201510251019A CN104849708B CN 104849708 B CN104849708 B CN 104849708B CN 201510251019 A CN201510251019 A CN 201510251019A CN 104849708 B CN104849708 B CN 104849708B
- Authority
- CN
- China
- Prior art keywords
- frequency domain
- data
- slow time
- apart
- range cell
- 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
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
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
一种基于频域多项式相位变换的高速机动目标参数估计方法。其是在天基雷达接收到的数据基础上对高速机动目标回波数据进行脉冲压缩处理,沿经过脉冲压缩处理后数据的快时间向进行快速傅里叶变换以得到距离频域‑慢时间数据,对距离频域‑慢时间数据进行延时,把经过延时的距离频域‑慢时间数据的共轭和未经过延时的距离频域‑慢时间数据进行相乘处理以实现频域多项式相位变换,接着沿经过频域多项式相位变换处理后数据的距离频域向进行逆快速傅里叶变换,再在此基础上沿着慢时间向进行快速傅里叶变换以得到距离单元‑多普勒数据,最后利用距离单元‑多普勒数据进行目标参数估计。本发明能在低信噪比和积累脉冲数有限的情况下获得较高精度的目标参数估计结果。
Description
技术领域
本发明属于雷达信号处理和空间高速机动目标参数估计技术领域,特别是涉及一种基于频域多项式相位变换的高速机动目标参数估计方法。
背景技术
随着世界各军事强国对于空间信息的重视程度不断加深,对空间高速机动目标的监视已成为各国急需发展的高精尖技术。空间高速机动目标主要包括卫星、空间碎片等。由于这些空间高速机动目标的运动速度较快且带有加速度,如空间碎片一般能够达到7000m/s的速度,加速度可以达到20g,因此在空间中运行的飞行器很容易受到其侵袭,从而造成极大伤害。由于空间高速机动目标的运动速度快,并且机动能力强,回波信号微弱,从而在低信噪比下对空间高速机动目标进行参数估计而为后续目标识别就会带来较大困难。
目前针对目标参数估计的方法有Radon-Fourier变换方法(RFT)、分数阶Fourier变换(FRFT)、Radon-FRFT变换方法等。其中Radon-Fourier变换方法只能对匀速直线运动目标进行能量积累,且只能对目标进行速度估计。分数阶Fourier变换(FRFT)方法是针对没有距离徙动的机动目标回波信号进行变换处理,它将机动目标回波的线性调频信号在时频域内进行旋转然后进行相应的Fourier变换,即完成相参积累,但是该方法在脉冲数有限和信号能量分散到不同距离单元时很难获得较好的参数估计效果。同样Radon-FRFT变换方法在沿着机动目标参数轨迹取出数据时只能将整数单元的数据取出,由于机动目标距离走动量往往不是整数个距离单元,而是会出现走动分数个距离单元,故Radon-FRFT方法不能够较好地将机动目标数据取出并进行相应的分数阶Fourier变换处理,同样也不能直接用于空间高速机动目标的参数估计。
发明内容
为了解决上述问题,本发明的目的在于提供一种基于频域多项式相位变换的高速机动目标参数估计方法。
为了达到上述目的,本发明提供的基于频域多项式相位变换的高速机动目标参数估计方法包括按顺序进行的下列步骤:
1)将天基雷达接收到的高速机动目标回波数据进行脉冲压缩处理;
2)沿上述经过脉冲压缩处理后数据的距离单元向进行快速傅里叶变换处理,得到距离频域-慢时间数据并留一个相同的备份数据,接着沿备份的距离频域-慢时间数据的慢时间向进行延时处理;
3)将上述未经过延时和经过延时的距离频域-慢时间数据的共轭进行相乘以得到带有延时量的距离频域-慢时间数据并留一个备份数据;然后沿带有延时量的距离频域-慢时间数据的距离频域向进行逆快速傅里叶变换以得到经过频域多项式相位变换后的距离单元-慢时间数据;
4)沿步骤3)中得到的经过频域多项式相位变换后得到的距离单元-慢时间数据的慢时间向进行快速傅里叶变换,以得到用于目标加速度估计的距离单元-多普勒数据;接着将得到的距离单元-多普勒数据中的最大值求出,并把最大值所对应的距离单元和多普勒频率求出,利用求出的多普勒频率值估计得到目标加速度值;
5)利用步骤4)中得到的最大值所在的距离单元求出目标的速度范围,在速度范围内进行搜索,根据搜索结果和加速度估计值构造距离频域补偿项,并将补偿项与步骤3)中得到的带有延时量的距离频域-慢时间备份数据相乘;将乘积先沿距离频域进行逆快速傅里叶变换,再沿慢时间向进行快速傅里叶变换以得到用于速度估计的距离单元-多普勒数据;
6)将步骤5)中得到的距离单元-多普勒数据中第一个距离单元内数据的最大值求出,然后将最大值对应的速度搜索值取出并一同放到速度估计信息矩阵中;重复步骤5)至步骤6)操作直到目标速度范围搜索完为止;最后将速度估计信息矩阵存储的所有最大值中的峰值求出,并将峰值对应的速度搜索值作为目标速度的估计结果。
在步骤2)中,所述的沿经过脉冲压缩处理后数据的距离单元向进行快速傅里叶变换处理,得到距离频域-慢时间数据并留一个相同的备份数据,接着沿备份的距离频域-慢时间数据的慢时间向进行延时处理的方法是:
设天基雷达接收到的经过脉冲压缩处理后的信号表达式为:
其中A1为回波信号复幅度,p(·)为回波包络,t'为快时间,tm=mTr为慢时间,c为光速,为信号波长,M为相干积累脉冲数,m=1,...,M,Tr为脉冲重复周期,fc为载波频率,为空间高速机动目标与天基雷达之间的距离,R0为目标的初始距离,v0为目标运动速度,a0为目标加速度;沿式(1)中的快时间向进行快速傅里叶变换得到:
f为距离频域,为调频率,TP为脉冲时宽;留一个与S(f,tm)相同的备份数据,记为S'(f,tm),沿S'(f,tm)的慢时间向进行延时,延时量为τ,可以得到:
在步骤3)中,所述的将未经过延时和经过延时的距离频域-慢时间数据的共轭进行相乘以得到带有延时量的距离频域-慢时间数据并留一个备份数据;然后沿带有延时量的距离频域-慢时间数据的距离频域向进行逆快速傅里叶变换以得到经过频域多项式相位变换后的距离单元-慢时间数据的方法是:
将式(2)与式(3)的共轭相乘以实现频域多项式相位变换:
其中“*”表示共轭处理,接着留一个和X(f,tm;τ)同样的备份数据并记为X'(f,tm;τ);
沿X(f,tm;τ)的距离频域进行逆快速傅里叶变换以得到经过频域多项式相位变换处理后的距离单元-慢时间数据:
其中A2是经过频域多项式相位变换后数据的复幅度。
在步骤4)中,所述的沿步骤3)中得到的经过频域多项式相位变换后得到的距离单元-慢时间数据的慢时间向进行快速傅里叶变换,以得到用于目标加速度估计的距离单元-多普勒数据;接着将得到的距离单元-多普勒数据中的最大值求出,并把最大值所对应的距离单元和多普勒频率求出,利用求出的多普勒频率值估计得到目标加速度值的方法是:
先沿步骤3)中得到的s'r(t',tm)的慢时间进行快速傅里叶变换处理可以得到用于目标加速度估计的距离单元-多普勒数据,把得到的距离单元-多普勒数据中的最大值求出来,并且相应地得到最大值所在距离单元和多普勒单元位置,根据多普勒单元位置得到相应的多普勒频率估计值,从而由得到目标加速度的估计值,其中是估计得到的多普勒频率。
在步骤5)中,所述的利用步骤4)中得到的最大值所在的距离单元求出目标的速度范围,在速度范围内进行搜索,根据搜索结果和加速度估计值构造距离频域补偿项,并将补偿项与步骤3)中得到的带有延时量的距离频域-慢时间备份数据相乘;将乘积先沿距离频域进行逆快速傅里叶变换,再沿慢时间向进行快速傅里叶变换以得到用于速度估计的距离单元-多普勒数据的方法是:
根据步骤4)中得到的最大值所在距离单元并根据得到目标速度的粗估计值,其中△r为天基雷达的距离分辨率,再根据这个粗估计值能够得到目标速度所在的范围;在得到的目标速度范围内进行搜索,步长为将速度范围内的搜索值个数记为N,利用搜索得到的速度值和步骤4)中得到的加速度估计值构造两个相位项:和其中v表示速度搜索值,tm=mTr,m=1,...,M,用这两个相位项一同乘以步骤3)中得到的备份数据X'(f,tm;τ),将乘积记为X”(f,tm;τ),沿X”(f,tm;τ)的距离频域维进行逆快速傅里叶变换,然后再沿慢时间向进行快速傅里叶变换以得到用于估计目标速度的距离单元-多普勒数据。
在步骤6)中,所述的将步骤5)中得到的距离单元-多普勒数据中第一个距离单元内数据的最大值求出,然后将最大值对应的速度搜索值取出并一同放到速度估计信息矩阵中;重复步骤5)至步骤6)操作直到目标速度范围搜索完为止;最后将速度估计信息矩阵存储的所有最大值中的峰值求出,并将峰值对应的速度搜索值作为目标速度的估计结果的方法是:
将步骤5)中得到的X”(f,tm;τ)将应用于估计目标速度的距离单元-多普勒数据第一个距离单元内的最大值求出,并将最大值对应的速度搜索值记录下来,将它们放入速度估计信息矩阵G2×N中,其中第一行放最大值,第二行放最大值对应的速度搜索值;重复步骤5)至步骤6),直到目标速度范围搜索并经过相应处理完为止;最后把G2×N中第一行中所有最大值中的峰值求出来,并把峰值对应的速度搜索值作为目标速度的估计结果。
本发明提供的基于频域多项式相位变换的高速机动目标参数估计方法是在天基雷达接收到的数据基础上对机动目标回波数据进行脉冲压缩处理,沿经过脉冲压缩处理后的快时间向进行快速傅里叶变换以得到距离频域-慢时间数据,对距离频域-慢时间数据进行延时,把经过延时的距离频域-慢时间数据的共轭和未经过延时的距离频域-慢时间数据进行相乘处理,接着沿经过频域多项式相位变换后数据的距离频域向进行逆快速傅里叶变换,然后再在此基础上沿着慢时间向进行快速傅里叶变换以得到距离单元-多普勒数据,最后利用距离单元-多普勒数据进行目标检测参数估计。本发明方法处理过程中使用的都是快速傅里叶变换和逆快速傅里叶变换,在对目标速度估计中也只进行速度参数的一维搜索。将仿真结果与常规RFT、FRFT和Radon-FRFT方法进行对比可知,本发明方法能在低信噪比和积累脉冲数有限的情况下获得较高精度的目标参数估计结果,从而验证了本发明方法的有效性。
附图说明
图1为本发明提供的基于频域多项式相位变换的高速机动目标参数估计方法流程图。
图2为经过解调和脉冲压缩处理后的天基雷达接收数据分布图。
图3为对接收到的数据进行延时后的结果图。
图4为对接收到的回波数据进行频域多项式相位变换后的结果图。
图5为应用本发明方法、RFRFT方法、FRFT方法、RFT方法对目标速度进行估计的均方根误差随信噪比变化曲线图。
图6为应用本发明方法、RFRFT方法、FRFT方法对目标加速度进行估计的均方根误差随信噪比变化曲线图。
具体实施方式
下面结合附图和具体实施例对本发明提供的基于频域多项式相位变换的高速机动目标参数估计方法进行详细说明。
如图1所示,本发明提供的基于频域多项式相位变换的高速机动目标参数估计方法包括按顺序进行的下列步骤:
1)将天基雷达接收到的高速机动目标回波数据进行脉冲压缩处理;
设天基雷达接收到的未经脉冲压缩处理的距离频域-慢时间域信号形式为:
其中A1为信号幅度,Bs为信号带宽,TP为脉冲信号宽度,c为光速,为调频率,fc表示信号的载波频率,f表示距离频域(即本文中所指的频域),tm=mTr为慢时间,且m=1,...M,为空间高速机动目标与天基雷达间的距离,且有:
对式(1)进行匹配滤波,即乘以的共轭,沿相乘后的距离频域-慢时间回波数据的距离频域向进行逆快速傅里叶变换后得到经过脉冲压缩处理后的回波信号为:
式(3)中p(·)为回波包络,t'为快时间。
2)沿上述经过脉冲压缩处理后数据的距离单元向进行快速傅里叶变换处理,得到距离频域-慢时间数据并留一个相同的备份数据,接着沿备份的距离频域-慢时间数据的慢时间向进行延时处理;
在步骤2)中,对步骤1)中式(3)的距离单元向进行快速傅里叶变换后得到的距离频域-慢时间数据表达式为:
对距离频域-慢时间数据留一个与S(f,tm)相同的备份数据,记为S'(f,tm);沿S'(f,tm)的慢时间向进行时延(时延量τ为0.01秒)后得到:
3)将上述未经过延时和经过延时的距离频域-慢时间数据的共轭进行相乘以得到带有延时量的距离频域-慢时间数据并留一个备份数据;然后沿带有延时量的距离频域-慢时间数据的距离频域向进行逆快速傅里叶变换以得到经过频域多项式相位变换后的距离单元-慢时间数据;
在步骤3)中,将未经延时和经过延时的距离频域-慢时间数据的共轭进行相乘处理以实现频域多项式相位变换处理,即将式(4)和式(5)的共轭进行相乘:
其中“*”表示共轭处理,留一个和X(f,tm;τ)同样的备份数据并记为X'(f,tm;τ)。
沿X(f,tm;τ)的距离频域进行逆快速傅里叶变换以得到经过频域多项式相位变换处理后的距离单元-慢时间数据:
其中A2是经过频域多项式相位变换后数据的复幅度。
4)沿步骤3)中得到的经过频域多项式相位变换后得到的距离单元-慢时间数据的慢时间向进行快速傅里叶变换,以得到用于目标加速度估计的距离单元-多普勒数据;接着将得到的距离单元-多普勒数据中的最大值求出,并把最大值所对应的距离单元和多普勒频率求出,利用求出的多普勒频率值估计得到目标加速度值;
在步骤4)中,先沿步骤3)中得到的距离单元-慢时间数据s'r(t',tm)的慢时间向进行快速傅里叶变换处理以得到用于目标加速度估计的距离单元-多普勒数据,然后将距离单元-多普勒数据域中的最大值求出来,并且相应地得到最大值所在距离单元和多普勒单元位置,根据多普勒单元位置得到相应的多普勒频率估计值,从而由(是估计得到的多普勒频率)得到目标加速度的估计值。
5)利用步骤4)中得到的最大值所在的距离单元求出目标的速度范围,在速度范围内进行搜索,根据搜索结果和加速度估计值构造距离频域补偿项,并将补偿项与步骤3)中得到的带有延时量的距离频域-慢时间备份数据相乘;将乘积先沿距离频域进行逆快速傅里叶变换,再沿慢时间向进行快速傅里叶变换以得到用于速度估计的距离单元-多普勒数据;
在步骤5)中,记步骤4)中得到的最大值所在的距离单元为则根据(△r是雷达的距离分辨率)得到目标速度的粗估计值,根据这个粗估计值能够得到目标速度所在的范围。在得到的目标速度范围内进行搜索(步长为将速度范围内的搜索值个数记为N),利用搜索得到的速度值和步骤4)中得到的加速度估计值构造两个相位项:和(其中v表示速度搜索值,tm=mTr,m=1,...,M),用这两个相位项一同乘以步骤3)中得到的备份数据X'(f,tm;τ),将乘积记为X”(f,tm;τ),沿X”(f,tm;τ)的距离频域维进行逆快速傅里叶变换,然后再沿慢时间向进行快速傅里叶变换以得到X”(f,tm;τ)对应的用于估计目标速度的距离单元-多普勒数据(如果两个相位项:和中的参数和目标参数一致时,即为和则用这两个相位项对X'(f,tm;τ)进行补偿后,目标回波数据就位于第一个距离单元内,此时对第一个距离单元内的数据进行快速傅里叶变换能够获得对目标能量的最大积累,故可用这一性质对目标速度进行估计)。
6)将步骤5)中得到的距离单元-多普勒数据中第一个距离单元内数据的最大值求出,然后将最大值对应的速度搜索值取出并一同放到速度估计信息矩阵中;重复步骤5)至步骤6)操作直到目标速度范围搜索完为止;最后将速度估计信息矩阵存储的所有最大值中的峰值求出,并将峰值对应的速度搜索值作为目标速度的估计结果。
在步骤6)中,将步骤5)中得到的X”(f,tm;τ)将应用于估计目标速度的距离单元-多普勒数据第一个距离单元内的最大值求出,并将最大值对应的速度搜索值记录下来,将它们放入速度估计信息矩阵G2×N中,其中第一行放最大值,第二行放最大值对应的速度搜索值;重复步骤5)至步骤6),直到目标速度范围搜索并经过相应处理完为止,此时G2×N中的第一行中存储了每个速度搜索值并完成相应处理后得到的距离单元-慢时间数据中第一个距离单元内的最大值,第二行存储了每个最大值对应的速度搜索值。最后把G2×N中第一行中所有最大值中的峰值求出来,并把峰值对应的速度搜索值作为目标速度的估计结果。
仿真结果及分析:
本发明提供的基于频域多项式相位变换的高速机动目标参数估计方法的效果可以通过以下仿真实验进一步说明。仿真参数设置:天基雷达发射信号波长为0.3m,脉冲重复频率为1500Hz,信号带宽为30MHz,载频为1GHz,输入信噪比为-10dB。目标参数设置为:目标初始在第200个距离单元中,速度v0=3000m/s,加速度a0=180m/s2。
1、回波数据分布图
图2为经过解调和脉冲压缩处理后的天基雷达接收数据分布图。从图2可以看出,经过脉冲压缩处理后,回波信号在距离单元-慢时间内会出现距离徙动现象,回波能量会分散在不同的距离单元内,从而不利于对目标参数进行估计。
2、接收数据进行延时处理
图3为对接收到的数据进行延时后的结果图,其中延时量为0.01秒。从图中可以看出接收数据经过延时后仍然存在距离徙动现象。
3、对回波数据进行频域多项式相位变换后的结果
图4为对接收到的回波数据进行频域多项式相位变换后的结果图。从图中可以看出,经过频域多项式相位变换后,回波数据校正到了一个距离单元内,此时能够将校正到一个距离单元内的回波数据进行能量积累以用于后续的目标参数估计。
4、目标参数估计性能对比
图5为采用本发明方法、RFRFT方法、FRFT方法、RFT方法对目标速度进行估计的均方根误差随信噪比变化曲线图。图6为采用本发明方法、RFRFT方法和FRFT方法对目标加速度进行估计的均方根误差随信噪比变化曲线图。从图5可以看出,各个方法对目标速度估计的均方根误差随信噪比增加都呈现出逐渐降低的趋势,但本发明方法能够在低信噪比情况下迅速向CRB界靠拢,并且能够更加接近CRB界,对于速度的估计性能优于其他方法。从图6也能够看出本发明方法对目标加速度的估计均方根误差在低信噪比下能迅速向CRB界靠拢,且能够更接近加速度估计的CRB界,估计性能优于RFRFT方法和FRFT方法。
Claims (6)
1.一种基于频域多项式相位变换的高速机动目标参数估计方法,其特征在于,所述的方法包括按顺序进行下列的步骤:
1)将天基雷达接收到的高速机动目标回波数据进行脉冲压缩处理;
2)沿上述经过脉冲压缩处理后数据的距离单元向进行快速傅里叶变换处理,得到距离频域-慢时间数据并留一个相同的备份数据,接着沿备份的距离频域-慢时间数据的慢时间向进行延时处理;
3)将上述未经过延时和经过延时的距离频域-慢时间数据的共轭进行相乘以得到带有延时量的距离频域-慢时间数据并留一个备份数据;然后沿带有延时量的距离频域-慢时间数据的距离频域向进行逆快速傅里叶变换以得到经过频域多项式相位变换后的距离单元-慢时间数据;
4)沿步骤3)中得到的经过频域多项式相位变换后得到的距离单元-慢时间数据的慢时间向进行快速傅里叶变换,以得到用于目标加速度估计的距离单元-多普勒数据;接着将得到的距离单元-多普勒数据中的最大值求出,并把最大值所对应的距离单元和多普勒频率求出,利用求出的多普勒频率值估计得到目标加速度值;
5)利用步骤4)中得到的最大值所在的距离单元求出目标的速度范围,在速度范围内进行搜索,根据搜索结果和加速度估计值构造距离频域补偿项,并将补偿项与步骤3)中得到的带有延时量的距离频域-慢时间备份数据相乘;将乘积先沿距离频域进行逆快速傅里叶变换,再沿慢时间向进行快速傅里叶变换以得到用于速度估计的距离单元-多普勒数据;
6)将步骤5)中得到的距离单元-多普勒数据中第一个距离单元内数据的最大值求出,然后将最大值对应的速度搜索值取出并一同放到速度估计信息矩阵中;重复步骤5)至步骤6)操作直到目标速度范围搜索完为止;最后将速度估计信息矩阵存储的所有最大值中的峰值求出,并将峰值对应的速度搜索值作为目标速度的估计结果。
2.根据权利要求1所述的基于频域多项式相位变换的高速机动目标参数估计方法,其特征在于:在步骤2)中,所述的沿经过脉冲压缩处理后数据的距离单元向进行快速傅里叶变换处理,得到距离频域-慢时间数据并留一个相同的备份数据,接着沿备份的距离频域-慢时间数据的慢时间向进行延时处理的方法是:
设天基雷达接收到的经过脉冲压缩处理后的信号表达式为:
其中A1为回波信号复幅度,p(·)为回波包络,t'为快时间,tm=mTr为慢时间,c为光速,为信号波长,M为相干积累脉冲数,m=1,...,M,Tr为脉冲重复周期,fc为载波频率,为空间高速机动目标与天基雷达之间的距离,R0为目标的初始距离,v0为目标运动速度,a0为目标加速度;沿式(1)中的快时间向进行快速傅里叶变换得到:
f为距离频域,为调频率,TP为脉冲时宽;留一个与S(f,tm)相同的备份数据,记为S'(f,tm),沿S'(f,tm)的慢时间向进行延时,延时量为τ,可以得到:
3.根据权利要求1所述的基于频域多项式相位变换的高速机动目标参数估计方法,其特征在于:在步骤3)中,所述的将未经过延时和经过延时的距离频域-慢时间数据的共轭进行相乘以得到带有延时量的距离频域-慢时间数据并留一个备份数据;然后沿带有延时量的距离频域-慢时间数据的距离频域向进行逆快速傅里叶变换以得到经过频域多项式相位变换后的距离单元-慢时间数据的方法是:
将式(2)与式(3)的共轭相乘以实现频域多项式相位变换:
其中“*”表示共轭处理,接着留一个和X(f,tm;τ)同样的备份数据并记为X'(f,tm;τ);
沿X(f,tm;τ)的距离频域进行逆快速傅里叶变换以得到经过频域多项式相位变换处理后的距离单元-慢时间数据:
其中A2是经过频域多项式相位变换后数据的复幅度。
4.根据权利要求1所述的基于频域多项式相位变换的高速机动目标参数估计方法,其特征在于:在步骤4)中,所述的沿步骤3)中得到的经过频域多项式相位变换后得到的距离单元-慢时间数据的慢时间向进行快速傅里叶变换,以得到用于目标加速度估计的距离单元-多普勒数据;接着将得到的距离单元-多普勒数据中的最大值求出,并把最大值所对应的距离单元和多普勒频率求出,利用求出的多普勒频率值估计得到目标加速度值的方法是:
先沿步骤3)中得到的s'r(t',tm)的慢时间进行快速傅里叶变换处理可以得到用于目标加速度估计的距离单元-多普勒数据,把得到的距离单元-多普勒数据中的最大值求出来,并且相应地得到最大值所在距离单元和多普勒单元位置,根据多普勒单元位置得到相应的多普勒频率估计值,从而由得到目标加速度的估计值,其中是估计得到的多普勒频率。
5.根据权利要求1所述的基于频域多项式相位变换的高速机动目标参数估计方法,其特征在于:在步骤5)中,所述的利用步骤4)中得到的最大值所在的距离单元求出目标的速度范围,在速度范围内进行搜索,根据搜索结果和加速度估计值构造距离频域补偿项,并将补偿项与步骤3)中得到的带有延时量的距离频域-慢时间备份数据相乘;将乘积先沿距离频域进行逆快速傅里叶变换,再沿慢时间向进行快速傅里叶变换以得到用于速度估计的距离单元-多普勒数据的方法是:
根据步骤4)中得到的最大值所在距离单元并根据得到目标速度的粗估计值,其中Δr为天基雷达的距离分辨率,再根据这个粗估计值能够得到目标速度所在的范围;在得到的目标速度范围内进行搜索,步长为将速度范围内的搜索值个数记为N,利用搜索得到的速度值和步骤4)中得到的加速度估计值构造两个相位项:和其中v表示速度搜索值,tm=mTr,m=1,...,M,用这两个相位项一同乘以步骤3)中得到的备份数据X'(f,tm;τ),将乘积记为X”(f,tm;τ),沿X”(f,tm;τ)的距离频域维进行逆快速傅里叶变换,然后再沿慢时间向进行快速傅里叶变换以得到用于估计目标速度的距离单元-多普勒数据。
6.根据权利要求1所述的一种基于频域多项式相位变换的高速机动目标参数估计方法,其特征在于:在步骤6)中,所述的将步骤5)中得到的距离单元-多普勒数据中第一个距离单元内数据的最大值求出,然后将最大值对应的速度搜索值取出并一同放到速度估计信息矩阵中;重复步骤5)至步骤6)操作直到目标速度范围搜索完为止;最后将速度估计信息矩阵存储的所有最大值中的峰值求出,并将峰值对应的速度搜索值作为目标速度的估计结果的方法是:
将步骤5)中得到的X”(f,tm;τ)将应用于估计目标速度的距离单元-多普勒数据第一个距离单元内的最大值求出,并将最大值对应的速度搜索值记录下来,将它们放入速度估计信息矩阵G2×N中,其中第一行放最大值,第二行放最大值对应的速度搜索值;重复步骤5)至步骤6),直到目标速度范围搜索并经过相应处理完为止;最后把G2×N中第一行中所有最大值中的峰值求出来,并把峰值对应的速度搜索值作为目标速度的估计结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510251019.8A CN104849708B (zh) | 2015-05-18 | 2015-05-18 | 基于频域多项式相位变换的高速机动目标参数估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510251019.8A CN104849708B (zh) | 2015-05-18 | 2015-05-18 | 基于频域多项式相位变换的高速机动目标参数估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104849708A CN104849708A (zh) | 2015-08-19 |
CN104849708B true CN104849708B (zh) | 2017-03-08 |
Family
ID=53849492
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510251019.8A Expired - Fee Related CN104849708B (zh) | 2015-05-18 | 2015-05-18 | 基于频域多项式相位变换的高速机动目标参数估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104849708B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108226929B (zh) * | 2018-01-12 | 2019-12-31 | 北京航空航天大学 | 一种正侧视sar慢速目标的检测方法及检测系统 |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104880698B (zh) * | 2015-06-05 | 2017-03-08 | 中国民航大学 | 基于距离频域多项式相位变换的空间机动目标检测方法 |
CN106154241B (zh) * | 2016-06-21 | 2018-12-28 | 大连大学 | 脉冲噪声环境下的韧性平行因子分析算法 |
CN106154242B (zh) * | 2016-06-21 | 2018-11-20 | 大连大学 | 基于分数低阶类相关熵的目标参数联合估计算法 |
CN107132534B (zh) * | 2017-06-27 | 2020-02-18 | 西安电子科技大学 | 一种高速雷达目标频域检测的优化方法 |
US11899095B2 (en) | 2020-09-17 | 2024-02-13 | Texas Instruments Incorporated | Doppler processing in frequency-modulated continuous wave radar systems using dither |
CN113884999A (zh) * | 2021-09-27 | 2022-01-04 | 中国人民解放军32802部队 | 一种外辐射源雷达的跨距离单元补偿方法 |
CN116482631A (zh) * | 2023-06-26 | 2023-07-25 | 四川天中星航空科技有限公司 | 一种基于雷达回波信号的目标模拟方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2251849A2 (de) * | 2009-05-09 | 2010-11-17 | EADS Deutschland GmbH | Verfahren zur verbesserten Erkennung von leitungsartigen Objekten |
CN103412295A (zh) * | 2013-08-30 | 2013-11-27 | 西安电子科技大学 | 基于回波精确模型的高速机动弱目标检测方法 |
CN104200072A (zh) * | 2014-08-18 | 2014-12-10 | 中国民航大学 | 基于改进三次相位变换的高机动目标检测和参数估计方法 |
-
2015
- 2015-05-18 CN CN201510251019.8A patent/CN104849708B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2251849A2 (de) * | 2009-05-09 | 2010-11-17 | EADS Deutschland GmbH | Verfahren zur verbesserten Erkennung von leitungsartigen Objekten |
CN103412295A (zh) * | 2013-08-30 | 2013-11-27 | 西安电子科技大学 | 基于回波精确模型的高速机动弱目标检测方法 |
CN104200072A (zh) * | 2014-08-18 | 2014-12-10 | 中国民航大学 | 基于改进三次相位变换的高机动目标检测和参数估计方法 |
Non-Patent Citations (3)
Title |
---|
Cross-range scaling method of inverse synthetic aperture radar image based on discrete polynomial-phase transform;Lei Liu et al.;《IET Radar,Sonar & Navigation》;20150305;第333-341页 * |
一种天波雷达机动目标运动补偿算法;游伟 等;《现代雷达》;20131115;第35卷(第11期);第26-29页 * |
基于离散多项式相位变换和分数阶傅里叶变换的加速目标检测算法;庞存锁;《电子学报》;20120115(第1期);第184-188页 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108226929B (zh) * | 2018-01-12 | 2019-12-31 | 北京航空航天大学 | 一种正侧视sar慢速目标的检测方法及检测系统 |
Also Published As
Publication number | Publication date |
---|---|
CN104849708A (zh) | 2015-08-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104849708B (zh) | 基于频域多项式相位变换的高速机动目标参数估计方法 | |
CN104076351B (zh) | 一种用于高速高机动目标的相参积累检测方法 | |
CN104237883B (zh) | 一种采用稀疏表示的机载雷达空时自适应处理方法 | |
CN103744068B (zh) | 双通道调频连续波sar系统的动目标检测成像方法 | |
CN106872954B (zh) | 一种高超声速平台杂波抑制和运动目标成像方法 | |
CN105137410B (zh) | 基于ofdm的高分辨率雷达通信一体化的波形优化方法 | |
CN103399310A (zh) | 基于相位差分Radon-Lv分布的雷达微弱动目标检测方法 | |
CN105738879A (zh) | 基于稀疏恢复的雷达杂波空时自适应预滤波方法 | |
CN106772304B (zh) | 基于空域多级分解的机载mimo雷达后多普勒自适应处理方法 | |
CN104977571B (zh) | 基于俯仰频率分集stap的距离模糊杂波抑制方法 | |
CN108008386B (zh) | 一种基于单快拍music算法的距离向处理方法 | |
CN104597434B (zh) | 改进包络移位补偿和分数阶傅里叶变换的多帧相参tbd方法 | |
CN104502898B (zh) | 将修正rft和修正mdcft相结合的机动目标参数估计方法 | |
Li et al. | ISAR imaging of nonuniformly rotating target based on the multicomponent CPS model under low SNR environment | |
CN104597435B (zh) | 基于修正频域补偿和分数阶傅里叶变换的多帧相参tbd方法 | |
CN104793194B (zh) | 基于改进的自适应多脉冲压缩的距离‑多普勒估计方法 | |
CN102621536B (zh) | 基于relax的空中多机动目标检测与参数估计方法 | |
Fang et al. | Radar maneuvering target detection based on two steps scaling and fractional Fourier transform | |
CN104950295A (zh) | 一种基于相关函数与尺度变换的高速机动目标检测方法 | |
CN107367715A (zh) | 基于稀疏表示的杂波抑制方法 | |
CN104345301A (zh) | 机载mimo雷达非自适应杂波预滤波空时二维对消方法 | |
CN104375128B (zh) | 一种基于互相关函数的高机动目标快速积累检测方法 | |
CN103983959A (zh) | 基于数据重构的sar系统运动目标径向速度估计方法 | |
CN111045002A (zh) | 基于TRT和SNuFFT的机动目标相参积累方法 | |
CN103760540B (zh) | 基于重构信号和1-范数的动目标检测与参数估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170308 Termination date: 20170518 |
|
CF01 | Termination of patent right due to non-payment of annual fee |