CN106908771A - 微波暗室中脉冲体制雷达目标探测模拟方法 - Google Patents

微波暗室中脉冲体制雷达目标探测模拟方法 Download PDF

Info

Publication number
CN106908771A
CN106908771A CN201710081095.8A CN201710081095A CN106908771A CN 106908771 A CN106908771 A CN 106908771A CN 201710081095 A CN201710081095 A CN 201710081095A CN 106908771 A CN106908771 A CN 106908771A
Authority
CN
China
Prior art keywords
signal
target
receiving
pulse
transmitting
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
CN201710081095.8A
Other languages
English (en)
Other versions
CN106908771B (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.)
National University of Defense Technology
Original Assignee
National University of Defense 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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN201710081095.8A priority Critical patent/CN106908771B/zh
Publication of CN106908771A publication Critical patent/CN106908771A/zh
Application granted granted Critical
Publication of CN106908771B publication Critical patent/CN106908771B/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
    • G01S7/285Receivers
    • G01S7/292Extracting wanted echo-signals
    • G01S7/2923Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
    • 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
    • G01S13/00Systems 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/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • 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/36Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
    • 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/40Means for monitoring or calibrating
    • G01S7/4052Means for monitoring or calibrating by simulation of echoes

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明一种微波暗室中脉冲体制雷达目标探测模拟方法:第一步,利用间歇收发方式,获取目标回波信号:第二步,目标回波恢复:第三步,脉冲压缩获取目标信息;第四步,目标回波能量补偿信息重构。本发明的有益效果包括:第一,解决了微波暗室中雷达收发分时模式的收发信号耦合问题与收发同时模式中收发信号的遮挡问题。第二,解决了目标回波时域波形的恢复问题。第三,目标信息获取较为精确。第四,目标回波能量补偿与信息重构结果精确。

Description

微波暗室中脉冲体制雷达目标探测模拟方法
【技术领域】
本发明一种微波暗室中脉冲体制雷达目标探测模拟方法,属于雷达电子对抗的仿真领域,具体涉及在微波暗室中模拟脉冲体制雷达信号的收发处理技术,特别是在微波暗室中实现脉冲体制雷达目标探测全过程的模拟与信号恢复方法。
【背景技术】
在微波暗室内进行雷达仿真具有可操作性强、可重复性高、保密性好、成本经济等优点。由于微波暗室空间有限,在静态雷达目标特性测量方面,冲激脉冲或连续波扫频激励的方式成为主要的测量手段,并提供了大量的仿真数据。但是,在实际的雷达电子对抗中,脉冲信号常被作为探测信号用于目标探测。冲激脉冲、扫频波形等难以在暗室中作为等效替代信号,模拟雷达电子对抗的实际过程。因此,在微波暗室中采用与真实雷达探测一致的脉冲信号作为内场辐射式仿真激励源,模拟雷达目标探测的全过程,是实现暗室中电子对抗仿真需要解决的基本问题。
通常,微波暗室环境有限,雷达与目标距离较近,若使用脉冲信号进行目标探测,在收发同时模式中,发射信号与接收信号将会在接收天线处产生互耦;而收发分时情况下,目标回波将会在发射信号未完全辐射前返回,难以实现回波的有效接收。
【发明内容】
本发明要解决的技术问题是:提供一种微波暗室中脉冲体制雷达目标探测模拟方法,在微波暗室中利用间歇收发方法,实现脉冲雷达信号无耦合收发,根据间歇收发控制信号的参数,对目标回波信号进行恢复和能量补偿,实现脉冲信号在暗室中进行目标探测的全过程仿真。
本发明一种微波暗室中脉冲体制雷达目标探测模拟方法,采取的技术方案如下:
第一步,利用间歇收发方式,获取目标回波信号
脉冲信号的收发过程等价于用方波信号p(t)与发射脉冲信号s0(t)相乘,即有s1(t)=s0(t)·p(t),而回波的目标调制过程相当于将发射信号与目标响应进行卷积。设雷达目标时域冲激响应函数为h(t),则间歇收发后目标回波信号可等价表示为s2(t)=h(t)*(s0(t)·p(t)),*表示卷积。
第二步,目标回波恢复
根据发射信号带宽等参数,设计低通滤波器,对目标回波信号s2(t)进行低通滤波,获取间歇收发之后的目标回波。
第三步,脉冲压缩获取目标信息
脉冲压缩,获取间歇收发后目标信息sh(t),然后根据间歇收发参数决定是否采用开窗的方法提取目标信息。
第四步,目标回波能量补偿信息重构
根据间歇收发参数,得到脉冲压缩后信号的能量补偿值,通过能量补偿,重构目标探测信息。
其中,所述第一步间歇收发方式下,收发交替工作过程可以抽象为一方波信号p(t)对雷达脉冲信号进行控制,p(t)可表示为
其中τ为发射通道工作时段,Ts为间歇收发周期,δ(·)为单位冲击函数,n表示第n个冲击脉冲函数。
根据傅氏变换性质,p(t)的频域形式为
其中,fs=1/Ts,sinc(x)=sin(πx)/(πx)。可以发现,间歇收发控制信号的频谱为p(t)各时域谐波的组合,且主峰处幅值由收发参数τfs决定,并决定目标信息重构的能量补偿值。
通过精确的时延和波形控制可以保证实际发射信号的相位连续性和幅度一致性,从而确保收发过程等效为方波信号p(t)与雷达发射脉冲信号s0(t)的乘积。
目标散射过程可等效为一线性系统与激励信号相卷积。设目标与雷达相对距离为R,散射强度为σ0,径向运动速度为v,目标在雷达视线方向投影长度为L。雷达目标时域冲激响应函数为h(t)=σ0δ(t-Δt)。假设雷达发射脉冲信号为一线性调频脉冲信号(上调频)。从而,间歇收发条件下,目标回波信号去载频后可表示为:
其中,A为回波幅度,可由回波功率得到,Pt为发射功率,G为天线收发增益,λ为波长,Δt=2(R-vt)/C,C为电磁波传播速度。rect(·)为矩形窗函数,Tp为脉宽,u(t)=exp(jπγt2)为复包络信号,为单位虚数,f0为中心频率,γ为线性调频率。可以发现,目标回波信号由一系列子脉冲信号构成,各段子脉冲宽度与发射通道工作时段τ相等,子脉冲带宽BΔ=γτ。
其中,所述第二步目标回波恢复,具体方法如下:s2(t)对应频域形式为
其中sinc(·)为辛克函数,U(f)为u(t)的频谱。
根据线性调频LFM信号带宽,设计低通滤波器截止频率刚好覆盖信号带宽,从而间歇收发之后,位于原始LFM信号带宽之外的频谱分量将被有效滤除,进一步可以得到信号时域波形。考虑W(f)为低通滤波器的频谱,理想情况下,滤波器的频谱满足如下条件
从而,经过低通滤波之后,可得
S3(f)=W(f)S2(f) (6)
进一步,可以获得经过低通滤波之后信号的时域形式,如图5所示。
其中,所述第三步脉冲压缩获取目标信息,具体方法如下:在接收机中,根据发射信号可得幅度归一化匹配滤波器频率响应为
其中f是频率分量,B是雷达信号带宽。
根据匹配滤波器特点,W(f)在LFM信号带宽内为1,带宽外为0,与H(f)基本相同,因此相乘后可以忽略W(f)。从而,目标回波信号经匹配滤波并取模后输出形式可以表示为
其中|·|表示取模运算,IFFT表示逆傅里叶变换。对|sh(t)|的主峰位置(n=0)进行时域开窗处理可滤出主峰,完成目标信息重构。需要条件分析如下:
|sh(t)|相邻两个辛克函数尖峰相距
要保证匹配滤波后相邻峰值不发生混叠,则需满足
ΔR>L (10)
此外要满足目标回波信号不遮挡和完全接收条件,还须有
其中表示信号遍历目标耗时。
从而,约束条件可写为
一方面Ts越小(fs越大),则相邻两个峰值距离越远,越利于目标信息重构,一旦Ts足够小,在距离门内将只能观察到实际目标处的尖峰;另一方面Ts越大则信号辐射能量越小,越不易恢复目标信息。
根据上述分析,当间歇收发频率小于信号带宽时,需要在匹配滤波之后,通过合理的开窗,提取脉冲压缩后的主峰,重构目标信息。开窗的目的主要是为了消除n≠0的峰值。在理想情况下,开窗类似于用矩形通断截取脉压主峰,通断时间对应的距离宽要小于相邻峰的距离,从而
因此,截断之后间歇收发回波的匹配滤波输出为
其中,回波延时Δt可以根据加窗前式(4)中的主峰位置估计得到。这样,就可以利用“间歇”收发方法解决空间有限条件下的长脉冲信号目标探测问题。此外,根据式(9),当间歇收发周期Ts与完整信号带宽相当时,相邻虚假峰将无法在距离门内观测到。此时,不采用时域截断即可有效获取目标信息。
其中,所述第四步目标回波能量补偿与信息重构,方法如下:
根据间歇控制信号p(t)的频域形式以及所得间歇收发回波匹配滤波输出,可以发现,对于实际目标信息而言,间歇收发匹配滤波输出的主峰(n=0)反映了真实目标信息。但是由于发射信号的通断,损失了部分信号能量,因此所得目标信息需要进行能量补偿。根据式(8),取n=0处的脉压输出信息,并对信号幅度补偿1/τfs,即可精确的恢复得到目标回波的匹配滤波输出,从而实现暗室中脉冲雷达信号的目标探测模拟过程。恢复后的脉压输出表达式为
本发明的有益效果主要包括:
第一,解决了微波暗室中雷达收发分时模式的收发信号耦合问题与收发同时模式中收发信号的遮挡问题。基于间歇收发方法实现收发通道交替工作,等效于用收发分时方式从时间域对收发通道的耦合信号解耦;根据雷达与目标相对距离、目标尺寸等条件,精确控制收发通道的交替工作时间,将原始雷达发射脉冲信号分成多个短时子脉冲信号进行收发,从而解决了信号遮挡的问题。
第二,解决了目标回波时域波形的恢复问题。通过调制间歇收发参数,使得间歇频率与信号带宽相当,通过低通滤波能够有效获取目标实际回波信号,从而解决了间歇收发后时域波形与理想时域波形的差异问题。
第三,目标信息获取较为精确。理想间歇收发等价于用矩形信号对雷达发射脉冲信号进行“截断”,属于线性作用过程,通过后续信号处理手段可精确去除“截断”效应,保证了目标信息重构结果可靠。根据间歇收发参数,可以决定是否采用开窗的方法获取目标信息,所得目标信息真实可靠。
第四,目标回波能量补偿与信息重构结果精确。根据间歇收发参数,设置回波能量补偿值,能够精确得到补偿后的目标回波信息,完成微波暗室中目标探测的等效模拟。
【附图说明】
图1是微波暗室与间歇收发工作示意图。
图2是本发明方法实现流程图。
图3(a)是原始目标回波时域波形。
图3(b)是原始目标回波频域波形。
图4(a)是τ=0.1us,Ts=0.2us的间歇收发时域波形。
图4(b)是τ=0.1us,Ts=0.2us的间歇收发频域波形。
图5(a)是低通滤波后目标回波时域波形。
图5(b)是低通滤波后目标回波频域波形。
图6是重构目标信息与原始目标信息对比示意图。
图7是重构与原始目标信息的信噪比(SNR)对比示意图。
【具体实施方式】
下面结合附图对本发明作进一步的说明。本发明一种微波暗室中脉冲体制雷达目标探测模拟方法,如图2所示,步骤如下:
第一步,利用间歇收发方法,获取目标回波信号。
间歇收发方式下,收发交替工作过程可以抽象为一方波信号p(t)对雷达脉冲信号进行控制,p(t)可表示为
其中τ为发射通道工作时段,Ts为间歇收发周期,δ(·)为单位冲击函数,n表示第n个冲击脉冲函数。
利用傅氏变换性质,p(t)的频域形式为
其中,fs=1/Ts,sinc(x)=sin(πx)/(πx)。可以发现,间歇收发控制信号的频谱为p(t)各时域谐波的组合,且主峰处幅值由收发参数τfs决定,并决定目标信息重构的能量补偿值。
通过精确的时延和波形控制可以保证实际发射信号的相位连续性和幅度一致性,从而确保收发过程等效为方波信号p(t)与雷达发射脉冲信号s0(t)的乘积,如图1所示。
目标散射过程可等效为一线性系统与激励信号相卷积。设目标与雷达相对距离为R,散射强度为σ0,径向运动速度为v,目标在雷达视线方向投影长度为L。雷达目标时域冲激响应为h(t)=σ0δ(t-Δt)。假设雷达发射脉冲信号为一线性调频脉冲信号(上调频),从而,间歇收发条件下,目标回波信号去载频后可表示为:
其中,A为回波幅度,可由回波功率得到,Pt为发射功率,G为天线收发增益,λ为波长,Δt=2(R-vt)/C,C为电磁波传播速度。rect(·)为矩形窗函数,Tp为脉宽,u(t)=exp(jπγt2)为复包络信号,为单位虚数,f0为中心频率,γ为线性调频率。可以发现,目标回波信号由一系列子脉冲信号构成,各段子脉冲宽度与发射通道工作时段τ相等,子脉冲带宽BΔ=γτ。
第二步,低通滤波恢复目标回波。
s2(t)对应频域形式为
其中sinc(·)为辛克函数,U(f)为u(t)的频谱。
根据LFM信号带宽,设计低通滤波器截止频率刚好覆盖信号带宽,从而间歇收发之后,位于原始LFM信号带宽之外的频谱分量将被有效滤除,进一步可以得到信号时域波形。考虑W(f)为低通滤波器的频谱,理想情况下,滤波器的频谱满足如下条件
从而,经过低通滤波之后,可得
S3(f)=W(f)S2(f) (21)
进一步,可以获得经过低通滤波之后信号的时域形式,如图5所示。
第三步,脉冲压缩获取目标信息。
接收到的目标回波信号理论上等价于将全脉冲发射条件下的目标回波信号进行“截断”,匹配滤波后可以通过开窗手段去除“截断”效应,实现目标信息重构。在接收机中,根据发射信号可得幅度归一化匹配滤波器频率响应为
其中f是频率分量,B是雷达信号带宽。
根据匹配滤波器特点,W(f)在LFM信号带宽内为1,带宽外为0,与H(f)基本相同,因此相乘后可以忽略W(f)。从而,目标回波信号经匹配滤波并取模后输出形式可以表示为
其中|·|表示取模运算,IFFT表示逆傅里叶变换。对|sh(t)|的主峰位置(n=0)进行时域开窗处理可滤出主峰,完成目标信息重构。需要条件分析如下:
|sh(t)|相邻两个辛克函数尖峰相距
要保证匹配滤波后相邻峰值不发生混叠,则需满足
ΔR>L (25)
此外要满足目标回波信号不遮挡和完全接收条件,还须有
其中表示信号遍历目标耗时。
从而,约束条件可写为
一方面Ts越小(fs越大),则相邻两个峰值距离越远,越利于目标信息重构,一旦Ts足够小,在距离门内将只能观察到实际目标处的尖峰;另一方面Ts越大则信号辐射能量越小,越不易恢复目标信息。
根据上述分析,当间歇收发频率小于信号带宽时,需要在匹配滤波之后,通过合理的开窗,提取脉冲压缩后的主峰,重构目标信息。开窗的目的主要是为了消除n≠0的峰值。在理想情况下,开窗类似于用矩形通断截取脉压主峰,通断时间对应的距离宽要小于相邻峰的距离,从而
因此,截断之后间歇收发回波的匹配滤波输出为
其中,回波延时Δt可以根据加窗前式(19)中的主峰位置估计得到。这样,就可以利用“间歇”收发方法解决空间有限条件下的长脉冲信号目标探测问题。此外,根据式(24),当间歇收发周期Ts与完整信号带宽相当时,相邻虚假峰将无法在距离门内观测到,此时,不采用时域截断即可有效获取目标信息。
第四步,目标回波能量补偿与信息重构。
根据间歇控制信号p(t)的频域形式以及所得间歇收发回波匹配滤波输出,可以发现,对于实际目标信息而言,间歇收发匹配滤波输出的主峰(n=0)反映了真实目标信息。但是由于发射信号的通断,损失了部分信号能量,因此所得目标信息需要进行能量补偿。根据式(23),取n=0处的脉压输出信息,并对信号幅度补偿1/τfs,即可精确的恢复得到目标回波的匹配滤波输出,从而实现暗室中脉冲雷达信号的目标探测模拟过程。恢复后的脉压输出表达式为
假设目标是点目标,目标与雷达距离为R=45m,静区反射率电平为-40dB。按照图1所示情景,将雷达收发控制系统与目标置于暗室,使两者的距离为R。雷达长脉冲信号脉宽为Tp=100us,带宽为B=5MHz,波长为0.3m,发射功率为1W,天线收发增益为30dB,接收机信噪比为10dB。目标散射截面积为σ0=0.1m2,速度为v=1m/s,匹配滤波后目标等效位于45m处。根据图2所示处理流程,设定间歇收发参数τ=0.1us,Ts=0.2us。首先,图3(a)、(b)给出完整LFM脉冲信号的时频域波形。然后,得到的间歇收发回波时频域波形如图4(a)、(b)所示。经过低通滤波,得到滤波后的间歇收发时频域波形如图5(a)、(b)所示,其时频域特性与完整信号的时频域波形基本一致,只有幅度上存在差异。由于间歇收发方法中fs=5MHz,与带宽相当,因此,匹配滤波后将观察不到相邻虚假尖峰。进一步,对所得脉压输出进行能量补偿即可获取目标信息,如图6所示。由于存在多普勒,脉压输出峰值位置偏离真实位置45m,但间歇收发方法与完整信号的脉压输出结果基本一致。在不同信噪比条件下,对比间歇收发方法与完整回波信噪比的估计值,得到图7。可以发现,两者信噪比基本一致,从而验证了间歇收发方法的有效性。

Claims (5)

1.一种微波暗室中脉冲体制雷达目标探测模拟方法,其特征在于:该方法步骤如下:
第一步,利用间歇收发方式,获取目标回波信号
脉冲信号的收发过程等价于用方波信号p(t)与发射脉冲信号s0(t)相乘,即有s1(t)=s0(t)·p(t),而回波的目标调制过程相当于将发射信号与目标响应进行卷积;设雷达目标时域冲激响应函数为h(t),则间歇收发后目标回波信号可等价表示为s2(t)=h(t)*(s0(t)·p(t)),*表示卷积;
第二步,目标回波恢复
根据发射信号带宽等参数,设计低通滤波器,对目标回波信号s2(t)进行低通滤波,获取间歇收发之后的目标回波;
第三步,脉冲压缩获取目标信息
脉冲压缩,获取间歇收发后目标信息sh(t),然后根据间歇收发参数决定是否采用开窗的方法提取目标信息;
第四步,目标回波能量补偿信息重构
根据间歇收发参数,得到脉冲压缩后信号的能量补偿值,通过能量补偿,重构目标探测信息。
2.根据权利要求1所述的微波暗室中脉冲体制雷达目标探测模拟方法,其特征在于:第一步所述间歇收发方式下,收发交替工作过程可以抽象为一方波信号p(t)对雷达脉冲信号进行控制,p(t)可表示为
p ( t ) = r e c t ( t τ ) * Σ n = - ∞ ∞ δ ( t - nT s ) - - - ( 1 )
其中τ为发射通道工作时段,Ts为间歇收发周期,δ(·)为单位冲击函数,n表示第n个冲击脉冲函数;
根据傅氏变换性质,p(t)的频域形式为
P ( f ) = τf s Σ n = - ∞ n = + ∞ sin c ( nf s τ ) δ ( f - nf s ) - - - ( 2 )
其中,fs=1/Ts,sinc(x)=sin(πx)/(πx);可以发现,间歇收发控制信号的频谱为p(t)各时域谐波的组合,且主峰处幅值由收发参数τfs决定,并决定目标信息重构的能量补偿值;
通过精确的时延和波形控制可以保证实际发射信号的相位连续性和幅度一致性,从而确保收发过程等效为方波信号p(t)与雷达发射脉冲信号s0(t)的乘积;
目标散射过程可等效为一线性系统与激励信号相卷积;设目标与雷达相对距离为R,散射强度为σ0,径向运动速度为v,目标在雷达视线方向投影长度为L;雷达目标时域冲激响应函数为h(t)=σ0δ(t-Δt);假设雷达发射脉冲信号为一线性调频脉冲信号;从而,间歇收发条件下,目标回波信号去载频后可表示为:
s 2 ( t ) = A u ( t - Δ t ) exp ( - j 2 πf 0 Δ t ) [ r e c t ( t τ ) * Σ n = - ∞ ∞ δ ( t - nT s ) ] - - - ( 3 )
其中,A为回波幅度,可由回波功率得到,Pt为发射功率,G为天线收发增益,λ为波长,Δt=2(R-vt)/C,C为电磁波传播速度;rect(·)为矩形窗函数,Tp为脉宽,u(t)=exp(jπγt2)为复包络信号,为单位虚数,f0为中心频率,γ为线性调频率;可以发现,目标回波信号由一系列子脉冲信号构成,各段子脉冲宽度与发射通道工作时段τ相等,子脉冲带宽BΔ=γτ。
3.根据权利要求1所述的微波暗室中脉冲体制雷达目标探测模拟方法,其特征在于:第二步所述目标回波恢复,具体方法如下:s2(t)对应频域形式为
S 2 ( f ) = Aτf s Σ n = - ∞ n = + ∞ sin c ( nf s τ ) U ( f - nf s ) exp ( - j 2 π ( f + f 0 - nf s ) Δ t ) - - - ( 4 )
其中sinc(·)为辛克函数,U(f)为u(t)的频谱;
根据线性调频LFM信号带宽,设计低通滤波器截止频率刚好覆盖信号带宽,从而间歇收发之后,位于原始LFM信号带宽之外的频谱分量将被有效滤除,进一步可以得到信号时域波形;考虑W(f)为低通滤波器的频谱,理想情况下,滤波器的频谱满足如下条件
W ( f ) = 1 , f < = B 0 , f > B - - - ( 5 )
从而,经过低通滤波之后,可得
S3(f)=W(f)S2(f) (6)
进一步,可以获得经过低通滤波之后信号的时域形式。
4.根据权利要求1所述的微波暗室中脉冲体制雷达目标探测模拟方法,其特征在于:第三步所述脉冲压缩获取目标信息,具体方法如下:在接收机中,根据发射信号可得幅度归一化匹配滤波器频率响应为
H ( f ) = r e c t ( f B ) &CenterDot; exp ( j &pi; f 2 &gamma; - j &pi; 4 ) - - - ( 7 )
其中f是频率分量,B是雷达信号带宽;
根据匹配滤波器特点,W(f)在LFM信号带宽内为1,带宽外为0,与H(f)基本相同,因此相乘后可以忽略W(f);从而,目标回波信号经匹配滤波并取模后输出形式可以表示为
| s h ( t ) | = | I F F T &lsqb; S 3 ( f ) H ( f ) &rsqb; | = | AB&tau;f s &Sigma; n = - &infin; n = + &infin; sin c ( nf s &tau; ) sin c ( B ( t + nf s &gamma; - &Delta; t ) ) | - - - ( 8 )
其中|·|表示取模运算,IFFT表示逆傅里叶变换;对|sh(t)|的主峰位置(n=0)进行时域开窗处理可滤出主峰,完成目标信息重构;需要条件分析如下:
|sh(t)|相邻两个辛克函数尖峰相距
&Delta; R = cf s 2 &gamma; - - - ( 9 )
要保证匹配滤波后相邻峰值不发生混叠,则需满足
ΔR>L (10)
此外要满足目标回波信号不遮挡和完全接收条件,还须有
&tau; &le; &Delta; t &tau; + &Delta; t + &Delta; &tau; &le; T s - - - ( 11 )
其中表示信号遍历目标耗时;
从而,约束条件可写为
&tau; &le; 2 R C &tau; + 2 ( R + L ) C &le; T s < CT p 2 B L - - - ( 12 )
一方面Ts越小(fs越大),则相邻两个峰值距离越远,越利于目标信息重构,一旦Ts足够小,在距离门内将只能观察到实际目标处的尖峰;另一方面Ts越大则信号辐射能量越小,越不易恢复目标信息;
根据上述分析,当间歇收发频率小于信号带宽时,需要在匹配滤波之后,通过合理的开窗,提取脉冲压缩后的主峰,重构目标信息;开窗的目的主要是为了消除n≠0的峰值;在理想情况下,开窗类似于用矩形通断截取脉压主峰,通断时间对应的距离宽要小于相邻峰的距离,从而
因此,截断之后间歇收发回波的匹配滤波输出为
其中,回波延时Δt可以根据加窗前式(4)中的主峰位置估计得到;这样,就可以利用“间歇”收发方法解决空间有限条件下的长脉冲信号目标探测问题;此外,根据式(9),当间歇收发周期Ts与完整信号带宽相当时,相邻虚假峰将无法在距离门内观测到;此时,不采用时域截断即可有效获取目标信息。
5.根据权利要求1所述的微波暗室中脉冲体制雷达目标探测模拟方法,其特征在于:第四步所述目标回波能量补偿与信息重构,方法如下:
根据间歇控制信号p(t)的频域形式以及所得间歇收发回波匹配滤波输出,可以发现,对于实际目标信息而言,间歇收发匹配滤波输出的主峰(n=0)反映了真实目标信息;但是由于发射信号的通断,损失了部分信号能量,因此所得目标信息需要进行能量补偿;根据式(8),取n=0处的脉压输出信息,并对信号幅度补偿1/τfs,即可精确的恢复得到目标回波的匹配滤波输出,从而实现暗室中脉冲雷达信号的目标探测模拟过程;恢复后的脉压输出表达式为
CN201710081095.8A 2017-02-15 2017-02-15 微波暗室中脉冲体制雷达目标探测模拟方法 Active CN106908771B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710081095.8A CN106908771B (zh) 2017-02-15 2017-02-15 微波暗室中脉冲体制雷达目标探测模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710081095.8A CN106908771B (zh) 2017-02-15 2017-02-15 微波暗室中脉冲体制雷达目标探测模拟方法

Publications (2)

Publication Number Publication Date
CN106908771A true CN106908771A (zh) 2017-06-30
CN106908771B CN106908771B (zh) 2019-10-11

Family

ID=59208660

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710081095.8A Active CN106908771B (zh) 2017-02-15 2017-02-15 微波暗室中脉冲体制雷达目标探测模拟方法

Country Status (1)

Country Link
CN (1) CN106908771B (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108169746A (zh) * 2017-12-21 2018-06-15 南京理工大学 线性调频脉冲半主动雷达导引头信号处理方法
CN108880666A (zh) * 2018-03-26 2018-11-23 北京大学 一种基于微波光子技术的串行通信分析仪及其波形重构方法
CN110018449A (zh) * 2017-08-31 2019-07-16 成都玖锦科技有限公司 一种利用包络信息的信号合成方法
CN110095757A (zh) * 2018-01-29 2019-08-06 株式会社万都 用于补偿雷达的天线的回波损耗的装置和方法及雷达设备
CN111123256A (zh) * 2019-12-06 2020-05-08 中国人民解放军国防科技大学 微波暗室中脉冲雷达进动目标微动特征提取方法
CN111712725A (zh) * 2018-12-26 2020-09-25 北京航迹科技有限公司 激光雷达测距的多脉冲融合分析
CN112731308A (zh) * 2020-12-21 2021-04-30 北京机电工程研究所 一种自适应低频有源对消雷达隐身实现方法
WO2021189287A1 (zh) * 2020-03-25 2021-09-30 深圳市汇顶科技股份有限公司 基于飞行时间的测距方法和测距系统
CN113791395A (zh) * 2021-08-04 2021-12-14 中国人民解放军国防科技大学 一种基于正交间歇采样调制雷达信号的同时极化测量方法
CN114200418A (zh) * 2021-11-10 2022-03-18 中国人民解放军国防科技大学 一种脉内多载频信号收发设计及回波重构方法
CN117031444A (zh) * 2023-10-10 2023-11-10 北京航空航天大学 一种高精度快速相位式激光测距方法
CN118191775A (zh) * 2024-03-14 2024-06-14 中国工程物理研究院应用电子学研究所 一种高功率微波后门强耦合参数提取方法及装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102323570A (zh) * 2011-05-24 2012-01-18 中国人民解放军国防科学技术大学 一种雷达目标回波信号模拟器的幅相特性估计方法
CN105137404A (zh) * 2015-08-13 2015-12-09 北京理工大学 基于预脉冲间处理的雷达压缩采样方法及系统
CN105866755A (zh) * 2016-05-30 2016-08-17 中国人民解放军国防科学技术大学 微波暗室内脉冲体制雷达目标回波信息重构方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102323570A (zh) * 2011-05-24 2012-01-18 中国人民解放军国防科学技术大学 一种雷达目标回波信号模拟器的幅相特性估计方法
CN105137404A (zh) * 2015-08-13 2015-12-09 北京理工大学 基于预脉冲间处理的雷达压缩采样方法及系统
CN105866755A (zh) * 2016-05-30 2016-08-17 中国人民解放军国防科学技术大学 微波暗室内脉冲体制雷达目标回波信息重构方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
黄光才等: "对双基地雷达的间歇采样转发干扰效果分析", 《航天电子对抗》 *

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110018449A (zh) * 2017-08-31 2019-07-16 成都玖锦科技有限公司 一种利用包络信息的信号合成方法
CN110018449B (zh) * 2017-08-31 2022-11-22 成都玖锦科技有限公司 一种利用包络信息的信号合成方法
CN108169746B (zh) * 2017-12-21 2021-09-21 南京理工大学 线性调频脉冲半主动雷达导引头信号处理方法
CN108169746A (zh) * 2017-12-21 2018-06-15 南京理工大学 线性调频脉冲半主动雷达导引头信号处理方法
CN110095757A (zh) * 2018-01-29 2019-08-06 株式会社万都 用于补偿雷达的天线的回波损耗的装置和方法及雷达设备
CN110095757B (zh) * 2018-01-29 2023-12-22 汉拿科锐动电子股份公司 用于补偿雷达的天线的回波损耗的装置和方法及雷达设备
CN108880666A (zh) * 2018-03-26 2018-11-23 北京大学 一种基于微波光子技术的串行通信分析仪及其波形重构方法
CN108880666B (zh) * 2018-03-26 2020-10-09 北京大学 基于微波光子技术的串行通信分析仪及其波形重构方法
CN111712725A (zh) * 2018-12-26 2020-09-25 北京航迹科技有限公司 激光雷达测距的多脉冲融合分析
CN111712725B (zh) * 2018-12-26 2023-11-17 广州沃芽来得灵科技有限公司 激光雷达测距的多脉冲融合分析
CN111123256A (zh) * 2019-12-06 2020-05-08 中国人民解放军国防科技大学 微波暗室中脉冲雷达进动目标微动特征提取方法
WO2021189287A1 (zh) * 2020-03-25 2021-09-30 深圳市汇顶科技股份有限公司 基于飞行时间的测距方法和测距系统
CN112731308A (zh) * 2020-12-21 2021-04-30 北京机电工程研究所 一种自适应低频有源对消雷达隐身实现方法
CN113791395B (zh) * 2021-08-04 2023-11-17 中国人民解放军国防科技大学 一种基于正交间歇采样调制雷达信号的同时极化测量方法
CN113791395A (zh) * 2021-08-04 2021-12-14 中国人民解放军国防科技大学 一种基于正交间歇采样调制雷达信号的同时极化测量方法
CN114200418A (zh) * 2021-11-10 2022-03-18 中国人民解放军国防科技大学 一种脉内多载频信号收发设计及回波重构方法
CN117031444A (zh) * 2023-10-10 2023-11-10 北京航空航天大学 一种高精度快速相位式激光测距方法
CN117031444B (zh) * 2023-10-10 2023-12-12 北京航空航天大学 一种高精度快速相位式激光测距方法
CN118191775A (zh) * 2024-03-14 2024-06-14 中国工程物理研究院应用电子学研究所 一种高功率微波后门强耦合参数提取方法及装置

Also Published As

Publication number Publication date
CN106908771B (zh) 2019-10-11

Similar Documents

Publication Publication Date Title
CN106908771A (zh) 微波暗室中脉冲体制雷达目标探测模拟方法
Smith et al. Multistatic micro-Doppler radar signatures of personnel targets
CN105866755B (zh) 微波暗室内脉冲体制雷达目标回波信息重构方法
CN103412295B (zh) 基于回波精确模型的高速机动弱目标检测方法
Jakabosky et al. Waveform design and receive processing for nonrecurrent nonlinear FMCW radar
CN107192991B (zh) 一种射频辐射式仿真宽带脉冲雷达成像等效模拟方法
CN106644030A (zh) 一种基于多普勒雷达的非接触式振动测量方法
CN106707258A (zh) 一种非高斯背景下微动目标多参数估计方法
CN101872014A (zh) 基于改进的cosgo-cfar的目标信号检测方法
CN104160296A (zh) 目标检测方法
CN106383340A (zh) 一种随机脉冲初始相位雷达的速度假目标识别方法
CN109975775B (zh) 一种雷达回波半实测数据仿真方法
WO2014201795A1 (zh) 一种多扫频雷达信号的产生及压缩方法
CN109782250A (zh) 基于有限新息率采样的雷达目标参数提取方法
CN106291489A (zh) 适用于多种发射信号波形的合成孔径雷达回波仿真方法
CN107271980B (zh) 一种对间歇调制信号的分段匹配滤波处理方法
CN106680791A (zh) 基于宽带扫频数据的雷达回波仿真方法
Higgins et al. Passive multistatic radar experiment using WiMAX signals of opportunity. Part 1: Signal processing
CN104062641B (zh) 一种ofdm雷达距离像处理的新方法
CN106019244A (zh) 基于正负斜率lfm波形正交极化雷达的抗移频干扰方法
CN104237864B (zh) 基于匹配模糊函数的距离扩展目标检测方法
CN103412289A (zh) 基于模拟sinc核采样的雷达高分辨距离成像方法
Liu et al. Transmitter power comparison between incoherent noise and coherent deceptive jamming against synthetic aperture radar
Adjrad et al. Experimental vital signs estimation using commercially available IR-UWB radar
CN108801660B (zh) 一种高速列车紧急制动性能的实验室检测装置及方法

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