MATLAB直接序列扩频仿真：PN码同步与误码率验证

简介面向通信工程、电子信息类专业学生及扩频通信入门者的MATLAB仿真设计文档围绕直接序列扩频DSSS系统的建模与性能评估展开可用于课程设计、实验报告撰写与原理验证。内容以BPSK调制为主线梳理直扩收发链路中信源、伪随机码生成、扩频调制、相关解扩与解调各环节并对m序列、GOLD码、WALSH码、OVSF码的自相关、互相关特性及在IS-95、TD-SCDMA等标准中的适用情况逐一说明同时给出信源速率2kbps、扩频因子64、高斯信道条件下的仿真参数设置与部分程序代码结合误码率随信噪比变化的曲线讨论处理增益以及抗同频干扰、抗衰落、抗多径的能力。资源包仅含1个doc文档约320KB篇幅紧凑公式、框图和代码并排呈现便于边读边在MATLAB中复现。已有143人学习下载适合希望从原理推导过渡到代码实现的读者作为参考底稿。1. 一条 31 位 PN 码把 -20 dB 的信噪比重新利用起来把接收机丢进 -20 dB 信噪比的信道里普通 BPSK 波形基本被噪声淹没判决输出接近随机换成本地一段 31 位 PN 码做相关相关输出端却重新出现清晰的符号峰。直接序列扩频做的就是这么一件事用码片率换处理增益把一条窄带信号摊到远大于信息带宽的频带上接收端靠同一段码把相位对齐的分量叠加回来。用 MATLAB 做这套系统的仿真设计重点不在于画几张好看的频谱图而在于能定量说清扩频增益从哪里来、同步差多少会崩、误码率曲线能不能和理论线对上。这篇内容面向正在做链路级仿真验证的工程师、要复现一条完整收发链路的学生也适合想把扩频模块放进自己仿真框架、却卡在同步和误码统计上的人。接下来按「算清楚 → 搭起来 → 同步 → 排错 → 验证」的顺序走一遍。2. 处理增益 15 dB 到底换来了什么链路预算与模块切分2.1 10log10(N) 换的是抗干扰容限不是抗白噪声扩频因子 N 31 时处理增益 Gp 10log10(31) ≈ 14.9 dB。很多人第一反应是「信噪比改善了 15 dB误码率曲线应该整体左移 15 dB」这是最常见的误解。在纯 AWGN 信道下扩频前后的误码率曲线在 Eb/N0 坐标上是重合的扩频把信号能量摊到更宽的频带上但每比特的总能量没变噪声单边功率谱密度 N0 也没变解扩时把 31 个码片相干累加信号幅度涨 31 倍、噪声功率涨 31 倍信噪比回到原来的水平。真正被改善的是窄带干扰容限和截获概率。干扰容限的常用估算是 Mj Gp − [S/N]o − Ls其中 [S/N]o 是解调所需的最小输出信噪比Ls 是实现损耗。取 Gp 14.9 dB、[S/N]o 10 dB、Ls 2 dB能扛住约 3 dB 的干扰余量。这个数看起来不大但意义在于一个占不到总带宽 1/30 的窄带干扰落进解扩器后只贡献 1/31 的相关输出而有用信号是全增益通过的。扩频因子 N处理增益 Gp占用带宽相对 Rb典型用途3114.9 dB31 倍教学链路、短帧控制信道6318.0 dB63 倍低速率遥测12721.0 dB127 倍抗窄带干扰场景25524.1 dB255 倍低信噪比捕获选 N 的时候不要只盯着 Gp。N 越大捕获搜索空间越大要搜 N 个码相位、同步建立越慢、仿真采样点越多。做链路验证时 N 31 已经足够把扩频、同步、解扩的全套流程跑通用 5 阶 m 序列一个周期刚好 31 个码片长度短、周期整、便于对齐是最省事的起点。等同步逻辑稳定了再往 127、255 扩改动只是寄存器阶数和本原多项式。2.2 五个模块的信号口径比特级、码片级、采样点级仿真跑不起来八成是信号长度口径混了。直接序列扩频链路上同时存在三种「速率」比特率 Rb、码片率 Rc Rb × N、采样率 fs Rc × sps。每个模块吃进去和吐出来的数组长度不一样reshape 一旦对不上就报错或者算出错的误码率。模块输入口径输出口径数组长度信源无比特级0/1 或 ±1Nsym扩频比特级码片级±1Nsym × N码片成型码片级采样点级浮点Nsym × N × sps载波调制采样点级采样点级带通Nsym × N × sps信道采样点级采样点级同上相干解调采样点级采样点级基带同上解扩与积分采样点级符号级浮点相关值Nsym判决与统计符号级比特级Nsym关键在于每符号采样点数 Ns fs / Rb N × sps 必须严格取整。取 Rb 50 bps、N 31、sps 8得到 Rc 1550 chip/s、fs 12400 Hz、Ns 248全是整数。如果随手把 fs 设成 12000 HzNs 240 但 N × sps 248两边对不上后面 reshape(z, Ns, Nsym) 会直接报维度错误或者更糟——不报错但每个符号窗跨到了下一个符号误码率曲线莫名其妙地永远 0.5。我的习惯是在脚本最前面把所有派生量一次算完并打一行断言Rb 50; N 31; sps 8; Rc Rb*N; fs Rc*sps; Ns fs/Rb; assert(abs(Ns - N*sps) 1e-9, 采样率与扩频因子不匹配); assert(abs(Ns - round(Ns)) 1e-9, 每符号采样点必须是整数);这两行断言挡掉的调试时间比后面所有排错加起来都值。2.3 载波频率与采样率的三个约束载波不是随便设的。第一条约束是奈奎斯特fc 必须明显小于 fs/2最好留出码片主瓣宽度的余量。码片率 1550 chip/s 的 sinc 主瓣零点在 1550 Hzfc 取 2 kHz 时上边带延伸到约 3.5 kHzfs/2 6.2 kHz余量够。第二条约束是一个符号周期内载波整数周期符号周期 20 msfc 2 kHz 正好 40 个周期积分清零时载波分量完全抵消不会在相关输出上留下残余纹波。第三条是码片与采样点的整数比sps 取 8 意味着每个码片 8 个样点码片边界严格落在样点格子上。sps 的取值需要权衡。取 1每码片一个样点跑得飞快但波形图看不出扩频形状滤波和延迟锁定的亚码片精度也无从体现。取 8 是常用的折中单个码片波形能看出矩形形状延迟锁定环的超前滞后间隔 0.5 码片等于 4 个样点够用仿真速度也不至于太慢。sps 提到 32 以上对误码率结果没有实质改善只是让频谱图更干净代价是仿真时间翻倍。还有一点常被忽略直接序列扩频信号是恒包络的 BPSK成型滤波如果用根升余弦接收端必须匹配滤波否则码片间干扰会吃掉 1~2 dB。教学链路里我一般先用矩形码片保持把同步和误码统计跑对再换成根升余弦滚降 0.35 观察性能损失这样出问题时能分清是同步的问题还是成型的问题。3. 从 m 序列到判决输出三段能在本地直接跑的 MATLAB 骨架3.1 m 序列生成与自相关自检扩频码选 5 阶 m 序列本原多项式 f(x) x^5 x^2 1周期 31。移位寄存器抽头在 5 和 2反馈为两者的异或初相只要非全零即可不同初相得到的是同一序列的循环移位。function pn mseq5(Nc) % 5 阶 m 序列f(x)x^5x^21周期 31 reg [1 0 1 0 1]; % 任一非全零初相 pn zeros(1, Nc); for k 1:Nc pn(k) reg(end); % 当前输出取最后一级 fb xor(reg(end), reg(end-3)); % 抽头 5 与 2 reg [fb, reg(1:end-1)]; % 右移并注入反馈 end pn 2*pn - 1; % 单极性转双极性0--1, 1-1 end生成完别急着往下走先做一次自相关自检。双极性 m 序列的周期自相关在零延迟处为 N其余延迟处为 −1这个「尖峰 平坦底部」的形状是后面捕获能工作的全部依据。用xcorr(pn,unbiased)画一下主峰 31、旁瓣约 −1如果旁瓣出现接近 1/3 主峰的值多半是抽头写错了或者寄存器方向反了此时生成的序列周期不是 31。再补一步互相关检查把circshift(pn, 1)和pn做相关结果应该是 −1 量级。这一步模拟的就是捕获时本地码相位差一个码片的情况能直观看到「错一个码片相关值就掉到噪声底」——这正是同步必须做码片级搜索的原因。3.2 扩频、载波调制与带宽观察扩频就是把每个符号乘以一整周期 PN用kron一步完成码片成型用重复采样实现注意矩阵展开方向。Rb 50; N 31; sps 8; fc 2e3; Rc Rb*N; fs Rc*sps; Ns fs/Rb; Nsym 2000; rng(7); % 固定随机种子结果可复现 bits randi([0 1], 1, Nsym); sym 2*bits - 1; pn mseq5(N); spread kron(sym, pn); % 比特级 - 码片级 wave reshape(repmat(spread, sps, 1), 1, []); % 码片级 - 采样点级 t (0:numel(wave)-1) / fs; tx wave .* cos(2*pi*fc*t); % BPSK 载波调制kron(sym, pn)展开后长度是 Nsym × N第 i 个符号占据第 (i−1)N1 到 iN 个码片。repmat(spread, sps, 1)生成 sps 行矩阵再reshape(..., 1, [])按列优先展开效果是每个码片值重复 sps 次长度变成 Nsym × Ns。带宽验证用pwelch比fft稳因为扩频信号功率谱密度低、动态范围大直接 fft 不加窗会泄漏得看不出形状。取 4096 点、汉宁窗、50% 重叠能看到主瓣宽度约在 ±Rc 附近、带内功率谱基本平坦的矩形轮廓。如果频谱看起来还是一根细线先查 sps 是不是取了 1码片波形根本没画出来再查 FFT 点数是不是小于一个符号周期的采样点数 248。3.3 相干解扩、积分清零与判决接收端做三件事乘本地载波恢复基带、乘本地 PN 完成解扩、按符号周期积分清零。三步可以合成一次逐点相乘再 reshape 求和写起来短也避开了滤波器引入的群延迟。EbN0_dB 6; Eb sum(tx.^2) / Nsym; % 每比特离散能量 N0 Eb / 10^(EbN0_dB/10); rx tx sqrt(N0*fs/2) * randn(size(tx)); % 实信号双边带噪声 ref reshape(repmat(kron(ones(1,Nsym), pn), sps, 1), 1, []); % 本地 PN 模板 z 2 * rx .* cos(2*pi*fc*t) .* ref; % 下变频 解扩 y sum(reshape(z, Ns, Nsym), 1); % 积分清零得到相关值 dec y 0; % 硬判决 ber sum(bits ~ dec) / Nsym;噪声方差取 N0·fs/2 是因为实信号双边功率谱密度是 N0/2带宽 fs 内的总噪声功率为 N0·fs/2。这个口径必须从头到尾一致否则误码率曲线会整体平移几个 dB。2*rx.*cos(...)里的 2 是为了把下变频后的幅度恢复到 1cos² 展开有 1/2 的直流分量乘ref后相关峰的理论值是 Ns这里不归一化也没关系判决只看符号。reshape(z, Ns, Nsym)是按列切分第 j 列正好是第 j 个符号的 Ns 个采样点sum(...,1)沿列求和就是积分清零。这一步同时完成了低通滤波的功能——两倍频分量在整符号周期积分下严格为零前提是 2.3 节说的整数周期约束。如果环境里没有 Communications Toolboxbiterr用不了就用上面这行sum(bits ~ dec)/Nsym结果一样。跑一个点确认 ber 在 1e-2 量级链路就通了接下来才是有意义的扫点。4. 码相位差半个码片就掉 3 dB捕获与跟踪的仿真做法4.1 滑动相关捕获搜索窗、门限与虚警捕获要解决的问题是接收信号里的 PN 相位未知得在 N 个可能的码片偏移里找出正确的那一个。最直接的做法是滑动相关——本地 PN 依次循环移位 0 到 N−1 个码片每个相位与该符号窗内的基带样点做相关取最大者。rxc 2 * rx(1:Ns) .* cos(2*pi*fc*t(1:Ns)); % 取一个符号周期剥载波 base rxc(sps/2 : sps : end); % 每码片取中点一个样点 corr zeros(1, N); for d 0:N-1 corr(d1) abs(sum(base .* circshift(pn, d))); end [peak, dhat] max(corr); % dhat-1 即估计的码相位取码片中点的原因很实际矩形码片在边界处有跳变样点落在边界上会把相邻码片的部分能量算进来相关值抖动大。中点采样是最省事的近似sps 8 时中点样点离两侧边界各 4 个样点足够干净。门限设置是捕获里最容易翻车的一环。理论上无噪声时正确相位相关值为 31、错误相位为 −1门限放 15 就能分得很开。但 Eb/N0 0 dB 时相关值本身就是一个随机变量正确相位的相关输出均值 31、标准差随噪声增大错误相位的输出均值 −1、标准差相近。门限卡太紧会漏检正确相位没超门限重新搜索卡太松会虚警错误相位超门限锁到错的位置上后面误码率永远是 0.5。工程上的做法是用虚警概率反推门限错误相位的相关输出近似高斯均值 −1、方差 σ²门限设为 −1 4σ 时单次虚警概率约 3e-5N 31 次搜索的总体虚警概率仍在千分之一量级。仿真里可以直接把噪声方差算出来代进去比手动试门限靠谱得多。4.2 延迟锁定环鉴相曲线与环路带宽捕获只保证对齐到一个码片以内剩下的小数倍码片偏差要靠跟踪。延迟锁定环用两路相关器本地码分别超前和滞后半个码片两者相关值之差就是鉴相误差符号指示调整方向。tau 0.5; % 超前滞后间隔单位码片 early circshift(pn, 1); % 近似超前 1 码片 late circshift(pn, -1); % 近似滞后 1 码片 E sum(base .* early); L sum(base .* late); err E - L; % 鉴相输出过零即对齐circshift只能做到整数码片移位真实的半码片超前滞后在采样点级实现sps 8 时半码片是 4 个样点直接对采样点索引做 ±4 的偏移即可精度比整数码片移位高一个量级。上面的写法是为了让鉴相曲线能在码片级快速验证——把本地码相位人为从 −1 码片扫到 1 码片画 err 随相位差的曲线应该看到一条过零的 S 曲线零点在相位差为 0 处斜率与相关峰斜率成正比。环路带宽决定跟踪速度和抖动。带宽取大能跟上快速的多普勒变化但稳态抖动大解扩相关值会有几个 dB 的随机衰减带宽取小抖动小但捕获后收敛慢仿真里要跑更长的前导序列。链路级仿真常用归一化带宽 BnT 在 0.005 到 0.02 之间T 是码片周期。这个参数没有普适最优值取决于你的信道是静态还是有时变先把静态信道跑对再引入时变。4.3 载波不同步与码不同步现象不一样两种不同步在误码率曲线上表现不同分清了能省一半调试时间。载波相位有残余偏差时相关输出按 cos(Δθ) 衰减Δθ 30° 时损失约 0.3 dBΔθ 60° 时损失约 1.2 dB误码率曲线整体右移但形状不变且在低信噪比下几乎看不出异常。码相位偏差时相关峰按三角波衰减差 0.5 码片损失约 3 dB差 1 个码片直接掉到旁瓣、相关值从 31 掉到 −1误码率瞬间冲到 0.5。判断方法很直接把发射端和接收端的本地码相位强制设成已知只加噪声跑一遍曲线正常再放开码同步如果曲线崩到 0.5问题在码相位而不是载波。反过来如果曲线只是整体平移几个 dB形状没变优先怀疑载波残余相位和频偏。5. 避坑直接序列扩频仿真里最容易出错的 5 类场景5.1 误码率稳定在 0.5 附近不动现象扫 Eb/N0 从 0 dB 到 10 dB误码率一直在 0.45~0.55 之间曲线横着走看起来像信道里根本没有有用信号。原因本地 PN 和发送端 PN 不是同一段或者相位差了一个码片以上。常见来源有三个接收端重新调用了一次mseq5但寄存器初相和发送端不同得到的是循环移位后的序列circshift的方向写反kron展开后 reshape 的维度顺序把码片顺序打乱了。这三种情况下相关输出都是噪声底判决完全随机。解决先在无噪声条件下跑一遍把y的前 5 个相关值打印出来看是不是 Ns 量级。如果是几百的小数而不是 248 附近问题确定在码不对齐。然后把接收端 PN 直接赋成发送端的pn变量不重新生成曲线立刻恢复正常就说明是初相不一致。养成习惯发送端和接收端共用同一个pn数组只允许在信道上引入延迟不要在接收端重新生成序列。5.2 低信噪比区间的曲线比理论还漂亮现象Eb/N0 0 dB 时仿真误码率 0.06而理论 BPSK 是 0.078看起来「仿真比理论还好」随着信噪比升高两条线又慢慢贴合。原因噪声功率归一化口径错了。Eb sum(tx.^2)/Nsym算出来的是离散采样后的每比特能量如果后面用awgn(tx, EbN0_dB, measured)又让它自己测一遍功率两个口径打架或者把每比特能量和每符号能量搞混扩频后符号能量是比特能量的 N 倍按符号能量算 N0 会让噪声偏小曲线整体左移约 15 dB。解决全程手动加噪不要混用awgn。固定用N0 Eb / 10^(EbN0_dB/10)、sigma sqrt(N0*fs/2)这一套Eb 用调制后信号的离散能量除以比特数。验证方法是把扩频因子 N 从 31 改成 1不扩频跑出来的曲线应该和 BPSK 理论线0.5*erfc(sqrt(10^(EbN0/10)))重合重合了说明能量口径对了再把 N 改回 31。5.3 频谱图上看不出扩频还是一根细线现象pwelch画出来的功率谱是窄窄一根主瓣宽度远小于 Rc看不出被摊开的形状。原因最常见是画频谱时用的信号是比特级或码片级数据没有经过重复采样等于把一个码片当一个样点采样率变成 Rc频谱自然只有一条。其次是 FFT 点数取得太少小于一个符号周期的 248 个样点频率分辨率不够。第三种是没加窗扩频信号本身是宽带的不加窗泄漏会把旁瓣抬得很高掩盖主瓣结构。解决画频谱前先assert(numel(x) Nsym*Ns)确认输入是采样点级信号。pwelch用 4096 点、汉宁窗、50% 重叠fs参数显式传进去。如果主瓣宽度看起来还是不对把 sps 临时提到 16 再看主瓣轮廓会更清楚——这一步只是诊断用不影响误码率结论。5.4 同步环锁定后相关值抖得厉害现象捕获和跟踪都跑通了但解扩后的相关峰值在理论值附近上下跳 30% 甚至更多误码率比预期高 2~3 dB。原因环路带宽取太大或者鉴相增益没有归一化。延迟锁定环的鉴相输出幅度与相关峰斜率成正比相关峰幅度又随信噪比变化如果不除以相关峰幅度做归一化低信噪比时环路增益自动变小、高信噪比时自动变大等效带宽随信噪比漂移。另一个来源是超前滞后间隔取得太宽比如 1 个码片鉴相曲线在零点附近斜率小、线性区窄。解决把鉴相输出除以abs(peak)做幅度归一化超前滞后间隔固定在 0.5 码片。仿真初始化时给环路 200 个符号的收敛期统计误码时从第 201 个符号开始避免把收敛过程中的错误算进去。这个「预热期」的坑很隐蔽不扣掉的话低信噪比下的误码率会偏高一个数量级。5.5 高信噪比端曲线出现平台现象Eb/N0 超过 8 dB 后仿真曲线不再下降停在 1e-4 附近打平而理论线继续往下走。原因蒙特卡洛样本数不够。2000 个比特在 8 dB 时理论误码率约 2e-4期望错误比特数只有 0.4 个统计上完全不可信。另一个来源是判决门限固定为 0而实际相关输出有直流偏置积分清零不彻底、或者载波周期不是整数残留了两倍频分量。解决按「至少 100 个错误比特」定样本数Eb/N0 8 dB 时需要约 50 万比特。分段仿真、每段独立加噪比一次性生成大数组省内存。门限偏置可以在无噪声条件下测一次相关输出的均值判决时减去这个均值再和 0 比较代码只多一行。6. 把仿真曲线压到 erfc 理论线上可复现的验证流程6.1 扫点流程与样本数下限验证的核心是扫一组 Eb/N0每个点独立生成噪声、独立统计误码和0.5*erfc(sqrt(10^(EbN0/10)))对照。注意这条理论线是 BPSK 在 AWGN 下的结果扩频不改变它——2.1 节说过处理增益换的不是抗白噪声能力。如果你的扩频曲线明显偏离这条线先查能量口径别急着怀疑扩频逻辑。Eb/N0 (dB)理论 BER建议比特数期望错误数07.9e-25,00039523.8e-210,00038041.3e-230,00039062.4e-3100,00024081.9e-4600,000114扫点循环写成函数输入 EbN0 和比特数输出误码率外层用arrayfun或普通 for 循环都行。关键是把随机种子显式固定同一个点跑两次结果必须完全一致否则曲线抖动你分不清是统计噪声还是逻辑错误。6.2 三条排查路径曲线整体左移 15 dB 左右基本可以确定是把扩频增益重复算进了噪声——要么噪声按符号能量算的要么 Eb 用了扩频前的比特能量。曲线整体右移 3 dB 以内优先怀疑载波残余相位或者成型滤波没做匹配滤波。曲线形状对但低信噪比端翘得比理论高检查同步预热期有没有扣掉以及捕获阶段有没有漏检导致整帧丢弃。我自己的习惯是每改一次参数先跑 N 1 的退化版本和 BPSK 理论线对一遍对上了再打开扩频。这个「退化对照」看起来多花十分钟实际省下的时间是它的好几倍——扩频链路里同时有码相位、载波相位、能量口径、采样率四个变量在动不把变量一个个锁死血泪经验就是只能靠猜。最后一个建议把参数区集中放在脚本最前面包括 Rb、N、sps、fc、EbN0、随机种子全部写死不要分散在函数里。仿真跑不顺的时候能一眼看完全部输入的人永远比翻十个文件找默认值的人先找到问题。希望帮到你。本文还有配套的精品资源点击获取