简介一份基于Matlab实现的体制雷达信号仿真源码包面向雷达信号处理方向的科研人员、高校学生及算法工程师解决难以在真实设备上验证多种雷达信号特性的问题。资源包共8个文件全部为m脚本压缩后仅5KB轻量易用覆盖单体制、频率分集体制、重频参差与抖动、PRI跳变与滑变、脉组PRI变化及双脉冲等典型信号类型。已有403人学习下载适合用于教学演示、课程设计或算法验证。通过运行脚本读者可直观对比不同脉冲重复间隔策略下的信号波形与频谱差异理解杂波抑制、抗干扰及距离/速度分辨等核心概念代码结构简洁、模块划分清晰既可作为快速上手的仿真模板也能在此基础上扩展复杂雷达系统实验为后续研究和工程应用提供坚实基础。1. 体制雷达信号仿真PRI变化为什么是抗干扰的第一课我最早接触“体制雷达信号仿真”时以为就是把一个脉冲串画出来直到被要求复现重频参差和重频抖动对盲速的改善才发现脉冲重复间隔PRI的调制方式直接决定了雷达能不能在强杂波里把动目标捞出来。这份源码包把单体制、频率分集、重频参差、重频抖动、PRI跳变、PRI滑变、脉组PRI变化和双脉冲八类体制雷达信号摔在你面前每个都有可直接跑的Matlab实现。它适合两类人一类是做雷达系统课程设计、需要交波形图和频谱图的学生另一类是刚接手雷达信号处理、想快速验证PRI变化对距离-速度模糊影响的工程师。别指望它替代完整雷达系统仿真它解决的是“波形怎么生成、参数怎么设置、后续FFT和模糊函数怎么接”这一段正好是多数人卡住的地方。2. 单体制与频率分集先把两种基线信号在Matlab里跑通2.1 单体制雷达信号一个相干脉冲串从参数到波形单体制雷达信号是所有PRI变化样式的地基。它指雷达以固定脉冲重复间隔发射相同载频的脉冲串每个脉冲内部是矩形包络加一个中频或射频载波。仿真时最常见的选择是生成复数基带信号也就是把载频移到复数指数上这样采样率不用拉太高后续匹配滤波也方便。fs 100e6; % 采样率 100 MHz T 1e-6; % 脉宽 1 us PRI 100e-6; % 脉冲重复间隔 100 us fc 10e6; % 中频载频 10 MHz N 8; % 脉冲个数 t_all 0:1/fs:N*PRI-1/fs; s zeros(1, length(t_all)); for i 1:N % 当前脉冲在整段时间的起止索引 idx round((i-1)*PRI*fs)1 : round(i*PRI*fs); t_loc t_all(idx) - (i-1)*PRI; % 脉冲内相对时间 s(idx) rectpuls(t_loc - T/2, T) .* exp(1j*2*pi*fc*t_all(idx)); end这段代码的关键在于指数里用的是全局时间t_all而不是局部时间。如果写成exp(1j*2*pi*fc*t_loc)每个脉冲起始相位都会清零脉冲串的相位就断了后续做相干积累会看到频谱出现奇异尖峰。用全局时间才能保证脉冲之间的相位连续这也是相干雷达仿真和普通“画波形”的本质区别。参数说明fs决定了时间步进至少要满足采样定理。对于中频载频10 MHz采样率100 MHz意味着每个载频周期采10个点波形平滑度足够。T与PRI的比值是占空比这里为1%接近真实雷达脉冲宽度越宽发射能量越大距离分辨率越差做课程设计时这两个参数要按实际需求折中。2.2 频率分集体制让载频在脉冲间跳变频率分集是为了抗干扰和改善目标起伏检测。它和单体制的区别只有一个每个脉冲的载频不再是固定值而是从一组频点里按顺序取。实现方式几乎不用改上面的骨架只需要把fc换成向量并循环取值。fc_set [10e6, 10.5e6, 11e6, 11.5e6]; % 4个频率点 N 8; s zeros(1, length(t_all)); for i 1:N idx round((i-1)*PRI*fs)1 : round(i*PRI*fs); t_loc t_all(idx) - (i-1)*PRI; fc_i fc_set(mod(i-1, length(fc_set)) 1); % 循环取频点 s(idx) rectpuls(t_loc - T/2, T) .* exp(1j*2*pi*fc_i*t_all(idx)); end注意两点一是mod取余保证索引不越界二是频率分集后每个脉冲的载频不同事后做包络对齐的测速处理时要从s里分离出各频点对应的距离门再做FFT否则多个载频的相位混在一起测速值会乱。这个在后续章节做避坑单独讲。2.3 两种波形的典型参数对照参数单体制频率分集载频固定单值频点集循环PRI固定可固定也可变抗干扰能力弱较强处理复杂度低中典型应用常规搜索抗有源干扰频率分集在工程上常与重频参差叠加使用即在跳频的同时PRI也按固定节拍变化。先单独跑通单频段再叠加跳频调试时定位问题会快很多。3. 重频参差、抖动与PRI跳变三种抗盲速波形的实现差异3.1 重频参差用固定小PRI拼出大周期重频参差的本质是让PRI在几个固定值之间循环比如一个“参差周期”内先发100微秒间隔再发120微秒再发110微秒然后重复。它的作用是消除盲速——单个PRI下的盲速对应多普勒频率等于PRI倒数而用多个不同PRI盲速点不会重合相当于把盲速推远了。PRI_set [100e-6, 120e-6, 110e-6]; % 参差PRI序列 N_cycle 100; % 重复多少个参差周期 s []; t_start 0; for k 1:N_cycle for m 1:length(PRI_set) PRI PRI_set(m); n round(PRI * fs); t t_start : 1/fs : t_start PRI - 1/fs; s [s, rectpuls(t - T/2 - t_start, T) .* exp(1j*2*pi*fc*(t - t_start t_start))]; t_start t_start PRI; end end实现时最容易犯的错是用for叠加造成数组反复拼接数据量大时极慢。正确做法是先计算总长度预分配零向量再填值。上面的写法为了可读性用了拼接但如果你跑10万个脉冲建议改成预分配。另外参差PRI不是随机抖动它有严格的周期结构。处理端会把一个完整的参差周期作为相干处理时间而不是像固定PRI那样直接对连续脉冲做FFT。这决定了后续匹配滤波的参考时间轴必须与发射的时间轴一致。3.2 重频抖动中心PRI加随机扰动重频抖动与参差不同它的PRI是围绕一个中心值随机变化的。目的是破坏干扰机对信号规律的预测也用于降低目标多普勒模糊的固定性。PRI0 100e-6; % 中心PRI dither 0.2; % 抖动深度20% rng(42); % 固定随机种子方便复现 N 1000; PRI_seq PRI0 * (1 dither * (2*rand(1,N) - 1)); % 保证最小间隔防止脉宽重叠 PRI_seq max(PRI_seq, T*5);抖动深度是核心参数。工程实物中抖动深度一般不超过20%30%太大则脉冲回波在距离维上散开积累增益下降。抖动数据要保留一份PRI_seq存成.mat文件因为后端处理时要用发射的真实PRI顺序去构造匹配滤波参考随机数种子一旦变整个仿真结果就无法复现。3.3 PRI跳变伪随机顺序的PRI序列重频跳变和抖动的区别在于抖动是连续分布跳变是从有限的几个PRI值里按伪随机顺序挑选。比如固定PRI为100、110、120微秒每来一个脉冲随机选一个但长期统计上三个值出现次数接近。PRI_pool [100e-6, 110e-6, 120e-6]; rng(7); idx randi([1, 3], 1, N); PRI_seq PRI_pool(idx);这里的randi生成的是均匀随机索引短期序列可能出现连续两个相同PRI。工程上会采用“不重复抽取”的策略把三个索引做成一个序列perms或 shuffle 后循环播放保证每三个脉冲内三个PRI各出现一次。否则连续相同PRI容易让后端误判为固定PRI。3.4 三种PRI变化方式的适用场景对比类型随机性来源对消盲速效果处理复杂度适用场景重频参差无固定循环好盲速被推远中常规MTI重频抖动连续随机残差噪声高低抗欺骗干扰PRI跳变离散伪随机介于两者之间中搜索与跟踪切换选型时我把重频参差作为默认首选因为它的波形设计确定性强后端好处理抖动和跳变则更多用在对抗特定干扰的场景。三者不是互斥的现代雷达经常把参差和抖动混合但仿真时一定要分开验证再谈组合。4. 滑变、脉组PRI与双脉冲更复杂的体制雷达波形怎么搭4.1 PRI滑变PRI按线性斜率连续变化PRI滑变滑变重复频率的PRI不固定而是按照某个斜率连续递增或递减。它的特点是无级变化不像参差那样阶跃常用于多普勒解模糊与目标距离跟踪同时需要的情况。PRI_start 100e-6; slope 1e-9; % 每个脉冲增加1 ns N 500; PRI_seq PRI_start (0:N-1) * slope; % 计算发射时刻 tx_time cumsum([0, PRI_seq(1:end-1)]);cumsum把PRI序列累加成发射时刻数组这一步是通用的。很多初学者只生成PRI_seq却不生成tx_time后面匹配滤波时根本不知道每个脉冲对应什么时刻导致画出的时域波形和脉冲实际位置对不上。PRI滑变有个坑斜率过大会让远端PRI变得很大相干处理时间内脉冲数减少速度分辨率变差。所以斜率选择要看你在乎的是距离无模糊还是速度无模糊两者是矛盾的。4.2 脉组PRI变化多组波形时间分集脉组PRI变化也叫重频参差组间变化意思是把整个工作时间分成若干组组内PRI固定组与组之间PRI不同。比如第一组用100微秒发50个脉冲第二组用120微秒发50个脉冲第三组用110微秒发50个脉冲。它的好处是每组内部仍是固定PRI可以沿用成熟的FFT测速流程组间PRI变化又提供了多普勒解模糊所需的多组观测。group_PRI [100e-6, 120e-6, 110e-6]; pulses_per_group 50; s []; tx_time []; for g 1:length(group_PRI) PRI group_PRI(g); for k 1:pulses_per_group n round(PRI * fs); t (k-1)*PRI : 1/fs : k*PRI - 1/fs; s [s, rectpuls(t - T/2, T) .* exp(1j*2*pi*fc*t)]; tx_time [tx_time, (k-1)*PRI]; end end注意这里每组内部我重新开始了相位参考也就是每组第一个脉冲相位从0开始。这样做符合部分真实雷达的分组处理逻辑但组与组之间的相位是断裂的。如果你要研究组间相干就必须把tx_time累加到全局时间上用全局时间的复指数相位就连续了。两种处理在后续模糊函数图上差异明显。4.3 双脉冲信号两脉冲一组的特殊调制双脉冲体制不是简单发两个脉冲而是两脉冲构成一组组内两个脉冲间隔很短组间间隔较长。这样做的目的是用组内短间隔获得高的多普勒不模糊速度用组间间隔获得距离不模糊范围。PRI_inner 50e-6; % 组内间隔 PRI_outer 500e-6; % 组间间隔 N_group 100; s []; for g 1:N_group base_tx (g-1) * PRI_outer; % 第一个脉冲 s [s, rectpuls((0:1/fs:PRI_inner-1/fs) - T/2, T) .* exp(1j*2*pi*fc*((0:1/fs:PRI_inner-1/fs) base_tx))]; % 第二个脉冲组内延迟PRI_inner t2 PRI_inner : 1/fs : 2*PRI_inner-1/fs; s [s, rectpuls(t2 - PRI_inner - T/2, T) .* exp(1j*2*pi*fc*(t2 base_tx))]; end组内间隔取值要远大于脉宽否则两个脉冲会重叠。仿真时还要考虑发射脉冲占空比对雷达平均功率的影响同样的脉宽和峰值功率下双脉冲体制在PRI_outer周期内只发射两小段能量平均功率会比常规固定PRI低目标检测距离会缩短这是方案取舍时要付出的代价。5. 避坑与排查Matlab做体制雷达仿真最容易翻车的5个点5.1 采样率没跟上载频波形直接失真现象生成的信号时域图看起来正弦波变形成锯齿频谱图上在载频附近出现密密麻麻的镜像。原因采样率fs小于2倍载频或者载频高到采样点不足。我在第一次做频率分集时把载频设成30 MHz采样率只有20 MHz结果30 MHz信号被混叠到10 MHz后面匹配滤波全错。解决要么把载频降低到采样率的1/10以下例如fs100e6时fc10e6要么显著提高fs到至少5倍载频以上。我的习惯是直接设fs为载频的10倍仿真单脉冲时跑得也不慢波形与解析计算高度吻合。5.2 随机抖动不设种子两次结果对不上现象同一套代码跑两次PRI序列不同频谱图和模糊函数图完全不一样结论无法复现。原因rand和randi每次调用都会从系统时钟取种子没有人工固定。解决凡是涉及抖动的脚本第一行写rng(固定的整数)。我这边项目里固定用rng(2024)作为标定参数并把生成的PRI_seq连同波形一起保存到.mat。后续换电脑或换Matlab版本也能对得上结果。5.3 脉组切换时相位不连续频谱出现杂散现象脉组PRI波形的FFT频谱在目标多普勒附近出现额外的对称旁瓣看起来像多了一个假目标。原因每组起始相位被重置实际信号相当于对载波做了分段门控门控沿在频域产生扩散。解决统一用全局时间轴生成所有脉冲复指数写exp(1j*2*pi*fc*t_global)不要在组内重新累加。需要验证时对比局部时间版本和全局时间版本的频谱后者杂散会明显降低。5.4 频率分集后做FFT测速距离-速度模糊一团糟现象对频率分集信号直接做跨脉冲FFT峰值不在真实多普勒位置。原因不同载频的脉冲相位随距离变化率不同直接混在一起测速相当于对不同频率点做非相干积累模糊函数被线性的相位差污染。解决先按载频把脉冲分组每个频点单独做距离门匹配滤波得到各频点的幅度信息再对不同频点的幅度进行合成或检测而不是对原始复数回波做跨脉冲FFT。仿真里可以做两个版本对照一个混在一起测速一个分组处理你会立刻看到后者的速度峰清晰得多。5.5 把波形生成和参数配置写在一起改到崩溃现象一个脚本从头写到尾参数散落在各个for循环里想改PRI要翻好几处漏改一处结果全错。原因没有做参数与波形分离。Matlab脚本天然容易写成“直通式”参数、生成、画图混在一起。解决我一般会建一个参数结构体param所有仿真入口读同一份配置。param.fs 100e6; param.fc 10e6; param.T 1e-6; param.signal_type stagger; % 可选single/freq/stagger/dither/hop/slide/group/double param.PRI_set [100e-6, 120e-6, 110e-6]; param.N 1000;波形生成函数只接收param输出s和tx_time。这样换信号类型只需改signal_type跑完自动存param到.mat整条链路可回溯。6. 一个具体技巧用模糊函数验证你的PRI设计6.1 模糊函数的Matlab快速实现PRI变没变对光看时域波形说明不了问题模糊函数图才是一眼看清多普勒模糊结构的工具。它计算发射信号与自身经过时延和多普勒频移后的相干积分能直接展现距离-多普勒耦合情况。function [AF, tau, fd] ambiguity_function(s, fs, Fd_max) N length(s); tau -T:PRI/100:T; % 时延范围按脉冲宽度和PRI取 fd -Fd_max:1/(N/fs):Fd_max; AF zeros(length(tau), length(fd)); for i 1:length(tau) delay round(tau(i) * fs); s_shift zeros(1, N); if delay 0 s_shift(delay1:end) s(1:end-delay); else s_shift(1:enddelay) s(-delay1:end); end for j 1:length(fd) AF(i,j) abs(sum(s .* conj(s_shift) .* exp(1j*2*pi*fd(j)*(0:N-1)/fs))); end end end这个实现复杂度是O(N^2)只适合短序列验证。跑20个脉冲以内没问题超过100个建议用FFT近似算法。模糊图的纵轴时延、横轴多普勒固定PRI会出现一个突出的峰墙重频参差会把峰墙打碎成多个离散平台抖动则直接让旁瓣铺平。看到哪种形态就知道PRI设计是否达到意图。6.2 从模糊图读PRI变化的特征用重频参差信号跑出来的模糊图中你会发现主峰旁边每隔一段多普勒距离出现一个小峰这些小峰的间距正好对应参差周期的倒数。如果小峰高度接近主峰的一半说明参差周期设计不够远盲速改善并不彻底。重频抖动时模糊图不再是尖锐峰而是变成一个宽底的泡状说明单个脉冲的相参积累能力被随机PRI稀释了。从那以后我每次生成完体制雷达波形都不急着出频谱图而是强制自己先跑一遍模糊函数和距离维频谱检查确认PRI变化的特征符合理论值再继续搭后端处理。这么做能拦住八成以上的参数设计错误省下的时间远比写函数花的多。希望帮到你。这份源码里已经把八种波形的生成函数和配套demo整理好了拿回去把param里的参数按你的指标改一改就能直接看到每种体制的模糊图和频谱。本文还有配套的精品资源点击获取