comm金宝appSupportPackageCheck (“CST_WINNER2”）;

s = rng (10);％设置RNG种子可重复性％传输参数chanbw ='CBW80';％ 信道带宽麻木器= 3;％用户数量numstsvec = [3 1 4];每个用户的流数UserPos = [0 1 2];％用户职位mcsvec = [4 6 8];％MCS每用户：16QAM，64QAM，256QAMapepVec = [520 192 856];每个用户的有效载荷，以字节为单位chcodingvec = {'贝卡'那'ldpc'那'ldpc'};每个用户的频道编码%预编码和均衡参数precodingtype ='ZF';%预编码类型;ZF和MMSE的信噪比= 47;％snr在db中eqmethod ='ZF';％均衡方法%创建多用户VHT格式配置对象numtx = sum（numstsvec）;cfgvhtmu = wlanvhtconfig（'信道带宽'，chanbw，......'麻木商'numUsers,......'numtransmitantennas'numTx,......'groupid'2，......“NumSpaceTimeStreams”，numstsvec，......'userPositions'，userpos，......'MCS'，mcsvec，......“APEPLength”，apepvec，......'ChannelCoding'，chcodingvec）;

％VHT探测（NDP）配置，用于相同数量的流cfgVHTNDP = wlanVHTConfig ('信道带宽'，chanbw，......'麻木商'，1，......'numtransmitantennas'numTx,......'groupid'，0，......“NumSpaceTimeStreams”，sum（numstsvec），......'MCS'，0，......“APEPLength”, 0);

％生成NULL数据包，没有数据txNDPSig = wlanWaveformGenerator ([], cfgVHTNDP);NPDSigLen = size(txNDPSig, 1);

％为获胜者II通道设置布局参数aa = winner2.antennaArray（“台湾”、numTx 0.2);为了i = 1：numusers aa（i + 1）= winner2.antennaArray（'Ula'，numstsvec（i），0.05）;结尾Staidx = 2 :(麻木器+ 1）;apidx = {1};rndseed = 12;cfglayout = winner2.layoutparset（staidx，apidx，numusers，aa，[]，rndseed）;cfglayout.pairing = [（1，Numusers）; 2 :(麻木器+ 1）];对所有用户的一个接入点cfgLayout。numUsers ScenarioVector = 1 (1);所有链接的％A1场景cfglayout.propagconditionVector =零（1，Numusers）;%仿真结果为了我= 1：麻木器每个用户的％随机设置速度v = rand（3,1） -  0.5;cfglayout.stations（i + 1）.velocity = v / r $（v，'fro'）;结尾%设置WINNER II通道的模型参数cfgModel = winner2.wimparset;cfgModel。FixedPdpUsed ='是的';cfgmodel.fixedanglesused =.'是的';cfgmodel.intraclusterdsused ='不';cfgmodel.randomseed = 111;％可重复性％3用户的最大速度为1m / s。设置采样度％字段以确保采样率与通道带宽匹配。maxMSVelocity = max(cell2mat(cellfun(@(x)) norm(x，'fro'），......{cfglayout.stations.velocity}，“UniformOutput”，错误的）））;cfgmodel.uniformtimesampling ='是的';cfgmodel.sampleDenty = round（physconst（'LightSpeed'/......cfgmodel.centerfrequency / 2 /（maxmsvelocity / wlansamplervere（cfgvhtmu））））;％创建获胜者II通道系统对象winnerchan = comm.winner2channel（cfgmodel，cfglayout）;％调用info方法来检查一些派生信道参数chanInfo = info (WINNERChan)

chanDelay = chanInfo.ChannelFilterDelay;numPadZeros = 10;％设置modelconfig.numtimesample，以匹配输入信号的长度％避免警告winnerchan.modelconfig.numtimesamples = npdsiglen + numpadzeros;％声音所有用户的获胜者II频道chanoutndp = winnerchan（[txndpsig; zeros（numpadzeros，numtx）]）;％adgn.rxndpsig = cellfun（@ awgn，chanoutndp，......num2cell（SNR *那些（Numusers，1）），“UniformOutput”、假);

垫=细胞（麻木，1）;为了UIDX = 1：麻木器％基于每个用户的接收信号计算反馈矩阵垫{uIdx} = vhtCSIFeedback (rxNDPSig {uIdx} (chanDelay (uIdx) + 1:,:),......cfgVHTNDP、uIdx numSTSVec);结尾

％将每个用户CSI包装为矩阵numst = length（mat {1}）;%子载波数steeringmatrix = zeros（numst，sum（numstsvec），sum（numstsvec））;% Nst-by-Nt-by-Nsts为了uIdx = 1:numUsers stsIdx = sum(numSTSVec(1:uIdx-1))+(1:numSTSVec(uIdx));steeringMatrix (:,:, stsIdx) =垫{uIdx};% Nst-by-Nt-by-Nsts结尾％零强制或MMSE预编码解决方案如果比较字符串(precodingType'ZF'）delta = 0;％零强制别的delta =（numtx /（10 ^（snr / 10）））*眼睛（numtx）;%的患者结尾为了我= 1：numst％信道反演预编码h =挤压（steeringmatrix（i，：，:)）;steeringmatrix（我，：，:) = h /（h'* h + delta）;结尾%根据转向矩阵设置空间映射cfgvhtmu.spatialmapping =.“自定义”;cfgvhtmu.spatialmappingmatrix = perfute（steeringmatrix，[1 3 2]）;

%创建数据序列，每个用户一个Txdatabits = Cell（Numusers，1）;psdudatabits = cell（麻木器，1）;为了UIDX = 1：麻木器%为每个用户生成有效负载txdatabits {uidx} = randi（[0 1]，cfgvhtmu.apeplenth（uidx）* 8,1，'int8'）;用零的％pad有效载荷形成psdu(txDataBits psduDataBits {uIdx} = {uIdx};......零（（cfgvhtmu.psdulength（uidx）-cfgvhtmu.apeplength（uidx））* 8,1，'int8'));结尾

txSig = wlanWaveformGenerator (psduDataBits cfgVHTMU);

释放(WINNERChan);%设置ModelConfig。NumTimeSamples来匹配NDP +的总长度％波形和填充零WINNERChan.ModelConfig.NumTimeSamples =......winnerchan.modelconfig.numtimsamples +长度（txsig）+ numpadzeros;％通过获胜者II通道发送所有用户，10个全零％样本附加以考虑通道过滤器延迟chanOut = WINNERChan ([txNDPSig;0 (numPadZeros numTx);......TXSIG;零（numpadzeros，numtx）]）;％提取每个用户的波形输出chanout = cellfun（@（x）x（npdsiglen + numpadzeros + 1：结束，:)，chanout，“UniformOutput”、假);％adgn.rxSig = cellfun (@awgn chanOut,......num2cell（SNR *那些（Numusers，1）），“UniformOutput”、假);

％从配置中获取现场索引，在接收器时已知IND = WLANFIELDINDINDICES（CFGVHTMU）;%单用户接收者恢复有效负载位rxdatabits = cell（Numusers，1）;缩放器=零（麻木器，1）;spaxes = gobjects（总和（numstsvec），1）;hfig = figure（'名称'那'每流均等符号星座'）;为了UIDX = 1：NUMUSERS RXNSIG = RXSIG {UIDX}（Chandelay（UIDX）+1：结束，:);%用户时空流stsU = numSTSVec (uIdx);VHT字段中的估计噪声功率％lltf = rxNSig (ind.LLTF (1): ind.LLTF (2):);demodLLTF = wlanLLTFDemodulate (lltf chanBW);据nVar = helperNoiseEstimate (demodLLTF chanBW,总和(numSTSVec));%基于VHT-LTF进行信道估计rxvhtltf = rxnsig（ind.vhtltf（1）：Ind.vhtltf（2），:);Demodvhtltf = wlanvhtltfdemodulate（rxvhtltf，chanbw，numstsvec）;chanest = wlanvhtltfchannelestimate（demodvhtltf，chanbw，numstsvec）;％恢复VHT数据字段中的信息位rxvhtdata = rxnsig（ind.vhtdata（1）：Ind.vhtdata（2），:);[rxdatabits {uidx}，〜，eqsym] = wlanvhtdatarecover（rxvhtdata，......据nVar成龙,cfgVHTMU uIdx,......'均衡 - 方法'，eqmethod，“PilotPhaseTracking”那'没有任何'那......“LDPCDecodingMethod”那'分层BP'那'themplyldpciterationcount'6）;％绘制每个用户所有流的均等符号缩放器（UIDX）= CEIL（MAX（ABS（ABS（eqsym（:)）; impt（eqsym（:)））））））））;为了i = 1：stsu子图（NumUsers，Max（NumStsvec），（UIDX-1）* MAX（NumStsvec）+ i）;绘图（重塑（eqsym（：：：，i），[]，1），'。'）;轴正方形spax（sum（[0 numstsvec（1 :( uidx-1）]）+ i）= gca;％储存轴手柄标题（['用户'num2str（UIDX）”,流”num2str（i）]）;网格在;结尾结尾所有子尺度和比例图的％尺度轴为了i = 1：NUMER（SPAX）XLIM（SPAX（I），[ -  MAX（Scaler）Max（Scaler）]）;ylim（spaxes（i），[ -  max（scaler）max（scaler）]）;结尾pos = get（hfig，'位置'）;套装（HFIG，'位置'，[POS（1）* 0.7 POS（2）* 0.7 1.3 * POS（3）1.3 * POS（4）]）;

％将恢复位与每个用户Apepleencth信息位进行比较BER = INF（1，NUMUSERS）;为了uIdx = 1:numUsers idx = (1:cfgVHTMU.APEPLength(uIdx)*8).';[~,误码率(uIdx)] = biterr (txDataBits {uIdx} (idx) rxDataBits {uIdx} (idx));disp ([用户的误码率'num2str（UIDX）'：'num2str (ber (uIdx))));结尾RNG;％恢复rng状态