简介：这份资源围绕雷达发射信号、雷达回波信号与脉冲压缩三大核心概念展开，面向雷达信号处理初学者、通信电子相关专业学生及需要快速搭建仿真验证环境的工程人员，帮助理解匹配滤波、回波分析与脉冲压缩的实现思路。压缩包共6个文件，均为m脚本文件，整体约8KB，涵盖主程序入口、回波生成、发射信号构造、匹配滤波与脉冲压缩等模块，便于在MATLAB中直接运行与修改参数。已有296人学习下载，适合作为课程设计或自学练手的参考素材。读者可借助源码观察发射信号参数对探测性能的影响，理解回波信号中距离、速度信息的提取过程，并掌握匹配滤波器与发射信号保持一致的设计要点，从而把抽象理论落到可运行的代码上，为后续雷达系统仿真与性能优化打下基础。
1. 从一组 MATLAB 源码看雷达回波仿真：脉冲压缩到底在算什么
如果你正在做雷达信号处理的课程设计、算法预研，或者单纯想跑通一次从发射到脉压的完整链路，这套rader2源码包值得拆开看看。它把雷达发射信号生成、目标回波模拟、匹配滤波脉冲压缩这几件事拆成了几个独立的.m文件，入口是main1105.m，配套rect.m、signal_tran1101.m、Signal_transmit.m、Signal_match.m、Signal_echo.m。说白了，这就是一个用 MATLAB 搭起来的最小雷达信号处理闭环：你改几个参数，就能看到发射波形长什么样、回波长什么样、脉压之后主瓣有多窄、旁瓣有多高。
它适合谁？一是刚接触雷达信号处理、需要一份能跑通的参考实现来对照课本公式的人；二是已经懂 LFM、匹配滤波原理，但想快速验证不同脉宽、带宽、采样率组合下脉压效果的人。不适合指望它直接上硬件或者做工程级实时处理的人——这是仿真验证代码，不是产品级 DSP 工程。下面我按“先跑通、再拆解、最后避坑”的顺序，把这几个文件的关系和参数逻辑讲清楚。
2. 发射、回波、脉压三段式：文件依赖与调用链拆解
2.1 入口脚本 main1105.m 在做什么
拿到源码包，第一件事不是急着改参数，而是先理清调用关系。main1105.m是主入口，它负责定义雷达系统的基本参数，然后依次调用发射信号生成、回波模拟、匹配滤波三个环节。常见做法是：主脚本里集中放参数，子函数只负责计算，这样调参时不用满文件找变量。
一个典型的入口脚本结构大致是这样：
% main1105.m 主入口：定义参数并串联发射-回波-脉压流程 clear; clc; close all; %% 雷达系统参数 c = 3e8; % 光速 (m/s) fc = 10e9; % 载频 10GHz (X波段) Tp = 10e-6; % 脉冲宽度 10us B = 5e6; % 信号带宽 5MHz fs = 20e6; % 采样率 20MHz，满足 fs > 2B R_target = 3000; % 目标距离 3km PRF = 1000; % 脉冲重复频率 1kHz %% 生成发射信号 [st, t] = signal_tran1101(Tp, B, fs, fc); %% 生成目标回波 sr = Signal_echo(st, t, R_target, c, fs); %% 匹配滤波脉冲压缩 [pc, t_pc] = Signal_match(st, sr, fs); %% 绘图 figure; subplot(3,1,1); plot(t*1e6, real(st)); title('发射信号实部'); xlabel('时间/us'); subplot(3,1,2); plot(t*1e6, real(sr)); title('回波信号实部'); xlabel('时间/us'); subplot(3,1,3); plot(t_pc*1e6, abs(pc)); title('脉压输出'); xlabel('时间/us');这段代码的逻辑很直白：先定参数，再按顺序调三个函数，最后画三张图。参数里最需要留意的是采样率fs和带宽B的关系——fs必须大于2B，否则回波频谱会混叠，脉压结果直接失真。Tp和B的乘积就是时宽带宽积，这个值越大，脉压后的压缩比越高，主瓣越窄，但计算量也越大。
2.2 发射信号生成：signal_tran1101.m 与 Signal_transmit.m 的分工
这两个文件名字很像，容易搞混。从命名习惯看，signal_tran1101.m更可能是带版本号的发射信号生成函数，而Signal_transmit.m可能是更通用的发射信号模块。实际使用时，先看main1105.m里调的是哪个，以调用为准。
发射信号的核心是线性调频（LFM）信号，数学表达是：
% signal_tran1101.m 线性调频信号生成 function [st, t] = signal_tran1101(Tp, B, fs, fc) % Tp: 脉冲宽度, B: 带宽, fs: 采样率, fc: 载频 N = round(Tp * fs); % 采样点数 t = (0:N-1) / fs; % 时间轴 K = B / Tp; % 调频斜率 st = exp(1j * pi * K * t.^2) .* exp(1j * 2 * pi * fc * t); end这里K = B/Tp是调频斜率，决定了频率在脉冲内线性变化的快慢。exp(1j*pi*K*t.^2)是基带 LFM 项，后面乘上exp(1j*2*pi*fc*t)是把信号搬移到载频。实际仿真里，如果只关心基带处理，载频项可以省掉，直接分析复包络，计算量更小。
rect.m大概率是矩形窗函数，用来限定脉冲的时间范围。有些实现会用rectpuls或者直接逻辑索引，效果一样。如果你发现发射信号首尾有突变导致的频谱泄漏，可以在rect.m里加一点升余弦滚降，但会牺牲一点分辨率。
2.3 回波模拟：Signal_echo.m 的距离延迟与衰减
回波信号本质就是发射信号延迟一段时间、幅度衰减一下。延迟由目标距离决定：tau = 2*R/c。Signal_echo.m要做的就是把这个延迟加到发射信号上。
% Signal_echo.m 目标回波生成 function sr = Signal_echo(st, t, R, c, fs) % st: 发射信号, t: 时间轴, R: 目标距离, c: 光速, fs: 采样率 tau = 2 * R / c; % 双程延迟 N_delay = round(tau * fs); % 延迟采样点数 sr = zeros(size(st)); if N_delay < length(st) sr(N_delay+1:end) = st(1:end-N_delay); end % 幅度衰减，与距离平方成反比 sr = sr / (R^2); end这段代码里有两个容易翻车的点。一是N_delay取整带来的距离量化误差，如果fs不够高，延迟点数取整后距离精度会下降。二是幅度衰减用1/R^2只是简化模型，实际雷达方程里还有天线增益、目标 RCS、系统损耗等因子，仿真里如果只关心脉压形状，衰减系数可以归一化掉。
2.4 匹配滤波：Signal_match.m 如何实现脉冲压缩
匹配滤波是脉冲压缩的核心。对 LFM 信号来说，匹配滤波器的冲激响应就是发射信号的时间反转共轭。在频域做就是发射信号频谱的共轭相乘。
% Signal_match.m 匹配滤波脉冲压缩 function [pc, t_pc] = Signal_match(st, sr, fs) % st: 发射参考信号, sr: 回波信号, fs: 采样率 N = length(st) + length(sr) - 1; Nfft = 2^nextpow2(N); % 取2的幂次加速FFT St = fft(st, Nfft); % 发射信号频谱 Sr = fft(sr, Nfft); % 回波信号频谱 H = conj(St); % 匹配滤波器频响 Pc = Sr .* H; % 频域相乘 pc = ifft(Pc, Nfft); % 反变换回时域 t_pc = (0:Nfft-1) / fs; endconj(St)就是匹配滤波器的频域形式。Nfft取 2 的幂次是为了 FFT 加速，这是工程上的常规操作。脉压输出的峰值位置对应目标距离，主瓣宽度约等于1/B，旁瓣电平对于矩形窗 LFM 大约是 -13.2dB。如果你看到旁瓣高得离谱，先检查Nfft是不是太小导致循环卷积混入了线性卷积的结果。
3. 参数怎么设：脉宽、带宽、采样率的联动关系
3.1 时宽带宽积决定压缩比
脉冲压缩的收益来自时宽带宽积D = Tp * B。发射一个宽脉冲保证能量，脉压后等效成窄脉冲保证分辨率。压缩后的等效脉宽约1/B，所以距离分辨率ΔR = c/(2B)。
举个例子：B = 5MHz时，ΔR = 3e8/(2*5e6) = 30m。如果你想要 15m 的分辨率，带宽得翻倍到 10MHz。但带宽增大后，采样率fs也得跟着提，否则混叠。常见做法是fs = 2.5B ~ 4B，留够过渡带。
| 参数 | 符号 | 典型取值 | 影响 |
|---|---|---|---|
| 脉冲宽度 | Tp | 10~100us | 越大能量越强，但盲区越大 |
| 带宽 | B | 1~20MHz | 决定距离分辨率 |
| 采样率 | fs | 2.5B~4B | 不足会混叠，过高浪费算力 |
| 调频斜率 | K | B/Tp | 影响频谱展宽 |
| 载频 | fc | 1~10GHz | 决定波长和穿透性 |
3.2 采样率不够会怎样
采样率低于2B时，回波频谱发生混叠，脉压输出的主瓣会展宽、旁瓣会抬高，严重时目标峰值位置都会偏移。我一般会在main1105.m里加一句检查：
% 采样率合规检查 if fs < 2 * B warning('采样率 %.2f MHz 低于 2B=%.2f MHz，存在混叠风险', fs/1e6, 2*B/1e6); end这个检查不解决根本问题，但能让你在跑出奇怪结果之前先意识到参数不匹配。
3.3 目标距离与延迟点数
目标距离R对应的延迟tau = 2R/c。如果tau超过一个脉冲重复周期1/PRF，就会出现距离模糊——第二个脉冲的回波还没收完，第一个脉冲的回波已经跑到了下一个周期。仿真里如果只发一个脉冲，这个问题不明显，但做多脉冲积累时就会暴露。
N_delay = round(tau * fs)的取整误差对应距离误差ΔR_q = c/(2*fs)。fs = 20MHz时，ΔR_q = 7.5m。如果目标距离精度要求高于这个值，要么提高采样率，要么在回波生成时用插值而不是取整。
4. 避坑与排查：跑不通、结果不对时先看这几条
4.1 现象：脉压输出没有明显峰值
原因通常是回波和参考信号的时间轴没对齐，或者匹配滤波时用了错误的参考信号。Signal_match.m里的st必须是和回波对应的发射信号副本，不能是另一个参数下生成的。
解决：在main1105.m里确认st只生成一次，然后同时传给Signal_echo和Signal_match。如果两个函数内部各自重新生成了发射信号，参数稍有出入就会导致失配。
4.2 现象：脉压主瓣很宽，分辨率不达标
先查带宽B是不是设小了。分辨率ΔR = c/(2B)，B = 1MHz时分辨率是 150m，B = 10MHz时才是 15m。另一个可能是Nfft太小，频域相乘变成了循环卷积，主瓣被展宽。
解决：把Nfft设为2^nextpow2(length(st)+length(sr)-1)，确保线性卷积不混叠。
4.3 现象：旁瓣电平异常高
矩形窗 LFM 的旁瓣理论值约 -13.2dB。如果实测远高于这个值，检查发射信号有没有加窗。有些实现会在rect.m里默认加汉明窗，这会把旁瓣压到 -40dB 以下，但主瓣会展宽约 1.5 倍。如果你要的是理论对照，就别加窗；如果要低旁瓣，就接受主瓣展宽。
4.4 现象：回波幅度为零或全为零
Signal_echo.m里N_delay如果大于等于length(st)，延迟后的信号就完全移出了时间窗，回波全零。这通常发生在目标距离设得太大、而时间轴长度不够的情况下。
解决：把时间轴长度设为Tp + 2*R_max/c，确保最远目标的回波也落在窗内。或者直接在main1105.m里根据最大探测距离反推需要的采样点数。
4.5 现象：FFT 之后频谱不对称
如果发射信号是实信号，频谱应该关于零频共轭对称。如果用了复指数载频项，频谱会搬移到fc处，这是正常的。但如果基带 LFM 的频谱不对称，检查K的符号和t的起始点。t从 0 开始还是从-Tp/2开始，会影响相位参考，但对脉压幅度没影响。
5. 进阶验证：用模糊函数和加窗对比确认脉压质量
跑通基本流程之后，我一般会做两件事来确认这套代码的脉压结果是否可信。第一是画模糊函数，第二是对比不同窗函数下的旁瓣表现。
模糊函数描述的是信号在距离-速度二维上的分辨能力。对 LFM 信号，模糊函数是斜椭圆，斜率和调频斜率相关。在 MATLAB 里可以用xcorr快速近似：
% 用自相关近似零多普勒切面的模糊函数 [af, lag] = xcorr(st, st); af = abs(af) / max(abs(af)); lag_us = lag / fs * 1e6; figure; plot(lag_us, 20*log10(af)); xlabel('延迟/us'); ylabel('幅度/dB'); title('LFM 自相关（零多普勒切面）'); grid on;这条曲线的主瓣宽度对应1/B，第一旁瓣应该在 -13.2dB 附近。如果对不上，说明发射信号生成有问题。
加窗对比更直接。在匹配滤波之前给参考信号乘一个汉明窗：
% 加窗匹配滤波对比 win = hamming(length(st))'; st_win = st .* win; [pc_win, ~] = Signal_match(st_win, sr, fs);然后和矩形窗的结果画在一起，看主瓣展宽了多少、旁瓣压了多少。汉明窗通常能把旁瓣压到 -40dB 以下，但主瓣展宽约 1.5 倍。这个取舍在雷达系统设计里很常见：要低旁瓣就得牺牲一点分辨率。
还有一个验证点是脉压增益。理论脉压增益是10*log10(D)，D = Tp*B。Tp = 10us、B = 5MHz时，D = 50，增益约 17dB。你可以在脉压前后分别算信号峰值和噪声均方根，看信噪比提升了多少，和理论值对一下。如果差太多，检查回波里有没有加噪声、匹配滤波器是不是真的匹配。
从那以后我每次拿到新的雷达仿真代码，都强制先跑一遍单目标、无噪声、矩形窗的基准 case，确认主瓣宽度和旁瓣电平和理论对得上，再去改参数做复杂场景。这一步花不了几分钟，但能省掉后面大量“结果不对但不知道哪不对”的排查时间。希望帮到你。
本文还有配套的精品资源，点击获取
