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.PreHESpatialMapping = FALSE;的％空间映射预HE字段cfgHEBase.GuardInterval = 0.8;％保护间隔持续时间cfgHEBase。HELTFType = 4;％HE-LTF压缩模式cfgHEBase.ChannelCoding =“方法”;%信道编码cfgHEBase。MCS = 3;％调制和编码方案cfgHEBase.SpatialMapping =“自定义”;％自定义的波束赋形

cfgNDP = cfgHEBase;cfgNDP.APEPLength = 0;NDP％没有数据cfgNDP。NumSpaceTimeStreams = NumTxAnts;％有关反馈矩阵计算cfgNDP。SpatialMapping ='直接';每个文件携带一个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 =“没有”;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（'HE-数据'，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）...'和= '的位元数num2str（NumBitsPsi）]）;情况下1 NumBitsPsi = 4;为psi的比特数NumBitsPhi = 6;位％号披disp (“端至端仿真与压缩波束成形量化与”);DISP（['的位元数= 'num2str（NumBitsPhi）...'和= '的位元数num2str（NumBitsPsi）]）;否则disp (“端至端仿真与非压缩波束成形”);结束%parfor isnr = 1:numSNR %使用“parfor”来加速模拟对于ISNR = 1：numSNR在每次迭代中设置随机的子流索引，以确保每个迭代使用一组可重复的随机数流= RandStream ('combRecursive'，'种子'，100）;stream.Substream = ISNR;RandStream.setGlobalStream（流）;为每个模拟的信噪比点创建一个AWGN通道实例awgnChannel = comm.AWGNChannel;awgnChannel.NoiseMethod =“信噪比”;awgnChannel。SignalPower = 1 / tgaxChannel.NumReceiveAntennas;%说明噪声能量在零，所以信噪比是定义每％的有效子载波awgnChannel。信噪比=信噪比(isnr) -10 * log10 (ofdmInfo.FFTLength / ofdmInfo.NumTones);为每个SNR点创建一个配置对象实例％模拟。这将使使用PARFORcfgHE = cfgHEBase;％循环来模拟多个包numPacketErrors = 0;numPkt = 1;数据包传输的索引%而numPacketErrors < = maxNumErrors & & numPkt < = maxNumPackets％空数据分组传输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（codeBookSize（IBF））找出波束形成反馈矩阵的量子化角度angidx = bfCompressQuantize（V（：，1：NumSTS，:)，NumBitsPhi，NumBitsPsi）;%从量子化角度计算转向矩阵％波束形成器：％假定在发送量化的角度零延迟%从beamformee到beamformer，转向矩阵为从量子化角度计算得到的，在％波束形成器的数据传输。[~,数控,Nr] =大小(V (1,1: NumSTS:));V = bfDecompress (angidx、Nr、数控、NumBitsPhi NumBitsPsi);结束steeringMat = V (: 1: NumSTS:);%波束形成数据传输psduLength = getPSDULength (cfgHE);％以字节为单位的长度PSDUtxPSDU =兰迪（[0 1]，psduLength * 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）），:);demodSym = wlanHEDemodulate（RXDATA，'HE-数据'，cfgHE）;％试验阶段跟踪％平均单流导频估计超过码元（第二尺寸）pilotEstTrack =平均（pilotEst，2）;demodSym = heCommonPhaseErrorTracking（demodSym，pilotEstTrack，cfgHE）;估计噪声功率据nvar = heNoiseEstimate (demodSym (ofdmInfo.PilotIndices,:,:), pilotEstTrack, cfgHE);从解调的符号和信道％提取数据子载波％估算demodDataSym = demodSym（ofdmInfo.DataIndices，：，:);chanEstData = CHANEST（ofdmInfo.DataIndices，：，:);%均衡和STBC联合[eqDataSym, csi] = heEqualizeCombine (demodDataSym chanEstData,据nvar, cfgHE);％恢复数据rxPSDU = wlanHEDataBitRecover（eqDataSym，nVarEst，CSI cfgHE，“LDPCDecodingMethod”，'分层-BP');％确定是否有任何位是错误的，即一个数据包错误packetError = ~ isequal (txPSDU rxPSDU);numPacketErrors = numPacketErrors + packetError;numPkt = numPkt + 1;结束在信噪比点计算数据包误码率(PER)packetErrorRate（IBF，ISNR）= numPacketErrors /（numPkt-1）;DISP（['MCS'num2str（cfgHE.MCS）''...'SNR'num2str（SNR（ISNR））...“之后完成”num2str (numPkt-1)“包”...'PER：'num2str (packetErrorRate (ibf, isnr))));结束DISP（换行）;结束

数字;线型= [“k o”“B-S”“R *”];semilogy（SNR，packetErrorRate（1，:)，线型（1））;保持在;网格在;包含(“信噪比(dB)”);ylabel ('每');对于IBF = 2：numQuant semilogy（SNR，packetErrorRate（IBF，:)，线型（IBF））;结束dataStr = [字符串（[压缩的波束形成,换行符...'NumBitsPhi = 4，NumBitsPsi = 2'换行符)...串（[压缩的波束形成,换行符...'NumBitsPhi = 6，NumBitsPsi = 4'换行符)...“非压缩的波束成形”];传奇(dataStr);标题(sprintf ('802.11ax波束形成PER为通道%s， %s， %s',cfgHEBase.ChannelBandwidth tgaxChannel.DelayProfile cfgHEBase.ChannelCoding));