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位置：假定原点BS指定为[方位角;高度]度的％角度％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'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;％终止尾位数prm.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）;％发送天线阵列定义％数组位置和角度prm.postx = [0; 0; 0];％bs /发射阵列位置，[x; y; z]，米如果istxura.％均匀矩形阵列txarray = phased.partitionedArray（......'大批'，相位淘风器（[expactortx numsts]，0.5 * prm.lambda），......'子阵列选择'，（numsts，numtx），'子阵列'那'风俗'）;别的％均匀线性阵列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）），'子阵列'那'风俗'）;prm.posrxelem = getElementPosition（rxarray）/prm.lambda;别的如果numrx（Uidx）> 1％均匀线性阵列rxarray = phased.ula（numrx（uidx），......'ElementsPacing'，0.5 * prm.lambda，......'元素'，序列。异丙哒植物;prm.posrxelem = getElementPosition（rxarray）/prm.lambda;别的rxarray =阶段。异丙哒植物;prm.posrxelem = [0;0;0];％LCS.结尾结尾％手机职位[XRX，YRX，ZRX] = SPH2CART（DEG2RAD（PRM.MOBILEANGLES（1，UIDX）），......DEG2RAD（PRM.MOBILEANGLES（2，UIDX）），......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），......％账户用于路径损耗'噪音文件'，prm.nfig，'referenceTemperature'，290，......'采样率'，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空间，等间距％txang = [-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）;fbb =复杂（零（prm.numcarriers，numsts，numsts））;frf =复杂（零（prm.numcarriers，numsts，numtx））;为了carridx = 1：prm.numcarriers [fbb（carridx，：，:)，frf（carridx，：，:)] = Helperomptransmitweights（......erfut（hdp {1}（carridx，：：），[2 3]），numsts，numsts，at）;结尾v = fbb;％设置基带预编码器（FBB）平面通道的子载波中％FRF相同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};逐个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'那......'重量'，wts（：1），'繁殖'，prm.clight）;结尾prm.numsts = numsts;％恢复数据传输

％卷积编码器Encoder = Comm.convolutionAlencoder（......'trellisstructure'，Poly2trellis（7，[133 171 165]），......'terminationmethod'那'终止'）;％位为QAM符号映射modrqam = comm.rectangularqammodulator（......'调制顺序'，prm.modmode，'bitInpul'，真的，......'正常化方法'那'平均力量'）;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}）;mappedsym = modrqam（endodedbits）;％映射到图层：每个用户，每个符号，每个数据流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;结尾结尾％多天线飞行员飞行员= 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];％RF波束成形：将FRF应用于数字信号％每个天线元素都连接到每个数据流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;emodrqam = comm.rectangularqamdemodulator（......'调制顺序'，prm.modmode，'bitoup'，真的，......“决策方法”那'近似日志 - 似然比'那......'正常化方法'那'平均力量'那'方差'，nvar）;解码器= 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），......％账户用于路径损耗'噪音文件'，prm.nfig，'referenceTemperature'，290，......'采样率'，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（:) / sqrt（numtx）;rxllrbits = demodrqam（rxsymbs）;％在解码之前应用CSIrxllrtmp = Rehape（rxllrbits，prm.bitspersubcarier，[]，......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子平板（prm.numusers，max（numstsvec），（uidx-1）* max（numstsvec）+ i）;绘图（Rehape（RXEQ（：，：，i）/ sqrt（numtx），[]，1），'。'）;轴正方形XLIM（GCA，[ -  Scaler Scaler]）;ylim（GCA，[ - 缩放器缩放器]）;标题（['U'num2str（UIDX）'，DS'num2str（i）]）;网格在;结尾％计算并显示EVMEVM = comm.evm（'正常化'那'平均星座力量'那......'referenceignalsource'那'从参考星座估计'那......'ReferenceConstellation'，星座（Demodrqam））;RMSEVM = EVM（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状态