EbNo = 2:4:10;% Eb/No in dBSPS = 4;%每个符号的样本，必须大于1dataLen = 42;数据长度(以字节为单位)，包括报头、有效载荷和CRCsimMode = {“LE1M”，“LE2M”，“LE500K”，“LE125K”};% PHY模式考虑模拟

maxNumErrors = 100;Eb/No点最大误码数maxNumPackets = 10;每个Eb/No点的最大包数

numMode =数字(simMode);%模式数ber = 0 (numMode,length(EbNo));预分配用于存储误码率结果每=零(numMode，长度(EbNo));预分配用于存储PER结果bitsPerByte = 8;%每字节的比特数固定访问地址理想情况下，该访问地址值应满足蓝牙规范第2.1.2节中指定的%要求。accessAddress = [0 1 1 0 1 0 1 1 0 1 1 1 1 1 1 0 0 ..．.1 0 0 0 1 1 1 1 0 0 0 1]';为iMode = 1:numMode phyMode = simMode{iMode};%根据模式设置信噪比(SNR)点。对于编码的物理(LE500K和LE125K)，包括码率因子在信噪比计算中使用1/2速率FEC编码器。如果任何(strcmp (phyMode, {“LE1M”，“LE2M”})) snrVec = EbNo - 10*log10(sps);其他的codeRate = 1/2;snrVec = EbNo + 10*log10(codeRate) - 10*log10(sps);结束使用'parfor'来加快模拟速度为iSnr = 1:长度(snrVec)使用'for'来调试模拟每次迭代设置随机子流索引，以确保每个子流的索引迭代使用一组可重复的随机数流= RandStream(“combRecursive”，“种子”, 0);流。子流= iSnr;RandStream.setGlobalStream(流);创建一个错误率实例errorRate = com . errorRate (“样本”，“自定义”，.．.“CustomSamples”1: (dataLen * bitsPerByte-1));为损害创建和配置系统对象initImp = helperBLEImpairmentsInit(phyMode,sps);在结构中配置接收方参数rxCfg = struct(Mode=phyMode,SamplesPerSymbol=sps,ChannelIndex=37，.．.DFPacketType =“禁用”AccessAddress = AccessAddress);rxCfg。coarse频率补偿器= com . coarse频率补偿器(调制=“OQPSK”，.．.SampleRate = sps * (1 + strcmp (phyMode,“LE2M”e6)) * 1,.．.SamplesPerSymbol = 2 * sps,.．.FrequencyResolution = 30);rxCfg。PreambleDetector = com .PreambleDetector(检测=“第一”）;初始化错误计算参数[numErrs,perCnt] = deal(0);numPkt = 1;%循环模拟多个数据包而numErrs <= maxNumErrors && numPkt <= maxNumPackets%生成蓝牙LE波形txBits = randi([0 1]，dataLen*bitsPerByte,1，“int8”）;数据位生成百分比channelIndex = randi([0 39]，1,1);每个包的随机通道索引值tx波形= bleWaveformGenerator(txBits，“模式”phyMode,.．.“SamplesPerSymbol”sps,.．.“ChannelIndex”channelIndex,.．.“accessAddress”, accessAddress);定义射频损伤参数inittim .pfo. frequencyoffset = randsrc(1,1，-50e3:10:50e3);% In Hz，最大范围为+/- 150khzinittim .pfo. phaseoffset = randsrc(1,1，-10:5:10);%(度)Initoff = 0.15*sps;%静态定时偏移量步长= 20*1e-6;%定时漂移，最大范围为+/- 50ppminitImp。vdelay = (initoff:步长:initoff+步长*(长度(tx波形)-1))';可变定时偏移量initImp。Dc = 20;w.r.t最大振幅值的百分比%通过射频损伤传递生成的波形txImpairedWfm = helperBLEImpairmentsAddition(tx波形，inititimp);将传输的波形通过AWGN通道rx波形= awgn(txImpairedWfm,snrVec(iSnr));使用实际接收器恢复数据位rxCfg。ChannelIndex = ChannelIndex;[rxBits,recAccessAddress] = helperBLEPracticalReceiver(rx波形，rxCfg);确定BER和PER如果(长度(txBits) ==长度(rxBits)) errors = errorRate(txBits,rxBits);%输出累计错误数ber(iMode,iSnr) = errors(1);累积误码率%currentErrors =错误(2)-数字错误;%当前数据包错误数如果(currentErrors)检查当前数据包是否错误perCnt = perCnt + 1;增加PER计数结束numErrs =错误(2);累积误差百分比numPkt = numPkt + 1;结束结束per(iMode,iSnr) = perCnt/(numPkt-1);disp ([“模式”phyMode”、“.．.'模拟Eb/No = 'num2str (EbNo (iSnr)),“数据库””、“.．.的数量:, num2str (ber(将iSnr)),”、“.．.“每:”, num2str(每(将,iSnr))))结束结束

标记=牛*年代”；颜色=“bmcr”；dataStr = {0 (numMode,1)};为iMode = 1:numMode subplot(2,1,1)， semogy (EbNo,ber(iMode，:).'，[“- - -”标记(将)颜色(将)]);持有在；dataStr(iMode) = simMode(iMode);次要情节(2,1,2),semilogy (EbNo,每(将:)。”“- - -”标记(将)颜色(将)]);持有在；dataStr(iMode) = simMode(iMode);结束次要情节(2,1,1),网格在；包含(“Eb /不(dB)”）;ylabel (“方方面面”）;传奇(dataStr);标题(“射频损伤下蓝牙LE的误码率”）;次要情节(2,1,2)网格在；包含(“Eb /不(dB)”）;ylabel (“每”）;传奇(dataStr);标题(“射频损伤下蓝牙LE的PER”）;

[refBER,refPER, refbno] =交易(零(numMode,1));headerLen = 2;%头长度(以字节为单位)crcLen = 3;%以字节为单位的CRC长度payloadLen = dataLen-headerLen-crcLen;负载长度(以字节为单位)为iMode = 1:numMode [refBER(iMode)，refPER(iMode)， refbno (iMode)] =.．.helperBLEReferenceResults (simMode(将),payloadLen);disp ([“模式”simMode{将}”、“.．.'参考Eb/No = 'num2str (refEbNo(将)),“数据库””、“.．.' ber = 'num2str (refBER(将)),”、“.．.' per = 'num2str (refPER(将)),”、“.．.有效载荷长度为num2str (payloadLen),“字节”]）结束