numtxants = 4;%发送天线数numsts = 2;％空间流的数量numnrxants = 2;%接收天线数cfghebase = wlanhauconfig;cfghebase.channelbandwidth =“CBW20”；%信道带宽cfgHEBase。NumSpaceTimeStreams = NumSTS;％空间流的数量cfghebase.numtransmitantennas = numtxants;%发送天线数cfgHEBase.APEPLength=1e3；%有效负载长度，以字节为单位cfgHEBase。ExtendedRange = false;％不使用扩展范围格式cfgHEBase。Upper106ToneRU = false;％不要使用upper 106 tone rucfgHEBase.PreHESpatialMapping=false；% pre-HE字段空间映射cfgHEBase.GuardInterval=0.8；%保护间隔持续时间cfghebase.heltftype = 4;% HE-LTF压缩模式cfgHEBase.ChannelCoding='ldpc'；%信道编码cfgHEBase。MCS = 3;％调制和编码方案cfghebase.spatialMapping =.“自定义”；%自定义波束形成

cfgNDP = cfgHEBase;cfgNDP。A.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。NormalizeChannelOutputs = false;fs = wlanSampleRate (cfgHEBase);tgaxChannel。SampleRate = f;

codeBookSize = [0 1 Inf];

信噪比=10:2:18；

maxnumerrors = 10;％SNR点的最大分组错误数maxNumPackets = 100;％SNR点的最大数据包数

numQuant =元素个数(codeBookSize);numSNR =元素个数(信噪比);数量的SNR点数packetErrorRate = 0 (numQuant numSNR);％获取占用子载波指数和OFDM参数ofdmInfo = wlanHEOFDMInfo (“数据”，cfgHEBase）；％索引从PPDU中提取字段ind=WLANFIELDICES（cfgHEBase）；indSound=WLANFIELDICES（cfgNDP）；对于IBF = 1：NimQuant转变codeBookSize (ibf)％见P802.11AX / D4.1第9.4.1.64节情况下0磅/平方英寸=2；PSI的百分比数量NumBitsPhi = 4;%phi的位数disp ('使用压缩波束形成量化的端到端模拟')；disp(['phi ='的比特数num2str (NumBitsPhi)......'PSI ='的比特数num2str（NumBitsPsi）]；情况下1 numbitspsi = 4;PSI的百分比数量NumBitsPhi = 6;%phi的位数disp ('使用压缩波束形成量化的端到端模拟')；disp(['phi ='的比特数num2str (NumBitsPhi)......'PSI ='的比特数num2str（NumBitsPsi）]；否则disp (“非压缩波束形成的端到端仿真”）;终止%parfor isnr = 1:numSNR %使用“parfor”加速模拟对于isnr=1:numSNR%每次迭代设置随机子流索引，以确保每个％迭代使用可重复的随机数Stream = RANDSTREAM（“组合递归”,“种子”, 100);流。Substream = isnr;RandStream.setGlobalStream(流);％ocket用于空中的噪声能量，因此每次定义SNR%有源副载波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;零（50，cfgndp.numtransmitantennas）]];%通过衰落的室内TGax信道重置(tgaxChannel);用于不同实现的％重置通道rx=tgaxChannel（txPad）；%将波形通过AWGN通道rx=awgn（rx，packetSNR）；在波束形成时计算转向矩阵v = heuserbeamformingfeedback（rx，cfgndp，true）;如果isempty (V)%用户反馈失败，数据包错误numPacketErrors = numPacketErrors + 1;numPkt = numPkt + 1;继续；%转到下一个循环迭代终止如果~isinf（代码本大小（ibf））求波束形成反馈矩阵的量子化角度angidx=bfCompressQuantize（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]，psduLength*8,1）；%生成随机PSDUcfgHE。SpatialMappingMatrix = steeringMat;tx = wlanWaveformGenerator (txPSDU cfgHE);%添加尾随零以允许信道延迟txpad = [tx;零（50，cfghe.numtransmitantannas）];%通过衰落的室内TGax信道rx=tgaxChannel（txPad）；%将波形通过AWGN通道rx=awgn（rx，packetSNR）；％数据包检测并确定粗分组偏移量粗糙度PKTOFFSET=WLANPACKETDERTECT（rx，chanBW）；如果isempty（粗pktoffset）%如果为空没有检测到L-STF;包错误numPacketErrors = numPacketErrors + 1;numPkt = numPkt + 1;继续；%转到下一个循环迭代终止％提取L-STF并执行粗略频率偏移校正lstf = rx（coarsepktoffset +（ind.lstf（1）：Ind.lstf（2）），:);coarsfreqoff = wlancoarsecfoestimate（lstf，chanbw）;rx = helperfrequencyOffset（RX，FS，-CoarseFreqoff）;%提取非ht字段并确定精细数据包偏移量nonhtfields = rx（coarsepktoffset +（ind.lstf（1）：Ind.lsig（2）），:);finepktoffset = WlanyMboltimingEstime（非Hantfields，Chanbw）;确定最后的数据包偏移量pktOffset = coarsePktOffset + finePktOffset;％如果数据包检测到从预期延迟范围介绍%信道建模；分组错误如果pktOffset>50 numPacketErrors = numPacketErrors+1;numPkt = numPkt + 1;继续；%转到下一个循环迭代终止％提取L-LTF并执行细频偏移校正rxLLTF=rx（pktofset+（ind.LLTF（1）：ind.LLTF（2）），：；fineFreqOff=wlanFineCFOEstimate（rxLLTF，chanBW）；rx=helperfrequencycoffset（rx，fs，-fineFreqOff）；％HE-LTF解调和信道估计rxHELTF=rx（pktofset+（ind.HELTF（1）：ind.HELTF（2）），：）；heltfdem=wlanHEDemodulate（rxHELTF，“HE-LTF”, cfgHE);(陈,飞行员)= heLTFChannelEstimate (heltfDemod cfgHE);％数据解调rxData = rx (pktOffset + (ind.HEData (1): ind.HEData (2)),:);demodSym = wlanHEDemodulate (rxData,“数据”, cfgHE);%导频相位跟踪%符号上的单流试验估计平均值(2维)pilotEstTrack=平均值（pilotEst，2）；DemoSym=共相位误差跟踪（DemoSym，pilotEstTrack，cfgHE）；估计他领域的估计噪音力量据nvar = heNoiseEstimate (demodSym (ofdmInfo.PilotIndices,:,:), pilotEstTrack, cfgHE);从解调符号和信道中提取数据子载波%估计demodatasym=demodSym（ofdmInfo.dataindex，：，：）；chanEstData=chanEst（ofdmInfo.dataindex，：，：）；％均衡和STBC组合[eqdatasym，csi] = heequalizeCombine（Demoddatasym，Chanestdata，Nvarest，Cfghe）;%恢复数据rxPSDU=wlanHEDataBitRecover（eqDataSym、nVarEst、csi、cfgHE、，'ldpcdecodingmethod',“分层bp”）;%确定是否有任何位出错，即数据包错误packetError = ~ isequal (txPSDU rxPSDU);numPacketErrors = numPacketErrors + packetError;numPkt = numPkt + 1;终止%计算信噪比点的包错误率包装错误率（ibf，isnr）=numPacketErrors/（numPkt-1）；disp([“MCS”num2str（cfgHE.MCS）”、“......“信噪比”num2str(信噪比(isnr))......“完成后”num2str（numpkt-1）“包”,......“每:”num2str (packetErrorRate (ibf, isnr))));终止disp（新线）；终止

图;线型= [“k o”“b-s”“r - *”];semilogy(信噪比、packetErrorRate(1:),线型(1));持有在；网格在；Xlabel（“信噪比(dB)”）;ylabel（“每”）;对于ibf = 2:numQuant semilogy(snr,packetErrorRate(ibf，:)，lineTypes(ibf));终止dataStr=[string([压缩的波束形成,换行符......‘NumBitsPhi=4，NumBitsPsi=2’新队]）......字符串([压缩的波束形成,换行符......‘NumBitsPhi=6，NumBitsPsi=4’新队]）......“非压缩波束形成”];传奇（数据组）;标题（Sprintf（'802.11ax波束形成PER for Channel %s， %s， %s'，tgaxchannel.delayprofile，cfghebase.channelbandwidth，cfghebase.channelcoding））;