S = rng(67);为可重复性设置RNG状态

%多用户系统，每个用户有单个/多个流prm.numUsers= 4;%用户数量prm.numSTSVec= [3 2 1 2];%每个用户独立数据流的个数prm.numSTS= sum(prm.numSTSVec);%必须是2的幂人口、难民和移民事务局。numTx = prm.numSTS*8;% BS发射天线数(2的幂)人口、难民和移民事务局。numRx = prm.numSTSVec*4;每个用户接收天线数量% (any >= numSTSVec)每个用户具有相同的调制人口、难民和移民事务局。bitsPerSubCarrier = 4;% 2: qpsk, 4:16 qam, 6:64 qam, 8:256qamprm.numDataSymbols= 10;% OFDM数据符号数% MS位置:假设原点为BS以[方位角;仰角]度表示的角度。% az in range [-180 180]， el in range[-90 90]，例如[45;0]maxRange = 1000;%在BS 1000米范围内的所有MSs人口、难民和移民事务局。mobileRanges = randi([1 maxRange]，1,prm.numUsers);人口、难民和移民事务局。mobileAngles = [rand(1,prm.numUsers)*360-180;.．.兰特(prm.numUsers) * 180 - 90);人口、难民和移民事务局。Fc = 28e9;% 28 GHz系统人口、难民和移民事务局。chanSRate = 100e6;通道采样率，100 Msps人口、难民和移民事务局。ChanType =“散射”；通道选项:“散射”，“MIMO”人口、难民和移民事务局。NFig = 8;%噪音数值(增加至恶化，5-10分贝)prm.nRays= 500;用于Frf, Fbb分区的射线数

人口、难民和移民事务局。FFTLength = 256;人口、难民和移民事务局。CyclicPrefixLength = 64;人口、难民和移民事务局。numCarriers = 234;人口、难民和移民事务局。NullCarrierIndices = [1:7 129 256-5:256]';%卫兵和DC人口、难民和移民事务局。PilotCarrierIndices = [26 54 90 118 140 168 204 232]';nonDataIdx = [prm.NullCarrierIndices;prm.PilotCarrierIndices];人口、难民和移民事务局。CarriersLocations = setdiff((1:prm.FFTLength)'， sort(nonDataIdx));numSTS = prm.numSTS;numTx = prm.numTx;numRx = prm.numRx;numSTSVec = prm.numSTSVec; codeRate = 1/3;%每个用户相同的码率numTails = 6;%终止尾位数人口、难民和移民事务局。numFrmBits = numSTSVec.*(prm.numDataSymbols*prm.numCarriers* ..．.prm.bitsPerSubCarrier * codeRate) -numTails;人口、难民和移民事务局。modMode = 2^prm.bitsPerSubCarrier;调制阶数考虑信道滤波器延迟numPadSym = 3;到零opad的符号数量的%人口、难民和移民事务局。numPadZeros = numPadSym*(prm.FFTLength+prm.CyclicPrefixLength);

人口、难民和移民事务局。cLight = physconst(“光速”）;人口、难民和移民事务局。lambda = prm.cLight/prm.fc;获取发送和接收数组信息[isTxURA,expFactorTx,isRxURA,expFactorRx] = helperArrayInfo(prm,true);发射天线阵列定义%阵列位置和角度人口、难民和移民事务局。posTx = [0;0;0];% BS/发射阵列位置，[x;y;z]，米如果isTxURA均匀矩形阵列Txarray =相控。PartitionedArray (.．.“数组”,分阶段。([expFactorTx numSTS], 0.5 * prm.lambda),.．.“SubarraySelection”的(numSTS numTx),“SubarraySteering”，“自定义”）;其他的均匀线阵Txarray =相控。齿龈(numTx“ElementSpacing”0.5 * prm.lambda.．.“元素”,分阶段。IsotropicAntennaElement (“BackBaffled”、假));结束人口、难民和移民事务局。posTxElem = getElementPosition(txarray)/prm.lambda;spLoss = 0 (prm.numUsers,1);人口、难民和移民事务局。posRx = 0 (3,prm.numUsers);为uIdx = 1:prm.numUsers%接收数组如果isRxURA (uIdx)均匀矩形阵列Rxarray =相控。PartitionedArray (.．.“数组”, phased.URA ([expFactorRx (uIdx) numSTSVec (uIdx)],.．.0.5 * prm.lambda),“SubarraySelection”的(numSTSVec (uIdx),.．.numRx (uIdx)),“SubarraySteering”，“自定义”）;人口、难民和移民事务局。posRxElem = getElementPosition(rxarray)/prm.lambda;其他的如果numRx (uIdx) > 1均匀线阵rxarray = phase . ula (numRx(uIdx)，.．.“ElementSpacing”0.5 * prm.lambda.．.“元素”, phased.IsotropicAntennaElement);人口、难民和移民事务局。posRxElem = getElementPosition(rxarray)/prm.lambda;其他的rxarray = phase . isotropicantennaelement;人口、难民和移民事务局。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以发出所有通道preambleSig = helperGenPreamble(prm);通过信道传输前导prm.numSTS= numSTS;%保持相同的阵列配置通道[rxPreSig,chanDelay] = helperApplyMUChannel(preambleSig,prm,spLoss);通道状态信息反馈hDp = cell(prm.numUsers,1);prm.numSTS= numTx;%设置为numTx估计所有链接为uIdx = 1:prm.numUsers前端放大器增益和热噪声rxPreAmp =相控的。ReceiverPreamp (.．.“获得”spLoss (uIdx),.．.%考虑路径丢失“NoiseFigure”人口、难民和移民事务局。NFig,“ReferenceTemperature”, 290,.．.“SampleRate”, prm.chanSRate);rxPreSigAmp = rxPreAmp(rxPreSig{uIdx});所使用子载波的%比例功率rxPreSigAmp = rxPreSigAmp *(根号(prm。FFTLength -.．.长度(prm.NullCarrierIndices)) / prm.FFTLength);% OFDM解调rxOFDM = ofdmdemod(rxPreSigAmp(chanDelay(uIdx)+1:.．.结束(prm.numPadZeros-chanDelay (uIdx)):), prm.FFTLength,.．.prm.CyclicPrefixLength prm.CyclicPrefixLength,.．.prm.NullCarrierIndices prm.PilotCarrierIndices);%来自序言的信道估计% 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;rand(1,prm.nRays)*180-90];%随机At = steervec(prm.posTxElem,txang);AtExp =复杂(零(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 = permute(mean(Frf,1)，[2 3 1]);其他的多用户联合空间分复用[Fbb,mFrf] = helperJSDMTransmitWeights(hDp,prm);%多用户基带预编码将每个用户的CSI打包成一个矩阵(块对角线)steeringMatrix = 0 (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 = permute(steeringMatrix，[1 3 2]);结束传输阵列模式图如果isTxURA第一个子载波的URA元素响应模式(txarray、prm.fc 180:180 90:90,“类型”，“efield”，.．.“ElementWeights”, mFrf。*挤压(v (1::)),.．.“PropagationSpeed”, prm.cLight);其他的%齿龈第一个子载波的阵列响应wts = mFrf.'*挤压(v(1，:，:));模式(txarray、prm.fc 180:180 90:90,“类型”，“efield”，.．.“重量”出世(:1),“PropagationSpeed”, prm.cLight);结束prm.numSTS= numSTS;%恢复数据传输

%卷积编码器encoder = com . 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 =编码器(txDataBits{uIdx});位到QAM符号的映射mappedSym = qammod(encodedBits,prm.modMode，“InputType”，“一点”，.．.“UnitAveragePower”,真正的);映射到层:每个用户、每个符号、每个数据流stsIdx = sum(numSTSVec(1:(uIdx-1)))+(1:numSTSVec(uIdx));gridData(:，:，stsIdx) =重塑(mappedSym,prm.numCarriers，.．.prm.numDataSymbols numSTSVec (uIdx));结束对子载波应用预编码权重，假设完全反馈preData =复杂(零(prm.numCarriers,prm.numDataSymbols,numSTS));为symIdx = 1:prm.numDataSymbols为carrIdx = 1:prm。numCarriers Q = squeeze(v(carrIdx，:，:));normQ = Q *√(numTx)/norm(Q，“摇来摇去”）;preData(carrIdx,symIdx，:) = squeeze(gridData(carrIdx,symIdx，:))).'.．.* normQ;结束结束%多天线飞行员飞行员= 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))));生成带有反馈权重和数据前缀的序文preambleSigD = helperGenPreamble(prm,v);txSigSTS = [preambleSigD;txOFDM];% RF波束形成:将RF应用于数字信号每个天线元件连接到每个数据流txSig = txSigSTS*mFrf;

将空间定义的通道应用到传输信号[rxSig,chanDelay] = helperApplyMUChannel(txSig,prm,spLoss,preambleSig);

Hfig = figure(“名字”，“每个流的均衡符号星座”）;scFact = ((prm.FFTLength-length(prm.NullCarrierIndices)).．./ prm.FFTLength ^ 2) / numTx;nVar = noisepow(prm.chanSRate,prm.NFig,290)/scFact;解码器= com . viterbidecoder (“InputFormat”，“Unquantized”，.．.“TrellisStructure”，poly2trellis(7， [133 171 165])，.．.“TerminationMethod”，“终止”，“OutputDataType”，“双”）;为uIdx = 1:prm.numUsersstsU = numSTSVec(uIdx); stsIdx = sum(numSTSVec(1:(uIdx-1)))+(1:stsU);前端放大器增益和热噪声rxPreAmp =相控的。ReceiverPreamp (.．.“获得”spLoss (uIdx),.．.%考虑路径丢失“NoiseFigure”人口、难民和移民事务局。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.NullCarrierIndices prm.PilotCarrierIndices);来自映射序文的信道估计hD = helperMIMOChannelEstimate(rxOFDM(:，1:numSTS，:)，prm);% MIMO均衡为感兴趣的用户索引到流中[rxEq, CSI] = helperMIMOEqualize (rxOFDM (numSTS + 1::,:),高清(:,stsIdx:));%软解调rxSymbs = rxEq(:)/sqrt(numTx);rxLLRBits = qamdemod(rxSymbs,prm.modMode，“UnitAveragePower”,真的,.．.“OutputType”，“approxllr”，“NoiseVariance”据nVar);在解码之前应用CSIrxLLRtmp =重塑(rxLLRBits,prm.bitsPerSubCarrier，[]，.．.prm.numDataSymbols stsU);csitmp =重塑(CSI,1，[]，1,numSTSVec(uIdx));rxScaledLLR = rxLLRtmp.*csitmp;%软输入通道解码rxDecoded =解码器(rxScaledLLR(:));%已解码的接收位rxBits = rxDecoded(1:prm.numFrmBits(uIdx));为每个用户的所有流绘制均衡符号标量= ceil(max(abs([real(rxSymbs(:)));图像放大(rxSymbs (:))))));为i = 1:stsU subplot(prm。numUsers, max(numSTSVec)， (uIdx-1)*max(numSTSVec)+i);plot(重塑(rxEq(:，:，i)/sqrt(numTx)， []， 1)，“。”）;轴广场甘氨胆酸xlim(、(标量标量));甘氨胆酸ylim(、(标量标量));标题([“U”num2str (uIdx)'， ds 'num2str (i)));网格在；结束计算并显示EVMevm = com . evm“归一化”，“平均星座功率”，.．.“ReferenceSignalSource”，“从参考星座估计”，.．.“ReferenceConstellation”，.．.qammod ((0: prm.modMode-1)、prm.modMode“UnitAveragePower”1));rmsEVM = evm(rxSymbs);disp ([“用户”num2str (uIdx)]);disp ([' RMS evm (%) = 'num2str (rmsEVM)]);计算和显示误码率ber = com . errorrate;measures = ber(txDataBits{uIdx}，rxBits);流(' BER = %.5f;不。of Bits = %d;不。误差= %d\n'，.．.(1),(3)措施的人,措施(2));结束

rng(年代);%恢复RNG状态