CN107966686A - 一种基于线性调频的时间依赖性频控阵目标探测方法 - Google Patents

一种基于线性调频的时间依赖性频控阵目标探测方法 Download PDF

Info

Publication number
CN107966686A
CN107966686A CN201711119620.7A CN201711119620A CN107966686A CN 107966686 A CN107966686 A CN 107966686A CN 201711119620 A CN201711119620 A CN 201711119620A CN 107966686 A CN107966686 A CN 107966686A
Authority
CN
China
Prior art keywords
mrow
msub
angle
distance
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.)
Granted
Application number
CN201711119620.7A
Other languages
English (en)
Other versions
CN107966686B (zh
Inventor
欧阳缮
徐子安
谢跃雷
廖可非
黄东红
李晶晶
黄高见
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guilin University of Electronic Technology
Original Assignee
Guilin University of Electronic Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guilin University of Electronic Technology filed Critical Guilin University of Electronic Technology
Priority to CN201711119620.7A priority Critical patent/CN107966686B/zh
Publication of CN107966686A publication Critical patent/CN107966686A/zh
Application granted granted Critical
Publication of CN107966686B publication Critical patent/CN107966686B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28Details of pulse systems

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

本发明公开了一种基于线性调频的时间依赖性频控阵目标探测方法,其特征是,包括如下步骤:1)确定LFM‑FDA发射波束峰值时的距离、角度和频率增量与时间的关系;2)确定时间依赖性频率增量表达;3)得到时间依赖性频率增量;4)得到距离‑时间以及角度‑时间的关系;5)进行目标搜索。这种方法的优点是使得频控阵的主旁瓣比得到提高、能量汇聚效果更集中,且能使得频控阵波束能量持续汇聚在目标点所在的位置,使得能量在目标位置得到持续积累,提高了探测准确性。

Description

一种基于线性调频的时间依赖性频控阵目标探测方法
技术领域
本发明涉及阵列信号处理领域,具体涉是一种基于线性调频频控阵(LinearFrequency Modulation Frequency Diverse Array,LFM-FDA)的雷达探测技术,尤其是一种基于线性调频的时间依赖性频控阵目标探测方法。
背景技术
频控阵相较于相控阵而言,其相邻两个阵元之间存在一个远小于基准载频的频率增量。这使得它的波束方向图不仅与角度有关,且依赖于距离和时间,这在雷达目标探测成像中有很大的应用前景。
雷达定位技术广泛应用于航天、军事、气候探测等领域,有着很大的需求。在频控阵中使用线性调频技术,可以通过增加每个阵元的发射带宽来提高雷达定位的效果,使得主旁瓣比更大,进而使得雷达的能量汇聚更加集中。对于一个目标位置,将其在阵列的一次探测时间内实现持续照射,则可以大大提高目标探测准确性。在实际探测中,由于目标的位置并非预先知道的,所以需要在距离或角度进行扫描。
发明内容
本发明的目的是针对现有技术的不足,而提供一种基于线性调频的时间依赖性频控阵目标探测方法。这种方法的优点是使得频控阵的主旁瓣比得到提高、能量汇聚效果更集中,且能使得频控阵波束能量持续汇聚在目标点所在的位置,使得能量在目标位置得到持续积累,提高了探测准确性。
实现本发明目的的技术方案是:
一种基于线性调频的时间依赖性频控阵目标探测方法,包括如下步骤:
1)确定LFM-FDA发射波束峰值时的距离、角度和频率增量与时间的关系:通过LFM-FDA方向图峰值的表达式得出LFM-FDA发射波束峰值时的距离、角度和频率增量与时间的关系,假定线性调频频控阵的各阵元的发射频率依次线性增加,第m个阵元发射信号的载频fm为:
fm=f0+(m-1)Δf m=1,…,M,
第m个阵元到达远场参考点的距离为:
Rm=R0-(m-1)d sinθ m=1,...,M,
第m个阵元发射信号的表达式为:
由于线性调频频控阵每一次工作周期都是一个非常小的数值,因此其关于时间的二次项部分在计算其周期性时可以忽略,则LFM-FDA方向图峰值的表达式为:
其中f0为频控阵的基准载频,M为总的阵元个数,Δf为阵元间频率增量,R0为第一个阵元距离远场参考点的距离,d为相邻阵元的间距,θ为阵元到达远场的波束与法线方向的夹角,c为光速,B为每个阵元的带宽,T为阵元每次发射的脉冲持续时间;
2)确定时间依赖性频率增量表达式:结合步骤1)中方向图峰值的表达式,并将其写为与时间相关的表达形式:
n设置为1;
3)得到时间依赖性频率增量Δfp(t):对于目标点(θp,Rp),其中θp为目标的角度,Rp为目标的距离,结合步骤2)中时间依赖性LFM-FDA频率增量的表达式,在阵列的一次探测时间内得到目标点(θp,Rp)角度-距离下的时间依赖性频率增量Δfp(t)为:
n设置为1;
4)得到距离-时间以及角度-时间的仿真图:结合步骤3)得到的时间相关频率增量Δfp(t)与LFM-FDA频控阵信号,将整个仿真区域按距离和角度划分网格点,(θq,Rq)为其中一个距离-角度网格点,对LFM-FDA每个阵元的发射波束,针对该网格点进行时延补偿,再将各阵元经时延补偿后的发射波束进行累加,求得累加信号的模值并进行平方,在整个距离-角度网格中得到一个时间、距离和角度的关系函数:
将目标的角度θp带入关系函数,得到距离-时间关系函数:
再将目标的距离Rp带入关系函数,得到角度-时间关系函数:
依次遍历仿真范围内每个网格点,并在距离-角度维上进行仿真;
5)进行目标搜索:将阵列的一次探测工作时间划分为5-10个时间段,在每个时间段内对仿真范围内的距离和角度上进行步进,并遍历仿真中所需要的探测范围,完成在距离域和角度域的目标搜索过程,首先设置探测区域中一个感兴趣的距离R,使用LFM-FDA时间依赖性频控阵进行角度扫描,在第一个时间段内的初始角度值为θ1,相邻时间段的角度扫描步进值为Δθ,最后一个时间段的终止角度值为θ2,当一次探测时间结束时,相当于对此距离R上的角度范围θ12完成了扫描,再设置探测区域中一个感兴趣的角度θ,使用LFM-FDA时间依赖性频控阵进行距离扫描,在第一个时间段内的初始距离值为R1,相邻时间段的距离扫描步进值为ΔR,最后一个时间段的终止距离值为R2,当一次探测时间结束时,相当于对此角度θ上的距离范围R1-R2完成了扫描。
本发明要解决的技术问题是FDA阵列在融入线性调频技术后,通过设置与时间相关的频率增量,使得在一个时间段内,频控阵波束能量持续汇聚在某个目标区域的问题,分为两个部分内容,首先得出在一个时间段中LFM-FDA能量持续照射一个目标点所在的距离-角度的方法,之后提供一种在实际应用背景下,此种波束在一个感兴趣的角度上对需探测的距离区间进行扫描的过程,及在一个感兴趣的距离上对需探测的角度区间进行扫描的过程。由于频控阵以及线性调频频控阵的发射方向图是S形曲线,在目标的角度上无法持续聚焦到目标的距离,在目标的距离也无法持续聚焦到目标角度值,而时间依赖性频率增量可以实现对目标距离-角度位置的持续能量汇聚,在搜索探测中也能在某个感兴趣的距离上实现角度的步进搜索,在某个感兴趣的距离上对角度上实现步进搜索。
在搜索探测过程中,对于探测区域中感兴趣的角度值,可以使得波束在距离上实现步进搜索;对于探测区域中感兴趣的距离值,可以使得波束在角度上实现步进搜索。
这种方法的优点是使得频控阵的主旁瓣比得到提高、能量汇聚效果更集中,且能使得频控阵波束能量持续汇聚在目标点所在的位置,使得能量在目标位置得到持续积累,提高了探测准确性。
附图说明
图1为实施例中LFM-FDA阵列的结构示意图;
图2为实施例中LFM-FDA使用时间依赖性频偏后,距离和时间的关系示意图;
图3为实施例中LFM-FDA使用时间依赖性频偏后,角度和时间的关系示意图;
图4为实施例中LFM-FDA对于固定方位进行的距离搜索的图像;
图5为实施例中LFM-FDA对于固定距离进行的角度搜索的图像;
图6为实施例的方法流程示意图。
具体实施方式
下面结合附图和实施例对本发明内容做进一步的阐述,但不是对本发明的限定。
实施例:
参照图6,一种基于线性调频的时间依赖性频控阵目标探测方法,包括如下步骤:
1)确定LFM-FDA发射波束峰值时的距离、角度和频率增量与时间的关系:通过LFM-FDA方向图峰值的表达式得出LFM-FDA发射波束峰值时的距离、角度和频率增量与时间的关系,假定线性调频频控阵的各阵元的发射频率依次线性增加,第m个阵元发射信号的载频fm为:
fm=f0+(m-1)Δf m=1,…,M,
第m个阵元到达远场参考点的距离为:
Rm=R0-(m-1)d sinθ m=1,...,M,
第m个阵元发射信号的表达式为:
由于线性调频频控阵每一次工作周期都是一个非常小的数值,因此其关于时间的二次项部分在计算其周期性时可以忽略,则LFM-FDA方向图峰值的表达式为:
其中f0为频控阵的基准载频,M为总的阵元个数,Δf为阵元间频率增量,R0为第一个阵元距离远场参考点的距离,d为相邻阵元的间距,θ为阵元到达远场的波束与法线方向的夹角,c为光速,B为每个阵元的带宽,T为阵元每次发射的脉冲持续时间;
2)确定时间依赖性频率增量表达式:结合步骤1)中方向图峰值的表达式,并将其写为与时间相关的表达形式:
n设置为1;
3)得到时间依赖性频率增量Δfp(t):对于目标点(θp,Rp),其中θp为目标的角度,Rp为目标的距离,结合步骤2)中时间依赖性LFM-FDA频率增量的表达式,在阵列的一次探测时间内得到目标点(θp,Rp)角度-距离下的时间依赖性频率增量Δfp(t)为:
n设置为1;
4)得到距离-时间以及角度-时间的仿真图:结合步骤3)得到的时间相关频率增量Δfp(t)与LFM-FDA频控阵信号,将整个仿真区域按距离和角度划分网格点,(θq,Rq)为其中一个距离-角度网格点,对LFM-FDA每个阵元的发射波束,针对该网格点进行时延补偿,再将各阵元经时延补偿后的发射波束进行累加,求得累加信号的模值并进行平方,在整个距离-角度网格中得到一个时间、距离和角度的关系函数:
将目标的角度θp带入关系函数,得到目标位置的距离-时间关系函数:
再将目标的距离Rp带入关系函数,得到目标位置的角度-时间关系函数:
依次遍历仿真范围内每个网格点,并在距离-角度维上进行仿真;
5)进行目标搜索:将阵列的一次探测工作时间划分为5-10个时间段,在每个时间段内对仿真范围内的距离和角度上进行步进,并遍历仿真中所需要的探测范围,完成在距离域和角度域的目标搜索过程,首先设置探测区域中一个感兴趣的距离R,使用LFM-FDA时间依赖性频控阵进行角度扫描,在第一个时间段内的初始角度值为θ1,相邻时间段的角度扫描步进值为Δθ,最后一个时间段的终止角度值为θ2,当一次探测时间结束时,相当于对此距离R上的角度范围θ12完成了扫描,再设置探测区域中一个感兴趣的角度θ,使用LFM-FDA时间依赖性频控阵进行距离扫描,在第一个时间段内的初始距离值为R1,相邻时间段的距离扫描步进值为ΔR,最后一个时间段的终止距离值为R2,当一次探测时间结束时,相当于对此角度θ上的距离范围R1-R2完成了扫描。
本技术方案的效果可以由以下仿真结果进一步说明:
一、仿真条件:
线性调频频控阵的天线模型采用如图1所示的均匀线阵,其中阵列中的阵元数M为16,雷达信号的基准载频为f0为10GHz,阵元间距d为半波长,线性调频带宽B设置为1MHz,脉冲持续时间T为0.1ms,仿真中时间域t取(0-0.1)ms,成像距离范围:(0-1200)Km,距离域扫描间隔是120m;当阵列处于照射特定方位的模式时,成像角度范围:(-90-90)°,角度域扫描间隔为1.8°,在距离搜索过程中,起始距离为15Km,距离的步进值为5Km,终止距离为50Km;在角度搜索过程中,其实角度为-30°,角度的步进值为15°,终止角度为60°。
二、仿真内容与结果:
仿真1:设置一个点目标,相对于LFM-FDA阵列来说,其方位角为θp=30°,径向距离θp=300Km,采用本发明进行针对单个目标的定位,图2所示为角度为30°时,在阵元间使用时间依赖性频率增量Δfp(t)时距离和时间的关系图,在图中可以看出,当距离为300Km时,在整个0-0.1ms的时间范围内,能量曲线的斜率为0,即在此时间段内能量持续照射在这个距离上,而在其他有能量汇聚的距离区域,曲线的斜率不为0,因此能量无法持续汇聚在这些距离值上,这种特性提高了目标被探测到的概率;图3所示为距离300Km时,在阵元间使用与角度相关的时间依赖性频率增量Δfp(t)时角度和时间的关系图,在图中可以看出,当角度为30°时,在整个0-0.1ms的时间范围内,阵列的发射能量持续汇聚在此角度值上,因此在30°处斜率为0,综上所述,通过使用由目标的角度-距离二维信息所求得的时间依赖性频率增量,可以使阵列实现在一个时间段内对一个目标位置(θp,Rp)的持续照射。
仿真2:使用时间依赖性频率增量Δfp(t)在一定的距离范围和角度范围进行遍历搜索,先在感兴趣的角度上对距离扫描过程进行仿真,设置一个感兴趣的角度为30°,将整个扫描时间0.1ms划分为8段,每段历时0.125ms,起始扫描距离R1=15Km,扫描间隔ΔR=5Km,扫描的终止距离R2=50Km,仿真结果如图4所示,可见在每一个时间段内,能量主瓣只对一个距离值实现持续的照射,从15Km开始,以5Km为增量逐步增至50Km,可以在一个感兴趣的角度对整个距离范围进行依次探测,再在感兴趣的距离上对角度扫描过程进行仿真,设置感兴趣的距离为300Km,将时间区间像上文一样划分为8段,起始扫描角度为-30°,扫描间隔为15°,扫描的终止角度为60°,仿真结果如图5所示,可见在一个时间段内,能量主瓣只对一个角度值实现持续的照射,从-30°开始,以15°为增量逐步增至60°,可以在一个固定的距离上对整个方位范围进行依次探测,但当方位角较大时,其性能也将变差,仿真结果如图5所示。
综上所述,本技术方案在频控阵的基础上应用了线性调频技术,通过使用与时间依赖性频率增量Δfp(t)实现在阵列一次探测时间内对目标位置进行持续照射,对于其他区域不进行持续照射,使得目标区域照射能量得到累积,对目标的探测准确性有促进作用,之后又通过距离和角度步进的方式完成了在感兴趣的角度上对于所取距离范围的扫描及在感兴趣的距离上对所取角度范围的扫描。

Claims (1)

1.一种基于线性调频的时间依赖性频控阵目标探测方法,其特征是,包括如下步骤:
1)确定LFM-FDA发射波束峰值时的距离、角度和频率增量与时间的关系:通过LFM-FDA方向图峰值的表达式得出LFM-FDA发射波束峰值时的距离、角度和频率增量与时间的关系,假定线性调频频控阵的各阵元的发射频率依次线性增加,第m个阵元发射信号的载频fm为:
fm=f0+(m-1)Δf m=1,…,M,
第m个阵元到达远场参考点的距离为:
Rm=R0-(m-1)dsinθ m=1,...,M,
第m个阵元发射信号的表达式为:
<mrow> <msub> <mi>S</mi> <mi>m</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mi>j</mi> <mn>2</mn> <mi>&amp;pi;</mi> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mn>0</mn> </msub> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <mi>m</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> <mi>&amp;Delta;</mi> <mi>f</mi> </mrow> <mo>)</mo> <mi>t</mi> <mo>+</mo> <mi>j</mi> <mi>&amp;pi;</mi> <mfrac> <mi>B</mi> <mi>T</mi> </mfrac> <msup> <mi>t</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mo>,</mo> <mn>0</mn> <mo>&amp;le;</mo> <mi>t</mi> <mo>&amp;le;</mo> <mi>T</mi> </mrow>
则LFM-FDA方向图峰值的表达式为:
<mrow> <mn>2</mn> <mi>&amp;pi;</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>f</mi> <mn>0</mn> </msub> <mi>d</mi> <mi> </mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;theta;</mi> </mrow> <mi>c</mi> </mfrac> <mo>-</mo> <mfrac> <mrow> <msub> <mi>&amp;Delta;fR</mi> <mn>0</mn> </msub> </mrow> <mi>c</mi> </mfrac> <mo>+</mo> <mi>&amp;Delta;</mi> <mi>f</mi> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mn>2</mn> <mi>n</mi> <mi>&amp;pi;</mi> </mrow>
其中f0为频控阵的基准载频,M为总的阵元个数,Δf为阵元间频率增量,R0为第一个阵元距离远场参考点的距离,d为相邻阵元的间距,θ为阵元到达远场的波束与法线方向的夹角,c为光速,B为每个阵元的带宽,T为阵元每次发射的脉冲持续时间;
2)确定时间依赖性频率增量表达式:结合步骤1)中方向图峰值的表达式,并将其写为与时间相关的表达形式:
n设置为1;
3)得到时间依赖性频率增量Δfp(t):对于目标点(θp,Rp),其中θp为目标的角度,Rp为目标的距离,结合步骤2)中时间依赖性LFM-FDA频率增量的表达式,在阵列的一次探测时间内得到目标点(θp,Rp)角度-距离下的时间依赖性频率增量Δfp(t)为:
n设置为1;
4)得到距离-时间以及角度-时间的仿真图:结合步骤3)得到的时间相关频率增量Δfp(t)与LFM-FDA频控阵信号,将整个仿真区域按距离和角度划分网格点,(θq,Rq)为其中一个距离-角度网格点,对LFM-FDA每个阵元的发射波束,针对该网格点进行时延补偿,再将各阵元经时延补偿后的发射波束进行累加,求得累加信号的模值并进行平方,在整个距离-角度网格中得到一个时间、距离和角度的关系函数:
<mrow> <mi>P</mi> <mo>=</mo> <msup> <mrow> <mo>|</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <mi>exp</mi> <mo>(</mo> <mi>j</mi> <mn>2</mn> <mi>&amp;pi;</mi> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mn>0</mn> </msub> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <mi>m</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> <msub> <mi>&amp;Delta;f</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mi>q</mi> </msub> <mo>-</mo> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>d</mi> <mi> </mi> <msub> <mi>sin&amp;theta;</mi> <mi>q</mi> </msub> </mrow> <mi>c</mi> </mfrac> <mo>)</mo> </mrow> <mo>+</mo> <mi>j</mi> <mi>&amp;pi;</mi> <mfrac> <mi>B</mi> <mi>T</mi> </mfrac> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mi>q</mi> </msub> <mo>-</mo> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>d</mi> <mi> </mi> <msub> <mi>sin&amp;theta;</mi> <mi>q</mi> </msub> </mrow> <mi>c</mi> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <mo>)</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>,</mo> <mn>0</mn> <mo>&amp;le;</mo> <mi>t</mi> <mo>&amp;le;</mo> <mi>T</mi> </mrow>
将目标的角度θp带入关系函数,得到目标位置的距离-时间关系函数:
<mrow> <msub> <mi>P</mi> <mi>R</mi> </msub> <mo>=</mo> <msup> <mrow> <mo>|</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <mi>exp</mi> <mo>(</mo> <mi>j</mi> <mn>2</mn> <mi>&amp;pi;</mi> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mn>0</mn> </msub> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <mi>m</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> <msub> <mi>&amp;Delta;f</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mi>q</mi> </msub> <mo>-</mo> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>d</mi> <mi> </mi> <msub> <mi>sin&amp;theta;</mi> <mi>p</mi> </msub> </mrow> <mi>c</mi> </mfrac> <mo>)</mo> </mrow> <mo>+</mo> <mi>j</mi> <mi>&amp;pi;</mi> <mfrac> <mi>B</mi> <mi>T</mi> </mfrac> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mi>q</mi> </msub> <mo>-</mo> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>d</mi> <mi> </mi> <msub> <mi>sin&amp;theta;</mi> <mi>p</mi> </msub> </mrow> <mi>c</mi> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <mo>)</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>,</mo> <mn>0</mn> <mo>&amp;le;</mo> <mi>t</mi> <mo>&amp;le;</mo> <mi>T</mi> <mo>,</mo> </mrow>
再将目标的距离Rp带入关系函数,得到目标位置的角度-时间关系函数:
<mrow> <msub> <mi>P</mi> <mi>&amp;theta;</mi> </msub> <mo>=</mo> <msup> <mrow> <mo>|</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <mi>exp</mi> <mo>(</mo> <mi>j</mi> <mn>2</mn> <mi>&amp;pi;</mi> <mrow> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mn>0</mn> </msub> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <mi>m</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> <msub> <mi>&amp;Delta;f</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mi>p</mi> </msub> <mo>-</mo> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>d</mi> <mi> </mi> <msub> <mi>sin&amp;theta;</mi> <mi>q</mi> </msub> </mrow> <mi>c</mi> </mfrac> <mo>)</mo> </mrow> <mo>+</mo> <mi>j</mi> <mi>&amp;pi;</mi> <mfrac> <mi>B</mi> <mi>T</mi> </mfrac> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>R</mi> <mi>p</mi> </msub> <mo>-</mo> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>d</mi> <mi> </mi> <msub> <mi>sin&amp;theta;</mi> <mi>q</mi> </msub> </mrow> <mi>c</mi> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <mo>)</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>,</mo> <mn>0</mn> <mo>&amp;le;</mo> <mi>t</mi> <mo>&amp;le;</mo> <mi>T</mi> <mo>,</mo> </mrow>
依次遍历仿真范围内每个网格点,并在距离-角度维上进行仿真;
5)进行目标搜索:将阵列的一次探测工作时间划分为5-10个时间段,在每个时间段内对仿真范围内的距离和角度上进行步进,并遍历仿真中所需要的探测范围,完成在距离域和角度域的目标搜索过程,首先设置探测区域中一个感兴趣的距离R,使用LFM-FDA时间依赖性频控阵进行角度扫描,在第一个时间段内的初始角度值为θ1,相邻时间段的角度扫描步进值为Δθ,最后一个时间段的终止角度值为θ2,当一次探测时间结束时,相当于对此距离R上的角度范围θ12完成了扫描,再设置探测区域中一个感兴趣的角度θ,使用LFM-FDA时间依赖性频控阵进行距离扫描,在第一个时间段内的初始距离值为R1,相邻时间段的距离扫描步进值为ΔR,最后一个时间段的终止距离值为R2,当一次探测时间结束时,相当于对此角度θ上的距离范围R1-R2完成了扫描。
CN201711119620.7A 2017-11-14 2017-11-14 一种基于线性调频的时间依赖性频控阵目标探测方法 Active CN107966686B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711119620.7A CN107966686B (zh) 2017-11-14 2017-11-14 一种基于线性调频的时间依赖性频控阵目标探测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711119620.7A CN107966686B (zh) 2017-11-14 2017-11-14 一种基于线性调频的时间依赖性频控阵目标探测方法

Publications (2)

Publication Number Publication Date
CN107966686A true CN107966686A (zh) 2018-04-27
CN107966686B CN107966686B (zh) 2021-05-28

Family

ID=62000061

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711119620.7A Active CN107966686B (zh) 2017-11-14 2017-11-14 一种基于线性调频的时间依赖性频控阵目标探测方法

Country Status (1)

Country Link
CN (1) CN107966686B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108761414A (zh) * 2018-05-25 2018-11-06 电子科技大学 一种基于频控阵的s形干扰波束测试方法
CN109375213A (zh) * 2018-08-08 2019-02-22 西安电子科技大学 基于子阵划分的频率分集阵列信号处理方法
CN110109069A (zh) * 2019-05-14 2019-08-09 桂林电子科技大学 一种去时间依赖性频控阵点状干扰波束形成方法
CN110794400A (zh) * 2019-10-30 2020-02-14 南京航空航天大学 一种基于频控阵的地形匹配方法
CN111624614A (zh) * 2020-06-10 2020-09-04 哈尔滨工业大学 一种通过激光相干探测进行测距的方法及系统
CN113659341A (zh) * 2021-08-12 2021-11-16 阳光学院 一种多相位频控阵波束扫描方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090015474A1 (en) * 2005-12-20 2009-01-15 Wicks Michael C Method and apparatus for a frequency diverse array
KR101302624B1 (ko) * 2012-08-22 2013-09-03 국방과학연구소 펄스 내 주파수 변조 레이더 신호의 변조 형태 인식 방법 및 그 장치
CN103592635A (zh) * 2013-11-11 2014-02-19 电子科技大学 一种认知fda雷达的信号发射方法、装置及雷达
CN106788625A (zh) * 2016-12-21 2017-05-31 电子科技大学 一种基于线性递增频偏策略的频控阵安全通信方法
CN107356921A (zh) * 2017-08-11 2017-11-17 桂林电子科技大学 一种频率分集阵列雷达基于一次频偏目标定位的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090015474A1 (en) * 2005-12-20 2009-01-15 Wicks Michael C Method and apparatus for a frequency diverse array
KR101302624B1 (ko) * 2012-08-22 2013-09-03 국방과학연구소 펄스 내 주파수 변조 레이더 신호의 변조 형태 인식 방법 및 그 장치
CN103592635A (zh) * 2013-11-11 2014-02-19 电子科技大学 一种认知fda雷达的信号发射方法、装置及雷达
CN106788625A (zh) * 2016-12-21 2017-05-31 电子科技大学 一种基于线性递增频偏策略的频控阵安全通信方法
CN107356921A (zh) * 2017-08-11 2017-11-17 桂林电子科技大学 一种频率分集阵列雷达基于一次频偏目标定位的方法

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108761414A (zh) * 2018-05-25 2018-11-06 电子科技大学 一种基于频控阵的s形干扰波束测试方法
CN108761414B (zh) * 2018-05-25 2022-05-20 电子科技大学 一种基于频控阵的s形干扰波束测试方法
CN109375213A (zh) * 2018-08-08 2019-02-22 西安电子科技大学 基于子阵划分的频率分集阵列信号处理方法
CN110109069A (zh) * 2019-05-14 2019-08-09 桂林电子科技大学 一种去时间依赖性频控阵点状干扰波束形成方法
CN110109069B (zh) * 2019-05-14 2022-09-02 桂林电子科技大学 一种去时间依赖性频控阵点状干扰波束形成方法
CN110794400A (zh) * 2019-10-30 2020-02-14 南京航空航天大学 一种基于频控阵的地形匹配方法
CN111624614A (zh) * 2020-06-10 2020-09-04 哈尔滨工业大学 一种通过激光相干探测进行测距的方法及系统
CN111624614B (zh) * 2020-06-10 2022-11-01 哈尔滨工业大学 一种通过激光相干探测进行测距的方法及系统
CN113659341A (zh) * 2021-08-12 2021-11-16 阳光学院 一种多相位频控阵波束扫描方法
CN113659341B (zh) * 2021-08-12 2024-05-28 阳光学院 一种多相位频控阵波束扫描方法

Also Published As

Publication number Publication date
CN107966686B (zh) 2021-05-28

Similar Documents

Publication Publication Date Title
CN107966686A (zh) 一种基于线性调频的时间依赖性频控阵目标探测方法
CN108107431B (zh) 一种圆柱扫描sar三维成像快速实现方法
CN103869311B (zh) 实波束扫描雷达超分辨成像方法
US7532150B1 (en) Restoration of signal to noise and spatial aperture in squint angles range migration algorithm for SAR
Xiong et al. New applications of omega-K algorithm for SAR data processing using effective wavelength at high squint
CN106872974B (zh) 基于高超声速平台双通道雷达的高精度运动目标成像方法
CN106443671A (zh) 基于调频连续波的sar雷达动目标检测与成像方法
CN102749621B (zh) 一种双基地合成孔径雷达频域成像方法
CN102707269B (zh) 一种机载雷达距离走动校正方法
CN101344590B (zh) 一种扇扫sar的广域观测方法
CN108459321B (zh) 基于距离-方位圆模型的大斜视高分辨率sar成像方法
CN102778681B (zh) 固定发射站双基地前视合成孔径雷达成像方法
CN102819020B (zh) 一种俯冲模型方位向非线性变标的合成孔径雷达成像方法
CN104181531B (zh) 一种基于相控阵雷达的三维关联成像方法
CN108226891B (zh) 一种扫描雷达回波计算方法
CN105223572B (zh) 一种基于pfa算法的正前视双基sar成像处理方法
Liang et al. A high-order phase correction approach for focusing HS-SAR small-aperture data of high-speed moving platforms
CN105487074B (zh) 一种双基地合成孔径雷达数值距离多普勒成像方法
CN109814100A (zh) 基于子孔径参数估计的sar地面运动目标成像方法
CN103033811B (zh) 一种基于类单基等效的双基地合成孔径雷达成像方法
CN104330779A (zh) 机载合成孔径雷达运动误差补偿方法
CN102608598A (zh) 基于子空间投影的实孔径前视成像方法
CN107153191B (zh) 一种针对隐形飞机的双基isar成像检测方法
CN105676223B (zh) 一种非匀速太赫兹雷达平台sar成像方法
CN104062659B (zh) 基于压缩感知的滑动聚束sar三维成像方法

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