rng（211）；为可重复性设置RNG状态

prm.NCellID=1；%细胞IDPRM.FreqRange =.'fr1'；%频率范围:'FR1'或'FR2'人口、难民和移民事务局。CenterFreq = 3.5 e9;%赫兹prm.ssblockpattern =“案例B”；%例A / B / C / D / Eprm.ssbtransmited = [（1,8）零（1,0）];长度% 4/8或64prm.txarraysize = [8 8];%发送数组大小，[rows cols]人口、难民和移民事务局。TxAZlim = [-60 60];%发射方位扫描极限prm.TxELlim=[-90 0]；%传输高程扫描限制人口、难民和移民事务局。RxArraySize = [2 2];%接收数组大小，[rows cols]prm.rxazlim = [-180 180];％接受方位角扫描限制prm.rxellim = [0 90];%接收提升扫描限制prm.elevationsweep = false;%启用/禁用提升扫描prm.SNRdB=30；%信噪比,dBprm.rsrpmode ='ssswdmrs'；%{'SSSwDMRS'，'SSSonly'}

PRM = ValidateParams（PRM）;

txBurst。BlockPattern = prm.SSBlockPattern;txBurst。SSBTransmitted = prm.SSBTransmitted;txBurst。NCellID = prm.NCellID; txBurst.SSBPeriodicity = 20; txBurst.NFrame = 0; txBurst.Windowing = 0; txBurst.DisplayBurst = true;%假设载波与突发的子载波间距相同carrier = nrcarrierconfig（“NCellID”，prm.NCellID，“NFrame”, txBurst.NFrame);母舰。SubcarrierSpacing = prm.SCS;carrierInfo = nrOFDMInfo(载体);txBurst。SampleRate = carrierInfo.SampleRate;

c=物理常数(“光速”）;%传播速度λ=c/prm.CenterFreq；%波长prm.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，......“元素”相位各向同性天线元件('底板',真的));其他的％均匀线性阵列arrayTx=相控.ULA（prm.NumTx，......'ElementsPacing'，0.5 * lambda，......“元素”相位各向同性天线元件('底板',真的));结束%接收数组如果人口、难民和移民事务局。IsRxURA均匀矩形阵列arrayrx = phased.ura（prm.rxArraysize，0.5 * lambda，......“元素”，序列。异丙哒植物;其他的％均匀线性阵列arrayrx = phased.ula（prm.numrx，......'ElementsPacing'，0.5 * lambda，......“元素”，序列。异丙哒植物;结束％散射器位置prm.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（deg2rad（azAngInSph），deg2rad（elAngInSph），r）；prm.ScatPos=[x；y；z]+（prm.posTx+prm.posRx）/2；结束%配置通道channel=phased.ScatteringMIMOChannel；channel.PropagationSpeed=c；channel.CarrierFrequency=prm.CenterFreq；channel.SampleRate=txBurst.SampleRate；channel.SimulatedDirectPath=false；channel.ChannelResponseOutput=true；channel.Polarization='没有任何'；channel.transmitarray = arraytx;channel.transmitarrayposition = prm.postx;channel.receivearray = arrayrx;channel.receiveArrayPosition = PRM.POSRX;channel.scattererspecifiedsource =.“属性”；channel.scattererposition = prm.scatpos;channel.scatterercoeffdiod = soneg（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 =牛伯斯坦州（TXBURST）;disp（txburstinfo）;%产生突发波形和网格[burstWaveform, txBurstGrid] = hSSBurst (txBurst);

%发送端和接收端的波束数numBeams =总和(txBurst.SSBTransmitted);%发射波束的方位角和仰角，等间距人口、难民和移民事务局azBW =波束宽度(arrayTx。CenterFreq,“切”,“方位角”); elBW=波束宽度（阵列频率、prm.CenterFreq、，“切”,“高度”）;txbeamang = hgetbeamweepangles（numbeams，prm.txazlim，prm.txellim，......AZBW，ELBW，PRM.ELEVATIONSWEEV）;评估发射侧转向权重的％SteerVecTx =分阶段。SteeringVector (“传感器阵列”arrayTx,......“PropagationSpeed”c);%为每个单边带应用每个OFDM符号转向carrier.nsizegrid = txburstinfo.nrb;OFDMINFO = NROFDMINFO（载体）;gridsymlengths = repmat（Ofdminfo.symbollenghs，1，......尺寸（txburstgrid，2）/长度（Ofdminfo.symbollenge））;%在numTx上重复爆发以准备转向strTxWaveform = repmat (burstWaveform 1 prm.NumTx)。/√(prm.NumTx);对于ssb=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，prm.CenterFreq，“切”,“高度”）;rxbeamang = hgetbeamswealls（numbeams，prm.rxazlim，prm.rxellim，......AZBW，ELBW，PRM.ELEVATIONSWEEV）;评估接收侧转向权重的％steervecrox = phased.steringvector（“传感器阵列”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=nrpssindex；pbchdmrsRef=nrPBCHDMRS（carrier.NCellID，txBurstInfo.ibar_SSB（1））；pbchDMRSInd=nrPBCHDMRSIndices（carrier.NCellID）；pssGrid=零（[240 4]）；pssGrid（pssInd）=pssRef；pssGrid（pBchdmrsrefgrid）=pbchdmrsRef；refGrid=零（[12*carrier.NSizeGrid.SymbolsPerSlot]）；refGrid（txBurstInfo.OccupiedSubcarriers，......txburstinfo.occupiedsymbols（1，:)）= pssgrid;%循环所有接收波束rsrp = zeros（numbeams，numbeams）;对于rIdx = 1: numBeams%具有路径损耗的衰落信道txwave = [strtxwaveform;零（maxchdelay，尺寸（strtxwaveform，2））];fadwave =通道（TxWave）;％接收收益，以补偿路径损耗fadWaveG = fadWave * rxGain;％添加WGN.噪音= N0 *复杂(randn(大小(fadWaveG)), randn(大小(fadWaveG)));rx波形= fadWaveG +噪声;%为被引导的方向生成权重WR = STEVECRX（PRM.CONERERFREQ，RXBeamang（：，RIDX））;％每次接收元素施加权重如果比较字符串(人口、难民和移民事务局。FreqRange,'fr1'）strrxwaveform = rxwaveform。*（wr'）;其他的%对于FR2，合并来自天线元件的信号strrxwaveform = rxwaveform *结合（wr）;结束%正确的时机偏移量=nrTimingEstimate（载波，......strRxWaveform (1: ofdmInfo.SampleRate * 1 e - 3:), refGrid * wR (1) ');如果offset > maxChDelay offset = 0;结束strRxWaveformS = strRxWaveform(1 +抵消:最终,);％OFDM解调rxgrid = nrofdmdemodulate（载体，strrxwaveforms）;%rxGrid中所有SSB上的环路（传输端）对于tidx = 1：numbeams％获取每个SSB网格rxssbgrid = rxgrid（txburstinfo.occupiedsubcarriers，......txburstinfo.occupiedsymbols（tidx，:)，:);进行测量，存储每个接收和发送波束RSRP（RIDX，TIDX）= MeasureSB（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）......'（方位角：'num2str（txBeamAng（1，txBeamID））的海拔高度:......Num2str（TXBeamang（2，TXBeamid））“)”13.“收到#”num2str (rxBeamID)......'（方位角：'rxBeamID num2str (rxBeamAng (1))的海拔高度:......num2str（rxBeamAng（2，rxBeamID））“)”]）;％显示最终光束对模式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 =.“选定接收阵列响应模式”；wR=SteerVecRx（prm.CenterFreq，rxBeamAng（：，rxBeamID））；模式（arrayRx，prm.CenterFreq，“PropagationSpeed”c“重量”，wr）;%用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（prm.freqrange，'fr1')如果人口、难民和移民事务局。CenterFreq > 7.125e9 || prm。CenterFreq < 410e6错误(['指定的中心频率超出FR1',......'频率范围（410 MHz  -  7.125 GHz）。]）;结束如果strcmpi(人口、难民和移民事务局。SSBlockPattern,'案例d'）||......strcmpi(人口、难民和移民事务局。SSBlockPattern,“案例e”） 错误（[“选定FR1频率的ssblock模式无效”......'范围。SSBlockPattern必须是'Case A'或'......''案例b''或''case c''for fr1。']）;结束如果~ (((prm.SSBTransmitted)长度= = 4)| |......（长度（prm.SSB传输）==8）错误([' ssbtransmit必须是长度为4或8的矢量',......“用于FR1频率范围。”]）;结束如果(人口、难民和移民事务局。CenterFreq <= 3e9) && (length(prm. ssbtransmit)~=4) error([' ssbtransmit必须是长度为4的矢量'......“中心频率小于或等于3GHz。”]）;结束如果（prm.CenterFreq>3e9）和（长度（prm.SSB传输）~=8）错误([' ssbtransmit必须是长度为8的矢量',......“中心频率大于3GHz且小于”,......“或等于7.125GHz。”]）;结束其他的％'fr2'如果prm.centerfreq> 52.6e9 ||prm.centerfreq <24.25e9错误（['指定的中心频率超出FR2',......频率范围(24.25 GHz - 52.6 GHz)]）;结束如果~ (strcmpi(人口、难民和移民事务局。SSBlockPattern,'案例d'）||......strcmpi(人口、难民和移民事务局。SSBlockPattern,“案例e”))错误([“选定FR2频率的ssblock模式无效”......的范围内。SSBlockPattern必须是“Case D”或“......“Case E”代表FR2。”]）;结束如果长度（prm.sbtransmited）〜= 64错误（['ssbtransmited必须是'的矢量长度为64,......FR2频率范围。]）;结束结束人口、难民和移民事务局。NumTx = prod(prm.TxArraySize); prm.NumRx = prod(prm.RxArraySize);如果prm.NumTx==1 | | prm.NumRx==1错误(['发送或接收天线元素的数量必须是',......'大于1.']）;结束人口、难民和移民事务局。IsTxURA = (prm.TxArraySize(1)>1) && (prm.TxArraySize(2)>1);人口、难民和移民事务局。IsRxURA = (prm.RxArraySize(1)>1) && (prm.RxArraySize(2)>1);如果~ (strcmpi(人口、难民和移民事务局。RSRPMode,'ssonly'）||......strcmpi（prm.RSRPMode，'ssswdmrs'） ） 错误（['无效的RSRP测量模式。指定“,......“‘SSSonly’或‘SSSwDMRS’作为模式。”]）;结束%根据SSBlockPattern选择SCS转变低(prm.SSBlockPattern)情况下“案例a”scs = 15;情况下{“案例b”,“案例c”} scs = 30;情况下“d”scs = 120;情况下“e”scs=240；结束prm.scs = scs;结束功能rsrp=measureSSB（rxSSBGrid，模式，NCellID）％基于SSS计算电源（RSRP）的参考信号，如果%已选择，也可选择PBCH DM-RS。sssInd = nrSSSIndices;％SSS指数numRx=大小（rxSSBGrid，3）；rsrpSSS=零（numRx，1）；对于rxidx = 1：numrx提取每个rx元素的信号rxssbgridperrx = rxssbgrid（：：，rxidx）;rxsss = rxssbgridperrx（sssind）;对所有RS的I％的平均功率贡献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的I％的平均功率贡献rsrpDMRS (rxIdx) = (rxPBCHDMRS意思。*连词(rxPBCHDMRS));结束结束转变低(模式)情况下'ssonly'％只有SSSrsrp = max (rsrpSSS);%最大接收单元数情况下'ssswdmrs'％SSS和PBCH-DMRS，每个RS的计计rsrp = max ((rsrpSSS * 127 + rsrpDMRS * 144) / 271);%最大接收单元数结束结束