numPRACHSlots = 10;%在每个信噪比下模拟的PRACH槽数SNRdB = [-21， -16， -11， -6， -1];%信噪比范围，单位为dBfoffset = 400.0;%频率偏移，单位为HztimeErrorTolerance = 2.55;%时间容错，以微秒为单位阈值= [];%检测阈值prachEnabled = true;%启用PRACH传输。为了模拟虚警测试，禁用PRACH传输。

载体= nrCarrierConfig;母舰。SubcarrierSpacing = 15;母舰。NSizeGrid = 25;计算ofdm相关信息ofdmInfo = nrOFDMInfo(载体);

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

频道= nrTDLChannel;通道。DelayProfile =“TDLC300”；%延迟概要通道。MaximumDopplerShift = 100.0;%以Hz为单位的最大多普勒频移通道。SampleRate = ofdmInfo.SampleRate;%输入信号采样率，单位为Hz通道。MIMOCorrelation =“低”；% MIMO相关通道。TransmissionDirection =“上行”；%上行传输通道。NumReceiveAntennas = 2;%接收天线数通道。NormalizePathGains = true;规范化延迟配置文件功率通道。种子= 42;%通道种子。为不同的通道实现改变这一点通道。NormalizeChannelOutputs = true;接收天线归一化获取通道特征信息channelInfo =信息(渠道);

初始化每个信噪比下存储检测概率的变量pDetection = 0(大小(SNRdB));存储生成PRACH波形所需的配置参数waveconfig。NumSubframes = prach.SubframesPerPRACHSlot;waveconfig。窗口= [];waveconfig。运营商=载体;waveconfig.PRACH.Config =开环;waveconfig.PRACH.Enable = prachEnabled;临时变量'prach_init'， 'waveconfig_init'， 'ofdmInfo_init'，%和'channelInfo_init'用于创建临时变量'prach'， 'waveconfig'， 'ofdmInfo'和'channelInfo'中的信噪比循环%以创建独立的实例，以防并行模拟prach_init =开环;waveconfig_init = waveconfig;ofdmInfo_init = ofdmInfo;channelInfo_init = channelInfo;为snrIdx = 1:元素个数(SNRdB)%注释用于并行计算% parfor snrIdx = 1: numl (SNRdB) % uncomment用于并行计算要减少总模拟时间，可以执行此循环%通过使用并行计算工具箱并行。注释掉for语句语句，并取消注释'parfor'语句。如果并行计算% Toolbox(TM)未安装，“parfor”默认为普通的“for”语句在命令窗口中显示进度timeNow = char (datetime (“现在”，“格式”，“HH: mm: ss”));流([timeNow':模拟信噪比= %+5.1f dB…'), SNRdB (snrIdx));将随机数生成器设置为默认值rng (“默认”）;初始化此信噪比点的变量，初始化所需%的变量使用并行计算工具箱开环= prach_init;waveconfig = waveconfig_init;ofdmInfo = ofdmInfo_init;channelInfo = channelInfo_init;重置通道，使每个信噪比点将经历相同的%通道实现重置(渠道);%归一化噪声功率以考虑抽样率，即aOFDM调制中使用的IFFT大小的%函数。定义了信噪比每个接收天线的每个载波资源单元的%。信噪比= 10 ^ (SNRdB (snrIdx) / 10);N0 = 1 /√(2.0 * channel.NumReceiveAntennas *双(ofdmInfo.Nfft) *信噪比);%检测到的前导计数detectedCount = 0;为每个PRACH槽循环numActivePRACHSlots = 0;为nSlot = 0: numPRACHSlots-1为当前槽位生成PRACH波形开环。NPRACHSlot = nSlot;waveconfig.PRACH.Config.NPRACHSlot = nSlot;[波形,~,winfo] = hNRPRACHWaveformGenerator (waveconfig);根据TS 38.141-1图设置PRACH定时偏移(微秒)% 8.4.1.4.2-2和图8.4.1.4.2-3如果开环。上帝抵抗军= = 839%长序文，值如图8.4.1.4.2-2所示baseOffset = ((winfo.WaveformResources.PRACH.Resources.PRACHSymbolsInfo.NumCyclicShifts / 2) / prach.LRA) / prach.SubcarrierSpacing * 1 e3;%(微秒)timinggoffset = baseOffset + mod(nSlot,10)/10;%(微秒)其他的短序言，值如图8.4.1.4.2-3所示baseOffset = 0;%(微秒)timinggoffset = baseOffset + mod(nSlot,9)/10;%(微秒)结束sampleDelay = fix(timinggoffset / 1e6 * ofdmInfo.SampleRate);%生成发射波形(sampleDelay txwave = [0, 1);波形);通过通道模型传递数据。的末尾加上0%传输波形以刷新通道内容。这些零%考虑到通道中引入的任何延迟。这是一个混合%的多路径时延和实现时延。这个值可能%的变化取决于采样率，延迟配置和延迟%的传播rxwave =通道([txwave;0 (channelInfo。MaximumChannelDelay、大小(txwave, 2)));%添加噪声noise = N0*complex(randn(size(rxwave))， randn(size(rxwave)));Rxwave = Rxwave + noise;如果PRACH未激活则跳过此槽位。。前进通道到后，跳过此槽的检测。%，确保通道始终与当前槽%。如果PRACH在这个槽位是不活动的，接收器不应该预期任何PRACH传输，因此甚至不应该尝试%检测一个PRACH。跳过未激活插槽的检测是%在执行一致性测试时尤其重要。如果% PRACH是不活动的，在内部计算的参考波形% |nrPRACHDetect|函数为空。这就导致了一个空%相关性，从而到一个空的检测到的序言。这个空% preamble导致检测值不正确%的概率。如果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(1i*2*pi*foffset*t)， 1, size(Rxwave, 2);% PRACH检测所有细胞前导索引[detected, offset] = nrPRACHDetect(carrier, prach, rxwave，“DetectionThreshold”阈值);测试前导检测如果(长度(发现)= = 1)如果~ prachEnabled对于虚警测试，检测到的任何前导都是错误的。detectedCount = detectedCount + 1;其他的测试正确的前导检测如果(= = prach.PreambleIndex检测)计算定时估计误差trueOffset = timingOffset / 1 e6;% (s)measuredOffset =补偿(1)/ channel.SampleRate;timingerror = abs (measuredOffset-trueOffset);测试可接受的计时错误如果(timingerror<=timeErrorTolerance/1e6) detectedCount = detectedCount + 1;%检测序言结束结束结束结束结束%的nSlot循环计算此信噪比的最终检测概率pDetection (snrIdx) = detectedCount / numActivePRACHSlots;%显示此信噪比的检测概率流(“探测概率:% d % % \ n 'pDetection (snrIdx) * 100);结束信噪比环路的%

Plot detection概率图(“名字”，“探测概率”）;情节(SNRdB pDetection,“这”，“线宽”2,“MarkerSize”7);标题([“检测概率”num2str (numPRACHSlots)“开槽(s)”]);包含(“信噪比(dB)”）;ylabel (“探测概率”）;网格在；持有在；Plot目标概率如果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))