s = rng (67);为可重复性设置RNG状态

％多用户系统，每个用户单个/多个流prm.numUsers= 4;%用户数量prm.numSTSVec = [3 2 1 2];每个用户的独立数据流数量prm.numsts = sum（prm.numstsvec）;%一定是2的幂人口、难民和移民事务局。numTx = prm.numSTS * 8;% BS发射天线数(2的功率)prm.numrx = prm.numstsvec * 4;%每个用户接收天线数量(任意>= numSTSVec)%每个用户具有相同的调制prm.bitspersubcarier = 4;％2：QPSK，4：16QAM，6：64QAM，8：256QAMprm.numdatasymbols = 10;% OFDM数据符号数％MS的位置：在假定原点BS％角指定为[方位角;仰角]度％az范围[-180180]，El在范围内[-9090]，例如，[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人口、难民和移民事务局。ChanType ='散射';%通道选项:“散射”，“MIMO”人口、难民和移民事务局。NFig = 8;％噪声系数（增加到恶化，5-10 dB）prm.nrays = 500;% Frf, Fbb分区的射线数

人口、难民和移民事务局。FFTLength = 256;人口、难民和移民事务局。CyclicPrefixLength = 64;人口、难民和移民事务局。numCarriers = 234;人口、难民和移民事务局。NullCarrierIndices = [1:7 129 256-5:256]';防护和直流百分比人口、难民和移民事务局。PilotCarrierIndices = [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;%终止尾比特数prm.numFrmBits = numSTSVec。*（prm.numDataSymbols * prm.numCarriers *.．.prm.bitsPerSubCarrier * codeRate) -numTails;人口、难民和移民事务局。modMode = 2 ^ prm.bitsPerSubCarrier;%调制顺序频道滤波器延迟的％帐户numPadSym = 3;百分数的符号到零prm.numpadzeros = numpadsym *（prm.fftlength + prm.cyclicprefixlength）;

人口、难民和移民事务局。cLight = physconst (“光速”）;人口、难民和移民事务局。λ= prm.cLight / prm.fc;%获取发送和接收数组信息[Istxura，Expactortx，Isrxura，Expfactorrx] = HelperarrayInfo（PRM，True）;%发射天线阵列定义%数组位置和角度人口、难民和移民事务局。posTx = (0, 0, 0);% BS/发射阵列位置，[x;y;z]，米如果isTxURA均匀矩形阵列txarray =分阶段。PartitionedArray (.．.“数组”,分阶段。([expFactorTx numSTS], 0.5 * prm.lambda),.．.'子阵列选择'，（numsts，numtx），'SubarraySteering'那“自定义”）;其他的％均匀线性阵列txarray = phased.ula（numtx，“ElementSpacing”0.5 * prm.lambda.．.'元素'，序列。异丙啶anteNneElement（“BackBaffled”，错误的））;结尾prm.postxelem = getElementPosition（txArray）/prm.lambda;sploss = zeros（prm.numusers，1）;prm.posrx = zeros（3，prm.numusers）;为了uIdx = 1: prm.numUsers％接收阵列如果isRxURA (uIdx)均匀矩形阵列rxarray =分阶段。PartitionedArray (.．.“数组”，相算诅咒（[expfactorrx（Uidx）numstsvec（Uidx）]，.．.0.5 * prm.lambda),'子阵列选择'，那些（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 =阶段。异丙哒植物;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.NUMSTS.= numSTS;%保持通道相同的阵列配置[rxPreSig, chanDelay] = helperApplyMUChannel(人口、难民和移民事务局,preambleSig spLoss);％信道状态信息反馈黄芪丹参滴丸=细胞(prm.numUsers, 1);PRM.NUMSTS.= numTx;％设置为numTx估计所有链接为了uIdx = 1: prm.numUsers%前端放大器增益和热噪声rxPreAmp =分阶段。ReceiverPreamp (.．.“获得”，spLoss（uIdx），.．.%表示路径损失“噪声系数”人口、难民和移民事务局。NFig,'referenceTemperature', 290,.．.'采样率', prm.chanSRate);rxpresigamp = rxpreamp（rxpresig {uidx}）;使用的副载波的%比例功率rxPreSigAmp = rxPreSigAmp * (sqrt(prm. prm. prm))FFTLength -.．.长度(prm.NullCarrierIndices)) / prm.FFTLength);% OFDM解调rxOFDM = ofdmdemod (rxPreSigAmp (chanDelay (uIdx) + 1:.．.结束(prm.numPadZeros-chanDelay (uIdx)):), prm.FFTLength,.．.prm.CyclicPrefixLength prm.CyclicPrefixLength,.．.PRM.NULLCARRIDIDICES，PRM.PILOTCARRIDIDICES）;信道估计从序言% numCarr, numTx, numRx黄芪丹参滴丸{uIdx} = helperMIMOChannelEstimate (rxOFDM (:, 1: numTx,:),人口、难民和移民事务局);结尾

％计算发射侧的混合体重如果prm.numUsers = = 1%单用户OMP% Spread rays in [az;el]=[-180:180;-90:90] 3D空间，等距％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:人口、难民和移民事务局。numCarriers AtExp(carrIdx，:，:) = At;所有子运营商相同结尾％正交匹配追求混合体重润扬悬索桥(Fbb) = omphybweights(黄芪丹参滴丸{1},numSTS, numSTS, AtExp);v = Fbb;%设置基带预编码器(Fbb)％FRF是跨越子载波平坦信道相同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-逐NSTS逐NSTS结尾v = permute(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'(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:人口、难民和移民事务局。numCarriers Q = squeeze(v(carrIdx，:，:));normQ = Q * sqrt(numTx)/norm(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.nullcarrierindes）））））;％生成与反馈的权重，并在前面加上数据序言preambleSigD = helperGenPreamble(人口、难民和移民事务局,v);txSigSTS = [preambleSigD; txOFDM];％RF波束形成：应用FRF到数字信号%每个天线元件连接到每个数据流txsig = txsigsts * mfrf;

%对发送信号应用一个空间定义的通道[rxSig, chanDelay] = helperApplyMUChannel(人口、难民和移民事务局,txSig spLoss, preambleSig);

hfig = figure（'姓名'那“每个流均衡码元星座”）;scFact = ((prm.FFTLength-length (prm.NullCarrierIndices)).．./ prm.FFTLength ^ 2) / numTx;据nVar = noisepow (prm.chanSRate prm.NFig 290) / scFact;解码器= comm.ViterbiDecoder (“InputFormat”那“Unquantized”那.．.'TrellisStructure'，poly2trellis（7，[133 171 165]），.．.“TerminationMethod”那“终止”那“OutputDataType”那“双”）;为了uIdx = 1: prm.numUsersstsU = numSTSVec(uIdx); stsIdx = sum(numSTSVec(1:(uIdx-1)))+(1:stsU);%前端放大器增益和热噪声rxPreAmp =分阶段。ReceiverPreamp (.．.“获得”，spLoss（uIdx），.．.%表示路径损失“噪声系数”人口、难民和移民事务局。NFig,'referenceTemperature', 290,.．.'采样率', 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.NULLCARRIDIDICES，PRM.PILOTCARRIDIDICES）;%从映射的序言进行信道估计高清= helperMIMOChannelEstimate (rxOFDM (: 1: numSTS,:),人口、难民和移民事务局);%分配均衡%索引为感兴趣的用户流[rxeq，csi] = helpermimoequalize（rxofdm（：，numsts + 1：结束，:)，hd（：，stsidx，:)）;%软解调rxSymbs = RXEQ（：）/ SQRT（numTx）;rxLLRBits = qamdemod（rxSymbs，prm.modMode，“UnitAveragePower”,真的,.．.'OutputType'那'近似'那'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）/ SQRT（numTx），[]，1），'。'）;轴广场甘氨胆酸xlim(、(标量标量));甘氨胆酸ylim(、(标量标量));标题([“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）]）;％计算并显示的误码率数量= comm.ErrorRate;措施=误码率(txDataBits {uIdx}, rxBits);流(' BER = %.5f;不。= %d;不。错误= %d\n'那.．.措施（1），测量（3），测量（2））;结尾

RNG;恢复RNG状态