设置传输参数numSTAs =4；%上行链路sta个数稳重的=[1 19 2 4]；% STA关联ID，分配给每个与AP关联的STAmultiplexingFlag =[0 1 0 0]；%每个STA在NFRP触发框中有信号feedbackStatus =[1 0 10 0]；每个STA发出的资源请求信息startingAID =0；% NFRP型触发框内有信号设置仿真参数numPackets =50；%要模拟的数据包数snrRange =2:2:6；信噪比% (dB)chanBW =“CBW20”；通道带宽%numTx =1；%假设所有sta的发射天线数量相同numRx =4；%接收(AP)天线数

返回发送sta的索引。计算所有sta的RUToneSetIndex和起始时空流。[txSTAIndex,ruToneSetIndexPerSTA,startingSTSNumPerSTA] = heTBNDPMappingParams(chanBW,numSTAs,multiplexingFlag,startingAID,STAID);numTxSTAs = numel(txSTAIndex);%计划发送的STAscfgSTA = cell(1,numTxSTAs);生成配置对象并为所有sta设置反馈状态属性cfgBase = wlanHETBConfig(“ChannelBandwidth”chanBW,“NumTransmitAntennas”numTx,“SpatialMapping”，“傅里叶”）;cfgNDP = getNDPFeedbackConfiguration(cfgBase);如果numSTAs ~ =元素个数(feedbackStatus)错误(FeedbackStatus中的元素数量必须等于sta的数量）;结束为u = 1: numTxSTAs cfgNDP。RUToneSetIndex = ruToneSetIndexPerSTA(u);cfgNDP。年代tartingSpaceTimeStream = startingSTSNumPerSTA(u); cfgNDP.FeedbackStatus = feedbackStatus(txSTAIndex(u)); cfgSTA{u} = cfgNDP;结束

delayProfile =“模型”；% TGax通道多路径延迟配置文件tgaxBase = wlanTGaxChannel;tgaxBase。DelayProfile = DelayProfile;tgaxBase。年代ampleRate = wlanSampleRate(cfgSTA{1}); tgaxBase.TransmissionDirection =“上行”；tgaxBase。TransmitReceiveDistance = 10;tgaxBase。ChannelBandwidth = chanBW;tgaxBase。numreceiveantenna = numRx;tgaxBase。NormalizeChannelOutputs = false;tgaxBase。PathGainsOutputPort = true;

单元格数组存储通道对象，每个STA一个。tgax = cell(1,numTxSTAs);为u=1:numTxSTAs tgax{u} =克隆(tgaxBase);tgax{你}。numtransmitantenna = numTx;tgax{你}。UserIndex = u;结束chInfo = info(tgaxBase);chFilterCoefficients = chInfo.ChannelFilterCoefficients;%信道滤波器系数

%处理信噪比点ofdmInfo = wlanHEOFDMInfo(“HE-LTF”, cfgNDP);numSNR = nummel (snrRange);信噪比百分比misdetectionProbability = 0 (numTxSTAs,numSNR);ind = wlanFieldIndices(cfgNDP);%所有sta相同%parfor isnr=1:numSNR %使用'parfor'加快模拟速度为isnr = 1: numSNR每次迭代设置随机子流索引，以确保每个子流的索引迭代使用一组可重复的随机数流= RandStream(“combRecursive”，“种子”, 0);流。年代ubstream = isnr; RandStream.setGlobalStream(stream); rxFeedbackStatus = zeros(numPackets,numTxSTAs); chDelay = zeros(1,numTxSTAs);为pktIdx=1:numPackets rx波形= 0;为每个包随机生成信道实现通道的UserIndex属性。这假设了所有的sta%拥有相同数量的发射天线。chPermutations = randperm(numTxSTAs);为u = 1: numTxSTAs为每个STA生成HE TB反馈NDP波形txSTA = wlanWaveformGenerator([]，cfgSTA{u});%通过随机TGax通道传递波形channelIdx = chPermutations(u);重置(tgax {channelIdx});%新渠道实现[rxSTA,h] = tgax{channelIdx}([txSTA;0(50、大小(txSTA 2))));%执行完美的信道延迟估计以找到起始点%数据包chDelay(u) = channelDelay(h,chFilterCoefficients);将所有sta的上行波形合并为一个波形rx波形= rx波形+rxSTA;结束%稍后使用最大通道延迟同步%作为包的开始pktOffset = max(chDelay(u));包开始索引通过AWGN通道传递波形。考虑噪音%未使用子载波能量。snrVal = snrRange(isnr)-10*log10(ofdmInfo.FFTLength/ofdmInfo.NumTones);rx波形= awgn(rx波形，snrVal);%上行链路处理(AP)rxHELTF = rx波形(pktOffset+(ind.HELTF(1):ind.HELTF(2))，:);heltfDemod = wlanHEDemodulate(rxHELTF，“HE-LTF”, cfgNDP.GuardInterval chanBW cfgNDP.HELTFType);恢复所有sta的反馈状态为u=1:numTxSTAs rxFeedbackStatus(pktIdx,u) = wlanHETBNDPFeedbackStatus(heltfDemod,cfgSTA{u});结束结束每个STA的漏检概率%misdetectionProbability(:，isnr) = 1-sum(rxFeedbackStatus==feedbackStatus(txSTAIndex))/numPackets;disp ([“信噪比”num2str (snrRange (isnr))“完成”]);结束

标记=“牛* sd ^ v > < ph值+ '；numMarkers =数字(标记);为u = 1: numTxSTAs semilogy (snrRange misdetectionProbability (u,:), (“- - -”标记(mod (u-1 numMarkers) + 1)]);持有在；结束包含(“信噪比(dB)”）;ylabel (“漏检概率”）;dataStr = arrayfun(@(x)sprintf(“沉静的% d ', x)的(txSTAIndex),“UniformOutput”、假);标题(sprintf ('%s，通道%s， %d-by-%d', delayProfile chanBW numTx numRx));传奇(dataStr);网格在；