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职位:以本科学历为基础%指定为[方位角;仰角]度的角度％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 mspsprm.chantype.='散射';％频道选项：'散射'，'mimo'prm.nfig = 8;％噪声系数（增加到恶化，5-10 dB）prm.nrays = 500;FBB分区的FRF的百分比数量

人口、难民和移民事务局。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;%终止尾比特数人口、难民和移民事务局。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）;%发射天线阵列定义%数组位置和角度prm.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”那“自定义”）;人口、难民和移民事务局。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 =阶段。异丙哒植物;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.NUMSTS.= numSTS;％保持相同的阵列配置[rxPreSig, chanDelay] = helperApplyMUChannel(人口、难民和移民事务局,preambleSig spLoss);％信道状态信息反馈HDP = CELL（PRM.NUMUSERS，1）;prm.numsts = numtx;%设置为numTx以估计所有链接为了uIdx = 1: prm.numUsers%前端放大器增益和热噪声rxPreAmp =分阶段。ReceiverPreamp (......'获得'spLoss (uIdx),......%表示路径损失“NoiseFigure”人口、难民和移民事务局。NFig,“ReferenceTemperature”，290，......“SampleRate”, 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.NullCarrierIndices prm.PilotCarrierIndices);信道估计从序言％numcarr，numtx，numrxHDP {UIDX} = HelperMimochanneleStimate（Rxofdm（：，1：numtx，:)，prm）;结尾

％计算发射侧的混合体重如果prm.numUsers = = 1%单用户OMP% Spread rays in [az;el]=[-180:180;-90:90] 3D空间，等距% tang = [-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，size（at，1），尺寸（at，2））））;为了carrIdx = 1:人口、难民和移民事务局。numCarriers AtExp(carrIdx，:，:) = At;所有子运营商相同结尾％正交匹配追求混合体重润扬悬索桥(Fbb) = omphybweights(黄芪丹参滴丸{1},numSTS, numSTS, AtExp);v = Fbb;％设置基带预编码器（FBB）%频响在平坦信道的子载波之间是相同的MFRF =润滑（平均值（FRF，1），[2 3 1]）;其他的多用户联合空间分割多路复用(Fbb mFrf] = helperJSDMTransmitWeights(黄芪丹参滴丸,人口、难民和移民事务局);多用户基带预编码将每个用户的CSI打包成一个矩阵(块对角线)steeringMatrix = 0 (prm.numCarriers总和(numSTSVec)和(numSTSVec));为了UIDX = 1：prm.numusers stsidx = sum（numstsvec（1：uidx-1））+（1：numstsvec（Uidx））;steeringmatrix（：，stsidx，stsidx）= fbb {uidx};% Nst-by-Nsts-by-Nsts结尾v = permute(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;%恢复数据传输

％卷积编码器Encoder = Comm.convolutionAlencoder（......“TrellisStructure”(133 171 165), poly2trellis(7日),......'terminationmethod'那“终止”）;txdatabits = cell（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 = complex（zeros（prm.numcarriers，prm.numdatasymbols，numsts））;为了symIdx = 1: prm.numDataSymbols为了carrIdx = 1:人口、难民和移民事务局。numCarriers Q = squeeze(v(carrIdx，:，:));normQ = Q * sqrt(numTx)/norm(Q，'fro'）;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）））））;%生成带有反馈权重的前导和数据的前置PreeAblesigd = HelpergenProdeamble（PRM，V）;txsigsts = [前期行程指数; txofdm];%射频波束形成:对数字信号施加频响%每个天线元件连接到每个数据流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),......%表示路径损失“NoiseFigure”人口、难民和移民事务局。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);从映射前导码的％信道估计高清= helperMIMOChannelEstimate (rxOFDM (: 1: numSTS,:),人口、难民和移民事务局);%分配均衡%索引为感兴趣的用户流[rxeq，csi] = helpermimoequalize（rxofdm（：，numsts + 1：结束，:)，hd（：，stsidx，:)）;%软解调rxSymbs = rxEq:) /√(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 subplot(prm.;numUsers马克斯(numSTSVec), (uIdx-1) *马克斯(numSTSVec) + i);情节(重塑(rxEq (:,:, i) /√(numTx), [], 1),'。'）;轴广场XLIM（GCA，[ -  Scaler Scaler]）;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状态