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;%前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 =“没有”；fs = wlanSampleRate (cfgHEBase);tgaxChannel。SampleRate = f;

codeBookSize = [0 1 Inf];

信噪比= 10:2:18;

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

numQuant =元素个数(codeBookSize);numSNR =元素个数(信噪比);%信噪比点数packetErrorRate = 0 (numQuant numSNR);%获取已占用子载波指标和OFDM参数ofdmInfo=wlanHEOFDMInfo(“他需要数据”, cfgHEBase);%从PPDU中提取字段的索引印第安纳州= 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）．..' and Number of Bits for psi = 'num2str (NumBitsPsi)]);情况下1 NumBitsPsi = 4;% psi的位数NumBitsPhi = 6;%的比特数disp (压缩波束形成量化的端到端仿真）;disp ([' =的比特数'num2str（NumBitsPhi）．..' and Number of Bits for psi = 'num2str (NumBitsPsi)]);否则disp (“非压缩波束形成的端到端模拟”）;结束%parfor isnr = 1:numSNR %使用“parfor”加速模拟为isnr = 1: numSNR%每次迭代设置随机子流索引，以确保每个%迭代使用一组可重复的随机数流= RandStream (“combRecursive”，“种子”,100); 子流=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配置对象实例%模拟。这将允许使用parforcfgHE = cfgHEBase;% Loop模拟多个数据包numPacketErrors = 0;numPkt = 1;%发送的数据包索引虽然numPacketErrors < = maxNumErrors & & numPkt < = maxNumPackets%空数据包传输tx = wlanWaveformGenerator ([], cfgNDP);添加尾随零以允许通道延迟txPad = [tx;0 (50, cfgNDP.NumTransmitAntennas)];%通过衰落的室内TGax信道复位（TGAX通道）；为不同的实现重置通道rx = tgaxChannel (txPad);%将波形通过AWGN通道rx = awgnChannel (rx);%在beamformee处计算转向矩阵V = heUserBeamformingFeedback (rx, cfgNDP,真);如果isempty (V)%用户反馈失败，数据包错误numPacketErrors = numPacketErrors + 1;numPkt = numPkt + 1;继续；%进入下一个循环迭代结束如果~ isinf (codeBookSize (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;0 (50, cfgHE.NumTransmitAntennas)];%通过衰落的室内TGax信道rx = tgaxChannel (txPad);%将波形通过AWGN通道rx = awgnChannel (rx);包检测并确定粗包偏移量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=粗pktOffset+精pktOffset；%如果数据包被检测出，且从%通道建模;包错误如果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（heltfdemm，cfgHE）；%数据解调rxData=rx（pktofset+（ind.HEData（1）：ind.HEData（2）），：；demosym=wlanHEDemodulate（rxData，“他需要数据”, cfgHE);导频相位跟踪%符号上的单流试验估计平均值(2维)pilotEstTrack=平均值（pilotEst，2）；DemoSym=共相错误跟踪（DemoSym，PilotTestTrack，cfgHE）；估算HE领域的噪声功率据nvar = heNoiseEstimate (demodSym (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;结束%计算信噪比点的包错误率packetErrorRate (ibf isnr) = numPacketErrors / (numPkt-1);disp ([“主持人”num2str (cfgHE.MCS)”、“．..“信噪比”num2str(信噪比(isnr))．..“完成后”num2str (numPkt-1)“包，”．..“PER:”num2str (packetErrorRate (ibf, isnr))));结束disp(换行符);结束

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