简介：扩频通信系统实用仿真技术资源包，面向通信工程、电子信息类学生、科研人员及工程师，系统覆盖直接序列扩频（DSSS）与跳频扩频（FHSS）的核心原理，并围绕码序列生成、调制解调、信道与干扰建模、同步机制、性能评估等关键环节提供可运行的仿真支撑。压缩包内共68个文件，其中37个mat数据文件用于存储波形与仿真结果，18个m脚本实现调制、解扩及误码率分析等算法逻辑，13个mdl模型展示Simulink系统框图，总计14.81MB，轻量易用，便于本地复现与二次开发。资源可帮助读者从基础概念验证进阶到完整系统级仿真，结合具体代码与模型快速理解扩频通信的抗干扰与多址接入特性，适合课程设计、毕业设计及实际项目预研。目前已有420人学习下载，是快速上手扩频仿真实践的实用参考资料。

1. 扩频通信系统仿真：为什么直接序列扩频和跳频扩频的代码不能通用

扩频通信的核心是用一个高速伪随机码把窄带信号在频域“摊开”，换取抗干扰、抗截获和码分多址能力。我拆这份“扩频通信系统实用仿真技术”资源时发现，它的价值不在公式推导，而在把直接序列扩频和跳频扩频两条链路的仿真代码分开组织，从m序列生成一路做到误码率评估。直接序列扩频靠码片速率扩展带宽，跳频扩频靠载波在预设频点上快速跳变躲避干扰，两者的发射机结构、同步方案和信道建模方式几乎完全不同，拿一套代码通吃必然翻车。对正在做通信课程设计，或者需要用MATLAB搭物理层验证平台的工程师来说，这份资源能省掉大部分搭框架的时间。下面按两条链路分别拆解，重点放在能跑的代码和参数边界上。

2. 直接序列扩频仿真链路：从m序列发生器到BPSK解扩与判决

2.1 为什么DSSS仿真第一步必须做PN码发生器

DSSS的扩频操作本质上是把每个信息符号乘以一个高速伪随机码。这个码的码片速率远高于信息速率，信号带宽因此被扩展到码片速率量级。仿真中最常用的PN码是m序列，由线性反馈移位寄存器生成，序列长度由寄存器级数n决定，取2^n - 1。m序列有良好的自相关性，这个特性决定了接收端能通过相关运算从噪声中把信号提出来。生成m序列时必须选对反馈抽头，抽头对应本原多项式，选错就退化成普通随机序列，自相关旁瓣抬高，解扩性能直接恶化。

手写LFSR比直接用MATLAB通信工具箱的comm.PNSequence更直观，因为你能看到每个时钟沿的移位和反馈过程。初版仿真我建议手写，调试时才好查状态。

function seq = generate_mseq(n, taps) % n: 寄存器级数，序列长度 = 2^n - 1 % taps: 反馈抽头位置，例如 [7 6] 对应本原多项式 x^7 + x^6 + 1 L = 2^n - 1; reg = ones(1, n); % 初始状态全部置1，不能全零 seq = zeros(1, L); for k = 1:L seq(k) = reg(n); % 输出最高位 fb = mod(sum(reg(taps)), 2); % 反馈位 = 抽头位异或 reg = [fb reg(1:n-1)]; % 左移一位 end seq = 2*seq - 1; % 0/1转成+1/-1，便于与BPSK符号直接相乘 end

这段代码的逻辑：每个时钟周期先取出最高位作为当前码片，然后把指定抽头位置的值做模2求和得到反馈位，最后移位寄存器左移一位，反馈位补到最低位。把输出从0/1映射到+1/-1，是为了后续扩频乘法可以直接用符号运算。taps参数是生成m序列的关键，常用本原多项式如下表：

仿真中n=7最常用，码长127对应扩频增益10log10(127)约21dB。这个增益直接决定了系统能承受多大的干扰和噪声。如果只是跑通流程，n=3的7位码就够，但画BER曲线时码长太短，扩频增益小，和理论BPSK的差异看不出明显规律。实际项目里选码长还要考虑多径时延扩展，时延扩展超过码片周期时，码间干扰会压过扩频增益。

2.2 扩频与BPSK调制：信息符号和PN码片的逐点乘法

发射端流程是信息比特先映射成BPSK符号，0映射为-1、1映射为+1，然后每个符号与整个PN码序列逐点相乘。每来一个信息比特，就输出一个完整的PN码周期，等效于把符号速率从Rb提升到Rb乘N。N就是PN码长度，也就是扩频因子。这个乘法的物理意义是，信号的频谱被PN码的码片速率展宽，功率密度随之下降，信号变得“淹没”在噪声和干扰里。

data = randi([0 1], 1000, 1); % 1000个随机信息比特 data_bpsk = 2*data - 1; % BPSK映射：0->-1, 1->+1 pn = generate_mseq(7, [7 6]); % 127位m序列 N = length(pn); % 扩频因子 spread = zeros(length(data_bpsk)*N, 1); for i = 1:length(data_bpsk) spread((i-1)*N+1 : i*N) = data_bpsk(i) * pn; end

逻辑说明：循环体内把第i个信息符号乘以整个PN码序列，结果写入输出向量中以i为中心的N个位置。乘法在+1/-1域进行，符号相同时输出码片保持原样，符号相反时码片整体翻转。这正是DSSS“用PN码给信号打烙印”的过程。

参数说明：这个实现里每个码片只取一个采样点，是最简形式。后续如果加脉冲成形或匹配滤波器，要先确定过采样倍数，比如每个码片采8个点。我的习惯是先完成扩频，再对码片序列做上采样。先上采样再扩频会让PN码片边沿叠加得很混乱，出了问题不好排查。

如果要在Simulink里搭同样链路，可以用“PN Sequence Generator”模块配合码片时钟，再用“Multiply”模块完成扩频。手写代码的优势是能精确控制每一级信号维度，Simulink的优势是模块化修改方便。初学阶段先在MATLAB脚本里跑通，再迁移到Simulink，排查错位问题会容易很多。

2.3 理想同步下的解扩与判决：先验证扩频增益本身

解扩是扩频的逆过程。接收端用本地PN码与接收信号逐点相乘再累加。当本地码片序列与接收信号中的PN码完全对齐时，累加结果等于信息符号能量加上被平摊的噪声；未对齐时累加结果接近零。这个差距就是解扩能工作、窃听者无法解调的根本原因。

received = spread + 0.08*randn(size(spread)); % 加高斯噪声 despread_symbol = zeros(1000, 1); for i = 1:1000 seg = received((i-1)*N+1 : i*N); % 截出当前符号的N个码片 corr = sum(seg .* pn); % 与本地PN码相关累加 despread_symbol(i) = corr; % 保存相关值供判决 end received_bits = double(despread_symbol > 0); % 大于0判1 ber = sum(received_bits ~= data) / 1000; % 误码率

逻辑说明：seg .* pn把接收码片逐个和PN码相乘，噪声是随机量，相乘后仍是随机量；信号分量由于相位对齐，累加后幅度为N。相关后的信噪比相比码片级信噪比提升了N倍，这就是扩频增益的来源。这里先做理想同步，即假设本地PN码在时间和码相位上完全对齐。这是验证扩频机制的第一步，同步模块后面再补。

参数说明：噪声标准差0.08对应的码片级信噪比约10log10(1/0.08^2)约22dB，解扩后符号级信噪比又提升21dB，理论误码率应非常低。若要画完整BER曲线，就要在循环里遍历Eb/N0，每个信噪比下至少跑几百个符号取平均，否则曲线抖动大。系统性的BER扫描框架放到第4章展开。

2.4 多址接入与调制选项：DSSS的CDMA能力怎么体现

DSSS的多址能力来自PN码族。不同用户分配不同的PN码，接收端用本用户的码解扩，其他用户的码因为互相关性低而被抑制。仿真多用户时，只需要生成两组不同抽头或不同码相位的PN码，把两路扩频信号相加后再解扩，观察目标用户误码率随干扰用户功率的变化。这就是CDMA的基本模型。

pn_user1 = generate_mseq(7, [7 6]); % 用户1的码 pn_user2 = generate_mseq(7, [7 3]); % 用户2的码，抽头不同 spread2 = zeros(size(spread)); for i = 1:1000 spread2((i-1)*N+1 : i*N) = data_bpsk(i) * pn_user2; end rx_cdma = spread + 0.3*spread2; % 两用户叠加，干扰功率比为0.3 despread_user1 = sum(rx_cdma(1:N) .* pn_user1); % 第一符号解扩

逻辑说明：用户2的扩频信号对用户1来说相当于宽谱干扰，解扩后残留的能量取决于两个PN码的互相关。抽头差异越大，互相关越低，多址干扰越弱。这个示例把多址干扰简化成功率加权叠加，适合观察码间隔离度的量级。

调制方面，摘要里提到DSSS可以和GMSK、FSK结合。仿真时扩频和调制是两级串接，扩频后的码片序列作为调制器输入。BPSK最简单，解扩前不需要额外处理载波相位；QPSK频谱效率更高，但解扩前要区分I/Q两路。GMSK的相位连续特性会和码片序列的矩形边沿相互作用，需要把调制器放在扩频之后，并保证码片速率与GMSK的BT参数匹配。第一版建议先用BPSK跑通，再替换调制模块。

3. 跳频扩频仿真：跳频图案参数、同步窗口设置与多径分集的三步仿真法

3.1 跳频三要素：频率集合、跳频速率、跳频图案

FHSS与DSSS的根本区别在于信号瞬时带宽不扩，载波频率在预设的频带内按伪随机图案跳变。接收端必须知道同一时刻的载波频率才能正确解调。仿真的三个核心参数是频率集合、跳频速率和跳频图案。频率集合定义跳变范围，通常等间隔分布；跳频速率决定每个频点的驻留时间；跳频图案决定跳变顺序，真实系统由PN码发生器驱动频率合成器实现。

fc_base = 900e6; % 起始频率 freq_spacing = 25e3; % 频点间隔，需大于信号带宽 num_freqs = 16; % 频点数 freq_set = fc_base + (0:num_freqs-1)*freq_spacing; % 用PN码控制跳变序列，这里先用均匀随机近似 hop_indices = randi(num_freqs, 10000, 1); tx_freq = freq_set(hop_indices);

逻辑说明：randi生成的随机索引等效于伪随机跳频图案，后续可以替换成由m序列发生器驱动的索引选择器。仿真初版用均匀随机就能跑通流程，重点是先把频率集合的定义和切换机制理清楚。

参数说明：freq_spacing必须大于信号单边带宽，否则相邻频点频谱重叠，解调时邻道干扰严重。25kHz间隔适合带宽小于25kHz的窄带信号。跳频速率在慢跳频系统中一般取符号速率的几分之一，比如符号速率10kbps、每频点驻留100个符号，跳频速率就是100跳/秒。

3.2 跳频同步：理想同步与同步头两种方案

跳频同步比DSSS码同步更直观也更麻烦。接收端必须在每个跳频时刻准确切换到发射端当前频点，差一个频点整个符号全错。工程上常用同步头方案：每帧开头发送一段已知频点的同步序列，接收端在同步序列期间扫描频率集合，找到相关峰后锁定跳频时钟。仿真时可以先做理想同步，把收发端的hop_indices设为完全相同，验证解调链路本身没有错。

sync_frame = [1 1 1 -1 -1 1 -1]; % 已知同步序列（BPSK） corr_out = zeros(num_freqs, 1); for f = 1:num_freqs % 用第f个频点下变频后与sync_frame相关 % 相关值最大的f就是当前发射频点 corr_out(f) = abs(sum(rx_downconverted(f).*sync_frame)); end [~, detected_freq] = max(corr_out);

逻辑说明：接收端对每个候选频点做一次下变频和相关运算，同步序列与本地参考的相关峰只在真实频点处出现，选择相关值最大的频点作为当前跳频位置。这种方法在低信噪比下需要更长的同步序列来保证检测概率，同步序列越长，帧开销越大。

参数说明：频点扫描需要接收机能快速切换本振，扫描时间受限于接收机架构。仿真中理想同步没有这部分开销，但要知道真实系统里跳频驻留时间必须大于同步捕获时间加上频率合成器锁定时间。慢跳频系统里通常预留10~50us的保护间隔，保护间隔内不传输有效数据。

3.3 慢跳频在频率选择性衰落下的分集效果

跳频抗多径的机理是频率分集。当两个跳频频点的间隔大于信道相干带宽时，它们经历的衰落相互独立。某个频点深衰落时，跳到另一个频点就能避开。仿真时给每个频点生成独立的信道增益系数，观察多次跳频平均后的误码率。

ber_hops = zeros(5000, 1); for hop = 1:5000 ch_gain = (randn + 1i*randn)/sqrt(2); % 瑞利衰落，每个频点独立 rx_sig = ch_gain * tx_sig + ... 10^(-snr_dB/20) * (randn(size(tx_sig)) + 1i*randn(size(tx_sig)))/sqrt(2); % 解调并判决当前跳频周期的符号 % ber_hops(hop) = 当前跳频周期的误码数 / 符号数 end avg_ber = sum(ber_hops) / 5000;

逻辑说明：每个跳频周期生成一个独立信道增益，等效于该频点经历一次独立的瑞利衰落实现。多次跳频平均后的误码率体现了频率分集的效果。如果所有频点共享同一个信道增益，就没有分集增益，系统退化为窄带瑞利衰落。

参数说明：snr_dB这里表示符号信噪比。瑞利信道下BPSK误码率比AWGN信道差很多，但有了跳频分集后，高信噪比区的BER会明显改善。要量化这个增益，可以对比单频点下5000次独立衰落平均的BER和跳频下5000个频点平均的BER，差异就是分集带来的收益。

4. 信道与干扰模型：AWGN加噪标定、窄带干扰与多径衰落的建模细节

4.1 AWGN加噪的标定：Eb/N0、码片能量与噪声方差的换算

扩频仿真的噪声标定是翻车重灾区。最稳妥的做法是先固定信号能量，再从目标Eb/N0反推噪声功率。BPSK信号的符号能量Eb在归一化后为1，每个码片能量为Ec=Eb/N。接收端的噪声功率谱密度N0由Eb/N0决定，复基带噪声的实部和虚部方差各为N0/2。

EbNo_dB = 6; EbNo = 10^(EbNo_dB/10); Eb = 1; N0 = Eb / EbNo; noise_var = N0 / 2; % 复噪声每维的方差 noise = sqrt(noise_var) * (randn(size(spread)) + 1i*randn(size(spread))); rx_with_noise = spread + noise;

逻辑说明：这里的关键是N0是双边功率谱密度，复基带噪声的总功率就是N0，每维度各占一半。如果误用sqrt(N0)作为每维标准差，噪声功率会翻倍，BER曲线右移约3dB。建议每改一次信噪比，都把信号和噪声的实际功率打印出来核对。

参数说明：如果信号在扩频后还经过了脉冲成形或功率归一化，Eb就不再等于1。更通用的做法是先计算spread的均方能量P_sig，再令噪声方差为P_sig乘以N再除以2倍的Eb/N0，即使信号形式变了也不容易错。

4.2 窄带干扰：单音干扰在解扩后的功率摊平效应

窄带干扰是DSSS最直观的抗干扰对象。仿真时在信号上叠加强度远高于信号功率的单音正弦波，经过解扩和相关运算后，单音干扰的能量被PN码摊开到整个扩频带宽，等效干扰功率下降N倍。这就是DSSS抗窄带干扰的数学体现。

fs = 127e3; % 采样率，对应127个码片每毫秒 t = (0:length(spread)-1)' / fs; sjr_dB = -20; % 干扰功率比信号高20dB interf_power = 10^(sjr_dB/10) * var(spread); f_interf = 30e3; % 干扰频率 interf = sqrt(2*interf_power) * sin(2*pi*f_interf*t); rx_interf = spread + interf + noise;

逻辑说明：解扩后，单音干扰被PN码序列调制，频谱从单根谱线变成宽谱，相关累加器对宽谱干扰只有噪声级的响应。SJR=-20dB、扩频增益21dB时，解扩后等效信干比约为1dB，信号仍能解调。把干扰频率放在码片速率附近，让干扰落在信号带内，能更明显观察到BER变化。

参数说明：干扰频率f_interf和采样率fs不能成整数倍关系。如果fs恰好是f_interf的整数倍，采样点固定在某一个相位上，正弦波退化成常数，干扰模型失真。fs=127e3、f_interf=30e3就避开了这个问题。

4.3 多径衰落：抽头延迟线模型和频率选择性信道的出现

多径对DSSS的影响表现为码片级的时间分散。当多径时延超过一个码片周期时，不同路径在接收端叠加出频率选择性衰落。抽头延迟线是仿真最常用的多径模型，每个抽头的延迟和增益由功率时延谱决定。

path_delays = [0 2 5]; % 以码片为单位的时延 path_gains = [1 0.5 0.2]; % 各径增益，主径为1 rx_mp = zeros(size(spread)); for p = 1:length(path_delays) d = path_delays(p); scaled = path_gains(p) * [zeros(d,1); spread(1:end-d)]; rx_mp = rx_mp + scaled; end

逻辑说明：每条径就是原码片序列延迟d个码片再乘上增益。第二径和第三径在主径之后到达叠加，产生符号间干扰。DSSS接收端用相关解扩时，主径相关峰两侧会出现副峰，副峰高度由多径增益决定。当延迟超过PN码周期时，前后码周期的自相关旁瓣会干扰主峰，单用户仿真一般让总时延小于码长。

参数说明：path_delays建议取小于PN码周期的值，比如码长127时取0、2、5。若需要研究超延迟扩展场景，PN码长要相应加长，否则扩频增益不足以分辨多径。真实系统里还可以用RAKE接收机把多径能量合并起来，仿真时对每条径单独解扩再做最大比合并，这是直接序列扩频在移动通信里的经典做法。

5. 扩频仿真避坑：码同步、噪声标定、增益计算的五个高频雷区

5.1 解扩后BER不降反升

现象：同样的信噪比下，解扩后的误码率比未扩频还高，或者BER曲线完全乱掉。

原因：本地PN码与接收PN码的相位没有对齐。相差一个码片，相关值会迅速掉到峰值以下，信号能量无法累积。

解决：先做相关峰搜索。把接收信号与本地PN码的循环移位版本逐一相关，找到最大相关峰对应的移位量，再按这个移位量对齐本地码。仿真里可以先人为设置一个已知延迟，验证搜索算法正确后再接入噪声信道。我一般会写一个函数输出相关峰位置和峰值幅度，跑无噪声信号时峰值应该非常尖锐。

5.2 扩频增益与理论值差20dB

现象：理论增益21dB，仿真结果只有1dB，或者完全没有增益。

原因：把过采样倍数当成了扩频因子。有些仿真为了波形精度，每个码片采8个点，误把8倍上采样也算进带宽扩展，导致N算错。

解决：扩频因子的定义只有一个——信息符号持续时间除以码片持续时间。上采样只是仿真精度手段，不改变物理带宽关系。代码里把N固定写成length(pn)，凡是涉及采样率的环节单独用一个变量oversample表示，两条线分开，就不容易混。出现这个坑的时候，检查代码里凡是出现N的地方是不是真的都指向同一个变量。

5.3 AWGN信道BER曲线右移3dB

现象：理论BPSK曲线在6dB处误码率约10^-3，仿真却在9dB处才到这个值，整条曲线往右挪了3dB。

原因：复噪声的方差标定错了。很多人用sqrt(N0)作为复噪声每维标准差，实际每维应该是sqrt(N0/2)。两维合起来才是总功率N0，多算的3dB就是噪声功率翻倍。

解决：统一用复基带模型。实部虚部各生成方差N0/2的独立高斯噪声，再合成为复噪声。在噪声相加前后分别计算信号功率和噪声功率，用打印值验证比值等于目标Eb/N0。这个检查每次都要做，改采样率或改信号功率之后尤其容易复发。

5.4 跳频仿真每换频点后连续误码

现象：某个跳频周期开头连续几十个符号全错，其余时间段正常。

原因：仿真里频点瞬时切换，没有建模频率合成器的锁定时间。真实系统中切换瞬间新频点的载波相位和幅度尚未稳定，不能承载数据。

解决：在跳频图案中每个频点驻留期的开头预留一段保护间隔，长度大于锁定时间。仿真时把保护间隔对应的采样点从数据中删掉，接收端也不对这些采样点做判决。锁定时间取10us，符号速率10ksps时对应0.1个符号，但在高速系统中可能吃掉好几个符号，设计跳频速率时要把这个开销算进去。

5.5 衰落信道BER曲线一直抖动不收敛

现象：加大循环次数，BER还在10^-2和10^-3之间反复横跳。

原因：Monte Carlo次数不够，或者每次循环只用一个信道实现。瑞利衰落的深衰落事件发生概率低，少数深衰落样本会剧烈拉高误码率。

解决：把循环次数提高到1万以上，并且每次循环都重新生成信道系数。观察BER随循环次数的收敛曲线，如果到末尾仍然跳跃，继续加次数。还可以固定随机种子分块运行，对比不同种子结果是否接近。收敛以后再做参数调整，否则改一个参数看不出真实效果。

6. 进阶验证：用BER曲线反推扩频增益并校准整个仿真链路

仿真的最终输出不是一张BER图，而是一组能证明“链路机理正确”的对比数据。我最常用的校准方法是把三条曲线叠在同一张图上：理想同步DSSS下的BER、理论BPSK BER、FHSS在AWGN下的BER。横轴统一用解扩后的Eb/N0。理论上DSSS曲线应该和BPSK曲线重合，因为解扩恢复的是符号能量，扩频不改变符号级信噪比；FHSS在无干扰时和BPSK也应接近，因为跳频本身不提供信噪比增益，只提供频率分集。如果DSSS曲线明显左移或右移，说明解扩前后的能量标定不一致；如果FHSS曲线偏离，多半是频点切换处符号丢失没处理好。

验证方法还有一个微妙之处：横轴如果改用码片信噪比Ec/N0，DSSS曲线会比BPSK左移10log10(N)dB，这个左移量就是扩频增益。所以同一套代码只要改一下横轴定义，就能同时验证“扩频增益存在”和“系统链路正确”两个结论。初学者最容易把这两个坐标系混在一起解释，最后得到的结论自己都圆不回来。建议在画图脚本里把横轴变量名写得非常明确，要么是ebno_dB，要么是ecno_dB，一个字符都不能含糊。

从那以后，我每次拆完一个扩频仿真包，都强制跑一遍“无干扰通道零误码、AWGN曲线与理论对齐、干扰场景增益可解释”三步检查，过了这三关，后面加多径、加多用户才有意义。这个习惯帮我少排了无数个定位不明的错。这个资源里的DSSS和FHSS脚本已经把这三步检查所需的链路搭好了，你下载后直接替换成自己的速率和码长参数就能复现，希望帮到你。

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