2026/9/8 9:16:14
FEATURED · 精选文章 
简介一份面向通信工程专业学生、研究人员及算法开发者的MATLAB仿真资源聚焦MFSK/FH高速跳频通信系统的建模与性能分析。该系统充分利用MFSK在有限带宽内传输更多信息、抗多径衰落能力强的特点以及FH快速跳变载波带来的抗干扰、抗截获和隐蔽性优势可服务于军事抗干扰通信、无线局域网等应用场景。仿真模型覆盖完整的信号发射、传输和接收流程支持自定义跳频图案、调制阶数与信道参数方便对比不同设置下的误码率表现也可深入研究跳频图案对通信安全性的提升作用。压缩包为zip格式共包含2个文件一个可直接运行的MATLAB脚本.m与一份说明文档.md整体仅4KB体积紧凑却便于快速查阅。目前已有92人学习下载。通过运行主程序可看到MFSK调制解调和FH跳频控制的实现细节结合说明文档可快速掌握参数调整与结果分析方法适合作为扩频通信课程设计、毕业设计或科研预研的速查参考。 接到这个选题时我第一反应是写MFSK/FH仿真的文章不少但大多数都停在贴代码跑通一版的层面真正把为什么要这么设计参数误码率曲线背后反映什么问题讲透的很少。这次我把完整做下来的经验和踩过的坑整理出来希望能给正在做通信仿真或者准备做毕业设计、课题预研的朋友一些参考。从仿真目标反推系统参数MFSK/FH链路该怎么拆很多人拿到高速跳频通信系统仿真这个题目第一件事就是打开MATLAB开始写代码。我劝你先停下来花半天时间把链路框图和数据流画清楚。因为MFSK和FH虽然原理都不复杂但合在一起之后采样率、符号速率、跳速、频率间隔这些参数会互相牵制参数设不合理后面跑出来的曲线自己都解释不了。为什么MFSK和FH是天生一对MFSK多进制频移键控是一种频率调制方式M进制意味着每个符号携带log2(M)比特信息用M个不同频率中的一个来承载。FH跳频则是把整个信号按照跳频图案在很宽的频带内搬移两者结合后信号在任一时刻既携带MFSK的调制频率信息又位于跳频图案指定的载波频率上。从抗干扰的角度看MFSK给系统带来了频率分集的效果某个频率点被干扰时其他频率点仍然可以正常工作FH则提供了处理增益干扰机不知道下一跳跳到哪个频率只能全频带压制代价极高。更关键的是MFSK是非相干检测友好的调制方式而跳频通信中接收机很难保持本振与发射机载波的严格相位相干MFSK天然适合这种场景。这套组合还有个工程上的好处实现门槛低。相比于需要精确相位跟踪的PSK类调制MFSK的调制解调用模拟器件或者DSP都能很方便实现这也是很多战术电台选择它的原因。确定仿真参数的关键考虑参数设计是整个仿真的起点也是我踩坑最多的地方。先给出一套完整可用的参数后面所有分析和代码都基于这套配置数据速率Rb20 kbps调制阶数M4即每个符号携带2比特符号速率Rs10 ksymbol/s频率间隔Δf25 kHz近似等于符号速率的设计值跳频速率2000 hop/s每跳承载5个符号跳频频率集规模Nhop64个频点频率间隔5 MHz总跳频带宽320 MHz采样率fs10 MHz单个跳频频点内看10 MHz采样足够容纳25 kHz的MFSK信号并留出余量这套参数下单跳内要处理5个MFSK符号跳频图案每毫秒更新一次。仿真时有两种建模方式一种是基带等效建模把所有频点搬移到基带处理速度快但不好观察频谱另一种是带通建模生成真实的射频波形直观但计算量大。我的建议是先用基带建模验证算法再抽出一小段做带通波形验证兼顾效率和可信度。MFSK调制的MATLAB实现为什么非相干解调更适合跳频MFSK的调制在MATLAB里其实就一句话根据输入的比特映射到对应的正弦频率。但解调方案的选择直接决定了接收机的结构和仿真代码的复杂度。频率映射与调制实现4FSK系统里每2个比特映射为一个符号频率。我采用的映射关系是00 - f0 100 kHz01 - f1 125 kHz10 - f2 150 kHz11 - f3 175 kHz注意这里的频率是相对跳频载波的偏置信号实际发射频率是跳频载波加上这个偏置。关键点在于频率间隔的取值。MFSK的非相干解调要求频率间隔不小于符号速率即Δf ≥ Rs否则相邻频率的频谱会发生严重重叠误码率会显著恶化。实践中我取Δf 2.5Rs留出一定保护带这样即使用矩形脉冲成型相邻频点间的串扰也可以忽略。调制代码核心部分如下% 参数定义 M 4; % 调制阶数 Rs 10000; % 符号速率 10 ksymbol/s T_symbol 1/Rs; % 符号周期 fs 10e6; % 采样率 10 MHz samples_per_symbol fs/Rs; % 每个符号的采样点数 % 频率映射表 freq_table [100e3, 125e3, 150e3, 175e3]; % 比特到符号索引的映射 Gray编码 bit_to_sym [0 1 3 2]; % 对输入比特流分组并调制 bits randi([0 1], 1, 10000); symbols_idx bit_to_sym(2*bits(1:2:end) bits(2:2:end) 1); % 生成MFSK调制信号单跳 t (0:samples_per_symbol-1)/fs; tx_signal []; for k 1:length(symbols_idx) tx_signal [tx_signal, cos(2*pi*freq_table(symbols_idx(k))*t)]; end这里用了一个细节比特到符号的映射采用了Gray编码映射00对应001对应111对应210对应3。Gray编码的好处是相邻符号之间只有1个比特差异当噪声导致解调器误判到相邻频率时只会错1个比特而不是2个。这个映射在高信噪比区域能让误码率直接减半属于投入产出比极高的小技巧。非相干解调的判定逻辑之所以弃用相干解调原因前面说过跳频系统里收发双方的本地振荡器频率难以完全同相加上快速跳频带来的频率切换瞬态做载波相位恢复非常痛苦。非相干解调则完全不在乎相位直接用能量检测。实现方式也简单接收端对每个符号周期做FFT然后找幅度最大的频率点对应的索引反查映射表得到比特。% 非相干解调 rx_symbols zeros(1, length(symbols_idx)); for k 1:length(symbols_idx) seg rx_signal((k-1)*samples_per_symbol1 : k*samples_per_symbol); % 加窗后FFT减少频谱泄漏 win hann(length(seg)); spectrum abs(fft(seg .* win)); % 只需要关注M个频点附近的幅度 [~, max_idx] max(spectrum(1:length(spectrum)/2)); % 找到对应的频率索引 [~, rx_freq_idx] min(abs(freq_table - max_idx*fs/length(seg))); rx_symbols(k) rx_freq_idx; end实际做的时候我加了汉宁窗后再做FFT因为矩形窗在频点非整数周期截断时会产生严重频谱泄漏干扰判决结果。跳频图案与频率合成器仿真里最容易拖慢速度的环节FH的仿真核心有两块伪随机序列的生成以及跳频图案到实际发射频率的映射。这两块做好了整个系统才有跳的灵魂不然就只是个MFSK仿真。伪随机序列与跳频图案生成跳频图案通常由一个伪随机序列控制频率合成器产生。仿真中用m序列最长线性反馈移位寄存器序列来生成其自相关性好、实现简单、具有伪随机性。m序列生成器的核心是反馈多项式。我用的本原多项式是x^8 x^6 x^5 x^4 1即8级移位寄存器周期255足够支撑64个频点的跳频图案而不重复。% 生成m序列 reg ones(1, 8); % 寄存器初始状态 seq_len 255; m_seq zeros(1, seq_len); for n 1:seq_len m_seq(n) reg(8); feedback xor(reg(8), xor(reg(6), xor(reg(5), reg(4)))); reg [feedback, reg(1:7)]; end % 映射到跳频频率每跳需要6比特控制64个频点 % 串并转换取6比特转10进制 hop_pattern zeros(1, floor(seq_len/6)); for h 1:length(hop_pattern) bits6 m_seq((h-1)*61 : h*6); hop_pattern(h) bi2de(bits6, left-msb); end % 通过异或加扰方式增加复杂度这里简化直接用映射 hop_freqs hop_freq_table(hop_pattern 1);注意一个细节m序列的周期255不能被6整除所以最后会浪费3个比特。实际系统中通常会用多个m序列组合生成更长的跳频图案仿真阶段不用追求极致的随机性能但至少保证跳频图案在仿真时长内不出现明显周期。跳频信号拼接的工程细节MFSK信号和跳频图案的时序配合是仿真中容易出问题的地方。跳频速率2000 hop/s意味着每0.5毫秒切换一次载波频率每跳内正好承载5个MFSK符号。发射信号的表达式是s(t)cos(2π(f_hop(n)f_mfsksymbol)·t)但这里有个隐含问题频率切换瞬间相位不连续。理想跳频系统允许相位不连续实际也能工作因为接收端用的是非相干解调。但在FFT分析跳频信号时会看到很宽的频谱扩散。仿真中我建议在每跳的边界处留出几个采样点的保护间隔更贴近工程中频率合成器建立时间。% 拼接带跳频的发射信号 tx_hop_signal []; hop_duration 0.5e-3; % 每跳0.5ms hop_samples hop_duration * fs; for h 1:length(hop_pattern) % 本跳的MFSK基带信号 start_sym (h-1)*5 1; end_sym h*5; hop_baseband tx_signal((start_sym-1)*samples_per_symbol1 : end_sym*samples_per_symbol); % 搬移到跳频频率 t_hop (0:hop_samples-1)/fs; hop_rf hop_baseband .* cos(2*pi*hop_freqs(h)*t_hop); tx_hop_signal [tx_hop_signal, hop_rf]; end这段代码跑起来你会发现一个问题循环拼接很慢数据长了之后内存消耗也大。更好的做法是预分配数组或者用矩阵化操作。跳频仿真尤其要注意性能优化我后面单独说。完整链路跑通后的性能分析误码率曲线能告诉我们什么系统搭完之后接下来才是重头戏把信道效应加进去跑不同信噪比下的误码率画出BER曲线分析系统性能。这一步做得好不好直接决定这篇仿真有没有说服力。高斯信道下的基准性能首先在AWGN加性高斯白噪声信道下测试系统的本征性能。这是所有性能分析的基准线任何一个跳频系统的误码率都不会优于同参数MFSK在AWGN下的理论值。4FSK非相干解调的理论误码率为P_e ∑_{k1}^{M-1} (-1)^{k1} C(M-1,k)/(k1) · exp(-k·Eb/(k1)·N0)我实测下来仿真曲线和理论值在低信噪比区域贴合得很好高信噪比区域也就是在Eb/N0大于14dB以后由于仿真数据量不足开始出现抖动。要解决这个问题至少要保证每个信噪比点统计到200个以上的错误比特低误码率区域建议统计到500个错误比特再停下来。干扰场景下的抗干扰增益跳频系统最核心的卖点是抗干扰。仿真中我设置了两种典型的干扰场景单频点连续干扰和宽带部分频带干扰。单频点干扰的模型是在某个固定频率上叠加强干扰信号。由于跳频图案的作用这个干扰只在系统跳到被干扰频点的那一刻产生影响。假设跳频有64个频点干扰占用1个频点那么系统只有1/64的时间被干扰到这就是处理增益。64个频点对应18dB的处理增益实际仿真跑下来的结果和这个理论值吻合。部分频带干扰模型更贴近实际情况干扰机不知道跳频图案但知道系统的工作频段于是压制其中一部分频带。比如干扰25%的频带那么系统约有四分之一的跳频点被干扰干扰频带内的符号会大量出错但依然有75%的频点能正常工作。仿真结果显示在干扰条件下MFSK/FH系统的误码率并不是线性恶化而是呈现明显的平层效应——即当干扰足够强时误码率不再随干扰功率显著变化而是被未被干扰的频点比例所决定。% 部分频带干扰仿真 interference_ratio 0.25; % 干扰25%频带 interfered_freqs randperm(length(hop_freq_table), round(interference_ratio * length(hop_freq_table))); % 在接收信号上叠加干扰 for h 1:length(hop_pattern) if ismember(hop_freqs(h), interfered_freqs) t_hop (0:hop_samples-1)/fs; rx_hop_signal((h-1)*hop_samples1 : h*hop_samples) ... rx_hop_signal((h-1)*hop_samples1 : h*hop_samples) ... A_jammer * cos(2*pi*hop_freqs(h)*t_hop rand*2*pi); end end同步误差对性能的影响很多人在仿真里默认收发双方完全同步这是最理想化的假设也是仿真结果和实测差距最大的地方。我专门做了一组对比当接收机的跳频图案相对发射机存在时间偏移时误码率如何变化。结果很有意思时间偏移低于一个MFSK符号周期的20%时性能几乎无损偏移达到半个符号周期误码率开始明显上升超过一个符号周期系统彻底无法工作。这主要是因为跳频同步误差会导致接收端解调窗口错位每个符号周期的FFT窗口内包含了两个不同符号的波形能量检测判错。这组数据说明了一个工程结论跳频系统的同步精度要求大约在符号周期的20%以内这也为接收机同步算法的设计提供了数字化指标。跳频仿真中那些只有跑过才懂的坑代码调通只算完成了一半真正耗时间的往往是那些看起来不起眼却又致命的细节。我把这次仿真中遇到的几个典型问题列出来希望能帮你少走弯路。频谱泄漏FFT频率分辨率不够怎么办MFSK解调用的是FFT找最大能量点但当符号周期不是FFT频率分辨率的整数倍时频谱能量会泄漏到相邻频点导致判决出错率高。我实测发现当采样率10MHz、FFT点数取1000时频率分辨率正好是10kHz而频点间隔是25kHz理论上足够区分。但矩形窗的旁瓣衰减只有13dB遇到噪声稍大就可能在旁瓣位置出现虚假峰值。解决方案有两个一是换窗函数汉宁窗旁瓣衰减可达31dB代价是主瓣变宽二是增加FFT点数用零填充到2048点能有效降低栅栏效应带来的幅度误差。实际我两个方法都用了误码率在高信噪比区域明显改善。跳变沿的毛刺忽略它会让误码率多一个数量级跳频频率切换不是瞬时完成的仿真中如果简单地在每跳边界处硬切换频率会在时域上产生一个相位突变造成很宽的频谱扩散。虽然在纯数学仿真中不影响MFSK解调但一旦加入带限信道滤波器毛刺就会变成真正的干扰。我在每跳边界的前后各留了20个采样点作为保护间隔将频率切换放在保护间隔内完成。这样每个跳频周期内有效数据只占96%损失可以忽略但信号频谱干净得多。仿真速度矩阵化是唯一出路跳频仿真最大的痛点是速度。如果坚持用for循环逐跳拼接信号仿真到100万比特数据量时会非常慢。我的做法是尽可能用MATLAB矩阵化操作代替循环。一个很实用的技巧是先生成完整的MFSK基带信号再用一个频率切换矩阵对信号进行分块搬移把逐跳乘载波变成向量点乘。此外如果只是看误码率趋势完全可以在基带等效域完成仿真也就是不显式生成射频频率而是给每个频点分配一个等效噪声和干扰功率速度可以快上数十倍。我建议先跑基带等效仿真摸清参数再跑带通仿真验证几个关键点。误码率统计的置信度问题还有一个容易忽略的问题误码率统计的样本量。很多人跑到感觉曲线光滑就停了但低误码率区域需要极大的数据量才能保证置信度。我的经验是每一个信噪比点错误比特数至少统计200个如果误码率低于1e-4统计到500个错误比特再停。这么做代价是仿真时间显著增加但换来的是可以放心写进报告里的结果。个人做这套仿真下来的总体感受是MFSK/FH系统的仿真难度不在单个模块而在模块之间的时序配合和参数的耦合关系。建议动手写代码前先把各模块的接口数据格式明确好——调制输出是什么形式、跳频叠加后变成什么形式、接收端先解跳还是先解调这些一旦定清楚代码写起来会顺很多。跳频通信在军事通信、应急通信、抗干扰数据链领域至今仍是不可替代的技术方案这套仿真做完之后你再去看实际的跳频电台指标会发现很多参数都能对应上。如果后续想做深一点还可以加进来信道编码、分集接收、自适应跳频这些方向都是在这套链路上做扩展。本文还有配套的精品资源点击获取
MFSK/FH跳频通信系统仿真：参数设计、误码率与工程实现
简介一份面向通信工程专业学生、研究人员及算法开发者的MATLAB仿真资源聚焦MFSK/FH高速跳频通信系统的建模与性能分析。该系统充分利用MFSK在有限带宽内传输更多信息、抗多径衰落能力强的特点以及FH快速跳变载波带来的抗干扰、抗截获和隐蔽性优势可服务于军事抗干扰通信、无线局域网等应用场景。仿真模型覆盖完整的信号发射、传输和接收流程支持自定义跳频图案、调制阶数与信道参数方便对比不同设置下的误码率表现也可深入研究跳频图案对通信安全性的提升作用。压缩包为zip格式共包含2个文件一个可直接运行的MATLAB脚本.m与一份说明文档.md整体仅4KB体积紧凑却便于快速查阅。目前已有92人学习下载。通过运行主程序可看到MFSK调制解调和FH跳频控制的实现细节结合说明文档可快速掌握参数调整与结果分析方法适合作为扩频通信课程设计、毕业设计或科研预研的速查参考。 接到这个选题时我第一反应是写MFSK/FH仿真的文章不少但大多数都停在贴代码跑通一版的层面真正把为什么要这么设计参数误码率曲线背后反映什么问题讲透的很少。这次我把完整做下来的经验和踩过的坑整理出来希望能给正在做通信仿真或者准备做毕业设计、课题预研的朋友一些参考。从仿真目标反推系统参数MFSK/FH链路该怎么拆很多人拿到高速跳频通信系统仿真这个题目第一件事就是打开MATLAB开始写代码。我劝你先停下来花半天时间把链路框图和数据流画清楚。因为MFSK和FH虽然原理都不复杂但合在一起之后采样率、符号速率、跳速、频率间隔这些参数会互相牵制参数设不合理后面跑出来的曲线自己都解释不了。为什么MFSK和FH是天生一对MFSK多进制频移键控是一种频率调制方式M进制意味着每个符号携带log2(M)比特信息用M个不同频率中的一个来承载。FH跳频则是把整个信号按照跳频图案在很宽的频带内搬移两者结合后信号在任一时刻既携带MFSK的调制频率信息又位于跳频图案指定的载波频率上。从抗干扰的角度看MFSK给系统带来了频率分集的效果某个频率点被干扰时其他频率点仍然可以正常工作FH则提供了处理增益干扰机不知道下一跳跳到哪个频率只能全频带压制代价极高。更关键的是MFSK是非相干检测友好的调制方式而跳频通信中接收机很难保持本振与发射机载波的严格相位相干MFSK天然适合这种场景。这套组合还有个工程上的好处实现门槛低。相比于需要精确相位跟踪的PSK类调制MFSK的调制解调用模拟器件或者DSP都能很方便实现这也是很多战术电台选择它的原因。确定仿真参数的关键考虑参数设计是整个仿真的起点也是我踩坑最多的地方。先给出一套完整可用的参数后面所有分析和代码都基于这套配置数据速率Rb20 kbps调制阶数M4即每个符号携带2比特符号速率Rs10 ksymbol/s频率间隔Δf25 kHz近似等于符号速率的设计值跳频速率2000 hop/s每跳承载5个符号跳频频率集规模Nhop64个频点频率间隔5 MHz总跳频带宽320 MHz采样率fs10 MHz单个跳频频点内看10 MHz采样足够容纳25 kHz的MFSK信号并留出余量这套参数下单跳内要处理5个MFSK符号跳频图案每毫秒更新一次。仿真时有两种建模方式一种是基带等效建模把所有频点搬移到基带处理速度快但不好观察频谱另一种是带通建模生成真实的射频波形直观但计算量大。我的建议是先用基带建模验证算法再抽出一小段做带通波形验证兼顾效率和可信度。MFSK调制的MATLAB实现为什么非相干解调更适合跳频MFSK的调制在MATLAB里其实就一句话根据输入的比特映射到对应的正弦频率。但解调方案的选择直接决定了接收机的结构和仿真代码的复杂度。频率映射与调制实现4FSK系统里每2个比特映射为一个符号频率。我采用的映射关系是00 - f0 100 kHz01 - f1 125 kHz10 - f2 150 kHz11 - f3 175 kHz注意这里的频率是相对跳频载波的偏置信号实际发射频率是跳频载波加上这个偏置。关键点在于频率间隔的取值。MFSK的非相干解调要求频率间隔不小于符号速率即Δf ≥ Rs否则相邻频率的频谱会发生严重重叠误码率会显著恶化。实践中我取Δf 2.5Rs留出一定保护带这样即使用矩形脉冲成型相邻频点间的串扰也可以忽略。调制代码核心部分如下% 参数定义 M 4; % 调制阶数 Rs 10000; % 符号速率 10 ksymbol/s T_symbol 1/Rs; % 符号周期 fs 10e6; % 采样率 10 MHz samples_per_symbol fs/Rs; % 每个符号的采样点数 % 频率映射表 freq_table [100e3, 125e3, 150e3, 175e3]; % 比特到符号索引的映射 Gray编码 bit_to_sym [0 1 3 2]; % 对输入比特流分组并调制 bits randi([0 1], 1, 10000); symbols_idx bit_to_sym(2*bits(1:2:end) bits(2:2:end) 1); % 生成MFSK调制信号单跳 t (0:samples_per_symbol-1)/fs; tx_signal []; for k 1:length(symbols_idx) tx_signal [tx_signal, cos(2*pi*freq_table(symbols_idx(k))*t)]; end这里用了一个细节比特到符号的映射采用了Gray编码映射00对应001对应111对应210对应3。Gray编码的好处是相邻符号之间只有1个比特差异当噪声导致解调器误判到相邻频率时只会错1个比特而不是2个。这个映射在高信噪比区域能让误码率直接减半属于投入产出比极高的小技巧。非相干解调的判定逻辑之所以弃用相干解调原因前面说过跳频系统里收发双方的本地振荡器频率难以完全同相加上快速跳频带来的频率切换瞬态做载波相位恢复非常痛苦。非相干解调则完全不在乎相位直接用能量检测。实现方式也简单接收端对每个符号周期做FFT然后找幅度最大的频率点对应的索引反查映射表得到比特。% 非相干解调 rx_symbols zeros(1, length(symbols_idx)); for k 1:length(symbols_idx) seg rx_signal((k-1)*samples_per_symbol1 : k*samples_per_symbol); % 加窗后FFT减少频谱泄漏 win hann(length(seg)); spectrum abs(fft(seg .* win)); % 只需要关注M个频点附近的幅度 [~, max_idx] max(spectrum(1:length(spectrum)/2)); % 找到对应的频率索引 [~, rx_freq_idx] min(abs(freq_table - max_idx*fs/length(seg))); rx_symbols(k) rx_freq_idx; end实际做的时候我加了汉宁窗后再做FFT因为矩形窗在频点非整数周期截断时会产生严重频谱泄漏干扰判决结果。跳频图案与频率合成器仿真里最容易拖慢速度的环节FH的仿真核心有两块伪随机序列的生成以及跳频图案到实际发射频率的映射。这两块做好了整个系统才有跳的灵魂不然就只是个MFSK仿真。伪随机序列与跳频图案生成跳频图案通常由一个伪随机序列控制频率合成器产生。仿真中用m序列最长线性反馈移位寄存器序列来生成其自相关性好、实现简单、具有伪随机性。m序列生成器的核心是反馈多项式。我用的本原多项式是x^8 x^6 x^5 x^4 1即8级移位寄存器周期255足够支撑64个频点的跳频图案而不重复。% 生成m序列 reg ones(1, 8); % 寄存器初始状态 seq_len 255; m_seq zeros(1, seq_len); for n 1:seq_len m_seq(n) reg(8); feedback xor(reg(8), xor(reg(6), xor(reg(5), reg(4)))); reg [feedback, reg(1:7)]; end % 映射到跳频频率每跳需要6比特控制64个频点 % 串并转换取6比特转10进制 hop_pattern zeros(1, floor(seq_len/6)); for h 1:length(hop_pattern) bits6 m_seq((h-1)*61 : h*6); hop_pattern(h) bi2de(bits6, left-msb); end % 通过异或加扰方式增加复杂度这里简化直接用映射 hop_freqs hop_freq_table(hop_pattern 1);注意一个细节m序列的周期255不能被6整除所以最后会浪费3个比特。实际系统中通常会用多个m序列组合生成更长的跳频图案仿真阶段不用追求极致的随机性能但至少保证跳频图案在仿真时长内不出现明显周期。跳频信号拼接的工程细节MFSK信号和跳频图案的时序配合是仿真中容易出问题的地方。跳频速率2000 hop/s意味着每0.5毫秒切换一次载波频率每跳内正好承载5个MFSK符号。发射信号的表达式是s(t)cos(2π(f_hop(n)f_mfsksymbol)·t)但这里有个隐含问题频率切换瞬间相位不连续。理想跳频系统允许相位不连续实际也能工作因为接收端用的是非相干解调。但在FFT分析跳频信号时会看到很宽的频谱扩散。仿真中我建议在每跳的边界处留出几个采样点的保护间隔更贴近工程中频率合成器建立时间。% 拼接带跳频的发射信号 tx_hop_signal []; hop_duration 0.5e-3; % 每跳0.5ms hop_samples hop_duration * fs; for h 1:length(hop_pattern) % 本跳的MFSK基带信号 start_sym (h-1)*5 1; end_sym h*5; hop_baseband tx_signal((start_sym-1)*samples_per_symbol1 : end_sym*samples_per_symbol); % 搬移到跳频频率 t_hop (0:hop_samples-1)/fs; hop_rf hop_baseband .* cos(2*pi*hop_freqs(h)*t_hop); tx_hop_signal [tx_hop_signal, hop_rf]; end这段代码跑起来你会发现一个问题循环拼接很慢数据长了之后内存消耗也大。更好的做法是预分配数组或者用矩阵化操作。跳频仿真尤其要注意性能优化我后面单独说。完整链路跑通后的性能分析误码率曲线能告诉我们什么系统搭完之后接下来才是重头戏把信道效应加进去跑不同信噪比下的误码率画出BER曲线分析系统性能。这一步做得好不好直接决定这篇仿真有没有说服力。高斯信道下的基准性能首先在AWGN加性高斯白噪声信道下测试系统的本征性能。这是所有性能分析的基准线任何一个跳频系统的误码率都不会优于同参数MFSK在AWGN下的理论值。4FSK非相干解调的理论误码率为P_e ∑_{k1}^{M-1} (-1)^{k1} C(M-1,k)/(k1) · exp(-k·Eb/(k1)·N0)我实测下来仿真曲线和理论值在低信噪比区域贴合得很好高信噪比区域也就是在Eb/N0大于14dB以后由于仿真数据量不足开始出现抖动。要解决这个问题至少要保证每个信噪比点统计到200个以上的错误比特低误码率区域建议统计到500个错误比特再停下来。干扰场景下的抗干扰增益跳频系统最核心的卖点是抗干扰。仿真中我设置了两种典型的干扰场景单频点连续干扰和宽带部分频带干扰。单频点干扰的模型是在某个固定频率上叠加强干扰信号。由于跳频图案的作用这个干扰只在系统跳到被干扰频点的那一刻产生影响。假设跳频有64个频点干扰占用1个频点那么系统只有1/64的时间被干扰到这就是处理增益。64个频点对应18dB的处理增益实际仿真跑下来的结果和这个理论值吻合。部分频带干扰模型更贴近实际情况干扰机不知道跳频图案但知道系统的工作频段于是压制其中一部分频带。比如干扰25%的频带那么系统约有四分之一的跳频点被干扰干扰频带内的符号会大量出错但依然有75%的频点能正常工作。仿真结果显示在干扰条件下MFSK/FH系统的误码率并不是线性恶化而是呈现明显的平层效应——即当干扰足够强时误码率不再随干扰功率显著变化而是被未被干扰的频点比例所决定。% 部分频带干扰仿真 interference_ratio 0.25; % 干扰25%频带 interfered_freqs randperm(length(hop_freq_table), round(interference_ratio * length(hop_freq_table))); % 在接收信号上叠加干扰 for h 1:length(hop_pattern) if ismember(hop_freqs(h), interfered_freqs) t_hop (0:hop_samples-1)/fs; rx_hop_signal((h-1)*hop_samples1 : h*hop_samples) ... rx_hop_signal((h-1)*hop_samples1 : h*hop_samples) ... A_jammer * cos(2*pi*hop_freqs(h)*t_hop rand*2*pi); end end同步误差对性能的影响很多人在仿真里默认收发双方完全同步这是最理想化的假设也是仿真结果和实测差距最大的地方。我专门做了一组对比当接收机的跳频图案相对发射机存在时间偏移时误码率如何变化。结果很有意思时间偏移低于一个MFSK符号周期的20%时性能几乎无损偏移达到半个符号周期误码率开始明显上升超过一个符号周期系统彻底无法工作。这主要是因为跳频同步误差会导致接收端解调窗口错位每个符号周期的FFT窗口内包含了两个不同符号的波形能量检测判错。这组数据说明了一个工程结论跳频系统的同步精度要求大约在符号周期的20%以内这也为接收机同步算法的设计提供了数字化指标。跳频仿真中那些只有跑过才懂的坑代码调通只算完成了一半真正耗时间的往往是那些看起来不起眼却又致命的细节。我把这次仿真中遇到的几个典型问题列出来希望能帮你少走弯路。频谱泄漏FFT频率分辨率不够怎么办MFSK解调用的是FFT找最大能量点但当符号周期不是FFT频率分辨率的整数倍时频谱能量会泄漏到相邻频点导致判决出错率高。我实测发现当采样率10MHz、FFT点数取1000时频率分辨率正好是10kHz而频点间隔是25kHz理论上足够区分。但矩形窗的旁瓣衰减只有13dB遇到噪声稍大就可能在旁瓣位置出现虚假峰值。解决方案有两个一是换窗函数汉宁窗旁瓣衰减可达31dB代价是主瓣变宽二是增加FFT点数用零填充到2048点能有效降低栅栏效应带来的幅度误差。实际我两个方法都用了误码率在高信噪比区域明显改善。跳变沿的毛刺忽略它会让误码率多一个数量级跳频频率切换不是瞬时完成的仿真中如果简单地在每跳边界处硬切换频率会在时域上产生一个相位突变造成很宽的频谱扩散。虽然在纯数学仿真中不影响MFSK解调但一旦加入带限信道滤波器毛刺就会变成真正的干扰。我在每跳边界的前后各留了20个采样点作为保护间隔将频率切换放在保护间隔内完成。这样每个跳频周期内有效数据只占96%损失可以忽略但信号频谱干净得多。仿真速度矩阵化是唯一出路跳频仿真最大的痛点是速度。如果坚持用for循环逐跳拼接信号仿真到100万比特数据量时会非常慢。我的做法是尽可能用MATLAB矩阵化操作代替循环。一个很实用的技巧是先生成完整的MFSK基带信号再用一个频率切换矩阵对信号进行分块搬移把逐跳乘载波变成向量点乘。此外如果只是看误码率趋势完全可以在基带等效域完成仿真也就是不显式生成射频频率而是给每个频点分配一个等效噪声和干扰功率速度可以快上数十倍。我建议先跑基带等效仿真摸清参数再跑带通仿真验证几个关键点。误码率统计的置信度问题还有一个容易忽略的问题误码率统计的样本量。很多人跑到感觉曲线光滑就停了但低误码率区域需要极大的数据量才能保证置信度。我的经验是每一个信噪比点错误比特数至少统计200个如果误码率低于1e-4统计到500个错误比特再停。这么做代价是仿真时间显著增加但换来的是可以放心写进报告里的结果。个人做这套仿真下来的总体感受是MFSK/FH系统的仿真难度不在单个模块而在模块之间的时序配合和参数的耦合关系。建议动手写代码前先把各模块的接口数据格式明确好——调制输出是什么形式、跳频叠加后变成什么形式、接收端先解跳还是先解调这些一旦定清楚代码写起来会顺很多。跳频通信在军事通信、应急通信、抗干扰数据链领域至今仍是不可替代的技术方案这套仿真做完之后你再去看实际的跳频电台指标会发现很多参数都能对应上。如果后续想做深一点还可以加进来信道编码、分集接收、自适应跳频这些方向都是在这套链路上做扩展。本文还有配套的精品资源点击获取 
