numPRACHSlots = 10;%的开环名额来模拟在每个信噪比SNRdB = [-21、-16、-11、6、1];在dB %信噪比范围foffset = 400.0;%在赫兹频率偏移timeErrorTolerance = 2.55;%在微秒时间误差公差阈值= [];%检测阈值prachEnabled = true;%使开环传输。模拟假警报测试,禁用开环传输。

载体= nrCarrierConfig;母舰。SubcarrierSpacing = 15;母舰。NSizeGrid = 25;%计算OFDM-related信息ofdmInfo = nrOFDMInfo(载体);

%设置开环配置开环= nrPRACHConfig;开环。FrequencyRange =“FR1”;%频率范围开环。DuplexMode =“FDD”;%频分双工(FDD)开环。ConfigurationIndex = 27个;% 0配置指数格式开环。SubcarrierSpacing = 1.25;%副载波间距开环。SequenceIndex = 22;%逻辑序列索引开环。PreambleIndex = 32;%序言指数开环。RestrictedSet =“UnrestrictedSet”;%正常模式开环。FrequencyStart = 0;%的频率位置%定义的价值ZeroCorrelationZone表存储在使用nc% nrPRACHConfig对象开关prach.Format情况下{' 0 ',' 1 ',' 2 '}ncsTable = nrPRACHConfig.Tables.NCSFormat012;ncsTableCol =(字符串(ncsTable.Properties.VariableNames) = = prach.RestrictedSet);情况下“3”ncsTable = nrPRACHConfig.Tables.NCSFormat3;ncsTableCol =(字符串(ncsTable.Properties.VariableNames) = = prach.RestrictedSet);否则ncsTable = nrPRACHConfig.Tables.NCSFormatABC;ncsTableCol =包含(string (ncsTable.Properties.VariableNames) num2str (prach.LRA));结束nc = 13;zeroCorrelationZone = ncsTable.ZeroCorrelationZone (ncsTable {:, ncsTableCol} = = nc);开环。ZeroCorrelationZone = ZeroCorrelationZone;%循环移位指数

频道= nrTDLChannel;通道。DelayProfile =“TDLC300”;%延迟概要通道。MaximumDopplerShift = 100.0;%最大赫兹的多普勒频移通道。SampleRate = ofdmInfo.SampleRate;在赫兹%输入信号采样率通道。MIMOCorrelation =“低”;% MIMO相关通道。TransmissionDirection =“上行”;%上行传输通道。NumReceiveAntennas = 2;%的接收天线通道。NormalizePathGains = true;%正常化延迟概要文件的权力通道。种子= 42;%通道种子。改变这种对不同通道实现通道。NormalizeChannelOutputs = true;%接收天线的正常化%获得信道特性的信息channelInfo =信息(渠道);

%初始化变量存储在每个信噪比检测概率pDetection = 0(大小(SNRdB));%存储配置参数需要生成开波形waveconfig。NumSubframes = prach.SubframesPerPRACHSlot;waveconfig。窗口= [];waveconfig。运营商=载体;waveconfig.PRACH。配置=开环;waveconfig.PRACH。使= prachEnabled;%的临时变量‘prach_init’,‘waveconfig_init’,‘ofdmInfo_init’,%和“channelInfo_init”是用来创建临时变量%‘开’,‘waveconfig’,‘ofdmInfo’,‘channelInfo信噪比内循环%创建独立的并行仿真实例prach_init =开环;waveconfig_init = waveconfig;ofdmInfo_init = ofdmInfo;channelInfo_init = channelInfo;为snrIdx = 1:元素个数(SNRdB)%的评论进行并行计算% parfor snrIdx = 1:元素个数为并行计算(SNRdB) %取消%减少总模拟时间,你可以执行这个循环%并行使用并行计算工具箱。注释掉”的%声明和取消“parfor”声明。如果并行计算%的工具箱(TM)没有安装,“parfor”默认为正常”的声明%在命令窗口中显示进展timeNow = char (datetime (“现在”,“格式”,“HH: mm: ss”));流([timeNow':模拟信噪比= % + 5.1 f dB……”),SNRdB (snrIdx));%设置随机数发生器设置为默认值rng (“默认”);%初始化变量的信噪比,需要初始化%的变量在使用并行计算工具箱开环= prach_init;waveconfig = waveconfig_init;ofdmInfo = ofdmInfo_init;channelInfo = channelInfo_init;%复位频道,这样每个信噪比点会经历相同的%通道实现重置(渠道);%正常化噪声功率占采样率,这是一个%传输线大小的函数用于OFDM调制。信噪比的定义% /载波资源元素对于每一个接收天线。信噪比= 10 ^ (SNRdB (snrIdx) / 10);N0 = 1 /√(2.0 * channel.NumReceiveAntennas *双(ofdmInfo.Nfft) *信噪比);%检测序言计数detectedCount = 0;%为每个开槽循环numActivePRACHSlots = 0;为nSlot = 0: numPRACHSlots-1%为当前生成开波形槽开环。NPRACHSlot = nSlot;waveconfig.PRACH.Config。NPRACHSlot = nSlot;[波形,~,winfo] = hNRPRACHWaveformGenerator (waveconfig);%设置是开在微秒时间抵消/ TS 38.141 - 1图% 8.4.1.4.2-2和图8.4.1.4.2-3如果prach.LRA = = 839%长序言,如图8.4.1.4.2-2值baseOffset = ((winfo.WaveformResources.PRACH.Resources.PRACHSymbolsInfo.NumCyclicShifts / 2) / prach.LRA) / prach.SubcarrierSpacing * 1 e3;%(微秒)timingOffset = baseOffset +国防部(nSlot 10) / 10;%(微秒)其他的如图8.4.1.4.2-3 %短暂的序言,值baseOffset = 0;%(微秒)timingOffset = baseOffset +国防部(nSlot 9) / 10;%(微秒)结束sampleDelay =修复(timingOffset / 1 e6 * ofdmInfo.SampleRate);%产生传输波形(sampleDelay txwave = [0, 1);波形);%通过通道模型传递数据。附加零结束的时候%传播波形冲洗频道内容。这些零%考虑推迟引入的通道。这是一个组合%的多路延迟和实现延迟。这个值可能变化百分比取决于采样率,延迟和延迟%的传播rxwave =通道([txwave;0 (channelInfo。MaximumChannelDelay、大小(txwave, 2)));%添加噪声噪音= N0 *复杂(randn(大小(rxwave)), randn(大小(rxwave)));rxwave = rxwave +噪声;%跳过这个位置,如果是开是不活跃的。%跳过槽后推进的检测通道%确保频道总是同步的当前槽%。%如果在这个位置开环是不活跃的,接收方不应%指望开环传输,因此不应该%检测开环。跳过一个不活跃的插槽的检测%执行一致性测试时特别重要。如果%是开是不活跃的,参考波形计算内部%的| nrPRACHDetect |函数是空的。这导致了一个空%的相关性,从而检测到一个空的序言。这个空%序言导致不正确的值的检测%的概率。如果isempty (winfo.WaveformResources.PRACH.Resources.PRACHSymbols)继续;结束numActivePRACHSlots = numActivePRACHSlots + 1;%去除通道滤波器的实现延迟rxwave = rxwave ((channelInfo。ChannelFilterDelay + 1):最终,);%应用频率偏移t =((0:大小(rxwave, 1) 1) / channel.SampleRate)。”;rxwave = rxwave。* repmat (exp(1 * 2 *π* foffset * t), 1,大小(rxwave, 2));%所有细胞中开环检测指标(检测、补偿)= nrPRACHDetect (rxwave载体,开环,“DetectionThreshold”阈值);%测试中检测如果(长度(发现)= = 1)如果~ prachEnabled%的假警报测试,任何序言发现是错的detectedCount = detectedCount + 1;其他的%检查正确的序言如果(= = prach.PreambleIndex检测)%计算定时估计误差trueOffset = timingOffset / 1 e6;% (s)measuredOffset =补偿(1)/ channel.SampleRate;timingerror = abs (measuredOffset-trueOffset);%测试可接受的计时误差如果(timingerror < = timeErrorTolerance / 1 e6) detectedCount = detectedCount + 1;%检测序言结束结束结束结束结束% nSlot循环%计算最终的检测概率的信噪比pDetection (snrIdx) = detectedCount / numActivePRACHSlots;%显示这个信噪比的检测概率流(“探测概率:% d % % \ n 'pDetection (snrIdx) * 100);结束%的信噪比循环

%的阴谋探测概率图(“名字”,“探测概率”);情节(SNRdB pDetection,“这”,“线宽”2,“MarkerSize”7);标题([的探测概率的num2str (numPRACHSlots)“开槽(s)”]);包含(“信噪比(dB)”);ylabel (“探测概率”);网格在;持有在;%绘制目标概率如果prachEnabled%为错过的检测试验,检测概率应> = 99%pTarget = 99;其他的%假警报测试,检测概率应< 0.1%pTarget = 0.1;% #好< UNRCH >结束图(-6.0,pTarget / 100,“处方”,“线宽”2,“MarkerSize”7);传奇(的仿真结果,(“目标”num2str (pTarget)“%概率”),“位置”,“最佳”);minP = 0;如果(~ isnan (min (pDetection))) minP = min ([pDetection (:);pTarget]);结束轴([SNRdB (1) -0.1 SNRdB(结束)+ 0.1 minp - 0.05 - 1.05))