numSubframes = 10;在每个信噪比下模拟的1毫秒子帧数SNRdB = [-21， -16， -11， -6， -1];%信噪比范围，单位为dBFoffset = 400.0;%频率偏移，单位为HztimeErrorTolerance = 2.55;时间容错率(微秒)

carrier = nrCarrierConfig;母舰。SubcarrierSpacing = 15;母舰。NSizeGrid = 25;

设置PRACH配置prach = nrPRACHConfig;开环。FrequencyRange =“FR1”；%频率范围开环。DuplexMode =“FDD”；频分双工(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.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));结束NCS = 13;zeroCorrelationZone (ncsTable{:，ncsTableCol}==NCS);开环。ZeroCorrelationZone = ZeroCorrelationZone;%循环位移指数计算这个PRACH配置的ofdm相关信息Windowing = [];ofdmInfo = nrPRACHOFDMInfo(载体，路径，“窗口”,窗口);

channel = nrTDLChannel;通道。DelayProfile =“TDL-C”；%延迟配置文件通道。DelaySpread = 300e-9;%延迟扩散，以秒为单位通道。MaximumDopplerShift = 100.0;最大多普勒频移(以Hz为单位)通道。SampleRate = ofdmInfo.SampleRate;输入信号采样率(Hz)通道。MIMOCorrelation =“低”；% MIMO相关性通道。TransmissionDirection =“上行”；%上行传输通道。numtransmitantenna = 1;%发射天线个数通道。numreceiveantenna = 2;%接收天线个数通道。normalizepathgain = true;正常化延迟配置文件功率通道。种子= 42;%通道种子。为不同的渠道实现更改此选项通道。NormalizeChannelOutputs = true;接收天线归一化

初始化变量存储在每个信噪比下的检测概率pDetection =零(大小(SNRdB));获取通道多路径延迟采样的最大数量%的组件。这是从最大的通道路径计算出来的%延迟和信道滤波器的实现延迟。这个例子%要求刷新信道滤波器以获得接收到的信号。channelInfo =信息(通道);maxChDelay = ceil(max(channelinfo . pathdelayed *channel.SampleRate)) + channelInfo.ChannelFilterDelay;%模拟周期内PRACH插槽总数numPRACHSlots = floor(numSubframes / prach.SubframesPerPRACHSlot);存储生成PRACH波形所需的配置参数waveconfig。NumSubframes = prach.SubframesPerPRACHSlot;waveconfig。开窗=开窗;waveconfig。carrier =承运人;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)%注释用于并行计算% parfor snrIdx = 1:数值(SNRdB) % uncomment用于并行计算要减少总的模拟时间，您可以执行此循环%并行，使用并行计算工具箱。注释掉for并取消'parfor'语句的注释。如果并行计算% Toolbox(TM)未安装，“parfor”默认为普通的“for”语句将随机数生成器设置设置为默认值rng (“默认”）;初始化此信噪比点的变量，初始化所需%的变量时，使用并行计算工具箱Prach = prach_init;Waveconfig = waveconfig_init;ofdmInfo = ofdmInfo_init;channelInfo = channelInfo_init;重置信道，使每个信噪比点都有相同的体验%信道实现重置(渠道);将噪声功率归一化以考虑采样率，即a用于OFDM调制的IFFT大小的%函数。信噪比被定义每个接收天线的每个资源元素的%。信噪比= 10^(SNRdB(snrIdx)/10);N0 = 1/√(2.0*channel. numreceiveantenna *double(ofdmInfo.Nfft)*信噪比);检测到的前导数detectedCount = 0;每个PRACH槽的%循环numActivePRACHSlots = 0;为nSlot = 0:numPRACHSlots-1 prach。NPRACHSlot = nSlot;为当前槽位生成PRACH波形waveconfig.PRACH.Config.NPRACHSlot = nSlot;[波形，~，winfo] = hNRPRACHWaveformGenerator(波形配置);如果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如果开环。上帝抵抗军= = 839%长序文，取值如图8.4.1.4.2-2所示baseOffset = ((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;%(微秒)timinggoffset = baseOffset + mod(nSlot,9)/10;%(微秒)结束sampleDelay = fix(timegoffset / 1e6 * ofdmInfo.SampleRate);生成发射波形txwave = [0 (sampleDelay,1);波形(1:(end-sampleDelay)));通过通道模型传递数据。属性的末尾附加0%传输波形以刷新通道内容。这些0%，考虑通道中引入的任何延迟。这是一种混合多路径时延和实现时延的百分比。这个值可能%的变化取决于采样率，延迟配置文件和延迟%的传播Rxwave =通道([txwave;0 (maxChDelay、尺寸(txwave 2))));%添加噪音噪声= N0*complex(randn(size(rxwave))， randn(size(rxwave)));Rxwave = Rxwave +噪声;移除通道建模的实现延迟rxwave = rxwave((channelInfo. rxwave。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] = hPRACHDetect(carrier, prach, rxwave， (0:63).');测试前导检测如果(长度(发现)= = 1)测试正确的前导检测如果(= = prach.PreambleIndex检测)计算时间估计误差trueOffset = timingOffset/1e6;% (s)measuredOffset = offset (1)/channel.SampleRate;timingerror = abs(measuredooffset - trueoffset);测试可接受的计时错误如果(timingerror<=timeErrorTolerance/1e6) detectedCount = detectedCount + 1;%检测到前导其他的disp (“定时误差”）;结束其他的disp (“检测到不正确的序言”）;结束其他的disp (“检测到多个或零序文”）;结束结束%的nSlot循环计算该信噪比的最终检测概率pDetection(snrIdx) = detectedCount/numActivePRACHSlots;结束信噪比环路的%

hPRACHDetectionResults(SNRdB, numSubframes, pDetection);