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”人口、难民和移民事务局。NFig = 8;％噪声系数（增加恶化，5-10分贝）prm.nRays = 500;用于ff, 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];人口、难民和移民事务局。CarriersLocations = setdiff((1:prm.FFTLength)'， sort(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）;％发送天线阵列定义％阵列位置和角度人口、难民和移民事务局。posTx = (0, 0, 0);传送阵列位置，[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);％信道状态信息反馈黄芪丹参滴丸=细胞(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, numRx黄芪丹参滴丸{uIdx} = helperMIMOChannelEstimate (rxOFDM (:, 1: numTx,:),人口、难民和移民事务局);结束

％计算在发送侧的混合权重如果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 =复杂(0 (prm.numCarriers、尺寸(,1),大小(,2)));对于carrIdx = 1:人口、难民和移民事务局。AtExp(carrIdx，:，:) = At;所有子运营商相同结束％正交匹配追踪混合的权重润扬悬索桥(Fbb) = omphybweights(黄芪丹参滴丸{1},numSTS, numSTS, AtExp);v = Fbb;设置基带预编码器(Fbb)对于平坦信道，子载波间的频响相同mFrf = permute(mean(Frf,1)，[2 3 1]);其他的%多用户联合空间分多路复用(Fbb mFrf] = helperJSDMTransmitWeights(黄芪丹参滴丸,人口、难民和移民事务局);%多用户基带预编码%将每个用户CSI打包成矩阵(块对角线)steeringMatrix = 0 (prm.numCarriers总和(numSTSVec)和(numSTSVec));对于uIdx = 1: prm.numUsersstsIdx = sum(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 =细胞(人口、难民和移民事务局。numUsers, 1);gridData =复杂(0 (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 =复杂(0 (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”，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),…%考虑路径损失“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子图(prm)。numUsers马克斯(numSTSVec), (uIdx-1) *马克斯(numSTSVec) + i);情节(重塑(rxEq (:,:, i) /√(numTx), [], 1),'');轴广场甘氨胆酸xlim(、(标量标量));甘氨胆酸ylim(、(标量标量));标题([“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）]）;%计算和显示误码率数量= comm.ErrorRate;措施=误码率(txDataBits {uIdx}, rxBits);流(' BER = %.5f;不。的位= %d;不。误差= %d\n'，…(1),(3)措施的人,措施(2));结束

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