NumTxAnts = 4;%发射天线个数NumSTS = 2;%时空流个数NumRxAnts = 2;%接收天线个数cfgHEBase = wlanHESUConfig;cfgHEBase。ChannelBandwidth =“CBW20”；通道带宽%cfgHEBase。NumSpaceTimeStreams = NumSTS;%时空流个数cfgHEBase。numtransmitantenna = NumTxAnts;%发射天线个数cfgHEBase。APEPLength = 1e3;负载长度(以字节为单位)cfgHEBase。ExtendedRange = false;不要使用扩展范围格式cfgHEBase。Upper106ToneRU = false;不使用106调以上的俄服cfgHEBase。PreHESpatialMapping = false;pre-HE字段的空间映射cfgHEBase。GuardInterval = 3.2;防护间隔时间%cfgHEBase。HELTFType = 4;HE-LTF压缩模式cfgHEBase。ChannelCoding =“方法”；%信道编码cfgHEBase。MCS = 3;调制和编码方案cfgHEBase。SpatialMapping =“自定义”；%自定义波束形成

cfgNDP = cfgHEBase;cfgNDP。APEPLength = 0;% NDP没有数据cfgNDP。NumSpaceTimeStreams = NumTxAnts;用于反馈矩阵计算cfgNDP。SpatialMapping =“直接”；每个TxAnt携带一个STS

创建和配置TGax通道chanBW = cfgHEBase.ChannelBandwidth;tgaxChannel = wlanTGaxChannel;tgaxChannel。DelayProfile =“b型”；tgaxChannel。numtransmitantenna = NumTxAnts;tgaxChannel。numreceiveantenna = NumRxAnts;tgaxChannel。TransmitReceiveDistance = 5;%距离(以米为单位)tgaxChannel。ChannelBandwidth = chanBW;tgaxChannel。LargeScaleFadingEffect =“没有”；tgaxChannel。NormalizeChannelOutputs = false;fs = wlanSampleRate(cfgHEBase);tgaxChannel。SampleRate = fs;

codeBookSize = [0 1 Inf];

信噪比= 10:2:18;

maxNumErrors = 10;在一个信噪比点上的最大错误包数maxNumPackets = 100;一个信噪比点的最大包数

numQuant = numel(codeBookSize);numSNR = numel(snr);信噪比百分比packetErrorRate = 0 (numQuant,numSNR);获取占用子载波指数和OFDM参数ofdmInfo = wlanHEOFDMInfo(“数据”, cfgHEBase);%索引从PPDU中提取字段ind = wlanFieldIndices(cfgHEBase);indSound = wlanFieldIndices(cfgNDP);为ibf = 1:numQuant开关codeBookSize (ibf)参见P802.11ax/D4.1章节9.4.1.64情况下0 NumBitsPsi = 2;% psi的位数NumBitsPhi = 4;的比特数disp (带压缩波束形成量化的端到端模拟）;disp ([=的比特数num2str (NumBitsPhi).．.'和psi = '的位数num2str (NumBitsPsi)]);情况下1 NumBitsPsi = 4;% psi的位数NumBitsPhi = 6;的比特数disp (带压缩波束形成量化的端到端模拟）;disp ([=的比特数num2str (NumBitsPhi).．.'和psi = '的位数num2str (NumBitsPsi)]);否则disp (采用非压缩波束形成的端到端模拟）;结束%parfor isnr = 1:numSNR %使用'parfor'加快模拟速度为isnr = 1:numSNR每次迭代设置随机子流索引，以确保每个子流的索引迭代使用一组可重复的随机数流= RandStream(“combRecursive”，“种子”, 100);流。Substream = isnr;RandStream.setGlobalStream(流);%考虑噪声能量为空，因此信噪比定义为有源子载波百分比packetSNR = snr(isnr)-10*log10(ofdmInfo.FFTLength/ofdmInfo.NumTones);为每个信噪比点创建一个HE配置对象实例%的模拟。这将允许使用parforcfgHE = cfgHEBase;%循环模拟多个数据包numPacketErrors = 0;numPkt = 1;%发送的报文索引而numPacketErrors<=maxNumErrors && numPkt<=maxNumPackets%空数据包传输tx = wlanWaveformGenerator([]，cfgNDP);添加尾随零以允许信道延迟txPad = [tx;0 (50, cfgNDP.NumTransmitAntennas)];通过一个衰落的室内TGax通道重置(tgaxChannel);为不同的实现重置通道rx = txxchannel (txPad);通过AWGN通道传递波形rx = awgn(rx,packetSNR);计算波束范围内的转向矩阵V = heUserBeamformingFeedback(rx,cfgNDP,true);如果isempty (V)%用户反馈失败，数据包错误numPacketErrors = numPacketErrors+1;numPkt = numPkt+1;继续；%转到下一次循环迭代结束如果~ isinf (codeBookSize (ibf))求波束形成反馈矩阵的量子化角度angidx = bfCompressQuantize(V(:，1:NumSTS，:)，NumBitsPhi,NumBitsPsi);从的量子化角度计算转向矩阵% beamformer:%假设传输量子化角度的延迟为零%从波束发生器到波束形成器，转向矩阵为%由量子化角度计算，用于%波束形成器的数据传输。[~，Nc,Nr] = size(V(1,1:NumSTS，:));V = bfDecompress(angidx,Nr,Nc,NumBitsPhi,NumBitsPsi);结束steeringMat = V(:，1:NumSTS，:);%波束形成的数据传输psduLength = getPSDULength(cfgHE);% PSDU长度，单位为字节txPSDU = randi([0 1]，psduLength*8,1);%生成随机PSDUcfgHE。SpatialMappingMatrix = steeringMat;tx = wlanWaveformGenerator(txPSDU,cfgHE);添加尾随零以允许信道延迟txPad = [tx;0 (50, cfgHE.NumTransmitAntennas)];通过一个衰落的室内TGax通道rx = txxchannel (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))，:);粗freqoff = wlan粗secfoestimate (lstf,chanBW);rx = frequencyOffset(rx,fs，-粗freqoff);提取非ht字段并确定细包偏移量nonhtfields = rx(coarsePktOffset+(ind.LSTF(1):ind.LSIG(2))，:);finePktOffset = wlanSymbolTimingEstimate(nonhtfields,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 = frequencyOffset(rx,fs，-fineFreqOff);HE-LTF解调和信道估计rxHELTF = rx(pktOffset+(ind.HELTF(1):ind.HELTF(2))，:);heltfDemod = wlanHEDemodulate(rxHELTF，“HE-LTF”, cfgHE);[chanEst,pilotEst] = heLTFChannelEstimate(heltfDemod,cfgHE);数据解调rxData = rx(pktOffset+(ind.HEData(1):ind.HEData(2))，:);demodSym = wlanHEDemodulate(rxData，“数据”, cfgHE);导频相位跟踪对符号的平均单流导频估计(第2维)pilotEstTrack = mean(pilotEst,2);demodSym = wlanHETrackPilotError(demodSym,pilotEstTrack,cfgHE，“数据”）;估计HE场的噪声功率nVarEst = heNoiseEstimate(demodSym(ofdmInfo.PilotIndices，:，:)，pilotEstTrack,cfgHE);从解调的符号和信道中提取数据子载波%的估计demodDataSym = demodSym(ofdinfo . dataindices，:，:);chanEstData = chanEst(ofdinfo . dataindices，:，:);%平衡和STBC结合[eqDataSym,csi] = heEqualizeCombine(demodDataSym,chanEstData,nVarEst,cfgHE);%恢复数据rxPSDU = wlanHEDataBitRecover(eqDataSym,nVarEst,csi,cfgHE，“LDPCDecodingMethod”，“norm-min-sum”）;确定是否有任何位错误，即数据包错误。packetError = ~isequal(txPSDU,rxPSDU);numPacketErrors = numPacketErrors+packetError;numPkt = numPkt+1;结束计算信噪比点的包错误率(PER)packeterrors (ibf,isnr) = numPacketErrors/(numPkt-1);disp ([“主持人”num2str (cfgHE.MCS)”、“.．.“信噪比”num2str(信噪比(isnr)).．.“完成后”num2str (numPkt-1)“包”,.．.“每:”num2str (packetErrorRate (ibf, isnr))));结束disp(换行符);结束

图;lineTypes = [“k o”“台球”“r - *”];semilogy(信噪比、packetErrorRate(1:),线型(1));持有在；网格在；包含(“信噪比(dB)”）;ylabel (“每”）;为ibf = 2:numQuant semilogy(snr,packetErrorRate(ibf，:)，lineTypes(ibf));结束dataStr = [string([“压缩波束形成”换行符.．.' numbitspi = 4, NumBitsPsi = 2'换行符).．.字符串([“压缩波束形成”换行符.．.' numbitspi = 6, NumBitsPsi = 4'换行符).．.“没有波束形成”];传奇(dataStr);标题(sprintf (用于信道%s， %s， %s的802.11ax波束形成PER ', cfgHEBase.ChannelBandwidth tgaxChannel.DelayProfile cfgHEBase.ChannelCoding));