cfgVHT = wlanVHTConfig;创建包配置cfgVHT。ChannelBandwidth =“CBW80”；% 80 MHzcfgVHT。numtransmitantenna = 1;一个发射天线cfgVHT。NumSpaceTimeStreams = 1;%一个时空流cfgVHT。STBC = false;%没有STBC，所以是一个空间流cfgVHT。MCS = 8;调制:256-QAMcfgVHT。APEPLength = 3000;% A-MPDU eof前填充长度(以字节为单位)

Osf = 3;% 3x过采样因子numPackets = 20;生成20个包idleTime = 10e-6;包之间的空闲时间% 10微秒

%创建随机数据;PSDULength以字节为单位savedState = rng(0);%设置随机状态data = randi([0 1]，cfgVHT.PSDULength*8*numPackets,1);%生成多包波形tx波形= wlanWaveformGenerator(数据，cfgVHT，.．.“OversamplingFactor”osf,“NumPackets”numPackets,“IdleTime”, idleTime);fs = wlanSampleRate(cfgVHT);波形的基带采样率

modeliq失衡= true;%设置为true添加智商增益和相位不平衡如果modeliqbalance iqGaindB = 1;%智商增益不平衡的dB，指定范围为[-1 1]iqPhaseDeg = 1;% IQ相位不平衡程度，从范围[-2]指定iqGainLin = db2mag(iqGaindB);将增益从dB转换为线性值tx波形= real(tx波形)+ 1i*imag(tx波形)*iqGainLin*exp(1j*iqPhaseDeg*pi/180);为[5]的方程1所示结束

p饱和度= 25;功率放大器的饱和功率% dBmhpaBackoff = 13;%功率放大器回退，单位为dB创建和配置一个无记忆非线性模型放大器非线性;非线性。方法=“拉普模式”；非线性。平滑度= 3;% p参数非线性。线性增益= -hpaBackoff;非线性。OutputSaturationLevel = db2mag(pSaturation-30);将模型应用于每个发射天线tx波形=非线性(tx波形);

thNoise = com . thermalnoise (“NoiseMethod”，“噪声图”，“SampleRate”fs * osf,“NoiseFigure”6);为i = 1:cfgVHT。numtransmitantenna tx波形(:，i) = thNoise(tx波形(:，i));结束

stop = 40;阻带衰减%ofdmInfo = wlanVHTOFDMInfo(“VHT-Data”, cfgVHT);OFDM参数SCS = fs/ofdmInfo.FFTLength;子载波间距%txbw = max(abs(ofdmInfo.ActiveFrequencyIndices))*2*SCS;已占用带宽百分比[L,M] =鼠(osf);maxLM = max([L M]);R = (fs-txbw)/fs;TW = 2*R/maxLM;%过渡宽度b = designMultirateFIR(L,M,TW, stopp);

firinterp = dsp.FIRRateConverter(M,L,b);rx波形= firinterp(tx波形);

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

获取访问时域数据包中每个字段的索引ind = wlanFieldIndices(cfgVHT);rxWaveformLength = size(rxwavelength,1);pktLength = double(ind.VHTData(2));定义我们可以检测到的最小数据长度;L-STF的长度%样本minPktLen = double(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;%如果没有检测到数据包或在波形边界外偏移，则停止如果isempty(pktOffset) || (pktOffset<0) ||.．.((pktOffset + ind.LSIG (2)) > rxWaveformLength)打破；结束提取非ht字段并进行粗频偏校正允许可靠的符号计时nonht = rx波形(pktOffset+(ind.LSTF(1):ind.LSIG(2))，:);coarfreqoff = wlanCoarseCFOEstimate(nonht,cfgVHT.ChannelBandwidth);nonht = helperFrequencyOffset(nonht,fs，-粗糙freqoff);确定L-LTF的预期开始和实际开始之间的偏移量L-LTF的%lltfOffset = wlanSymbolTimingEstimate(nonht,cfgVHT.ChannelBandwidth);确定数据包偏移量pktOffset = pktOffset+lltfOffset;%如果偏移量没有波形边界，则跳过采样并继续%在剩余波形中搜索如果(pktOffset<0) || ((pktOffset+pktLength)>rxWaveformLength) searchOffset = pktOffset+double(ind.LSTF(2))+1;继续；结束%定时同步完成;提取检测到的数据包rxPacket = rx波形(pktOffset+(1:pktLength)，:);pktNum = pktNum+1;disp ([“包”num2str (pktNum)' at index: 'num2str (pktOffset + 1)]);对提取的数据包进行粗频率校正rxPacket = helperFrequencyOffset(rxPacket,fs，-粗糙freqoff);对提取的数据包执行精细的频偏校正lltf = rxPacket(ind.LLTF(1):ind.LLTF(2)，:);提取L-LTFfineFreqOff = wlanFineCFOEstimate(lltf,cfgVHT.ChannelBandwidth);rxPacket = helperFrequencyOffset(rxPacket,fs，-fineFreqOff);提取VHT-LTF样本，解调并进行信道估计vhtltf = rxPacket(ind.VHTLTF(1):ind.VHTLTF(2)，:);vhtltfDemod = wlanVHTLTFDemodulate(vhtltf,cfgVHT);获得单个流通道估计chanEstSSPilots = vhtSingleStreamChannelEstimate(vhtltfDemod,cfgVHT);%渠道估计chanEst = whtltfchannelestimate (vhtltfDemod,cfgVHT);执行IQ增益和相位不平衡估计(增益、相位、α,β,γ,dataRot] =.．.helperIQImbalanceEstimate(成龙、cfgVHT);流(测量IQ增益和相位不平衡:%2.2f dB， %2.2f deg\n、增益、相位);在信道上执行IQ增益和相位不平衡校正%的估计chanEst = chanEst。/(alphaEst + betaEst.*gamma);为[5]的公式29所示谱平坦度测量vhtTxSpectralFlatnessMeasurement(成龙、cfgVHT pktNum hSF);从波形中提取VHT数据样本vhtdata = rxPacket(ind.VHTData(1):ind.VHTData(2)，:);估计VHT数据场的噪声功率noisevvht = vhtNoiseEstimate(vhtdata,chanEstSSPilots,cfgVHT);提取VHT数据样本并进行OFDM解调、均衡%和相位跟踪[~，~，eqSym] = wlanVHTDataRecover(vhtdata,chanEst, noisevvht,cfgVHT，.．.“EqualizationMethod”，“ZF”，“PilotPhaseTracking”，“PreEQ”）;%采用零强制算法进行均衡对VHT数据进行IQ增益和相位不平衡校正eqSym = eqSym.*dataRot;%数据子载波上的载波旋转eqSym =((连词(α)* eqSym) -(β*连词(eqSym(结束:1:1 ,:,:))))/(( abs(α)^ 2)——(abs(β)^ 2));为[5]的公式30所示为包计算所有空间流的RMS EVMrmsEVM(pktNum) = EVMPerPkt(eqSym);流(' RMS EVM: %2.2f%%， %2.2fdB\n'rmsEVM (pktNum), 20 * log10 (rmsEVM (pktNum) / 100));计算每个子载波和包的空间流的RMS EVMevmPerSC = evmPerSC (eqSym);% Nst-by-1-by-Nss绘制每个子载波和均衡星座的RMS EVMvhtTxEVMConstellationPlots (eqSym、evmPerSC cfgVHT、pktNum hCon, hEVM);存储波形内每个包的偏移量pktOffsetStore(pktNum) = pktOffset;增量波形偏移量并搜索数据包的剩余波形searchOffset = pktOffset+pktLength+minPktLen;结束

如果pktNum > 0流(' %d数据包的平均EVM: %2.2f%%， %2.2fdB\n'，.．.pktNum,意味着(rmsEVM (1: pktNum)), 20 * log10(平均(rmsEVM (1: pktNum)) / 100));其他的disp (“未检测到完整的数据包”）;结束

用于访问时域数据包中每个字段的索引ind = wlanFieldIndices(cfgVHT，“OversamplingFactor”(osf);startIdx = ind.VHTData(1);%数据起始endIdx = ind.VHTData(2);数据结束idleNSamps = idleTime*fs*osf;%空闲时间样本perPktLength = endIdx+idleNSamps;idx = 0 (endIdx-startIdx+1,numPackets);为i = 1:numPackets% tx波形中数据包的开始，考虑过滤器延迟pktOffset = (i-1)*perPktLength;tx波形中非ht数据的%指数idx(:，i) = pktOffset+(startIdx:endIdx);结束gatedVHTData = tx波形(idx(:)，:);

如果pktNum > 0 helperSpectralMaskTest (fs, gatedVHTData osf);结束