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）人口、难民和移民事务局。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]最大范围=1000；%在BS 1000米范围内的所有MSsprm.mobileRanges=randi（[1最大范围]，1，prm.numUsers）；prm.mobileAngles=[rand（1，prm.numUsers）*360-180；......兰特（1，百万美元）*180-90]；prm.fc=28e9；%28GHz系统人口、难民和移民事务局。chanSRate = 100 e6;%通道采样率，100 Msps圆型=“散射”；％频道选项：'散射'，'mimo'prm.nfig = 8;%噪声值(增加至恶化，5-10 dB)prm.nRays= 500;FBB分区的FRF的百分比数量

prm.fft长度=256；prm.CyclicPrefixLength=64；prm.numCarriers=234；prm.NullCarrierIndices=[1:7129256-5:256]；%警卫和DCprm.PilotCarrierDices=[26 54 90 118 140 168 204 232]；Nondaidx=[prm.NullCarrierDices；prm.PilotCarrierDices]；prm.CarrierLocations=setdiff（（1:prm.FFTLength）’，sort（Nondaidx））；numSTS=prm.numSTS；numTx=prm.numTx；numRx=prm.numRx；numSTSVec=prm.numSTSVec；codeRate=1/3；%每个用户的代码速率相同数值=6；%终端尾位数人口、难民和移民事务局。numFrmBits = numSTSVec。* (prm.numDataSymbols * prm.numCarriers *......prm.bitsPerSubCarrier*编码率）-numails；prm.modMode=2^prm.bitsPerSubCarrier；%调制顺序%通道滤波器延迟numPadSym=3；%要归零的符号数人口、难民和移民事务局。numPadZeros = numPadSym * (prm.FFTLength + prm.CyclicPrefixLength);

prm.cLight=physconst(“光速”)；prm.lambda=prm.cLight/prm.fc；%获取发送和接收阵列信息[isTxURA, expFactorTx isRxURA expFactorRx] = helperArrayInfo(人口、难民和移民事务局,真实);%发射天线阵列定义%阵列位置和角度prm.postx = [0; 0; 0];％bs /发射阵列位置，[x; y; z]，米如果isTxURA％均匀矩形阵列txarray=相控分区阵列(......“数组”，分阶段.URA（[expFactorTx 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.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);人口、难民和移民事务局。posRxElem = getElementPosition (rxarray) / prm.lambda;其他的rxarray = phased.IsotropicAntennaElement;人口、难民和移民事务局。posRxElem = [0;0;0);%LCS结束结束%移动位置[xRx，yRx，zRx]=sph2cart（脱胶剂（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]=助手应用多通道（前置码、prm、spLoss）；%通道状态信息反馈HDP = CELL（PRM.NUMUSERS，1）;prm.numsts = numtx;%设置为numTx以估计所有链接为uIdx=1:prm.numUsers%前端放大器增益和热噪声RxStreamp=相控接收前置放大器(......'获得'spLoss (uIdx),......%路径损耗的原因“NoiseFigure”，prm.NFig，“参考温度”，290，......“SampleRate”, prm.chanSRate);rxPreSigAmp = rxPreAmp (rxPreSig {uIdx});%所用子载波的标度功率rxPreSigAmp = rxPreSigAmp * (sqrt(prm. prm. prm))FFTLength -......长度（prm.nullcarrierDices））/prm.FFTLength）；%正交频分复用解调rxOFDM=ofdmdemod（rxPreSigAmp（chanDelay（uIdx）+1:......完（努帕泽罗斯-钱德雷公墓（uIdx）），：），公墓长度，......prm.循环喷射长度，prm.循环喷射长度，......prm.无载货方、prm.飞行员载货方）；信道估计从序言％numcarr，numtx，numrxHDP {UIDX} = HelperMimochanneleStimate（Rxofdm（：，1：numtx，:)，prm）;结束

%计算发送侧的混合权值如果prm.numUsers==1%单用户OMP%在[az；el]=[-180:180；-90:90]三维空间中传播光线，等间距% tang = [-180:360/prm.nRays:180;90:180 / prm.nRays: 90);txang=[兰特（1，公升）*360-180；兰特（1，公升）*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，AtExp）；v=Fbb；％设置基带预编码器（FBB）%频响在平坦信道的子载波之间是相同的MFRF =润滑（平均值（FRF，1），[2 3 1]）;其他的%多用户联合空分复用[Fbb，mFrf]=助手JSDMTransmitweights（hDp，prm）；%多用户基带预编码将每个用户的CSI打包成一个矩阵(块对角线)steeringMatrix=零（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=排列（steeringMatrix[1 3 2]）；结束%发送阵列模式图如果isTxURA第一个子载波的％URA元素响应图案（TxArray，PRM.FC，-180：180，-90：90，“类型”，“efield”，......“元素权重”挤压（v（1，：，：），......“传播速度”，prm.cLight）；其他的%齿龈%第一副载波的阵列响应挤压（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，“输入类型”，“一点”，......“单位平均功率”，真的）;%映射到图层：每个用户、每个符号、每个数据流stsIdx=sum（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:prm.numCarriers Q=挤压（v（carrIdx，：，：）；normQ=Q*sqrt（numTx）/normQ（Q，'fro'）;preData (carrIdx symIdx:) =挤压(gridData (carrIdx、symIdx:))。......* normQ;结束结束% Multi-antenna飞行员pilots=helperGenPilots（prm.numDataSymbols，numSTS）；ofdm调制数据txOFDM = ofdmmod (preData prm.FFTLength prm.CyclicPrefixLength,......prm.空载体、prm.引航员载体、引航员）；%所用子载波的标度功率txofdm = txofdm *（prm.fftlength /......√(prm.FFTLength-length prm.NullCarrierIndices))));%生成带有反馈权重的前导和数据的前置PreeAblesigd = HelpergenProdeamble（PRM，V）;txsigsts = [前期行程指数; txofdm];%射频波束形成:对数字信号施加频响%每个天线单元连接到每个数据流txSig = txSigSTS * mFrf;

%将空间定义的信道应用于发射信号[rxSig，chanDelay]=帮助器应用多通道（txSig，prm，spLoss，前置码SIG）；

hfig =图(“名字”，“每个流的均衡符号星座”)；scFact=（（prm.fft长度（prm.nullcarrierndices））....../prm.fftlelength^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）；%前端放大器增益和热噪声RxStreamp=相控接收前置放大器(......'获得'spLoss (uIdx),......%路径损耗的原因“NoiseFigure”，prm.NFig，“参考温度”，290，......“SampleRate”，prm.chanSRate）；rxSigAmp=rxStreamp（rxSig{uIdx}）；%已占用子载波功率比例rxSigAmp=rxSigAmp*（sqrt（prm.FFT长度（prm.NullCarrierIndices））....../prm.fftlength）;%正交频分复用解调rxOFDM = ofdmdemod (rxSigAmp (chanDelay (uIdx) + 1:......完（努帕泽罗斯-钱德雷公墓（uIdx）），：），公墓长度，......prm.循环喷射长度，prm.循环喷射长度，......prm.无载货方、prm.飞行员载货方）；从映射前导码的％信道估计hD=帮助终端信道估计（rxOFDM（：，1:numSTS，：），prm）；%MIMO均衡%为感兴趣的用户索引到流中[rxEq, CSI] = helperMIMOEqualize (rxOFDM (numSTS + 1::,:),高清(:,stsIdx:));%软解调rxSymbs = rxEq:) /√(numTx);rxLLRBits = qamdemod (rxSymbs prm.modMode,“单位平均功率”，真的，......“OutputType”，“approxllr”，“NoiseVariance”，nvar）;%解码前应用CSIrxLLRtmp =重塑(rxLLRBits prm.bitsPerSubCarrier, [],......prm.numDataSymbols，stsU）；csitmp=重塑（CSI，1，[]，1，numSTSVec（uIdx））；rxScaledLLR=rxLLRtmp.*csitmp；％软输入通道解码rxdecode=解码器（rxScaledLLR（：）；％解码接收位rxBits = rxDecoded (1: prm.numFrmBits (uIdx));%为每个用户的所有流绘制均衡符号scaler=ceil（max（abs（[real（rxSymbs（：）；imag（rxSymbs（：）]））））；为i = 1:stsU subplot(prm.;numUsers马克斯(numSTSVec), (uIdx-1) *马克斯(numSTSVec) + i);情节(重塑(rxEq (:,:, i) /√(numTx), [], 1),“。”）;轴广场XLIM（GCA，[ -  Scaler Scaler]）;ylim（GCA，[ - 缩放器缩放器]）;标题（[“你”num2str (uIdx)”,DS”num2str（i）]；网格在；结束%计算并显示EVM维生素与= comm.EVM (“归一化”，的平均星座力量”，......“ReferenceSignalSource”，“根据参考星座估计”，......“ReferenceConstellation”，......qammod（（0:prm.modMode-1）’，prm.modMode，“单位平均功率”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状态