脉冲雷达系统设计与Matlab仿真【附仿真代码】

【摘要】

  • 参考书目:《雷达系统设计Matlab仿真》
  • 其他参考资料:屈卫《雷达系统与原理》视频课
  • 本文主题:雷达分类、作用距离、分辨率、雷达方程、脉冲积累、噪声系数
  • 实验:设计一个相对完整的脉冲雷达系统,实现对目标的搜索和探测,并用matlab进行仿真

文章目录

  • 1 雷达分类
  • 2 作用距离
  • 3 距离分辨率
  • 4 多普勒频移
  • 5 雷达方程
  • 7 脉冲积累
  • 8 雷达损失
  • 9 雷达设计实例与Matlab仿真
    • 9.2 需求描述
    • 9.3 雷达设计
    • 9.4 matlab仿真程序
  • 附录1:脉冲的信号宽度


1 雷达分类

  • 按照平台:地基、记载、天基、舰载雷达系统。
  • 按照任务和功能:气象、捕获和搜索、跟踪、边扫描边跟踪、获利控制、预警、超视距、地形跟随、地形回避。
  • 相控阵雷达使用相控阵天线,通常称为多功能(多模式)雷达。
  • 按工作频段:
波段符号 频率/GHz 说明 图片
HF高频 0.003-0.03 使用电离层反射的电磁波检测地平线以下的目标
VHF甚高频 0.03-0.3 甚远距离预警雷达(EWR),波长大、测量距离远,雷达孔径大! 对空远距离跟踪战略轰炸机、弹道导弹等。此类地基雷达发射功率一般在10MW以上。 脉冲雷达系统设计与Matlab仿真【附仿真代码】_第1张图片
UHf超高频 0.3-1(基本与P波段0.23-1.0重合) 同VHF,典型型号:美军AN/FPS-115,工作频率420MHz-450MHz,探测距离一般为4800公里。 脉冲雷达系统设计与Matlab仿真【附仿真代码】_第2张图片
L波段 1.0-2.0 主要用于地基和舰载雷达, 远程 军事和空中交通控制搜索行动
S波段 2.0-4.0 大部分地基和舰载 中程 雷达。典型型号: “ 宙 斯 盾 ” 舰 载 无源相控阵雷达AN/SPY-1D,3.1-3.5GHz,对空探测距离约400km。 脉冲雷达系统设计与Matlab仿真【附仿真代码】_第3张图片
C波段 4.0-8.0 大半部分气象雷达,中程搜索和火力控制军用雷达以及测量雷达。典型型号:“爱国者”导弹制导雷达AN/MPQ-53,对 1 m 2 1\rm m^2 1m2RCS的探测距离约为160km。 脉冲雷达系统设计与Matlab仿真【附仿真代码】_第4张图片
X波段 8.0-12.5 用于天线大小受物理限制的雷达系统,包括大部分军用多模式雷达。典型型号:萨德系统配套的AN/TPY-2,X波段(9.5GMz)。 脉冲雷达系统设计与Matlab仿真【附仿真代码】_第5张图片
Ku波段 12.5-18.0 大气和气象衰减严重,多用于近程,如警用交通雷达、近程地形回避雷达,地形跟随
K波段 18.0-26.5 同Ku
Ka波段 26.5-40 同Ku
MMW毫米波 通常>34.0 用于非常近距离的射频搜索和实验雷达

2 作用距离

脉冲雷达系统设计与Matlab仿真【附仿真代码】_第6张图片
图1.1是一个简化的脉冲雷达框图,收发开关允许天线在发射和接收模式之间切换。目标距离天线的距离用 R R R表示,信号从发出到经过目标反射再回到天线经过的时间记为 Δ t \Delta t Δt,则
R = c Δ t 2 R=\frac{c\Delta t}{2} R=2cΔt其中 c = 3 × 1 0 8 m / s c=3\times10^8\rm m/s c=3×108m/s为电磁波在真空中的传播速度。
脉冲雷达系统设计与Matlab仿真【附仿真代码】_第7张图片
图1.2中的几个名词:
脉冲间周期(IPP)也称为脉冲重复间隔(PRI),一般用符合 T T T表示,它的倒数称为脉冲重复频率(PRF),用 f r f_r fr表示,则:
f r = 1 P R I = 1 T f_r=\frac{1}{PRI}=\frac{1}{T} fr=PRI1=T1
脉冲宽度记为 τ \tau τ,每个PRI期间,雷达只发射 τ \tau τ秒能量,其余时间侦听目标反射信号,站空因子 d t = τ / T d_t=\tau/T dt=τ/T,雷达平均功率
P a v = P t × d t P_{av}=P_t\times d_t Pav=Pt×dt
其中 P t P_t Pt表示雷达发射峰值功率,脉冲能量 E p = P t τ = P a v T = P a v / f r E_p=P_t\tau=P_{av}T=P_{av}/f_r Ep=Ptτ=PavT=Pav/fr
脉冲雷达系统设计与Matlab仿真【附仿真代码】_第8张图片
脉冲测距,距离可能存在模糊。如图1.3,第一个脉冲发出后,先后遇到了远近两个目标后产生两个回波,第二个回波还没到达接收机,雷达又发出了第二个脉冲,此时第二个回波的时延应该是 T + Δ t T+\Delta t T+Δt,因此 R 2 = c ( T + Δ t ) / 2 R_2=c(T+\Delta t)/2 R2=c(T+Δt)/2

3 距离分辨率

雷达分辨率一般用 Δ R \Delta R ΔR表示,描述了雷达区分位置非常靠近的不同目标的能力。很明显
Δ R = c τ 2 = c 2 B \Delta R=\frac{c\tau}{2}=\frac{c}{2B} ΔR=2cτ=2Bc
其中 B = 1 / τ B=1/\tau B=1/τ,表示带宽。

4 多普勒频移

目标运动的相对速度为 v v v,脉冲
靠近的目标: f d = 2 v λ f_d=\frac{2v}{\lambda} fd=λ2v,
远去的目标: f d = − 2 v λ fd=-\frac{2v}{\lambda} fd=λ2v

5 雷达方程

对全向天线来说,设其峰值发射功率为 P t P_t Pt,距离天线 R R R处的功率密度为
P D = P t 4 π R 2 (1) P_D=\frac{P_t}{4\pi R^2}\tag{1} PD=4πR2Pt(1)
方向性天线可以将能量定向到某个方向,假设定向天线的增益为 G G G,其有效孔径为 A e A_e Ae,二者的关系为 G = 4 π A e λ 2 G=\frac{4\pi A_e}{\lambda^2} G=λ24πAe目标的雷达反射截面积(RCS)为 σ \sigma σ,则天线接收的回波信号功率为
P r = P t G σ ( 4 π R 2 ) 2 A e P_r=\frac{P_tG\sigma}{(4\pi R^2)^2}A_e Pr=(4πR2)2PtGσAe
单基地雷达,收发用同一个天线,结合公式(1),则
P r = P t G 2 λ 2 σ ( 4 π ) 3 R 4 , P r = P t A 2 σ 4 π λ 2 R 4 P_r=\frac{P_tG^2\lambda^2\sigma}{(4\pi)^3 R^4}, P_r=\frac{P_tA^2\sigma}{4\pi\lambda^2 R^4} Pr=(4π)3R4PtG2λ2σ,Pr=4πλ2R4PtA2σ
对应的最大雷达探测距离为
R max ⁡ = ( P t G 2 λ 2 σ ( 4 π ) 3 S min ⁡ ) 1 4 , R max ⁡ = ( P t A 2 σ ( 4 π ) 3 λ 2 S min ⁡ ) 1 4 R_{\max}=\left(\frac{P_tG^2\lambda^2\sigma}{(4\pi)^3 S_{\min}}\right)^{\frac{1}{4}}, R_{\max}=\left(\frac{P_tA^2\sigma}{(4\pi)^3 \lambda^2 S_{\min}}\right)^{\frac{1}{4}} Rmax=((4π)3SminPtG2λ2σ)41,Rmax=((4π)3λ2SminPtA2σ)41
其中 S min ⁡ S_{\min} Smin表示最小可检测信号功率。定义噪声系数(参考这篇文章)
F = S N R i S N R o = S i / N i S o / N o F=\frac{SNR_i}{SNR_o}=\frac{S_i/N_i}{S_o/N_o} F=SNRoSNRi=So/NoSi/Ni
且输入噪声为热噪声
N i = k T B N_i=kTB Ni=kTB即玻尔兹曼常数、温度、带宽。因此
S min ⁡ = k T B F ( S N R ) o min ⁡ S_{\min}=kTBF(SNR)_{o\min} Smin=kTBF(SNR)omin
进一步,引入雷达系统损耗,记为 L L L,则
S N R o = P t G 2 λ 2 σ ( 4 π ) 3 R 4 k T B F L (1.56) SNR_o=\frac{P_tG^2\lambda^2\sigma}{(4\pi)^3R^4kTBFL}\tag{1.56} SNRo=(4π)3R4kTBFLPtG2λ2σ(1.56)

一部雷达部署后,(1.56)里面的 G G G, λ \lambda λ, k k k, T T T, B B B, F F F, L L L什么的都是确定不变的。接收机输出端的信噪比(SNR)和目标的RCS( σ \sigma σ)、距离( R R R),峰值功率( P t P_t Pt)之间遵循什么样的规律呢?

根据式(1.56),图1.12(a)画出了峰值功率功率为 P t P_t Pt=1.5MW,不同大小的目标,雷达接收机输出信噪比SNR随目标距离 R R R的变化的曲线。可以 看出,目标越大,在相同距离时,接收机输出信噪比越高。另一方面,假设检测目标所需的最小信噪比为10dB,则三种大小的目标分别对应的探测距离为60km、108km、181km。

图1.12(b)画出了 σ = 0.1 m 2 \sigma=0.1\rm m^2 σ=0.1m2,峰值功率功率分别为 P t P_t Pt=2.16MW、1.5MW、0.6MW时,雷达接收机输出信噪比SNR随目标距离 R R R的变化的曲线。可以 看出,功率越高,看得越远。

脉冲雷达系统设计与Matlab仿真【附仿真代码】_第9张图片 脉冲雷达系统设计与Matlab仿真【附仿真代码】_第10张图片
图1.12(a)不同RCS下,SNR相对检测距离的曲线 图1.12(b)不同峰值功率下,SNR相对检测距离的曲线

点赞、收藏、并评论“仿真代码”获取以上曲线图的仿真代码。

7 脉冲积累

当一个目标在单次扫描期间位于雷达波束内时,它可能反射几个脉冲C。通过把一个给定目标在单次打描期间反射的所有脉冲的回波相加,雷达的灵敏度 (SNR) 会增加。

单个脉冲的信噪比记为 S N R 1 SNR_1 SNR1 n p n_p np个脉冲进行相干积累后的信噪比
S N R C I = n p S N R 1 SNR_{CI}=n_pSNR_1 SNRCI=npSNR1
非相干积累后的信噪比
S N R N C I = n p S N R 1 2 1 + S N R 1 SNR_{NCI}=n_p\frac{SNR_1^2}{1+SNR_1} SNRNCI=np1+SNR1SNR12

8 雷达损失

发射和接收损失(波导损失),1-2dB。

天线方向图损失,由于天线没有正对着目标。
脉冲雷达系统设计与Matlab仿真【附仿真代码】_第11张图片

此外还有大气损失、跌价损失、处理损失等。

9 雷达设计实例与Matlab仿真

9.2 需求描述

脉冲雷达系统设计与Matlab仿真【附仿真代码】_第12张图片

9.3 雷达设计

脉冲雷达系统设计与Matlab仿真【附仿真代码】_第13张图片
脉冲雷达系统设计与Matlab仿真【附仿真代码】_第14张图片
脉冲雷达系统设计与Matlab仿真【附仿真代码】_第15张图片
脉冲雷达系统设计与Matlab仿真【附仿真代码】_第16张图片
脉冲雷达系统设计与Matlab仿真【附仿真代码】_第17张图片

9.4 matlab仿真程序

点赞、收藏、并评论“仿真代码”获取以上曲线图的仿真代码。


附录1:脉冲的信号宽度

关于脉冲宽度和它的带宽之间的关系,可以看下面的推导(樊昌信《通信原理》(第六版):

脉冲雷达系统设计与Matlab仿真【附仿真代码】_第18张图片

你可能感兴趣的:(雷达与通信信号处理,雷达原理,雷达方程,脉冲积累,噪声系数,雷达Matlab仿真,1024程序员节)