简介适用于雷达电子战与干扰技术研究的MATLAB脚本围绕切片干扰Chopping and Interleaving算法提供完整实现帮助工程师和研究人员直观理解通过切割、重排雷达信号来破坏敌方探测与跟踪的机制。压缩包内仅含1个m文件大小1KB属于轻量级源码便于直接运行查看效果并在此基础上修改参数或扩展功能。目前已有1051人学习下载适合作为高校雷达课程、电子战培训或抗干扰算法设计的教学示例与参考起点。脚本内部涵盖信号生成、切片分割、交错组合、干扰效果评估等关键环节其中信号生成部分支持脉冲序列、频率调制、相位编码等常见雷达波形切片与交织模块可按预设规则或随机方式重组信号片段读者可对照代码梳理雷达干扰的典型流程也可将其嵌入SystemVue等系统仿真环境进行硬件在环验证进一步评估误检率、漏检率等干扰性能指标。1. ChoppingAndInterleaving一种面向线性调频雷达的切片重排干扰在工程仿真里复制一个线性调频脉冲然后整体转发脉压匹配滤波后在固定延迟处只会出现一个又高又窄的单峰稍微有点经验的雷达操作员都能把它从真实回波里区分出来。ChoppingAndInterleaving后文统一写 CI换了一种思路把截获到的雷达脉冲在时间轴上切成若干段再按新的顺序拼接发射。切片宽度和重排间隔设计得当的话匹配滤波输出端会出现一串幅度相近、位置可控的假目标并且这串假目标还能在多次积累后保持稳定干扰效果比单峰转发更贴合工程需要。下面直接从信号结构拆起用 MATLAB 从基带脉冲开始把 CI 的完整链路写一遍给出一套能改参数复现的脚本再讨论切片宽度、重排间隔和相位连续性问题。适合正在做雷达信号级仿真、DRFM 干扰建模或电子对抗效果评估的工程师参考。2. CI 干扰的信号结构切片宽度、重排间隔与脉压输出2.1 从截获脉冲到切片序列雷达发射的线性调频脉冲在基带可以写成 s(t)rect(t/T)·exp(jπKt²)其中 T 是脉宽KB/T 是调频斜率。干扰机截获信号后利用 DRFM 得到基带复制品这一步通常还会叠加上一个与距离对应的延迟。CI 的第一步是把这段截获脉冲按固定切片宽度 τ 切成 Nround(T/τ) 段每段保留自己的幅度和相位段内不做调制只改变它在时间轴上的排列位置。切片后的信号模型可以写作 x(t)Σ Aᵢ·s(t−tᵢ)其中 tᵢt₀i·τ 是第 i 片的原始起点。这里的 t₀ 是干扰机转发相对雷达发射时刻的固有延迟。切片动作本身在 MATLAB 里不需要真的做 FFT 之类的频率处理只需按时间索引截取复基带样点。下面这段代码生成一段基带 LFM并演示按索引切片的思路fs 200e6; % 采样率 200 MHz T 10e-6; % 脉宽 10 us B 20e6; % 带宽 20 MHz K B / T; % 调频斜率 t (-T/2 : 1/fs : T/2 - 1/fs).; % 基带时间轴 s exp(1j * pi * K * t.^2); % 复基带 LFM 脉冲 tau 1e-6; % 切片宽度 1 us Nc round(T / tau); % 切片数量 rise (0:Nc-1). * tau; % 每片相对脉冲起点的起始时间 for n 1:Nc idx (t rise(n) - T/2) (t rise(n) - T/2 tau); seg s(idx); % 第 n 片包含幅度和相位 % 后续按重排顺序把 seg 放到新位置即可 end这段代码里K B / T决定了信号的瞬时频率变化率是后面匹配滤波峰值位置计算的基础。tau是切片宽度它直接影响假目标之间的最小间隔切得越窄重排后单位时间内的假目标越多但单片的时宽也越短脉压后的单峰增益会下降。rise保存的是每片在原始脉冲内部的时间偏移重排时真正要移动的正是这些时间偏移量。2.2 重排公式与匹配滤波器输出中的主峰位置CI 的关键操作在“重排”这两个字上。设切片原始编号为 i重排后它被放到编号 p(i) 的位置其中 p 是一个从 0 到 N-1 的排列。再引入重排间隔 d它表示重排后相邻切片中心的时间距离。第 i 片在新序列中的起始时间变成 t₀p(i)·d那么它与原始位置的差值为Δᵢ p(i)·d − i·τ匹配滤波器的输出是干扰信号与发射脉冲共轭反转的卷积。由于每片内部仍然是完整的 LFM 片段脉压后该片会在时延轴上形成一个主峰峰值位置就落在 Δᵢ 对应的时延上。也就是说N 片切片经过重排后理论上有 N 个主峰主峰之间的相对间距由 Δᵢ 的差值决定。实际仿真里这些主峰之间还会出现一些低幅度的交叉项它们来自不同切片之间的互相关幅度通常低于主峰 10 dB 以上不影响假目标群的判读但当切片段不连续、相位突变明显时交叉项会抬高基底表现为脉压输出整体抬升。从频域看切片操作等效于用一组矩形窗去截断信号重排则是对截断产物做时移线性叠加。时域切片会带来频谱旁瓣这在匹配滤波后表现为假目标峰两侧的裙边。所以 CI 不仅仅是在时间轴上“打乱”它的频谱结构已经和原始 LFM 不同这就是干扰机发射信号和雷达发射信号在脉压后不重合的根本原因。2.3 三个关键设计参数与一张参考表做 CI 仿真时最先定下来的应该是切片宽度 τ、重排间隔 d 和重排序列 p。它们分别对应“切多碎”“拉开多远”“按什么顺序放”。先给一组我常用的参数范围后面章节再深入分析参数含义常用范围对脉压输出的影响τ切片宽度T/50 ~ T/10决定假目标间隔下限τ 越小假目标越密d重排间隔τ ~ 2τ控制假目标在时延轴上的散布程度p重排序列伪随机排列决定假目标群的分布形状N切片数量T/τ假目标数量上限N 增加旁瓣基底升高JSR干信比10 ~ 20 dB干扰峰相对噪声底的高度一个常见的误区是把 τ 当作干扰带宽来理解。切片宽度影响的是匹配滤波后主峰的时域位置和宽度而不是干扰信号的频谱范围。干扰信号的带宽在切片瞬间就已经被展宽后面的重排只是重新分配这些频谱分量在时间上的落点。理解这一点就能解释为什么 CI 对 LFM 雷达效果明显LFM 的匹配滤波依赖完整的时频线性关系切片重排刚好破坏了这种关系让原本应该压缩成单峰的能量分散到多个时延上。3. MATLAB 实现切片干扰最小复现脚本与参数表3.1 可以直接运行的基带仿真脚本这一节给出一段完整的基带级 CI 仿真脚本逻辑是按照“生成 LFM → 模拟截获延迟 → 切片 → 随机重排 → 匹配滤波”的顺序往下走。代码里对切片重排部分做了比较清楚的索引处理可以直接复制到 MATLAB 里运行也可以改参数后作为仿真模板。% cinj_demo.m CI 切片重排干扰最小复现脚本 fs 200e6; T 10e-6; B 20e6; K B / T; pad round(5e-6 * fs); % 前后填充避免重排越界 t (-T/2 : 1/fs : T/2 - 1/fs).; s exp(1j * pi * K * t.^2); s_delay [zeros(round(2e-6*fs), 1); s]; % 模拟 2 us 截获延迟 Nbuf length(s_delay) 2 * pad; buf zeros(Nbuf, 1); buf(pad1 : padlength(s_delay)) s_delay; tau 1e-6; % 切片宽度 1 us d_hop 1.2e-6; % 重排间隔 1.2 us Nc round(T / tau); perm randperm(Nc); % 伪随机重排序列 jsr_db 15; % 干信比 15 dB jsr_amp 10^(jsr_db / 20); seg_len round(tau * fs); x_int zeros(Nbuf, 1); for n 1:Nc % 提取第 n 片切片 seg_start pad round( (n-1) * tau * fs ); seg buf(seg_start : seg_start seg_len - 1); % 按重排序号放到新位置 out_start pad round( (perm(n)-1) * d_hop * fs ); x_int(out_start : out_start seg_len - 1) ... x_int(out_start : out_start seg_len - 1) jsr_amp * seg; end % 匹配滤波 h conj(s(end:-1:1)); y conv(x_int, h, same); y y(pad1 : end-pad); % 截掉填充区 t_us (0:length(y)-1) / fs * 1e6; figure; plot(t_us, 20*log10(abs(y) eps)); xlabel(时间 / us); ylabel(相对幅度 / dB); grid on;代码里的seg_start是根据切片编号 n 和切片宽度 τ 计算的原始位置out_start则是重排后的新位置。用randperm(Nc)生成随机序列每运行一次得到不同的假目标分布实际工程里可以固定随机种子来对比不同重排序列的干扰效果。jsr_amp控制干扰信号的幅度改成 1 就是 0 dB 干信比匹配滤波后干扰峰和噪声底的关系会随之变化。运行后会在时间轴上看到一串间隔不等的峰峰数量接近 Nc10主峰之间的间距与d_hop和 τ 的差值有关。如果d_hop远大于 τ假目标群会明显拉散峰值个数不变但分布变稀疏。3.2 核心函数与参数说明上面的脚本适合一次跑完但做参数扫描时更建议把切片重排封装成独立函数。下面是我常用的接口输入为原始脉冲、采样率、切片宽度、重排间隔、重排序列和干信比输出为基带干扰信号function x_int cinj_interleave(s, fs, tau, d_hop, perm, jsr_db) % CI 切片干扰基带信号生成 % s : 截获到的单脉冲基带信号列向量 % fs : 采样率 % tau : 切片宽度单位秒 % d_hop : 重排间隔单位秒 % perm : 切片重排顺序1 到 round(T/tau) 的排列 % jsr_db : 干信比单位 dB T length(s) / fs; Nc length(perm); pad round(0.5e-6 * fs); seg_len round(tau * fs); Nbuf length(s) 2 * pad round(Nc * d_hop * fs); buf zeros(Nbuf, 1); buf(pad1 : padlength(s)) s.; x_int zeros(Nbuf, 1); jsr_amp 10^(jsr_db / 20); for n 1:Nc seg_start pad round((n-1) * tau * fs); seg buf(seg_start : seg_start seg_len - 1); out_start pad round((perm(n)-1) * d_hop * fs); x_int(out_start : out_start seg_len - 1) ... x_int(out_start : out_start seg_len - 1) jsr_amp * seg; end % 截取有效范围长度与输入一致 x_int x_int(1:length(s) round(Nc * d_hop * fs)); end函数内部用的是复基带信号实信号需要先做希尔伯特变换或正交下变频再传入。pad的作用是给重排后的信号留出前后裕量避免切片移动到负数索引。返回值不再是单个脉冲的长度需要考虑重排后的时延扩展通常取输入长度加上Nc * d_hop对应的采样点数。这样设计后调用方可以方便地对比不同 τ 和 d 组合的输出而不需要每次都重写截取逻辑。3.3 用图观察干扰效果仿真脚本跑完画图时重点看三件事。一是主峰数量是否等于切片数 Nc如果少于 Nc说明有切片落到了匹配滤波窗口之外需要加大缓冲区的预留长度。二是主峰间隔是否与d_hop - tau接近CI 在理想情况下相邻假目标的时延差主要取决于这两个参数间隔完全随机会让假目标群看起来更像噪声压制。三是基底电平交叉项会把脉压旁瓣抬高观察20*log10(abs(y))里主峰之间的最低电平如果高于噪声底 20 dB说明切片间的相位不连续已经明显影响了干扰质量。4. 参数边界与对抗从 CI 到工程可用的干扰样式4.1 切片宽度与码元对齐相位编码雷达的约束CI 最常用的对象是 LFM 雷达但实际工程中会遇到相位编码雷达比如 Barker 码或 m 序列调相。相位编码信号的匹配滤波依赖码元序列的完整性如果切片宽度 τ 没有对齐码元边界重排后的干扰信号在解码时会产生大量随机相位噪声假目标峰难以形成。解决办法是把 τ 取成码元宽度 Tc 的整数倍即 τM·TcM 为正整数。这样每个切片内部都包含完整的码元序列重排后在新的位置上仍然能产生相关峰。用 MATLAB 验证时可以先生成一段二相编码信号再按照 M·Tc 切片、重排、做匹配滤波。对比 τ 等于 3 倍码元和 3.5 倍码元两种情况前者脉压后主峰清晰后者峰值会下降 6 dB 以上且旁瓣升高。这一点对做干扰效果仿真很重要很多初学仿真的人直接用 LFM 的参数去干扰相位编码雷达得到的假目标群要么杂乱无章要么能量分散到无法检测。4.2 重排间隔的选择稀疏假目标与压制区的折中重排间隔 d 决定假目标在时延轴上的稀疏程度。d 取 τ 时切片重排后正好前后相接假目标群连续而密集适合制造一个看起来像扩展目标的压制区。d 取 2τ 时切片之间留出空白时隙匹配滤波后主峰之间出现明显间隔适合做多个可分辨假目标。实际操作中 d 最好不要小于 τ否则相邻切片会重叠重叠部分的复数叠加会破坏切片内部的相位结构导致主峰幅度抖动。如果目标是欺骗式干扰要求假目标位置精确可控应该用固定 d 加上固定重排序列。如果目标是压制式干扰侧重让假目标群铺满一段距离区间伪随机 d 的抖动会更好用。我一般在参数扫描时把 d 限制在 τ 到 1.5τ 之间先跑一次快速仿真看假目标分布再根据雷达的距离分辨单元调整。4.3 MTI 雷达与干扰重排顺序的耦合动目标显示MTI雷达利用相邻脉冲对消来滤除静止目标。CI 干扰如果每个脉冲都用相同的重排序列和相同的 d假目标群在多普勒维度上表现为固定回波会被 MTI 当作静止目标抑制掉。应对做法是每个脉冲周期切换重排序列但这种切换会带来假目标群的位置抖动在多次积累后降低峰值稳定性。工程上常用的做法是在一段时间内固定重排序列保持假目标群的相干性积累数个脉冲后整体切换到另一组序列。这样 MTI 对消器在切换瞬间会产生残差之后的驻留窗内仍然能看到较稳定的假目标。在 MATLAB 里做 MTI 对抗仿真可以构造一列脉冲重复间隔 PRI每个 PRI 内生成一段 CI 干扰信号然后做两脉冲对消。对比固定 perm 和逐脉冲切换 perm 两种模式的输出观察对消后剩余峰的幅度。固定 perm 的剩余峰会明显衰减逐脉冲切换则会在对消输出端产生宽带抖动这就是 CI 与 MTI 对抗的基本权衡。4.4 一个边界案例表参数组合脉压输出特征适用场景τ 接近 Td≈0单峰加高旁瓣近似转发式干扰距离欺骗制造单个假目标τT/20dτ固定 perm间隔约 τ 的密集假目标群距离门拖引扩展目标模拟τT/50d2τ随机 perm稀疏且随机分布压制式干扰迷惑检测器τM·Tc固定 permdτ相位编码雷达多个相关峰针对相位编码雷达τT/10dτ逐脉冲切换 perm脉压峰数多但驻留间不稳定对抗 MTI兼顾距离欺骗这个表的核心信息是CI 不是一组固定参数而是切片、重排两个自由度组合出来的干扰样式族。调整 τ 可以控制假目标的最小间隔调整 d 控制散布范围调整 perm 控制分布形状三者解耦独立设置。5. 用 MATLAB 验证干扰效果峰位检测、旁瓣代价与驻留稳定性5.1 峰位与假目标计数验证 CI 干扰是否达到预期第一件事是统计脉压输出里的主峰位置和数量。MATLAB 里用findpeaks加高度门限就能完成y_env abs(y); thr max(y_env) * 0.35; % 主峰高度门限按最大峰值的 35% 取 [pks, locs] findpeaks(y_env, MinPeakHeight, thr, MinPeakDistance, 20); peak_pos_us locs / fs * 1e6; % 转成微秒 % 理论峰位原始切片起点加上 Δi 对应的时延偏移 theory_pos (perm. * d_hop - (0:Nc-1). * tau) * 1e6 ... 2e6; % 加上固有延迟 2 usMinPeakDistance按采样点数设置避免同一个主峰的旁瓣被重复计数。把实测peak_pos_us和theory_pos做差如果最大偏差在 1~2 个采样点以内说明重排逻辑正确。偏差大的话基本可以确定是切片对齐或缓冲索引问题优先检查seg_start和out_start的计算。5.2 多 PRF 驻留下的一致性判断实际雷达会切换脉冲重复频率。CI 干扰如果只针对单 PRF 设计切到新 PRF 后假目标群可能明显变形。验证方法是在多个 PRF 下分别跑同一组 CI 参数记录每次的峰位集合再计算峰位集合的 Jaccard 相似度prf_list [1000, 1200, 1500]; % 不同 PRF单位 Hz peak_sets {}; for prf prf_list PRI 1 / prf; % 这里重新生成一个 PRI 长度的回波加入 CI 干扰 % 略去与 3.1 相同的信号生成部分得到 y_env [~, locs] findpeaks(y_env, MinPeakHeight, thr, ... MinPeakDistance, 20); peak_sets{end1} sort(locs); end inter length(intersect(peak_sets{1}, peak_sets{2})); union_len length(union(peak_sets{1}, peak_sets{2})); jaccard inter / union_len;Jaccard 相似度高于 0.8说明干扰样式对不同 PRF 稳定低于 0.5 则说明重排序列与 PRF 耦合严重实战中需要把 perm 序列设计成与 PRI 无关的版本或者采用固定切片宽度但可变重排间隔的策略。5.3 一个进阶技巧用复包络斜率快速检查切片拼接处的相位不连续CI 干扰在切片拼接处经常出现相位跳变跳变幅度过大会抬高脉压旁瓣。检查方法不一定要上频谱分析直接看脉压输出信号的复包络相位差分dphase angle(y(2:end) .* conj(y(1:end-1))); spikes find(abs(dphase) 0.3); % 相位差超过 0.3 rad 的位置spikes出现的位置通常对应切片边界。如果边界数量过多且分布在整个脉冲宽度内说明需要给每段切片做初始相位补偿即在重排前把每片的起点相位对齐到同一参考相位。补偿后的干扰信号脉压旁瓣基底一般能降低 6 dB 以上。反过来如果某些切片在重排后相位连续脉压峰值会异常突出这时可以主动引入随机相位扰动来掩盖这种结构让干扰样式更接近噪声特性。本文还有配套的精品资源点击获取
基于MATLAB的CI切片重排干扰仿真与参数分析
简介适用于雷达电子战与干扰技术研究的MATLAB脚本围绕切片干扰Chopping and Interleaving算法提供完整实现帮助工程师和研究人员直观理解通过切割、重排雷达信号来破坏敌方探测与跟踪的机制。压缩包内仅含1个m文件大小1KB属于轻量级源码便于直接运行查看效果并在此基础上修改参数或扩展功能。目前已有1051人学习下载适合作为高校雷达课程、电子战培训或抗干扰算法设计的教学示例与参考起点。脚本内部涵盖信号生成、切片分割、交错组合、干扰效果评估等关键环节其中信号生成部分支持脉冲序列、频率调制、相位编码等常见雷达波形切片与交织模块可按预设规则或随机方式重组信号片段读者可对照代码梳理雷达干扰的典型流程也可将其嵌入SystemVue等系统仿真环境进行硬件在环验证进一步评估误检率、漏检率等干扰性能指标。1. ChoppingAndInterleaving一种面向线性调频雷达的切片重排干扰在工程仿真里复制一个线性调频脉冲然后整体转发脉压匹配滤波后在固定延迟处只会出现一个又高又窄的单峰稍微有点经验的雷达操作员都能把它从真实回波里区分出来。ChoppingAndInterleaving后文统一写 CI换了一种思路把截获到的雷达脉冲在时间轴上切成若干段再按新的顺序拼接发射。切片宽度和重排间隔设计得当的话匹配滤波输出端会出现一串幅度相近、位置可控的假目标并且这串假目标还能在多次积累后保持稳定干扰效果比单峰转发更贴合工程需要。下面直接从信号结构拆起用 MATLAB 从基带脉冲开始把 CI 的完整链路写一遍给出一套能改参数复现的脚本再讨论切片宽度、重排间隔和相位连续性问题。适合正在做雷达信号级仿真、DRFM 干扰建模或电子对抗效果评估的工程师参考。2. CI 干扰的信号结构切片宽度、重排间隔与脉压输出2.1 从截获脉冲到切片序列雷达发射的线性调频脉冲在基带可以写成 s(t)rect(t/T)·exp(jπKt²)其中 T 是脉宽KB/T 是调频斜率。干扰机截获信号后利用 DRFM 得到基带复制品这一步通常还会叠加上一个与距离对应的延迟。CI 的第一步是把这段截获脉冲按固定切片宽度 τ 切成 Nround(T/τ) 段每段保留自己的幅度和相位段内不做调制只改变它在时间轴上的排列位置。切片后的信号模型可以写作 x(t)Σ Aᵢ·s(t−tᵢ)其中 tᵢt₀i·τ 是第 i 片的原始起点。这里的 t₀ 是干扰机转发相对雷达发射时刻的固有延迟。切片动作本身在 MATLAB 里不需要真的做 FFT 之类的频率处理只需按时间索引截取复基带样点。下面这段代码生成一段基带 LFM并演示按索引切片的思路fs 200e6; % 采样率 200 MHz T 10e-6; % 脉宽 10 us B 20e6; % 带宽 20 MHz K B / T; % 调频斜率 t (-T/2 : 1/fs : T/2 - 1/fs).; % 基带时间轴 s exp(1j * pi * K * t.^2); % 复基带 LFM 脉冲 tau 1e-6; % 切片宽度 1 us Nc round(T / tau); % 切片数量 rise (0:Nc-1). * tau; % 每片相对脉冲起点的起始时间 for n 1:Nc idx (t rise(n) - T/2) (t rise(n) - T/2 tau); seg s(idx); % 第 n 片包含幅度和相位 % 后续按重排顺序把 seg 放到新位置即可 end这段代码里K B / T决定了信号的瞬时频率变化率是后面匹配滤波峰值位置计算的基础。tau是切片宽度它直接影响假目标之间的最小间隔切得越窄重排后单位时间内的假目标越多但单片的时宽也越短脉压后的单峰增益会下降。rise保存的是每片在原始脉冲内部的时间偏移重排时真正要移动的正是这些时间偏移量。2.2 重排公式与匹配滤波器输出中的主峰位置CI 的关键操作在“重排”这两个字上。设切片原始编号为 i重排后它被放到编号 p(i) 的位置其中 p 是一个从 0 到 N-1 的排列。再引入重排间隔 d它表示重排后相邻切片中心的时间距离。第 i 片在新序列中的起始时间变成 t₀p(i)·d那么它与原始位置的差值为Δᵢ p(i)·d − i·τ匹配滤波器的输出是干扰信号与发射脉冲共轭反转的卷积。由于每片内部仍然是完整的 LFM 片段脉压后该片会在时延轴上形成一个主峰峰值位置就落在 Δᵢ 对应的时延上。也就是说N 片切片经过重排后理论上有 N 个主峰主峰之间的相对间距由 Δᵢ 的差值决定。实际仿真里这些主峰之间还会出现一些低幅度的交叉项它们来自不同切片之间的互相关幅度通常低于主峰 10 dB 以上不影响假目标群的判读但当切片段不连续、相位突变明显时交叉项会抬高基底表现为脉压输出整体抬升。从频域看切片操作等效于用一组矩形窗去截断信号重排则是对截断产物做时移线性叠加。时域切片会带来频谱旁瓣这在匹配滤波后表现为假目标峰两侧的裙边。所以 CI 不仅仅是在时间轴上“打乱”它的频谱结构已经和原始 LFM 不同这就是干扰机发射信号和雷达发射信号在脉压后不重合的根本原因。2.3 三个关键设计参数与一张参考表做 CI 仿真时最先定下来的应该是切片宽度 τ、重排间隔 d 和重排序列 p。它们分别对应“切多碎”“拉开多远”“按什么顺序放”。先给一组我常用的参数范围后面章节再深入分析参数含义常用范围对脉压输出的影响τ切片宽度T/50 ~ T/10决定假目标间隔下限τ 越小假目标越密d重排间隔τ ~ 2τ控制假目标在时延轴上的散布程度p重排序列伪随机排列决定假目标群的分布形状N切片数量T/τ假目标数量上限N 增加旁瓣基底升高JSR干信比10 ~ 20 dB干扰峰相对噪声底的高度一个常见的误区是把 τ 当作干扰带宽来理解。切片宽度影响的是匹配滤波后主峰的时域位置和宽度而不是干扰信号的频谱范围。干扰信号的带宽在切片瞬间就已经被展宽后面的重排只是重新分配这些频谱分量在时间上的落点。理解这一点就能解释为什么 CI 对 LFM 雷达效果明显LFM 的匹配滤波依赖完整的时频线性关系切片重排刚好破坏了这种关系让原本应该压缩成单峰的能量分散到多个时延上。3. MATLAB 实现切片干扰最小复现脚本与参数表3.1 可以直接运行的基带仿真脚本这一节给出一段完整的基带级 CI 仿真脚本逻辑是按照“生成 LFM → 模拟截获延迟 → 切片 → 随机重排 → 匹配滤波”的顺序往下走。代码里对切片重排部分做了比较清楚的索引处理可以直接复制到 MATLAB 里运行也可以改参数后作为仿真模板。% cinj_demo.m CI 切片重排干扰最小复现脚本 fs 200e6; T 10e-6; B 20e6; K B / T; pad round(5e-6 * fs); % 前后填充避免重排越界 t (-T/2 : 1/fs : T/2 - 1/fs).; s exp(1j * pi * K * t.^2); s_delay [zeros(round(2e-6*fs), 1); s]; % 模拟 2 us 截获延迟 Nbuf length(s_delay) 2 * pad; buf zeros(Nbuf, 1); buf(pad1 : padlength(s_delay)) s_delay; tau 1e-6; % 切片宽度 1 us d_hop 1.2e-6; % 重排间隔 1.2 us Nc round(T / tau); perm randperm(Nc); % 伪随机重排序列 jsr_db 15; % 干信比 15 dB jsr_amp 10^(jsr_db / 20); seg_len round(tau * fs); x_int zeros(Nbuf, 1); for n 1:Nc % 提取第 n 片切片 seg_start pad round( (n-1) * tau * fs ); seg buf(seg_start : seg_start seg_len - 1); % 按重排序号放到新位置 out_start pad round( (perm(n)-1) * d_hop * fs ); x_int(out_start : out_start seg_len - 1) ... x_int(out_start : out_start seg_len - 1) jsr_amp * seg; end % 匹配滤波 h conj(s(end:-1:1)); y conv(x_int, h, same); y y(pad1 : end-pad); % 截掉填充区 t_us (0:length(y)-1) / fs * 1e6; figure; plot(t_us, 20*log10(abs(y) eps)); xlabel(时间 / us); ylabel(相对幅度 / dB); grid on;代码里的seg_start是根据切片编号 n 和切片宽度 τ 计算的原始位置out_start则是重排后的新位置。用randperm(Nc)生成随机序列每运行一次得到不同的假目标分布实际工程里可以固定随机种子来对比不同重排序列的干扰效果。jsr_amp控制干扰信号的幅度改成 1 就是 0 dB 干信比匹配滤波后干扰峰和噪声底的关系会随之变化。运行后会在时间轴上看到一串间隔不等的峰峰数量接近 Nc10主峰之间的间距与d_hop和 τ 的差值有关。如果d_hop远大于 τ假目标群会明显拉散峰值个数不变但分布变稀疏。3.2 核心函数与参数说明上面的脚本适合一次跑完但做参数扫描时更建议把切片重排封装成独立函数。下面是我常用的接口输入为原始脉冲、采样率、切片宽度、重排间隔、重排序列和干信比输出为基带干扰信号function x_int cinj_interleave(s, fs, tau, d_hop, perm, jsr_db) % CI 切片干扰基带信号生成 % s : 截获到的单脉冲基带信号列向量 % fs : 采样率 % tau : 切片宽度单位秒 % d_hop : 重排间隔单位秒 % perm : 切片重排顺序1 到 round(T/tau) 的排列 % jsr_db : 干信比单位 dB T length(s) / fs; Nc length(perm); pad round(0.5e-6 * fs); seg_len round(tau * fs); Nbuf length(s) 2 * pad round(Nc * d_hop * fs); buf zeros(Nbuf, 1); buf(pad1 : padlength(s)) s.; x_int zeros(Nbuf, 1); jsr_amp 10^(jsr_db / 20); for n 1:Nc seg_start pad round((n-1) * tau * fs); seg buf(seg_start : seg_start seg_len - 1); out_start pad round((perm(n)-1) * d_hop * fs); x_int(out_start : out_start seg_len - 1) ... x_int(out_start : out_start seg_len - 1) jsr_amp * seg; end % 截取有效范围长度与输入一致 x_int x_int(1:length(s) round(Nc * d_hop * fs)); end函数内部用的是复基带信号实信号需要先做希尔伯特变换或正交下变频再传入。pad的作用是给重排后的信号留出前后裕量避免切片移动到负数索引。返回值不再是单个脉冲的长度需要考虑重排后的时延扩展通常取输入长度加上Nc * d_hop对应的采样点数。这样设计后调用方可以方便地对比不同 τ 和 d 组合的输出而不需要每次都重写截取逻辑。3.3 用图观察干扰效果仿真脚本跑完画图时重点看三件事。一是主峰数量是否等于切片数 Nc如果少于 Nc说明有切片落到了匹配滤波窗口之外需要加大缓冲区的预留长度。二是主峰间隔是否与d_hop - tau接近CI 在理想情况下相邻假目标的时延差主要取决于这两个参数间隔完全随机会让假目标群看起来更像噪声压制。三是基底电平交叉项会把脉压旁瓣抬高观察20*log10(abs(y))里主峰之间的最低电平如果高于噪声底 20 dB说明切片间的相位不连续已经明显影响了干扰质量。4. 参数边界与对抗从 CI 到工程可用的干扰样式4.1 切片宽度与码元对齐相位编码雷达的约束CI 最常用的对象是 LFM 雷达但实际工程中会遇到相位编码雷达比如 Barker 码或 m 序列调相。相位编码信号的匹配滤波依赖码元序列的完整性如果切片宽度 τ 没有对齐码元边界重排后的干扰信号在解码时会产生大量随机相位噪声假目标峰难以形成。解决办法是把 τ 取成码元宽度 Tc 的整数倍即 τM·TcM 为正整数。这样每个切片内部都包含完整的码元序列重排后在新的位置上仍然能产生相关峰。用 MATLAB 验证时可以先生成一段二相编码信号再按照 M·Tc 切片、重排、做匹配滤波。对比 τ 等于 3 倍码元和 3.5 倍码元两种情况前者脉压后主峰清晰后者峰值会下降 6 dB 以上且旁瓣升高。这一点对做干扰效果仿真很重要很多初学仿真的人直接用 LFM 的参数去干扰相位编码雷达得到的假目标群要么杂乱无章要么能量分散到无法检测。4.2 重排间隔的选择稀疏假目标与压制区的折中重排间隔 d 决定假目标在时延轴上的稀疏程度。d 取 τ 时切片重排后正好前后相接假目标群连续而密集适合制造一个看起来像扩展目标的压制区。d 取 2τ 时切片之间留出空白时隙匹配滤波后主峰之间出现明显间隔适合做多个可分辨假目标。实际操作中 d 最好不要小于 τ否则相邻切片会重叠重叠部分的复数叠加会破坏切片内部的相位结构导致主峰幅度抖动。如果目标是欺骗式干扰要求假目标位置精确可控应该用固定 d 加上固定重排序列。如果目标是压制式干扰侧重让假目标群铺满一段距离区间伪随机 d 的抖动会更好用。我一般在参数扫描时把 d 限制在 τ 到 1.5τ 之间先跑一次快速仿真看假目标分布再根据雷达的距离分辨单元调整。4.3 MTI 雷达与干扰重排顺序的耦合动目标显示MTI雷达利用相邻脉冲对消来滤除静止目标。CI 干扰如果每个脉冲都用相同的重排序列和相同的 d假目标群在多普勒维度上表现为固定回波会被 MTI 当作静止目标抑制掉。应对做法是每个脉冲周期切换重排序列但这种切换会带来假目标群的位置抖动在多次积累后降低峰值稳定性。工程上常用的做法是在一段时间内固定重排序列保持假目标群的相干性积累数个脉冲后整体切换到另一组序列。这样 MTI 对消器在切换瞬间会产生残差之后的驻留窗内仍然能看到较稳定的假目标。在 MATLAB 里做 MTI 对抗仿真可以构造一列脉冲重复间隔 PRI每个 PRI 内生成一段 CI 干扰信号然后做两脉冲对消。对比固定 perm 和逐脉冲切换 perm 两种模式的输出观察对消后剩余峰的幅度。固定 perm 的剩余峰会明显衰减逐脉冲切换则会在对消输出端产生宽带抖动这就是 CI 与 MTI 对抗的基本权衡。4.4 一个边界案例表参数组合脉压输出特征适用场景τ 接近 Td≈0单峰加高旁瓣近似转发式干扰距离欺骗制造单个假目标τT/20dτ固定 perm间隔约 τ 的密集假目标群距离门拖引扩展目标模拟τT/50d2τ随机 perm稀疏且随机分布压制式干扰迷惑检测器τM·Tc固定 permdτ相位编码雷达多个相关峰针对相位编码雷达τT/10dτ逐脉冲切换 perm脉压峰数多但驻留间不稳定对抗 MTI兼顾距离欺骗这个表的核心信息是CI 不是一组固定参数而是切片、重排两个自由度组合出来的干扰样式族。调整 τ 可以控制假目标的最小间隔调整 d 控制散布范围调整 perm 控制分布形状三者解耦独立设置。5. 用 MATLAB 验证干扰效果峰位检测、旁瓣代价与驻留稳定性5.1 峰位与假目标计数验证 CI 干扰是否达到预期第一件事是统计脉压输出里的主峰位置和数量。MATLAB 里用findpeaks加高度门限就能完成y_env abs(y); thr max(y_env) * 0.35; % 主峰高度门限按最大峰值的 35% 取 [pks, locs] findpeaks(y_env, MinPeakHeight, thr, MinPeakDistance, 20); peak_pos_us locs / fs * 1e6; % 转成微秒 % 理论峰位原始切片起点加上 Δi 对应的时延偏移 theory_pos (perm. * d_hop - (0:Nc-1). * tau) * 1e6 ... 2e6; % 加上固有延迟 2 usMinPeakDistance按采样点数设置避免同一个主峰的旁瓣被重复计数。把实测peak_pos_us和theory_pos做差如果最大偏差在 1~2 个采样点以内说明重排逻辑正确。偏差大的话基本可以确定是切片对齐或缓冲索引问题优先检查seg_start和out_start的计算。5.2 多 PRF 驻留下的一致性判断实际雷达会切换脉冲重复频率。CI 干扰如果只针对单 PRF 设计切到新 PRF 后假目标群可能明显变形。验证方法是在多个 PRF 下分别跑同一组 CI 参数记录每次的峰位集合再计算峰位集合的 Jaccard 相似度prf_list [1000, 1200, 1500]; % 不同 PRF单位 Hz peak_sets {}; for prf prf_list PRI 1 / prf; % 这里重新生成一个 PRI 长度的回波加入 CI 干扰 % 略去与 3.1 相同的信号生成部分得到 y_env [~, locs] findpeaks(y_env, MinPeakHeight, thr, ... MinPeakDistance, 20); peak_sets{end1} sort(locs); end inter length(intersect(peak_sets{1}, peak_sets{2})); union_len length(union(peak_sets{1}, peak_sets{2})); jaccard inter / union_len;Jaccard 相似度高于 0.8说明干扰样式对不同 PRF 稳定低于 0.5 则说明重排序列与 PRF 耦合严重实战中需要把 perm 序列设计成与 PRI 无关的版本或者采用固定切片宽度但可变重排间隔的策略。5.3 一个进阶技巧用复包络斜率快速检查切片拼接处的相位不连续CI 干扰在切片拼接处经常出现相位跳变跳变幅度过大会抬高脉压旁瓣。检查方法不一定要上频谱分析直接看脉压输出信号的复包络相位差分dphase angle(y(2:end) .* conj(y(1:end-1))); spikes find(abs(dphase) 0.3); % 相位差超过 0.3 rad 的位置spikes出现的位置通常对应切片边界。如果边界数量过多且分布在整个脉冲宽度内说明需要给每段切片做初始相位补偿即在重排前把每片的起点相位对齐到同一参考相位。补偿后的干扰信号脉压旁瓣基底一般能降低 6 dB 以上。反过来如果某些切片在重排后相位连续脉压峰值会异常突出这时可以主动引入随机相位扰动来掩盖这种结构让干扰样式更接近噪声特性。本文还有配套的精品资源点击获取