NumTxAnts = 4;百分比的发射天线数量numsts = 2;%时空流的数量NumRxAnts = 2;%接收天线数cfgHEBase = wlanHESUConfig;cfgHEBase。ChannelBandwidth =“CBW20”；%通道带宽cfghebase.numspacetimestreams = numsts;%时空流的数量cfgHEBase。NumTransmitAntennas = NumTxAnts;百分比的发射天线数量cfgHEBase。一个PEPLength = 1e3;%有效负载长度，以字节为单位cfgHEBase。ExtendedRange = false;%不要使用扩展范围格式cfgHEBase。Upper106ToneRU = false;%不要使用高于106音RUcfgHEBase。PreHESpatialMapping = false;% pre-HE字段空间映射cfgHEBase。GuardInterval = 0.8;保护间隔时间cfgHEBase。HELTFType = 4;% HE-LTF压缩模式cfgHEBase。ChannelCoding =“方法”；％信道编码cfghebase.mcs = 3;%调制和编码方案cfgHEBase。SpatialMapping =“自定义”；%自定义波束形成

cfgNDP = cfgHEBase;cfgNDP。一个PEPLength = 0;% NDP没有数据cfgNDP。NumSpaceTimeStreams = NumTxAnts;%用于反馈矩阵计算cfgndp.spatialMapping =.“直接”；％每个Txant携带STS

％创建并配置TGAX通道chanBW = cfgHEBase.ChannelBandwidth;tgaxChannel = wlanTGaxChannel;tgaxChannel。DelayProfile =“b型”；tgaxChannel。NumTransmitAntennas = NumTxAnts;tgaxChannel。NumReceiveAntennas = NumRxAnts;tgaxChannel。TransmitReceiveDistance = 5;% NLOS的距离以米为单位tgaxChannel。ChannelBandwidth = chanBW;tgaxChannel。LargeScaleFadingEffect ='没有任何'；tgaxchannel.normalizechanneloutpuls = false;FS = WLANSAMPLEDE（CFGHEBASE）;tgaxchannel.samplere = fs;

codebooksize = [0 1 inf];

信噪比= 10:2:18;

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

nimitant = numel（码本大小）;numsnr = numel（snr）;数量的SNR点数packeterrorrate = zeros（nimquant，numsnr）;%获取已占用子载波指标和OFDM参数ofdmInfo = wlanHEOFDMInfo (“数据”, cfgHEBase);%从PPDU中提取字段的索引印第安纳州= wlanFieldIndices (cfgHEBase);indSound = wlanFieldIndices (cfgNDP);为IBF = 1：NimQuant开关码本（IBF）%见P802.11ax/D4.1章节9.4.1.64案件0 NumBitsPsi = 2;PSI的百分比数量NumBitsPhi = 4;%的比特数DISP（压缩波束形成量化的端到端仿真）;disp (['phi ='的比特数num2str (NumBitsPhi)．．.' and Number of Bits for psi = 'num2str (NumBitsPsi)]);案件1 NumBitsPsi = 4;PSI的百分比数量numbitsphi = 6;%的比特数DISP（压缩波束形成量化的端到端仿真）;disp (['phi ='的比特数num2str (NumBitsPhi)．．.' and Number of Bits for psi = 'num2str (NumBitsPsi)]);除此以外DISP（“非压缩波束形成的端到端仿真”）;结束％partsnr = 1：numsnr％使用'par for'加速模拟为isnr = 1: numSNR％迭代设置随机子流索引，以确保每个%迭代使用一组可重复的随机数流= RandStream (“combRecursive”，“种子”, 100);流。Substream = isnr;RandStream.setGlobalStream(流);%考虑噪声能量为零，因此信噪比定义为%活跃副载波packetSNR =信噪比(isnr) -10 * log10 (ofdmInfo.FFTLength / ofdmInfo.NumTones);%为每个SNR点创建一个HE配置对象实例%的模拟。这将允许使用parforcfgHE = cfgHEBase;% Loop模拟多个数据包numPacketErrors = 0;numPkt = 1;传输的数据包％索引而numPacketErrors < = maxNumErrors & & numPkt < = maxNumPackets%空数据包传输tx = wlanWaveformGenerator ([], cfgNDP);添加尾随零以允许通道延迟txPad = [tx;0 (50, cfgNDP.NumTransmitAntennas)];%通过衰落的室内TGax信道重置(tgaxChannel);为不同的实现重置通道rx = tgaxChannel (txPad);％通过AWGN通道传递波形rx = awgn (rx, packetSNR);在波束形成时计算转向矩阵V = heUserBeamformingFeedback (rx, cfgNDP,真);如果isempty（v）%用户反馈失败，数据包错误numpacketerrors = numpacketerrors + 1;numPkt = numPkt + 1;继续；%进入下一个循环迭代结束如果~ isinf (codeBookSize (ibf))％找到波束形成反馈矩阵的量化角度angidx = bfCompressQuantize (V (: 1: NumSTS,:), NumBitsPhi, NumBitsPsi);％从量化角度计算转向矩阵% beamformer:%假设发射量化角的延迟为零从波束形成器到波束形成器，转向矩阵为％由量化角度计算并用于%波束形成器数据传输。[〜，nc，nr] =大小（v（1,1：numsts，:)）;v = bfdecpress（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;零（50，cfghe.numtransmitantannas）];%通过衰落的室内TGax信道rx = tgaxChannel (txPad);％通过AWGN通道传递波形rx = awgn (rx, packetSNR);包检测并确定粗包偏移量chanBW coarsePktOffset = wlanPacketDetect (rx);如果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 = 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 = 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);导频相位跟踪%符号上的单流试验估计平均值(2维)pilotEstTrack =意味着(飞行员,2);demodSym = heCommonPhaseErrorTracking (demodSym pilotEstTrack cfgHE);估算HE领域的噪声功率nvarest = henoiseestimate（remodsym（Ofdminfo.pilotindices，：:)，pilotesttrack，cfghe）;从解调符号和信道中提取数据子载波%的估计demodDataSym = demodSym (ofdmInfo.DataIndices,:,);chanEstData =成龙(ofdmInfo.DataIndices,:,);％均衡和STBC组合[eqDataSym, csi] = heEqualizeCombine (demodDataSym chanEstData,据nvar, cfgHE);%恢复数据据nvar rxPSDU = wlanHEDataBitRecover (eqDataSym, csi, cfgHE,“LDPCDecodingMethod”，“layered-bp”）;%确定是否有任何位错误，即数据包错误Packeterror =〜Isequal（TXPSDU，RXPSDU）;numpacketerrors = numpacketerrors + packeterror;numPkt = numPkt + 1;结束％计算SNR点的数据包错误率（每个）packetErrorRate (ibf isnr) = numPacketErrors / (numPkt-1);disp ([“主持人”num2str (cfgHE.MCS)”、“．．.“信噪比”num2str(信噪比(isnr))．．.“完成后”num2str (numPkt-1)“包”,．．.“每:”num2str（packeterrorrate（IBF，ISNR））]）;结束disp(换行符);结束

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