s = rng (67);％设置重复性的RNG状态

％多用户系统，每个用户单个/多个流prm.numusers = 4;％用户数量PRM.NUMSTSVEC.= [3 2 1 2];每个用户的独立数据流数量prm.numsts = sum（prm.numstsvec）;％必须是2的力量prm.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;所有MSS在BS的1000米范围内prm.mobileranges = randi（[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 mspsprm.chantype =.'散射';%通道选项:“散射”，“MIMO”人口、难民和移民事务局。NFig = 8;％噪声系数（增加到恶化，5-10 dB）prm.nrays = 500;% Frf, Fbb分区的射线数

prm.fftlength = 256;prm.cyclicprefixlength = 64;prm.numcarriers = 234;prm.nullcarrierindices = [1：7 129 256-5：256]'％守卫和直流prm.pilotcarrierindices = [26 54 90 118 140 168 204 232]';nondataidx = [prm.nullcarrierindices;prm.pilotcarrierindices];PRM.CarriersLocations = setDiff（（1：prm.fftlength）'，sort（nondataidx））;numsts = prm.numsts;numtx = prm.numtx;numrx = prm.numrx;numstsvec = prm.numstsvec;分数= 1/3;每个用户相同的代码率numtails = 6;％终止尾位数人口、难民和移民事务局。numFrmBits = numSTSVec。* (prm.numDataSymbols * prm.numCarriers *．..prm.bitspersubcarier * coderate） -  numtails;prm.modmode = 2 ^ prm.bitspersubcarier;％调制顺序频道滤波器延迟的％帐户numpadsym = 3;零的符号数量prm.numpadzeros = numpadsym *（prm.fftlength + prm.cyclicprefixlength）;

prm.clight = physconst（'LightSpeed'）;Prm.lambda = prm.clight / prm.fc;％获取发送和接收数组信息[Istxura，Expactortx，Isrxura，Expfactorrx] = HelperarrayInfo（PRM，True）;％发送天线阵列定义％数组位置和角度人口、难民和移民事务局。posTx = (0, 0, 0);% BS/发射阵列位置，[x;y;z]，米如果isTxURA均匀矩形阵列txarray = phased.partitionedArray（．..'大批'，相位淘风器（[expactortx numsts]，0.5 * prm.lambda），．..'子阵列选择'，（numsts，numtx），“SubarraySteering”那'风俗'）;别的％均匀线性阵列txarray = phased.ula（numtx，“ElementSpacing”0.5 * prm.lambda．..'元素'，序列。异丙啶anteNneElement（'底板'、假));结尾prm.postxelem = getElementPosition（txArray）/prm.lambda;sploss = zeros（prm.numusers，1）;prm.posrx = zeros（3，prm.numusers）;为了UIDX = 1：PRM.Numusers％接收阵列如果Isrxura（UIDX）均匀矩形阵列rxarray = phased.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）），．..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声音所有频道PreeAblesig = HelpergenProdeamble（PRM）;％在通道上传输序言prm.numsts = numsts;%保持通道相同的阵列配置[rxpresig，chandelay] = helperapplymuchannel（preamblesig，prm，sploss）;％信道状态信息反馈黄芪丹参滴丸=细胞(prm.numUsers, 1);prm.numsts = numtx;%设置为numTx以估计所有链接为了UIDX = 1：PRM.Numusers％前端放大器增益和热噪声rxpreamp = phased.receiverProamp（．..“获得”spLoss (uIdx),．..％账户用于路径损耗“NoiseFigure”，prm.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.NULLCARRIDIDICES，PRM.PILOTCARRIDIDICES）;信道估计从序言% numCarr, numTx, numRx黄芪丹参滴丸{uIdx} = helperMIMOChannelEstimate (rxOFDM (:, 1: numTx,:),人口、难民和移民事务局);结尾

％计算发射侧的混合体重如果prm.numusers == 1％单用户OMP[az; el] = [ -  180：180; -90：90] 3d空间，等间距% tang = [-180:360/prm.nRays:180;90:180 / prm.nRays: 90);txang = [rand（1，prm.nrays）* 360-180;兰特（1，prm.nrays）* 180-90];％ 随机的在= steervec (prm.posTxElem txang);AtExp =复杂(0 (prm.numCarriers、尺寸(,1),大小(,2)));为了carrIdx = 1:人口、难民和移民事务局。numCarriers AtExp(carrIdx，:，:) = At;所有子运营商相同结尾％正交匹配追求混合体重[FBB，FRF] = Omphybweights（HDP {1}，Numsts，Numsts，Atexch）;v = fbb;%设置基带预编码器(Fbb)%频响在平坦信道的子载波之间是相同的mFrf = permute(mean(Frf,1)，[2 3 1]);别的％多用户联合空间分割多路复用[FBB，MFRF] = HelperjsdMtransmitweights（HDP，PRM）;％多用户基带预编码将每个用户的CSI打包成一个矩阵(块对角线)steeringmatrix = zeros（prm.numcarriers，sum（numstsvec），sum（numstsvec））;为了UIDX = 1：PRM.NumusersstsIdx = sum(numSTSVec(1:uIdx-1))+(1:numSTSVec(uIdx)); steeringMatrix(:,stsIdx,stsIdx) = Fbb{uIdx};% Nst-by-Nsts-by-Nsts结尾v = erfute（steeringmatrix，[1 3 2]）;结尾%发送阵列模式图如果isTxURA% URA元素对第一副载波的响应模式(txarray、prm.fc 180:180 90:90,“类型”那“efield”那．..'Elementweights'，MFRF'*挤压（v（1，：，:)），．..'繁殖'，prm.clight）;别的％ula针对第一子载波的％数组响应WTS = MFRF。'*挤压（v（1，：，:)）;模式(txarray、prm.fc 180:180 90:90,“类型”那“efield”那．..'重量'出世(:1),'繁殖'，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} = randi（[0,1]，prm.numfrmbits（uidx），1）;EncodedBits = encoder（TxDatabits {UIDX}）;％位为QAM符号映射mappedsym = qammod（encodedbits，prm.modmode，'输入类型'那“一点”那．..'onemaveragepower',真正的);％映射到图层：每个用户，每个符号，每个数据流stsidx = sum（numstsvec（1 :( uidx-1）））+（1：numstsvec（Uidx））;griddata（：，：，stsidx）=重塑（映射yms，prm.numcarriers，．..prm.numDataSymbols numSTSVec (uIdx));结尾%对子载波应用预编码权，假设是完美反馈preData =复杂(0 (prm.numCarriers, prm.numDataSymbols numSTS));为了Symidx = 1：prm.numdatasymbols为了carridx = 1：prm.numcarriers q =挤压（v（carridx，：，:)）;organq = 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）））））;%生成带有反馈权重的前导和数据的前置PreeAblesigd = HelpergenProdeamble（PRM，V）;txsigsts = [前期行程指数; txofdm];%射频波束形成:对数字信号施加频响％每个天线元素都连接到每个数据流txsig = txsigsts * mfrf;

％将空间定义的通道应用于发射信号[rxsig，chandelay] = helperapplymuchannel（txsig，prm，sploss，preeamblesig）;

hfig = figure（'姓名'那“每个流的均衡符号星座”）;scfact =（（prm.fftlength-length-（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 = sum（numstsvec（1 :( uidx-1）））+（1：stsu）;％前端放大器增益和热噪声rxpreamp = phased.receiverProamp（．..“获得”spLoss (uIdx),．..％账户用于路径损耗“NoiseFigure”，prm.nfig，'referenceTemperature', 290,．..“SampleRate”，prm.chansrate）;rxsigamp = rxpreamp（rxsig {uidx}）;占用子运营商的％级别功率rxsigamp = rxsigamp *（sqrt（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）;%从映射的序言进行信道估计HD = Helpermimochannelestimate（Rxofdm（：，1：numsts，:)，prm）;％MIMO均衡为您感兴趣的用户的物流索引[rxeq，csi] = helpermimoequalize（rxofdm（：，numsts + 1：结束，:)，hd（：，stsidx，:)）;%软解调rxSymbs = rxEq:) /√(numTx);rxLLRBits = qamdemod (rxSymbs prm.modMode,'onemaveragepower',真的,．..'OutputType'那'近似'那“NoiseVariance”据nVar);％在解码之前应用CSIrxLLRtmp =重塑(rxLLRBits prm.bitsPerSubCarrier, [],．..Prm.numdatasymbols，Stsu）;csitmp = REPAPE（CSI，1，[]，1，NumStsvec（UIDX））;rxscaledllr = rxllrtmp。* csitmp;软输入信道译码rxdecoded =解码器（rxscaledllr（:)）;%解码接收位rxbits = rxdecoded（1：prm.numfrmbits（uidx））;％绘制每个用户所有流的均等符号缩放者= ceil（max（abs（abs（rxsymbs（rxsymbs（rxsymbs）; imag（rxsymbs（:)）））））;为了i = 1:stsU subplot(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，'onemaveragepower'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状态