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）人口、难民和移民事务局。numRx = prm.numSTSVec * 4;每个用户的接收天线数量（任何> = numstsvec）％每个用户具有相同的调制人口、难民和移民事务局。bitsPerSubCarrier = 4;% 2: qpsk, 4: 16qam, 6: 64qam, 8: 256qamprm.numDataSymbols= 10;％数量的OFDM数据符号% MS职位:以本科学历为基础%指定为[方位角;仰角]度的角度% az在范围[-180 180]，el在范围[-90 90]，例如[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系统人口、难民和移民事务局。chanSRate = 100 e6;%通道采样率，100 Mspsprm.chantype =.“散射”；％频道选项：'散射'，'mimo'prm.nfig = 8;%噪声值(增加至恶化，5-10 dB)prm.nRays= 500;FBB分区的FRF的百分比数量

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;零的符号数量人口、难民和移民事务局。numPadZeros = numPadSym * (prm.FFTLength + prm.CyclicPrefixLength);

prm.clight = physconst（'LightSpeed'）;Prm.lambda = prm.clight / prm.fc;％获取发送和接收数组信息[isTxURA, expFactorTx isRxURA expFactorRx] = helperArrayInfo(人口、难民和移民事务局,真实);％发送天线阵列定义％数组位置和角度prm.postx = [0; 0; 0];％bs /发射阵列位置，[x; y; z]，米如果isTxURA％均匀矩形阵列txarray = phased.partitionedArray（......'大批'，相位淘风器（[expactortx numsts]，0.5 * prm.lambda），......“SubarraySelection”的(numSTS numTx),“SubarraySteering”，'风俗'）;别的均匀线性阵列txarray =分阶段。齿龈(numTx'ElementsPacing'，0.5 * prm.lambda，......“元素”,分阶段。IsotropicAntennaElement ('底板'、假));结束人口、难民和移民事务局。posTxElem = getElementPosition (txarray) / prm.lambda;spLoss = 0 (prm.numUsers, 1);人口、难民和移民事务局。prm.numUsers posRx = 0(3日);为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;人口、难民和移民事务局。posRxElem = [0;0;0);％LCS.结束结束％手机职位(xRx yRx zRx] = sph2cart(函数(prm.mobileAngles (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）;%通道状态信息反馈HDP = CELL（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.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，numrxHDP {UIDX} = HelperMimochanneleStimate（Rxofdm（：，1：numtx，:)，prm）;结束

%计算发送侧的混合权值如果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 =复合物（零（prm.numcarriers，size（at，1），尺寸（at，2））））;为carrIdx = 1:人口、难民和移民事务局。numCarriers AtExp(carrIdx，:，:) = At;所有子运营商相同结束%正交匹配追求混合权重[FBB，FRF] = Omphybweights（HDP {1}，Numsts，Numsts，Atexch）;v = fbb;％设置基带预编码器（FBB）%频响在平坦信道的子载波之间是相同的MFRF =润滑（平均值（FRF，1），[2 3 1]）;别的％多用户联合空间分割多路复用[FBB，MFRF] = HelperjsdMtransmitweights（HDP，PRM）;多用户基带预编码将每个用户的CSI打包成一个矩阵(块对角线)steeringmatrix = zeros（prm.numcarriers，sum（numstsvec），sum（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 = erfute（steeringmatrix，[1 3 2]）;结束%发送阵列模式图如果isTxURA第一个子载波的％URA元素响应图案（TxArray，PRM.FC，-180：180，-90：90，“类型”，“efield”，......'Elementweights'，MFRF'*挤压（v（1，：，:)），......'繁殖'，prm.clight）;别的%齿龈针对第一子载波的％数组响应WTS = MFRF。'*挤压（v（1，：，:)）;图案（TxArray，PRM.FC，-180：180，-90：90，“类型”，“efield”，......“重量”出世(:1),'繁殖'，prm.clight）;结束prm.numsts = numsts;％恢复数据传输

%卷积编码器Encoder = Comm.convolutionAlencoder（......“TrellisStructure”，Poly2trellis（7，[133 171 165]），......'terminationmethod'，“终止”）;txdatabits = cell（prm.numusers，1）;griddata =复杂（零（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 = complex（zeros（prm.numcarriers，prm.numdatasymbols，numsts））;为Symidx = 1：prm.numdatasymbols为carridx = 1：prm.numcarriers q =挤压（v（carridx，：，:)）;organq = q * sqrt（numtx）/ norm（q，'fro'）;preData (carrIdx symIdx:) =挤压(gridData (carrIdx、symIdx:))。......* normQ;结束结束% Multi-antenna飞行员飞行员= HelperGenPilots（PRM.Numdatasymbols，Numsts）;ofdm调制数据txOFDM = ofdmmod (preData prm.FFTLength prm.CyclicPrefixLength,......PRM.NULLCARRIERINDICES，PRM.PILOTCARRIERINDICES，飞行员）;二手分载流子的％级别功率txofdm = txofdm *（prm.fftlength /......√(prm.FFTLength-length prm.NullCarrierIndices))));%生成带有反馈权重的前导和数据的前置PreeAblesigd = HelpergenProdeamble（PRM，V）;txsigsts = [前期行程指数; txofdm];%射频波束形成:对数字信号施加频响％每个天线元素都连接到每个数据流txSig = txSigSTS * mFrf;

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

hfig =图(“名字”，“每个流的均衡符号星座”）;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::,:),高清(:,stsIdx:));％软解调rxSymbs = rxEq:) /√(numTx);rxLLRBits = qamdemod (rxSymbs prm.modMode,'onemaveragepower'，真的，......“OutputType”，“approxllr”，“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（GCA，[ -  Scaler Scaler]）;ylim（GCA，[ - 缩放器缩放器]）;标题（['U'num2str (uIdx)”,DS”num2str（i）]）;网格在；结束%计算并显示EVM维生素与= comm.EVM (“归一化”，的平均星座力量”，......“ReferenceSignalSource”，'从参考星座估计'，......“ReferenceConstellation”，......qammod（（0：prm.modmode-1）'，prm.modmode，'onemaveragepower'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状态