MATLAB扩频通信系统仿真：DSSS与FHSS从链路搭建到误码率分析

简介：这份资源是一篇基于MATLAB的扩频通信系统仿真毕业论文，面向通信工程、电子信息类专业的学生及需要完成相关毕业设计的学习者，帮助解决扩频通信理论建模与仿真验证的完整方案需求。压缩包内仅含1个PDF文件，大小约3.58MB，内容涵盖直接序列扩频、跳频扩频等技术的理论基础，以及利用SIMULINK搭建扩频、解扩与干扰处理模块的详细设计过程，并给出仿真参数设定与结果分析。论文重点研究了扩频增益与输出端信噪比、误码率之间的关系，验证了增大扩频增益可提升系统抗干扰能力的结论，同时包含摘要、关键词、正文、参考文献等规范章节，结构完整。目前已有267人学习下载，适合作为毕业设计写作模板、仿真建模参考与通信系统课程学习的辅助材料。

1. 扩频通信系统仿真：从毕业设计选题到能跑通的 MATLAB 工程

很多通信工程专业的同学拿到「基于 MATLAB 的扩频通信系统仿真」这个题目时，第一反应是去搜 matlab 下载安装教程，然后打开软件对着空白脚本发呆。这个题目看起来简单——扩频嘛，不就是把信号带宽展宽再解扩回来——但真正动手做的时候，你会发现涉及的知识点横跨信号生成、伪随机序列设计、调制解调、扩频解扩、信道建模、误码率分析、频谱分析等多个环节。它不是一个「跑一个函数就出图」的题目，而是一个需要你从系统层面理解信号流向、从工程层面搭建仿真链路的完整项目。

这篇文章面向正在做这个毕业设计、或者想用 MATLAB 把扩频通信系统跑通的同学。我会从系统架构讲起，把直接序列扩频（DSSS）和跳频扩频（FHSS）两条主流路线的仿真搭建过程拆开，给出可复现的代码框架和参数设置逻辑，最后把我在调试过程中踩过的坑和验证方法一并交代。读完你应该能独立搭出一套完整的扩频通信仿真系统，并且知道每个参数为什么这样设、出了问题去哪里找原因。

2. 扩频通信仿真的系统架构与 MATLAB 选型：先想清楚信号怎么流

2.1 扩频通信的核心机制与仿真映射关系

扩频通信的本质是用一个高速的伪随机码（PN 码）去调制一个低速的窄带信息信号，使得发射信号的带宽远大于信息本身所需的带宽。接收端用同一个 PN 码做相关解扩，把信号恢复成窄带，同时把非相关干扰和噪声扩散到整个带宽上。这个过程带来的处理增益 $G_p = B_{ss} / B_{info}$，直接决定了系统抗干扰能力的上限。

在 MATLAB 里做仿真，你需要把这个物理过程映射成离散时间信号处理链路。具体来说，信息比特流 → 符号映射 → 基带成形 → 与 PN 码相乘（扩频）→ 载波调制 → 信道叠加噪声/干扰 → 载波解调 → 与同步 PN 码相乘（解扩）→ 积分判决 → 误码统计。每一步在 MATLAB 里都有对应的函数或矩阵运算方式，关键在于采样率的设置要贯穿始终、各环节的维度要对齐。

我一般建议在仿真开始前先画一张信号流图（纸上画就行），标注每个节点的采样率、信号维度、数据类型。这张图后面会帮你省下大量调试时间，因为扩频仿真里最常见的翻车就是某个环节的采样率没对齐，导致解扩后星座图完全散掉。

2.2 DSSS 与 FHSS 两条路线的仿真复杂度对比

直接序列扩频（DSSS）和跳频扩频（FHSS）是这个题目最常选的两种方案。DSSS 的实现逻辑更直观：信息比特与 PN 码直接相乘，在基带完成扩频，然后统一上载波。FHSS 则需要在每个跳频周期内切换载波频率，涉及频率合成器的建模和跳频图案的生成。

从仿真复杂度来看，DSSS 的代码量更少、调试更直接，适合作为毕业论文的核心仿真方案。FHSS 的亮点在于跳频图案设计和频率切换时序，如果你论文里想体现「抗窄带干扰」或「抗截获」的特性分析，FHSS 的仿真结果会更有说服力。我的建议是：主仿真用 DSSS 跑通全链路和误码率曲线，再用 FHSS 做一个对比实验或补充分析，这样论文的仿真工作量和技术深度都够。

在 MATLAB 版本选择上，R2020b 之后的版本对通信工具箱（Communications Toolbox）的函数命名做了统一，比如comm.PNSequence系统对象在后续版本中逐步被函数式写法替代。如果你用的是 matlab 2023b 或更新版本，建议直接用函数式写法（如pnsequence配合矩阵操作），避免系统对象在新版本中的兼容性问题。另外，如果你在 Linux 下安装 MATLAB，注意检查通信工具箱是否完整安装，部分精简安装包会缺少berawgn、comm.ErrorRate等关键函数。

2.3 仿真参数体系：采样率、扩频因子与处理增益的联动设置

参数设置是扩频仿真的骨架。以下是我在多个项目中反复验证过的一套参数框架，你可以直接作为起点：

这里的关键联动关系是：扩频因子决定处理增益，处理增益决定误码率曲线的水平位置。如果你把 SF 从 31 改成 127，在相同 Eb/N0 下误码率会明显下降，但仿真时间也会成倍增加。我一般先用 SF=31 快速验证链路正确性，确认无误后再用 SF=63 或 127 跑最终的性能曲线。

采样率的设置有个容易忽略的点：扩频后的信号带宽是码片速率的两倍（BPSK 调制下），如果你的采样率不够，频谱会混叠，解扩后的眼图会糊成一团。经验法则是采样率至少取码片速率的 4 倍，如果后面要加成型滤波器（如根升余弦），最好取 8 倍。

3. 用 MATLAB 搭建 DSSS 仿真链路：从 PN 码生成到误码率曲线

3.1 m 序列与 Gold 序列的生成代码与优选逻辑

PN 码是扩频系统的核心。最常用的是 m 序列，由线性反馈移位寄存器（LFSR）生成。MATLAB 里可以用comm.PNSequence系统对象，也可以自己写 LFSR。我倾向于自己写，因为这样你能完全控制反馈抽头，也方便在论文里解释原理。

% m序列生成函数：基于LFSR，指定阶数和反馈抽头 function seq = generate_m_sequence(n, taps, seed) % n: 寄存器阶数，taps: 反馈抽头位置（如[6 1]表示第6和第1位异或） % seed: 初始状态（非零向量） N = 2^n - 1; % m序列周期 reg = seed(:)'; % 寄存器初始状态 seq = zeros(1, N); % 预分配输出 for k = 1:N seq(k) = reg(end); % 输出最右位 fb = mod(sum(reg(taps)), 2); % 反馈位：抽头异或 reg = [fb, reg(1:end-1)]; % 移位并注入反馈 end seq = 2*seq - 1; % 映射到 +1/-1，方便后续BPSK调制 end

这段代码的逻辑是：每次输出寄存器最后一位，同时根据抽头位置计算反馈位，将反馈位插入寄存器最前面，其余位右移。参数taps决定了序列的周期和自相关特性，常用的本原多项式抽头可以在通信教材的附录里查到。seed不能全零，否则序列会锁死。输出映射到 ±1 是为了后续直接与 BPSK 符号相乘完成扩频，省去额外的电平转换步骤。

Gold 序列适合多用户场景，由一对优选 m 序列异或生成。如果你论文里要做多用户扩频（CDMA 方向），Gold 序列的互相关特性比纯 m 序列更好，代码上只需要生成两个 m 序列然后按位异或、循环移位即可。

3.2 扩频调制与解扩的完整代码实现

下面是一条完整的 DSSS-BPSK 基带仿真链路，包含扩频、加噪、解扩和误码统计：

% DSSS-BPSK 基带仿真主脚本 clear; close all; % ---- 参数设置 ---- SF = 31; % 扩频因子 Nbits = 50000; % 信息比特数 EbN0_dB = 0:2:10; % Eb/N0 扫描范围 n = 5; % m序列阶数，2^5-1=31 taps = [5 2]; % 本原多项式抽头 seed = [1 0 1 0 1]; % 初始状态 % ---- 生成PN码 ---- pn = generate_m_sequence(n, taps, seed); % 长度31的±1序列 pn = repmat(pn, 1, ceil(Nbits*SF/length(pn))); pn = pn(1:Nbits*SF); % 裁剪到所需长度 % ---- 生成信息比特并扩频 ---- bits = randi([0 1], 1, Nbits); symbols = 2*bits - 1; % BPSK映射：0→-1, 1→+1 spread_signal = zeros(1, Nbits*SF); for k = 1:Nbits spread_signal((k-1)*SF+1 : k*SF) = symbols(k) * pn((k-1)*SF+1 : k*SF); end % ---- AWGN信道与误码率计算 ---- ber = zeros(size(EbN0_dB)); for idx = 1:length(EbN0_dB) EbN0 = 10^(EbN0_dB(idx)/10); % 扩频后每比特能量分散到SF个码片上，噪声方差按码片SNR计算 SNR_chip = EbN0 / SF; % 码片信噪比 noise_var = 1 / (2 * SNR_chip); % 实信号噪声方差 noise = sqrt(noise_var) * randn(1, length(spread_signal)); rx_signal = spread_signal + noise; % ---- 解扩：与同步PN码相乘并积分 ---- decision = zeros(1, Nbits); for k = 1:Nbits segment = rx_signal((k-1)*SF+1 : k*SF) .* pn((k-1)*SF+1 : k*SF); decision(k) = sum(segment); % 积分判决 end % ---- 误码统计 ---- decoded_bits = decision > 0; % 正判为1，负判为0 ber(idx) = sum(decoded_bits ~= bits) / Nbits; end % ---- 绘制误码率曲线 ---- semilogy(EbN0_dB, ber, 'bo-', 'LineWidth', 1.5); hold on; theory = berawgn(EbN0_dB, 'psk', 2, 'nondiff'); % BPSK理论误码率 semilogy(EbN0_dB, theory, 'r--', 'LineWidth', 1.5); xlabel('Eb/N0 (dB)'); ylabel('误码率 (BER)'); legend('DSSS仿真', 'BPSK理论'); grid on; title(['DSSS-BPSK误码率曲线 (SF=' num2str(SF) ')']);

这段代码的核心逻辑分三步：扩频时每个信息符号乘以一整段 PN 码，解扩时用同一段 PN 码做相关运算。注意噪声方差的计算——扩频后信号能量分散到 SF 个码片上，所以码片信噪比是 $E_b/N_0$ 除以 SF。这个细节如果搞错，你的误码率曲线会和理论值差出好几个数量级，这是扩频仿真里最经典的翻车点之一。

参数方面，SF决定了处理增益，也决定了仿真耗时。Nbits建议不低于 10000，否则高信噪比下的误码率统计会不稳定（可能一个误码都没有，曲线出现断点）。taps的选择影响 PN 码的自相关旁瓣，如果你要做同步算法仿真，选旁瓣更低的抽头组合会有帮助。

4. 跳频扩频仿真与干扰场景建模：让论文的仿真结果有对比度

4.1 FHSS 跳频图案生成与频率切换时序

FHSS 的仿真思路和 DSSS 不同：信息信号本身不扩频，而是载波频率按照跳频图案快速切换。每个跳频周期内传输一段窄带信号，接收端按相同图案同步切换本地振荡器频率完成解跳。

% FHSS跳频图案生成与信号合成 Nhop = 100; % 跳频时隙数 Nfreq = 8; % 可用频点数 hop_rate = 1000; % 跳频速率（跳/秒） symbol_rate = 100; % 每跳内符号速率 samples_per_hop = hop_rate / symbol_rate; % 生成跳频图案：这里用随机跳频，也可改为正交拉丁方等确定性图案 hop_pattern = randi([1 Nfreq], 1, Nhop); % 频率集合（归一化频率） freq_set = linspace(0.05, 0.4, Nfreq); % 合成跳频信号 fs = 8000; % 采样率 t_hop = (0:samples_per_hop-1) / fs; fhss_signal = []; for k = 1:Nhop fc = freq_set(hop_pattern(k)); tone = cos(2*pi*fc*t_hop); % 每跳内的载波 data = randi([0 1], 1, samples_per_hop); modulated = (2*data - 1) .* tone; % BPSK调制 fhss_signal = [fhss_signal, modulated]; end

跳频图案的设计是 FHSS 论文的加分项。随机跳频实现简单但抗干扰性能一般，如果你用混沌序列或 RS 码构造跳频图案，可以在论文里做性能对比。频率集合的选择要注意避开直流和奈奎斯特频率，否则频谱分析时会出现虚假分量。

4.2 窄带干扰与部分频带干扰下的性能对比

扩频系统最核心的优势就是抗干扰。在仿真里加入干扰源，观察误码率曲线的变化，是论文里体现「扩频通信系统」价值的必要环节。常见的干扰类型有两种：窄带干扰（单频或窄带噪声）和部分频带干扰（占据一定比例带宽的噪声）。

% 在DSSS仿真基础上加入窄带干扰 JSR_dB = 10; % 干信比 JSR = 10^(JSR_dB/10); interf_freq = 0.15; % 干扰频率（归一化） interf = sqrt(JSR * SF / 2) * cos(2*pi*interf_freq*(0:length(spread_signal)-1)); rx_signal = spread_signal + noise + interf;

加入干扰后重新跑误码率曲线，你会看到：没有扩频处理时，窄带干扰会直接恶化误码率；而扩频系统通过解扩把干扰功率扩散到整个带宽，误码率曲线只会略微抬升。这个对比图放在论文里非常有说服力。注意干扰功率的归一化方式要和信号功率对齐，否则干信比的物理意义就不对了。

部分频带干扰的建模稍微复杂一些：需要生成一段带限噪声，中心频率对准信号带宽的某个部分。可以用带通滤波器对白噪声滤波来实现，滤波器的带宽占比就是干扰带宽占比。仿真时建议扫描不同的干扰带宽占比（如 10%、30%、50%），画出误码率随干扰带宽变化的曲线，这样论文的分析深度会明显提升。

5. 扩频仿真避坑与排查：那些让误码率曲线跑偏的细节

5.1 误码率曲线不收敛或与理论值偏差过大

现象：仿真跑出来的误码率曲线在高 Eb/N0 处突然变平，或者整体比理论值高一个数量级。

原因：最常见的原因是噪声方差计算错误。扩频后每个码片的能量是 $E_b / SF$，如果你按信息比特的 $E_b/N_0$ 直接算噪声方差，噪声功率会被低估 SF 倍，导致误码率偏高。另一个原因是 PN 码没有对齐——扩频和解扩用的 PN 码序列起始位置不一致，相关运算结果接近零。

解决：检查噪声方差公式，确认用的是码片信噪比 $E_c/N_0 = (E_b/N_0) / SF$。PN 码对齐问题可以通过在解扩前打印前 10 个码片的相关值来排查，如果相关值在零附近随机波动，说明码没对齐。

5.2 扩频后频谱混叠导致解扩失败

现象：解扩后的星座图完全散开，误码率接近 0.5，但信噪比设得很高。

原因：采样率不足。扩频信号的带宽是码片速率的两倍，如果采样率低于码片速率的 2 倍，频谱会发生混叠，解扩时的相关运算失去意义。

解决：把采样率提高到码片速率的 4 倍以上。如果你用了根升余弦滤波器，采样率建议取 8 倍。在 MATLAB 里可以用pwelch函数检查扩频前后的功率谱密度，确认频谱没有混叠。

5.3 MATLAB 版本差异导致的函数兼容性问题

现象：代码在实验室电脑上能跑，换到自己笔记本上报错「未定义函数或变量」。

原因：MATLAB 不同版本对通信工具箱函数的命名和支持范围有差异。比如comm.PNSequence在 R2020b 之后逐步被pnsequence函数替代，berawgn的参数格式在旧版本中也有所不同。另外，如果你在 matlab 2023 中打开旧代码，中文注释可能出现乱码，这是因为编码格式不匹配。

解决：统一使用函数式写法而非系统对象，减少版本依赖。中文注释乱码问题可以在 MATLAB 首选项里把编码设为 UTF-8，或者用feature('DefaultCharacterSet', 'UTF-8')在脚本开头强制指定。如果实验室和自用电脑版本差距大，建议在论文里注明仿真所用的 MATLAB 版本号。

5.4 高信噪比下误码率为零导致曲线断点

现象：误码率曲线在高 Eb/N0 处突然中断，semilogy报错或曲线不连续。

原因：仿真比特数不够。当理论误码率低于 $1/N_{bits}$ 时，仿真中可能一个误码都不出现，计算出的误码率为零，对数坐标无法显示。

解决：要么增加信息比特数（比如从 10000 增加到 100000），要么在高信噪比区域改用误符号率或误帧率来统计。另一个技巧是在semilogy之前把零值替换为一个很小的非零值（如1/Nbits），这样曲线能连续显示，但要在论文里注明这是统计下限。

5.5 扩频因子与处理增益的混淆

现象：论文里写「处理增益 15 dB」，但仿真结果和理论计算对不上。

原因：处理增益的定义是 $G_p = 10\log_{10}(SF)$，但有些人误用了码片速率与信息速率的比值再取对数，如果中间经过了成形滤波或过采样，这个比值会被放大，导致处理增益被高估。

解决：处理增益只跟扩频因子有关，跟采样率无关。SF=31 对应约 14.9 dB，SF=63 对应约 18 dB，SF=127 对应约 21 dB。在论文里写处理增益时，直接引用扩频因子计算即可，不要引入采样率。

6. 从仿真到论文：用误码率曲线和频谱图讲清楚扩频的价值

做完整套仿真之后，你手里应该有几张关键图：DSSS 的误码率曲线（含理论对比）、FHSS 的跳频图案时频图、加入窄带干扰前后的误码率对比、扩频前后的功率谱密度对比。这几张图就是你论文仿真章节的核心素材。

我一般会这样组织论文的仿真部分：先给系统框图，说明信号流程和参数设置；然后放扩频前后的频谱对比，直观展示带宽展宽效果；接着给无干扰下的误码率曲线，验证链路正确性；最后给干扰场景下的性能对比，说明扩频系统的抗干扰增益。每一张图都要配一段分析，说清楚「为什么这条曲线是这样的」「和理论值的偏差来自哪里」「参数变化会带来什么影响」。

有一个技巧值得单独说：在画误码率曲线时，把不同扩频因子（SF=31、63、127）的曲线画在同一张图上，可以非常直观地展示处理增益对性能的提升。这张图放在论文里，审稿老师一眼就能看出你做了实质性的仿真工作，而不是随便跑了一个 demo。

另外，如果你想让仿真结果更有工程说服力，可以加入同步误差的影响分析。在实际系统中，PN 码同步不可能做到完全对齐，存在一定的相位偏差。你可以在解扩前给 PN 码加入一个可控的延迟（比如偏移 1 个、2 个码片），观察误码率的变化。这个实验能说明扩频系统对同步精度的要求，也是论文里体现「工程思维」的加分项。

我自己做这类仿真最大的教训是：不要一上来就追求复杂的信道模型和多用户场景，先把单用户 AWGN 信道下的链路跑通、误码率曲线和理论值对上，再逐步加干扰、加多径、加同步误差。每一步只改一个变量，确认结果符合预期后再进入下一步。这样即使出了问题，你也能快速定位是哪个环节引入的。希望帮到你。

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