s=rng（21）；%为重复性设置RNG种子%传输参数chanBW=“CBW80”;%信道带宽numUsers=3；%活动用户数numSTSAll=[3 1 2]；%4个用户的流数userPos=[0123]；%最多4个用户的用户位置mcsVec=[4 6 2]；%最多4个用户的MCSapepVec=[15120819254006000]；%4个用户的有效负载（字节）chCodingVec={“密件抄送”,“自民党”,“自民党”,“密件抄送”};%4个用户的信道编码%信道和接收机参数香奈儿=“A型”;%TGac衰落信道模型预编码类型=“ZF”;%预编码类型；ZF或MMSE信噪比=38；%信噪比（dB）eqMethod=“ZF”;%均衡法%适当地创建多用户VHT格式配置对象%为活动用户建立向量值索引如果（numUsers==1）groupID=0；其他的groupID=2；终止numSTSVec=numSTSAll（1:numUsers）；numTx=sum（numSTSVec）；cfgVHTMU=wlanVHTConfig(“信道带宽”，chanBW，...“裸体者”，缪斯家，...“Numtransmitatenenas”，numTx，...“GroupID”，groupID，...“NumSpaceTimeStreams”，numSTSVec，...“用户位置”，userPos（1:numUsers），...“MCS”，mcsVec（1:numUsers），...“APEPLength”，apepVec（1:numUsers），...“信道编码”，chCodingVec（1:numUsers））；

%VHT探测（NDP）配置，用于相同数量的流cfgVHTNDP=wlanVHTConfig(“信道带宽”，chanBW，...“裸体者”1....“Numtransmitatenenas”，numTx，...“GroupID”, 0,...“NumSpaceTimeStreams”，总和（numSTSVec），...“MCS”, 0,...“APEPLength”, 0);

%生成无数据的空数据包txNDPSig=wlanWaveformGenerator（[]，cfgVHTNDP）；

%创建三个独立的频道TGAC=细胞（numUsers，1）；chanSeeds=[111122223344444]；%为最多4个用户选择uIndex=[105 2 1]；%为最多4个用户选择钱德雷=零（numUsers，1）；对于uIdx=1:numUsers TGAC{uIdx}=wlanTGacChannel(...“信道带宽”，cfgVHTMU.ChannelBandwidth，...“延迟配置文件”，chanMdl，...“用户索引”，uIndex（uIdx），...“Numtransmitatenenas”，numTx，...“numreceiveantens”，numSTSVec（uIdx），...“随机流”,“带种子的mt19937ar”,...“种子”，chanSeeds（uIdx），...“采样器”，wlanSampleRate（cfgVHTMU），...“传输接收距离”，5）；chanInfo=info（TGAC{uIdx}）；chanDelay（uIdx）=chanInfo.ChannelFilterDelay；终止

%附加零以允许信道滤波器延迟txNDPSig=[txNDPSig；零（10，numTx）]；%为所有传输流的每个用户设置独立通道rxNDPSig=单元（数值，1）；对于uIdx=1:numUsers rxNDPChan=TGAC{uIdx}（txNDPSig）；%每个接收器添加WGNrxNDPSig{uIdx}=awgn（rxNDPChan，snr）；终止

mat=单元（numUsers，1）；对于uIdx=1:numUsers%根据每个用户收到的信号计算反馈矩阵mat{uIdx}=vhtsif反馈（rxNDPSig{uIdx}（chanDelay（uIdx）+1:end，：），...cfgVHTNDP、uIdx、numSTSVec）；终止

%将每个用户的CSI打包到矩阵中numST=长度（mat{1}）；%子载波数steeringMatrix=0（numST，sum（numSTSVec），sum（numSTSVec））；%由新界北至新界北至新界北对于uIdx=1:numUsers stsIdx=sum（numSTSVec（1:uIdx-1））+（1:numSTSVec（uIdx））；steeringMatrix（：，：，stsIdx）=mat{uIdx}；%由新界北至新界北至新界北终止%迫零或MMSE预编码解决方案如果strcmp（预编码型，“ZF”)δ=0；%迫零其他的delta=（numTx/（10^（snr/10））*眼（numTx）；%彩信终止对于i=1:numST%信道反转预编码h=挤压（steeringMatrix（i，：，：）；steeringMatrix（i，：，：）=h/（h'*h+delta）；终止%基于转向矩阵设置空间映射空间映射=“习俗”; cfgVHTMU.SpatialMappingMatrix=permute（steeringMatrix[1 3 2]）；

%创建数据序列，每个用户一个txDataBits=单元（numuers，1）；psduDataBits=单元（numuers，1）；对于uIdx=1:numUsers%为每个用户生成有效负载txDataBits{uIdx}=randi（[01]，cfgVHTMU.APEPLength（uIdx）*8,1，“int8”);%用零填充有效负载以形成PSDUpsduDataBits{uIdx}=[txDataBits{uIdx}；...零（（cfgVHTMU.PSDULength（uIdx）-cfgVHTMU.APEPLength（uIdx））*8,1，“int8”)];终止

%生成多用户VHT波形txSig=WLANWAVEORMG发生器（psduDataBits，cfgVHTMU）；%通过每用户衰落信道传输rxSig=单元（numuers，1）；对于uIdx=1:numUsers%附加零以允许信道滤波器延迟rxSig{uIdx}=TGAC{uIdx}（[txSig；零（10，numTx）]）；终止

%从配置中获取字段索引，假设接收器已知ind=WLANFIELDICES（cfgVHTMU）；%单用户接收机恢复有效负载位rxDataBits=单元（numuers，1）；定标器=零（定标器，1）；spAxes=gobjects（总和（numSTSVec），1）；hfig=图(“姓名”,“每流均衡符号星座”);对于uIdx=1:numUsers%每个接收器添加WGNrxNSig=awgn（rxSig{uIdx}，snr）；rxNSig=rxNSig（chanDelay（uIdx）+1:end，：）；%用户时空流stsU=numSTSVec（uIdx）；%基于VHT-LTF执行信道估计rxVHTLTF=rxNSig（ind.VHTLTF（1）：ind.VHTLTF（2），：）；demovhtltf=wlanvhtltf解调（rxvhttf，chanBW，numSTSVec）；chanEst=wlanvhtltf信道估计（demovhtltf，chanBW，numSTSVec）；%获得单流信道估计chanEstSSPilots=vhtSingleStreamChannelEstimate（解调VHTLTF、cfgVHTMU）；%从波形中提取VHT数据样本rxVHTData=rxNSig（ind.VHTData（1）：ind.VHTData（2），：）；%估计VHT数据场中的噪声功率nVar=vhtNoiseEstimate（rxVHTData、chanEstSSPilots、cfgVHTMU）；%恢复VHT数据字段中的信息位[rxDataBits{uIdx}，~，eqsym]=wlanVHTDataRecover（rxVHTData，...chanEst、nVar、cfgVHTMU、uIdx、，“均衡法”，eqMethod，...“PilotPhaseTracking”,“没有”,“LDPCDecoding方法”,“分层bp”);%为每个用户的所有流绘制均衡符号scaler（uIdx）=ceil（max（abs（[real（eqsym（：））；imag（eqsym（：）]）））；对于i=1:stsU子图（numUsers，最大值（numSTSVec），（uIdx-1）*最大值（numSTSVec）+i）；绘图（重塑（等式SYM（：，：，i），[]，1），'.')；轴线广场SPAX（总和（[0 numSTSVec（1:（uIdx-1）））+i）=gca；%存储轴手柄头衔([“用户”num2str（uIdx）'溪流'num2str（i）]；网格在…上;终止终止%所有子地块和比例图形的比例轴对于i=1:numel（spAxes）xlim（spAxes（i），[-max（scaler）max（scaler）]；ylim（spAxes（i），[-max（scaler）max（scaler）]；终止pos=获取（hfig，“位置”)；设置（hfig，“位置”，[pos（1）*0.7 pos（2）*0.71.3*pos（3）1.3*pos（4）]；

%将恢复的位与每用户APEPLength信息位进行比较ber=inf（1，数值）；对于uIdx=1:numUsers idx=（1:cfgVHTMU.APEPLength（uIdx）*8）。；[~，ber（uIdx）]=biterr（txDataBits{uIdx}（idx），RXDDATABITS{uIdx}（idx））；disp(['用户的位错误率'num2str（uIdx）': 'num2str（ber（uIdx））]；终止rng（s）；%恢复RNG状态