s=rng（67）；%为重复性设置RNG状态

%每个用户具有单个/多个流的多用户系统prm.numUsers=4；％用户数量prm.numSTSVec=[3 2 1 2]；%每个用户的独立数据流数prm.numSTS=总和（prm.numSTSVec）；%必须是2的幂prm.numtx = prm.numsts * 8;%基站发射天线数量（2个功率）prm.numRx=prm.numSTSVec*4；%每个用户的接收天线数（任意>=numSTSVec）%每个用户具有相同的调制prm.bitsPerSubCarrier=4；%2:QPSK，4:16QAM，6:64QAM，8:256QAMprm.numDataSymbols=10；%OFDM数据符号数%MS位置：假设原点为BS%指定为[方位角；仰角]度的角度%az在[-180]范围内，el在[-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；％28 GHz系统prm.chanSRate=100e6；%信道采样率，100毫秒/秒圆型=“分散”;%信道选项：“散射”、“MIMO”prm.NFig=8；%噪声系数（增加至恶化，5-10 dB）prm.N天数=500；%Frf、Fbb分区的光线数

prm.fftlength = 256;prm.cyclicprefixlength = 64;prm.numcarriers = 234;prm.nullcarrierindices = [1：7 129 256-5：256]'％守卫和直流prm.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；%每个用户的代码速率相同numtails = 6;％终止尾位数prm.numFrmBits=numSTSVec.*（prm.numDataSymbols*prm.numCarriers*...prm.bitsPerSubCarrier*编码率）-numails；prm.modMode=2^prm.bitsPerSubCarrier；%调制顺序%信道滤波器延迟的原因numpadsym = 3;零的符号数量prm.numPadZeros=numPadSym*（prm.FFTLENGHT+prm.CYCLICPREFIXLENGHT）；

prm.clight = physconst（'LightSpeed');Prm.lambda = prm.clight / prm.fc;％获取发送和接收数组信息[isTxURA，expFactorTx，isRxURA，expFactorRx]=helperArrayInfo（prm，true）；％发送天线阵列定义％数组位置和角度prm.posTx=[0；0；0]；%基站/发射阵列位置[x；y；z]，米如果伊斯特塞拉%均匀矩形阵列txarray=相控分区阵列(...'大批'，分阶段.URA（[expFactorTx numSTS]，0.5*prm.lambda），...“子阵列选择”，一个（numSTS，numTx），“亚阵列转向”,“习俗”);其他的%均匀线阵txarray=相控.ULA（numTx，“元素间距”，0.5*prm.lambda，...“元素”相位各向同性天线元件(“背靠背”，假）；终止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.URA（[expFactorRx（uIdx）numSTSVec（uIdx）]，...0.5*prm.lambda），“子阵列选择”，一个（numSTSVec（uIdx），...numrx（Uidx）），“亚阵列转向”,“习俗”)；prm.posRxElem=getElementPosition（rxarray）/prm.lambda；其他的如果numRx（uIdx）>1%均匀线阵rxarray=相控.ULA（numRx（uIdx），...“元素间距”，0.5*prm.lambda，...“元素”，相位各向同性天线元素）；prm.posRxElem=getElementPosition（rxarray）/prm.lambda；其他的rxarray=相位各向同性天线元素；prm.posRxElem=[0；0；0]；%LCS终止终止%移动位置[xRx，yRx，zRx]=sph2cart（脱胶剂（prm.mobileAngles（1，uIdx）），...deg2rad（prm.mobileAngles（2，uIdx）），...prm.mobileRanges（uIdx））；prm.posRx（：，uIdx）=[xRx；yRx；zRx]；[toRxRange，toRxAng]=射程角（prm.posTx，prm.posRx（：，uIdx））；spLoss（uIdx）=fspl（toRxRange，prm.lambda）；终止

%生成前导信号prm.numsts = numtx;%设置为numTx可发出所有频道的声音PreeAblesig = HelpergenProdeamble（PRM）;%通过信道传输前导码prm.numsts = numsts;%为通道保留相同的阵列配置[rxPreSig，chanDelay]=助手应用多通道（前置码、prm、spLoss）；%信道状态信息反馈hDp=单元（prm.numUsers，1）；prm.numSTS=numTx；%设置为numTx以估计所有链接对于uIdx=1:prm.numUsers%前端放大器增益和热噪声rxpreamp = phased.receiverProamp（...“收益”，spLoss（uIdx），...%路径损耗的原因“噪音图”，prm.nfig，“参考温度”,290,...“采样器”，prm.chanSRate）；rxPreSigAmp=rxPreAmp（rxPreSig{uIdx}）；%所用子载波的标度功率rxPreSigAmp=rxPreSigAmp*（sqrt（prm.fft长度-...长度（prm.nullcarrierindices）/ prm.fftlength）;%正交频分复用解调rxOFDM=ofdmdemod（rxPreSigAmp（chanDelay（uIdx）+1:...完（努帕泽罗斯-钱德雷公墓（uIdx）），：），公墓长度，...prm.cyclicprefixlength，prm.cyclicprefixlength，...prm.无载货方、prm.飞行员载货方）；%基于前导码的信道估计%numCarr、numTx、numRxhDp{uIdx}=helpermmochannelestimate（rxOFDM（：，1:numTx，：），prm）；终止

%计算传输侧的混合权重如果prm.numUsers==1%单用户OMP[az; el] = [ -  180：180; -90：90] 3d空间，等间距%txang=[-180:360/prm.天数：180；-90:180/prm.天数：90]；txang=[兰特（1，公升）*360-180；兰特（1，公升）*180-90]；%随机的At=steervec（prm.posTxElem，txang）；AtExp=complex（零（prm.numCarriers，大小（At，1），大小（At，2））；对于carrIdx=1:prm.numCarriers AtExp（carrIdx，：，：）=At；%所有子载波相同终止%正交匹配追踪混合权重[Fbb，Frf]=omphybweights（hDp{1}，numSTS，numSTS，AtExp）；v=Fbb；%设置基带预编码器（Fbb）%对于平坦信道，Frf在子载波上是相同的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 Nst by Nst终止v = erfute（steeringmatrix，[1 3 2]）;终止%发射阵列模式图如果伊斯特塞拉%第一子载波的URA元素响应模式（txarray，prm.fc，-180:180，-90:90，“类型”,“埃菲尔德”,...'Elementweights'挤压（v（1，：，：），...'繁殖'，prm.clight）;其他的%乌拉%第一副载波的阵列响应WTS = MFRF。'*挤压（v（1，：，:)）;模式（txarray，prm.fc，-180:180，-90:90，“类型”,“埃菲尔德”,...“重量”，wts（：，1），'繁殖'，prm.clight）;终止prm.numsts = numsts;％恢复数据传输

%卷积编码器编码器=公共卷积编码器(...“网格结构”，poly2trellis（7[133 171 165]），...“终止方法”,“终止”)；txDataBits=单元（prm.numesers，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=复数（零（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（预数据、prm.FFT长度、prm.CyclicPrefixLength、，...PRM.NULLCARRIERINDICES，PRM.PILOTCARRIERINDICES，飞行员）;%所用子载波的标度功率txOFDM=txOFDM*（prm.fft长度/...sqrt（（prm.fft长度（prm.nullcarrierndices））；%生成带有反馈权重的前导码并预先发送到数据preambleSigD=helperGenPreamble（prm，v）；txSigSTS=[preambleSigD；txOFDM]；%射频波束形成：将Frf应用于数字信号%每个天线单元连接到每个数据流txSig=TXSIGST*mFrf；

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

hfig=图(“姓名”,“每个流的均衡符号星座”);scfact =（（prm.fftlength-length-（prm.nullcarrierindices））.../prm.fftlelength^2）/numTx；nVar=noisepow（prm.chanSRate，prm.NFig，290）/scFact；解码器=comm.ViterbiDecoder('InputFormat','缺乏尺寸',...“网格结构”，poly2trellis（7[133 171 165]），...“终止方法”,“终止”,“OutputDataType”,“双人”);对于UIDX = 1：prm.numusers stsu = numstsvec（Uidx）;stsidx = sum（numstsvec（1 :( uidx-1）））+（1：stsu）;%前端放大器增益和热噪声rxpreamp = phased.receiverProamp（...“收益”，spLoss（uIdx），...%路径损耗的原因“噪音图”，prm.nfig，“参考温度”,290,...“采样器”，prm.chansrate）;rxsigamp = rxpreamp（rxsig {uidx}）;%占用子载波的缩放功率rxSigAmp=rxSigAmp*（sqrt（prm.FFT长度（prm.NullCarrierIndices））.../prm.fft长度）；%正交频分复用解调rxOFDM=ofdmdemo（rxSigAmp（chanDelay（uIdx）+1:...完（努帕泽罗斯-钱德雷公墓（uIdx）），：），公墓长度，...prm.cyclicprefixlength，prm.cyclicprefixlength，...prm.无载货方、prm.飞行员载货方）；%基于映射前导的信道估计HD = Helpermimochannelestimate（Rxofdm（：，1：numsts，:)，prm）;%MIMO均衡%为感兴趣的用户索引到流中[rxEq，CSI]=helperminoequalize（rxOFDM（：，numSTS+1:end，：），hD（：，stsIdx，：）；%软解调rxSymbs=rxEq（：）/sqrt（numTx）；rxLLRBits=qamdemod（rxSymbs，prm.modMode，'onemaveragepower'符合事实的...“输出类型”,“approxllr”,“噪音回避”，nVar）；%解码前应用CSIrxLLRtmp=重塑（rxLLRBits，prm.bitsPerSubCarrier，[]，...Prm.numdatasymbols，Stsu）;csitmp = REPAPE（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子地块（prm.numUsers，最大值（numSTSVec），（uIdx-1）*最大值（numSTSVec）+i）；绘图（重塑（rxEq（：，：，i）/sqrt（numTx），[]，1），'.')；轴线正方形xlim（gca，[-scaler-scaler]）；ylim（gca，[-scaler-scaler]）；title([“你”num2str（uIdx）“，DS”num2str（i）]）;网格在…上;终止%计算并显示EVMevm=通信evm(“正常化”,“平均星座功率”,...“引用信号源”,“根据参考星座估计”,...“参考告知”,...qammod（（0:prm.modMode-1）’，prm.modMode，'onemaveragepower',1)); rmsEVM=evm（rxSymbs）；disp([“用户”num2str（uIdx）]；disp(['RMS EVM（%）='num2str（rmsEVM）]；%计算并显示误码率ber=comm.ErrorRate；measures=ber（txDataBits{uIdx}，rxBits）；fprintf('BER=%.5f；位数=%d；错误数=%d\n',...措施（1）、措施（3）、措施（2））；终止

rng（s）；％恢复rng状态