cfgVHT = wlanVHTConfig;%创建报文配置信道带宽='CBW80'；%80兆赫cfgvht.numtransmitantennas = 1;％一个发射天线cfgvht.numspacetimestreams = 1;%一个时空流cfgvht.stbc = false;%没有STBC所以一个空间流cfgVHT。MCS = 8;%调制:256 - qamcfgVHT.APEPLength=3000；％A-MPDU长度预换页垫以字节为单位

osf=3；%3倍过采样因子numPackets=20；％生成20个数据包idleTime = 10 e-6;%数据包之间10微秒的空闲时间

%创建随机数据；pModule长度以字节为单位savedState = rng (0);%设置随机状态数据= randi([0 1]，cfgVHT.PSDULength*8*numPackets,1);%生成多数据包波形txWaveform = wlanWaveformGenerator(数据、cfgVHT...“过度采样物料”，OSF，“NumPackets”，numPackets，“空闲时间”，空闲时间）;FS = WLANSAMPLEDE（CFGVHT）;波形基带采样率%

modelIQImbalance = true;％设置为true以添加IQ增益和相位不平衡如果modeliq失衡iqGaindB = 1;% IQ增益不平衡在dB中，从范围[-1 1]指定iqPhaseDeg = 1;%IQ相位不平衡度，从范围[-2]指定iqgainlin = db2mag（iqgaindb）;%将增益从dB转换为线性值txWaveform = real(txWaveform) + 1i*imag(txWaveform)*iqGainLin*exp(1j*iqPhaseDeg*pi/180); / /输出％在[5]的公式-1中规定结束

PsAturation = 25;%以dBm为单位的功率放大器饱和功率hpaBackoff = 13;DB中％功率放大器退避%创建并配置无记忆非线性以模拟放大器非线性= comm.MemorylessNonlinearity;非线性。方法=“拉普模式”非线性，平滑度=3；% p参数非线性。LinearGain = -hpaBackoff;非线性。OutputSaturationLevel = db2mag (pSaturation-30);%将模型应用于每个发射天线TxWaveForm =非线性（TxWaveForm）;

thNoise = comm.ThermalNoise (“噪音法”，'噪音'，“SampleRate”fs * osf,“NoiseFigure”6);为i = 1：cfgvht.numtransmitantennas txwaveform（：，i）= thnoise（txwaveform（：，i））;结束

aStop=40；%阻带衰减ofdmInfo=wlanVHTOFDMInfo(“甚高频数据”, cfgVHT);% OFDM参数SCS=fs/ofdmInfo.fft长度；%副载波间距txbw = max (abs (ofdmInfo.ActiveFrequencyIndices)) * 2 * SCS;%占用带宽[L，M]=大鼠（osf）；maxLM=最大值（[L M]）；R=（fs-txbw）/fs；TW=2*R/maxLM；%过渡宽度B = DesignMultirateFIR（L，M，TW，ASTOP）;

firinterp = dsp.FIRRateConverter (M, L, b);rxWaveform = firinterp (txWaveform);

%设置EVM测量值[EVMPerPkt, EVMPerSC] = vhtEVMSetup (cfgVHT);

％获取用于访问时域数据包中的每个字段的索引印第安纳州= wlanFieldIndices (cfgVHT);rxWaveformLength =大小(rxWaveform, 1);pktLength =双(ind.VHTData (2));％定义我们可以检测到的最小数据长度;L-STF的长度%样品minPktLen =双(ind.LSTF (2) -ind.LSTF (1) + 1;%设置测量图[据hSF, hCon, hEVM] = vhtTxSetupPlots (cfgVHT);rmsEVM = 0 (numPackets, 1);pktOffsetStore = 0 (numPackets, 1);rng (savedState);恢复随机状态

pktnum = 0;searchOffset = 0;％在第一个样本开始（没有偏移量）尽管(searchOffset + minPktLen) < = rxWaveformLength％数据包检测pktOffset=wlanPacketDetect（rx波形，cfgVHT.ChannelBandwidth）；%从波形开始的数据包偏移量pktOffset = searchOffset + pktOffset;%如果未检测到数据包或偏移量超出波形范围，则停止如果is空(pktOffset) || (pktOffset<0) ||...（（pktoffset + Ind.lsig（2））> rxwaveFormLength）休息；结束%提取非ht字段并进行粗频偏移校正%允许可靠的符号定时nonht = rxwaveform（pktoffset +（ind.lstf（1）：Ind.lsig（2）），:);coarsefreqoff = wlancoarsecfoeStimate（nonht，cfgvht.channelbandwidth）;nonht = helperfrequencyOffset（非HU，FS，-CoarseFreqoff）;%确定L-LTF预期启动和实际启动之间的偏移量%的L-LTFlltfOffset=wlanSymbolTimingEstimate（非HT，cfgVHT.信道带宽）；%确定数据包偏移量pktOffset = pktOffset + lltfOffset;%如果偏移量没有波形边界，跳过采样并继续%在剩余波形内搜索如果(pktOffset<0) || ((pktOffset+pktLength)>rxWaveformLength) searchOffset = pktOffset+double(ind.LSTF(2))+1;持续；结束%定时同步完成；提取检测到的数据包rxpacket = rxwaveform（pktoffset +（1：pktlength），:);pktnum = pktnum + 1;DISP（[“包”num2str（pktnum）'在索引：'num2str（pktoffset + 1）]）;%对提取的数据包应用粗频率校正rxPacket = helperFrequencyOffset (fs, rxPacket -coarsefreqOff);%对提取的数据包进行精细的频率偏移校正lltf=rxPacket（ind.lltf（1）：ind.lltf（2），：）；%提取L-LTFFINEFREKOFF=WLANFINECFOFESTIMATE（lltf，cfgVHT.ChannelBandwidth）；rxPacket=helperFrequencyOffset（rxPacket，fs，-FINEFREKOFF）；提取VHT-LTF样本，解调并进行信道估计vhtltf = rxPacket (ind.VHTLTF (1): ind.VHTLTF (2):);vhtltfDemod = wlanVHTLTFDemodulate (vhtltf cfgVHT);获得单流信道估计chanEstSSPilots=vhtSingleStreamChannelEstimate（VHTLTFDEM，cfgVHT）；%信道估计chanEst=wlanVHTLTFChannelEstimate（VHTLTFDEMO，cfgVHT）；％执行IQ增益和相位不平衡估计(增益、相位、α,β,γ,dataRot] =...Helperiqimbalance（Chanest，CFGVHT）;fprintf（测量的IQ增益和相位不平衡:%2.2f dB， %2.2f deg\n、增益、相位);在信道上执行IQ增益和相位不平衡校正%的估计chanEst=chanEst./（alphaEst+betaEst.*伽马）；%如[5]中等式-29所述%光谱平坦度测量VHTTXSpectralflatnessmeasurement（Chanest，CFGVHT，PKTNUM，HSF）;％从波形中提取VHT数据样本vhtdata=rxPacket（ind.vhtdata（1）：ind.vhtdata（2），：）；%估计VHT数据场中的噪声功率noiseVarVHT = vhtNoiseEstimate (vhtdata chanEstSSPilots cfgVHT);％提取VHT数据样本并执行OFDM解调，均衡%相位跟踪[~, ~, eqSym] = wlanVHTDataRecover (vhtdata,成龙、noiseVarVHT cfgVHT,...“EqualizationMethod”，“ZF”，'先行者屏蔽'，“PreEQ”）;%使用迫零算法进行均衡%对VHT数据执行IQ增益和相位不平衡校正eqSym = eqSym。* dataRot;%数据子载波上的载波旋转eqSym =((连词(α)* eqSym) -(β*连词(eqSym(结束:1:1 ,:,:))))/(( abs(α)^ 2)——(abs(β)^ 2));%如[5]中等式-30所述在数据包的所有空间流中计算RMS EVMrmsEVM（pktNum）=EVMPerPkt（eqSym）；fprintf('RMS EVM:%2.2f%%，%2.2fdB\n'rmsEVM (pktNum), 20 * log10 (rmsEVM (pktNum) / 100));％计算RMS EVM每个子载波和数据包的空间流evmPerSC = evmPerSC (eqSym);% Nst-by-1-by-Nss%绘制每个子载波和均衡星座的RMS EVMVHTTXEVM星座图（等式SYM、evmPerSC、cfgVHT、pktNum、hCon、hEVM）；％存储波形内每个分组的偏移量pktOffsetStore (pktNum) = pktOffset;%增加波形偏移量并搜索数据包的剩余波形searchOffset = pktOffset + pktLength + minPktLen;结束

如果pktNum>0 fprintf('用于%d个数据包的平均EVM:%2.2f%%，%2.2fdB\n'，...PKTNUM，平均值（RMSEVM（1：PKTNUM）），20 * LOG10（平均值（RMSEVM（1：PKTNUM））/ 100））;别的disp (“未检测到完整的数据包”）;结束

％索引用于访问时域数据包中的每个字段印第安纳州= wlanFieldIndices (cfgVHT,“过度采样物料”，OSF）;startidx = Ind.Vhtdata（1）;%数据的开始endIdx=ind.VHTData（2）；%数据结束idleNSamps = idleTime * fs * osf;％空闲时间样本PERPKTLENGHT=endIdx+idleNSamps；idx=0（endIdx startIdx+1，numPackets）；为i = 1：numpacketsTXWAVEFORM中数据包的％开始，占滤波器延迟pktOffset =(张)* perPktLength;TxWaveForm中的非HT数据的％指标idx（：，i）=pktOffset+（startIdx:endIdx）；结束gatedVHTData=tx波形（idx（：），：）；

如果pktNum>0 helperSpectralMaskTest（gatedVHTData、fs、osf）；结束