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

%返回发送sta的索引。计算RUToneSetIndex并为所有sta启动空间时间流。[txSTAIndex, ruToneSetIndexPerSTA startingSTSNumPerSTA] = heTBNDPMappingParams (chanBW、numSTAs multiplexingFlag, startingAID,沉静的);numTxSTAs =元素个数(txSTAIndex);预定发送的STAscfgSTA =细胞(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。NumReceiveAntennas = numRx;tgaxBase。NormalizeChannelOutputs = false;tgaxBase。PathGainsOutputPort = true;

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

%处理信噪比点ofdmInfo = wlanHEOFDMInfo (“HE-LTF”, cfgNDP);numSNR =元素个数(snrRange);%信噪比点数misdetectionProbability = 0 (numTxSTAs numSNR);印第安纳州= wlanFieldIndices (cfgNDP);所有STAs都一样%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{你});%通过随机的TGax通道通过波形channelIdx = chPermutations (u);重置(tgax {channelIdx});%新通道实现[rxSTA h] = tgax {channelIdx} ([txSTA;0(50、大小(txSTA 2))));执行完美的信道延迟估计以找到开始%的包chDelay (u) = channelDelay (h, chFilterCoefficients);%合并上行波形从所有sta成一个波形rxWaveform = rxWaveform + rxSTA;结束%使用最大通道延迟稍后进行同步%在所有通道之间作为数据包的开始pktOffset = max (chDelay (u));%数据包启动索引%将波形通过AWGN通道。占噪音未使用子载波能量%。snrVal = snrRange (isnr) -10 * log10 (ofdmInfo.FFTLength / ofdmInfo.NumTones);rxWaveform = awgn (rxWaveform snrVal);%上行处理(AP)rxHELTF = rxWaveform (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 (“misdetection概率”）;dataStr = arrayfun (@ (x) sprintf (“沉静的% d ', x)的(txSTAIndex),“UniformOutput”、假);标题(sprintf ('%s，通道%s， %d-by-%d', delayProfile chanBW numTx numRx));传奇(dataStr);网格在；