NumTxAnts = 4;%发送天线数NumSTS = 2;%时空流的数量NumRxAnts = 2;%接收天线数cfgHEBase = wlanHESUConfig;cfgHEBase。ChannelBandwidth ='CBW20'；％ 信道带宽cfgHEBase。NumSpaceTimeStreams = NumSTS;%时空流的数量cfgHEBase。NumTransmitAntennas = NumTxAnts;%发送天线数cfghebase.apeplength = 1e3;％有效载荷长度以字节为单位cfgHEBase。ExtendedRange = false;%不要使用扩展范围格式cfgHEBase。Upper106ToneRU = false;%不要使用高于106音RUcfghebase.prehepatialmapping = false;预先田地的％空间映射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 = numnrxants;tgaxchannel.transmitreceivedistance = 5;% NLOS的距离以米为单位TGAXCHANNEL.CHANNELBANDWIDTH = CHANBW;tgaxchannel.largescalefadingeffect =“没有”；fs = wlanSampleRate (cfgHEBase);tgaxChannel。SampleRate = f;

codeBookSize = [0 1 Inf];

SNR = 10：2：18;

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

numQuant =元素个数(codeBookSize);numSNR =元素个数(信噪比);%信噪比点数packetErrorRate = 0 (numQuant numSNR);%获取已占用子载波指标和OFDM参数Ofdminfo = wlanheofdminfo（'他是数据'，cfghebase）;%从PPDU中提取字段的索引IND = WLANFIELDINDIDICES（CFGHEBASE）;Indsound = WlanFieldIndices（CFGNDP）;为了ibf = 1: numQuant开关codeBookSize (ibf)%见P802.11ax/D4.1章节9.4.1.64情况下0 numbitspsi = 2;% psi的位数NumBitsPhi = 4;PHI的百分比数量disp ('具有压缩波束成形量化的端到端模拟'）;DISP（[' =的比特数'num2str（numbitsphi）．..' and Number of Bits for psi = 'num2str（numbitspsi）]）;情况下1 NumBitsPsi = 4;% psi的位数NumBitsPhi = 6;PHI的百分比数量disp ('具有压缩波束成形量化的端到端模拟'）;DISP（[' =的比特数'num2str（numbitsphi）．..' and Number of Bits for psi = 'num2str（numbitspsi）]）;否则disp ('具有非压缩波束形成的端到端模拟'）;结尾%parfor isnr = 1:numSNR %使用“parfor”加速模拟为了ISNR = 1：NUMSNR%每次迭代设置随机子流索引，以确保每个%迭代使用一组可重复的随机数流= RandStream ('幸福父母'那'种子'，100）;Stream.substream = ISNR;randstream.setGlobalstream（Stream）;%创建一个AWGN信道实例每个SNR点模拟Awgnchannel = comm.awgnchannel;awgnchannel.noisemethod =“信噪比”；awgnChannel。SignalPower = 1 / tgaxChannel.NumReceiveAntennas;%考虑噪声能量为零，因此信噪比定义为％活动子载波awgnChannel。信噪比=信噪比(isnr) -10 * log10 (ofdmInfo.FFTLength / ofdmInfo.NumTones);%为每个SNR点创建一个HE配置对象实例％模拟。这将启用可用cfghe = cfghebase;% Loop模拟多个数据包numPacketErrors = 0;numPkt = 1;%发送的数据包索引尽管numPacketErrors < = maxNumErrors & & numPkt < = maxNumPackets％null数据包传输tx = wlanWaveformGenerator ([], cfgNDP);％添加尾随零以允许通道延迟txPad = [tx;0 (50, cfgNDP.NumTransmitAntennas)];％通过褪色室内TGAX通道重置（TGAXCHANNEL）;为不同的实现重置通道rx = tgaxchannel（txpad）;%将波形通过AWGN通道rx = awgnchannel（rx）;％计算波束上的转向矩阵V = heUserBeamformingFeedback (rx, cfgNDP,真);如果isempty (V)％用户反馈失败，数据包错误numPacketErrors = numPacketErrors + 1;numpkt = numpkt + 1;继续；％转到下一个循环迭代结尾如果〜ISINF（码本（IBF））求波束形成反馈矩阵的量子化角度angidx = BfChageSquantize（V（：，1：numsts，:)，numbitsphi，numbitspsi）;从点的量子化角度计算转向矩阵% beamformer:在发送量化角度时，％假设零延迟从波束形成器到波束形成器，转向矩阵为%由量子化角度计算，用于波束形成器的％数据传输。[~,数控,Nr] =大小(V (1,1: NumSTS:));V = bfDecompress (angidx、Nr、数控、NumBitsPhi NumBitsPsi);结尾steeringMat = V (: 1: NumSTS:);%波束形成的数据传输psduLength = getPSDULength (cfgHE);％PSDU长度为字节TXPSDU = RANDI（[0 1]，PSDULENG * 8,1）;％生成随机psducfghe.spatialmappingmatrix = steeringmat;Tx = WlanWaveFormGenerator（TXPSDU，CFGHE）;％添加尾随零以允许通道延迟txPad = [tx;0 (50, cfgHE.NumTransmitAntennas)];％通过褪色室内TGAX通道rx = tgaxchannel（txpad）;%将波形通过AWGN通道rx = awgnchannel（rx）;％数据包检测并确定粗分组偏移量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 = 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）;[Chanest，Pilotest] = HeltfchanneLeStimate（HeltFdemod，CFGHE）;%数据解调rxdata = rx（pktoffset +（ind.hedata（1）：Ind.hedata（2）），:);emodsym = wlanhedemodulate（rxdata，'他是数据'，cfghe）;％导频阶段跟踪%符号上的单流试验估计平均值(2维)pilotesttrack =卑鄙（Pilotest，2）;Demodsym = HeCommOnphaseErrortracking（Demodsym，Pilotesttrack，Cfghe）;估算HE领域的噪声功率据nvar = heNoiseEstimate (demodSym (ofdmInfo.PilotIndices,:,:), pilotEstTrack, cfgHE);从解调符号和信道中提取数据子载波％ 估计emoddatasym = remodsym（Ofdminfo.dataindices，：，:);chanestdata = chanest（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;结尾%计算信噪比点的包错误率Packeterrorrate（IBF，ISNR）= NumPacketerrors /（NumPKT-1）;DISP（['MCS'num2str（cfghe.mcs）”、“．..'snr'num2str(信噪比(isnr))．..'完成后'num2str (numPkt-1)'包，'．..'每：'num2str (packetErrorRate (ibf, isnr))));结尾DISP（甜点）;结尾

数字;linetypes = [“k o”“b-s”“r- *”];半径（SNR，Packeterrorrate（1，:)，Lineetypes（1））;抓住在；网格在；包含(“信噪比(dB)”）;ylabel ('每'）;为了IBF = 2：Nimquant半径（SNR，Packeterrorrate（IBF，:)，LineTypes（IBF））;结尾data口= [string（[压缩的波束形成,换行符．..'numbitsphi = 4，numbitspsi = 2'换行符)．..字符串([压缩的波束形成,换行符．..'numbitsphi = 6，numbitspsi = 4'换行符)．..“非压缩波束形成”];传奇(dataStr);标题(sprintf ('802.11ax波束形成PER for Channel %s， %s， %s', cfgHEBase.ChannelBandwidth tgaxChannel.DelayProfile cfgHEBase.ChannelCoding));