numSubframes=10；%在每个信噪比下要模拟的1ms子帧数SNRdB=[-21，-16，-11，-6，-1]；%信噪比范围，单位为dBfoffset = 400.0;%频率偏移量(Hz)TimeRortolerance = 2.55;microSeconds中的％误差容差

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

%设置PRACH配置Prach = nrprachconfig;Prach.frequencyRange =“FR1”;%频率范围prach.duplemmode =.'FDD';%频分双工（FDD）prach.ConfigurationIndex=27；%格式0的配置索引prach.SubcarrierSpacing=1.25；％子载波间距开环。SequenceIndex = 22;％逻辑序列索引开环。PreambleIndex = 32;%序言索引prach.RestrictedSet =.“UnrestrictedSet”;%正常模式开环。FrequencyStart = 0;%频率定位中存储的NCS表定义ZeroCorrelationZone的值% nrPRACHConfig对象开关开环。格式案例{'0'那'1'那'2'}ncsTable=nrPRACHConfig.Tables.NCSFormat012；ncsTableCol=（字符串（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(载体、开环“开窗”，开窗）；

channel = nrtdlchannel;channel.delayprofile =“TDL-C”;%延迟概要channel.delayspread = 300e-9;％延迟在几秒钟内传播channel.maximumdopplershift = 100.0;Hz中的％最大多普勒偏移通道。SampleRate = ofdmInfo.SampleRate;Hz中的％输入信号采样率信道相关性=“低的”;%MIMO相关channel.TransmissionDirection =.“上行”;％上行链路传输channel.NumTransmitAntennas=1；%发射天线的数量通道。NumReceiveAntennas = 2;%接收天线的数量channel.NormalizePathGains=true；％标准化延迟配置文件channel.seed = 42;%通道种子。为不同的通道实现更改此参数channel.normalizechanneloutpuls = true;％标准化为接收天线

%初始化存储每个SNR下检测概率的变量pDetection = 0(大小(SNRdB));％通过通道多径获得最大延迟样本数％ 成分。这是从最大的通道路径计算的%延迟和通道滤波器的实现延迟。榜样%这需要冲洗通道滤波器以获得接收信号。channelInfo=info（频道）；maxChDelay=ceil（最大值（channelInfo.PathDelays*channel.SampleRate））+channelInfo.ChannelFilterDelay；仿真期间PRACH插槽总数NumFrachSlots=地板（numSubframes/prach.SubFramesPermachSlots）；%存储生成PRACH波形所需的配置参数waveconfig.numsubframes = prach.subframesperprachslot;waveconfig.windowing =窗口;waveconfig.carriers =载体;waveconfig.prach.config = prach;％临时变量'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)%并行计算的注释对于并行计算的％pardx = 1：numel（snRDB）％取消注释％要减少总模拟时间，可以执行此循环%使用并行计算工具箱进行并行计算。注释掉“for”%语句并取消对“parfor”语句的注释。如果% Toolbox(TM)未安装，'parfor'默认为普通的'for'语句％将随机数生成器设置设置为默认值rng(“默认”）;初始化所需的SNR点的％初始化变量%使用并行计算工具箱时的变量数量prach = prach_init;waveconfig = waveconfig_init;OFDMINFO = OFDMINFO_INIT;channelinfo = channelinfo_init;％重置通道，以便每个SNR点都会遇到相同的％渠道实现重置（频道）;%标准化噪声功率以考虑采样率，这是一个OFDM调制中使用的IFFT尺寸的效果。SNR定义每个接收天线的每个资源元素%。信噪比= 10 ^ (SNRdB (snrIdx) / 10);N0 = 1 /√(2.0 * channel.NumReceiveAntennas *双(ofdmInfo.Nfft) *信噪比);%检测到的前导码计数检测到= 0;%每个PRACH插槽的循环numActivePRACHSlots = 0;为了nSlot=0:numPRACHSlots-1 prach.nprachlots=nSlot；%产生电流槽的PRACH波形waveconfig.PRACH.Config.NPRACHSlot=nSlot；[waveform，~，winfo]=HNRPachWaveOrmGenerator（waveconfig）；%如果PRACH处于非活动状态，请跳过此插槽如果（isempty（winfo.WaveformResources.PRACH））继续;结尾numActivePRACHSlots = numActivePRACHSlots + 1;根据TS 38.141-1数字设置PRACH时序偏移的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-3所示baseOffset = 0;%(微秒)定时偏移=基准偏移+模（nSlot，9）/10；%(微秒)结尾采样延迟=固定（定时偏移/1e6*ofdmInfo.SampleRate）；%产生发射波形(sampleDelay txwave = [0, 1);波形(1:(end-sampleDelay)));％通过频道模型传递数据。在末尾附加零%发送波形以刷新频道内容。这些零取零百分比频道引入的任何延迟。这是一个混合多路径延迟和实现延迟的百分比。这个值可能%根据采样率、延迟曲线和延迟进行更改%散布rxwave =通道([txwave;0 (maxChDelay、尺寸(txwave 2))));%增加噪音噪声= N0*complex(randn(size(rxwave))， randn(size(rxwave))));Rxwave = Rxwave +噪声;%消除通道建模的实现延迟rxwave=rxwave（（channelInfo.ChannelFilterDelay+1）：结束，：）；％施加频率偏移量t=（（0:size（rxwave，1）-1）/channel.SampleRate）；rxwave=rxwave.*repmat（exp（1i*2*pi*foffset*t），1，大小（rxwave，2））；所有小区前导码索引的％PRACH检测[detected, offset] = hPRACHDetect(carrier, prach, rxwave， (0:63).');序言检测的％测试如果（长度（检测到）=1）测试正确的序言检测如果（检测到==prach.PreambleIndex）%计算定时估计误差TrueOffset = TimingOffset / 1E6;%(s)测量OFFSET =偏移（1）/channel.SAMPLEDE;timingError = ABS（测量OFFSET-TRUEOFFSET）;％测试可接受定时错误如果（timingerror<=时间误差公差/1e6）检测计数=检测计数+1；%检测序言别的disp('定时错误'）;结尾别的disp('检测到序言不正确'）;结尾别的disp('检测到多个或零序言'）;结尾结尾%nSlot环路的设计%计算该信噪比的最终检测概率pdetection（snridx）=检测到的count / numactiveprachslots;结尾回路信噪比的%

hprachdetectionResults（SNRDB，NumSubframes，PDetection）;