简介雷达基础知识又名雷达原理是一份面向电子工程、通信工程及雷达相关专业初学者的幻灯片课件也适用于教师备课与自习巩固。课件以雷达定义和系统组成为切入点完整介绍发射机、接收机、天线、信号处理等核心模块并借助组成框图帮助理解整体工作流程。在测距方面详细讲解电磁波直线传播与反射测距原理给出回波延迟时间求距离的计算公式进一步推导最大不模糊距离、最小测量距离并通过实例演示脉冲重复频率与距离计算使抽象概念变得直观。在分辨力方面系统阐述距离分辨力和角度分辨力的定义、判定准则结合天线方向图说明方位角与俯仰角的获取方法同时配备战术雷达横向距离计算例题强化公式应用能力。资源共一个幻灯片文件大小约一点四兆字节幻灯片结构完整、公式与图例并重便于按章节系统学习雷达参数体系已有四千一百一十九人学习下载是入门雷达原理的高性价比参考资料。1. 雷达原理的第一课雷达在“看”什么把一辆静止的卡车和一辆移动中的自行车同时放在雷达面前回波强度差出几十倍但不少雷达算法最后跟踪的往往是自行车而不是卡车。这不是算法 bug而是雷达的工作原理决定的雷达并不像摄像头那样“看见”物体的轮廓它测量的是回波相对发射信号的延迟、频率和相位变化。延迟给出距离频率偏移给出径向速度不同接收通道之间的相位差给出角度。这三个物理量就是雷达原理里所有公式展开的原点。所以学雷达基础不能用“发射电磁波、接收回波”这九个字一笔带过。真正要建立的是从物理量到数据量的换算框架目标距离 1km对应电磁波往返约 6.7μs目标以 10m/s 径向运动在 10GHz 频段产生约 667Hz 的多普勒频移。这些数字是后面一切工程判断的标尺。这篇内容适合做自动驾驶感知、气象观测、安防监视、无人机避障以及正在转行做无线电和信号处理的工程师。雷达原理不只是一门电磁场课它本质上是“用有限的时间、带宽和功率从噪声里把想要的信息榨出来”的系统工程这恰恰是 IT 背景的人最容易上手的角度。2. 雷达测距原理与雷达方程功率怎么决定最远能看到哪2.1 测距的起点电磁波往返时间脉冲雷达的工作方式说起来很简单每隔一段固定时间发射一个短脉冲然后等待回波。如果脉冲发射后经过时间 t 收到回波目标距离就是 R c·t/2系数 1/2 是因为电磁波走了个来回。光速 c 取 299792458 m/s1μs 的时延对应 150m 距离。这里的核心结论是测距本质上是在测时间。时间测量的精度、时钟源的稳定性、接收机对回波到达时刻的判决方式直接决定距离精度。工程上通常不用回波前沿做判决因为前沿容易被噪声扰动而是用匹配滤波输出后的峰值位置。峰值判决更稳健但会引入至多半个采样间隔的固定偏置这在后文的脉冲压缩里还会再提到。2.2 雷达方程里的每个参数都在回答工程问题雷达方程是雷达原理里最不该跳过的一个公式它把发射、传播、目标、接收四个环节压缩成一行Pr Pt·G²·λ²·σ / ((4π)³·R⁴)逐项拆开看Pt 是峰值发射功率G 是天线增益收发共用一副天线所以出现 G²λ 是波长决定天线尺寸与大气损耗σ 是目标雷达截面积 RCS单位是平方米它不等于目标的物理面积R⁴ 是双向空间衰减距离翻一倍回波功率掉到原来的 1/16。这四个变量里R 是雷达无法控制的其余全是设计自由度。工程上更常用 dB 形式做预算。一个典型算例X 波段 10GHzλ0.03m峰值功率 10W天线增益 30dB目标 RCS 为 1m²距离 5km。套公式算下来接收功率约 -111dBm而 1MHz 带宽、3dB 噪声系数下的热噪声底约 -111dBm。也就是说单脉冲 SNR 接近 0dB刚好卡在检测边缘。这个数字说明一个工程事实雷达要看得远不能只靠提高峰值功率更要靠相干积累把多个脉冲的能量叠起来。下表是一组 X 波段脉冲雷达的常见参数后面章节的讨论都以它为基础参数符号典型值说明峰值功率Pt10W车载/ 1kW地面监视受器件与散热限制中心频率f010GHzX 波段λ3cm脉宽τ1μs决定基础距离分辨信号带宽B50~200MHz脉内调频后决定实际分辨天线增益G30dB收发共用脉冲重复频率PRF1~10kHz决定模糊边界噪声系数F3dB接收机前端质量积累脉冲数N64~128相干积累增益写个小脚本算一下 SNR 随距离的变化比手算快得多import numpy as np c 299_792_458.0 f0 10e9 lam c / f0 Pt 10.0 # 峰值功率单位 W G 10 ** (30 / 10) # 天线增益 30 dB sigma 1.0 # RCS单位 m^2 kB 1.380649e-23 T0 290.0 # 参考噪声温度单位 K B 1e6 # 匹配滤波带宽 1 MHz NF 10 ** (3 / 10) # 噪声系数 3 dB R np.linspace(1e3, 20e3, 200) Pr Pt * G**2 * lam**2 * sigma / ((4 * np.pi)**3 * R**4) Pn kB * T0 * B * NF SNR 10 * np.log10(Pr / Pn) for r in [3000, 5000, 10000]: pr Pt * G**2 * lam**2 * sigma / ((4 * np.pi)**3 * r**4) print(fR {r/1000:.0f} km, SNR {10*np.log10(pr / Pn):.1f} dB)这段代码把雷达方程拆成线性计算再用 log 转成 SNR。运行结果里 5km 处 SNR 在 0dB 附近与手算一致。脚本里的可调参数是 Pt、G、NF 和 B把 NF 从 3dB 降到 1.5dBSNR 能改善 1.5dB想靠提高功率补回 3dB得把 Pt 翻倍。所以接收机噪声系数往往是雷达链路里最值得花钱的地方。2.3 距离分辨率关键不在脉宽在带宽测距精度和距离分辨率是两个概念但常被混淆。距离分辨率指雷达能把两个相距多近的目标区分开的能力。对单载频矩形脉冲分辨率约等于 δR c·τ/21μs 脉宽对应 150m。也就是说两辆车相距 150m 以内在普通脉冲雷达上看就是一个目标。但现代雷达几乎都采用脉内线性调频LFM发射脉冲被展宽到 1μs但内部频率在 200MHz 带宽内扫描。接收端做匹配滤波后等效脉冲宽度被压缩到约 1/B距离分辨率变成 δR c/(2B)200MHz 带宽对应 0.75m。这就是雷达原理里最常被新手绕晕的一点分辨率由带宽决定不是由脉宽决定。压窄脉宽能提高分辨率但峰值功率要求急剧上升展宽带宽再压缩成本在信号处理链路而不是发射机。现代雷达系统几乎都选择后者这也是下一章要展开的波形设计问题。3. 脉冲雷达的波形与模糊重频、盲速和距离模糊怎么选3.1 脉冲重复频率决定了能看多远、测多快脉冲雷达每发射一个脉冲要在下一个脉冲发射前收到上一脉冲的回波否则距离就模糊了。因此最大不模糊距离是 R_unamb c/(2·PRF)。与此同时对回波做多普勒处理时脉冲采样定理限制了不模糊速度范围最大不模糊速度是 v_unamb λ·PRF/2实际可测范围是 ±v_unamb/2。这两个公式放在一起就是一个经典矛盾PRF 高看得近但测速范围大PRF 低看得远但速度容易模糊。以第 2 章那组 X 波段参数为例PRF最大不模糊距离最大不模糊速度±典型场景1kHz150km±7.5m/s远程警戒5kHz30km±37.5m/s中程监视20kHz7.5km±150m/s高速目标跟踪注意 1kHz 那一行±7.5m/s 的测速范围连正常行驶的汽车都覆盖不了。所以远程雷达不是单纯选低 PRF 就完事速度模糊必须配合解模糊算法一起处理。3.2 重频参差与双重频解模糊的常见工程做法既然低 PRF 测距清晰、高 PRF 测速清晰常见的工程做法就是“两个都不够就两个都发”。典型方案是重频参差交替使用两个相近但不相同的 PRF各自得到带模糊的测量值再通过余数定理思想在候选距离里找一致解。下面给出一个双重频解距离模糊的简化实现def resolve_range(r1, r2, prf1, prf2, max_r): c 299_792_458.0 ramb1 c / (2 * prf1) ramb2 c / (2 * prf2) tol min(ramb1, ramb2) * 0.25 candidates [] for n1 in range(int(max_r / ramb1) 1): d1 r1 n1 * ramb1 for n2 in range(int(max_r / ramb2) 1): d2 r2 n2 * ramb2 if abs(d1 - d2) tol: candidates.append((d1 d2) / 2) return candidates # 例真实距离 82.3kmPRF11kHz(不模糊150km)PRF21.1kHz # 两个雷达各自量到的表观距离需要由前端给出 print(resolve_range(82.3e3, 82.3e3, 1000, 1100, 300e3))这里的核心思想是遍历模糊重数 n1、n2把两个带模糊的距离各自外推找出互相吻合的候选值。tol 的取值要大于距离测量误差但又不能大到把多个候选都吞进来。实际系统里往往加入第三个 PRF 来消除余数定理的二次模糊三重频的正确解概率会显著高于双重频。3.3 距离-多普勒耦合的边界模糊函数不是考试概念LFM 信号有一个绕不开的性质距离-多普勒耦合。目标的多普勒频移 f_d 会让压缩后的峰值在时延轴上偏移一个量偏移量约等于 f_d 除以调频斜率 k单位 Hz/s。目标速度越快测距偏差越大。对 200MHz/1ms 的扫频参数10m/s 目标大约造成 1m 量级的距离偏移这在近距高精度场景里已经不能忽略。模糊函数描述的是“波形对距离和速度联合分辨的能力边界”它回答的问题是给定一个波形在距离-速度二维平面上能分辨到多细。工程上选波形本质上就是选模糊函数的形状LFM 是倾斜的刀刃形相位编码是图钉形各自适合不同的场景。理解到这一层再看后面 MTD 和 FMCW 的设计就能明白每个参数为什么出现在那个位置。4. 从回波到点迹脉冲压缩、MTD 与 CFAR 的实现链路4.1 匹配滤波与脉冲压缩把 1μs 变成 1.5m回波进接收机后第一个关键步骤是匹配滤波。匹配滤波的理论结论是在白色高斯噪声背景下使输出 SNR 最大的线性滤波器是发射信号的共轭时间翻转。对 LFM 信号做匹配滤波输出波形被压缩成宽度约 1/B 的窄峰同时获得 B·τ 倍的 SNR 增益。1μs 脉宽、200MHz 带宽的 LFM压缩比是 200即 23dB 的积累增益距离分辨率从 150m 提升到 0.75m。这一步实现起来并不神秘频域做一次复数相乘加 IFFT 即可import numpy as np fs 200e6 T 1e-6 B 200e6 N int(fs * T) t np.arange(N) / fs k B / T lfm_tx np.exp(1j * 2 * np.pi * (k / 2) * t**2) # 发射信号 delay 100 # 回波延迟 100 个采样点 lfm_rx np.zeros(N, dtypecomplex) lfm_rx[delay:] lfm_tx[:N - delay] # 模拟回波 H np.conj(np.fft.fft(lfm_tx)) # 匹配滤波器 pulse_comp np.fft.ifft(np.fft.fft(lfm_rx) * H) # 频域滤波 peak np.argmax(np.abs(pulse_comp)) print(f峰值位置: {peak}, 理论延迟: {delay}, 分辨率: {c / (2*B):.2f} m)注意匹配滤波器用的是发射信号的共轭谱回波经过它之后输出峰值的位置就对应目标时延。这里的关键参数是采样率 fs 和带宽 B峰值定位精度受采样率限制如果 fs 只取 200MHz逐点峰值会有半个采样点的量化误差工程上常对峰值附近做插值如抛物线插值或 CZT来细化。4.2 MTI/MTD用多普勒滤波器组把地杂波滤掉雷达环境里最讨厌的不是噪声而是地物杂波。地面、建筑物、海面的回波强度往往比目标高出几个数量级而且多普勒频率集中在零速附近。对付它们有两级手段MTI动目标显示和 MTD动目标检测。MTI 用脉冲对消器最简单的两脉冲对消就是把相邻两个脉冲的回波相减。静止杂波在两次回波里几乎不变相减后被消除运动目标因为相位在脉冲间变化相减后仍有残余。三脉冲对消的陷波更宽对慢速杂波的抑制更深但也会压低慢速目标的响应。MTD 更进一步把同一距离单元上一连串脉冲的回波做 N 点 FFT等效于一组窄带多普勒滤波器每个滤波器对应一个速度通道。64 点 FFT、PRF 10kHz 时多普勒分辨率为 156.25Hz换算成 X 波段速度分辨约 2.34m/s。静止杂波落在 0 号通道目标落在哪个通道就大致知道速度。速度分辨越细需要的脉冲数越多但脉冲数增加意味着处理时间变长数据率下降。系统设计里这是个需要权衡的指标。MTD 滤波器组的另一个好处是相参积累一个目标如果能量集中在一个多普勒通道里FFT 等效于把 N 个脉冲的 SNR 叠加获得 N 倍的功率增益。128 点 FFT 能带来 21dB 增益这正是第 2 章算出的“单脉冲 SNR 0dB 也能检测”的底气。4.3 恒虚警检测 CFAR阈值为什么不写死经过脉冲压缩和 MTD 之后数据变成一张距离-多普勒图接下来要做的就是从图里找出“亮”的位置。最直观的方法是设一个固定阈值超过它就是目标。但噪声和杂波水平会随时间、环境变化固定阈值要么在噪声高时灌入大量虚警要么在噪声低时漏掉弱目标。CFAR 的思路是让阈值跟着背景走对待检测单元周围的参考单元求平均得到局部噪声估计再乘以一个系数作为阈值。最常用的是单元平均 CFARCA-CFAR阈值系数 α 与参考单元总数 N 和期望虚警率 Pfa 的关系是α N·(Pfa^(-1/N) − 1)这个公式的直观含义是参考窗越长噪声估计越稳α 可以取得越小Pfa 要求越严格α 必须越大。一个实用的 CA-CFAR 实现如下def ca_cfar_1d(amp, ref16, guard4, pfa1e-6): n len(amp) k_total 2 * ref alpha k_total * (pfa ** (-1.0 / k_total) - 1.0) threshold np.zeros(n) for i in range(n): left max(0, i - ref - guard) right min(n, i ref guard 1) win np.concatenate([amp[left:i - guard], amp[i guard 1:right]]) if len(win) 0: threshold[i] alpha * float(np.mean(win)) else: threshold[i] np.inf return threshold # 用法amp 为某一距离-多普勒切片的幅度 # det np.where(amp threshold)[0]这里 ref16 表示左右各取 16 个参考单元guard4 是保护单元。保护单元的作用是防止目标自身能量泄漏进参考窗把强目标的旁瓣误当成噪声抬高阈值。CA-CFAR 的局限在于参考窗内如果有另一个强目标阈值会被抬高弱目标被“遮蔽”。工程上的改良包括 GO-CFAR取左右窗较大值和 OS-CFAR对参考单元排序取第 k 个值前者保强目标后者抗多目标具体选型看场景。注意 CFAR 的输入应当是幅度而不是 dB 值。dB 是对数压缩会把弱小目标的对比度压扁影响检测效果。代码里 ref 和 pfa 都是可调项实际调测时先固定 pfa1e-6再根据虚警率实测值微调 ref 长度这是常规调参顺序。5. 连续波与 FMCW近距高精度测距测速的另一条路5.1 FMCW 原理差拍频率里同时藏着距离和速度脉冲雷达要发大功率短脉冲对器件要求高近距离场景里还有距离盲区。FMCW调频连续波换了个思路连续发射频率随时间线性变化的信号回波与发射信号混频后得到差拍频率这个频率里同时包含距离和速度信息。对于静止目标差拍频率 f_b 2·B·R/(c·Tc)B 是扫频带宽Tc 是扫频周期。目标运动时多普勒频移 f_d 2v/λ 会叠加到差拍频率上。锯齿波只能测到 f_b 和 f_d 的和或差没法分离三角波用上升沿和下降沿两次扫频加加减减就能解出来。三角波测速测距公式如下v (f_bu − f_bd)·λ/4R (f_bu f_bd)·c·Tc/(8B)。一个典型算例Tc1ms、B200MHz、载频 10GHz目标距离 100m 时距离对应的差拍是 133.3kHz目标速度 5m/s 时多普勒偏移约 333Hz。可见多普勒量级远小于距离项锯齿波当作纯测距用一般够但高精度测速或远距离时必须用三角波分离。5.2 一个可复现的 FMCW 距离-多普勒二维 FFT 仿真FMCW 的标准处理流程是对每个 chirp 做距离维 FFT得到“距离-慢时间”矩阵再沿慢时间维做 FFT得到距离-多普勒二维谱。下面这段代码生成两个目标的基带信号跑完整的二维 FFT 处理import numpy as np c 299_792_458.0 fc 10e9 lam c / fc B 200e6 Tc 1e-3 fs 4e6 # ADC 采样率 4 MHz Ns int(Tc * fs) # 单个 chirp 的采样点数 M 256 # chirp 数量 rng np.random.default_rng(0) t np.arange(Ns) / fs s np.zeros((M, Ns), dtypecomplex) # 两个目标(距离 m, 速度 m/s) targets [(50.0, 5.0), (120.0, -3.0)] for m in range(M): tm m * Tc for R, v in targets: fb 2 * B * R / (c * Tc) 2 * v / lam s[m] np.exp(2j * np.pi * fb * t 1j * 2 * np.pi * 2 * v / lam * tm) s 0.05 * (rng.normal(sizes.shape) 1j * rng.normal(sizes.shape)) win_r np.hanning(Ns) win_d np.hanning(M) S np.fft.fftshift(np.fft.fft(s * win_r, axis1), axes1) S np.fft.fftshift(np.fft.fft(S * win_d[:, None], axis0), axes0) S_db 20 * np.log10(np.abs(S) 1e-9) idx np.unravel_index(np.argmax(S_db), S_db.shape) fb_est (idx[1] - Ns / 2) / Ns * fs fd_est (idx[0] - M / 2) / M * (1 / Tc) R_est fb_est * c * Tc / (2 * B) v_est fd_est * lam / 2 print(f峰值: 距离 {R_est:.2f} m, 速度 {v_est:.2f} m/s)这段代码里第一个累加项 2BR/(cTc) 是距离产生的差拍第二个累加项 2v/λ 是多普勒。距离维 FFT 和速度维 FFT 各加了汉宁窗用来压低旁瓣。运行后峰值位置会落在 50m/5m/s 附近误差来自窗函数展宽和栅栏效应。想提高精度可以对二维峰值做 CZT 细化或者增大 chirp 数量 M。5.3 FMCW 的参数选择带宽、扫频时间、采样率怎么搭配FMCW 的参数设计有三个约束捋清楚就基本入门了。第一距离分辨率只由扫频带宽决定δR c/(2B)200MHz 给 0.75m1GHz 给 0.15m第二最大不模糊距离由 ADC 采样率和扫频时间共同决定差拍频率最高只能到 fs/2第三最大不模糊速度由 chirp 重复率 1/Tc 决定v_unamb λ/(4Tc)考虑正负范围。参数近距高精度远距离监视扫频带宽 B200MHz~1GHz50~100MHz扫频时间 Tc100~500μs1~5msADC 采样率 fs4~20MHz1~5MHzchirp 数量 M128~25664~128距离分辨率0.15~0.75m1.5~3mTc 拉长能让差拍频率分辨率变细距离测量更准但会压缩最大不模糊速度。所以车载雷达普遍把 Tc 控制在 100μs 量级宁可牺牲一点距离精度也要保住 30m/s 以上的测速范围。这个权衡与脉冲雷达里 PRF 的选择是同一个问题的两面。6. 雷达参数自检与体制选型的落地技巧6.1 用雷达方程做灵敏度闭环数字能对上才算开机雷达样机第一次上电先别急着看目标先验证接收链路是否工作在理论值附近。做法是取一个 RCS 已知的标准目标用雷达方程反推。工程上常用三面角反射器作为标准目标它的 RCS 理论值为 σ 4πa⁴/(3λ²)a 是角反射器的棱边长。这个公式值得记住外场实验时经常要现场算import math a 0.3 lam 0.03 sigma 4 * math.pi * a**4 / (3 * lam**2) print(fRCS {sigma:.2f} m^2 {10*math.log10(sigma):.1f} dBsm)a0.3m 的角反射器在 X 波段的理论 RCS 大约是 34m²15dBsm。把它放在已知距离处测量回波 SNR和雷达方程算出的理论值对比差值超过 3dB 基本能断定链路有问题。常见原因依次是天线增益与标称不符、接收机噪声系数劣化、发射功率实际值低于设定值。这个闭环验证方法比任何仿真都更能说明“雷达到底有没有调好”。6.2 三种体制和频段怎么选型雷达体制没有绝对优劣只看场景约束。脉冲多普勒雷达适合远距离、高低速目标同时存在的场景但需要处理距离和速度模糊FMCW 在近距高精度上有优势连续波体制峰值功率低、收发隔离是难点相控阵则是在波束灵活性上碾压前两者代价是成本和散热。维度脉冲多普勒FMCW相控阵作用距离远数十km中近数百m~数km远距离分辨率依赖带宽可达亚米级依赖带宽可达厘米级同左测速能力需解速度模糊天然测速范围受 Tc 限制与波形相关同时多目标强中强波束扫描机械/电子机械电扫成本与复杂度中低高选型给出三个自问作用距离是否超过 10km距离和速度分辨率是否要同时拉满成本与体积容不容得下多通道第一个问题回答“是”就直接排除掉 FMCW第二个问题指向脉冲多普勒加参差重频第三个问题是数字相控阵的入场券。6.3 外场验证的三个小技巧射频注入是外场快速验证接收链路的常用手段把发射信号经过 40~60dB 衰减器和一段已知延迟线直接馈入接收机替代天线收发。一段 10m 的电缆对应约 33ns 时延接收机测出来的距离应该严格等于这段时延换算值差多少就说明射频前端哪里偏了。这个方法能隔离天线问题快速定位链路失效环节。多径虚假目标在外场很常见特征是距离维出现间隔固定的双峰。若一个远离目标的回波峰和真实峰之间的距离恰好在某个角度上对应地面反射路径差基本就是镜面多径。判别的办法很直接把天线高度改变 10cm真实目标峰值移动很小多径假峰会移动两倍量因为反射路径对高度变化更敏感。最后一个常被忽略的检查是天线旁瓣。把波束指向偏开 3dB 波束宽度真实目标的回波 SNR 应该下降约 6dB发射和接收各 3dB。如果实测只降了 1~2dB说明旁瓣电平异常后续所有角度测量都会受到旁瓣干扰。用这个方法给天线做个体检比看方向图数据更贴近实际工作状态。本文还有配套的精品资源点击获取 简介雷达基础知识又名雷达原理是一份面向电子工程、通信工程及雷达相关专业初学者的幻灯片课件也适用于教师备课与自习巩固。课件以雷达定义和系统组成为切入点完整介绍发射机、接收机、天线、信号处理等核心模块并借助组成框图帮助理解整体工作流程。在测距方面详细讲解电磁波直线传播与反射测距原理给出回波延迟时间求距离的计算公式进一步推导最大不模糊距离、最小测量距离并通过实例演示脉冲重复频率与距离计算使抽象概念变得直观。在分辨力方面系统阐述距离分辨力和角度分辨力的定义、判定准则结合天线方向图说明方位角与俯仰角的获取方法同时配备战术雷达横向距离计算例题强化公式应用能力。资源共一个幻灯片文件大小约一点四兆字节幻灯片结构完整、公式与图例并重便于按章节系统学习雷达参数体系已有四千一百一十九人学习下载是入门雷达原理的高性价比参考资料。1. 雷达原理的第一课雷达在“看”什么把一辆静止的卡车和一辆移动中的自行车同时放在雷达面前回波强度差出几十倍但不少雷达算法最后跟踪的往往是自行车而不是卡车。这不是算法 bug而是雷达的工作原理决定的雷达并不像摄像头那样“看见”物体的轮廓它测量的是回波相对发射信号的延迟、频率和相位变化。延迟给出距离频率偏移给出径向速度不同接收通道之间的相位差给出角度。这三个物理量就是雷达原理里所有公式展开的原点。所以学雷达基础不能用“发射电磁波、接收回波”这九个字一笔带过。真正要建立的是从物理量到数据量的换算框架目标距离 1km对应电磁波往返约 6.7μs目标以 10m/s 径向运动在 10GHz 频段产生约 667Hz 的多普勒频移。这些数字是后面一切工程判断的标尺。这篇内容适合做自动驾驶感知、气象观测、安防监视、无人机避障以及正在转行做无线电和信号处理的工程师。雷达原理不只是一门电磁场课它本质上是“用有限的时间、带宽和功率从噪声里把想要的信息榨出来”的系统工程这恰恰是 IT 背景的人最容易上手的角度。2. 雷达测距原理与雷达方程功率怎么决定最远能看到哪2.1 测距的起点电磁波往返时间脉冲雷达的工作方式说起来很简单每隔一段固定时间发射一个短脉冲然后等待回波。如果脉冲发射后经过时间 t 收到回波目标距离就是 R c·t/2系数 1/2 是因为电磁波走了个来回。光速 c 取 299792458 m/s1μs 的时延对应 150m 距离。这里的核心结论是测距本质上是在测时间。时间测量的精度、时钟源的稳定性、接收机对回波到达时刻的判决方式直接决定距离精度。工程上通常不用回波前沿做判决因为前沿容易被噪声扰动而是用匹配滤波输出后的峰值位置。峰值判决更稳健但会引入至多半个采样间隔的固定偏置这在后文的脉冲压缩里还会再提到。2.2 雷达方程里的每个参数都在回答工程问题雷达方程是雷达原理里最不该跳过的一个公式它把发射、传播、目标、接收四个环节压缩成一行Pr Pt·G²·λ²·σ / ((4π)³·R⁴)逐项拆开看Pt 是峰值发射功率G 是天线增益收发共用一副天线所以出现 G²λ 是波长决定天线尺寸与大气损耗σ 是目标雷达截面积 RCS单位是平方米它不等于目标的物理面积R⁴ 是双向空间衰减距离翻一倍回波功率掉到原来的 1/16。这四个变量里R 是雷达无法控制的其余全是设计自由度。工程上更常用 dB 形式做预算。一个典型算例X 波段 10GHzλ0.03m峰值功率 10W天线增益 30dB目标 RCS 为 1m²距离 5km。套公式算下来接收功率约 -111dBm而 1MHz 带宽、3dB 噪声系数下的热噪声底约 -111dBm。也就是说单脉冲 SNR 接近 0dB刚好卡在检测边缘。这个数字说明一个工程事实雷达要看得远不能只靠提高峰值功率更要靠相干积累把多个脉冲的能量叠起来。下表是一组 X 波段脉冲雷达的常见参数后面章节的讨论都以它为基础参数符号典型值说明峰值功率Pt10W车载/ 1kW地面监视受器件与散热限制中心频率f010GHzX 波段λ3cm脉宽τ1μs决定基础距离分辨信号带宽B50~200MHz脉内调频后决定实际分辨天线增益G30dB收发共用脉冲重复频率PRF1~10kHz决定模糊边界噪声系数F3dB接收机前端质量积累脉冲数N64~128相干积累增益写个小脚本算一下 SNR 随距离的变化比手算快得多import numpy as np c 299_792_458.0 f0 10e9 lam c / f0 Pt 10.0 # 峰值功率单位 W G 10 ** (30 / 10) # 天线增益 30 dB sigma 1.0 # RCS单位 m^2 kB 1.380649e-23 T0 290.0 # 参考噪声温度单位 K B 1e6 # 匹配滤波带宽 1 MHz NF 10 ** (3 / 10) # 噪声系数 3 dB R np.linspace(1e3, 20e3, 200) Pr Pt * G**2 * lam**2 * sigma / ((4 * np.pi)**3 * R**4) Pn kB * T0 * B * NF SNR 10 * np.log10(Pr / Pn) for r in [3000, 5000, 10000]: pr Pt * G**2 * lam**2 * sigma / ((4 * np.pi)**3 * r**4) print(fR {r/1000:.0f} km, SNR {10*np.log10(pr / Pn):.1f} dB)这段代码把雷达方程拆成线性计算再用 log 转成 SNR。运行结果里 5km 处 SNR 在 0dB 附近与手算一致。脚本里的可调参数是 Pt、G、NF 和 B把 NF 从 3dB 降到 1.5dBSNR 能改善 1.5dB想靠提高功率补回 3dB得把 Pt 翻倍。所以接收机噪声系数往往是雷达链路里最值得花钱的地方。2.3 距离分辨率关键不在脉宽在带宽测距精度和距离分辨率是两个概念但常被混淆。距离分辨率指雷达能把两个相距多近的目标区分开的能力。对单载频矩形脉冲分辨率约等于 δR c·τ/21μs 脉宽对应 150m。也就是说两辆车相距 150m 以内在普通脉冲雷达上看就是一个目标。但现代雷达几乎都采用脉内线性调频LFM发射脉冲被展宽到 1μs但内部频率在 200MHz 带宽内扫描。接收端做匹配滤波后等效脉冲宽度被压缩到约 1/B距离分辨率变成 δR c/(2B)200MHz 带宽对应 0.75m。这就是雷达原理里最常被新手绕晕的一点分辨率由带宽决定不是由脉宽决定。压窄脉宽能提高分辨率但峰值功率要求急剧上升展宽带宽再压缩成本在信号处理链路而不是发射机。现代雷达系统几乎都选择后者这也是下一章要展开的波形设计问题。3. 脉冲雷达的波形与模糊重频、盲速和距离模糊怎么选3.1 脉冲重复频率决定了能看多远、测多快脉冲雷达每发射一个脉冲要在下一个脉冲发射前收到上一脉冲的回波否则距离就模糊了。因此最大不模糊距离是 R_unamb c/(2·PRF)。与此同时对回波做多普勒处理时脉冲采样定理限制了不模糊速度范围最大不模糊速度是 v_unamb λ·PRF/2实际可测范围是 ±v_unamb/2。这两个公式放在一起就是一个经典矛盾PRF 高看得近但测速范围大PRF 低看得远但速度容易模糊。以第 2 章那组 X 波段参数为例PRF最大不模糊距离最大不模糊速度±典型场景1kHz150km±7.5m/s远程警戒5kHz30km±37.5m/s中程监视20kHz7.5km±150m/s高速目标跟踪注意 1kHz 那一行±7.5m/s 的测速范围连正常行驶的汽车都覆盖不了。所以远程雷达不是单纯选低 PRF 就完事速度模糊必须配合解模糊算法一起处理。3.2 重频参差与双重频解模糊的常见工程做法既然低 PRF 测距清晰、高 PRF 测速清晰常见的工程做法就是“两个都不够就两个都发”。典型方案是重频参差交替使用两个相近但不相同的 PRF各自得到带模糊的测量值再通过余数定理思想在候选距离里找一致解。下面给出一个双重频解距离模糊的简化实现def resolve_range(r1, r2, prf1, prf2, max_r): c 299_792_458.0 ramb1 c / (2 * prf1) ramb2 c / (2 * prf2) tol min(ramb1, ramb2) * 0.25 candidates [] for n1 in range(int(max_r / ramb1) 1): d1 r1 n1 * ramb1 for n2 in range(int(max_r / ramb2) 1): d2 r2 n2 * ramb2 if abs(d1 - d2) tol: candidates.append((d1 d2) / 2) return candidates # 例真实距离 82.3kmPRF11kHz(不模糊150km)PRF21.1kHz # 两个雷达各自量到的表观距离需要由前端给出 print(resolve_range(82.3e3, 82.3e3, 1000, 1100, 300e3))这里的核心思想是遍历模糊重数 n1、n2把两个带模糊的距离各自外推找出互相吻合的候选值。tol 的取值要大于距离测量误差但又不能大到把多个候选都吞进来。实际系统里往往加入第三个 PRF 来消除余数定理的二次模糊三重频的正确解概率会显著高于双重频。3.3 距离-多普勒耦合的边界模糊函数不是考试概念LFM 信号有一个绕不开的性质距离-多普勒耦合。目标的多普勒频移 f_d 会让压缩后的峰值在时延轴上偏移一个量偏移量约等于 f_d 除以调频斜率 k单位 Hz/s。目标速度越快测距偏差越大。对 200MHz/1ms 的扫频参数10m/s 目标大约造成 1m 量级的距离偏移这在近距高精度场景里已经不能忽略。模糊函数描述的是“波形对距离和速度联合分辨的能力边界”它回答的问题是给定一个波形在距离-速度二维平面上能分辨到多细。工程上选波形本质上就是选模糊函数的形状LFM 是倾斜的刀刃形相位编码是图钉形各自适合不同的场景。理解到这一层再看后面 MTD 和 FMCW 的设计就能明白每个参数为什么出现在那个位置。4. 从回波到点迹脉冲压缩、MTD 与 CFAR 的实现链路4.1 匹配滤波与脉冲压缩把 1μs 变成 1.5m回波进接收机后第一个关键步骤是匹配滤波。匹配滤波的理论结论是在白色高斯噪声背景下使输出 SNR 最大的线性滤波器是发射信号的共轭时间翻转。对 LFM 信号做匹配滤波输出波形被压缩成宽度约 1/B 的窄峰同时获得 B·τ 倍的 SNR 增益。1μs 脉宽、200MHz 带宽的 LFM压缩比是 200即 23dB 的积累增益距离分辨率从 150m 提升到 0.75m。这一步实现起来并不神秘频域做一次复数相乘加 IFFT 即可import numpy as np fs 200e6 T 1e-6 B 200e6 N int(fs * T) t np.arange(N) / fs k B / T lfm_tx np.exp(1j * 2 * np.pi * (k / 2) * t**2) # 发射信号 delay 100 # 回波延迟 100 个采样点 lfm_rx np.zeros(N, dtypecomplex) lfm_rx[delay:] lfm_tx[:N - delay] # 模拟回波 H np.conj(np.fft.fft(lfm_tx)) # 匹配滤波器 pulse_comp np.fft.ifft(np.fft.fft(lfm_rx) * H) # 频域滤波 peak np.argmax(np.abs(pulse_comp)) print(f峰值位置: {peak}, 理论延迟: {delay}, 分辨率: {c / (2*B):.2f} m)注意匹配滤波器用的是发射信号的共轭谱回波经过它之后输出峰值的位置就对应目标时延。这里的关键参数是采样率 fs 和带宽 B峰值定位精度受采样率限制如果 fs 只取 200MHz逐点峰值会有半个采样点的量化误差工程上常对峰值附近做插值如抛物线插值或 CZT来细化。4.2 MTI/MTD用多普勒滤波器组把地杂波滤掉雷达环境里最讨厌的不是噪声而是地物杂波。地面、建筑物、海面的回波强度往往比目标高出几个数量级而且多普勒频率集中在零速附近。对付它们有两级手段MTI动目标显示和 MTD动目标检测。MTI 用脉冲对消器最简单的两脉冲对消就是把相邻两个脉冲的回波相减。静止杂波在两次回波里几乎不变相减后被消除运动目标因为相位在脉冲间变化相减后仍有残余。三脉冲对消的陷波更宽对慢速杂波的抑制更深但也会压低慢速目标的响应。MTD 更进一步把同一距离单元上一连串脉冲的回波做 N 点 FFT等效于一组窄带多普勒滤波器每个滤波器对应一个速度通道。64 点 FFT、PRF 10kHz 时多普勒分辨率为 156.25Hz换算成 X 波段速度分辨约 2.34m/s。静止杂波落在 0 号通道目标落在哪个通道就大致知道速度。速度分辨越细需要的脉冲数越多但脉冲数增加意味着处理时间变长数据率下降。系统设计里这是个需要权衡的指标。MTD 滤波器组的另一个好处是相参积累一个目标如果能量集中在一个多普勒通道里FFT 等效于把 N 个脉冲的 SNR 叠加获得 N 倍的功率增益。128 点 FFT 能带来 21dB 增益这正是第 2 章算出的“单脉冲 SNR 0dB 也能检测”的底气。4.3 恒虚警检测 CFAR阈值为什么不写死经过脉冲压缩和 MTD 之后数据变成一张距离-多普勒图接下来要做的就是从图里找出“亮”的位置。最直观的方法是设一个固定阈值超过它就是目标。但噪声和杂波水平会随时间、环境变化固定阈值要么在噪声高时灌入大量虚警要么在噪声低时漏掉弱目标。CFAR 的思路是让阈值跟着背景走对待检测单元周围的参考单元求平均得到局部噪声估计再乘以一个系数作为阈值。最常用的是单元平均 CFARCA-CFAR阈值系数 α 与参考单元总数 N 和期望虚警率 Pfa 的关系是α N·(Pfa^(-1/N) − 1)这个公式的直观含义是参考窗越长噪声估计越稳α 可以取得越小Pfa 要求越严格α 必须越大。一个实用的 CA-CFAR 实现如下def ca_cfar_1d(amp, ref16, guard4, pfa1e-6): n len(amp) k_total 2 * ref alpha k_total * (pfa ** (-1.0 / k_total) - 1.0) threshold np.zeros(n) for i in range(n): left max(0, i - ref - guard) right min(n, i ref guard 1) win np.concatenate([amp[left:i - guard], amp[i guard 1:right]]) if len(win) 0: threshold[i] alpha * float(np.mean(win)) else: threshold[i] np.inf return threshold # 用法amp 为某一距离-多普勒切片的幅度 # det np.where(amp threshold)[0]这里 ref16 表示左右各取 16 个参考单元guard4 是保护单元。保护单元的作用是防止目标自身能量泄漏进参考窗把强目标的旁瓣误当成噪声抬高阈值。CA-CFAR 的局限在于参考窗内如果有另一个强目标阈值会被抬高弱目标被“遮蔽”。工程上的改良包括 GO-CFAR取左右窗较大值和 OS-CFAR对参考单元排序取第 k 个值前者保强目标后者抗多目标具体选型看场景。注意 CFAR 的输入应当是幅度而不是 dB 值。dB 是对数压缩会把弱小目标的对比度压扁影响检测效果。代码里 ref 和 pfa 都是可调项实际调测时先固定 pfa1e-6再根据虚警率实测值微调 ref 长度这是常规调参顺序。5. 连续波与 FMCW近距高精度测距测速的另一条路5.1 FMCW 原理差拍频率里同时藏着距离和速度脉冲雷达要发大功率短脉冲对器件要求高近距离场景里还有距离盲区。FMCW调频连续波换了个思路连续发射频率随时间线性变化的信号回波与发射信号混频后得到差拍频率这个频率里同时包含距离和速度信息。对于静止目标差拍频率 f_b 2·B·R/(c·Tc)B 是扫频带宽Tc 是扫频周期。目标运动时多普勒频移 f_d 2v/λ 会叠加到差拍频率上。锯齿波只能测到 f_b 和 f_d 的和或差没法分离三角波用上升沿和下降沿两次扫频加加减减就能解出来。三角波测速测距公式如下v (f_bu − f_bd)·λ/4R (f_bu f_bd)·c·Tc/(8B)。一个典型算例Tc1ms、B200MHz、载频 10GHz目标距离 100m 时距离对应的差拍是 133.3kHz目标速度 5m/s 时多普勒偏移约 333Hz。可见多普勒量级远小于距离项锯齿波当作纯测距用一般够但高精度测速或远距离时必须用三角波分离。5.2 一个可复现的 FMCW 距离-多普勒二维 FFT 仿真FMCW 的标准处理流程是对每个 chirp 做距离维 FFT得到“距离-慢时间”矩阵再沿慢时间维做 FFT得到距离-多普勒二维谱。下面这段代码生成两个目标的基带信号跑完整的二维 FFT 处理import numpy as np c 299_792_458.0 fc 10e9 lam c / fc B 200e6 Tc 1e-3 fs 4e6 # ADC 采样率 4 MHz Ns int(Tc * fs) # 单个 chirp 的采样点数 M 256 # chirp 数量 rng np.random.default_rng(0) t np.arange(Ns) / fs s np.zeros((M, Ns), dtypecomplex) # 两个目标(距离 m, 速度 m/s) targets [(50.0, 5.0), (120.0, -3.0)] for m in range(M): tm m * Tc for R, v in targets: fb 2 * B * R / (c * Tc) 2 * v / lam s[m] np.exp(2j * np.pi * fb * t 1j * 2 * np.pi * 2 * v / lam * tm) s 0.05 * (rng.normal(sizes.shape) 1j * rng.normal(sizes.shape)) win_r np.hanning(Ns) win_d np.hanning(M) S np.fft.fftshift(np.fft.fft(s * win_r, axis1), axes1) S np.fft.fftshift(np.fft.fft(S * win_d[:, None], axis0), axes0) S_db 20 * np.log10(np.abs(S) 1e-9) idx np.unravel_index(np.argmax(S_db), S_db.shape) fb_est (idx[1] - Ns / 2) / Ns * fs fd_est (idx[0] - M / 2) / M * (1 / Tc) R_est fb_est * c * Tc / (2 * B) v_est fd_est * lam / 2 print(f峰值: 距离 {R_est:.2f} m, 速度 {v_est:.2f} m/s)这段代码里第一个累加项 2BR/(cTc) 是距离产生的差拍第二个累加项 2v/λ 是多普勒。距离维 FFT 和速度维 FFT 各加了汉宁窗用来压低旁瓣。运行后峰值位置会落在 50m/5m/s 附近误差来自窗函数展宽和栅栏效应。想提高精度可以对二维峰值做 CZT 细化或者增大 chirp 数量 M。5.3 FMCW 的参数选择带宽、扫频时间、采样率怎么搭配FMCW 的参数设计有三个约束捋清楚就基本入门了。第一距离分辨率只由扫频带宽决定δR c/(2B)200MHz 给 0.75m1GHz 给 0.15m第二最大不模糊距离由 ADC 采样率和扫频时间共同决定差拍频率最高只能到 fs/2第三最大不模糊速度由 chirp 重复率 1/Tc 决定v_unamb λ/(4Tc)考虑正负范围。参数近距高精度远距离监视扫频带宽 B200MHz~1GHz50~100MHz扫频时间 Tc100~500μs1~5msADC 采样率 fs4~20MHz1~5MHzchirp 数量 M128~25664~128距离分辨率0.15~0.75m1.5~3mTc 拉长能让差拍频率分辨率变细距离测量更准但会压缩最大不模糊速度。所以车载雷达普遍把 Tc 控制在 100μs 量级宁可牺牲一点距离精度也要保住 30m/s 以上的测速范围。这个权衡与脉冲雷达里 PRF 的选择是同一个问题的两面。6. 雷达参数自检与体制选型的落地技巧6.1 用雷达方程做灵敏度闭环数字能对上才算开机雷达样机第一次上电先别急着看目标先验证接收链路是否工作在理论值附近。做法是取一个 RCS 已知的标准目标用雷达方程反推。工程上常用三面角反射器作为标准目标它的 RCS 理论值为 σ 4πa⁴/(3λ²)a 是角反射器的棱边长。这个公式值得记住外场实验时经常要现场算import math a 0.3 lam 0.03 sigma 4 * math.pi * a**4 / (3 * lam**2) print(fRCS {sigma:.2f} m^2 {10*math.log10(sigma):.1f} dBsm)a0.3m 的角反射器在 X 波段的理论 RCS 大约是 34m²15dBsm。把它放在已知距离处测量回波 SNR和雷达方程算出的理论值对比差值超过 3dB 基本能断定链路有问题。常见原因依次是天线增益与标称不符、接收机噪声系数劣化、发射功率实际值低于设定值。这个闭环验证方法比任何仿真都更能说明“雷达到底有没有调好”。6.2 三种体制和频段怎么选型雷达体制没有绝对优劣只看场景约束。脉冲多普勒雷达适合远距离、高低速目标同时存在的场景但需要处理距离和速度模糊FMCW 在近距高精度上有优势连续波体制峰值功率低、收发隔离是难点相控阵则是在波束灵活性上碾压前两者代价是成本和散热。维度脉冲多普勒FMCW相控阵作用距离远数十km中近数百m~数km远距离分辨率依赖带宽可达亚米级依赖带宽可达厘米级同左测速能力需解速度模糊天然测速范围受 Tc 限制与波形相关同时多目标强中强波束扫描机械/电子机械电扫成本与复杂度中低高选型给出三个自问作用距离是否超过 10km距离和速度分辨率是否要同时拉满成本与体积容不容得下多通道第一个问题回答“是”就直接排除掉 FMCW第二个问题指向脉冲多普勒加参差重频第三个问题是数字相控阵的入场券。6.3 外场验证的三个小技巧射频注入是外场快速验证接收链路的常用手段把发射信号经过 40~60dB 衰减器和一段已知延迟线直接馈入接收机替代天线收发。一段 10m 的电缆对应约 33ns 时延接收机测出来的距离应该严格等于这段时延换算值差多少就说明射频前端哪里偏了。这个方法能隔离天线问题快速定位链路失效环节。多径虚假目标在外场很常见特征是距离维出现间隔固定的双峰。若一个远离目标的回波峰和真实峰之间的距离恰好在某个角度上对应地面反射路径差基本就是镜面多径。判别的办法很直接把天线高度改变 10cm真实目标峰值移动很小多径假峰会移动两倍量因为反射路径对高度变化更敏感。最后一个常被忽略的检查是天线旁瓣。把波束指向偏开 3dB 波束宽度真实目标的回波 SNR 应该下降约 6dB发射和接收各 3dB。如果实测只降了 1~2dB说明旁瓣电平异常后续所有角度测量都会受到旁瓣干扰。用这个方法给天线做个体检比看方向图数据更贴近实际工作状态。本文还有配套的精品资源点击获取