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位置:在原点假定为b指定为[方位角;仰角]度的%角度在范围[-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”人口、难民和移民事务局。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;终止尾比特的数量％人口、难民和移民事务局。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）;％发送天线阵列定义％阵列位置和角度人口、难民和移民事务局。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）），…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，PRM，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 * (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 = [兰特（1，prm.nRays）* 360-180;兰特（1，prm.nRays）* 180-90]。随机％在= steervec (prm.posTxElem txang);AtExp =复杂(0 (prm.numCarriers、尺寸(,1),大小(,2)));对于carrIdx = 1:人口、难民和移民事务局。= At;对所有副承运人相同结束％正交匹配追踪混合的权重[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.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”，poly2trellis（7，[133 171 165]），…“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，'输入类型'，“一点”，…'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 * (pr . fftlength /…SQRT（（prm.FFTLength长度（prm.NullCarrierIndices））））;%生成带有反馈权重的序文，并在数据前置preambleSigD = helperGenPreamble（PRM，V）;txSigSTS = [preambleSigD; txOFDM];% RF波束形成:对数字信号应用频响％每个天线元件被连接到每个数据流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”， [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 *（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:) /√(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) /√(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 = % 5;不。比特数= %d;不。错误数= %d\n'，…(1),(3)措施的人,措施(2));结束

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