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;％1000米以内BS的所有MSprm.mobileRanges =兰迪（[1 maxRange]，1，prm.numUsers）;prm.mobileAngles = [兰特（1，prm.numUsers）* 360-180;...兰特（1，prm.numUsers）* 180-90]。prm.fc = 28e9;％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]';％卫队和DCprm.PilotCarrierIndices = [26 54 90 118 140 168 204 232]';nonDataIdx = [prm.NullCarrierIndices;prm.PilotCarrierIndices];prm.CarriersLocations = setdiff（（1：prm.FFTLength）”，排序（nonDataIdx））;numSTS = prm.numSTS;numTx = prm.numTx;numRx = prm.numRx;numSTSVec = prm.numSTSVec;码率1/3 =;每个用户％相同的码速率numTails = 6;终止尾比特的数量％prm.numFrmBits = numSTSVec。*（prm.numDataSymbols * prm.numCarriers *...prm.bitsPerSubCarrier *码率）-numTails;prm.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（...“数组”，phased.URA（[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'，“自定义”）;prm.posRxElem = getElementPosition（rxarray）/prm.lambda;其他如果numRx（uIdx）> 1％均匀线性阵列rxarray = phased.ULA（numRx（uIdx），...'ElementSpacing'，0.5 * prm.lambda，...'元件'，phased.IsotropicAntennaElement）;prm.posRxElem = getElementPosition（rxarray）/prm.lambda;其他rxarray = phased.IsotropicAntennaElement;prm.posRxElem = [0;0;0];％LCS结束结束％移动位置[XRX，YRX，ZRX] = sph2cart（deg2rad（prm.mobileAngles（1，uIdx）），...deg2rad（prm.mobileAngles（2，uIdx）），...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，PRM，spLoss）;％信道状态信息反馈HDP =细胞（prm.numUsers，1）;prm.numSTS = numTx;％设置为numTx估计所有链接对于uIdx = 1：prm.numUsers％前端放大器的增益和热噪声rxPreAmp = phased.ReceiverPreamp（...'获得'，spLoss（uIdx），...％占路损“噪声系数”，prm.NFig，'ReferenceTemperature'，290，...'采样率'，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 = [兰特（1，prm.nRays）* 360-180;兰特（1，prm.nRays）* 180-90]。随机％在= steervec（prm.posTxElem，txang）;AtExp =络合物（零（prm.numCarriers，尺寸（在，1），大小（在，2）））;对于carrIdx = 1：prm.numCarriers AtExp（carrIdx，：，:) =在;％同样为所有子载波结束％正交匹配追踪混合的权重[FBB，FRF] = omphybweights（HDP {1}，numSTS，numSTS，AtExp）;V = FBB;％设置基带预编码（FBB）％FRF是跨越子载波平坦信道相同MFRF =置换（平均值（FRF，1），[2 3 1]）;其他％多用户联合空间分复用[FBB，MFRF] = helperJSDMTransmitWeights（HDP，PRM）;％的多用户预编码的基带％包中的每个用户的CSI成矩阵（方框对角线）steeringMatrix =零（prm.numCarriers，总和（numSTSVec），总和（numSTSVec））;对于uIdx = 1：prm.numUsers stsIdx =总和（numSTSVec（1：uIdx-1））+（1：numSTSVec（uIdx））;steeringMatrix（：，stsIdx，stsIdx）= {FBB uIdx};％NST-逐NSTS逐NSTS结束V =置换（steeringMatrix，[1 3 2]）;结束％发射阵列图案图如果isTxURA为第一子载波％URA元件响应图案（txarray，prm.fc，-180：180，-90：90，'类型'，“数字油田”，...'ElementWeights'，MFRF'*挤压（V（1，：，:)）...'PropagationSpeed'，prm.cLight）;其他％ULA为第一子载波％阵列响应WTS = MFRF'*挤压（V（1，：，:)）;图案（txarray，prm.fc，-180：180，-90：90，'类型'，“数字油田”，...“权重”中，WT（：，1），'PropagationSpeed'，prm.cLight）;结束prm.numSTS = numSTS;用于数据传输％恢复

％卷积编码器编码器= comm.ConvolutionalEncoder（...'TrellisStructure'，poly2trellis（7，[133 171 165]），...'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，'输入类型'，“位”，...'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（PRM，V）;txSigSTS = [preambleSigD; txOFDM];％RF波束形成：应用FRF到数字信号％每个天线元件被连接到每个数据流txSig = txSigSTS * MFRF;

％应用空间限定信道向发射信号[rxSig，chanDelay] = helperApplyMUChannel（txSig，PRM，spLoss，preambleSig）;

hfig =系数（'名称'，“每个流均衡码元星座”）;scFact =（（prm.FFTLength长度（prm.NullCarrierIndices））.../prm.FFTLength^2)/numTx;NVAR = noisepow（prm.chanSRate，prm.NFig，290）/ scFact;解码器= comm.ViterbiDecoder（'InputFormat'，“未量化”，...'TrellisStructure'，poly2trellis（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），...％占路损“噪声系数”，prm.NFig，'ReferenceTemperature'，290，...'采样率'，prm.chanSRate）;rxSigAmp = rxPreAmp（rxSig {uIdx}）;％比例功率占用的子载波rxSigAmp = rxSigAmp *（SQRT（prm.FFTLength长度（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）;％MIMO均衡％索引流感兴趣的用户[RXEQ，CSI] = helperMIMOEqualize（rxOFDM（：，numSTS + 1：端，:)，HD（：，stsIdx，:)）;％软解调rxSymbs = RXEQ（：）/ SQRT（numTx）;rxLLRBits = qamdemod（rxSymbs，prm.modMode，'UnitAveragePower'，真正，...“输出类型”，'approxllr'，'NoiseVariance'，NVAR）;％之前应用CSI到解码rxLLRtmp =重塑（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（:)）; IMAG（rxSymbs（:)）]）））;对于I = 1：stsU副区（prm.numUsers，最大值（numSTSVec），（uIdx-1）*最大（numSTSVec）+ I）;情节（重塑（RXEQ（：，：，I）/ SQRT（numTx），[]，1），''）;轴广场XLIM（GCA，[ - 定标器定标器]）;ylim（GCA，[ - 定标器定标器]）;标题（['U'num2str（uIdx）“S，D”num2str（I）]）;格上;结束％计算并显示EVMEVM = comm.EVM（'正常化'，“平均星座力量”，...'ReferenceSignalSource'，“从参考星座估计”，...'ReferenceConstellation'，...qammod（（0：prm.modMode-1）”，prm.modMode，'UnitAveragePower'，1））;rmsEVM = EVM（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状态