rng (211);％设置重复性的RNG状态

人口、难民和移民事务局。NCellID = 1;%细胞IDprm.freqrange ='fr1'；％频率范围：'fr1'或'fr2'人口、难民和移民事务局。CenterFreq = 3.5 e9;%赫兹人口、难民和移民事务局。SSBlockPattern =“案例B”；%例A / B / C / D / E人口、难民和移民事务局。ssbtransmit = [ones(1,8) zeros(1,0)];长度% 4/8或64prm.txarraysize = [8 8];%发送数组大小，[rows cols]prm.txazlim = [-60 60];发射方位扫描限制人口、难民和移民事务局。TxELlim = [-90 0];%发送仰角扫描限制prm.rxarraysize = [2 2];%接收数组大小，[rows cols]prm.rxazlim = [-180 180];％接受方位角扫描限制prm.rxellim = [0 90];%接收提升扫描限制prm.elevationsweep = false;%启用/禁用提升扫描人口、难民和移民事务局。SNRdB = 30;%信噪比,dB人口、难民和移民事务局。RSRPMode ='ssswdmrs'；%{‘SSSwDMRS’,‘SSSonly}

PRM = ValidateParams（PRM）;

txburst.blockpattern = prm.ssblockpattern;txburst.ssbtransmited = prm.ssbtransmited;txburst.ncellid = prm.ncellid;txburst.ssbperiodicity = 20;txburst.nframe = 0;txburst.windowing = 0;txburst.displayburst = true;%假设载波与突发的子载波间距相同载体= nrCarrierConfig ('ncellid'人口、难民和移民事务局。NCellID,“NFrame”, txBurst.NFrame);母舰。SubcarrierSpacing = prm.SCS;carrierInfo = nrOFDMInfo(载体);txBurst。SampleRate = carrierInfo.SampleRate;

c = physconst (“光速”）;%传播速度λ= c / prm.CenterFreq;%的波长人口、难民和移民事务局。posTx = (0, 0, 0);%发射阵列位置，[x;y;z]，米人口、难民和移民事务局。posRx = (100; 50 0);%接收数组位置，[x;y;z]，米torxrange = rangeangle（prm.postx，prm.posrx）;sploss = fspl（torxrange，lambda）;%空闲空间路径损失％发送阵列如果人口、难民和移民事务局。IsTxURA％均匀矩形阵列arraytx = phased.ura（prm.txArraysize，0.5 * lambda，......“元素”,分阶段。IsotropicAntennaElement (“BackBaffled”，真的））;别的％均匀线性阵列arrayTx = phased.ULA (prm.NumTx,......“ElementSpacing”0.5 *λ......“元素”,分阶段。IsotropicAntennaElement (“BackBaffled”，真的））;结尾%接收数组如果人口、难民和移民事务局。IsRxURA％均匀矩形阵列arrayRx = phased.URA (prm.RxArraySize, 0.5 *λ,......“元素”，序列。异丙哒植物;别的％均匀线性阵列arrayRx = phased.ULA (prm.NumRx,......“ElementSpacing”0.5 *λ......“元素”，序列。异丙哒植物;结尾％散射器位置人口、难民和移民事务局。FixedScatMode = true;如果人口、难民和移民事务局。FixedScatMode%固定单一散布器位置prm.scatpos = [50;80;0];别的％在随机位置生成散射体nscat = 10;%散射体数目azrange = -180：180;Elrange = -90：90;randazorder = randperm（长度（azrange））;randelorder = randperm（长度（elrange））;azanginsph = azrange（randazorder（1：nscat））;Elanginsph = Elrange（Randelorder（1：NSCAT））;r = 20;%半径[x, y, z] = sph2cart(函数(azAngInSph),函数(elAngInSph), r);人口、难民和移民事务局。= [x;y;z] + (prm。posTx + prm.posRx) / 2;结尾%配置通道频道= phased.ScatteringMIMOChannel;通道。PropagationSpeed = c;通道。CarrierFrequency = prm.CenterFreq;通道。SampleRate = txBurst.SampleRate;通道。SimulateDirectPath = false;通道。ChannelResponseOutputPort = true; channel.Polarization =“没有”；通道。TransmitArray = arrayTx;通道。TransmitArrayPosition = prm.posTx;通道。ReceiveArray = arrayRx;通道。ReceiveArrayPosition = prm.posRx;通道。ScattererSpecificationSource =“属性”；通道。ScattererPosition = prm.ScatPos;通道。ScattererCoefficient = 1(1、大小(prm.ScatPos 2));％获得最大频道延迟[〜，〜，tau] =通道（复杂（randn（txburst.samplege * 1e-3，prm.numtx），......randn (txBurst.SampleRate * 1 e - 3, prm.NumTx)));maxChDelay =装天花板(max(τ)* txBurst.SampleRate);

%创建和显示突发信息txBurstInfo = hSSBurstInfo (txBurst);disp (txBurstInfo);%产生突发波形和网格[BurstwaveForm，TxBurstGrid] =冰爆（TXBurst）;

%发送端和接收端波束数numbeams = sum（txburst.ssbtransmited）;发射波束的方位角和仰角，等间距人口、难民和移民事务局azBW =波束宽度(arrayTx。CenterFreq,“切”，“方位”）;人口、难民和移民事务局elBW =波束宽度(arrayTx。CenterFreq,“切”，'海拔'）;txbeamang = hgetbeamweepangles（numbeams，prm.txazlim，prm.txellim，......azBW、elBW prm.ElevationSweep);%用于评估传动侧转向重量steervectx = phased.steringvector（“SensorArray”arrayTx,......“PropagationSpeed”，C）;%为每个单边带应用每个OFDM符号转向carrier.nsizegrid = txburstinfo.nrb;OFDMINFO = NROFDMINFO（载体）;gridsymlengths = repmat（Ofdminfo.symbollenghs，1，......长度尺寸(txBurstGrid, 2) / (ofdmInfo.SymbolLengths));%在numTx上重复爆发以准备转向strtxwaveform = repmat（burstwaveform，1，prm.numtx）./ sqrt（prm.numtx）;为单边带= 1:长度(txBurstInfo.SSBIndex)从突发中提取SSB波形blockSymbols = txBurstInfo.OccupiedSymbols(单边带、:);startSSBInd =总和(gridSymLengths (1: blockSymbols (1) 1)) + 1;endSSBInd =总和(gridSymLengths (1: blockSymbols (4)));ssbWaveform = strTxWaveform (startSSBInd: endSSBInd, 1);％为转向方向产生重量wT = SteerVecTx (prm.CenterFreq txBeamAng(:,单边带));%将每个发送元素的权重应用到SSBstrTxWaveform (startSSBInd: endSSBInd:) = ssbWaveform。* (wT);结尾

接收方位角和仰角，等间距AZBW = BeamWidth（ArrayRx，PRM.Centerfreq，“切”，“方位”）;人口、难民和移民事务局elBW =波束宽度(arrayRx。CenterFreq,“切”，'海拔'）;rxbeamang = hgetbeamswealls（numbeams，prm.rxazlim，prm.rxellim，......azBW、elBW prm.ElevationSweep);评估接收侧转向权重的％SteerVecRx =分阶段。SteeringVector (“SensorArray”arrayRx,......“PropagationSpeed”，C）;% AWGN水平SNR = 10 ^（prm.snrdb / 20）;%转换为线性增益N0 = 1 /(√(2.0 * prm.NumRx *双(ofdmInfo.Nfft)) *信噪比);％噪声STD。开发。%以线性形式接收增益，以补偿路径损耗rxgain = 10 ^（sploss / 20）;生成一个参考网格用于定时校正%假设SSB在第一个槽位pssRef = nrPSS (carrier.NCellID);pssInd = nrPSSIndices;pbchdmrsRef = nrPBCHDMRS (carrier.NCellID txBurstInfo.ibar_SSB (1));pbchDMRSInd = nrPBCHDMRSIndices (carrier.NCellID);pssGrid = 0 ([240 4]);pssGrid (pssInd) = pssRef;pssGrid (pbchDMRSInd) = pbchdmrsRef;refGrid = 0([12 *运营商。NSizeGrid ofdmInfo.SymbolsPerSlot]); refGrid(txBurstInfo.OccupiedSubcarriers,......txBurstInfo.OccupiedSymbols (1)) = pssGrid;所有接收波束的％循环rsrp = 0 (numBeams numBeams);为RIDX = 1：numbeams衰落信道，带有路径损耗txWave = [strTxWaveform;0 (maxChDelay、尺寸(strTxWaveform 2)));fadWave =通道(txWave);%接收增益，以补偿路径损失fadWaveG = fadWave * rxGain;％添加WGN.噪声= N0 *复合物（Randn（大小（Fadwaveg）），Randn（大小（Fadwaveg）））;rxwaveform = fadwaveg +噪声;％为转向方向产生重量或者说是= SteerVecRx (prm.CenterFreq rxBeamAng (:, rIdx));％每次接收元素施加权重如果Strcmp（prm.freqrange，'fr1'）strrxwaveform = rxwaveform。*（wr'）;别的%对于FR2，合并来自天线元件的信号strRxWaveform = rxWaveform *连词(wR);结尾%的正确时机抵消= nrTimingEstimate(载体,......strrxwaveform（1：Ofdminfo.samplerate * 1E-3，:)，Refgrid * WR（1）'）;如果偏移> maxchdelay offset = 0;结尾strrxwaveforms = strrxwaveform（1 + offset：结束，:);% OFDM解调rxgrid = nrofdmdemodulate（载体，strrxwaveforms）;rxGrid中所有ssb的环路(发送端)为tIdx = 1: numBeams%得到每个SSB网格rxSSBGrid = rxGrid (txBurstInfo。OccupiedSubcarriers,......txburstinfo.occupiedsymbols（tidx，:)，:);％进行测量，每次接收存储，传输波束rsrp (rIdx tIdx) = measureSSB (rxSSBGrid、prm.RSRPMode txBurst.NCellID);结尾结尾

[m，i] = max（rsrp，[]，'全部'，'线性'）;%第一次出现是输出% I首先是列向下(接收)，然后是跨列(传输)[rxBeamID,txBeamID] = ind2sub([numBeams numBeams]，i(1));显示选定的光束对DISP（['与RSRP的选定梁配对：'num2str（10 * log10（rsrp（rxbeamid，......TXBeamid））+ 30）'dbm'，13.'传输＃'num2str (txBeamID)......'（方位角：'txBeamID num2str (txBeamAng (1))'，海拔：'......txBeamID num2str (txBeamAng (2))“)”13“收到#”num2str (rxBeamID)......'（方位角：'rxBeamID num2str (rxBeamAng (1))'，海拔：'......rxBeamID num2str (rxBeamAng (2))“)”]）;％显示最终光束对模式h =图(“位置”，图（[32 55 32 40]），'菜单栏'，'没有任何'）;H.Name =.“选定的传输阵列响应模式”；wT = SteerVecTx (prm.CenterFreq txBeamAng (:, txBeamID));模式(arrayTx,人口、难民和移民事务局。CenterFreq,“PropagationSpeed”，C，“重量”、wT);h =图(“位置”，图（[32 55 32 40]），'菜单栏'，'没有任何'）;H.Name =.“选定的接收阵列响应模式”；或者说是= SteerVecRx (prm.CenterFreq rxBeamAng (:, rxBeamID));模式(arrayRx,人口、难民和移民事务局。CenterFreq,“PropagationSpeed”，C，“重量”,或者说是);%用tx, rx，散射体和确定的波束绘制MIMO场景prmscene = struct（）;prmscene.txarraysize = prm.txArrassize;prmscene.rxArraysize = prm.rxArrassize;prmscene.txelempos = getElementPosition（arraytx）;％仪表prmscene.rxelempos = getElementPosition（arrayrx）;％仪表prmScene。txArrayPos = prm.posTx;prmScene。rxArrayPos = prm.posRx;prmScene。txAzAngles = 90:90;prmScene。rxAzAngles = [90:180 -179:-90];prmScene。scatPos = prm.ScatPos; prmScene.lambda = lambda; prmScene.arrayScaling = 1; hPlotSpatialMIMOScene(prmScene,wT,wR);如果人口、难民和移民事务局。ElevationSweep视图(2);结尾

函数PRM = ValidateParams（PRM）％验证用户指定的参数并返回更新的参数％％仅交叉依赖检查是参数一致性的。如果strcmpi(人口、难民和移民事务局。FreqRange,'fr1'）如果prm.centerfreq> 7.125e9 ||prm.centerfreq <410E6错误（['指定的中心频率超出FR1'，......'频率范围（410 MHz  -  7.125 GHz）。]）;结尾如果Strcmpi（prm.ssblockpattern，'案例d'）||......Strcmpi（prm.ssblockpattern，“案例e”） 错误（['选择FR1频率的SSBLOCKPattern无效'......的范围内。SSBlockPattern必须是“Case A”或“......“‘情况B’或FR1的情况C。”]）;结尾如果~ (((prm.SSBTransmitted)长度= = 4)| |......(长度(prm.SSBTransmitted) = = 8))错误([' ssbtransmit必须是长度为4或8的矢量'，......“用于FR1频率范围。”]）;结尾如果(人口、难民和移民事务局。CenterFreq <= 3e9) && (length(prm. ssbtransmit)~=4) error([' ssbtransmit必须是长度为4的矢量'......'中心频率小于或等于3GHz。]）;结尾如果(人口、难民和移民事务局。CenterFreq > 3e9) && (length(prm. ssbtransmission)~=8) error([' ssbtransmit必须是长度为8的矢量'，......“中心频率大于3GHz且小于”，......“或等于7.125GHz。”]）;结尾别的％'fr2'如果人口、难民和移民事务局。CenterFreq > 52.6e9 || prm。CenterFreq < 24.25e9 error(['指定的中心频率超出FR2'，......频率范围(24.25 GHz - 52.6 GHz)]）;结尾如果〜（strcmpi（prm.ssblockpattern，'案例d'）||......Strcmpi（prm.ssblockpattern，“案例e”)错误([“选定FR2频率的ssblock模式无效”......的范围内。ssblockpattern必须是''case d''或'......'''案例e''f​​or fr2。]）;结尾如果长度(prm.SSBTransmitted) ~ = 64错误(['ssbtransmited必须是'的矢量长度为64，......FR2频率范围。]）;结尾结尾prm.numtx = prod（prm.txArraysize）;prm.numrx = prod（prm.rxArraysize）;如果人口、难民和移民事务局。NumTx==1 || prm.NumRx==1 error([“发射或接收天线元件的数量必须为”，......'大于1.']）;结尾人口、难民和移民事务局。IsTxURA = (prm.TxArraySize(1)>1) && (prm.TxArraySize(2)>1);人口、难民和移民事务局。IsRxURA = (prm.RxArraySize(1)>1) && (prm.RxArraySize(2)>1);如果~ (strcmpi(人口、难民和移民事务局。RSRPMode,“SSSonly”）||......strcmpi(人口、难民和移民事务局。RSRPMode,'ssswdmrs')错误(['无效的RSRP测量模式。指定“，......''SSOSOLLY''或''SSSWDMRS''模式。']）;结尾％基于SSBLOCKPattern选择SCS开关低(prm.SSBlockPattern)情况下的情况下scs = 15;情况下｛“案例b”，'案例c'} sc = 30;情况下'案例d'scs = 120;情况下'案例e'scs = 240;结尾人口、难民和移民事务局。SCS = SCS;结尾函数NCellID rsrp = measureSSB (rxSSBGrid模式)%根据SSS和if计算参考信号接收功率(RSRP)%选择，也PBCH DM-RS。sssind = nrsssindices;% SSS指数numRx =大小(rxSSBGrid, 3);rsrpSSS = 0 (numRx, 1);为rxIdx = 1: numRx每个RX元素提取信号rxSSBGridperRx = rxSSBGrid (:,:, rxIdx);rxSSS = rxSSBGridperRx (sssInd);% RS的所有REs的平均功率贡献RSRPSSS（rxidx）=均值（rxsss。*结合（rxsss））;结尾如果Strcmpi（模式，'ssswdmrs'）pbchdmrsind = nrpbchdmrsindices（ncellid）;％PBCH DM-RS指数rsrpDMRS = 0 (numRx, 1);为rxIdx = 1: numRx每个RX元素提取信号rxSSBGridperRx = rxSSBGrid (:,:, rxIdx);rxPBCHDMRS = rxSSBGridperRx (pbchDMRSInd);% RS的所有REs的平均功率贡献rsrpDMRS (rxIdx) = (rxPBCHDMRS意思。*连词(rxPBCHDMRS));结尾结尾开关较低（模式）情况下'ssonly'％只有SSSrsrp = max (rsrpSSS);%最大接收元素情况下“ssswdmrs”％SSS和PBCH-DMRS，每个RS的计计rsrp = max ((rsrpSSS * 127 + rsrpDMRS * 144) / 271);%最大接收元素结尾结尾