cfgHE = wlanHESUConfig;cfgHE。ChannelBandwidth ='CBW20'；％ 信道带宽cfghe.numspacetimestreams = 2;%时空流数量cfgHE。NumTransmitAntennas = 2;%发射天线数cfgHE。一个PEPLength = 1e3;%有效负载长度(以字节为单位)cfghe.extendedRange = false;%不要使用扩展范围格式cfgHE。Upper106ToneRU = false;％不要使用upper 106 tone rucfgHE。PreHESpatialMapping = false;% pre-HE字段的空间映射cfgHE。GuardInterval = 0.8;％保护间隔持续时间cfgHE。HELTFType = 4;% HE-LTF压缩模式cfghe.channelcoding ='LDPC'；%信道编码cfghe.mcs = 3;调制和编码方案

%创建并配置TGax通道chanbw = cfghe.channelbandwidth;tgaxchannel = wlantgaxchannel;tgaxchannel.delayprofile =“b型”；tgaxchannel.numtransmitantennas = cfghe.numtransmitantennas;tgaxchannel.numreceiveantennas = 2;tgaxchannel.transmitreceivedistance = 5;% NLOS距离(米)TGAXCHANNEL.CHANNELBANDWIDTH = CHANBW;TGAXCHANNEL.LARGESCALEFADEFFECT =“没有”；tgaxChannel。NormalizeChannelOutputs = false;fs = wlanSampleRate (cfgHE);tgaxChannel。SampleRate = f;

信噪比= 10:5:35;

maxNumErrors = 10;％SNR点处的最大数据包错误数maxNumPackets = 100;％SNR点的最大数据包数

numSNR =元素个数(信噪比);％SNR点数numSNR packetErrorRate = 0 (1);％获得占用子载波指数和OFDM参数ofdmInfo = wlanHEOFDMInfo (“数据”, cfgHE);从PPDU中提取字段的％索引印第安纳州= wlanFieldIndices (cfgHE);%parfor isnr = 1:numSNR %使用'parfor'来加速模拟为ISNR = 1：NUMSNR设置每个迭代的随机子流索引，以确保每个％迭代使用可重复的随机数Stream = randstream（“combRecursive”那'种子', 99);流。Substream = isnr;RandStream.setGlobalStream(流);％ocket用于null中的噪声能量，因此每次定义SNR％活动子载波packetSNR =信噪比(isnr) -10 * log10 (ofdmInfo.FFTLength / ofdmInfo.NumTones);%循环，模拟多个数据包numpacketerrors = 0;numpkt = 1;%发送的报文索引而numpacketerrors <= maxnumerrors && numpkt <= maxnumpackets%生成随机PSDU的数据包psduLength = getPSDULength (cfgHE);% PSDU长度(字节)txpsdu = randi（[0 1]，psdulength * 8,1）;Tx = WlanWaveFormGenerator（TXPSDU，CFGHE）;％添加尾随零以允许通道延迟txPad = [tx;0 (50, cfgHE.NumTransmitAntennas)];%通过褪色室内TGax通道重置(tgaxChannel);%为不同的实现方式重置通道rx = tgaxChannel (txPad);%将波形通过AWGN通道rx = awgn（rx，packetsnr）;%检测和确定粗包偏移coarsepktoffset = wlanpacketdetect（rx，chanbw）;如果isempty (coarsePktOffset)%如果为空，没有检测到L-STF;数据包错误numPacketErrors = numPacketErrors + 1;numPkt = numPkt + 1;继续；％转到下一个循环迭代结束％提取L-STF并执行粗略频率偏移校正lstf = rx（coarsepktoffset +（Ind.lstf（1）：Ind.lstf（2）），:);coarsefreqoff = wlancoarsecfoestimate（LSTF，Chanbw）;rx = helperfrequencyoffset（rx，fs，-coarsefreqoff）;%提取非ht字段，确定细包偏移量nonhtfields = rx（coarsepktoffset +（ind.lstf（1）：Ind.lsig（2）），:);finepktoffset = Wlanymboltimingestimate（非Hantfields，Chanbw）;%决定最终数据包偏移量pktOffset = coarsePktOffset + finePktOffset;％如果数据包检测到来自的预期延迟范围％频道建模;数据包错误如果pktOffset>50 numPacketErrors = numPacketErrors+1;numPkt = numPkt + 1;继续；％转到下一个循环迭代结束％提取L-LTF并执行精细频率偏移校正rxLLTF = rx (pktOffset + (ind.LLTF (1): ind.LLTF (2)),:);fineFreqOff = wlanFineCFOEstimate (rxLLTF chanBW);rx = helperFrequencyOffset (rx, fs -fineFreqOff);％HE-LTF解调和信道估计rxheltf = rx（pktoffset +（Ind.heltf（1）：Ind.Heltf（2）），:);heltfdemod = wlanhedemodulate（rxheltf，'he-ltf', cfgHE);(陈,飞行员)= heLTFChannelEstimate (heltfDemod cfgHE);％数据解调rxData = rx (pktOffset + (ind.HEData (1): ind.HEData (2)),:);demodSym = wlanHEDemodulate (rxData,“数据”, cfgHE);％导频阶段跟踪Demodsym = HeCommOnphaseErrortracking（Demodsym，Chanest，Cfghe）;估计他领域的噪音力量nvarest = henoiseestimate（remodsym（Ofdminfo.pilotindices，：，:)，pilotest，cfghe）;从解调符号和信道中提取数据子载波％ 估计demodDataSym = demodSym (ofdmInfo.DataIndices,:,);chanEstData =成龙(ofdmInfo.DataIndices,:,);％均衡和STBC组合[eqDataSym, csi] = heEqualizeCombine (demodDataSym chanEstData,据nvar, cfgHE);％恢复数据rxpsdu = wlanhedatabitrocover（eqdatasym，nvarest，csi，cfghe，'ldpcdecodingmethod'那'分层BP'）;％确定是否有错误，即数据包错误packetError = ~ isequal (txPSDU rxPSDU);numPacketErrors = numPacketErrors + packetError;numPkt = numPkt + 1;结束%计算SNR点的包错误率packetErrorRate (isnr) = numPacketErrors / (numPkt-1);disp (['MCS'num2str (cfgHE.MCS)”、“......'snr'num2str(信噪比(isnr))......'后完成'num2str（numpkt-1）“包”,......“每:”num2str (packetErrorRate (isnr))));结束

数字;半机（SNR，Packeterrorrate，' -  *'）;持有在；网格在；包含(“信噪比(dB)”）;ylabel（“每”）;dataStr = arrayfun (@ (x) sprintf ('MCS％D'，x），cfghe.mcs，'统一输出'、假);传奇(dataStr);标题(sprintf (“每个人都是频道％s，％s，％s，psdulength：％d', cfgHE.ChannelBandwidth tgaxChannel.DelayProfile cfgHE.ChannelCoding cfgHE.APEPLength));