简介这是一套面向通信工程专业学生与科研人员的扩频通信系统实用仿真资源围绕直接序列扩频DSSS与跳频扩频FHSS两大核心体制涵盖PN码生成、GMSK/FSK调制、AWGN与多径信道建模、同步机制及BER性能评估等关键环节。压缩包共68个文件其中包含37个MAT数据文件.mat、18个MATLAB脚本.m和13个Simulink模型.mdl文件类型覆盖从底层仿真代码到可视化模型的不同需求整体大小14.81MB便于按模块调用学习。目前已有420人学习浏览适合具备一定通信原理基础、希望将理论落实到代码与仿真实验的读者。通过学习这些文件可以快速掌握扩频系统的模型搭建与参数调整方法理解抗干扰和多址接入的实现思路为课程设计或毕业设计提供可直接运行和二次开发的基础。1. 直接序列扩频仿真为什么说它是理解扩频通信最快的路径直接序列扩频DSSS的仿真表面上看就是“BPSK调制、加噪声、数误码”三段式但真正把扩频码、码片速率、同步和抗干扰串起来之后这套链路会暴露出一堆书上不写的问题同样的信噪比为什么扩频前后误码率差出好几个数量级PN码相位没对齐时为什么解出来的数据是全错的处理增益到底从哪里来又在哪里消失。这篇笔记按我平时搭扩频通信实用仿真链路的顺序展开先交代码的选择和基带模型再给一套能直接跑的MATLAB实现最后把我在参数匹配、同步、噪声换算上踩过的坑和排查方法写透。适合刚开始做扩频仿真的学生也适合需要在项目预案里快速评估直扩方案可行性的工程师。2. 直扩系统的基带等效模型从扩频码选型到参数匹配2.1 扩频码的选择m序列、Gold码与Walsh码在仿真中的差异直接序列扩频的仿真起点不是调制器而是扩频码。扩频码决定了序列的周期、自相关特性和多用户下的互相关水平这三者直接影响同步算法能不能抓到峰值、多用户仿真里会不会出现“自己干扰自己”。m序列由线性反馈移位寄存器LFSR生成周期为2^n-1自相关函数接近理想冲激旁瓣恒定为-1/周期这个特性让它在单用户测距和抗干扰仿真里非常好用。Gold码由一对优选m序列模2加生成周期相同互相关值被限制在三值集合内峰值有界所以两三个用户同频段跑的时候Gold码更接近工程习惯。Walsh码完全正交长度必须是2的幂同步条件下用户间干扰为零但它对相位偏差极其敏感一旦失去正交性互相关立刻恶化。从仿真选型的角度我一般遵循这样的优先级单用户验证抗干扰和误码率用长度31或63的m序列就够计算量小相关峰尖锐定位问题直观需要模拟多用户共信道换成Gold码只有做严格同步的正交信道模型比如前向链路结构才考虑Walsh码。三种码在仿真里的表现差异可以先用一张表记住码型生成方式典型周期自相关旁瓣互相关特征适合场景m序列LFSR本原多项式2^n-1-1/周期峰值尖锐可能较大需避开劣质多项式单用户抗干扰、同步捕获、测距Gold码优选m序列模2加2^n-1与m序列相当三值有界峰均比受控多用户CDMA仿真Walsh码Hadamard矩阵行向量2^n非冲激型同步时为零失步时恶化严格同步的正交信道选码时还有一个仿真层面的细节m序列的反馈多项式不是随便选。用同一个n值但不同的本原多项式生成序列的周期相同但自相关旁瓣和互相关峰值会不一样。仿真里如果要对比两种码的性能务必把多项式都列出来否则换了多项式结果变了会误以为是信道的问题。2.2 搭建基带等效模型发送端、信道与接收端的最小仿真框架很多第一次搭直扩仿真的人一上来就把载波频率、上变频、下变频全部写进代码结果定位一个bug要在频谱图和时域波形里翻半天。我的习惯是先把基带等效模型跑通所有信号都写成一维采样序列省去载波只在最后需要验证射频端行为时才把载波加回来。基带等效的好处是仿真速度更快Eb/N0可以直接用符号功率与噪声功率计算而且只要信道是线性的误码率结果与带载波版本完全一致。基带等效模型的最小框架包含三段发送端将数据比特按扩频因子sf扩展成码片序列再与PN码逐码片相乘完成扩频信道部分叠加复高斯白噪声可选加一个窄带干扰项接收端用本地PN码与接收序列做相关累加按比特周期判决统计错位数。这个框架看起来简单但它直接决定了后面所有参数怎么设。关键点在于采样率与码片速率的关系。最简配置是采样率等于码片速率一个码片恰好一个采样点解扩时不需要额外做抽取如果采样率是码片速率的整数倍解扩前就需要先做匹配滤波和抽取代码复杂度明显上升。初版仿真我强烈建议先用“每码片一个采样”的配置把功能验证完整后再升级到过采样模型。另外基带模型里要顺手把扰码加上——如果数据里出现长串的连续0或连续1PN码相乘后虽然不会改变统计特性但后续做同步捕获时相关峰会受到数据极性的调制导致某些相位上的峰被压平。很多仿真翻车就翻在忽略了这个小细节。2.3 关键仿真参数设置扩频因子、码速率、采样率的匹配关系扩频因子sf是直扩仿真最核心的参数每个数据比特对应sf个码片处理增益的理论值就是10lg(sf)。仿真中还要分清三个速率数据速率Rb、码片速率Rc和采样率fs。每码片一个采样时fs等于Rc一个数据比特在序列里占据sf个采样点误码率统计要以比特为单位逐段处理不能按采样点直接判。我设置参数的顺序是这样先定数据速率Rb比如1kbps按抗干扰需求或带宽约束定sf比如sf31那Rc就是31kchips/s初版仿真的采样率直接等于Rc后面要加成型滤波器的时候再把采样率提升到Rc的4倍或8倍。这里最容易踩的坑是把扩频因子加到数据速率上——实际上sf只改变码片速率不改变数据速率带宽扩展的倍数才是sf本身。给一组具体对照参数数值说明数据速率Rb1 kbps信息比特速率扩频因子sf31每比特码片数码片速率Rc31 kchips/sRb × sf处理增益14.9 dB10lg31采样率fs31 kHz每码片1采样初版用过采样率124 kHz4倍过采样成型滤波用噪声功率的设置同样依赖这套参数。在基带等效仿真里Eb等于单比特扩频序列的能量也就是sf因为每个码片幅度为1N0由Eb和Eb/N0反推实噪声方差是N0/2。后面章节的代码里会直接体现这套换算这里先记住结论噪声参数的设置必须以Eb/N0为入口而不是拍脑袋写一个SNR否则仿真结果和理论曲线对不上时根本不知道错在哪。3. 用MATLAB把直扩系统跑通发送、信道与解扩的完整实现3.1 发射端实现扩频码生成与数据调制先写一个m序列生成函数。m序列的本质是LFSR本原多项式决定反馈抽头位置初态不能全零否则序列卡死在零态不输出。下面是我最常用的一个版本function code mseq(taps, n) % 生成双极性m序列输出±1 % taps: 反馈抽头位置如 [5 2] 表示本原多项式 x^5 x^2 1 % n: 移位寄存器级数序列周期为 2^n - 1 reg ones(1, n); % 初态全1避免全零状态锁死 len 2^n - 1; code zeros(1, len); for k 1:len code(k) reg(end); % 最高位作为输出 fb mod(sum(reg(taps)), 2); % 反馈位抽头位置模2和 reg [fb, reg(1:end-1)]; % 寄存器右移一位 end code 2 * code - 1; % 0/1转成-1/1方便与数据符号相乘 end逻辑说明LFSR每一步先输出当前最高位再计算反馈位然后整体右移。反馈抽头必须与本原多项式对应比如n5时常用x^5x^21对应taps[5 2]换成其他本原多项式后序列周期仍是31但自相关旁瓣可能出现非最小值的退化所以仿真里固定抽头有利于结果复现。参数说明里最重要的一条是初态全零状态会让LFSR永远输出0初态置全1是标准的规避方式。拿到PN码之后发射端的扩频与数据调制可以合并写sf 31; % 扩频因子与m序列周期一致 Nbits 2000; % 仿真数据比特数 bits randi([0 1], 1, Nbits); % 随机数据比特 data 2*bits - 1; % BPSK符号映射0→-11→1 pn mseq([5 2], 5); % 周期31的m序列双极性 % 每个bit与整段PN码逐码片相乘完成扩频 tx_chips zeros(1, Nbits * sf); for k 1:Nbits tx_chips((k-1)*sf (1:sf)) data(k) * pn; end逻辑说明扩频就是把一个±1符号与一段PN码逐码片相乘得到长度为sf的扩频序列。这里data(k)*pn相当于把BPSK调制与扩频同时完成——在基带等效模型里不需要单独写调制器。参数说明sf取31而非任意数是为了让PN码周期与比特边界天然对齐如果sf取32需要把PN码截断或补齐功能上没问题但会破坏序列的周期完整性相关旁瓣特性会改变不建议初版仿真这么做。如果想确认发射频谱是否正确可以顺手画一个功率谱fs 31e3; % 采样率等于码片速率 figure; pwelch(tx_chips, 256, [], [], fs); % 观察主瓣与旁瓣 title(扩频后信号功率谱);这里的谱形状是sinx/x型的码片谱主瓣零点出现在码片速率整数倍处直观展示了扩频把信号带宽从数据速率展宽到码片速率的过程。3.2 信道干扰建模AWGN与窄带干扰注入信道模型这一段噪声功率的定义是重灾区。我用基带等效模型噪声直接加在码片级采样上方差按Eb/N0换算而不是猜SNR。噪声设置代码与参数含义一起说明EbN0_dB 6; % 目标Eb/N0单位dB EbN0 10^(EbN0_dB/10); % 线性值 Eb sum(tx_chips .^ 2) / Nbits; % 每bit能量基带下等于sf N0 Eb / EbN0; % 噪声功率谱密度 noise_var N0 / 2; % 实噪声方差 noise sqrt(noise_var) * randn(1, length(tx_chips)); rx tx_chips noise; % 接收序列逻辑说明Eb在这里恰好等于sf因为每个比特有sf个幅值为1的码片能量累加为sfN0由Eb/EbN0反推实高斯噪声的单边功率谱密度对应方差为N0/2。注意如果是复基带模型同相和正交支路各占N0/4噪声总方差才是N0/2复噪声与实噪声的系数不同这是仿真结果偏移最常见的原因之一。窄带干扰的注入方式是在码片序列上叠加一个正弦信号频率落在数据带宽内幅度相对信号幅度可调fi 0.1 * (1/sf); % 干扰频率设为数据带宽内的偏置 jam_amp 0.5; % 干信比相对信号幅度1 jamming jam_amp * sin(2*pi*fi*(1:length(rx))); rx_jam rx jamming;参数说明干扰频率fi设成0.1倍数据速率相当于一个落在主瓣内的单音干扰jam_amp为0.5时单音功率比信号低约6dB此时直扩系统解扩后仍能正常工作但误码率会有劣化。把jam_amp调大到1.5以上可以压出直扩系统的干扰容限边界观察误码率如何随干扰功率快速爬升。3.3 接收端解扩解调同步、解扩与误码率统计接收端在没有定时偏差时解扩就是一组相关累加rx_bits zeros(1, Nbits); for k 1:Nbits seg rx_jam((k-1)*sf (1:sf)); % 取第k个bit对应的sf个码片 y sum(seg .* pn); % 与本地PN码逐点相乘累加 rx_bits(k) y 0; % 判决正值判1负值判0 end err sum(rx_bits ~ bits); ber err / Nbits;逻辑说明解扩本质是相关运算。发送端每个bit与PN码相乘接收端再次与同样的PN码相乘再累加PN码的随机极性在累加中相互抵消只剩数据符号的贡献PN码未对齐时累加结果趋近于零判决输出以50%概率错误。参数说明里注意解扩窗口长度必须严格等于sf多取一个码片或少取一个码片都会破坏正交性。为扫出完整的误码率曲线把上面所有步骤包进一个Eb/N0循环EbN0_dB_list 0:2:12; ber_sim zeros(size(EbN0_dB_list)); for idx 1:length(EbN0_dB_list) EbN0_dB EbN0_dB_list(idx); % 按当前Eb/N0生成噪声并叠加解扩判决记录ber % 该点要循环多次取平均避免误码数太少导致抖动 ber_sim(idx) run_ber_point(tx_chips, pn, sf, EbN0_dB, 20); end我在仿真里习惯把“单点误码率计算”封装成函数主脚本只负责参数扫描和画图这样调试同步问题时不用反复跑整个循环。单点循环次数建议保证误码数大于100否则低误码率区的曲线抖动会让人误判性能。4. 直扩系统的同步与验证码相位、频偏和误码率曲线4.1 滑动相关法实现码同步直扩仿真最容易被忽略的模块是同步。发射端和接收端不可能天然共享同一段PN码相位接收端必须从采样序列里找到扩频码的起始位置。滑动相关法原理不复杂本地PN码与接收序列逐码片滑动计算每个候选相位下的相关值相关峰值所在相位就是码边界。delay_true 17; % 人为预设的码片时延 rx_delay [zeros(1, delay_true), rx]; % 模拟未知时延 corr_abs zeros(1, sf); for shift 0:sf-1 seg rx_delay(shift1 : shiftsf); corr_abs(shift1) abs(sum(seg .* pn)); end [~, best_shift] max(corr_abs); % 峰值相位逻辑说明每个候选相位对应一个相关幅度sf个相位扫描完后最大值位置就是码起始位置。这里与发射端使用同一段pn序列所以峰值会出现在delay_true对应的shift处。参数说明滑动窗口长度等于sf如果采样率是码片速率的整数倍这里的步进就不是1而是过采样倍数滑动窗口的搜索范围也要相应扩大。滑动相关只能粗同步到一个码片以内。后续如果要做精细跟踪常用超前-滞后环delay-locked loop维持码相位的对齐。仿真里一般把粗同步结果直接作为最佳相位进入解扩除非你要做动态多径信道否则不需要把跟踪环也写进去。4.2 载波频偏对解扩性能的影响及仿真观测如果仿真带载波频偏问题逃不掉。发射端上变频用cos(2πfct)接收端下变频也用cos(2πfct)但两边振荡器不可能完全一致残余频偏Δf会让解扩前的信号带一个旋转因子。基带等效模型里这等价于每个码片乘上exp(j·2π·Δf·t)。Δf与码片速率Rc的比值是决定影响程度的关键比值小于0.01时误码率损失很小比值达到0.1量级相关峰会因为相位旋转被抹平误码率直接到0.5。仿真里观测频偏影响的标准做法是固定Eb/N0扫描Δfdf_list [0 0.01 0.05 0.1 0.2] * (1/sf); % 相对数据速率的归一化频偏 for idx 1:length(df_list) df df_list(idx); t (0:length(tx_chips)-1) / (31e3); % 时间轴fs31kHz rot exp(1j*2*pi*df*t); % 频偏相位旋转 rx_df tx_chips .* rot noise; % 含频偏的接收序列 % 解扩、统计误码率记录到ber_df(idx) end参数说明df_list里的0.01、0.05是相对数据速率Rb的归一化频偏换算成实际频率要乘以Rb。比如Rb1kHz时0.05对应50Hz频偏这个量级的频偏在普通BPSK里影响不大但在直扩解扩里已经会看到明显的性能下降。扫参结果符合工程直觉直扩的扩频增益只对抗干扰有效对抗频偏的能力与未扩频BPSK相同并不存在额外优势。4.3 误码率曲线验证与理论值的对比方法仿真结果必须与理论曲线对比才能证明代码正确。直扩BPSK在AWGN信道下的理论误码率与未扩频BPSK完全一样Pb Q(sqrt(2·Eb/N0))。很多第一次做直扩仿真的人对这条结论很困惑——处理增益去哪了答案在于AWGN本身就是平坦白噪声扩频把信号带宽展宽、噪声功率同时增加比特级信噪比没有改善处理增益只有在对抗窄带干扰或窄带信道衰落时才体现出来。仿真时如果拿纯AWGN信道去验证“扩频增益”必然得不到想象中的性能提升这不是bug而是原理使然。画曲线时的注意事项项要求原因横轴Eb/N0单位dB与理论公式对齐不能用SNR纵轴log10坐标高信噪比区误码率跨度大对比点理论Q函数曲线偏差超过0.5dB要排查噪声换算最低误码数大于100否则统计抖动大低误码区失真验证时发现偏差先查噪声方差定义再查解扩窗口是否对齐。如果仿真点整体比理论好多半是噪声功率设小了如果整体差基本是扩频码相位偏了半个码片。5. 直扩仿真避坑指南5个最常见的翻车点和排查方法5.1 误码率高得离谱先检查扩频码是否对齐现象发射和接收代码看起来都没问题但误码率固定在0.5附近怎么调Eb/N0都没有反应。原因解扩用的PN码与接收序列中的PN码相位没对齐。直扩系统对码相位非常敏感哪怕只差一个码片相关累加结果也趋近于零判决输出基本是随机猜。解决先把发射端到接收端的码片延时人为设为0跑通后再加滑动相关同步模块。排查时打印解扩累加值的绝对值分布如果所有值都接近0基本可以断定是相位没对上。用一句话说直扩解扩不存在“半对齐”的模糊地带要不就全对要不就全错。5.2 采样率与码片速率不匹配频谱混叠导致结果失真现象误码率曲线在低Eb/N0处正常高Eb/N0处反而反弹出现类似地板的效应怎么加噪声都没用。原因采样率低于奈奎斯特要求或没有加成型滤波器码片波形在采样点之间混叠符号间干扰累积。初版用每码片一个采样没问题但只要引入滤波器或过采样采样率就必须是码片速率的整数倍。解决把采样率设置为码片速率的4倍并在收发两端加入根升余弦滤波器。滤波器滚降系数取0.22即可接收端滤波器要与发送端匹配否则会额外引入12dB损耗。仿真里常见做法是用rcosdesign生成滤波器系数发送端成型、接收端匹配滤波。5.3 处理增益计算口径不一致现象参考代码里写处理增益15dB自己仿真测出来只有3dB差别大得离谱。原因处理增益公式10lg(sf)中sf必须等于码片速率除以数据速率如果数据速率用的是信息比特速率而链路里还有纠错编码口径就变了需要区分信息速率与信道速率。解决仿真里明确标注Rb和Rc在代码注释里写明“处理增益基于未编码信息速率”。如果需要对比系统方案统一到同一个基准否则不同编码速率、不同扩频因子的方案之间没法横向比较。5.4 Eb/N0与SNR的换算雷区现象仿真曲线整体偏移2~3dB而且偏移量与扩频因子相关换成不同sf后偏移量跟着变。原因Eb/N0和SNR之间差一个带宽因子SNR Eb/N0 × (Rb/BW)。BW在过采样和滤波后会变化直接按SNR去设置噪声功率结果必然漂移而且漂移量随扩频因子改变。解决用Eb/N0作为噪声设置入口按3.2节的方差公式生成噪声不走SNR。如果一定要用SNR做接口务必把等效噪声带宽算准确并且在文档里记录换算过程。这个坑我踩过不止一次最后统一到Eb/N0口径后才消停。5.5 同步模块的“假锁定”现象现象滑动相关输出存在多个峰值取了最大值相位后误码率仍然很高数据恢复不出来。原因低信噪比或窄带干扰条件下相关输出存在多个局部极大值直接选最大值可能锁到旁瓣上尤其PN码周期较短如31时旁瓣水平相对主瓣不够低。解决把捕获门限设为理论主峰值的50%以上并要求峰值与次大值之比大于2才认定同步成功。仿真里可以统计同步成功概率作为第二个指标而不是只看误码率。多径信道下假锁定更容易出现这时候用Gold码这类互相关有界的序列配合多径搜索能有效降低假锁概率。6. 把仿真结果用起来一套验证技巧与日常检查习惯仿真跑通不等于结果可信我给自己定了一套“三重验证”流程每次直扩仿真都走一遍可以直接拿来当检查清单。第一重单点波形验证。固定一组参数把发射序列、含噪接收序列、解扩累加值三段画在同一张图上眼睛确认扩频结构存在、噪声幅度合理、解扩累加值明显分成正负两簇。这一步花两分钟能挡掉八成低级错误。第二重理论曲线对照。用AWGN信道跑出误码率曲线和理论BPSK曲线叠图。偏差超过0.5dB时先查噪声方差的定义再查解扩窗口长度。这里要强调一个原则直扩在纯AWGN下没有信噪比增益如果有人告诉你“扩频后误码率应该更好”直接用这条原理反驳别浪费时间调参数。第三重干扰容限实测。在信道里注入单音干扰逐渐加大干扰幅度画出误码率随干信比变化的曲线。正常情况下曲线会有一个明显的“膝盖”低于某个干信比时误码率保持在一个较低水平超过后快速恶化。这个膝盖位置与理论处理增益大约差3dB以内如果差得远说明同步或滤波环节还有问题。日常仿真我还要检查随机数种子和循环次数。误码率统计必须保证每个Eb/N0点的误码数超过100低误码率点需要更多比特随机数种子固定下来方便复现。做参数扫描时不要一次性跑所有点先用间隔4dB粗扫找到异常区间再局部细化能省下大量时间。这套流程坚持用下来直扩仿真的结果基本是可信的。我从第一次把PN码相位设错开始就养成了“先画波形、再对理论、最后抗干扰”的顺序希望你也能少走这些弯路。希望帮到你。本文还有配套的精品资源点击获取