loadFromFile = 0;% Set to 1 to load a captured waveformifloadFromFile% Load captured waveformrx = load(“capturedWaveformSIB1.mat”); rxWaveform = rx.waveform;% Configure receiver sample rate (samples/second)rxSampleRate = rx.sampleRate;% Symbol phase compensation frequency. Specify the carrier center%频率或设置为0禁用符号相位有限公司mpensationfPhaseComp = rx.fPhaseComp;% Carrier center frequency (Hz)% Set the minimum channel bandwidth for the NR band required to% configure CORESET 0 in FR1 (See TS 38.101-1 Table 5.3.5-1)minChannelBW = rx.minChannelBW;% 5, 10, 40 MHz% Configure necessary burst parameters at the receiver. The SSB pattern% can be 'Case A','Case B','Case C' for FR1 or 'Case D','Case E' for% FR2. The maximum number of blocks L_max can be 4 or 8 for FR1 and 64% for FR2.refBurst.BlockPattern = rx.ssbBlockPattern; refBurst.L_max = rx.L_max;else% Generate waveform containing SS burst and SIB1% Configure the cell identityconfig = struct(); config.NCellID = 102;% Configure an SS burstconfig.BlockPattern ='Case B';% FR1: 'Case A','Case B','Case C'. FR2: 'Case D','Case E'config.TransmittedBlocks = ones(1,8);% Bitmap of SS blocks transmittedconfig.SubcarrierSpacingCommon = 15;% SIB1 subcarrier spacing in kHz (15 or 30 for FR1. 60 or 120 for FR2)config.EnableSIB1 = 1;% Set to 0 to disable SIB1% Set the minimum channel bandwidth for the NR band required to% configure CORESET0 in FR1 (See TS 38.101-1 Table 5.3.5-1)config.MinChannelBW = 5;% 5, 10, 40 MHz% Configure and generate a waveform containing an SS burst and SIB1wavegenConfig = hSIB1WaveformConfiguration(config); [txWaveform,waveInfo] = nrWaveformGenerator(wavegenConfig); txOfdmInfo = waveInfo.ResourceGrids(1).Info;% Introduce a beamforming gain by boosting the SNR of one SSB and% associated SIB1 PDCCH and PDSCHssbIdx = 0;% Index of the SSB to boost (0-based)boost = 6;% SNR boost in dBtxWaveform = hSIB1Boost(txWaveform,wavegenConfig,waveInfo,ssbIdx,boost);% Add white Gaussian noise to the waveformrng('default');% Reset the random number generatorSNRdB = 20;% SNR for AWGNrxWaveform = awgn(txWaveform,SNRdB-boost,-10*log10(double(txOfdmInfo.Nfft)));% Configure receiver% Sample raterxSampleRate = txOfdmInfo.SampleRate;% Symbol phase compensation frequency (Hz). The function% nrWaveformGenerator does not apply symbol phase compensation to the% generated waveform.fPhaseComp = 0;% Carrier center frequency (Hz)% Minimum channel bandwidth (MHz)minChannelBW = config.MinChannelBW;% Configure necessary burst parameters at the receiverrefBurst.BlockPattern = config.BlockPattern; refBurst.L_max = numel(config.TransmittedBlocks);end% Get OFDM information from configured burst and receiver parametersnrbSSB = 20;scsSSB = hSSBurstSubcarrierSpacing(再保险fBurst.BlockPattern); rxOfdmInfo = nrOFDMInfo(nrbSSB,scsSSB,'SampleRate',rxSampleRate);% Display spectrogram of received waveformfigure; nfft = rxOfdmInfo.Nfft; spectrogram(rxWaveform(:,1),ones(nfft,1),0,nfft,'centered',rxSampleRate,'yaxis',“MinThreshold”,-130); title('Spectrogram of the Received Waveform')

disp(' -- Frequency correction and timing estimation --')% Specify the frequency offset search bandwidth in kHzsearchBW = 6*scsSSB; [rxWaveform,freqOffset,NID2] = hSSBurstFrequencyCorrect(rxWaveform,refBurst.BlockPattern,rxSampleRate,searchBW); disp([' Frequency offset: 'num2str(freqOffset,'%.0f')' Hz'])

% Create a reference grid for timing estimation using detected PSS. The PSS%放在第二个OFDM符号的引用ence grid to avoid the% special CP length of the first OFDM symbol.refGrid = zeros([nrbSSB*12 2]); refGrid(nrPSSIndices,2) = nrPSS(NID2);% Second OFDM symbol for correct CP length% Timing estimation. This is the timing offset to the OFDM symbol prior to% the detected SSB due to the content of the reference gridnSlot = 0; timingOffset = nrTimingEstimate(rxWaveform,nrbSSB,scsSSB,nSlot,refGrid,'SampleRate',rxSampleRate);% Synchronization, OFDM demodulation, and extraction of strongest SS blockrxGrid = nrOFDMDemodulate(rxWaveform(1+timingOffset:end,:),nrbSSB,scsSSB,nSlot,'SampleRate',rxSampleRate); rxGrid = rxGrid(:,2:5,:);% Display the timing offset in samples. As the symbol lengths are measured% in FFT samples, scale the symbol lengths to account for the receiver% sample rate.srRatio = rxSampleRate/(scsSSB*1e3*rxOfdmInfo.Nfft); firstSymbolLength = rxOfdmInfo.SymbolLengths(1)*srRatio; str = sprintf(' Time offset to synchronization block: %%.0f samples (%%.%.0ff ms) \n',floor(log10(rxSampleRate))-3); fprintf(str,timingOffset+firstSymbolLength,(timingOffset+firstSymbolLength)/rxSampleRate*1e3);

% Extract the received SSS symbols from the SS/PBCH blocksssIndices = nrSSSIndices; sssRx = nrExtractResources(sssIndices,rxGrid);% Correlate received SSS symbols with each possible SSS sequencesssEst = zeros(1,336);forNID1 = 0:335 ncellid = (3*NID1) + NID2; sssRef = nrSSS(ncellid); sssEst(NID1+1) = sum(abs(mean(sssRx .* conj(sssRef),1)).^2);end% Plot SSS correlationsfigure; stem(0:335,sssEst,'o'); title('SSS Correlations (Frequency Domain)'); xlabel('$N_{ID}^{(1)}$','Interpreter','latex'); ylabel('Magnitude'); axis([-1 336 0 max(sssEst)*1.1]);% Determine NID1 by finding the strongest correlationNID1 = find(sssEst==max(sssEst)) - 1;% Plot selected NID1holdon; plot(NID1,max(sssEst),'kx','LineWidth',2,'MarkerSize',8); legend(["correlations""$N_{ID}^{(1)}$ = "+ num2str(NID1)],'Interpreter','latex');% Form overall cell identity from estimated NID1 and NID2ncellid = (3*NID1) + NID2; disp([' Cell identity: 'num2str(ncellid)])

% Calculate PBCH DM-RS indicesdmrsIndices = nrPBCHDMRSIndices(ncellid);% Perform channel estimation using DM-RS symbols for each possible DM-RS% sequence and estimate the SNRdmrsEst = zeros(1,8);foribar_SSB = 0:7 refGrid = zeros([240 4]); refGrid(dmrsIndices) = nrPBCHDMRS(ncellid,ibar_SSB); [hest,nest] = nrChannelEstimate(rxGrid,refGrid,'AveragingWindow',[0 1]); dmrsEst(ibar_SSB+1) = 10*log10(mean(abs(hest(:).^2)) / nest);end% Plot PBCH DM-RS SNRsfigure; stem(0:7,dmrsEst,'o'); title('PBCH DM-RS SNR Estimates'); xlabel('$\overline{i}_{SSB}$','Interpreter','latex'); xticks(0:7); ylabel('Estimated SNR (dB)'); axis([-1 8 min(dmrsEst)-1 max(dmrsEst)+1]);% Record ibar_SSB for the highest SNRibar_SSB = find(dmrsEst==max(dmrsEst)) - 1;% Plot selected ibar_SSBholdon; plot(ibar_SSB,max(dmrsEst),'kx','LineWidth',2,'MarkerSize',8); legend(["SNRs""$\overline{i}_{SSB}$ = "+ num2str(ibar_SSB)],'Interpreter','latex');

refGrid = zeros([nrbSSB*12 4]); refGrid(dmrsIndices) = nrPBCHDMRS(ncellid,ibar_SSB); refGrid(sssIndices) = nrSSS(ncellid); [hest,nest,hestInfo] = nrChannelEstimate(rxGrid,refGrid,'AveragingWindow',[0 1]);

disp(' -- PBCH demodulation and BCH decoding -- ')% Extract the received PBCH symbols from the SS/PBCH block[pbchIndices,pbchIndicesInfo] = nrPBCHIndices(ncellid); pbchRx = nrExtractResources(pbchIndices,rxGrid);% Configure 'v' for PBCH scrambling according to TS 38.211 Section 7.3.3.1% 'v' is also the 2 LSBs of the SS/PBCH block index for L_max=4, or the 3% LSBs for L_max=8 or 64.ifrefBurst.L_max == 4 v = mod(ibar_SSB,4);elsev = ibar_SSB;endssbIndex = v;% PBCH equalization and CSI calculationpbchHest = nrExtractResources(pbchIndices,hest); [pbchEq,csi] = nrEqualizeMMSE(pbchRx,pbchHest,nest); Qm = pbchIndicesInfo.G / pbchIndicesInfo.Gd; csi = repmat(csi.',Qm,1); csi = reshape(csi,[],1);% Plot received PBCH constellation after equalizationfigure; plot(pbchEq,'o'); xlabel('In-Phase'); ylabel('Quadrature') title('Equalized PBCH Constellation'); m = max(abs([real(pbchEq(:)); imag(pbchEq(:))])) * 1.1; axis([-m m -m m]);% PBCH demodulationpbchBits = nrPBCHDecode(pbchEq,ncellid,v,nest);% Calculate RMS PBCH EVMpbchRef = nrPBCH(pbchBits<0,ncellid,v); evm = comm.EVM; pbchEVMrms = evm(pbchRef,pbchEq);% Display calculated EVMdisp([' PBCH RMS EVM: 'num2str(pbchEVMrms,'%0.3f')'%']);

% Apply CSIpbchBits = pbchBits .* csi;% Perform BCH decoding including rate recovery, polar decoding, and CRC% decoding. PBCH descrambling and separation of the BCH transport block% bits 'trblk' from 8 additional payload bits A...A+7 is also performed:% A ... A+3: 4 LSBs of system frame number% A+4: half frame number% A+5 ... A+7: for L_max=64, 3 MSBs of the SS/PBCH block index% for L_max=4 or 8, A+5 is the MSB of subcarrier offset k_SSBpolarListLength = 8; [~,crcBCH,trblk,sfn4lsb,nHalfFrame,msbidxoffset] =...nrBCHDecode(pbchBits,polarListLength,refBurst.L_max,ncellid);% Display the BCH CRCdisp([' BCH CRC: 'num2str(crcBCH)]);% Stop processing MIB and SIB1 if BCH was received with errorsifcrcBCH disp(' BCH CRC is not zero.');returnend% Use 'msbidxoffset' value to set bits of 'k_SSB' or 'ssbIndex', depending% on the number of SS/PBCH blocks in the burstif(refBurst.L_max==64) ssbIndex = ssbIndex + (bit2int(msbidxoffset,3) * 8); k_SSB = 0;elsek_SSB = msbidxoffset * 16;end% Displaying the SSB indexdisp([' SSB index: 'num2str(ssbIndex)]);

% Parse the last 23 decoded BCH transport block bits into a MIB message.% The BCH transport block 'trblk' is the RRC message BCCH-BCH-Message,% consisting of a leading 0 bit and 23 bits corresponding to the MIB. The% leading bit signals the message type transmitted (MIB or empty sequence).mib = fromBits(MIB,trblk(2:end));% Do not parse leading bit% Create set of subcarrier spacings signaled by the 7th bit of the decoded% MIB, the set is different for FR1 (L_max=4 or 8) and FR2 (L_max=64)if(refBurst.L_max==64) commonSCSs = [60 120];elsecommonSCSs = [15 30];endinitialSystemInfo =结构();initialSystemInfo。NFrame = mib.systemFrameNumber*2^4 + bit2int(sfn4lsb,4); initialSystemInfo.SubcarrierSpacingCommon = commonSCSs(mib.subCarrierSpacingCommon + 1); initialSystemInfo.k_SSB = k_SSB + mib.ssb_SubcarrierOffset; initialSystemInfo.DMRSTypeAPosition = 2 + mib.dmrs_TypeA_Position; initialSystemInfo.PDCCHConfigSIB1 = info(mib.pdcch_ConfigSIB1); initialSystemInfo.CellBarred = mib.cellBarred; initialSystemInfo.IntraFreqReselection = mib.intraFreqReselection;% Display the MIB structuredisp(' BCH/MIB Content:') disp(initialSystemInfo);% Check if a CORESET for Type0-PDCCH common search space (CSS) is present,% according to TS 38.213 Section 4.1if~isCORESET0Present(refBurst.BlockPattern,initialSystemInfo.k_SSB) fprintf('CORESET 0 is not present (k_SSB > k_SSB_max).\n');returnend

if(refBurst.L_max==64) scsKSSB = initialSystemInfo.SubcarrierSpacingCommon;elsescsKSSB = 15;endk_SSB = initialSystemInfo.k_SSB; kFreqShift = k_SSB*scsKSSB*1e3; rxWaveform = rxWaveform.*exp(1i*2*pi*kFreqShift*(0:length(rxWaveform)-1)'/rxSampleRate);% Adjust timing offset to the frame originframeOffset = hTimingOffsetToFrame(refBurst,timingOffset,ssbIndex,rxSampleRate);% If the frame offset is negative, the frame of interest is incomplete. Add% leading zeros to the waveform to align the wavefom to the frameifframeOffset < 0 rxWaveform = [zeros(-frameOffset,size(rxWaveform,2));rxWaveform];elserxWaveform = rxWaveform(1+frameOffset:end,:);end% Determine the OFDM demodulation bandwidth using CORESET 0 bandwidthcset0Idx = initialSystemInfo.PDCCHConfigSIB1.controlResourceSetZero; scsCommon = initialSystemInfo.SubcarrierSpacingCommon; scsPair = [scsSSB scsCommon]; [csetNRB,~,csetFreqOffset] = hCORESET0Resources(cset0Idx,scsPair,minChannelBW,k_SSB);% Minimum bandwidth in RB that includes CORESET 0 in received waveform.c0 = csetFreqOffset+10*scsSSB/scsCommon;% CORESET frequency offset from carrier centernrb = 2*max(c0,csetNRB-c0);% Minimum number of RB to cover CORESET 0ifrxSampleRate < nrb*12*scsCommon*1e3 disp(['SIB1 recovery cannot continue. CORESET 0 resources are beyond '...'the frequency limits of the received waveform for the sampling rate configured.']);return;end% OFDM demodulate received waveform with common subcarrier spacingnSlot = 0; rxGrid = nrOFDMDemodulate(rxWaveform, nrb, scsCommon, nSlot,...'SampleRate',rxSampleRate,'CarrierFrequency',fPhaseComp);% Display OFDM resource grid and highlight strongest SS blockfigure; imagesc(abs(rxGrid(:,:,1))); axisxyxlabel('OFDM symbol'); ylabel('Subcarrier'); numFrames = floor(length(rxWaveform)/rxSampleRate/10e-3); sfns = sprintf('(%d...%d)',initialSystemInfo.NFrame, initialSystemInfo.NFrame+numFrames-1); title(['Received Resource Grid. System Frame Number: 'sfns]); highlightSSBlock(refBurst,ssbIndex,nrb,scsPair,kFreqShift)

[csetNRB,csetDuration,csetOffset,csetPattern] = hCORESET0Resources(cset0Idx,scsPair,minChannelBW,k_SSB); ssIdx = initialSystemInfo.PDCCHConfigSIB1.searchSpaceZero; [ssSlot,ssFirstSym,isOccasion] = hPDCCH0MonitoringOccasions(ssIdx,ssbIndex,scsPair,csetPattern,csetDuration,initialSystemInfo.NFrame);% PDCCH monitoring occasions associated to different SS blocks can be in% different frames. If there are no monitoring occasions in this frame,% there must be one in the next one.slotsPerFrame = 10*scsCommon/15;if~isOccasion [ssSlot,ssFirstSym,isOccasion] = hPDCCH0MonitoringOccasions(ssIdx,ssbIndex,scsPair,csetPattern,csetDuration,initialSystemInfo.NFrame+1); ssSlot = ssSlot+slotsPerFrame;end% For FR1, UE monitors PDCCH in the Type0-PDCCH CSS over two consecutive% slots for CORESET pattern 1ifcsetPattern == 1 monSlotsPerPeriod = 2;elsemonSlotsPerPeriod = 1;end% Calculate 1-based subscripts of the subcarriers and OFDM symbols for the% slots containing the PDCCH0 associated to the detected SS block in this% and subsequent 2-frame blockscsetSubcarriers = 12*(nrb-20*scsSSB/scsCommon)/2 - csetOffset*12 + (1:csetNRB*12); numRxSym = size(rxGrid,2); symbolsPerSlot = 14; numRxSlots = ceil(numRxSym/symbolsPerSlot); monSlots = ssSlot + (0:monSlotsPerPeriod-1)' + (0:2*slotsPerFrame:(numRxSlots-ssSlot-1)); monSlots = monSlots(:)'; monSymbols = monSlots*symbolsPerSlot + (1:symbolsPerSlot)'; monSymbols = monSymbols(:)';%去除监测超过波形lim象征itsmonSymbols(monSymbols > numRxSym) = [];% Check if search space is beyond end of waveformifisempty(monSymbols) disp('Search space slot is beyond end of waveform.');return;end% Extract slots containing strongest PDCCH from the received gridrxMonSlotGrid = rxGrid(csetSubcarriers,monSymbols,:);

pdcch = hPDCCH0Configuration(ssbIndex,initialSystemInfo,scsPair,ncellid,minChannelBW);% Configure the carrier to span the BWP (CORESET 0)c0Carrier = nrCarrierConfig; c0Carrier.SubcarrierSpacing = initialSystemInfo.SubcarrierSpacingCommon; c0Carrier.NStartGrid = pdcch.NStartBWP; c0Carrier.NSizeGrid = pdcch.NSizeBWP; c0Carrier.NSlot = pdcch.SearchSpace.SlotPeriodAndOffset(2); c0Carrier.NFrame = initialSystemInfo.NFrame; c0Carrier.NCellID = ncellid;

% Specify DCI message with Format 1_0 scrambled with SI-RNTI (TS 38.212% Section 7.3.1.2.1)dci = DCIFormat1_0_SIRNTI(pdcch.NSizeBWP); disp(' -- Downlink control information message search in PDCCH -- '); siRNTI = 65535;% TS 38.321 Table 7.1-1dciCRC = true; mSlot = 0;% Loop over all monitoring slotswhile(mSlot < length(monSlots)) && dciCRC ~= 0 c0Carrier.NSlot = monSlots(mSlot+1);ifmonSlotsPerPeriod==2ifmod(mSlot,2) pdcch.SearchSpace.SlotPeriodAndOffset(2) = monSlots(2);elsepdcch.SearchSpace.SlotPeriodAndOffset(2) = monSlots(1);endend% Get PDCCH candidates according to TS 38.213 Section 10.1[pdcchInd,pdcchDmrsSym,pdcchDmrsInd] = nrPDCCHSpace(c0Carrier,pdcch); rxSlotGrid = rxMonSlotGrid(:,(1:symbolsPerSlot) + symbolsPerSlot*mSlot,:); rxSlotGrid = rxSlotGrid/max(abs(rxSlotGrid(:)));% Normalization of received RE magnitude% Loop over all supported aggregation levelsaLev = 1;while(aLev <= 5) && dciCRC ~= 0% Loop over all candidates at each aggregation level in SScIdx = 1; numCandidatesAL = pdcch.SearchSpace.NumCandidates(aLev);while(cIdx <= numCandidatesAL) && dciCRC ~= 0% Channel estimation using PDCCH DM-RS[hest,nVar,pdcchHestInfo] = nrChannelEstimate(rxSlotGrid,pdcchDmrsInd{aLev}(:,cIdx),pdcchDmrsSym{aLev}(:,cIdx));% Equalization and demodulation of PDCCH symbols[pdcchRxSym,pdcchHest] = nrExtractResources(pdcchInd{aLev}(:,cIdx),rxSlotGrid,hest); pdcchEqSym = nrEqualizeMMSE(pdcchRxSym,pdcchHest,nVar); dcicw = nrPDCCHDecode(pdcchEqSym,pdcch.DMRSScramblingID,pdcch.RNTI,nVar);% DCI message decodingpolarListLength = 8; [dcibits,dciCRC] = nrDCIDecode(dcicw,dci.Width,polarListLength,siRNTI);ifdciCRC == 0 disp([' Decoded PDCCH candidate #'num2str(cIdx)' at aggregation level 'num2str(2^(aLev-1))])endcIdx = cIdx + 1;endaLev = aLev+1;endmSlot = mSlot+1;endcIdx = cIdx-1; aLev = aLev-1; mSlot = mSlot-1; monSymbols = monSymbols(mSlot*symbolsPerSlot + (1:symbolsPerSlot));% Calculate RMS PDCCH EVMpdcchRef = nrPDCCH(double(dcicw<0),pdcch.DMRSScramblingID,pdcch.RNTI); evm = comm.EVM; pdcchEVMrms = evm(pdcchRef,pdcchEqSym);% Display calculated EVMdisp([' PDCCH RMS EVM: 'num2str(pdcchEVMrms,'%0.3f')'%']); disp([' PDCCH CRC: 'num2str(dciCRC)]);% Highlight CORESET 0/SS corresponding to strongest SSBbounding_box = @(y,x,h,w)rectangle('Position',[x+0.5 y-0.5 w h],'EdgeColor','r'); bounding_box(csetSubcarriers(1),monSymbols(1)+ssFirstSym-1,csetNRB*12,csetDuration); str = sprintf('CORESET 0/SS'); text(monSymbols(1)+ssFirstSym-7,csetSubcarriers(1)-20,0,str,'FontSize',10,'Color','w')ifdciCRC disp(' DCI decoding failed.');returnend% Plot received PDCCH constellation after equalizationfigure; plot(pdcchEqSym,'o'); xlabel('In-Phase'); ylabel('Quadrature') title('Equalized PDCCH Constellation'); m = max(abs([real(pdcchEqSym(:)); imag(pdcchEqSym(:))])) * 1.1; axis([-m m -m m]);% Display the OFDM grid of the slot containing strongest PDCCHfigure; imagesc(abs(rxSlotGrid(:,:,1))); axisxyxlabel('OFDM symbol'); ylabel('subcarrier'); title('Slot Containing Strongest PDCCH');% Highlight PDCCH in resource gridsubsPdcch = nrPDCCHSpace(c0Carrier,pdcch,'IndexStyle','Subs'); subsPdcch = double(subsPdcch{aLev}(:,:,cIdx)); x = min(subsPdcch(:,2))-1; X = max(subsPdcch(:,2))-x; y = min(subsPdcch(:,1)); Y = max(subsPdcch(:,1))-y+1; bounding_box(y,x,Y,X); str = sprintf(' PDCCH \n Aggregation Level: %d\n Candidate: %d',2.^(aLev-1),cIdx-1); text(x+X+1,y+Y/2,0,str,'FontSize',10,'Color','w')

disp(' -- PDSCH demodulation and DL-SCH decoding -- ')% Build DCI message structuredci = fromBits(dci,dcibits);% Get PDSCH configuration from cell ID, BCH information, and DCI[pdsch,K_0] = hSIB1PDSCHConfiguration(dci,pdcch.NSizeBWP,initialSystemInfo.DMRSTypeAPosition,csetPattern);% For CORESET pattern 2, the gNodeB can allocate PDSCH in the next slot,% which is indicated by the slot offset K_0 signaled by DCI. For more% information, see TS 38.214 Table 5.1.2.1.1-4.c0Carrier.NSlot = c0Carrier.NSlot+K_0; symbolOffset = symbolsPerSlot*(mSlot+K_0); monSymbols = monSymbols+symbolOffset; rxSlotGrid = rxGrid(csetSubcarriers,monSymbols,:); rxSlotGrid = rxSlotGrid/max(abs(rxSlotGrid(:)));% Normalization of received RE magnitudeifK_0 > 0% Display the OFDM grid of the slot containing associated PDSCHfigure; imagesc(abs(rxSlotGrid(:,:,1))); axisxyxlabel('OFDM symbol'); ylabel('subcarrier'); title('Slot Containing PDSCH (Slot Offset K_0 = 1)');end% PDSCH channel estimation and equalization using PDSCH DM-RSpdschDmrsIndices = nrPDSCHDMRSIndices(c0Carrier,pdsch); pdschDmrsSymbols = nrPDSCHDMRS(c0Carrier,pdsch);

mu = log2(scsCommon/15); bw = 2^mu*100;% Search bandwidth (kHz)freqStep = 2^mu;% Frequency step (kHz)freqSearch = -bw/2:freqStep:bw/2-freqStep; [~,fSearchIdx] = sort(abs(freqSearch));% Sort frequencies from centerfreqSearch = freqSearch(fSearchIdx);forfpc = fPhaseComp - kFreqShift + 1e3*freqSearch% OFDM demodulate received waveformnSlot = 0; rxGrid = nrOFDMDemodulate(rxWaveform, nrb, scsCommon, nSlot,...'SampleRate',rxSampleRate,'CarrierFrequency',fpc);% Extract monitoring slot from the received gridrxSlotGrid = rxGrid(csetSubcarriers,monSymbols,:); rxSlotGrid = rxSlotGrid/max(abs(rxSlotGrid(:)));% Normalization of received RE magnitude% Channel estimation and equalization of PDSCH symbols[hest,nVar,pdschHestInfo] = nrChannelEstimate(rxSlotGrid,pdschDmrsIndices,pdschDmrsSymbols); [pdschIndices,pdschIndicesInfo] = nrPDSCHIndices(c0Carrier,pdsch); [pdschRxSym,pdschHest] = nrExtractResources(pdschIndices,rxSlotGrid,hest); pdschEqSym = nrEqualizeMMSE(pdschRxSym,pdschHest,nVar);% PDSCH demodulationcw = nrPDSCHDecode(c0Carrier,pdsch,pdschEqSym,nVar);% Create and configure DL-SCH decoder with target code rate and% transport block sizedecodeDLSCH = nrDLSCHDecoder; Xoh_PDSCH = 0;% TS 38.214 Section 5.1.3.2tcr = hMCS(dci.ModulationCoding); NREPerPRB = pdschIndicesInfo.NREPerPRB; tbsLength = nrTBS(pdsch.Modulation,pdsch.NumLayers,length(pdsch.PRBSet),NREPerPRB,tcr,Xoh_PDSCH); decodeDLSCH.TransportBlockLength = tbsLength; decodeDLSCH.TargetCodeRate = tcr;% Decode DL-SCH[sib1bits,sib1CRC] = decodeDLSCH(cw,pdsch.Modulation,pdsch.NumLayers,dci.RedundancyVersion);ifsib1CRC == 0break;endend% Highlight PDSCH in resource gridsubsPdsch = double(nrPDSCHIndices(c0Carrier,pdsch,'IndexStyle','subscript')); x = min(subsPdsch(:,2))-1; X = max(subsPdsch(:,2))-x; y = min(subsPdsch(:,1)); Y = max(subsPdsch(:,1))-y+1; bounding_box(y,x,Y,X); str = sprintf('PDSCH (SIB1) \n Modulation: %s\n Code rate: %.2f',pdsch.Modulation,tcr); text(x+4,y+Y+60,0, str,'FontSize',10,'Color','w')% Plot received PDSCH constellation after equalizationfigure; plot(pdschEqSym,'o'); xlabel('In-Phase'); ylabel('Quadrature') title('Equalized PDSCH Constellation'); m = max(abs([real(pdschEqSym(:)); imag(pdschEqSym(:))])) * 1.1; axis([-m m -m m]);% Calculate RMS PDSCH EVM, including normalization of PDSCH symbols for any% offset between DM-RS and PDSCH powerpdschRef = nrPDSCH(c0Carrier,pdsch,double(cw{1}<0)); evm = comm.EVM; pdschEVMrms = evm(pdschRef,pdschEqSym/sqrt(var(pdschEqSym)));% Display PDSCH EVM and DL-SCH CRCdisp([' PDSCH RMS EVM: 'num2str(pdschEVMrms,'%0.3f')'%']); disp([' PDSCH CRC: 'num2str(sib1CRC)]);ifsib1CRC == 0 disp(' SIB1 decoding succeeded.');elsedisp(' SIB1 decoding failed.');end

functionpresent = isCORESET0Present(ssbBlockPattern,kSSB)switchssbBlockPatterncase{'Case A','Case B','Case C'}% FR1kssb_max = 23;case{'Case D','Case E'}% FR2kssb_max = 11;endif(kSSB <= kssb_max) present = true;elsepresent = false;endendfunctiontimingOffset = hTimingOffsetToFrame(burst,offset,ssbIdx,rxSampleRate)% As the symbol lengths are measured in FFT samples, scale the symbol% lengths to account for the receiver sample rate. Non-integer delays% are approximated at the end of the process.scs = hSSBurstSubcarrierSpacing(burst.BlockPattern); ofdmInfo = nrOFDMInfo(1,scs,'SampleRate',rxSampleRate);% smallest FFT size for SCS-SRsrRatio = rxSampleRate/(scs*1e3*ofdmInfo.Nfft); symbolLengths = ofdmInfo.SymbolLengths*srRatio;% Adjust timing offset to the start of the SS block. This step removes% the extra offset introduced in the reference grid during PSS search,% which contained the PSS in the second OFDM symbol.offset = offset + symbolLengths(1);% Timing offset is adjusted so that the received grid starts at the% frame head i.e. adjust the timing offset for the difference between% the first symbol of the strongest SSB, and the start of the frameburstStartSymbols = hSSBurstStartSymbols(burst.BlockPattern,burst.L_max);% Start symbols in SSB numerologyssbFirstSym = burstStartSymbols(ssbIdx+1);% 0-based% Adjust for whole subframessymbolsPerSubframe = length(symbolLengths); subframeOffset = floor(ssbFirstSym/symbolsPerSubframe); samplesPerSubframe = sum(symbolLengths); timingOffset = offset - (subframeOffset*samplesPerSubframe);% Adjust for remaining OFDM symbols and round offset if not integersymbolOffset = mod(ssbFirstSym,symbolsPerSubframe); timingOffset = round(timingOffset - sum(symbolLengths(1:symbolOffset)));endfunctionhighlightSSBlock(refBurst,ssbIndex,commonNRB,scs,kFreqShift) scsSSB = scs(1); scsCommon = scs(2);% Determine frequency origin of the SSB in common numerologybounding_box = @(y,x,h,w)rectangle('Position',[x+0.5 y-0.5 w h],'EdgeColor','r'); scsRatio = scsSSB/scsCommon; ssbFreqOrig = 12*(commonNRB-20*scsRatio)/2+1+kFreqShift/(scsCommon*1e3);% Determine time origin of the SSB in common numerologyssbStartSymbols = hSSBurstStartSymbols(refBurst.BlockPattern,refBurst.L_max); ssbHeadSymbol = ssbStartSymbols(ssbIndex+1)/scsRatio; ssbTailSymbol = floor((ssbStartSymbols(ssbIndex+1)+4)/scsRatio)-1; bounding_box(ssbFreqOrig,ssbHeadSymbol,240*scsRatio,ssbTailSymbol-ssbHeadSymbol+1); str = sprintf('Strongest \n SSB: %d',ssbIndex); text(ssbHeadSymbol,ssbFreqOrig-20,0, str,'FontSize',10,'Color','w')end