numsubframes = 10;%在每个信噪比下要模拟的1ms子帧数SNRdB = [-21， -16， -11， -6， -1];%信噪比范围，单位为dBfoffset = 400.0;%频率偏移量(Hz)timeErrorTolerance = 2.55;%容错时间，以微秒为单位

载体= nrCarrierConfig;母舰。SubcarrierSpacing = 15;母舰。NSizeGrid = 25;

％set prach配置开环= nrPRACHConfig;开环。FrequencyRange ='fr1'；％ 频率范围开环。DuplexMode =“FDD”；％频分双工（FDD）Prach.configurationIndex = 27;%格式0的配置索引Prach.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.RestrictedSet);案件“3”ncsTable = nrPRACHConfig.Tables.NCSFormat3;ncsTableCol = (string(ncsTable.Properties.VariableNames) == prach.RestrictedSet);除此以外ncsTable = nrPRACHConfig.Tables.NCSFormatABC;ncsTableCol = contains(string(ncsTable.Properties.VariableNames)， num2str(prach.LRA));结束nc = 13;zeroCorrelationZone = ncsTable.ZeroCorrelationZone (ncsTable {:, ncsTableCol} = = nc);开环。ZeroCorrelationZone = ZeroCorrelationZone;循环移位指数%计算该PRACH配置的ofdm相关信息窗口= [];ofdmInfo = nrPRACHOFDMInfo(载体、开环'窗口',窗口);

频道= nrTDLChannel;通道。DelayProfile =“TDL-C”；%延迟概要通道。DelaySpread = 300 e-9;%延迟传播以秒为单位通道。MaximumDopplerShift = 100.0;%最大多普勒频移(Hz)通道。SampleRate = ofdmInfo.SampleRate;%输入信号采样率(Hz)channel.mimocorrelation =.“低”；% MIMO相关通道。TransmissionDirection =“上行”；%上行传输通道。NumTransmitAntennas = 1;百分比的发射天线通道。NumReceiveAntennas = 2;%接收天线数channel.normalizepathgains = true;正规化延迟轮廓功率通道。种子= 42;%通道种子。为不同的通道实现更改此参数通道。NormalizeChannelOutputs = true;接收天线正常化

％初始化在每个SNR处存储检测概率的变量pDetection = 0(大小(SNRdB));%获取信道多路径的最大延迟采样数%的组件。这是从最大的通道路径计算出来的频率延迟和通道过滤器的实现延迟。这个例子％需要它刷新信道滤波器以获得接收的信号。ChannelInfo = Info（频道）;maxchdelay = ceil（max（channicinfo.pathdelays * channel.samplevere））+ channelinfo.channelfilterdelay;%仿真周期内PRACH槽位总数numPRACHSlots = floor(numSubframes / prach.SubframesPerPRACHSlot);％存储生成PRACH波形所需的配置参数waveconfig。NumSubframes = prach.SubframesPerPRACHSlot;waveconfig。窗口=窗口;waveconfig。运营商=载体;waveconfig.PRACH.Config =开环;%临时变量'prach_init'， 'waveconfig_init'， 'ofdmInfo_init'，%和'channelInfo_init'用于创建临时变量% 'prach'， 'waveconfig'， 'ofdmInfo'和'channelInfo'在信噪比环路中在并行模拟的情况下，％要创建独立实例Prach_init = Prach;waveconfig_init = waveconfig;OFDMINFO_INIT = OFDMINFO;ChannelInfo_init = ChannelInfo;为snrIdx = 1:元素个数(SNRdB)％评论并行计算% parfor snrIdx = 1:numel(SNRdB) % uncomment用于并行计算为了减少总的模拟时间，你可以执行这个循环通过使用并行计算工具箱并行的％并行。评论'for'%语句并取消“parfor”语句的注释。如果并行计算% Toolbox(TM)未安装，'parfor'默认为普通的'for'语句%将随机数生成器设置为默认值rng (“默认”）;%初始化这个SNR点的变量，需要初始化使用并行计算工具箱时变量的百分比开环= prach_init;waveconfig = waveconfig_init;ofdmInfo = ofdmInfo_init;channelInfo = channelInfo_init;重置信道，使每个信噪比点体验相同%通道实现重置(渠道);％归一化噪声功率，以考虑采样率，这是一个OFDM调制中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;%产生电流槽的PRACH波形waveconfig.prach.config.nprachslot = nslot;[波形，〜，winfo] = hnrprachwaveformgenerator（waveconfig）;%如果PRACH未激活，跳过该槽位如果（isempty（winfo.waveformresources.prach））继续；结束numActivePRACHSlots = numActivePRACHSlots + 1;%按照TS 38.141-1设置PRACH定时偏移量(微秒)% 8.4.1.4.2-2和图8.4.1.4.2-3如果Prach.LRA == 839.%长序言，值如图8.4.1.4.2-2baseoffset =（（winfo.waveformresources.prach.resources.prachsymbolsinfo.numcyclicshifts / 2）/prach.lra）/prach.subcarrierspacing*1e3;%(微秒)timinggoffset = baseOffset + mod(nSlot,10)/10;%(微秒)其他的％短前导码，值如图8.4.1.4.2-3baseOffset = 0;%(微秒)timinggoffset = baseOffset + mod(nSlot,9)/10;%(微秒)结束sampleDelay = fix(timinggoffset / 1e6 * ofdmInfo.SampleRate); / /修改时间％生成传输波形(sampleDelay txwave = [0, 1);波形(1:(end-sampleDelay)));%通过通道模型传递数据。的末尾加0％发送波形到刷新信道内容。这些零点取%考虑信道中引入的任何延迟。这是一种混合多路径延迟和实现延迟的百分比。这个值可能根据采样率，延迟配置文件和延迟，百分比变更％ 传播rxwave =通道([txwave;0 (maxChDelay、尺寸(txwave 2))));％添加噪声噪声= N0*complex(randn(size(rxwave))， randn(size(rxwave))));Rxwave = Rxwave +噪声;％删除了频道建模的实现延迟rxwave = rxwave（（channicinfo.channelfilterdelay + 1）：结束，:);应用频率偏移t =（（0：size（rxwave，1）-1）/channel.samplege）。';rxwave = rxwave。* repmat（exp（1i * 2 * pi * foffset * t），1，尺寸（rxwave，2））;% PRACH检测所有细胞前序指数[detected, offset] = hPRACHDetect(carrier, prach, rxwave， (0:63).');%前导检测测试如果（长度（检测到）== 1）测试正确的序言检测如果（检测到== Prach.ProdbleIndex）％计算定时估计错误trueOffset = timingOffset / 1 e6;% (s)measuredOffset =补偿(1)/ channel.SampleRate;timingerror = abs (measuredOffset-trueOffset);测试可接受的定时错误如果(timingerror<=timeErrorTolerance/1e6) detectedCount = detectedCount + 1;%检测序言其他的DISP（“定时误差”）;结束其他的DISP（“检测到错误的序言”）;结束其他的DISP（'检测到多个或零序言'）;结束结束nslot循环的％％计算此SNR的最终检测概率pDetection (snrIdx) = detectedCount / numActivePRACHSlots;结束回路信噪比的%

hPRACHDetectionResults (SNRdB numSubframes pDetection);