％创建格式配置对象cfgDMG = wlanDMGConfig;mcs =“9”;％MCS指定为字符串标量或向量串cfgDMG.PSDULength = 4096;％PSDULength字节

％获取采样率和指定的载波频率fs = wlanSampleRate (cfgDMG);fc = 60 e9;％创建TGay频道对象tgayChan = wlanTGayChannel;tgayChan.SampleRate = FS;tgayChan.CarrierFrequency = FC;tgayChan.Environment =“开放区域热点”;tgayChan.TransmitArray。Size = [4 4];tgayChan。TransmitArrayPosition = [0;0;6);％计tgayChan。TransmitArrayOrientation = [0;270;0);%度tgayChan.ReceiveArray.Size = [4 4];tgayChan.ReceiveArrayPosition = [6;6;1.5];％计tgayChan。ReceiveArrayOrientation = [90;0;0);%度tgayChan.BeamformingMethod =“最大功率雷”;tgayChan。NormalizeImpulseResponses = true;

showEnvironment（tgayChan）;

ppm = 20;驱动CFO的时钟精度(ppm)freqOffset = ppm * 1 e-6 *俱乐部;％载波频率偏移（Hz）的延迟= 500;%采样以延迟波形zeroPadding = 100;％添加尾随零，以允许信道延迟

:0.5:0.5 snrRanges = {-2.0,...% MCS 10.0：1.0：5.0，...% MCS 21.0：1.2：7.0，...% MCS 32.0:1.2:8.0,...% MCS 42.0:1.3:8.5,...% MCS 52.5:1.3:9.0,...% MCS 64.0:1.3:10.5,...% MCS 75.0:1.5:12.5,...% MCS 85.5：1.5：13.0，...% MCS 97.0:1.5:14.5,...% MCS 9.18.0：1.8：17.0，...％MCS 1010.0：2.0：20.0，...％MCS 1112.0：2.0：22.0，...％MCS 1212.0：2.0：22.0，...% MCS 12.114.0:2.0:24.0,...% MCS 12.216.0：2.5：28.5，...％MCS 12.317.0:2.5:29.5,...% MCS 12.417.0:2.5:29.5,...% MCS 12.520.0：2.5：32.5};% MCS 12.6

maxNumErrors = 10;在信噪比点的最大数据包错误数maxNumPackets = 100;在SNR点分组的％最大数量

numSNR = numel（snrRanges {1}）;%的信噪比点数如果~ isstring (mcs)错误(“MCS必须被指定为一个字符串标量或串矢量”);结束numMCS = numel（MCS）;MCS的数量％packetErrorRate = 0 (numMCS numSNR);进行下一代Ngi = 64;%标准中规定的固定GI长度(20.6.3.2.5)validMCS =串（排序（[1:12 9.1 12.1：0.1：12.6]））;对于imcs = 1:numMCS cfgDMG。MCS = MCS (imc);如果~ strcmp (phyType (cfgDMG),“SC”）错误（“本例仅支持DMG SC PHY仿真”万博1manbetx);结束印第安纳州= wlanFieldIndices (cfgDMG);信噪比= snrRanges {mcs (imc) = = validMCS};％SNR点从MCS模拟% parfor isnr = 1:numSNR %使用“parfor”来加速模拟对于ISNR = 1：numSNR使用'for'调试模拟每次迭代％设定随机子指数，以确保每个迭代使用一组可重复的随机数流= RandStream ('combRecursive'，'种子'10);流。Substream = isnr;RandStream.setGlobalStream(流);为每个模拟的信噪比点创建一个AWGN通道实例awgnChannel = comm.AWGNChannel;awgnChannel.NoiseMethod =“信噪比”;awgnChannel。信噪比=信噪比(isnr);设置仿真参数numPacketErrors = 0;numPkt = 1;传输数据包的%索引而numPacketErrors < = maxNumErrors & & numPkt < = maxNumPackets％生成波形包PSDU =兰迪（[0 1]，cfgDMG.PSDULength * 8,1）;txWaveform = wlanWaveformGenerator（PSDU，cfgDMG）;添加延迟和尾随零tx =[0(延迟,1);txWaveform;0 (zeroPadding 1)];通过TGay通道％发送。重置为一个信道%不同的实现每个包。重置(tgayChan);chanOut = tgayChan (tx);%添加噪声rx = awgnChannel (chanOut);%添加首席财务官rx = helperFrequencyOffset (rx, fs freqOffset);％数据包检测阈值= 0.03;％适合低信噪比pktStartOffset = dmgPacketDetect (rx 0阈值);如果isempty (pktStartOffset)%如果为空，未检测到STF;包错误numPacketErrors = numPacketErrors + 1;numPkt = numPkt + 1;继续;％进入下一个循环迭代结束％频率偏移估计和校正月31 = rx (pktStartOffset + (ind.DMGSTF (1): ind.DMGSTF (2)));fOffsetEst = dmgCFOEstimate（STF）;RX = helperFrequencyOffset（RX，FS，-fOffsetEst）;％符号定时和信道估计preamblefield = rx (pktStartOffset + 1: pktStartOffset + ind.DMGHeader (2):);[symbolTimingOffset,陈]= dmgTimingAndChannelEstimate (preamblefield);startOffset = pktStartOffset + symbolTimingOffset;%如果没有足够的样本来解码检测到的数据字段开始，%则假定同步错误和包错误如果（开始偏移+ ind.DMGData（2））>大小（RX，1）= numPacketErrors numPacketErrors + 1;numPkt = numPkt + 1;继续;％进入下一个循环迭代结束使用STF作为重复序列的%噪声估计月31 = rx (pktStartOffset + (ind.DMGSTF (1): ind.DMGSTF (2)));据nvar = dmgSTFNoiseEstimate (31);提取数据字段(忽略第一个GI)RXDATA = RX（开始偏移+（（ind.DMGData（1）+ NGI）：ind.DMGData（2）））;％线性频域均衡rxEqDataBlks = dmgSingleCarrierFDE（RXDATA，CHANEST，nVarEst）;％独有的字相位跟踪rxEqDataBlks = dmgUniqueWordPhaseTracking（rxEqDataBlks）;放弃所有方块的GIrxDataSym = rxEqDataBlks (1: end-Ngi,:);从DMG数据字段中恢复传输的PSDU据nvar dataDecode = wlanDMGDataBitRecover (rxDataSym, cfgDMG);％确定是否有任何位是错误的，即一个数据包错误packetError =任何（biterr（PSDU，dataDecode））;numPacketErrors = numPacketErrors + packetError;numPkt = numPkt + 1;结束在信噪比点计算数据包误码率(PER)packetErrorRate（的IMC，ISNR）= numPacketErrors /（numPkt-1）;DISP（加入（["主持人:"cfgDMG.MCS“信噪比”...num2str（SNR（ISNR））“之后完成的”...num2str (numPkt-1)“数据包，PER：”...num2str（packetErrorRate（的IMC，ISNR））]，“”));结束结束

标记='牛* SD ^ V> ;颜色=“bmcrgbrkymcrgbrkymc”;图;对于imcs = 1:numMCS semilogy(snrRanges{mcs(imcs)==validMCS}，packetErrorRate(imcs，:)。' - '标记物（的IMC）颜色（的IMC）]）;保持在;结束网格在;包含(“信噪比(dB)”);ylabel ('每');dataStr = arrayfun (@ (x) sprintf (“MCS % s”mcs, x)“UniformOutput”，假）;图例（dataStr）;标题（'PER为DMG SC-PHY，TGay通道');