简介这是一份基于MATLAB的雷达系统建模与仿真资料面向雷达信号处理初学者与研究者覆盖目标回波模拟、高频放大、混频、中频放大、相干检波、模数转换、脉冲压缩、MTI、积累和CFAR检测等完整信号处理链。压缩包共64个文件包含36个M脚本每个模块均配有可运行程序、对应ASV备份、11张JPG输出图形、两份DOC报告及FIG、LOG、DB等辅助文件整体约780KB结构清晰便于按模块对照学习。已有418人学习下载。借助这些脚本与图形读者可快速复现雷达收发链路的仿真结果理解各环节参数设置与算法原理也可在此基础上修改参数、扩展功能用于课程设计或科研验证。1. 为什么把雷达接收机拆成十个 MATLAB 模块拿到这套雷达系统建模与仿真的源码包时我先把几十个 .m 文件按接收链路顺序捋了一遍发现它覆盖的正是教材里雷达接收机的完整框图huibo.m 生成目标回波gaofang.m 做高频放大hunpin.m 混频zhongfang.m 中频放大xiangganjianbo.m 相干检波AD.m 模数转换maichongyasuo.m 脉冲压缩MTI.m 动目标指示jilei.m 多脉冲积累CFAR.m 恒虚警检测。这套东西适合两类人一是刚开始接触雷达信号处理需要把高放、混频、中放这些概念和代码一一对上的学生二是要做系统参数论证或课程设计的工程师想直接拿可运行的模板改参数。整个工程用 MATLAB 自带的信号处理函数就能跑通不依赖额外硬件每个模块输出的图形也都齐。下文按信号流逐级展开每级给出代码、参数依据和调参时容易踩的坑。2. 前端链路建模目标回波、增益分配与混频中放2.1 目标回波里该放哪些物理量huibo.m 要解决的第一件事是把目标的距离 R、径向速度 v、雷达截面积 σ 写进回波模型。若发射信号是线性调频脉冲理想的接收回波相对发射信号有三个变化延迟 τ 2R/c多普勒频移 fd 2v/λ幅度衰减正比于 R 的四次方σ 体现在回波功率上。很多初学者只做延迟不加多普勒结果后面的 MTI 和积累模块什么都看不到。包里还附带了 swerling.m这是 Swerling 起伏模型用来描述 σ 在帧间或脉冲间的随机性做蒙特卡洛统计时很有用定频仿真时可以先绕开。我建议回波生成只保留三个自由度延迟、多普勒、幅度。RCS 起伏在后续多帧统计时再叠加不要在一个脉冲里混入两种随机性否则出了问题很难定位。下面是简化后的回波生成核心代码% 参数声明 c 3e8; % 光速 m/s f0 10e9; % 载频 10 GHz lambda c / f0; % 波长 B 10e6; % 发射带宽 10 MHz T 50e-6; % 脉冲宽度 50 us fs 40e6; % 采样率 t 0 : 1/fs : T - 1/fs; s_tx exp(1j * pi * (B/T) * t.^2); % LFM 基带波形 R0 5000; v 30; % 目标距离 5 km径向速度 30 m/s tau 2 * R0 / c; % 回波延迟 fd 2 * v / lambda; % 多普勒频率 s_rx s_tx .* exp(-1j * 2 * pi * f0 * tau) ... .* exp(1j * 2 * pi * fd * t); % 延迟多普勒逻辑说明exp(-1j * 2 * pi * f0 * tau)是回波相对发射信号的固定相移它由目标距离决定这个相位在脉冲压缩前看不出作用但到 MTI 和相参积累时会直接影响相邻脉冲间的相位关系exp(1j * 2 * pi * fd * t)是目标运动引起的相位旋转这正是 MTI 和积累能区分动目标和静物回波的物理基础。发射信号用基带复包络表示而不是真的把 10 GHz 载频扔进采样是因为采样率不用跟着载频走这是雷达仿真里最常用的低通等效建模手法。提示如果严格按雷达方程算回波幅度需要给出发射功率、天线增益、噪声系数等一整套参数。仿真包里多数用归一化幅度再加高斯白噪声先把信噪比做成可调变量这对调试后端算法比追求绝对功率值更有价值。2.2 高放、混频、中放的增益与噪声分配高放 gaofang.m、混频 hunpin.m、中放 zhongfang.m 对应的是接收机的模拟前端。高放把天线收到的微弱回波放大到混频器能正常工作的电平混频把射频信号与本地振荡器信号相乘完成频率搬移中放再对中频信号做一次窄带放大和滤波。三级增益分配有一条基本约束高放增益不能太大否则强目标回波会把混频器推到饱和区产生交调失真但高放噪声系数又直接决定接收机整机灵敏度。级联噪声系数公式 F F1 (F2-1)/G1 (F3-1)/(G1·G2) 在这套仿真里的落地方式是每级放大后往信号里叠一份高斯白噪声噪声功率按该级带宽和增益系数折算而不是只做无噪声的纯放大。下面是一组常见的参数分配参考模块频率范围典型增益噪声系数带宽高放射频 9.4 GHz20 dB2 dB100 MHz混频射频转中频-6 dB8 dB100 MHz中放中频 30 MHz30 dB6 dB10 MHz相干检波基带 I/Q0 dB—10 MHz混频在数学上就是乘法s_if s_rf .* cos(2*pi*f_lo*t)MATLAB 里用点乘直接实现。注意混频输出不止有中频分量还叠加了和频分量实际仿真里要么在混频后加低通或带通滤波器要么直接用复信号让负频率分量不落在采样带宽内。如果看到 hunpin.m 的输出频谱出现两个对称峰值多半是漏了滤波这一级。中放的带宽要和发射信号带宽匹配带宽太窄会把 LFM 信号边缘削掉导致后续脉冲压缩的输出主瓣变宽。2.3 相干检波为什么必须保留 I/Q 两路中放输出到相干检波 xiangganjianbo.m这里最容易犯的错误是把检波理解成取包络。相干检波的工作原理是将中频信号分别与同相参考信号 cos(2πf_if·t) 和正交参考信号 -sin(2πf_if·t) 相乘经低通滤波后得到 I/Q 两路基带信号。这么做的意义有两个一是保留相位信息后续脉冲压缩的匹配滤波、MTI 的相位对消都依赖它二是信号变成复基带后采样率只需要覆盖信号带宽不需要扛着 30 MHz 的中频频率。% 相干检波核心I/Q 两路混频加低通 t 0 : 1/fs : N/fs - 1/fs; % N 为样本数 f_if 30e6; % 中频频率 x_i 2 * s_if .* cos(2 * pi * f_if * t); % 同相支路 x_q -2 * s_if .* sin(2 * pi * f_if * t); % 正交支路 % 设计低通滤波器通带 5 MHz 对应带宽一半 lpFilt designfilt(lowpassfir, PassbandFrequency, 5e6, ... StopbandFrequency, 8e6, SampleRate, fs); z filter(lpFilt, x_i 1j * x_q); % 合成复基带信号逻辑说明乘以 2 是工程惯例乘法器输出中差频项的系数只有 1/2乘 2 可以把幅度归一。低通滤波器去掉的是和频分量。PassbandFrequency 取 5 MHz对应 10 MHz 信号带宽的一半如果这里取得太大和频分量会残留基带信号会被污染。整条模拟前端调参的正确顺序是先固定噪声系数调整各级增益观察输出信噪比变化再回头看混频器和中放的带宽对脉冲沿的影响。信噪比异常时优先怀疑混频后的滤波而不是增益不够。3. ADC 量化与脉冲压缩距离分辨率的关键3.1 ADC 的采样率和量化位数怎么定AD.m 处理的是相干检波输出的模拟基带信号。两个参数决定这个环节的信息损失采样率和量化位数。基带信号带宽若为 B采样率取 fs ≥ B 就能满足复采样条件工程上通常会留 20% 裕量带宽 10 MHz 的信号采样率取 40 MHz 已经富余。量化位数决定信号动态范围仿真里默认取 12 位以上量化噪声远小于接收机噪声时可以忽略不需要刻意调低去测试量化效应。% 均匀量化模型 n_bits 12; % 量化位数 Vmax max(abs(z)); % 信号峰值 L 2^n_bits; % 量化电平数 delta 2 * Vmax / L; % 量化步长 z_q round(z / delta) * delta; % 量化后的复信号 % 量化信噪比理论值约 6.02*n_bits 1.76 dB sqnr_db 6.02 * n_bits 1.76;逻辑说明round(z / delta) * delta是均匀量化器的标准实现把连续幅值映射到最近的电平上。量化信噪比的理论公式 6.02n 1.76 dB 在这里可以作为参考12 位量化对应约 74 dB 的动态范围远超雷达接收链路中 -10 到 20 dB 的典型信噪比区间因此 ADC 不是系统瓶颈。如果发现量化后的仿真结果和理想情况差异明显先检查 Vmax 取值是否合理峰值取得过大或过小都会浪费量化区间。提示ADC 真实器件还存在采样抖动、孔径延迟和非线性这套仿真没有细化到器件级。做系统级论证时ADC 模型保持理想把注意力放在脉冲压缩、MTI、CFAR 这些算法链路上才是对的。3.2 脉冲压缩匹配滤波与 LFM 信号脉冲压缩是 maichongyasuo.m 的核心。宽脉冲保证足够的平均发射功率调制带宽决定距离分辨率接收端用匹配滤波把回波压缩成窄脉冲。LFM 信号的匹配滤波输出是一个 sinc 型脉冲主瓣宽度约 1/B因此距离分辨率 ΔR c/(2B)。这里有一个常见误解分辨率只由带宽决定和脉冲宽度无关带宽 10 MHz 对应分辨率约 15 m。匹配滤波有两种实现方式时域卷积和频域相乘。数据量大时频域实现更快因为一次 FFT 的复杂度是 O(N log N)而时域卷积是 O(N²)。代码如下% 脉冲压缩频域匹配滤波 Nfft 2^nextpow2(length(s_rx) length(s_tx) - 1); S_rx fft(s_rx, Nfft); S_ref fft(conj(fliplr(s_tx)), Nfft); % 匹配滤波器的频响 y ifft(S_rx .* S_ref, Nfft); y y(1 : length(s_rx) length(s_tx) - 1); % 归一化后取幅度 y_abs abs(y) / max(abs(y));逻辑说明匹配滤波器的冲激响应是发射信号的共轭时间反演对应代码里的conj(fliplr(s_tx))。频域相乘相当于时域卷积nextpow2保证 FFT 长度覆盖线性卷积输出长度避免循环卷积的混叠。输出主瓣宽度约 1/B若在频域给参考信号加窗比如汉明窗可以将距离旁瓣压低到 -40 dB 左右代价是主瓣展宽约 1.5 倍这个权衡在弱目标检测时经常要用到。表格里列出不同窗函数的典型旁瓣和主瓣展宽系数可以当作调参速查窗函数最高旁瓣电平主瓣展宽系数矩形窗不加窗-13.3 dB1.0汉宁窗-31.5 dB1.62汉明窗-42.7 dB1.47布莱克曼窗-58.1 dB1.73用窗的时机也很关键只在匹配滤波的频域参考信号上乘窗不要在回波域直接加窗否则会损失信噪比。4. MTI 对消与脉冲积累滤除静物回波4.1 一次对消与二次对消的传递函数分析MTI.m 实现的是动目标指示。它的思路建立在这样一个事实上静止地物回波在相邻脉冲间的相位保持不变而运动目标的相位按多普勒频率旋转。因此将相邻两个脉冲的回波相减静物回波被抵消运动目标保留下来这就是一次对消。在 z 域一次对消的传递函数是 H(z) 1 - z⁻¹频率响应在零多普勒处有一个零点。代码实现非常简洁% 一次对消相邻脉冲相减 mti_out z_pulse(2:end, :) - z_pulse(1:end-1, :); % 二次对消再对一次对消结果做一次差分 mti_out2 mti_out(2:end, :) - mti_out(1:end-1, :); % 频响对比多普勒频率轴 f_dop linspace(-PRF/2, PRF/2, 256); H1 abs(1 - exp(-1j * 2 * pi * f_dop / PRF)); H2 abs(1 - exp(-1j * 2 * pi * f_dop / PRF)).^2;逻辑说明z_pulse是脉冲维 × 距离维的二维矩阵z_pulse(2:end,:) - z_pulse(1:end-1,:)把相邻脉冲逐距离门相减。第一行代码做一次对消第二行做二次对消。频响表达式里 PRF 是脉冲重复频率零多普勒处的陷波越深对静止杂波的抑制越好。二次对消的陷波更宽适合处理有一定谱宽的地物杂波但代价是低速目标的响应也会被削掉一部分最小可检测速度会抬升。4.2 相参积累与非相参积累的取舍积累在 jilei.m 里实现作用是把多个脉冲的回波能量叠加提高信噪比。积累分两种相参积累在检波之前对复信号做非相参积累在检波之后对幅度做。相参积累理论上有 M 倍的 SNR 增益非相参积累只有大约 M 的 0.8 次方但相参积累要求目标在积累时间内不发生跨距离门走动并且多普勒频率必须落在滤波器通带内。常见做法是用 FFT 做相参积累把 M 个脉冲在距离维的复数据沿慢时间做 FFT得到距离-多普勒图% 相参积累沿慢时间维做 FFT M 64; % 积累脉冲数 rd_map fftshift(fft(z_matrix, M, 1), 1); % z_matrix: 脉冲×距离 % 频率轴换算 f_axis (-M/2 : M/2-1) * PRF / M; % 多普勒分辨率与积累时间的关系 df_dop PRF / M; v_res df_dop * lambda / 2; % 速度分辨率逻辑说明fft(..., M, 1)沿第一维脉冲维做 M 点 FFTfftshift把零多普勒移到频谱中心便于观察。多普勒分辨率等于 PRF/M速度分辨率再乘 λ/2。比如 PRF 为 1 kHz64 点积累多普勒分辨率约 15.6 HzX 波段速度分辨率约 0.23 m/s。积累脉冲数 M 不能贪多目标在积累时间内如果发生跨距离门走动能量会散到多个距离门里增益反而下降。遇到高速目标时要用包络走动补偿或者改用更短的积累时间。MTI 和积累的配合关系是先做 MTI 对消掉强静物杂波再做积累提高目标信噪比。顺序不能反如果先积累再做 MTI强杂波会把积累后的动态范围撑爆弱小目标就检测不出来了。5. CFAR 检测门限的工程设定与实测验证5.1 CA-CFAR 的参数依据CFAR.m 实现的是单元平均恒虚警检测。基本原理是检测单元两侧各取 N 个参考单元求平均作为杂波功率估计再乘以门限因子 α 得到检测门限。参考单元数和保护单元数是最关键的两个参数。保护单元紧邻检测单元用来挡住目标回波主瓣和旁瓣泄漏进参考单元导致门限被抬高。对脉压后的 LFM 信号主瓣大约占 2~3 个距离单元保护单元每侧取 2~4 个即可。门限因子 α 由虚警概率 Pfa 和参考单元数 N 决定关系式是 α N(Pfa^(-1/N) - 1)。它的推导基于参考单元噪声服从指数分布这一假设实际仿真里如果噪声是高斯分布经平方率检波后也近似成立。下面给出一段可运行的 CFAR 检测核心代码function det cfar_1d(x, N_ref, N_guard, Pfa) % x: 输入距离维数据, N_ref: 参考单元数(每侧), N_guard: 保护单元数(每侧) alpha N_ref * (Pfa^(-1/N_ref) - 1); det zeros(size(x)); len length(x); for k N_guard N_ref 1 : len - N_guard - N_ref ref_left x(k - N_guard - N_ref : k - N_guard - 1); ref_right x(k N_guard 1 : k N_guard N_ref); noise_est (sum(ref_left) sum(ref_right)) / (2 * N_ref); threshold alpha * noise_est; det(k) x(k) threshold; end end逻辑说明循环从头到尾滑动对每个距离单元 k取左侧和右侧各 N_ref 个参考单元求平均得到噪声功率估计noise_est门限是 α 乘以噪声估计。det(k) x(k) threshold输出二值判断结果。Pfa 取 1e-6 到 1e-4 之间是雷达仿真的常见设置Pfa 定得过低会大幅抬高门限导致检测概率下降过高则会出现大量虚警。参考单元数 N_ref 取 16 到 32 之间较稳妥太小则噪声估计方差大门限抖动明显太大则在非均匀杂波环境中响应迟缓。5.2 用蒙特卡洛验证虚警率与检测概率评价 CFAR 参数设得合不合理最直接的办法是跑蒙特卡洛仿真。对纯噪声输入统计虚警数对含目标输入统计检测概率两条曲线能直观看出门限因子和参考窗长度的权衡。% 蒙特卡洛验证虚警率统计 Nmc 1000; N_cells 1000; cnt 0; for i 1:Nmc noise abs(randn(1, N_cells) 1j*randn(1, N_cells)).^2 / 2; det cfar_1d(noise, 16, 4, 1e-4); cnt cnt sum(det); end pfa_measured cnt / (Nmc * N_cells); fprintf(理论 Pfa: 1e-4, 实测 Pfa: %e\n, pfa_measured);逻辑说明randn生成复高斯噪声取模平方除以 2 得到指数分布噪声功率样本这是平方率检波后的标准噪声模型。实测虚警率和理论值通常在一个数量级内吻合如果偏差过大优先检查保护单元是否过少或参考单元内的强杂波没有被清除。检测概率的仿真则是把已知信噪比的信号插进噪声统计检测到目标的次数占比通常希望单脉冲检测概率在 0.5 附近、积累后达到 0.9 以上这是雷达方程和信号处理增益共同决定的。值得注意的是CFAR 在所有距离单元上逐点滑窗处理时边缘位置的目标检测不到这是因为参考窗越界后无法有效估计背景功率。工程上可以在数据首尾补一段低电平噪声或者丢弃边缘若干个距离单元前者更适合仿真对比后者更适合实际系统因为真实系统里边缘距离单元往往对应近距离强杂波补噪声反而失真。本文还有配套的精品资源点击获取 雷达接收机MATLAB仿真：从回波生成到CFAR的模块化实现
简介这是一份基于MATLAB的雷达系统建模与仿真资料面向雷达信号处理初学者与研究者覆盖目标回波模拟、高频放大、混频、中频放大、相干检波、模数转换、脉冲压缩、MTI、积累和CFAR检测等完整信号处理链。压缩包共64个文件包含36个M脚本每个模块均配有可运行程序、对应ASV备份、11张JPG输出图形、两份DOC报告及FIG、LOG、DB等辅助文件整体约780KB结构清晰便于按模块对照学习。已有418人学习下载。借助这些脚本与图形读者可快速复现雷达收发链路的仿真结果理解各环节参数设置与算法原理也可在此基础上修改参数、扩展功能用于课程设计或科研验证。1. 为什么把雷达接收机拆成十个 MATLAB 模块拿到这套雷达系统建模与仿真的源码包时我先把几十个 .m 文件按接收链路顺序捋了一遍发现它覆盖的正是教材里雷达接收机的完整框图huibo.m 生成目标回波gaofang.m 做高频放大hunpin.m 混频zhongfang.m 中频放大xiangganjianbo.m 相干检波AD.m 模数转换maichongyasuo.m 脉冲压缩MTI.m 动目标指示jilei.m 多脉冲积累CFAR.m 恒虚警检测。这套东西适合两类人一是刚开始接触雷达信号处理需要把高放、混频、中放这些概念和代码一一对上的学生二是要做系统参数论证或课程设计的工程师想直接拿可运行的模板改参数。整个工程用 MATLAB 自带的信号处理函数就能跑通不依赖额外硬件每个模块输出的图形也都齐。下文按信号流逐级展开每级给出代码、参数依据和调参时容易踩的坑。2. 前端链路建模目标回波、增益分配与混频中放2.1 目标回波里该放哪些物理量huibo.m 要解决的第一件事是把目标的距离 R、径向速度 v、雷达截面积 σ 写进回波模型。若发射信号是线性调频脉冲理想的接收回波相对发射信号有三个变化延迟 τ 2R/c多普勒频移 fd 2v/λ幅度衰减正比于 R 的四次方σ 体现在回波功率上。很多初学者只做延迟不加多普勒结果后面的 MTI 和积累模块什么都看不到。包里还附带了 swerling.m这是 Swerling 起伏模型用来描述 σ 在帧间或脉冲间的随机性做蒙特卡洛统计时很有用定频仿真时可以先绕开。我建议回波生成只保留三个自由度延迟、多普勒、幅度。RCS 起伏在后续多帧统计时再叠加不要在一个脉冲里混入两种随机性否则出了问题很难定位。下面是简化后的回波生成核心代码% 参数声明 c 3e8; % 光速 m/s f0 10e9; % 载频 10 GHz lambda c / f0; % 波长 B 10e6; % 发射带宽 10 MHz T 50e-6; % 脉冲宽度 50 us fs 40e6; % 采样率 t 0 : 1/fs : T - 1/fs; s_tx exp(1j * pi * (B/T) * t.^2); % LFM 基带波形 R0 5000; v 30; % 目标距离 5 km径向速度 30 m/s tau 2 * R0 / c; % 回波延迟 fd 2 * v / lambda; % 多普勒频率 s_rx s_tx .* exp(-1j * 2 * pi * f0 * tau) ... .* exp(1j * 2 * pi * fd * t); % 延迟多普勒逻辑说明exp(-1j * 2 * pi * f0 * tau)是回波相对发射信号的固定相移它由目标距离决定这个相位在脉冲压缩前看不出作用但到 MTI 和相参积累时会直接影响相邻脉冲间的相位关系exp(1j * 2 * pi * fd * t)是目标运动引起的相位旋转这正是 MTI 和积累能区分动目标和静物回波的物理基础。发射信号用基带复包络表示而不是真的把 10 GHz 载频扔进采样是因为采样率不用跟着载频走这是雷达仿真里最常用的低通等效建模手法。提示如果严格按雷达方程算回波幅度需要给出发射功率、天线增益、噪声系数等一整套参数。仿真包里多数用归一化幅度再加高斯白噪声先把信噪比做成可调变量这对调试后端算法比追求绝对功率值更有价值。2.2 高放、混频、中放的增益与噪声分配高放 gaofang.m、混频 hunpin.m、中放 zhongfang.m 对应的是接收机的模拟前端。高放把天线收到的微弱回波放大到混频器能正常工作的电平混频把射频信号与本地振荡器信号相乘完成频率搬移中放再对中频信号做一次窄带放大和滤波。三级增益分配有一条基本约束高放增益不能太大否则强目标回波会把混频器推到饱和区产生交调失真但高放噪声系数又直接决定接收机整机灵敏度。级联噪声系数公式 F F1 (F2-1)/G1 (F3-1)/(G1·G2) 在这套仿真里的落地方式是每级放大后往信号里叠一份高斯白噪声噪声功率按该级带宽和增益系数折算而不是只做无噪声的纯放大。下面是一组常见的参数分配参考模块频率范围典型增益噪声系数带宽高放射频 9.4 GHz20 dB2 dB100 MHz混频射频转中频-6 dB8 dB100 MHz中放中频 30 MHz30 dB6 dB10 MHz相干检波基带 I/Q0 dB—10 MHz混频在数学上就是乘法s_if s_rf .* cos(2*pi*f_lo*t)MATLAB 里用点乘直接实现。注意混频输出不止有中频分量还叠加了和频分量实际仿真里要么在混频后加低通或带通滤波器要么直接用复信号让负频率分量不落在采样带宽内。如果看到 hunpin.m 的输出频谱出现两个对称峰值多半是漏了滤波这一级。中放的带宽要和发射信号带宽匹配带宽太窄会把 LFM 信号边缘削掉导致后续脉冲压缩的输出主瓣变宽。2.3 相干检波为什么必须保留 I/Q 两路中放输出到相干检波 xiangganjianbo.m这里最容易犯的错误是把检波理解成取包络。相干检波的工作原理是将中频信号分别与同相参考信号 cos(2πf_if·t) 和正交参考信号 -sin(2πf_if·t) 相乘经低通滤波后得到 I/Q 两路基带信号。这么做的意义有两个一是保留相位信息后续脉冲压缩的匹配滤波、MTI 的相位对消都依赖它二是信号变成复基带后采样率只需要覆盖信号带宽不需要扛着 30 MHz 的中频频率。% 相干检波核心I/Q 两路混频加低通 t 0 : 1/fs : N/fs - 1/fs; % N 为样本数 f_if 30e6; % 中频频率 x_i 2 * s_if .* cos(2 * pi * f_if * t); % 同相支路 x_q -2 * s_if .* sin(2 * pi * f_if * t); % 正交支路 % 设计低通滤波器通带 5 MHz 对应带宽一半 lpFilt designfilt(lowpassfir, PassbandFrequency, 5e6, ... StopbandFrequency, 8e6, SampleRate, fs); z filter(lpFilt, x_i 1j * x_q); % 合成复基带信号逻辑说明乘以 2 是工程惯例乘法器输出中差频项的系数只有 1/2乘 2 可以把幅度归一。低通滤波器去掉的是和频分量。PassbandFrequency 取 5 MHz对应 10 MHz 信号带宽的一半如果这里取得太大和频分量会残留基带信号会被污染。整条模拟前端调参的正确顺序是先固定噪声系数调整各级增益观察输出信噪比变化再回头看混频器和中放的带宽对脉冲沿的影响。信噪比异常时优先怀疑混频后的滤波而不是增益不够。3. ADC 量化与脉冲压缩距离分辨率的关键3.1 ADC 的采样率和量化位数怎么定AD.m 处理的是相干检波输出的模拟基带信号。两个参数决定这个环节的信息损失采样率和量化位数。基带信号带宽若为 B采样率取 fs ≥ B 就能满足复采样条件工程上通常会留 20% 裕量带宽 10 MHz 的信号采样率取 40 MHz 已经富余。量化位数决定信号动态范围仿真里默认取 12 位以上量化噪声远小于接收机噪声时可以忽略不需要刻意调低去测试量化效应。% 均匀量化模型 n_bits 12; % 量化位数 Vmax max(abs(z)); % 信号峰值 L 2^n_bits; % 量化电平数 delta 2 * Vmax / L; % 量化步长 z_q round(z / delta) * delta; % 量化后的复信号 % 量化信噪比理论值约 6.02*n_bits 1.76 dB sqnr_db 6.02 * n_bits 1.76;逻辑说明round(z / delta) * delta是均匀量化器的标准实现把连续幅值映射到最近的电平上。量化信噪比的理论公式 6.02n 1.76 dB 在这里可以作为参考12 位量化对应约 74 dB 的动态范围远超雷达接收链路中 -10 到 20 dB 的典型信噪比区间因此 ADC 不是系统瓶颈。如果发现量化后的仿真结果和理想情况差异明显先检查 Vmax 取值是否合理峰值取得过大或过小都会浪费量化区间。提示ADC 真实器件还存在采样抖动、孔径延迟和非线性这套仿真没有细化到器件级。做系统级论证时ADC 模型保持理想把注意力放在脉冲压缩、MTI、CFAR 这些算法链路上才是对的。3.2 脉冲压缩匹配滤波与 LFM 信号脉冲压缩是 maichongyasuo.m 的核心。宽脉冲保证足够的平均发射功率调制带宽决定距离分辨率接收端用匹配滤波把回波压缩成窄脉冲。LFM 信号的匹配滤波输出是一个 sinc 型脉冲主瓣宽度约 1/B因此距离分辨率 ΔR c/(2B)。这里有一个常见误解分辨率只由带宽决定和脉冲宽度无关带宽 10 MHz 对应分辨率约 15 m。匹配滤波有两种实现方式时域卷积和频域相乘。数据量大时频域实现更快因为一次 FFT 的复杂度是 O(N log N)而时域卷积是 O(N²)。代码如下% 脉冲压缩频域匹配滤波 Nfft 2^nextpow2(length(s_rx) length(s_tx) - 1); S_rx fft(s_rx, Nfft); S_ref fft(conj(fliplr(s_tx)), Nfft); % 匹配滤波器的频响 y ifft(S_rx .* S_ref, Nfft); y y(1 : length(s_rx) length(s_tx) - 1); % 归一化后取幅度 y_abs abs(y) / max(abs(y));逻辑说明匹配滤波器的冲激响应是发射信号的共轭时间反演对应代码里的conj(fliplr(s_tx))。频域相乘相当于时域卷积nextpow2保证 FFT 长度覆盖线性卷积输出长度避免循环卷积的混叠。输出主瓣宽度约 1/B若在频域给参考信号加窗比如汉明窗可以将距离旁瓣压低到 -40 dB 左右代价是主瓣展宽约 1.5 倍这个权衡在弱目标检测时经常要用到。表格里列出不同窗函数的典型旁瓣和主瓣展宽系数可以当作调参速查窗函数最高旁瓣电平主瓣展宽系数矩形窗不加窗-13.3 dB1.0汉宁窗-31.5 dB1.62汉明窗-42.7 dB1.47布莱克曼窗-58.1 dB1.73用窗的时机也很关键只在匹配滤波的频域参考信号上乘窗不要在回波域直接加窗否则会损失信噪比。4. MTI 对消与脉冲积累滤除静物回波4.1 一次对消与二次对消的传递函数分析MTI.m 实现的是动目标指示。它的思路建立在这样一个事实上静止地物回波在相邻脉冲间的相位保持不变而运动目标的相位按多普勒频率旋转。因此将相邻两个脉冲的回波相减静物回波被抵消运动目标保留下来这就是一次对消。在 z 域一次对消的传递函数是 H(z) 1 - z⁻¹频率响应在零多普勒处有一个零点。代码实现非常简洁% 一次对消相邻脉冲相减 mti_out z_pulse(2:end, :) - z_pulse(1:end-1, :); % 二次对消再对一次对消结果做一次差分 mti_out2 mti_out(2:end, :) - mti_out(1:end-1, :); % 频响对比多普勒频率轴 f_dop linspace(-PRF/2, PRF/2, 256); H1 abs(1 - exp(-1j * 2 * pi * f_dop / PRF)); H2 abs(1 - exp(-1j * 2 * pi * f_dop / PRF)).^2;逻辑说明z_pulse是脉冲维 × 距离维的二维矩阵z_pulse(2:end,:) - z_pulse(1:end-1,:)把相邻脉冲逐距离门相减。第一行代码做一次对消第二行做二次对消。频响表达式里 PRF 是脉冲重复频率零多普勒处的陷波越深对静止杂波的抑制越好。二次对消的陷波更宽适合处理有一定谱宽的地物杂波但代价是低速目标的响应也会被削掉一部分最小可检测速度会抬升。4.2 相参积累与非相参积累的取舍积累在 jilei.m 里实现作用是把多个脉冲的回波能量叠加提高信噪比。积累分两种相参积累在检波之前对复信号做非相参积累在检波之后对幅度做。相参积累理论上有 M 倍的 SNR 增益非相参积累只有大约 M 的 0.8 次方但相参积累要求目标在积累时间内不发生跨距离门走动并且多普勒频率必须落在滤波器通带内。常见做法是用 FFT 做相参积累把 M 个脉冲在距离维的复数据沿慢时间做 FFT得到距离-多普勒图% 相参积累沿慢时间维做 FFT M 64; % 积累脉冲数 rd_map fftshift(fft(z_matrix, M, 1), 1); % z_matrix: 脉冲×距离 % 频率轴换算 f_axis (-M/2 : M/2-1) * PRF / M; % 多普勒分辨率与积累时间的关系 df_dop PRF / M; v_res df_dop * lambda / 2; % 速度分辨率逻辑说明fft(..., M, 1)沿第一维脉冲维做 M 点 FFTfftshift把零多普勒移到频谱中心便于观察。多普勒分辨率等于 PRF/M速度分辨率再乘 λ/2。比如 PRF 为 1 kHz64 点积累多普勒分辨率约 15.6 HzX 波段速度分辨率约 0.23 m/s。积累脉冲数 M 不能贪多目标在积累时间内如果发生跨距离门走动能量会散到多个距离门里增益反而下降。遇到高速目标时要用包络走动补偿或者改用更短的积累时间。MTI 和积累的配合关系是先做 MTI 对消掉强静物杂波再做积累提高目标信噪比。顺序不能反如果先积累再做 MTI强杂波会把积累后的动态范围撑爆弱小目标就检测不出来了。5. CFAR 检测门限的工程设定与实测验证5.1 CA-CFAR 的参数依据CFAR.m 实现的是单元平均恒虚警检测。基本原理是检测单元两侧各取 N 个参考单元求平均作为杂波功率估计再乘以门限因子 α 得到检测门限。参考单元数和保护单元数是最关键的两个参数。保护单元紧邻检测单元用来挡住目标回波主瓣和旁瓣泄漏进参考单元导致门限被抬高。对脉压后的 LFM 信号主瓣大约占 2~3 个距离单元保护单元每侧取 2~4 个即可。门限因子 α 由虚警概率 Pfa 和参考单元数 N 决定关系式是 α N(Pfa^(-1/N) - 1)。它的推导基于参考单元噪声服从指数分布这一假设实际仿真里如果噪声是高斯分布经平方率检波后也近似成立。下面给出一段可运行的 CFAR 检测核心代码function det cfar_1d(x, N_ref, N_guard, Pfa) % x: 输入距离维数据, N_ref: 参考单元数(每侧), N_guard: 保护单元数(每侧) alpha N_ref * (Pfa^(-1/N_ref) - 1); det zeros(size(x)); len length(x); for k N_guard N_ref 1 : len - N_guard - N_ref ref_left x(k - N_guard - N_ref : k - N_guard - 1); ref_right x(k N_guard 1 : k N_guard N_ref); noise_est (sum(ref_left) sum(ref_right)) / (2 * N_ref); threshold alpha * noise_est; det(k) x(k) threshold; end end逻辑说明循环从头到尾滑动对每个距离单元 k取左侧和右侧各 N_ref 个参考单元求平均得到噪声功率估计noise_est门限是 α 乘以噪声估计。det(k) x(k) threshold输出二值判断结果。Pfa 取 1e-6 到 1e-4 之间是雷达仿真的常见设置Pfa 定得过低会大幅抬高门限导致检测概率下降过高则会出现大量虚警。参考单元数 N_ref 取 16 到 32 之间较稳妥太小则噪声估计方差大门限抖动明显太大则在非均匀杂波环境中响应迟缓。5.2 用蒙特卡洛验证虚警率与检测概率评价 CFAR 参数设得合不合理最直接的办法是跑蒙特卡洛仿真。对纯噪声输入统计虚警数对含目标输入统计检测概率两条曲线能直观看出门限因子和参考窗长度的权衡。% 蒙特卡洛验证虚警率统计 Nmc 1000; N_cells 1000; cnt 0; for i 1:Nmc noise abs(randn(1, N_cells) 1j*randn(1, N_cells)).^2 / 2; det cfar_1d(noise, 16, 4, 1e-4); cnt cnt sum(det); end pfa_measured cnt / (Nmc * N_cells); fprintf(理论 Pfa: 1e-4, 实测 Pfa: %e\n, pfa_measured);逻辑说明randn生成复高斯噪声取模平方除以 2 得到指数分布噪声功率样本这是平方率检波后的标准噪声模型。实测虚警率和理论值通常在一个数量级内吻合如果偏差过大优先检查保护单元是否过少或参考单元内的强杂波没有被清除。检测概率的仿真则是把已知信噪比的信号插进噪声统计检测到目标的次数占比通常希望单脉冲检测概率在 0.5 附近、积累后达到 0.9 以上这是雷达方程和信号处理增益共同决定的。值得注意的是CFAR 在所有距离单元上逐点滑窗处理时边缘位置的目标检测不到这是因为参考窗越界后无法有效估计背景功率。工程上可以在数据首尾补一段低电平噪声或者丢弃边缘若干个距离单元前者更适合仿真对比后者更适合实际系统因为真实系统里边缘距离单元往往对应近距离强杂波补噪声反而失真。本文还有配套的精品资源点击获取