简介本资源是一套基于MATLAB实现的IEEE 802.11b无线网络仿真系统面向通信工程、无线网络课程设计及科研入门的学习者聚焦Wi-Fi物理层与MAC层建模核心问题帮助理解DSSS扩频、CCK/BPSK/QPSK调制、瑞利衰落信道建模及BER性能评估等关键技术。压缩包共12个文件75KB含6个MATLAB脚本如test_level_1.m、ber_test.m、WiFi_init.m负责信号生成、误码率测试与系统初始化、4个Simulink模型文件WiFi.mdl与WiFi_lib.mdl等构建端到端链路仿真框架以及2个数据文件cck_codes.mat等预存标准CCK编码表。已有362人学习下载资源结构清晰、模块解耦合理提供可直接运行的完整仿真流程——从基带信号生成、信道损伤注入到接收端解扩解调与性能可视化配套代码注释充分便于教学演示、算法验证与二次开发。1. 我为什么又翻出802.11b做MATLAB仿真一个“老协议”的现代价值最近在整理手头的仿真代码时我又把802.11b的MATLAB仿真翻了出来——这套用DBPSK加Barker码实现的WiFi物理层链路虽然速率只有1Mbps却是我觉得最适合用来理解无线通信仿真的入门案例。很多人一上来就盯着WiFi 6/7的MIMO-OFDMA反而把最基本的扩频、调制和判决逻辑丢了。这篇文章就跟大家聊聊我如何用MATLAB把802.11b从发射到接收完整仿真出来以及在这个过程里踩过的坑。适合那些想入门无线网络仿真、又不想被复杂标准淹没的读者也适合需要快速验证2.4GHz链路性能的工程师。1.1 从WiFi物理层演进看802.11b的位置802.11b是1999年发布的WiFi标准工作在2.4GHz频段最高物理层速率只有11Mbps。拿今天的眼光看这个速率确实不够看但它有一个非常特殊的位置它是把一个无线比特从“如何调制”到“如何扩频”讲得最干净的协议。后边的802.11a/g开始用OFDM要处理子载波、导频、循环前缀、交织、信道估计复杂度一下子高了一个量级再往后WiFi 4/5/6的MIMO、波束成形、OFDMA更是刚入门就被术语淹没。802.11b不一样。它用的是直序扩频DSSS每个数据比特被一个11码片的Barker码扩频再用DBPSK或DQPSK调制。整个物理层链路可以精简成“比特 - 扩频 - 调制 - 过采样 - 加噪 - 匹配滤波 - 解扩 - 判决”这条流水线每一环都能在MATLAB里用几十行代码写清楚。对我个人来说做802.11b仿真不是为了追新而是为了把通信链路里最核心的概念——带宽、信噪比、扩频增益、相关接收——彻底搞明白。1.2 仿真能解决哪些实际工程问题很多人问我现在做WiFi仿真是不是都应该直接用网络仿真器跑NS-3我的回答是看你要回答什么问题。如果你关注的是TCP吞吐量、路由协议、节点调度那应该去用系统级仿真工具但如果你关注的是“某个距离下能不能听到信号”“2.4GHz多径环境下1Mbps和11Mbps速率谁更稳”那就绕不开物理层的调制解调和误码率这部分恰恰是MATLAB最擅长的。MATLAB仿真802.11b能解决几类实际问题链路预算验证给定发射功率、天线增益、路径损耗通过仿真得到接收端的误码率或丢包率比只算一条静态链路预算更接近真实。扩频增益的直观验证在AWGN信道下扩频前后的误码率理论值不应该改变但抗窄带干扰的能力会显著提升这个结论用仿真图展示出来比看公式有说服力得多。多径灵敏度的扫描2.4GHz室内环境时延扩展通常在50~100ns左右接近802.11b的码片周期约91ns所以多径是实际部署中躲不开的问题。用简单的两径信道模型就能看到误码率曲线出现地板效应这对接入点布放有直接参考价值。为协议级仿真做准备物理层的BER/PER曲线是MAC层CSMA/CA仿真的输入之一。把物理层跑稳了后面加退避算法、RTS/CTS才算有扎实的下层支撑。2. 动手前必须想清楚的802.11b物理层细节2.1 四种速率与调制/扩频方案的对应关系802.11b支持四种物理层速率很多教程只关注速率值没有真正说清速率为什么是1、2、5.5、11这四档。其实这四档对应的是完全不同的调制和编码组合我整理成了一张表速率调制方式扩频方案码片速率数据码元速率备注1 MbpsDBPSKDSSS 11位Barker码11 Mchip/s1 Msymbol/s每个bit直接扩频2 MbpsDQPSKDSSS 11位Barker码11 Mchip/s1 Msymbol/s每个符号携带2 bit5.5 MbpsDQPSKCCKCCK 8码片11 Mchip/s1.375 Msymbol/s每个符号4 bit11 MbpsDQPSKCCKCCK 8码片11 Mchip/s1.375 Msymbol/s每个符号8 bit这里最容易让新手犯迷糊的是“码片速率固定但数据速率不同”这件事。以1Mbps和2Mbps为例两者的射频码片速率都是11Mchip/s带宽都是约22MHz。区别在于1Mbps模式里1个数据bit先被扩展成11个Barker码片然后每个码片做DBPSK调制所以每秒传1M个bit、对应11M个码片2Mbps模式里2个bit映射成1个DQPSK符号再扩展成11个码片所以每秒传1M个符号、对应2Mbit数据、同样是11M个码片。换句话说扩频后的射频符号速率始终是1Msymbol/s数据速率靠调制阶数提高而不是靠提高带宽。2.2 仿真链路的基础参数做MATLAB仿真前我习惯先把参数表写在注释里避免到后面算信噪比时全部乱掉。以下是这次仿真的基础参数参数数值说明载波频率2.4 GHz802.11b工作频段码片速率11 Mchip/s扩频后的码片速率数据调制DBPSK / DQPSK本次重点做DBPSK扩频码Barker码长度11[1 1 1 -1 1 -1 -1 1 -1 -1 1]过采样率16每个码片用16个采样点表示仿真比特数1e5 ~ 1e6根据噪声环境动态调整信道模型AWGN / 两径衰落两径时延可配置为什么把过采样率设为16这是一个工程权衡。过采样率太低比如等于1仿真正弦波和多径延迟时根本表示不了时间上的精细移动太高比如等于64仿真时间一下子拉长噪声计算也容易出错。16是一个够用又不慢的默认值特别是加入多径信道时0.5个码片时延可以用8个采样点表达精度基本够了。2.3 Barker码自相关特性的价值很多人会忽略Barker码在接收端的作用。它其实是一个“自带同步辅助工具”的序列长度为11的Barker码具有一个非常好的性质当它和自己对齐时相关值最大达到11当它偏移一个或多个码片时相关值的绝对值不超过1。这意味着接收机只要做滑动相关就能找到扩频码的起点这就是802.11b能快速完成码片同步的底层原因。在MATLAB仿真里我们通常假设同步已经完成但如果你要写一个完整的接收机Barker码的自相关峰就是你的“定时度量”。这个细节在仿真里很容易被忽略真正做硬件或软件接收机时却至关重要。我建议大家仿真跑通后自己画一下接收信号和本地Barker码在不同偏移下的相关输出你会对“扩频系统抗干扰能力强”这句话有更直接的理解。3. MATLAB仿真链路搭建从比特流到基带波形3.1 发射端调制、扩频和过采样我这次没有用MATLAB Communications Toolbox里的现成调制器而是自己手写发射端。原因是802.11b的扩频细节很固定自己写反而更容易控制每一个中间信号方便中途画波形图检查。下面是核心代码% 发射端参数 numBits 1e5; % 统计比特数 sps 16; % 每码片采样数 barker [1 1 1 -1 1 -1 -1 1 -1 -1 1]; % 11位Barker码 % 随机比特 dataTx randi([0 1], numBits, 1); % DBPSK调制0 - 1, 1 - -1 modSymbols 1 - 2 * dataTx; % 扩频每个符号乘整段Barker码得到 numBits*11 的码片序列 chipSeq reshape(modSymbols * barker, [], 1); % 过采样每个码片重复 sps 次 waveform repelem(chipSeq, sps); % 归一化让每个数据比特的能量等于1 % 每个bit包含11个码片每个码片采样sps个点所以总共11*sps个采样点 waveform waveform / sqrt(11 * sps);最后一行归一化很关键。这里我把每个数据比特对应的采样点数固定为11 * sps把这些采样点上的总能量归一化为1。这样后面加噪声时噪声功率只和Eb/N0相关不会因为过采样率变化导致信噪比算错。这也是我在第5章要重点讲的坑。如果你追求更真实的谱形状可以在这里加一个升余弦脉冲成形滤波器替代简单的repelem。repelem生成的是矩形脉冲频谱旁瓣很高不符合802.11b的频谱模板。不过做误码率仿真时矩形脉冲对AWGN下的结果没有本质影响所以我一般先用它跑通链路再改成rcosdesign观察波形和频谱。3.2 AWGN信道与两径衰落信道信道部分我做了两层第一层是AWGN第二层是两径衰落。AWGN的噪声方差设置为EbNoVec 0:2:12; % Eb/N0单位dB for idx 1:length(EbNoVec) EbNoLin 10^(EbNoVec(idx)/10); noiseStd sqrt(1 / (2 * EbNoLin)); noise noiseStd * randn(size(waveform)); rxWaveform waveform noise; % ... 接收与统计BER ... end这里噪声方差用了1/(2*EbNoLin)原因很简单发射波形已经按每比特能量归一化为1所以Eb1那么Eb/N0 1/N0 EbNoLin得到N0 1/EbNoLin复基带噪声的实部和虚部各占N0/2因此单路噪声标准差是sqrt(N0/2)。如果只做实信号仿真写sqrt(1/(2*EbNoLin))就是对的。两径信道我一般这样模型化% 两径信道主径 延迟径 delayInSamples sps; % 延迟1个码片 attenuation 0.5; % 延迟径幅度 rxWaveform waveform attenuation * [zeros(delayInSamples,1); waveform(1:end-delayInSamples)]; rxWaveform rxWaveform noise;这个模型虽然简单但只要把延迟时间设置成接近室内环境的时延扩展就能看到频率选择性衰落对802.11b的影响。这里有一个容易犯的错误如果用zeros直接补头但接收端没有把延迟径造成的干扰纳入考虑BER会表现得很不合理。实际做的时候我会把主径信号和延迟径信号分别画出来确认相对时延和幅度都符合预期再继续。3.3 接收端相关解扩和判决接收端可以按照“匹配滤波 - 解扩 - 判决”的顺序做。由于发射端用的是矩形脉冲匹配滤波可以简化成对一个码片内的sp个采样点求平均。之后把11个码片和本地Barker码做内积得到的软值符号就是判决结果% 接收端处理 rxChipMatrix reshape(rxWaveform, sps, []); % 每一列是一个码片 rxChipAvg mean(rxChipMatrix, 1); % 码片级匹配滤波 % 多少个完整bit numBitsRec floor(length(rxChipAvg) / 11); rxSoft zeros(numBitsRec, 1); for k 1:numBitsRec segment rxChipAvg((k-1)*111 : k*11); rxSoft(k) sum(segment .* barker); % 相关解扩 end % 判决软值是负号 - 原始比特为1正号 - 原始比特为0 dataRx double(rxSoft 0); dataRx dataRx(1:numBits); % 截断多余部分 ber sum(dataTx ~ dataRx) / numBits;这段代码没有包含码片同步和载波同步原因是我们在理想AWGN仿真里默认已知起始位置。如果你要仿真真实突发传输还需要把前导码和SFD加进去用前导码做同步估计。我觉得初学阶段可以先跳过同步集中精力把能量和BER算对等把“信号处理链路”跑熟后再补同步环。4. 仿真结果分析DBPSK与DQPSK的误码率对比4.1 完整的仿真循环结构把上一章的代码拼起来套一个遍历Eb/N0的循环就是一个可用的误码率仿真脚本。为了节省时间我会在低Eb/N0时用更少的比特高Eb/N0时增加比特数否则在高信噪比下统计不到足够多的错误比特BER曲线会出现毛刺numBits 5e4; % 低SNR先用5e4 for idx 1:length(EbNoVec) if EbNoVec(idx) 8 numBits 5e5; % 高SNR增加比特数 end % 发射、加噪声、接收、统计BER end这样做的好处是在高Eb/N0区间仍然能统计到几十到几百个错误比特BER曲线更平滑。如果你只是画一个趋势图可以不用这么严谨但要拿曲线做工程参考这个细节就很重要了。4.2 结果曲线、理论值和我的解读在AWGN信道下DBPSK相干解调的理论误码率是BER 1/2 * erfc(sqrt(Eb/N0))这个公式可以放到MATLAB里直接画理论曲线再和仿真点比较。实测下来只要信噪比计算正确仿真点和理论曲线几乎重合。误差通常出现在Eb/N0高的时候原因是错误比特太少统计波动变大。我再用DQPSK做对比时会特意把差分编码的损耗也体现出来。802.11b里的2Mbps模式使用DQPSK接收端一般做差分检测损失比相干检测大约2~3dB。从这个角度看802.11b选择DBPSK/DQPSK不是随便定的而是为了在2.4GHz这个“噪声和干扰都不太友好”的频段里用简单的非相干或差分接收机稳定可靠地传信号。还有一个经常被初学者忽略的现象在AWGN信道下扩频并不会改变DBPSK的误码率曲线。很多学员做完DSSS仿真后发现BER和没扩频时一样就会怀疑自己代码是不是错了。其实不是扩频的真正收益在抗窄带干扰和多径分辨不在AWGN下的能量效率。你可以在接收端叠加一个单音干扰再对比有无扩频的BER就会看到显著差异。5. 从“能跑通”到“跑得对”仿真链路调试实录5.1 信噪比计算错误是最隐蔽的坑我在第一次写802.11b仿真时BER曲线比理论值差了很多花了两天时间排查最后发现是噪声方差算错了。很多教程在加AWGN时只考虑信号的平均功率为1然后写noiseStd sqrt(1/(2*10^(EbNo/10)))这在我前面第3章的归一化前提下是对的但前提是信号平均功率必须等于每比特能量。问题是如果做过采样后没有把所有采样点归一化到“每比特能量为1”那么实际信号采样幅度就不是1。比如我最初用repelem扩频后直接让Barker码片幅度保持1那每个比特包含11个码片、每码片16个采样每个采样幅度是1总能量就是11*16176而不是1。此时再用噪声方差1/(2*EbNoLin)去加噪声实际Eb/N0就被拉高了BER会显得比理论好好到不真实。解决方式就是我在发射端做的归一化waveform waveform / sqrt(11 * sps);这样每比特总能量严格为1之后再用统一公式加噪声就不会错。这个坑几乎每个人都会踩一次所以我不建议直接复制网上的信噪比公式强烈建议先用sum(waveform.^2)/numBits打印一下每比特能量确认它是1再跑BER。5.2 解扩时的码片极性问题802.11b的Barker码有两种极性定义有的资料说数据0发送Barker码本身数据1发送Barker码的反码也有的资料定义相反。只要发射端和接收端保持一致仿真结果是完全等价的。真正容易出错的是你在解扩时如果直接用原始Barker码做相关但发射端用了反码映射那么所有软值符号都会翻转误码率会接近100%。排查这个问题有个小技巧不要一上来就统计BER先取10个比特把发射端的调制符号、扩频码片和接收端解扩软值逐段打印出来。确认软值正负和原始比特的映射关系一致后再跑大规模仿真。这个“先手算10个bit”的习惯能帮你省下大量调代码的时间。5.3 多径衰落下的“伪BER”问题在第3.2节的两径信道中如果时延不是码片周期的整数倍MATLAB里要精确模拟非整数采样点延迟需要用到插值或频域移位。很多初学者图省事直接把延迟设为比如1.3个采样点然后用round取整。这样表面上代码没有报错实际上造成了一个“伪信道”延迟径和主径的相对相位被随机化了每次跑的BER差别很大甚至出现高SNR下BER不下降的假象。我的做法是要么把时延严格设为整数倍采样周期做一个“准”多径仿真要么使用interp1插值把信号移动到亚采样精度再用频域插入的方法实现真正的分数延迟。初期仿真建议先用整数倍时延至少保证结果可复现。等确认多径机理清楚了再引入更真实的信道抽头系数。5.4 仿真速度与内存的平衡如果numBits1e6、sps16扩频后波形长度是1e6*11*161.76e8个采样点一个double数组占1.4GB内存MATLAB很可能直接卡死。哪怕内存够每遍历一个Eb/N0点都跑一遍这么长序列脚本总时间也会让人崩溃。我常用三个策略分段处理把发射序列切成若干段每段numBitsSeg1e4逐段加噪声并解调最后把所有错误比特累加。降低过采样率如果只做AWGN仿真sps4甚至sps2都够只有需要展示脉冲波形或做多径时才用到16。用矩阵化替代循环接收端解扩那层循环其实可以改成矩阵乘但前提是11个码片在内存里连续排列当序列很长时直接矩阵乘可能比循环更慢所以要测试后决定。6. 往系统级走802.11b仿真还能玩出什么花6.1 把PER和MAC层CSMA/CA结合起来物理层BER仿真的下一步是把误码率转换成误包率PER。802.11b的MAC帧长度通常几百到一千多字节在MATLAB里可以用一个简单公式近似PER 1 - (1 - BER)^(frameBits)。有了不同信噪比下的PER值就可以搭建一个轻量级的系统级仿真来模拟CSMA/CA退避过程。我不是说MATLAB一定比NS-3适合系统级仿真而是如果你已经用MATLAB把物理层跑得很熟再往上叠加MAC层时数据接口是最顺的。你可以用一个时间推进循环让多个节点在不同时刻产生数据包用物理层查表方式决定这个包是否被正确接收。这样虽然不如NS-3严谨但做课程设计或技术方案预研足够了。6.2 多节点干扰与隐藏终端场景2.4GHz频段里真正影响WiFi性能的往往不是热噪声而是其他节点造成的同频干扰。你可以在接收端叠加一个干扰信号比如另一个调制后的802.11b信号观察BER随信干比的变化。这个实验能直观解释为什么“一个屋子里终端太多WiFi就会变卡”。隐藏终端场景还可以把CSMA/CA的RTS/CTS机制加进来两个节点同时向一个接入点发包物理层模型给碰撞一个丢包概率MAC层统计吞吐量。这个方向对做无线网络课程项目的人来说非常实用而且不需要改动物理层核心代码只要在上层加时间调度逻辑。6.3 工具箱与手写代码怎么选MATLAB的Communications Toolbox提供了很多现成函数比如dpskmod、dpskdemod、rcosdesign可以节省码代码时间。WLAN Toolbox则主要面向基于OFDM的802.11a/n/ac/ax对802.11b的DSSS/CCK模式支持相对有限。所以我的建议是802.11b物理层仿真完全值得手写一遍因为代码量不大却能让你彻底吃透扩频和调制的细节等后面做OFDM链路时再用WLAN Toolbox效率会高很多。如果你决定使用工具箱也别忘了检查它的参考波形是否带前导码和Signal字段。工具箱生成的波形通常是为了测试接收机算法不一定完全符合802.11b的时域帧结构。在这个问题上手写代码反而更容易控制。关于802.11b仿真我最想分享的经验是先把AWGN单用户链路跑通再把误差率公式和仿真结果对齐最后才叠加多径和干扰。很多项目做不好不是因为MATLAB不熟而是把信道、同步、信道估计这些因素一股脑堆进去最后分不清是协议问题还是代码问题排查成本非常高。这套“先简单后复杂、先单点后系统”的方法同样适用你后续做802.11a的OFDM仿真甚至延伸到5G NR链路。本文还有配套的精品资源点击获取
基于MATLAB的802.11b物理层仿真：从扩频调制到误码率分析
简介本资源是一套基于MATLAB实现的IEEE 802.11b无线网络仿真系统面向通信工程、无线网络课程设计及科研入门的学习者聚焦Wi-Fi物理层与MAC层建模核心问题帮助理解DSSS扩频、CCK/BPSK/QPSK调制、瑞利衰落信道建模及BER性能评估等关键技术。压缩包共12个文件75KB含6个MATLAB脚本如test_level_1.m、ber_test.m、WiFi_init.m负责信号生成、误码率测试与系统初始化、4个Simulink模型文件WiFi.mdl与WiFi_lib.mdl等构建端到端链路仿真框架以及2个数据文件cck_codes.mat等预存标准CCK编码表。已有362人学习下载资源结构清晰、模块解耦合理提供可直接运行的完整仿真流程——从基带信号生成、信道损伤注入到接收端解扩解调与性能可视化配套代码注释充分便于教学演示、算法验证与二次开发。1. 我为什么又翻出802.11b做MATLAB仿真一个“老协议”的现代价值最近在整理手头的仿真代码时我又把802.11b的MATLAB仿真翻了出来——这套用DBPSK加Barker码实现的WiFi物理层链路虽然速率只有1Mbps却是我觉得最适合用来理解无线通信仿真的入门案例。很多人一上来就盯着WiFi 6/7的MIMO-OFDMA反而把最基本的扩频、调制和判决逻辑丢了。这篇文章就跟大家聊聊我如何用MATLAB把802.11b从发射到接收完整仿真出来以及在这个过程里踩过的坑。适合那些想入门无线网络仿真、又不想被复杂标准淹没的读者也适合需要快速验证2.4GHz链路性能的工程师。1.1 从WiFi物理层演进看802.11b的位置802.11b是1999年发布的WiFi标准工作在2.4GHz频段最高物理层速率只有11Mbps。拿今天的眼光看这个速率确实不够看但它有一个非常特殊的位置它是把一个无线比特从“如何调制”到“如何扩频”讲得最干净的协议。后边的802.11a/g开始用OFDM要处理子载波、导频、循环前缀、交织、信道估计复杂度一下子高了一个量级再往后WiFi 4/5/6的MIMO、波束成形、OFDMA更是刚入门就被术语淹没。802.11b不一样。它用的是直序扩频DSSS每个数据比特被一个11码片的Barker码扩频再用DBPSK或DQPSK调制。整个物理层链路可以精简成“比特 - 扩频 - 调制 - 过采样 - 加噪 - 匹配滤波 - 解扩 - 判决”这条流水线每一环都能在MATLAB里用几十行代码写清楚。对我个人来说做802.11b仿真不是为了追新而是为了把通信链路里最核心的概念——带宽、信噪比、扩频增益、相关接收——彻底搞明白。1.2 仿真能解决哪些实际工程问题很多人问我现在做WiFi仿真是不是都应该直接用网络仿真器跑NS-3我的回答是看你要回答什么问题。如果你关注的是TCP吞吐量、路由协议、节点调度那应该去用系统级仿真工具但如果你关注的是“某个距离下能不能听到信号”“2.4GHz多径环境下1Mbps和11Mbps速率谁更稳”那就绕不开物理层的调制解调和误码率这部分恰恰是MATLAB最擅长的。MATLAB仿真802.11b能解决几类实际问题链路预算验证给定发射功率、天线增益、路径损耗通过仿真得到接收端的误码率或丢包率比只算一条静态链路预算更接近真实。扩频增益的直观验证在AWGN信道下扩频前后的误码率理论值不应该改变但抗窄带干扰的能力会显著提升这个结论用仿真图展示出来比看公式有说服力得多。多径灵敏度的扫描2.4GHz室内环境时延扩展通常在50~100ns左右接近802.11b的码片周期约91ns所以多径是实际部署中躲不开的问题。用简单的两径信道模型就能看到误码率曲线出现地板效应这对接入点布放有直接参考价值。为协议级仿真做准备物理层的BER/PER曲线是MAC层CSMA/CA仿真的输入之一。把物理层跑稳了后面加退避算法、RTS/CTS才算有扎实的下层支撑。2. 动手前必须想清楚的802.11b物理层细节2.1 四种速率与调制/扩频方案的对应关系802.11b支持四种物理层速率很多教程只关注速率值没有真正说清速率为什么是1、2、5.5、11这四档。其实这四档对应的是完全不同的调制和编码组合我整理成了一张表速率调制方式扩频方案码片速率数据码元速率备注1 MbpsDBPSKDSSS 11位Barker码11 Mchip/s1 Msymbol/s每个bit直接扩频2 MbpsDQPSKDSSS 11位Barker码11 Mchip/s1 Msymbol/s每个符号携带2 bit5.5 MbpsDQPSKCCKCCK 8码片11 Mchip/s1.375 Msymbol/s每个符号4 bit11 MbpsDQPSKCCKCCK 8码片11 Mchip/s1.375 Msymbol/s每个符号8 bit这里最容易让新手犯迷糊的是“码片速率固定但数据速率不同”这件事。以1Mbps和2Mbps为例两者的射频码片速率都是11Mchip/s带宽都是约22MHz。区别在于1Mbps模式里1个数据bit先被扩展成11个Barker码片然后每个码片做DBPSK调制所以每秒传1M个bit、对应11M个码片2Mbps模式里2个bit映射成1个DQPSK符号再扩展成11个码片所以每秒传1M个符号、对应2Mbit数据、同样是11M个码片。换句话说扩频后的射频符号速率始终是1Msymbol/s数据速率靠调制阶数提高而不是靠提高带宽。2.2 仿真链路的基础参数做MATLAB仿真前我习惯先把参数表写在注释里避免到后面算信噪比时全部乱掉。以下是这次仿真的基础参数参数数值说明载波频率2.4 GHz802.11b工作频段码片速率11 Mchip/s扩频后的码片速率数据调制DBPSK / DQPSK本次重点做DBPSK扩频码Barker码长度11[1 1 1 -1 1 -1 -1 1 -1 -1 1]过采样率16每个码片用16个采样点表示仿真比特数1e5 ~ 1e6根据噪声环境动态调整信道模型AWGN / 两径衰落两径时延可配置为什么把过采样率设为16这是一个工程权衡。过采样率太低比如等于1仿真正弦波和多径延迟时根本表示不了时间上的精细移动太高比如等于64仿真时间一下子拉长噪声计算也容易出错。16是一个够用又不慢的默认值特别是加入多径信道时0.5个码片时延可以用8个采样点表达精度基本够了。2.3 Barker码自相关特性的价值很多人会忽略Barker码在接收端的作用。它其实是一个“自带同步辅助工具”的序列长度为11的Barker码具有一个非常好的性质当它和自己对齐时相关值最大达到11当它偏移一个或多个码片时相关值的绝对值不超过1。这意味着接收机只要做滑动相关就能找到扩频码的起点这就是802.11b能快速完成码片同步的底层原因。在MATLAB仿真里我们通常假设同步已经完成但如果你要写一个完整的接收机Barker码的自相关峰就是你的“定时度量”。这个细节在仿真里很容易被忽略真正做硬件或软件接收机时却至关重要。我建议大家仿真跑通后自己画一下接收信号和本地Barker码在不同偏移下的相关输出你会对“扩频系统抗干扰能力强”这句话有更直接的理解。3. MATLAB仿真链路搭建从比特流到基带波形3.1 发射端调制、扩频和过采样我这次没有用MATLAB Communications Toolbox里的现成调制器而是自己手写发射端。原因是802.11b的扩频细节很固定自己写反而更容易控制每一个中间信号方便中途画波形图检查。下面是核心代码% 发射端参数 numBits 1e5; % 统计比特数 sps 16; % 每码片采样数 barker [1 1 1 -1 1 -1 -1 1 -1 -1 1]; % 11位Barker码 % 随机比特 dataTx randi([0 1], numBits, 1); % DBPSK调制0 - 1, 1 - -1 modSymbols 1 - 2 * dataTx; % 扩频每个符号乘整段Barker码得到 numBits*11 的码片序列 chipSeq reshape(modSymbols * barker, [], 1); % 过采样每个码片重复 sps 次 waveform repelem(chipSeq, sps); % 归一化让每个数据比特的能量等于1 % 每个bit包含11个码片每个码片采样sps个点所以总共11*sps个采样点 waveform waveform / sqrt(11 * sps);最后一行归一化很关键。这里我把每个数据比特对应的采样点数固定为11 * sps把这些采样点上的总能量归一化为1。这样后面加噪声时噪声功率只和Eb/N0相关不会因为过采样率变化导致信噪比算错。这也是我在第5章要重点讲的坑。如果你追求更真实的谱形状可以在这里加一个升余弦脉冲成形滤波器替代简单的repelem。repelem生成的是矩形脉冲频谱旁瓣很高不符合802.11b的频谱模板。不过做误码率仿真时矩形脉冲对AWGN下的结果没有本质影响所以我一般先用它跑通链路再改成rcosdesign观察波形和频谱。3.2 AWGN信道与两径衰落信道信道部分我做了两层第一层是AWGN第二层是两径衰落。AWGN的噪声方差设置为EbNoVec 0:2:12; % Eb/N0单位dB for idx 1:length(EbNoVec) EbNoLin 10^(EbNoVec(idx)/10); noiseStd sqrt(1 / (2 * EbNoLin)); noise noiseStd * randn(size(waveform)); rxWaveform waveform noise; % ... 接收与统计BER ... end这里噪声方差用了1/(2*EbNoLin)原因很简单发射波形已经按每比特能量归一化为1所以Eb1那么Eb/N0 1/N0 EbNoLin得到N0 1/EbNoLin复基带噪声的实部和虚部各占N0/2因此单路噪声标准差是sqrt(N0/2)。如果只做实信号仿真写sqrt(1/(2*EbNoLin))就是对的。两径信道我一般这样模型化% 两径信道主径 延迟径 delayInSamples sps; % 延迟1个码片 attenuation 0.5; % 延迟径幅度 rxWaveform waveform attenuation * [zeros(delayInSamples,1); waveform(1:end-delayInSamples)]; rxWaveform rxWaveform noise;这个模型虽然简单但只要把延迟时间设置成接近室内环境的时延扩展就能看到频率选择性衰落对802.11b的影响。这里有一个容易犯的错误如果用zeros直接补头但接收端没有把延迟径造成的干扰纳入考虑BER会表现得很不合理。实际做的时候我会把主径信号和延迟径信号分别画出来确认相对时延和幅度都符合预期再继续。3.3 接收端相关解扩和判决接收端可以按照“匹配滤波 - 解扩 - 判决”的顺序做。由于发射端用的是矩形脉冲匹配滤波可以简化成对一个码片内的sp个采样点求平均。之后把11个码片和本地Barker码做内积得到的软值符号就是判决结果% 接收端处理 rxChipMatrix reshape(rxWaveform, sps, []); % 每一列是一个码片 rxChipAvg mean(rxChipMatrix, 1); % 码片级匹配滤波 % 多少个完整bit numBitsRec floor(length(rxChipAvg) / 11); rxSoft zeros(numBitsRec, 1); for k 1:numBitsRec segment rxChipAvg((k-1)*111 : k*11); rxSoft(k) sum(segment .* barker); % 相关解扩 end % 判决软值是负号 - 原始比特为1正号 - 原始比特为0 dataRx double(rxSoft 0); dataRx dataRx(1:numBits); % 截断多余部分 ber sum(dataTx ~ dataRx) / numBits;这段代码没有包含码片同步和载波同步原因是我们在理想AWGN仿真里默认已知起始位置。如果你要仿真真实突发传输还需要把前导码和SFD加进去用前导码做同步估计。我觉得初学阶段可以先跳过同步集中精力把能量和BER算对等把“信号处理链路”跑熟后再补同步环。4. 仿真结果分析DBPSK与DQPSK的误码率对比4.1 完整的仿真循环结构把上一章的代码拼起来套一个遍历Eb/N0的循环就是一个可用的误码率仿真脚本。为了节省时间我会在低Eb/N0时用更少的比特高Eb/N0时增加比特数否则在高信噪比下统计不到足够多的错误比特BER曲线会出现毛刺numBits 5e4; % 低SNR先用5e4 for idx 1:length(EbNoVec) if EbNoVec(idx) 8 numBits 5e5; % 高SNR增加比特数 end % 发射、加噪声、接收、统计BER end这样做的好处是在高Eb/N0区间仍然能统计到几十到几百个错误比特BER曲线更平滑。如果你只是画一个趋势图可以不用这么严谨但要拿曲线做工程参考这个细节就很重要了。4.2 结果曲线、理论值和我的解读在AWGN信道下DBPSK相干解调的理论误码率是BER 1/2 * erfc(sqrt(Eb/N0))这个公式可以放到MATLAB里直接画理论曲线再和仿真点比较。实测下来只要信噪比计算正确仿真点和理论曲线几乎重合。误差通常出现在Eb/N0高的时候原因是错误比特太少统计波动变大。我再用DQPSK做对比时会特意把差分编码的损耗也体现出来。802.11b里的2Mbps模式使用DQPSK接收端一般做差分检测损失比相干检测大约2~3dB。从这个角度看802.11b选择DBPSK/DQPSK不是随便定的而是为了在2.4GHz这个“噪声和干扰都不太友好”的频段里用简单的非相干或差分接收机稳定可靠地传信号。还有一个经常被初学者忽略的现象在AWGN信道下扩频并不会改变DBPSK的误码率曲线。很多学员做完DSSS仿真后发现BER和没扩频时一样就会怀疑自己代码是不是错了。其实不是扩频的真正收益在抗窄带干扰和多径分辨不在AWGN下的能量效率。你可以在接收端叠加一个单音干扰再对比有无扩频的BER就会看到显著差异。5. 从“能跑通”到“跑得对”仿真链路调试实录5.1 信噪比计算错误是最隐蔽的坑我在第一次写802.11b仿真时BER曲线比理论值差了很多花了两天时间排查最后发现是噪声方差算错了。很多教程在加AWGN时只考虑信号的平均功率为1然后写noiseStd sqrt(1/(2*10^(EbNo/10)))这在我前面第3章的归一化前提下是对的但前提是信号平均功率必须等于每比特能量。问题是如果做过采样后没有把所有采样点归一化到“每比特能量为1”那么实际信号采样幅度就不是1。比如我最初用repelem扩频后直接让Barker码片幅度保持1那每个比特包含11个码片、每码片16个采样每个采样幅度是1总能量就是11*16176而不是1。此时再用噪声方差1/(2*EbNoLin)去加噪声实际Eb/N0就被拉高了BER会显得比理论好好到不真实。解决方式就是我在发射端做的归一化waveform waveform / sqrt(11 * sps);这样每比特总能量严格为1之后再用统一公式加噪声就不会错。这个坑几乎每个人都会踩一次所以我不建议直接复制网上的信噪比公式强烈建议先用sum(waveform.^2)/numBits打印一下每比特能量确认它是1再跑BER。5.2 解扩时的码片极性问题802.11b的Barker码有两种极性定义有的资料说数据0发送Barker码本身数据1发送Barker码的反码也有的资料定义相反。只要发射端和接收端保持一致仿真结果是完全等价的。真正容易出错的是你在解扩时如果直接用原始Barker码做相关但发射端用了反码映射那么所有软值符号都会翻转误码率会接近100%。排查这个问题有个小技巧不要一上来就统计BER先取10个比特把发射端的调制符号、扩频码片和接收端解扩软值逐段打印出来。确认软值正负和原始比特的映射关系一致后再跑大规模仿真。这个“先手算10个bit”的习惯能帮你省下大量调代码的时间。5.3 多径衰落下的“伪BER”问题在第3.2节的两径信道中如果时延不是码片周期的整数倍MATLAB里要精确模拟非整数采样点延迟需要用到插值或频域移位。很多初学者图省事直接把延迟设为比如1.3个采样点然后用round取整。这样表面上代码没有报错实际上造成了一个“伪信道”延迟径和主径的相对相位被随机化了每次跑的BER差别很大甚至出现高SNR下BER不下降的假象。我的做法是要么把时延严格设为整数倍采样周期做一个“准”多径仿真要么使用interp1插值把信号移动到亚采样精度再用频域插入的方法实现真正的分数延迟。初期仿真建议先用整数倍时延至少保证结果可复现。等确认多径机理清楚了再引入更真实的信道抽头系数。5.4 仿真速度与内存的平衡如果numBits1e6、sps16扩频后波形长度是1e6*11*161.76e8个采样点一个double数组占1.4GB内存MATLAB很可能直接卡死。哪怕内存够每遍历一个Eb/N0点都跑一遍这么长序列脚本总时间也会让人崩溃。我常用三个策略分段处理把发射序列切成若干段每段numBitsSeg1e4逐段加噪声并解调最后把所有错误比特累加。降低过采样率如果只做AWGN仿真sps4甚至sps2都够只有需要展示脉冲波形或做多径时才用到16。用矩阵化替代循环接收端解扩那层循环其实可以改成矩阵乘但前提是11个码片在内存里连续排列当序列很长时直接矩阵乘可能比循环更慢所以要测试后决定。6. 往系统级走802.11b仿真还能玩出什么花6.1 把PER和MAC层CSMA/CA结合起来物理层BER仿真的下一步是把误码率转换成误包率PER。802.11b的MAC帧长度通常几百到一千多字节在MATLAB里可以用一个简单公式近似PER 1 - (1 - BER)^(frameBits)。有了不同信噪比下的PER值就可以搭建一个轻量级的系统级仿真来模拟CSMA/CA退避过程。我不是说MATLAB一定比NS-3适合系统级仿真而是如果你已经用MATLAB把物理层跑得很熟再往上叠加MAC层时数据接口是最顺的。你可以用一个时间推进循环让多个节点在不同时刻产生数据包用物理层查表方式决定这个包是否被正确接收。这样虽然不如NS-3严谨但做课程设计或技术方案预研足够了。6.2 多节点干扰与隐藏终端场景2.4GHz频段里真正影响WiFi性能的往往不是热噪声而是其他节点造成的同频干扰。你可以在接收端叠加一个干扰信号比如另一个调制后的802.11b信号观察BER随信干比的变化。这个实验能直观解释为什么“一个屋子里终端太多WiFi就会变卡”。隐藏终端场景还可以把CSMA/CA的RTS/CTS机制加进来两个节点同时向一个接入点发包物理层模型给碰撞一个丢包概率MAC层统计吞吐量。这个方向对做无线网络课程项目的人来说非常实用而且不需要改动物理层核心代码只要在上层加时间调度逻辑。6.3 工具箱与手写代码怎么选MATLAB的Communications Toolbox提供了很多现成函数比如dpskmod、dpskdemod、rcosdesign可以节省码代码时间。WLAN Toolbox则主要面向基于OFDM的802.11a/n/ac/ax对802.11b的DSSS/CCK模式支持相对有限。所以我的建议是802.11b物理层仿真完全值得手写一遍因为代码量不大却能让你彻底吃透扩频和调制的细节等后面做OFDM链路时再用WLAN Toolbox效率会高很多。如果你决定使用工具箱也别忘了检查它的参考波形是否带前导码和Signal字段。工具箱生成的波形通常是为了测试接收机算法不一定完全符合802.11b的时域帧结构。在这个问题上手写代码反而更容易控制。关于802.11b仿真我最想分享的经验是先把AWGN单用户链路跑通再把误差率公式和仿真结果对齐最后才叠加多径和干扰。很多项目做不好不是因为MATLAB不熟而是把信道、同步、信道估计这些因素一股脑堆进去最后分不清是协议问题还是代码问题排查成本非常高。这套“先简单后复杂、先单点后系统”的方法同样适用你后续做802.11a的OFDM仿真甚至延伸到5G NR链路。本文还有配套的精品资源点击获取
