CN105301590B - 一种机动目标调频步进逆合成孔径成像方法 - Google Patents
一种机动目标调频步进逆合成孔径成像方法 Download PDFInfo
- Publication number
- CN105301590B CN105301590B CN201510738038.3A CN201510738038A CN105301590B CN 105301590 B CN105301590 B CN 105301590B CN 201510738038 A CN201510738038 A CN 201510738038A CN 105301590 B CN105301590 B CN 105301590B
- Authority
- CN
- China
- Prior art keywords
- mrow
- target
- msub
- mover
- 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.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 27
- 238000012545 processing Methods 0.000 claims abstract description 8
- 230000001133 acceleration Effects 0.000 claims description 33
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 238000005311 autocorrelation function Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 210000001367 artery Anatomy 0.000 claims 1
- 238000011156 evaluation Methods 0.000 claims 1
- 210000003462 vein Anatomy 0.000 claims 1
- 238000000034 method Methods 0.000 abstract description 19
- 230000001427 coherent effect Effects 0.000 abstract description 8
- 230000005012 migration Effects 0.000 abstract description 7
- 238000013508 migration Methods 0.000 abstract description 7
- 238000009825 accumulation Methods 0.000 abstract description 6
- 238000010200 validation analysis Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- JEIPFZHSYJVQDO-UHFFFAOYSA-N ferric oxide Chemical compound O=[Fe]O[Fe]=O JEIPFZHSYJVQDO-UHFFFAOYSA-N 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 238000003786 synthesis reaction Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 238000012935 Averaging Methods 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000000844 transformation 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
- 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
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9064—Inverse SAR [ISAR]
-
- 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
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- 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)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供了一种机动目标调频步进逆合成孔径成像方法，属于雷达信号处理和成像技术领域。该方法包括以下步骤：1)雷达回波经过去斜处理后得到粗距离像；2)从粗距离像信号FRAF变换谱图中估计目标运动参数；3)利用不同脉冲串的参数估计值拟合出精确值；4)构造相位补偿函数，补偿粗距离像信号；5)对补偿后的粗距离像信号在慢时间维FFT得到目标精距离像；6)对目标精距离像信号进行相位聚焦；7)进行转台成像，得到ISAR图像。该方法在低信噪比下可稳健估计目标运动参数，能补偿目标运动产生的距离和多普勒徙动，实现目标能量的相参积累，提升了调频步进雷达在低信噪比下对机动目标的成像能力，具有推广应用价值。
Description
技术领域
本发明涉及一种低信噪比下机动目标的调频步进雷达成像方法，该方法可以应用于调频步进雷达在低信噪比下对机动目标进行成像。
背景技术
调频步进信号是一种易于实现的距离高分辨信号，它发射一串载频线性跳变的窄带调频脉冲信号，经过对回波脉冲进行两次FFT处理的运算量即可以获得合成距离高分辨率。这种信号综合了步进频和线性调频脉冲信号的优点，已有较广应用，但是现有文献对调频步进雷达的分析不多。
脉冲积累成像是逆合成孔径(ISAR)雷达的常用方法，其相较于瞬时成像方法的好处在于能适用于低信噪比，而不利因素在于复杂的运动补偿机制。脉冲积累成像又分为相参和非相参两种方法，非相参方法计算量低，但相参方法能适用于低信噪比。
目标匀速运动时，可以通过包络相关、Keystone变换以及Hough变换等非相参方法，或Radon-傅里叶变换(RFT)等相参方法完成距离徙动校正。
目标机动运动时，回波信号具有高阶相位形式，除了距离徙动外还存在多普勒徙动。可以采用分数阶傅里叶变换(FrFT)、立方相位函数法(Cubic Phase Function,CPF)等方法进行多普勒徙动校正。
分数阶模糊函数(fractional ambiguity function,FRAF)的方法，综合了模糊函数和FRFT的优点，能灵活匹配高次相位信号，将该方法应用于机动目标成像将具有同时补偿距离和多普勒徙动的效果。
发明内容
本发明要解决的技术问题包括：1)在调频步进雷达体制下，机动目标的回波信号会出现距离徙动效应和多普勒徙动效应，导致目标能量跨越多个距离单元和多普勒单元，距离-多普勒成像会出现畸变，另外还会影响目标能量积累，在低信噪比下会出现没有足够信号能量，无法成像的情况；2)传统的相参累积技术方案在不能很好兼顾低信噪比和目标机动运动情况下的能量相参积累问题。
为了解决上述技术问题，本发明提供的一种机动目标调频步进逆合成孔径成像方法，包括如下步骤：
步骤1，将调频步进雷达回波信号经过本地的参考信号进行去斜处理，从而得到粗距离像信号；
步骤2，将选定的成像区域在距离上划分为各个较粗的距离单元，并将预警雷达提供的目标的速度和加速度作为FRAF参数离散区间取值的参考，将每个距离单元的时延作为信号瞬时自相关函数的时延参数离散区间参考；
步骤3，将粗距离像信号进行FRAF变换获得各个距离单元相应的FRAF谱图，比较各个FRAF谱图的峰值大小，将峰值最大的那幅FRAF谱图作为要求的参数估计谱；
步骤4，通过参数估计谱解出目标运动的初始速度和加速度参数，即得到一个脉冲串内目标的初始速度和加速度的估计值；
步骤5，重复执行步骤2至步骤4，从而获得不同脉冲串内目标的初始速度和加速度的估计值，再利用不同脉冲串内目标的初始速度和加速度的估计值拟合出目标精确的初始速度和加速度的估计值；
步骤6，根据目标精确的初始速度和加速度的估计值构造相位补偿函数，并将相位补偿函数和粗距离像信号相乘，从而补偿去除粗距离像信号中的高次相位项；
步骤7，对补偿后的粗距离像信号在慢时间域进行FFT，从而得到目标精距离像信号再对目标精距离像信号进行相位聚焦；
步骤8，对相位聚焦后的目标精距离像信号进行转台成像，得到目标逆合成孔径图像。
作为本发明的进一步限定方案，在步骤6中，构造的相位补偿函数为：
式中，和为目标的初始速度和加速度的精确估计值，m∈[0,M-1]，M为脉冲串内的脉冲数，l∈[0,L-1]，L为脉冲串的个数，△f为频率步进量，Tr为脉冲串重复周期，f0为载波频率，μ为子脉冲调频率，c为光速。
本发明的有益效果在于：
(1)采用FRAF变换能够在低信噪比下能补偿目标机动运动产生的距离和多普勒徙动，可以实现目标能量的相参积累；(2)利用对目标运动参数进行拟合，有效消除了误差，且根据调频步进雷达目标回波特征构造补偿函数，该方法能在低信噪比下得到目标运动精确参数，提升了调频步进雷达在恶劣背景下对机动目标的成像能力。
附图说明
图1为本发明的方法流程图；
图2为本发明的方法验证时的仿真飞机模型；
图3为本发明的方法验证时未经处理的精距离像；
图4为本发明的方法验证时运动补偿后的精距离像；
图5为本发明的方法验证时的FRAF谱图；
图6为本发明的方法验证时的ISAR像；
图7为本发明的方法在不同信噪比条件下对飞机目标的运动参数估计结果。
具体实施方式
为了便于解释机动目标的调频步进雷达信号模型，首先做如下推导：
调频步进雷达发射的每一簇线性调频(LFM)脉冲串可以表示为：
s(t)＝rect[(t-tm)/Tp]exp[jπμ(t-tm)2]exp[j2π(f0+mΔf)(t-tm)]
式中，t为全时间，tm＝mTr+lMTr为慢时间，m∈[0,M-1]，M为一簇子脉冲串的脉冲数，l∈[0,L-1]，L为子脉冲串的个数，△f为频率步进量，Tp为子脉冲宽度，Tr为子脉冲重复周期，f0为载波频率，μ为子脉冲调频率。
调频步进雷达接收到的目标回波为：
sr(t)＝rect[(t-tm-2R/c)/Tp]exp[jπμ(t-tm-2R/c)2]exp[j2π(f0+mΔf)(t-tm-2R/c)]
将该回波信号和本地参考信号进行去斜处理，可以得到目标回波的差频信号为：
式中，Rref为去斜参考距离。
将差频信号对快时间进行傅里叶变换就可以得到脉冲串的粗距离像信号：
由于目标在一个Tr内的运动可以忽略，所以在慢时间tm维考虑，具有加速度的机动目标与雷达之间的距离R为：
式中，R0为目标初始距离，v为目标速度，a为目标加速度。
将以上的距离R表达式和tm＝mTr+lMTr代入粗距离像信号公式中，由于距离合成只和mTr相关的相位项有关，所以只保留和mTr相关的低次相位项，则粗距离像信号的相位项部分剔除小值后可以近似为下式：
式中，R△＝R0-Rref，第一个指数项反映了目标的精距离像，第二个和第三个指数项在一个子脉冲串持续时间内为常数，但是在脉冲串间产生精距离像的走动，需要消除，后两个指数项是由于目标运动产生的调频项，同样需消除，可以看出上式是一个三次调频信号，如果能估计出二次和三次调频率，就可以解出目标的速度和加速度。
如图1所示，本发明的机动目标调频步进逆合成孔径成像方法，包括如下步骤：
步骤1，将调频步进雷达回波信号经过本地的参考信号进行去斜处理，从而得到粗距离像信号，具体地，在相参雷达接收端将接收到的回波信号经过发大、限幅处理后，进行距离向和方位向采样得到回波数据，再将雷达发射信号作为去斜参考信号，将参考信号与回波数据差频处理，得到差频信号：
式中，sr(t)是回波数据，是参考信号的复共轭，再将该差频信号存储为距离-方位二维数据矩阵，对该矩阵在距离向进行FFT，得到各个脉冲串的粗距离像信号：
将和tm＝mTr+lMTr代入粗距离像信号公式中，经过化简和忽略小项可以得到粗距离像信号的相位部分信号为：
步骤2，将选定的成像区域在距离上划分为各个较粗的距离单元，并将预警雷达提供的目标的速度和加速度作为FRAF参数离散区间取值的参考，将每个距离单元的时延作为信号瞬时自相关函数的时延参数离散区间参考，具体地，选定的预警雷达提供的成像区域，在距离上将其划分为若干较粗的距离单元，再由预警雷达提供目标的速度、加速度参考和离散参数(α,u)，再将每个距离单元的时延τ作为参数代入FRAF变换：
式中，Rf(tm,τ)为信号的瞬时自相关函数：
Kα(t,u)为分数阶变换核：
式中，
步骤3，将粗距离像信号进行FRAF变换获得各个距离单元相应的FRAF谱图，FRAF谱图与离散参数(α,u)相关，比较各个FRAF谱图的峰值大小，将峰值最大的那幅FRAF谱图作为要求的参数估计谱；
步骤4，通过参数估计谱解出目标运动的初始速度和加速度参数，即得到一个脉冲串内目标的初始速度和加速度的估计值，具体的，通过该参数估计谱得出离散参数α和u的估计值，根据公式12πa3τ+cotα＝0和u＝4πa2τsinα得出参数a2和a3，以a2＝2Δf(v/Tr+alM)/c和a3＝2μa(alMTr+v)/c2解出参数v和a，即得到了一个脉冲串内目标的初始速度和加速度的估计值和
步骤5，重复执行步骤2至步骤4，从而获得不同脉冲串内目标的初始速度和加速度的估计值，再利用不同脉冲串内目标的初始速度和加速度的估计值拟合出目标精确的初始速度和加速度的估计值，在典型的秒级成像时间内，目标的加速度是几乎恒定的，所以可以将不同脉冲串内的加速度估计值求平均得到更加精确的估计值，同时，由于加速度不变，不同脉冲串的初始速度估计值在成像时间内应该连成一直线，根据这一特征，可以将初始速度估计值进行拟合，可以得到更加精确的估计值；
步骤6，根据目标精确的初始速度和加速度的估计值构造相位补偿函数，并将相位补偿函数和粗距离像信号相乘，从而补偿去除粗距离像信号中的高次相位项，构造的相位补偿函数为：
式中，和为目标的初始速度和加速度的精确估计值，m∈[0,M-1]，M为脉冲串内的脉冲数，l∈[0,L-1]，L为脉冲串的个数，△f为频率步进量，Tr为脉冲串重复周期，f0为载波频率，μ为子脉冲调频率，c为光速；
步骤7，对补偿后的粗距离像信号在慢时间域进行FFT，从而得到目标精距离像信号再对目标精距离像信号进行相位聚焦；
步骤8，对相位聚焦后的目标精距离像信号进行转台成像，得到目标逆合成孔径图像。
如图2-6所示，为本发明算法的仿真结果，可以看出在没有噪声的仿真参数条件下可以精确成像。在仿真中，假设雷达发射信号的参数为：子脉冲宽度Tp＝1μs，载波频率f0＝10GHz，子脉冲重复周期Tr＝10μs，子脉冲串包含脉冲数M＝200，子脉冲带宽B＝4MHz，频率步进量△f＝4MHz，合成带宽800MHz，子脉冲串个数L＝500，即成像观测时间1s，目标初始距离为R0＝50km，初始速度v＝1000m/s，加速度a＝30m/s2，目标转动速度4°/s。
如图7所示，是在图2-6的参数条件下，在粗距离像信号中加入不同程度的噪声后，本发明方法的参数估计结果，从中可以看出本文方法的参数估计性能，在信噪比低至-20dB时是稳定有效的，信噪比在-25dB时运动参数估计出现较大误差。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉该技术的人在在本发明所揭露的技术范围内，可延伸到其它的修改、变化和应用，都应涵盖在本发明的包含范围之内。
Claims (2)
1.一种机动目标调频步进逆合成孔径成像方法，其特征在于，包括如下步骤：
步骤1，将调频步进雷达回波信号经过本地的参考信号进行去斜处理，从而得到粗距离像信号；
步骤2，将选定的成像区域在距离上划分为各个较粗的距离单元，并将预警雷达提供的目标的速度和加速度作为FRAF参数离散区间取值的参考，将每个距离单元的时延作为信号瞬时自相关函数的时延参数离散区间参考；
步骤3，将粗距离像信号进行FRAF变换获得各个距离单元相应的FRAF谱图，比较各个FRAF谱图的峰值大小，将峰值最大的那幅FRAF谱图作为要求的参数估计谱；
步骤4，通过参数估计谱解出目标运动的初始速度和加速度参数，即得到一个脉冲串内目标的初始速度和加速度的估计值；
步骤5，重复执行步骤2至步骤4，从而获得不同脉冲串内目标的初始速度和加速度的估计值，再利用不同脉冲串内目标的初始速度和加速度的估计值拟合出目标精确的初始速度和加速度的估计值；
步骤6，根据目标精确的初始速度和加速度的估计值构造相位补偿函数，并将相位补偿函数和粗距离像信号相乘，从而补偿去除粗距离像信号中的高次相位项；
步骤7，对补偿后的粗距离像信号在慢时间域进行FFT，从而得到目标精距离像信号，再对目标精距离像信号进行相位聚焦；
步骤8，对相位聚焦后的目标精距离像信号进行转台成像，得到目标逆合成孔径图像。
2.根据权利要求1所述的机动目标调频步进逆合成孔径成像方法，其特征在于，在步骤6中，构造的相位补偿函数为：
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Phi;</mi>
<mo>=</mo>
<mi>exp</mi>
<mrow>
<mo>{</mo>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>mT</mi>
<mi>r</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mn>2</mn>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mover>
<mi>v</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<mover>
<mi>a</mi>
<mo>~</mo>
</mover>
<msub>
<mi>lMT</mi>
<mi>r</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mi>c</mi>
</mrow>
<mo>}</mo>
</mrow>
<mi>exp</mi>
<mrow>
<mo>{</mo>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>mT</mi>
<mi>r</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>&Delta;</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mrow>
<mn>2</mn>
<mover>
<mi>v</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<mover>
<mi>a</mi>
<mo>~</mo>
</mover>
<msub>
<mi>lMT</mi>
<mi>r</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>lMT</mi>
<mi>r</mi>
</msub>
<mo>/</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>cT</mi>
<mi>r</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>}</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>&times;</mo>
<mi>exp</mi>
<mrow>
<mo>{</mo>
<mrow>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>mT</mi>
<mi>r</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&CenterDot;</mo>
<mn>2</mn>
<mi>&Delta;</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mrow>
<mover>
<mi>v</mi>
<mo>~</mo>
</mover>
<mo>/</mo>
<msub>
<mi>T</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<mover>
<mi>a</mi>
<mo>~</mo>
</mover>
<mi>l</mi>
<mi>M</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mi>c</mi>
</mrow>
<mo>}</mo>
</mrow>
<mi>exp</mi>
<mrow>
<mo>{</mo>
<mrow>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>mT</mi>
<mi>r</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
<mo>&CenterDot;</mo>
<mn>2</mn>
<mi>&mu;</mi>
<mover>
<mi>a</mi>
<mo>~</mo>
</mover>
<mrow>
<mo>(</mo>
<mrow>
<mover>
<mi>a</mi>
<mo>~</mo>
</mover>
<msub>
<mi>lMT</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<mover>
<mi>v</mi>
<mo>~</mo>
</mover>
</mrow>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
</mrow>
<mo>}</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中，和为目标的初始速度和加速度的精确估计值，m∈[0,M-1]，M为脉冲串内的脉冲数，l∈[0,L-1]，L为脉冲串的个数，△f为频率步进量，Tr为脉冲串重复周期，f0为载波频率，μ为子脉冲调频率，c为光速。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510738038.3A CN105301590B (zh) | 2015-11-03 | 2015-11-03 | 一种机动目标调频步进逆合成孔径成像方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510738038.3A CN105301590B (zh) | 2015-11-03 | 2015-11-03 | 一种机动目标调频步进逆合成孔径成像方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105301590A CN105301590A (zh) | 2016-02-03 |
| CN105301590B true CN105301590B (zh) | 2017-08-29 |
Family
ID=55199069
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510738038.3A Active CN105301590B (zh) | 2015-11-03 | 2015-11-03 | 一种机动目标调频步进逆合成孔径成像方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105301590B (zh) |
Families Citing this family (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106054154B (zh) * | 2016-05-27 | 2018-04-06 | 南京邮电大学 | 一种基于icpf的机动目标的步进频率信号高分辨率成像方法 |
| CN108072864B (zh) * | 2016-11-07 | 2021-03-26 | 中国科学院沈阳自动化研究所 | 一种基于变载频调频序列的多目标探测方法 |
| CN106814359B (zh) * | 2016-12-09 | 2019-06-04 | 北京无线电计量测试研究所 | 一种机场跑道异物检测雷达信号预处理方法 |
| CN106597440B (zh) * | 2016-12-12 | 2018-11-27 | 南京信息职业技术学院 | 一种调频步进雷达低信噪比成像方法 |
| CN107085212B (zh) * | 2017-04-01 | 2023-10-24 | 中国人民解放军空军工程大学 | 一种基于线性调频步进信号的自旋目标时变三维成像方法 |
| CN107085213A (zh) * | 2017-05-19 | 2017-08-22 | 中国人民解放军63892部队 | 基于随机调频步进波形设计的运动目标isar成像方法 |
| CN107861125B (zh) * | 2017-11-01 | 2021-09-03 | 河海大学 | 基于广义帕累托分布的贝叶斯压缩感知成像方法 |
| CN108594228B (zh) * | 2018-04-28 | 2022-02-22 | 西安电子科技大学 | 基于isar图像重聚焦的空间目标姿态估计方法 |
| CN108490417B (zh) * | 2018-05-22 | 2020-09-08 | 电子科技大学 | 一种精确的sar动目标参数估计方法 |
| CN108919221A (zh) * | 2018-07-17 | 2018-11-30 | 武汉大学 | 一种用于变加速运动目标的相参积累检测方法 |
| CN109001708B (zh) * | 2018-08-05 | 2020-09-04 | 中国人民解放军海军航空大学 | 基于分级积累检测的雷达机动目标快速精细化处理方法 |
| CN109709552B (zh) * | 2019-01-24 | 2022-04-05 | 南京信息职业技术学院 | 一种低信噪比isar成像运动补偿方法 |
| CN111999725B (zh) * | 2020-09-01 | 2023-07-11 | 中国电子科技集团公司第三十八研究所 | 窄带信号引导下的宽带多项式相位信号去斜方法及装置 |
| CN113985411A (zh) * | 2021-11-01 | 2022-01-28 | 湖南大学 | 一种基于高超声速平台随机跳频雷达子孔径合成成像方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102998672A (zh) * | 2012-11-27 | 2013-03-27 | 西安电子科技大学 | 基于相干化处理的步进频率isar成像方法 |
| CN105005045A (zh) * | 2015-06-26 | 2015-10-28 | 西安电子科技大学 | 基于信号预处理的高速目标isar步进频信号合成方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8138961B2 (en) * | 2009-03-24 | 2012-03-20 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Step frequency ISAR |
-
2015
- 2015-11-03 CN CN201510738038.3A patent/CN105301590B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102998672A (zh) * | 2012-11-27 | 2013-03-27 | 西安电子科技大学 | 基于相干化处理的步进频率isar成像方法 |
| CN105005045A (zh) * | 2015-06-26 | 2015-10-28 | 西安电子科技大学 | 基于信号预处理的高速目标isar步进频信号合成方法 |
Non-Patent Citations (3)
| Title |
|---|
| 匀加速运动对频率步进雷达目标一维距离像的影响及其运动补偿方法;崔应留 等;《弹箭与制导学报》;20141231;第34卷(第6期);135-138 * |
| 步进频率雷达目标的运动参数估计;刘峥 等;《电子学报》;20000331;第28卷(第3期);43-45 * |
| 步进频雷达目标ISAR成像运动补偿新方法;罗贤全 等;《宇航学报》;20070930;第28卷(第5期);1288-1294 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105301590A (zh) | 2016-02-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105301590B (zh) | 一种机动目标调频步进逆合成孔径成像方法 | |
| CN106970371B (zh) | 一种基于Keystone和匹配滤波的目标检测方法 | |
| CN106405541B (zh) | 全相参连续波多普勒雷达及其测距测速方法 | |
| CN104865569B (zh) | 一种基于单频网外辐射源雷达的飞机目标识别方法 | |
| CN103744068B (zh) | 双通道调频连续波sar系统的动目标检测成像方法 | |
| CN102426360B (zh) | 空中微动旋转目标的二维isar成像方法 | |
| US4450444A (en) | Stepped frequency radar target imaging | |
| US7450057B2 (en) | Signal processing for accelerating moving targets | |
| USH1720H (en) | Time frequency processor for radar imaging of moving targets | |
| CN107132534B (zh) | 一种高速雷达目标频域检测的优化方法 | |
| Song et al. | The velocity and range detection using the 2D-FFT scheme for automotive radars | |
| CN109613507B (zh) | 一种针对高阶机动目标雷达回波的检测方法 | |
| CN108051809A (zh) | 基于Radon变换的运动目标成像方法、装置及电子设备 | |
| Li et al. | Scaled Radon-Wigner transform imaging and scaling of maneuvering target | |
| CN111736128A (zh) | 基于skt-siaf-mscft的相参积累方法 | |
| CN106842166A (zh) | 一种适用于线性调频连续波雷达系统的解速度模糊方法 | |
| CN106597440B (zh) | 一种调频步进雷达低信噪比成像方法 | |
| CN103364783B (zh) | 基于单通道sar的运动目标径向速度无模糊估计方法 | |
| CN102540189A (zh) | 基于复数后向投影的自旋目标三维成像方法 | |
| CN110850384B (zh) | 一种基于扫频数据产生宽带去斜回波的方法 | |
| Wacks et al. | Passive synthetic aperture hitchhiker imaging of ground moving targets—Part 1: Image formation and velocity estimation | |
| CN1327242C (zh) | 逆合成孔径雷达机动多目标相对运动补偿方法 | |
| CN106291489A (zh) | 适用于多种发射信号波形的合成孔径雷达回波仿真方法 | |
| CN111856466A (zh) | 一种高效的复杂运动目标isar平动补偿方法 | |
| Karakasiliotis et al. | Two-dimensional ISAR model and image reconstruction with stepped frequency-modulated signal |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20221216 Address after: 102600 103-212, Floor 1, Building 6, Yard 1, Hongye East Road, Daxing District, Beijing Patentee after: Zhongji Zhonglian (Beijing) Energy Technology Co.,Ltd. Address before: Nanjing City, Jiangsu province 210046 City Xianlin University No. 99 Patentee before: Nanjing Vocational College of Information Technology |