S = RNG（67）;设置RNG状态为可重复性

每个用户单个/多个流％的多用户系统prm.numUsers= 4;用户的数量％prm.numSTSVec= [3 2 1 2];每个用户独立数据流的百分比prm.numSTS =总和（prm.numSTSVec）;%一定是2的幂prm.numTx = prm.numSTS * 8;BS的％数发射天线（2的幂）prm.numRx = prm.numSTSVec * 4;每个用户接收天线数量%(任意>= numSTSVec)每个用户有相同的调制prm.bitsPerSubCarrier = 4;% 2: QPSK, 4: 16QAM, 6: 64QAM, 8: 256QAMprm.numDataSymbols = 10;OFDM数据符号的数目% MS位置:假设BS在原点指定为[方位角;仰角]度的角度在范围[-180 180]，EL在范围[-90 90]，例如％AZ[45; 0]maxRange = 1000;离BS 1000米内的所有MSs人口、难民和移民事务局。mobileRanges = randi([1 maxRange]，1,prm.numUsers);人口、难民和移民事务局。mobileAngles =[兰德(prm.numUsers) * 360 - 180;...兰特(prm.numUsers) * 180 - 90);人口、难民和移民事务局。fc = 28 e9;％28 GHz的系统prm.chanSRate = 100e6;％频道采样率，100个Msps的prm.ChanType='散射';％的通道选项： '散射'， 'MIMO'prm.NFig = 8;％噪声系数（增加恶化，5-10分贝）prm.nRays = 500;光线为FRF，FBB分区的％数

prm.FFTLength = 256;prm.CyclicPrefixLength = 64;prm.numCarriers = 234;prm.NullCarrierIndices = [1：7 129 256-5：256]';％卫队和DC人口、难民和移民事务局。驾驶索引= [26 54 90 118 140 168 204 232];nonDataIdx = [prm.NullCarrierIndices;prm.PilotCarrierIndices];人口、难民和移民事务局。载波位置= setdiff((1:prm.FFTLength)'，排序(nonDataIdx);numSTS = prm.numSTS;numTx = prm.numTx;numRx = prm.numRx;numSTSVec = prm.numSTSVec; codeRate = 1/3;每个用户的代码率相同numTails = 6;终止尾比特的数量％人口、难民和移民事务局。numFrmBits = numSTSVec。* (prm.numDataSymbols * prm.numCarriers *...prm.bitsPerSubCarrier * codeRate) -numTails;人口、难民和移民事务局。modMode = 2 ^ prm.bitsPerSubCarrier;%调制顺序对信道滤波器延迟％帐户numPadSym = 3;符号ZEROPAD的％数prm.numPadZeros = numPadSym *（prm.FFTLength + prm.CyclicPrefixLength）;

prm.cLight = physconst（“光速”);prm.lambda = prm.cLight / prm.fc;％获取发射和接收阵列信息[isTxURA，expFactorTx，isRxURA，expFactorRx] = helperArrayInfo（PRM，TRUE）;％发送天线阵列定义％阵列位置和角度prm.posTx = [0; 0; 0];％BS /发射阵列的位置，[X; Y; ​​Z]，米如果isTxURA％均匀的矩形阵列txarray = phased.PartitionedArray (...“数组”,分阶段。([expFactorTx numSTS], 0.5 * prm.lambda),...'SubarraySelection'，酮（numSTS，numTx），“SubarraySteering”，“自定义”);其他的％均匀线性阵列txarray = phased.ULA（numTx，'ElementSpacing'，0.5 * prm.lambda，...'元件'，phased.IsotropicAntennaElement（“BackBaffled”、假));结束prm.posTxElem = getElementPosition（txarray）/prm.lambda;spLoss =零（prm.numUsers，1）;prm.posRx =零（3，prm.numUsers）;对于uIdx = 1: prm.numUsers％接收阵列如果isRxURA (uIdx)％均匀的矩形阵列rxarray = phased.PartitionedArray (...“数组”，phased.URA（[expFactorRx（uIdx）numSTSVec（uIdx）]...0.5 * prm.lambda),'SubarraySelection'，酮（numSTSVec（uIdx），...numRx（uIdx）），“SubarraySteering”，“自定义”);人口、难民和移民事务局。posRxElem = getElementPosition (rxarray) / prm.lambda;其他的如果numRx（uIdx）> 1％均匀线性阵列rxarray = phased.ULA（numRx（uIdx），...'ElementSpacing'，0.5 * prm.lambda，...'元件',phased.IsotropicAntennaElement);人口、难民和移民事务局。posRxElem = getElementPosition (rxarray) / prm.lambda;其他的rxarray = phased.IsotropicAntennaElement;prm.posRxElem = [0;0;0];% LCS结束结束%移动位置[XRX，YRX，ZRX] = sph2cart（deg2rad（prm.mobileAngles（1，uIdx）），...uIdx函数(prm.mobileAngles (2)),...prm.mobileRanges（uIdx））;prm.posRx（：，uIdx）= [XRX; YRX; ZRX];[toRxRange，toRxAng] = rangeangle（prm.posTx，prm.posRx（：，uIdx））;spLoss（uIdx）= fspl（toRxRange，prm.lambda）;结束

%产生前导信号prm.numSTS = numTx;％设置为numTx试探所有通道preambleSig = helperGenPreamble（PRM）;通过信道传送前导prm.numSTS = numSTS;％保留用于信道相同的阵列配置[rxPreSig, chanDelay] = helperApplyMUChannel(人口、难民和移民事务局,preambleSig spLoss);％信道状态信息反馈HDP =细胞（prm.numUsers，1）;prm.numSTS = numTx;%设置为numTx来估计所有链接对于uIdx = 1: prm.numUsers%前端放大器增益和热噪声rxPreAmp = phased.ReceiverPreamp（...'获得'spLoss (uIdx),...%考虑路径损耗“NoiseFigure”，prm.NFig，'ReferenceTemperature'，290，...“SampleRate”，prm.chanSRate）;rxPreSigAmp = rxPreAmp（rxPreSig {uIdx}）;使用过的副载波的比例功率rxPreSigAmp = rxPreSigAmp *（SQRT（prm.FFTLength  -...长度（prm.NullCarrierIndices））/ prm.FFTLength）;% OFDM解调rxOFDM = ofdmdemod (rxPreSigAmp (chanDelay (uIdx) + 1:...结束(prm.numPadZeros-chanDelay (uIdx)):), prm.FFTLength,...prm.CyclicPrefixLength，prm.CyclicPrefixLength，...prm.NullCarrierIndices，prm.PilotCarrierIndices）;从前导％的信道估计％numCarr，numTx，numRxHDP {uIdx} = helperMIMOChannelEstimate（rxOFDM（：，1：numTx，:)，PRM）;结束

％计算在发送侧的混合权重如果prm.numUsers = = 1%单用户OMP％传播光线在[AZ; EL] = [ -  180：180; -90：90]三维空间中，相等的间隔% txang = [-180:360/prm.nRays:180;90:180 / prm.nRays: 90);txang =[兰德(prm.nRays) * 360 - 180;兰德(prm.nRays) * 180 - 90);%随机在= steervec（prm.posTxElem，txang）;AtExp =络合物（零（prm.numCarriers，尺寸（在，1），大小（在，2）））;对于carrIdx = 1:人口、难民和移民事务局。numcarrier AtExp(carrIdx，:，:) = At;％同样为所有子载波结束％正交匹配追踪混合的权重润扬悬索桥(Fbb) = omphybweights(黄芪丹参滴丸{1},numSTS, numSTS, AtExp);v = Fbb;％设置基带预编码（FBB）对于平坦信道，在子载波之间的% Frf是相同的MFRF =置换（平均值（FRF，1），[2 3 1]）;其他的%多用户联合空间分割多路复用(Fbb mFrf] = helperJSDMTransmitWeights(黄芪丹参滴丸,人口、难民和移民事务局);%多用户基带预编码将每个用户CSI打包成一个矩阵(块对角)steeringMatrix = 0 (prm.numCarriers总和(numSTSVec)和(numSTSVec));对于uIdx = 1：prm.numUsers stsIdx =总和（numSTSVec（1：uIdx-1））+（1：numSTSVec（uIdx））;steeringMatrix（：，stsIdx，stsIdx）= {FBB uIdx};% Nst-by-Nsts-by-Nsts结束V =置换（steeringMatrix，[1 3 2]）;结束传输阵列模式图如果isTxURA为第一子载波％URA元件响应图案（txarray，prm.fc，-180：180，-90：90，“类型”，“efield”，...'ElementWeights',mFrf。*挤压(v (1::)),...'PropagationSpeed'，prm.cLight）;其他的％ULA第一副载波的阵列响应WTS = MFRF'*挤压（V（1，：，:)）;图案（txarray，prm.fc，-180：180，-90：90，“类型”，“efield”，...“权重”出世(:1),'PropagationSpeed'，prm.cLight）;结束prm.numSTS = numSTS;用于数据传输％恢复

％卷积编码器编码器= comm.ConvolutionalEncoder（...“TrellisStructure”(133 171 165),poly2trellis(7日),...'TerminationMethod'，“终止”);txDataBits =细胞（prm.numUsers，1）;的GridData =络合物（零（prm.numCarriers，prm.numDataSymbols，numSTS））;对于uIdx = 1: prm.numUsers%从每个用户的位生成映射符号txDataBits {uIdx} =兰迪([0,1],prm.numFrmBits (uIdx), 1);encodedBits =编码器(txDataBits {uIdx});%位到QAM符号映射mappedSym = qammod (encodedBits prm.modMode,“InputType”，“一点”，...'UnitAveragePower'，真正）;%映射到层:每个用户，每个符号，每个数据流stsIdx =总和（numSTSVec（1：（uIdx-1）））+（1：numSTSVec（uIdx））;的GridData（：，：，stsIdx）=重塑（mappedSym，prm.numCarriers，...prm.numDataSymbols，numSTSVec（uIdx））;结束假设反馈完美，对子载波应用预编码权preData =络合物（零（prm.numCarriers，prm.numDataSymbols，numSTS））;对于symIdx = 1: prm.numDataSymbols对于carrIdx = 1：prm.numCarriers Q =挤压（V（carrIdx，：，:)）;normQ = Q * SQRT（numTx）/常态（Q，“回回”);preData（carrIdx，symIdx，:) =挤压（的GridData（carrIdx，symIdx，:)）“。...* normQ;结束结束％的多天线的导频飞行员= helperGenPilots（prm.numDataSymbols，numSTS）;数据的％OFDM调制txOFDM = ofdmmod（preData，prm.FFTLength，prm.CyclicPrefixLength，...prm.NullCarrierIndices，prm.PilotCarrierIndices，飞行员）;使用过的副载波的比例功率txOFDM = txOFDM *（prm.FFTLength /...SQRT（（prm.FFTLength长度（prm.NullCarrierIndices））））;生成带有反馈权重的序言，并在数据前加上前缀preambleSigD = helperGenPreamble(人口、难民和移民事务局,v);txSigSTS = [preambleSigD; txOFDM];%射频波束形成:应用频响对数字信号每个天线元件连接到每个数据流txSig = txSigSTS * MFRF;

将空间定义的信道应用于发射信号[rxSig, chanDelay] = helperApplyMUChannel(人口、难民和移民事务局,txSig spLoss, preambleSig);

hfig =系数（'名称'，“均衡符号星座/流”);scFact =（（prm.FFTLength长度（prm.NullCarrierIndices））.../ prm.FFTLength ^ 2) / numTx;据nVar = noisepow (prm.chanSRate prm.NFig 290) / scFact;解码器= comm.ViterbiDecoder ('InputFormat'，“未量化”，...“TrellisStructure”，多边形格架(7，[133 171 165])，...'TerminationMethod'，“终止”，“OutputDataType”，“双”);对于uIdx = 1：prm.numUsers stsU = numSTSVec（uIdx）;stsIdx =总和（numSTSVec（1：（uIdx-1）））+（1：stsU）;%前端放大器增益和热噪声rxPreAmp = phased.ReceiverPreamp（...'获得'spLoss (uIdx),...%考虑路径损耗“NoiseFigure”，prm.NFig，'ReferenceTemperature'，290，...“SampleRate”，prm.chanSRate）;rxSigAmp = rxPreAmp（rxSig {uIdx}）;％比例功率占用的子载波rxSigAmp = rxSigAmp * (√prm.FFTLength-length (prm.NullCarrierIndices)).../prm.FFTLength）;% OFDM解调rxOFDM = ofdmdemod (rxSigAmp (chanDelay (uIdx) + 1:...结束(prm.numPadZeros-chanDelay (uIdx)):), prm.FFTLength,...prm.CyclicPrefixLength，prm.CyclicPrefixLength，...prm.NullCarrierIndices，prm.PilotCarrierIndices）;从映射前导％的信道估计HD = helperMIMOChannelEstimate（rxOFDM（：，1：numSTS，:)，PRM）;%分配均衡为感兴趣的用户将索引转换为流[RXEQ，CSI] = helperMIMOEqualize（rxOFDM（：，numSTS + 1：端，:)，HD（：，stsIdx，:)）;％软解调rxSymbs = rxEq:) /√(numTx);rxLLRBits = qamdemod (rxSymbs prm.modMode,'UnitAveragePower'，真正，...“输出类型”，'approxllr'，“NoiseVariance”，NVAR）;在解码之前应用CSIrxLLRtmp =重塑(rxLLRBits prm.bitsPerSubCarrier, [],...prm.numDataSymbols，stsU）;csitmp =重塑（CSI，1，[]，1，numSTSVec（uIdx））;rxScaledLLR = rxLLRtmp * csitmp。％软输入信道解码rxDecoded =解码器(rxScaledLLR (:));％解码接收到的比特rxBits = rxDecoded（1：prm.numFrmBits（uIdx））;％绘制均衡码为每个用户所有流定标器=装天花板(max (abs([真实(rxSymbs (:));图像放大(rxSymbs (:))))));对于i = 1:stsU subplot(prm)。numUsers马克斯(numSTSVec), (uIdx-1) *马克斯(numSTSVec) + i);情节(重塑(rxEq (:,:, i) /√(numTx), [], 1),'');轴广场XLIM（GCA，[ - 定标器定标器]）;ylim（GCA，[ - 定标器定标器]）;标题（[“U”num2str (uIdx)”,DS”num2str（I）]）;格上;结束％计算并显示EVMEVM = comm.EVM（“归一化”，“平均星座力量”，...“ReferenceSignalSource”，“参考星座估算”，...“ReferenceConstellation”，...qammod ((0: prm.modMode-1)、prm.modMode'UnitAveragePower'1));rmsEVM =维生素(rxSymbs);DISP（[“用户”num2str（uIdx）]）;DISP（['RMS EVM（％）='num2str（rmsEVM）]）;%计算和显示误码率BER = comm.ErrorRate;措施= BER（txDataBits {uIdx}，rxBits）;fprintf中（'BER =％.5f;位的数目=％d;错误的数目=％d \ n”，...(1),(3)措施的人,措施(2));结束

RNG（一个或多个）;％恢复RNG状态