简介面向合成孔径雷达SAR电子对抗领域的欺骗干扰功能MATLAB仿真代码贴合本科、硕士阶段物理应用方向教学科研使用。程序基于MATLAB 2019a设计涵盖欺骗干扰信号构造、回波模拟、SAR成像处理等关键环节并加入压制干扰对比程序便于从原理到实现系统理解干扰机制。压缩包共57个文件其中32个.m脚本与函数、20个.asv自动备份、2个Word运行说明、2个.mat数据和1幅bmp测试图整体大小约5.78MB代码与文档组织清晰。已有413人学习下载配合运行说明可快速复现典型欺骗干扰场景适合课程设计、毕业设计或雷达对抗课题的仿真参考。1. 为什么说SAR欺骗干扰是雷达对抗里性价比最高的方案SAR成像侦察现在几乎是星载和机载平台的标配图像分辨率能做到亚米级想在强侦察下隐藏真实部署已经很难了。压制干扰的思路是拿噪声把整个场景盖掉等效功率要高出好几个数量级才能压得住匹配处理增益而欺骗干扰只用一个很小的转发功率就能在最终的SAR图像里“造”出根本不存在的虚假目标和编队。这个资源包里的qipian.m做的正是这件事截获SAR信号参数通过延时、移频和相位调制把假目标直接注入到成像结果里。对我这种做物理仿真和雷达信号处理的人来说它把从信号生成、回波模拟到干扰效果验证的整条链路都串起来了适合本科、硕士阶段做SAR对抗课题时当脚手架也适合刚接触电子战仿真的工程师快速理解欺骗干扰的作用边界。2. 欺骗干扰的成像域机理延时、移频如何变成图上的假目标2.1 距离向与方位向的“作弊”入口SAR图像的距离向分辨率靠的是大带宽线性调频LFM信号做脉冲压缩方位向分辨率靠的是合成孔径期间的多普勒历史。欺骗干扰要在这两个维度分别做文章。先看距离向。雷达发射的LFM信号可以写成% LFM基带信号K为调频斜率 t -Tp/2:1/Fs:Tp/2; s_ref exp(1j*pi*K*t.^2); % 匹配滤波参考信号匹配滤波的本质是把线性调频信号的能量压成一个窄峰峰的位置对应回波延时。干扰机如果对截获到的信号整体延时τ再转发回去SAR端匹配滤波后峰的位置就会移动c*τ/2也就等效于在距离向上制造了一个偏移目标。更妙的是因为LFM信号的“延时-频率”耦合关系对信号做固定移频f_d等效于产生一个附加延时f_d/K同样能把假目标“推”到别的距离单元上去。这就是欺骗干扰在距离向的两种基本注入手段显式延时和移频。方位向的欺骗要稍微绕一点。SAR平台飞行时真实目标与雷达之间的斜距先缩短再拉长形成多普勒历史方位压缩就是利用这个历史聚焦的。干扰机逐脉冲地给转发信号叠加一个可控的相位增量相当于伪造了一个不同的多普勒历史SAR端方位压缩后就会在某个方位位置形成虚假目标。所以一套完整的欺骗干扰仿真必须同时具备两个自由度脉冲内调制延时/移频和脉冲间调制相位编码。2.2 欺骗干扰的功率优势匹配处理增益就是你的杠杆欺骗干扰的另一个关键参数是干信比。SAR接收端对信号做匹配滤波时宽带宽积B*Tp会带来等效的处理增益。压制干扰要有效噪声功率得淹没有用信号于是干扰机功率必须比真实回波大很多而欺骗干扰只要保证截获信号后转发的功率能压过接收机噪声底匹配输出后的假目标幅度就会和真目标相当。一边是面临B*Tp量级的处理增益劣势一边是享受同样的增益两边差的功率就是干扰机的“杠杆”。从工程上说典型星载SAR的时宽带宽积动辄上千这意味着压制干扰要付出30dB以上的额外功率代价。下面这组对比在实际评估时很有参考价值。维度压制干扰欺骗干扰作用效果图像噪化、目标淹没图像出现虚假目标/编队所需功率高需对抗匹配增益低利用匹配增益信号知识依赖只需大致载频和带宽需要精确的LFM参数与PRF对抗检测难度易被识别为有源干扰假目标融入场景难以判别主要失效条件高增益天线和窄带滤波参数测量误差和算法失配2.3 欺骗干扰的“算力”瓶颈在哪里欺骗干扰在工程上最大的难点不是发射功率而是对SAR信号的实时参数测量。你要在极短的脉冲宽度内估计出载频、调频斜率、脉冲宽度和脉冲重复频率然后马上决定延时量、移频量和相位编码序列再转发出去。这个包里gonglvpu.m就是用来做功率谱分析和参数粗测的MAICHONG.m负责生成典型的SAR脉冲模板。真正到工程实现这些模块通常要跑在本振同步和FPGA高速信号链上MATLAB这份代码的价值在于先把算法链路验证通把误差敏感度摸清楚再往硬件移植就快得多。跑这套代码时我一般会先不加欺骗干扰只用回波模拟与成像模块看一个点目标的聚焦效果确认SAR成像链路本身干净再去叠加欺骗干扰模块看假目标是否按要求出现。这个顺序能省掉大量联调时间。3. MATLAB模块逐段拆解脉冲生成、功率谱与欺骗参数求解3.1 从MAICHONG.m看SAR脉冲模板设计SAR发射信号通常是宽带LFM脉冲脉冲宽度Tp决定平均功率带宽B决定距离向分辨率。MAICHONG.m本质上就是按照给定的B、Tp、采样率Fs和载频fc构造一个复数基带或中频脉冲。下面这段是从原工程思路里抽出来的教学版保持功能一致变量名直接可见% MAICHONG.m 教学版构造SAR发射脉冲模板 function s MAICHONG(B, Tp, Fs, fc) % B 信号带宽单位Hz % Tp 脉冲宽度单位s % Fs 采样率单位Hz % fc 载频单位Hz用于后续中频搬移基带仿真可置0 K B / Tp; % 调频斜率 N round(Tp * Fs); % 单脉冲采样点数 t (0:N-1) / Fs - Tp/2; % 以脉冲中心为零点的时间轴 s_b exp(1j * pi * K * t.^2); % 基带LFM信号 s s_b .* exp(1j * 2 * pi * fc * t); % 搬移到载频便于后续模拟下变频 end调用时比如B 30e6; Tp 10e-6; Fs 60e6; s MAICHONG(B, Tp, Fs, 0);就生成一个30MHz带宽、10微秒脉宽的LFM基带脉冲。注意Fs 60e6是带宽的两倍这是为了留出过采样余量避免匹配滤波后峰值旁瓣在采样点上失真。如果Fs只设成和带宽相同脉冲压缩后会看到峰值有细微的栅瓣效应幅度误差不大但做定量分析时会影响旁瓣电平读数。3.2 功率谱分析gonglvpu.m的截获参数测量思路干扰机截获到SAR信号后第一步是测带宽和中心频率。gonglvpu.m用FFT求功率谱从谱宽估计带宽。工程上更稳的做法是CZT线性调频Z变换做局部细化不过对仿真级验证来说FFT加矩形窗或汉明窗已经够用。% gonglvpu.m 教学版功率谱估计与带宽测量 function [B_est, fc_est] gonglvpu(x, Fs) N length(x); w hamming(N).; % 加窗抑制频谱泄漏 X fftshift(fft(x .* w, N)); f_axis (-N/2:N/2-1) * Fs / N; P abs(X).^2 / (N * Fs); % 功率谱密度 % 找到峰值附近功率比峰值低3dB的两个频点 [Pmax, idx] max(P); th 10^(0.3) * Pmax; % -3dB门限计算要注意单位 idx_range find(P Pmax / 2); % 半功率点 B_est f_axis(max(idx_range)) - f_axis(min(idx_range)); fc_est (f_axis(max(idx_range)) f_axis(min(idx_range))) / 2; end半功率点的选择直接对应3dB带宽这是衡量LFM信号带宽的常用口径。要注意如果输入信号里混有强干扰或多次谐波FFT谱会变宽估计出来的带宽会偏大。这时候可以先用带通滤波器把感兴趣的频段框出来再做谱分析。gonglvpu.m后接的juchi.m则是从距离维上计算目标回波的延时为欺骗干扰的延时量设定提供基准。3.3 欺骗参数求解把假目标坐标翻译成延时和移频假目标的位置由你在距离向和方位向的设定值决定。假设场景中存在一个真实点目标你要在它旁边ΔR米处放一个假目标那么对应的距离向延时就是τ 2*ΔR/c。如果要通过移频实现同样的偏移则需要的移频量为f_d K*τ。从工程上讲显式延时的好处是直观、对LFM的调频斜率不敏感缺点是需要大的延时线移频的好处是硬件实现简单一个混频器就搞定但前提是你对K的估计要准。TIAOPIN.m和TIAOXIANG.m分别对应调频和调相两种干扰调制样式适合用来对比同一假目标在不同调制方式下的成像差异。下面给出一组常用的参数设定参考参数典型取值作用带宽 B30 MHz决定距离向分辨率脉宽 Tp10 μs决定作用距离与平均功率调频斜率 K3×10^12 Hz/s决定移频-距离换算比例PRF1000 Hz决定方位向无模糊范围移频量 f_d0.55 MHz控制假目标距离偏移量相位编码步进02π控制假目标方位向位置参数表的逻辑链条是这样的你先定距离向偏移算出延时和移频再定方位向偏移算出每个脉冲上该叠的相位增量最后把这些参数输入干扰调制模块。我在调试时习惯先把f_d设成B/2以内的值超过这个范围会引起距离向频谱混叠假目标会出现成对虚影。4. 全链路复现回波模拟、欺骗叠加与SAR成像验证4.1 回波模拟huibo.m与huiboshuju.mat回波模拟是欺骗干扰的靶场。huibo.m负责生成不含干扰的真实回波huiboshuju.mat则是把回波数据存下来的结果。做干扰仿真前先有干净回波才能在成像图上分清哪个是真实目标、哪个是假目标。回波生成的核心是计算每个方位采样时刻目标到雷达的斜距并把发射脉冲按延时叠加进接收向量。% huibo.m 教学版生成单个点目标的SAR回波 function echo huibo(s_ref, R0, V, X_target, Tp, Fs, PRF, Na) % R0 最近斜距(m) % V 平台速度(m/s) % X_target 目标方位向坐标(m) % Na 方位向脉冲数 c 3e8; K 30e6 / Tp; % 从带宽和脉宽换算调频斜率 t_fast -Tp/2:1/Fs:Tp/2; % 距离向快时间 echo zeros(Na, length(t_fast)); for ia 1:Na t_az (ia - Na/2) / PRF; % 方位向慢时间 R sqrt(R0^2 (V*t_az - X_target)^2); % 瞬时斜距 tau 2*R/c; % 双程延时 idx round(tau*Fs) 1; if idx length(t_fast) echo(ia, :) s_ref .* exp(-1j*2*pi*fc*tau); % 简化的接收信号 end end end这段代码里的R用勾股定理计算没有做抛物线近似好处是在正侧视大斜视角时依然能保持较高的相位精度。慢时间轴以Na/2为中心对应合成孔径中心时刻。实际工程中为了节省内存不会把整个回波矩阵都存下来而是按方位块处理huiboshuju.mat里存的其实是一帧已经截断、去斜后的回波数据块。4.2 欺骗信号生成与叠加qipian.m的主流程qipian.m是整个工程的发动机。它从截获信号里提取参数生成调制后的欺骗脉冲再和真实回波相加。下面这个版本是经过我简化重写的教学主干完整工程里还挂了.asv自动保存的中间版本跑的时候直接认run.m即可% qipian.m 教学版欺骗干扰主流程 % 假设已经截获并测得了信号参数fs, B, Tp, fc, PRF s_jam_all zeros(size(echo)); % 预分配干扰矩阵 for ia 1:Na % 距离向移频 延时 t -Tp/2:1/Fs:Tp/2; s_jam exp(1j*pi*K*t.^2); % 复制真实脉冲 s_jam s_jam .* exp(1j*2*pi*fd*t); % 移频fd实现距离向偏移 tau_deceive 2*R_target/c; % 假目标设定的距离延时 n_offset round(tau_deceive * Fs); s_jam circshift(s_jam, n_offset, 2); % 循环移位模拟延时 % 方位向叠加相位步进伪造多普勒历史 phase_az exp(1j * 2*pi * ia * fdaz / PRF); s_jam_all(ia, :) phase_az * s_jam; % 干信比控制按所需JSR放大干扰幅度 s_jam_all(ia, :) s_jam_all(ia, :) * sqrt(10^(JSR/10)); end echo_jam echo s_jam_all; % 真实回波上面叠假目标这里circshift是循环移位注意它会把信号尾部卷到头部来。对于远距离目标延时远大于脉宽的场景这种卷绕会造成距离模糊假象正确做法是开一个长度Na×Ns的大矩阵用n_offset直接索引写入。上面代码为了保证可读性用了循环移位实际工程我改成线性写入。JSR干信比一般设成-10到10dB之间就能看到明显的欺骗效果太高会连假目标旁瓣一起变亮反而暴露干扰痕迹。4.3 成像验证距离压缩、距离徙动校正与方位压缩欺骗干扰成不成功最终要在SAR图像上看。最基本的验证链路是先距离压缩再做距离徙动校正RCMC最后方位压缩。下面给出一个不含插值的简化RD算法流程。% 距离压缩匹配滤波 echo_rc ifft(fft(echo_jam, Nf, 2) .* conj(fft(s_ref, Nf, 2)), Nf, 2); echo_rc echo_rc(:, 1:length(t_fast)); % 距离徙动校正按最近斜距计算徙动曲线逐距离门搬移简化最近邻 for ia 1:Na R_true sqrt(R0^2 (V*(ia-Na/2)/PRF)^2); dr 2*(R_true - R0)/c * Fs; % 需要校正的采样点偏移 shift round(dr); if shift ~ 0 echo_rc(ia, :) circshift(echo_rc(ia, :), -shift, 2); end end % 方位压缩对每个距离门做方位向FFT img fftshift(fft(echo_rc, Na, 1), 1); img abs(img).^2;这段代码的用意是把回波先压缩成距离-慢时间矩阵再沿方位向压缩得到聚焦图像。干扰是否生效看两个特征第一假目标是否出现在你设定的距离和方位坐标上第二假目标的聚焦质量是否和真目标一致。如果假目标出现但明显散焦基本可以判断是移频量或相位编码步进设错了如果假目标出现在镜像位置则是方位向频谱混叠需要把PRF调大或降低相位编码变化速率。这个包的lena.bmp不是拿来当目标图像的而是用来模拟SAR场景后向散射系数分布把图像灰度值映射成地面各点的反射率再逐点生成回波。这样你能直观地看到欺骗干扰在一个真实纹理场景里的视觉效果比单纯点目标有说服力得多。5. 参数边界与工程排错五个典型坑与新版MATLAB迁移5.1 假目标位置偏移的定量校验调试时最常问的问题是“我设了ΔR和ΔX为什么图上假目标不在预期位置”。原因多半是单位换算里漏了双程因子。距离向延时换算用τ 2ΔR/c移频换算用ΔR c*f_d/(2K)。第二项要特别注意K的正负上扫频和下扫频的LFM符号不同移频产生的距离偏移方向正好相反。在qipian.m里改干扰前先打印一组换算结果和成像图上读出的峰值位置对一下通常十分钟内能定位问题。方位向的偏移量则按ΔX V_fake * t_az算这里的V_fake是你通过相位步进伪装的等效速度不是平台速度。5.2 多普勒模糊与成对回波当干扰相位步进对应的多普勒频率超过PRF/2时就会产生方位向混叠假目标会出现在真实的镜像位置形成成对目标。解决思路是限制最大相位变化率或者提高PRF重频。工程上另一个常见现象是假目标在方位向“拉丝”这通常是因为相位编码序列存在固定偏差导致多普勒中心偏移未完全补偿试试在方位压缩前对慢时间信号做去斜处理。5.3.asv文件和新版MATLAB的兼容性包里那些.asv是MATLAB自动保存的草稿不是必须的运行文件。拿到工程后先看运行说明.doc再直接运行run.m。如果用的不是R2019a而是新版比如R2023b或R2026b需要注意新版对circshift和矩阵索引规则没有重大改动但fftshift的默认维度行为在R2021a后更严格了代码里最好明确写出fftshift(X, 1)这样的维度参数否则二维矩阵时可能得到和旧版不一致的显示结果。另外新版MATLAB的plot配色和坐标轴属性默认值变了不影响数据计算但脚本里若有依赖旧版坐标轴颜色的后续叠加绘图要检查一遍。5.4 与压制干扰的对比判断信号“是否被识破”从接收端视角回看这套代码可以额外做个实验对同一份真实SAR数据分别加噪声压制干扰和欺骗干扰再用CFAR检测统计虚警率和漏警率。欺骗干扰产生的假目标如果功率调得和真实目标接近CFAR的检测门限完全无法区分压制干扰则直接把背景噪声抬高CFAR门限上移真实目标反而可能被淹掉并触发全局告警——从隐蔽对抗的角度讲欺骗干扰明显更难被识破这也是它被频繁用于SAR对抗仿真教研的原因。如果你想在这个包里继续深入可以试着把净输出改成多个假目标沿同一方向同步移动观察SAR成像结果会不会出现“幽灵编队”效应。本文还有配套的精品资源点击获取 SAR欺骗干扰原理与MATLAB仿真：从延时移频到假目标成像
简介面向合成孔径雷达SAR电子对抗领域的欺骗干扰功能MATLAB仿真代码贴合本科、硕士阶段物理应用方向教学科研使用。程序基于MATLAB 2019a设计涵盖欺骗干扰信号构造、回波模拟、SAR成像处理等关键环节并加入压制干扰对比程序便于从原理到实现系统理解干扰机制。压缩包共57个文件其中32个.m脚本与函数、20个.asv自动备份、2个Word运行说明、2个.mat数据和1幅bmp测试图整体大小约5.78MB代码与文档组织清晰。已有413人学习下载配合运行说明可快速复现典型欺骗干扰场景适合课程设计、毕业设计或雷达对抗课题的仿真参考。1. 为什么说SAR欺骗干扰是雷达对抗里性价比最高的方案SAR成像侦察现在几乎是星载和机载平台的标配图像分辨率能做到亚米级想在强侦察下隐藏真实部署已经很难了。压制干扰的思路是拿噪声把整个场景盖掉等效功率要高出好几个数量级才能压得住匹配处理增益而欺骗干扰只用一个很小的转发功率就能在最终的SAR图像里“造”出根本不存在的虚假目标和编队。这个资源包里的qipian.m做的正是这件事截获SAR信号参数通过延时、移频和相位调制把假目标直接注入到成像结果里。对我这种做物理仿真和雷达信号处理的人来说它把从信号生成、回波模拟到干扰效果验证的整条链路都串起来了适合本科、硕士阶段做SAR对抗课题时当脚手架也适合刚接触电子战仿真的工程师快速理解欺骗干扰的作用边界。2. 欺骗干扰的成像域机理延时、移频如何变成图上的假目标2.1 距离向与方位向的“作弊”入口SAR图像的距离向分辨率靠的是大带宽线性调频LFM信号做脉冲压缩方位向分辨率靠的是合成孔径期间的多普勒历史。欺骗干扰要在这两个维度分别做文章。先看距离向。雷达发射的LFM信号可以写成% LFM基带信号K为调频斜率 t -Tp/2:1/Fs:Tp/2; s_ref exp(1j*pi*K*t.^2); % 匹配滤波参考信号匹配滤波的本质是把线性调频信号的能量压成一个窄峰峰的位置对应回波延时。干扰机如果对截获到的信号整体延时τ再转发回去SAR端匹配滤波后峰的位置就会移动c*τ/2也就等效于在距离向上制造了一个偏移目标。更妙的是因为LFM信号的“延时-频率”耦合关系对信号做固定移频f_d等效于产生一个附加延时f_d/K同样能把假目标“推”到别的距离单元上去。这就是欺骗干扰在距离向的两种基本注入手段显式延时和移频。方位向的欺骗要稍微绕一点。SAR平台飞行时真实目标与雷达之间的斜距先缩短再拉长形成多普勒历史方位压缩就是利用这个历史聚焦的。干扰机逐脉冲地给转发信号叠加一个可控的相位增量相当于伪造了一个不同的多普勒历史SAR端方位压缩后就会在某个方位位置形成虚假目标。所以一套完整的欺骗干扰仿真必须同时具备两个自由度脉冲内调制延时/移频和脉冲间调制相位编码。2.2 欺骗干扰的功率优势匹配处理增益就是你的杠杆欺骗干扰的另一个关键参数是干信比。SAR接收端对信号做匹配滤波时宽带宽积B*Tp会带来等效的处理增益。压制干扰要有效噪声功率得淹没有用信号于是干扰机功率必须比真实回波大很多而欺骗干扰只要保证截获信号后转发的功率能压过接收机噪声底匹配输出后的假目标幅度就会和真目标相当。一边是面临B*Tp量级的处理增益劣势一边是享受同样的增益两边差的功率就是干扰机的“杠杆”。从工程上说典型星载SAR的时宽带宽积动辄上千这意味着压制干扰要付出30dB以上的额外功率代价。下面这组对比在实际评估时很有参考价值。维度压制干扰欺骗干扰作用效果图像噪化、目标淹没图像出现虚假目标/编队所需功率高需对抗匹配增益低利用匹配增益信号知识依赖只需大致载频和带宽需要精确的LFM参数与PRF对抗检测难度易被识别为有源干扰假目标融入场景难以判别主要失效条件高增益天线和窄带滤波参数测量误差和算法失配2.3 欺骗干扰的“算力”瓶颈在哪里欺骗干扰在工程上最大的难点不是发射功率而是对SAR信号的实时参数测量。你要在极短的脉冲宽度内估计出载频、调频斜率、脉冲宽度和脉冲重复频率然后马上决定延时量、移频量和相位编码序列再转发出去。这个包里gonglvpu.m就是用来做功率谱分析和参数粗测的MAICHONG.m负责生成典型的SAR脉冲模板。真正到工程实现这些模块通常要跑在本振同步和FPGA高速信号链上MATLAB这份代码的价值在于先把算法链路验证通把误差敏感度摸清楚再往硬件移植就快得多。跑这套代码时我一般会先不加欺骗干扰只用回波模拟与成像模块看一个点目标的聚焦效果确认SAR成像链路本身干净再去叠加欺骗干扰模块看假目标是否按要求出现。这个顺序能省掉大量联调时间。3. MATLAB模块逐段拆解脉冲生成、功率谱与欺骗参数求解3.1 从MAICHONG.m看SAR脉冲模板设计SAR发射信号通常是宽带LFM脉冲脉冲宽度Tp决定平均功率带宽B决定距离向分辨率。MAICHONG.m本质上就是按照给定的B、Tp、采样率Fs和载频fc构造一个复数基带或中频脉冲。下面这段是从原工程思路里抽出来的教学版保持功能一致变量名直接可见% MAICHONG.m 教学版构造SAR发射脉冲模板 function s MAICHONG(B, Tp, Fs, fc) % B 信号带宽单位Hz % Tp 脉冲宽度单位s % Fs 采样率单位Hz % fc 载频单位Hz用于后续中频搬移基带仿真可置0 K B / Tp; % 调频斜率 N round(Tp * Fs); % 单脉冲采样点数 t (0:N-1) / Fs - Tp/2; % 以脉冲中心为零点的时间轴 s_b exp(1j * pi * K * t.^2); % 基带LFM信号 s s_b .* exp(1j * 2 * pi * fc * t); % 搬移到载频便于后续模拟下变频 end调用时比如B 30e6; Tp 10e-6; Fs 60e6; s MAICHONG(B, Tp, Fs, 0);就生成一个30MHz带宽、10微秒脉宽的LFM基带脉冲。注意Fs 60e6是带宽的两倍这是为了留出过采样余量避免匹配滤波后峰值旁瓣在采样点上失真。如果Fs只设成和带宽相同脉冲压缩后会看到峰值有细微的栅瓣效应幅度误差不大但做定量分析时会影响旁瓣电平读数。3.2 功率谱分析gonglvpu.m的截获参数测量思路干扰机截获到SAR信号后第一步是测带宽和中心频率。gonglvpu.m用FFT求功率谱从谱宽估计带宽。工程上更稳的做法是CZT线性调频Z变换做局部细化不过对仿真级验证来说FFT加矩形窗或汉明窗已经够用。% gonglvpu.m 教学版功率谱估计与带宽测量 function [B_est, fc_est] gonglvpu(x, Fs) N length(x); w hamming(N).; % 加窗抑制频谱泄漏 X fftshift(fft(x .* w, N)); f_axis (-N/2:N/2-1) * Fs / N; P abs(X).^2 / (N * Fs); % 功率谱密度 % 找到峰值附近功率比峰值低3dB的两个频点 [Pmax, idx] max(P); th 10^(0.3) * Pmax; % -3dB门限计算要注意单位 idx_range find(P Pmax / 2); % 半功率点 B_est f_axis(max(idx_range)) - f_axis(min(idx_range)); fc_est (f_axis(max(idx_range)) f_axis(min(idx_range))) / 2; end半功率点的选择直接对应3dB带宽这是衡量LFM信号带宽的常用口径。要注意如果输入信号里混有强干扰或多次谐波FFT谱会变宽估计出来的带宽会偏大。这时候可以先用带通滤波器把感兴趣的频段框出来再做谱分析。gonglvpu.m后接的juchi.m则是从距离维上计算目标回波的延时为欺骗干扰的延时量设定提供基准。3.3 欺骗参数求解把假目标坐标翻译成延时和移频假目标的位置由你在距离向和方位向的设定值决定。假设场景中存在一个真实点目标你要在它旁边ΔR米处放一个假目标那么对应的距离向延时就是τ 2*ΔR/c。如果要通过移频实现同样的偏移则需要的移频量为f_d K*τ。从工程上讲显式延时的好处是直观、对LFM的调频斜率不敏感缺点是需要大的延时线移频的好处是硬件实现简单一个混频器就搞定但前提是你对K的估计要准。TIAOPIN.m和TIAOXIANG.m分别对应调频和调相两种干扰调制样式适合用来对比同一假目标在不同调制方式下的成像差异。下面给出一组常用的参数设定参考参数典型取值作用带宽 B30 MHz决定距离向分辨率脉宽 Tp10 μs决定作用距离与平均功率调频斜率 K3×10^12 Hz/s决定移频-距离换算比例PRF1000 Hz决定方位向无模糊范围移频量 f_d0.55 MHz控制假目标距离偏移量相位编码步进02π控制假目标方位向位置参数表的逻辑链条是这样的你先定距离向偏移算出延时和移频再定方位向偏移算出每个脉冲上该叠的相位增量最后把这些参数输入干扰调制模块。我在调试时习惯先把f_d设成B/2以内的值超过这个范围会引起距离向频谱混叠假目标会出现成对虚影。4. 全链路复现回波模拟、欺骗叠加与SAR成像验证4.1 回波模拟huibo.m与huiboshuju.mat回波模拟是欺骗干扰的靶场。huibo.m负责生成不含干扰的真实回波huiboshuju.mat则是把回波数据存下来的结果。做干扰仿真前先有干净回波才能在成像图上分清哪个是真实目标、哪个是假目标。回波生成的核心是计算每个方位采样时刻目标到雷达的斜距并把发射脉冲按延时叠加进接收向量。% huibo.m 教学版生成单个点目标的SAR回波 function echo huibo(s_ref, R0, V, X_target, Tp, Fs, PRF, Na) % R0 最近斜距(m) % V 平台速度(m/s) % X_target 目标方位向坐标(m) % Na 方位向脉冲数 c 3e8; K 30e6 / Tp; % 从带宽和脉宽换算调频斜率 t_fast -Tp/2:1/Fs:Tp/2; % 距离向快时间 echo zeros(Na, length(t_fast)); for ia 1:Na t_az (ia - Na/2) / PRF; % 方位向慢时间 R sqrt(R0^2 (V*t_az - X_target)^2); % 瞬时斜距 tau 2*R/c; % 双程延时 idx round(tau*Fs) 1; if idx length(t_fast) echo(ia, :) s_ref .* exp(-1j*2*pi*fc*tau); % 简化的接收信号 end end end这段代码里的R用勾股定理计算没有做抛物线近似好处是在正侧视大斜视角时依然能保持较高的相位精度。慢时间轴以Na/2为中心对应合成孔径中心时刻。实际工程中为了节省内存不会把整个回波矩阵都存下来而是按方位块处理huiboshuju.mat里存的其实是一帧已经截断、去斜后的回波数据块。4.2 欺骗信号生成与叠加qipian.m的主流程qipian.m是整个工程的发动机。它从截获信号里提取参数生成调制后的欺骗脉冲再和真实回波相加。下面这个版本是经过我简化重写的教学主干完整工程里还挂了.asv自动保存的中间版本跑的时候直接认run.m即可% qipian.m 教学版欺骗干扰主流程 % 假设已经截获并测得了信号参数fs, B, Tp, fc, PRF s_jam_all zeros(size(echo)); % 预分配干扰矩阵 for ia 1:Na % 距离向移频 延时 t -Tp/2:1/Fs:Tp/2; s_jam exp(1j*pi*K*t.^2); % 复制真实脉冲 s_jam s_jam .* exp(1j*2*pi*fd*t); % 移频fd实现距离向偏移 tau_deceive 2*R_target/c; % 假目标设定的距离延时 n_offset round(tau_deceive * Fs); s_jam circshift(s_jam, n_offset, 2); % 循环移位模拟延时 % 方位向叠加相位步进伪造多普勒历史 phase_az exp(1j * 2*pi * ia * fdaz / PRF); s_jam_all(ia, :) phase_az * s_jam; % 干信比控制按所需JSR放大干扰幅度 s_jam_all(ia, :) s_jam_all(ia, :) * sqrt(10^(JSR/10)); end echo_jam echo s_jam_all; % 真实回波上面叠假目标这里circshift是循环移位注意它会把信号尾部卷到头部来。对于远距离目标延时远大于脉宽的场景这种卷绕会造成距离模糊假象正确做法是开一个长度Na×Ns的大矩阵用n_offset直接索引写入。上面代码为了保证可读性用了循环移位实际工程我改成线性写入。JSR干信比一般设成-10到10dB之间就能看到明显的欺骗效果太高会连假目标旁瓣一起变亮反而暴露干扰痕迹。4.3 成像验证距离压缩、距离徙动校正与方位压缩欺骗干扰成不成功最终要在SAR图像上看。最基本的验证链路是先距离压缩再做距离徙动校正RCMC最后方位压缩。下面给出一个不含插值的简化RD算法流程。% 距离压缩匹配滤波 echo_rc ifft(fft(echo_jam, Nf, 2) .* conj(fft(s_ref, Nf, 2)), Nf, 2); echo_rc echo_rc(:, 1:length(t_fast)); % 距离徙动校正按最近斜距计算徙动曲线逐距离门搬移简化最近邻 for ia 1:Na R_true sqrt(R0^2 (V*(ia-Na/2)/PRF)^2); dr 2*(R_true - R0)/c * Fs; % 需要校正的采样点偏移 shift round(dr); if shift ~ 0 echo_rc(ia, :) circshift(echo_rc(ia, :), -shift, 2); end end % 方位压缩对每个距离门做方位向FFT img fftshift(fft(echo_rc, Na, 1), 1); img abs(img).^2;这段代码的用意是把回波先压缩成距离-慢时间矩阵再沿方位向压缩得到聚焦图像。干扰是否生效看两个特征第一假目标是否出现在你设定的距离和方位坐标上第二假目标的聚焦质量是否和真目标一致。如果假目标出现但明显散焦基本可以判断是移频量或相位编码步进设错了如果假目标出现在镜像位置则是方位向频谱混叠需要把PRF调大或降低相位编码变化速率。这个包的lena.bmp不是拿来当目标图像的而是用来模拟SAR场景后向散射系数分布把图像灰度值映射成地面各点的反射率再逐点生成回波。这样你能直观地看到欺骗干扰在一个真实纹理场景里的视觉效果比单纯点目标有说服力得多。5. 参数边界与工程排错五个典型坑与新版MATLAB迁移5.1 假目标位置偏移的定量校验调试时最常问的问题是“我设了ΔR和ΔX为什么图上假目标不在预期位置”。原因多半是单位换算里漏了双程因子。距离向延时换算用τ 2ΔR/c移频换算用ΔR c*f_d/(2K)。第二项要特别注意K的正负上扫频和下扫频的LFM符号不同移频产生的距离偏移方向正好相反。在qipian.m里改干扰前先打印一组换算结果和成像图上读出的峰值位置对一下通常十分钟内能定位问题。方位向的偏移量则按ΔX V_fake * t_az算这里的V_fake是你通过相位步进伪装的等效速度不是平台速度。5.2 多普勒模糊与成对回波当干扰相位步进对应的多普勒频率超过PRF/2时就会产生方位向混叠假目标会出现在真实的镜像位置形成成对目标。解决思路是限制最大相位变化率或者提高PRF重频。工程上另一个常见现象是假目标在方位向“拉丝”这通常是因为相位编码序列存在固定偏差导致多普勒中心偏移未完全补偿试试在方位压缩前对慢时间信号做去斜处理。5.3.asv文件和新版MATLAB的兼容性包里那些.asv是MATLAB自动保存的草稿不是必须的运行文件。拿到工程后先看运行说明.doc再直接运行run.m。如果用的不是R2019a而是新版比如R2023b或R2026b需要注意新版对circshift和矩阵索引规则没有重大改动但fftshift的默认维度行为在R2021a后更严格了代码里最好明确写出fftshift(X, 1)这样的维度参数否则二维矩阵时可能得到和旧版不一致的显示结果。另外新版MATLAB的plot配色和坐标轴属性默认值变了不影响数据计算但脚本里若有依赖旧版坐标轴颜色的后续叠加绘图要检查一遍。5.4 与压制干扰的对比判断信号“是否被识破”从接收端视角回看这套代码可以额外做个实验对同一份真实SAR数据分别加噪声压制干扰和欺骗干扰再用CFAR检测统计虚警率和漏警率。欺骗干扰产生的假目标如果功率调得和真实目标接近CFAR的检测门限完全无法区分压制干扰则直接把背景噪声抬高CFAR门限上移真实目标反而可能被淹掉并触发全局告警——从隐蔽对抗的角度讲欺骗干扰明显更难被识破这也是它被频繁用于SAR对抗仿真教研的原因。如果你想在这个包里继续深入可以试着把净输出改成多个假目标沿同一方向同步移动观察SAR成像结果会不会出现“幽灵编队”效应。本文还有配套的精品资源点击获取