简介：这套 MATLAB 源码以雷达波形调频率捷变技术为核心，面向通信、电子信息工程、数学等专业的学生，以及需要快速验证抗干扰算法的研究人员，解决多目标干扰场景下雷达波形设计与频率捷变仿真的问题。代码采用参数化编程，关键参数可灵活修改，注释细致，并附有可直接运行的案例数据，能够帮助新手快速掌握从波形生成到抗干扰性能评估的完整流程。压缩包共 15 个文件，包含 9 个 m 脚本、3 个 asv 备份文件、2 个 mat 数据文件及 1 张示意图，整体约 276KB，结构清晰，适合课程设计、期末大作业或毕业设计使用。目前已有 150 人学习下载。通过该资源，使用者可以直观了解频率捷变如何削弱多目标干扰，并基于提供的脚本继续扩展不同调制方式、观察参数变化对信噪比和抗干扰性能的影响，为后续深入研究雷达电子对抗技术打下基础。

1. 雷达波形调频率捷变抗多目标干扰：rar 里的 matlab 仿真链路到底能干什么

雷达波形调频率捷变抗多目标干扰的 matlab 代码，解压下来以后不是只看一张频谱图就完事，而是一条从波形参数到干扰评估的完整仿真链路。多目标干扰一进来，固定载频脉冲串基本白给。频率捷变的核心是让每个脉冲的载频按伪随机规律切换，干扰机很难提前猜出下一个点频，多假目标没法在每个脉冲上保持同相积累，脉压后再做多普勒处理，真目标进得去、假目标散得开。

这份 rar 适合两类人：刚接抗干扰课题、需要快速出验证波形的同学，和做电子对抗仿真的工程师。代码本身并不复杂，甚至可以说有点“原始”：一排放射脉冲、一段回波、一次匹配滤波和二维 FFT。但边界条件很多，频率集间隔、相位补偿、统计次数，全都踩一遍才会有手感。我建议把它当成一个能跑通的基准，先不动算法，只改随机种子和观测时长做对比，后面再谈复杂干扰。

2. 频率捷变不止是改载频：从波形参数到跳频序列

当你把 rar 解压完，第一次搜索代码里哪个变量是“载频”，会发现它多半不是单个数字，而是一列数。这就是频率捷变的第一印象。但“改载频”和“做频率捷变”之间的差距，比想象中大。关键不在频点多，而在干扰机对下一脉冲频率的预测难度。

2.1 固定载频为什么怕多目标干扰

多假目标干扰的典型做法是 DRFM 转发：干扰机截获雷达信号后存储脉内波形，再按预设的多个延迟时刻重新发射。只要载频确定，干扰机测一次频率就能在所有脉冲上保持同相转发，造假目标的效率非常高。接收端做匹配滤波时，每一个假目标都对应一个相关峰，做 MTD 之后这些峰还会同步积累，看起来就像一群真目标。

问题本质不是脉冲压缩做得不够好，而是干扰机拥有“复制并重放”的能力。频率捷变打破的恰恰是复制条件。当每个脉冲的载频都变化，干扰机要么提前不知道频点，要么只能复现上一个脉冲的波形。复制波形与当前脉冲载频不匹配，假目标经过匹配滤波后不会形成规则峰值，多普勒维上也会散开。这就是频率捷变抗多假目标干扰的基本原理。

2.2 每脉冲跳频与脉组跳频的工程取舍

工程实现上有两种常见方式。第一种是 PRI 级跳频，每个脉冲都换频，对抗性最强，但对频率合成器的稳定时间、接收机本振的跟踪速度都有要求，脉压之后还要做相位补偿。第二种是 CPI 级跳频，一个 CPI 内所有脉冲用同一个频点，下一个 CPI 再跳走，处理简单，但干扰机有一个 CPI 的适应时间，抗干扰能力弱一档。

我一般会建议先用脉组跳频把流程跑通，再切到每脉冲跳频看收益差异。代码里留一个模式开关最省事：

agile_mode = 'pulse'; % 'pulse'：每个脉冲都跳；'cpi'：一个CPI内同频

当agile_mode等于'pulse'时，每个脉冲从频率集里独立抽取；等于'cpi'时，只在 CPI 边界重新抽一次。后面讲到相位补偿时会看到，脉组跳频几乎不需要补偿，每脉冲跳频则必须补偿。

2.3 频率集设计的三条硬约束与可复现步骤

选频率集有三个硬约束，缺一个都会在仿真里翻车。

第一，跳频间隔不小于信号带宽。带宽 10MHz，跳频间隔至少 10MHz，实际留一倍余量取 20MHz 更稳妥。第二，跳频覆盖范围要控制在收发系统瞬时带宽内，否则有的频点接收机根本处理不了。第三，伪随机序列周期不能短于观察时间，否则一串脉冲里重复频点太多，干扰机又有机会学习规律。

复现步骤很简单：先定频率集，再锁随机种子，最后看跳频序列。下面的代码是解除 rar 之后我通常会先跑的第一段：

rng(2025); % 固定随机种子，保证同学之间可复现 fc_set = (9.3e9:20e6:10.1e9); % 41 个频点，间隔 20MHz freq_seq = fc_set(randi(length(fc_set), 64, 1)); plot(freq_seq / 1e9, '.-'); xlabel('脉冲序号'); ylabel('载频 (GHz)'); title('跳频序列');

这里randi(length(fc_set), 64, 1)生成 64 个随机索引，对应一个 CPI 内 64 个脉冲。把64改成128就是 128 脉冲。执行结果会看到频率在 9.3GHz 到 10.1GHz 之间乱跳，没有规律。这就是后面所有仿真的源头：随机种子不变，跳变序列就完全一致，任何人复现都是同一个结果。

3. matlab 代码拆解：回波生成、匹配滤波与二维 FFT 的具体实现

这一章把 rar 里的主流程拆成三段：参数定义、回波与干扰生成、匹配滤波和多普勒处理。只想改参数，看 3.1 和 3.4；要加新干扰样式，在 3.2 的回波叠加处下手。

3.1 主脚本里的参数定义与跳频序列生成

仿真最开始是一段参数定义，频率捷变的“频率”就在这里产生。我习惯把所有能改的物理量先集中放在文件头部，避免后面为找一个常数翻半天。

c = 3e8; B = 10e6; % 线性调频带宽 10MHz Tp = 10e-6; % 脉宽 10us fs = 20e6; % 基带采样率 20MHz prf = 10e3; % 脉冲重复频率 10kHz npulse = 64; % 一个CPI内脉冲数 fc_set = 9.3e9:20e6:10.1e9; rng(2025); freq_seq = fc_set(randi(length(fc_set), npulse, 1)); Ns = round(Tp * fs); % 每个脉冲的采样点数 t = (0:Ns - 1) / fs; mu = B / Tp; % LFM调频率 tx = exp(1j * pi * mu * t.^2); % 发射基带信号

采样率fs只要大于带宽两倍就够，20MHz 对应带宽 10MHz 是两倍过采样。freq_seq是跳频序列，每个脉冲的载频都存在这里。tx是发射信号的复基带形式，后面匹配滤波也要用它。

3.2 目标回波与延迟转发干扰建模

接下来是主体循环。这段代码把真目标、复制式干扰和热噪声压进同一个脉冲里，再逐脉冲做匹配滤波。注意干扰机只能复现前一个脉冲，所以它的相位项用的是freq_seq(k-1)。

R1 = 100; tau1 = 2 * R1 / c; D1 = round(tau1 * fs); R2 = 185; tau2 = 2 * R2 / c; D2 = round(tau2 * fs); fd = 120; pri = 1 / prf; data = zeros(Ns, npulse); for k = 1:npulse f0 = freq_seq(k); % 当前脉冲的真实载频 rx = zeros(1, Ns); % 真实目标：延迟D1，载频相位随 f0 变化 if D1 < Ns rx(D1 + 1:end) = tx(1:end - D1) * exp(-1j * 2 * pi * f0 * tau1) * exp(1j * 2 * pi * fd * k * pri); end % 复制式干扰：只能复现前一个脉冲，因此用上一脉冲载频 if k > 1 && D2 < Ns rx(D2 + 1:end) = rx(D2 + 1:end) + tx(1:end - D2) * exp(-1j * 2 * pi * freq_seq(k - 1) * tau2); end rx = rx + 0.02 * (randn(1, Ns) + 1j * randn(1, Ns)); % 热噪声 % 匹配滤波 y = conv(rx, conj(fliplr(tx)), 'same'); % 关键补偿：以真实目标参考距离 D1 对应的 tau1 做相位补偿 y = y * exp(1j * 2 * pi * f0 * tau1); data(:, k) = y; end

逻辑说明：真目标回波里，exp(-1j*2*pi*f0*tau1)是载频相位，exp(1j*2*pi*fd*k*pri)是多普勒相位。干扰机延迟转发，但它没有跟上当前跳变，只能使用上一脉冲载频，所以它的跨脉冲相位是乱的。匹配滤波用conv(rx, conj(fliplr(tx)), 'same')，本质是把接收信号和发射基带信号的镜像卷积在一起，峰值出现在目标延迟处。

最后一行y * exp(1j*2*pi*f0*tau1)是频率捷变必须做的相位补偿。如果不补偿，真目标的多普勒谱会被捷变相位调制掉，看起来也像假目标。补偿以参考距离tau1为准，工程上对应“距离门级相位补偿”。

3.3 距离-多普勒二维 FFT

逐脉冲处理完之后，data是距离门乘脉冲矩阵。对每一行沿着脉冲维做 FFT，就得到距离-多普勒图。

RD = fftshift(fft(data, npulse, 2), 2); rd_amp = abs(RD); imagesc(0:Ns - 1, (-npulse / 2:npulse / 2 - 1) * prf / npulse, 20 * log10(rd_amp + 1e-12)); xlabel('距离门'); ylabel('多普勒 (Hz)'); colorbar;

这里fft(data, npulse, 2)沿第二维做 64 点 FFT，fftshift把零频移到中间。多普勒分辨率是prf/npulse，也就是 10000/64，约 156Hz。目标多普勒 120Hz 会落在靠近中心的谱线附近，假目标因为相位不连续，能量不会集中到一根谱线上。

3.4 参数修改的四个入口

做实验时主要改四处。第一是fc_set，它决定频率覆盖范围和频点数量；第二是npulse，它决定多普勒分辨率和积累时间；第三是fd，目标速度对应多普勒；第四是干扰幅度系数 0.02，它决定信干比基线。每次只动一个变量，不要同时改，否则无法判断结果变化来自哪个因素。

4. 避坑与排查：频率捷变仿真里常见的 5 个翻车点

这一章是解包这些代码最想写的血泪经验。频率捷变看着只有一行跳频，实际跑起来坑非常多，下面五条每一条我都踩过。

4.1 跳频间隔设得太小：距离维出现频谱重叠

现象：把fc_set的间隔设成 5MHz，信号带宽也是 10MHz，跑出来的脉压旁瓣明显抬高，目标峰周围多出一圈干扰状波纹。

原因：相邻频点间隔小于信号带宽，两个频点对应的频谱重叠，脉压输出互相污染。频率捷变本来想让干扰机猜不到，结果频谱先把自己人糊在一起。

解决：跳频间隔至少取信号带宽的 1 到 2 倍。10MHz 带宽就用 20MHz 间隔，宁可频点数少一点，也不要让频谱重叠。

4.2 不补偿载频相位：真目标在多普勒维被“捷变”掉了

现象：真目标在距离维的峰还在，但跨脉冲做 FFT 后多普勒峰值很弱，甚至找不到目标。第一次遇到时我以为目标参数写错了，查了半天才发现是跳频相位没处理。

原因：每个脉冲的载频不同，目标回波的初始相位也不同，这个相位不是多普勒带来的，而是exp(-j*2*pi*fc*tau)。不做补偿，它会在多普勒维形成随机调制。

解决：在匹配滤波后乘exp(+j*2*pi*fc*tau_ref)，tau_ref是参考距离对应的时延。先粗测目标距离，再把相位补偿补到那一格距离门上。补偿距离选错了，目标仍然会散开。

4.3 固定载频与捷变载频对比不公：干扰电平设置不一致

现象：跑完对比，捷变载频抗干扰效果很好，但换个干扰功率数字，效果又变差了，得不出稳定结论。

原因：两种模式里干扰幅度、噪声功率、目标相位都没有严格对齐。比较抗干扰收益必须先固定发射功率、干扰功率和目标散射系数，只改变载频是否跳变。

解决：把干扰幅度、噪声功率和目标幅度写成同一个变量，用蒙特卡洛跑至少 50 次，取平均信干比，不要只看一次仿真图像。统计口径一致，结论才有说服力。

4.4 解压后中文注释乱码：编码问题和 rar 解密是两回事

现象：rar 解开以后，.m文件里的中文注释全部变成乱码，代码本身能跑，但注释没法读。

原因：压缩包内文件的编码和 MATLAB 默认编辑器编码不一致，常见是源文件用 UTF-8，而 MATLAB 在旧环境下按本地编码打开。

解决：先用记事本或 VS Code 打开文件另存为 UTF-8 带 BOM，再放回 MATLAB。如果 rar 本身设了密码，不要硬猜，先找来源方要口令；确实忘了可以使用 Advanced RAR Password Recovery 这类工具尝试，但对长口令成功率很低，别抱太大希望。重点还是先把编码理顺。

4.5 循环写得太重：64 个脉冲能跑，512 个脉冲卡死

现象：把npulse改成 512 以后，脚本跑了十几分钟还没出结果，CPU 占用很高。

原因：逐脉冲做conv，每个脉冲都重新生成一段接收向量，512 个脉冲就是 512 次全数组卷积，MATLAB 在循环里还会不断重新分配内存。

解决：先把匹配滤波器换到频域，一次生成H = fft(conj(fliplr(tx)), Nfft)，再用ifft(fft(rx, Nfft) .* H)实现匹配滤波。更快的做法是把所有脉冲的回波排成矩阵，一次做二维 FFT 滤波。速度提升非常明显，结果和时域卷积一致。

5. 验证与扩展：把抗干扰收益量化成一个 dB 数字

5.1 蒙特卡洛 SINR 评估

频率捷变好不好，不能只看距离-多普勒图里假目标少了。我习惯用信干噪比 SINR 来量化，把抗干扰收益变成能写进报告的数字。

mc = 50; sinr_fix = zeros(1, mc); sinr_jump = zeros(1, mc); for m = 1:mc % 固定载频：所有脉冲 fc 都取 9.5GHz % 频率捷变：freq_seq 每个脉冲随机抽取 % 两种模式使用同一目标和同一干扰幅度 % 得到 RD 图后统计目标峰值与全局中值之比 signal_peak = max(max(abs(RD(D1-2:D1+2, :)))); noise_floor = median(abs(RD(:))); sinr = 20 * log10(signal_peak / noise_floor); end % 比较 mean(sinr_jump) 与 mean(sinr_fix)

这段代码里D1-2:D1+2是目标所在距离门附近取峰值，避免旁瓣被当成目标。RD固定载频版本和捷变版本要在同一套参数下各跑 50 次，最后比较均值。捷变带来的 SINR 增益就是两个均值之差。

5.2 与 LFM、自适应跳频组合

频率捷变可以单独用，也可以和脉内线性调频组合。LFM 提供距离分辨率，跳频提供抗干扰能力，两者不冲突。更进阶的做法是自适应跳频：先对每个频点上的干扰能量做估计，下一轮只选干扰能量最低的几个频点。这样不再是盲目随机跳，而是“避开拥挤频段”的智能跳变。

仿真里验证自适应跳频，可以把每根频点的干噪比记录成向量，排序后取后面一半作为新的fc_set。这套代码里只需要改频率集生成逻辑，不需要动匹配滤波和 MTD 部分。论文出图时，我会把距离-多普勒图导出成 EPS 矢量图再放进文档，比位图清晰得多。

频率捷变仿真最怕的就是“看起来有效，但说不清有效多少”。从那以后我每次拿到一个雷达仿真 rar，都强制先把固定载频的基线 SINR 统计出来，再用同一套统计口径去跑所有改进算法。这条流程救过我很多次，希望帮到你。

本文还有配套的精品资源，点击获取