В этом примере показано, как использовать луч, развертывающийся и в передатчике (gNB) и в приемнике (UE) концы 5G система NR. Используя блоки сигнала синхронизации (SSB), этот пример иллюстрирует некоторые процедуры управления луча, используемые во время начального доступа. Чтобы выполнить развертку луча, пример использует несколько компонентов от Phased Array System Toolbox™.
Поддержка волны миллиметра (mmWave) частоты требует направленных ссылок, которые привели к спецификации процедур управления луча по начальному доступу в NR. Управление лучом является набором (физического) Слоя 1 и Слоя 2 (среднее управление доступом) процедуры, чтобы получить и обеспечить набор ссылок пары луча (луч, используемый в gNB, соединенном с лучом, используемым в UE). Процедуры управления луча применяются и для нисходящей и для восходящей передачи и приема [1], [2]. Эти процедуры включают:
Развертка луча
Излучите измерение
Излучите определение
Излучите создание отчетов
Излучите восстановление
Этот пример фокусируется на первоначальных процедурах доступа для неактивных пользователей, когда связь устанавливается между оборудованием пользователя (UE) и узлом сети доступа (gNB). На физическом уровне, с помощью блоков сигнала синхронизации (SSB), переданных как пакет в нисходящем направлении (gNB к UE), пример подсвечивает, и передайте/получите точку (TRP) развертка луча и луч UE, развертывающийся, чтобы установить ссылку пары луча. Среди нескольких процедур управления луча TR 38.802 задает эту развертку двойного конца как процедуру p-1 [1].
После того, как соединенный, та же ссылка пары луча может использоваться для последующих передач. При необходимости лучи далее усовершенствованы с помощью CSI-RS (для нисходящего канала) и SRS (для восходящего канала). В случае отказа луча могут быть восстановлены эти парные ссылки. Для примера улучшения пары луча смотрите, что Улучшение Луча Конца Передачи Нисходящего канала NR Использует CSI-RS (5G Toolbox).
Этот пример генерирует пакет сигнала синхронизации NR, beamforms каждый из SSBs в пакете, чтобы развернуться и по азимуту и по направлениям вертикального изменения, передает этот сигнал beamformed по пространственному каналу рассеивания и обрабатывает этот полученный сигнал по нескольким лучам получать-конца. Пример измеряется, мощность приемника опорного сигнала (RSRP) для каждой передачи - получают пары луча (в двойном цикле), и определяет ссылку пары луча с максимальным RSRP. Эта ссылка пары луча таким образом показывает лучшее парное лучом в передаче, и получите концы для симулированного пространственного сценария. Этот рисунок показывает основные шаги обработки с единицами управления лучом, подсвеченными в цвете.
rng(211); % Set RNG state for repeatability
Задайте системные параметры для примера. Измените эти параметры, чтобы исследовать их удар на систему.
prm.NCellID = 1; % Cell ID prm.FreqRange = 'FR1'; % Frequency range: 'FR1' or 'FR2' prm.CenterFreq = 3.5e9; % Hz prm.SSBlockPattern = 'Case B'; % Case A/B/C/D/E prm.SSBTransmitted = [ones(1,8) zeros(1,0)]; % 4/8 or 64 in length prm.TxArraySize = [8 8]; % Transmit array size, [rows cols] prm.TxAZlim = [-60 60]; % Transmit azimuthal sweep limits prm.TxELlim = [-90 0]; % Transmit elevation sweep limits prm.RxArraySize = [2 2]; % Receive array size, [rows cols] prm.RxAZlim = [-180 180]; % Receive azimuthal sweep limits prm.RxELlim = [0 90]; % Receive elevation sweep limits prm.ElevationSweep = false; % Enable/disable elevation sweep prm.SNRdB = 30; % SNR, dB prm.RSRPMode = 'SSSwDMRS'; % {'SSSwDMRS', 'SSSonly'}
Пример использует эти параметры:
ID ячейки для сценария отдельной ячейки с одним BS и UE
Частотный диапазон, как строка, чтобы определять операцию FR1 или FR2
Центральная частота, в Гц, зависящем от частотного диапазона
Состав блока сигнала синхронизации как один из Случая A/B/C для FR1 и Случая D/E для FR2. Это также выбирает интервал поднесущей.
Переданный SSBs в шаблоне, как бинарный вектор из длины 4 или 8 для FR1 и длины 64 для FR2. Количество переданного SSBs определяет номер лучей и в передаче, и получите концы.
Передайте размер массивов, как двухэлементный вектор-строка, задающий количество антенных элементов в строках и столбцах массива передачи, соответственно. Универсальный прямоугольный массив (URA) используется, когда оба значения больше того.
Передайте азимутальные пределы развертки в градусах, чтобы задать запуск и конечные углы азимута для развертки
Передайте пределы развертки вертикального изменения в градусах, чтобы задать запуск и окончание углов возвышения для развертки
Получите размер массивов, как двухэлементный вектор-строка, задающий количество антенных элементов в строках и столбцах получить массива, соответственно. Универсальный прямоугольный массив (URA) используется, когда оба значения больше того.
Получите азимутальные пределы развертки в градусах, чтобы задать запуск и конечные углы азимута для развертки
Получите пределы развертки вертикального изменения в градусах, чтобы задать запуск и окончание углов возвышения для развертки
Включите или отключите развертку вертикального изменения и для передачи и получите концы. Включите развертку вертикального изменения для FR2 и/или URAs
Отношение сигнал-шум в дБ
Режим Measurement для SSB, чтобы задать использование только вторичных сигналов синхронизации ('SSSonly') или использование PBCH DM-RS наряду со вторичными сигналами синхронизации ('SSSwDMRS')
prm = validateParams(prm);
Настройте пакетные параметры сигнала синхронизации при помощи заданных системных параметров. Для начального доступа, набор периодичность SSB к 20 мс.
txBurst.BlockPattern = prm.SSBlockPattern; txBurst.SSBTransmitted = prm.SSBTransmitted; txBurst.NCellID = prm.NCellID; txBurst.SSBPeriodicity = 20; txBurst.NFrame = 0; txBurst.Windowing = 0; txBurst.DisplayBurst = true; % Assume same subcarrier spacing for carrier as the burst carrier = nrCarrierConfig('NCellID',prm.NCellID,'NFrame',txBurst.NFrame); carrier.SubcarrierSpacing = prm.SCS; carrierInfo = nrOFDMInfo(carrier); txBurst.SampleRate = carrierInfo.SampleRate;
Обратитесь к учебным Блокам Сигнала Синхронизации и Пакетам (5G Toolbox) для получения дополнительной информации о блоках сигнала синхронизации и пакетах.
Сконфигурируйте пространственное рассеивание канал MIMO channel. Эта модель канала применяет потери при распространении в свободном пространстве и, опционально, другие атмосферные затухания к входу. Задайте местоположения для BS и UE как [x,y,z] координаты в декартовой системе. В зависимости от заданных размеров массивов используйте или универсальные линейные матрицы (ULA) или универсальные прямоугольные массивы (URA). Используйте изотропные антенные элементы для массивов.
c = physconst('LightSpeed'); % Propagation speed lambda = c/prm.CenterFreq; % Wavelength prm.posTx = [0;0;0]; % Transmit array position, [x;y;z], meters prm.posRx = [100;50;0]; % Receive array position, [x;y;z], meters toRxRange = rangeangle(prm.posTx,prm.posRx); spLoss = fspl(toRxRange,lambda); % Free space path loss % Transmit array if prm.IsTxURA % Uniform rectangular array arrayTx = phased.URA(prm.TxArraySize,0.5*lambda, ... 'Element',phased.IsotropicAntennaElement('BackBaffled',true)); else % Uniform linear array arrayTx = phased.ULA(prm.NumTx, ... 'ElementSpacing',0.5*lambda, ... 'Element',phased.IsotropicAntennaElement('BackBaffled',true)); end % Receive array if prm.IsRxURA % Uniform rectangular array arrayRx = phased.URA(prm.RxArraySize,0.5*lambda, ... 'Element',phased.IsotropicAntennaElement); else % Uniform linear array arrayRx = phased.ULA(prm.NumRx, ... 'ElementSpacing',0.5*lambda, ... 'Element',phased.IsotropicAntennaElement); end % Scatterer locations prm.FixedScatMode = true; if prm.FixedScatMode % Fixed single scatterer location prm.ScatPos = [50; 80; 0]; else % Generate scatterers at random positions Nscat = 10; % Number of scatterers azRange = -180:180; elRange = -90:90; randAzOrder = randperm(length(azRange)); randElOrder = randperm(length(elRange)); azAngInSph = azRange(randAzOrder(1:Nscat)); elAngInSph = elRange(randElOrder(1:Nscat)); r = 20; % radius [x,y,z] = sph2cart(deg2rad(azAngInSph),deg2rad(elAngInSph),r); prm.ScatPos = [x;y;z] + (prm.posTx + prm.posRx)/2; end % Configure channel channel = phased.ScatteringMIMOChannel; channel.PropagationSpeed = c; channel.CarrierFrequency = prm.CenterFreq; channel.SampleRate = txBurst.SampleRate; channel.SimulateDirectPath = false; channel.ChannelResponseOutputPort = true; channel.Polarization = 'None'; channel.TransmitArray = arrayTx; channel.TransmitArrayPosition = prm.posTx; channel.ReceiveArray = arrayRx; channel.ReceiveArrayPosition = prm.posRx; channel.ScattererSpecificationSource = 'Property'; channel.ScattererPosition = prm.ScatPos; channel.ScattererCoefficient = ones(1,size(prm.ScatPos,2)); % Get maximum channel delay [~,~,tau] = channel(complex(randn(txBurst.SampleRate*1e-3,prm.NumTx), ... randn(txBurst.SampleRate*1e-3,prm.NumTx))); maxChDelay = ceil(max(tau)*txBurst.SampleRate);
Создайте пакетную форму волны SS [3] путем вызова hSSBurst функция помощника. Сгенерированная форма волны еще не beamformed.
% Create and display burst information txBurstInfo = hSSBurstInfo(txBurst); disp(txBurstInfo); % Generate burst waveform and grid [burstWaveform,txBurstGrid] = hSSBurst(txBurst);
SubcarrierSpacing: 30
NCRB_SSB: -20
k_SSB: 0
FrequencyOffsetSSB: 0
MIB: [24x1 double]
L: 8
SSBIndex: [0 1 2 3 4 5 6 7]
i_SSB: [0 1 2 3 4 5 6 7]
ibar_SSB: [0 1 2 3 4 5 6 7]
SampleRate: 30720000
Nfft: 1024
NRB: 72
CyclicPrefix: 'Normal'
OccupiedSubcarriers: [240x1 double]
OccupiedSymbols: [8x4 double]
Windowing: 0
Чтобы достигнуть развертки луча TRP, beamform, каждый из SS блокируется в сгенерированном пакете с помощью аналога beamforming. На основе количества блоков SS в пакете и указанных диапазонах развертки, определите и азимут и направления вертикального изменения для различных лучей. Затем beamform индивидуум блокируется в пакете к каждому из этих направлений.
% Number of beams at both transmit and receive ends numBeams = sum(txBurst.SSBTransmitted); % Transmit beam angles in azimuth and elevation, equi-spaced azBW = beamwidth(arrayTx,prm.CenterFreq,'Cut','Azimuth'); elBW = beamwidth(arrayTx,prm.CenterFreq,'Cut','Elevation'); txBeamAng = hGetBeamSweepAngles(numBeams,prm.TxAZlim,prm.TxELlim, ... azBW,elBW,prm.ElevationSweep); % For evaluating transmit-side steering weights SteerVecTx = phased.SteeringVector('SensorArray',arrayTx, ... 'PropagationSpeed',c); % Apply steering per OFDM symbol for each SSB carrier.NSizeGrid = txBurstInfo.NRB; ofdmInfo = nrOFDMInfo(carrier); gridSymLengths = repmat(ofdmInfo.SymbolLengths,1, ... size(txBurstGrid,2)/length(ofdmInfo.SymbolLengths)); % repeat burst over numTx to prepare for steering strTxWaveform = repmat(burstWaveform,1,prm.NumTx)./sqrt(prm.NumTx); for ssb = 1:length(txBurstInfo.SSBIndex) % Extract SSB waveform from burst blockSymbols = txBurstInfo.OccupiedSymbols(ssb,:); startSSBInd = sum(gridSymLengths(1:blockSymbols(1)-1))+1; endSSBInd = sum(gridSymLengths(1:blockSymbols(4))); ssbWaveform = strTxWaveform(startSSBInd:endSSBInd,1); % Generate weights for steered direction wT = SteerVecTx(prm.CenterFreq,txBeamAng(:,ssb)); % Apply weights per transmit element to SSB strTxWaveform(startSSBInd:endSSBInd,:) = ssbWaveform.*(wT'); end
beamformed разрываются, форма волны затем передается по пространственно осведомленному каналу рассеивания.
Для развертки луча получать-конца разрываются переданные beamformed, форма волны получена последовательно по каждому, получают луч. Для N передайте лучи и M получите лучи в процедуре p-1, каждом N лучами является переданный M времена от gNB так, чтобы каждый луч передачи был получен по M получите лучи.
Пример принимает оба N и M быть равным количеству SSBs в пакете. Для простоты пример генерирует только один пакет, но подражать пакетному приему по воздуху M времена, приемник обрабатывает этот один пакет M \times.
Этот рисунок показывает основанную на луче схему для разверток и в gNB и в UE для N = M = 4, в азимутальной плоскости. Схема показывает время, потраченное для двойной развертки, где каждый интервал в gNB соответствует SSB, и каждый интервал в UE соответствует пакету SS. Для изображенного сценария, S3 лучей и U2 подсвечены как выбранная парная лучом ссылка умозрительно. Пример реализует двойную развертку по длительности времени N*M моменты времени.
Получить обработка переданного пакета включает
Приложение пространственно осведомленного исчезающего канала
Получите усиление, чтобы компенсировать вызванную потерю пути и AWGN
Получать-конец beamforming
Синхронизация коррекции
Демодуляция OFDM
Извлечение известной сетки SSB
Измерение RSRP на основе заданного режима измерения
Обработка повторяет эти шаги для каждого из получить лучей, затем выбирает лучшее парное лучом на основе полного набора сделанных измерений.
Чтобы подсветить развертку луча, пример принимает известную информацию SSB в приемнике. Для получения дополнительной информации о восстановлении обработка смотрите Поиск Ячейки NR и MIB и Восстановление SIB1 (5G Toolbox).
Для нерабочего режима измерение SS-RSRP используйте или только вторичные сигналы синхронизации (SSS) или физический канал телевизионного вещания (PBCH) опорные сигналы демодуляции (DM-RS) в дополнение к SSS (Раздел 5.1.1. из [4]). Задайте это RSRPMode параметр примера. Для FR2 измерение RSRP основано на объединенном сигнале от антенных элементов, в то время как измерение на антенный элемент для FR1.
% Receive beam angles in azimuth and elevation, equi-spaced azBW = beamwidth(arrayRx,prm.CenterFreq,'Cut','Azimuth'); elBW = beamwidth(arrayRx,prm.CenterFreq,'Cut','Elevation'); rxBeamAng = hGetBeamSweepAngles(numBeams,prm.RxAZlim,prm.RxELlim, ... azBW,elBW,prm.ElevationSweep); % For evaluating receive-side steering weights SteerVecRx = phased.SteeringVector('SensorArray',arrayRx, ... 'PropagationSpeed',c); % AWGN level SNR = 10^(prm.SNRdB/20); % Convert to linear gain N0 = 1/(sqrt(2.0*prm.NumRx*double(ofdmInfo.Nfft))*SNR); % Noise Std. Dev. % Receive gain in linear terms, to compensate for the path loss rxGain = 10^(spLoss/20); % Generate a reference grid for timing correction % assumes an SSB in first slot pssRef = nrPSS(carrier.NCellID); pssInd = nrPSSIndices; pbchdmrsRef = nrPBCHDMRS(carrier.NCellID,txBurstInfo.ibar_SSB(1)); pbchDMRSInd = nrPBCHDMRSIndices(carrier.NCellID); pssGrid = zeros([240 4]); pssGrid(pssInd) = pssRef; pssGrid(pbchDMRSInd) = pbchdmrsRef; refGrid = zeros([12*carrier.NSizeGrid ofdmInfo.SymbolsPerSlot]); refGrid(txBurstInfo.OccupiedSubcarriers, ... txBurstInfo.OccupiedSymbols(1,:)) = pssGrid; % Loop over all receive beams rsrp = zeros(numBeams,numBeams); for rIdx = 1:numBeams % Fading channel, with path loss txWave = [strTxWaveform; zeros(maxChDelay,size(strTxWaveform,2))]; fadWave = channel(txWave); % Receive gain, to compensate for the path loss fadWaveG = fadWave*rxGain; % Add WGN noise = N0*complex(randn(size(fadWaveG)),randn(size(fadWaveG))); rxWaveform = fadWaveG + noise; % Generate weights for steered direction wR = SteerVecRx(prm.CenterFreq,rxBeamAng(:,rIdx)); % Apply weights per receive element if strcmp(prm.FreqRange, 'FR1') strRxWaveform = rxWaveform.*(wR'); else % for FR2, combine signal from antenna elements strRxWaveform = rxWaveform*conj(wR); end % Correct timing offset = nrTimingEstimate(carrier, ... strRxWaveform(1:ofdmInfo.SampleRate*1e-3,:),refGrid*wR(1)'); if offset > maxChDelay offset = 0; end strRxWaveformS = strRxWaveform(1+offset:end,:); % OFDM Demodulate rxGrid = nrOFDMDemodulate(carrier,strRxWaveformS); % Loop over all SSBs in rxGrid (transmit end) for tIdx = 1:numBeams % Get each SSB grid rxSSBGrid = rxGrid(txBurstInfo.OccupiedSubcarriers, ... txBurstInfo.OccupiedSymbols(tIdx,:),:); % Make measurements, store per receive, transmit beam rsrp(rIdx,tIdx) = measureSSB(rxSSBGrid,prm.RSRPMode,txBurst.NCellID); end end
После того, как развертка двойного конца и измерения завершены, определяют лучшую парную лучом ссылку на основе измерения RSRP.
[m,i] = max(rsrp,[],'all','linear'); % First occurence is output % i is column-down first (for receive), then across columns (for transmit) [rxBeamID,txBeamID] = ind2sub([numBeams numBeams],i(1)); % Display the selected beam pair disp(['Selected Beam pair with RSRP: ' num2str(10*log10(rsrp(rxBeamID, ... txBeamID))+30) ' dBm', 13 ' Transmit #' num2str(txBeamID) ... ' (Azimuth: ' num2str(txBeamAng(1,txBeamID)) ', Elevation: ' ... num2str(txBeamAng(2,txBeamID)) ')' 13 ' Receive #' num2str(rxBeamID) ... ' (Azimuth: ' num2str(rxBeamAng(1,rxBeamID)) ', Elevation: ' ... num2str(rxBeamAng(2,rxBeamID)) ')' ]); % Display final beam pair patterns h = figure('Position',figposition([32 55 32 40]),'MenuBar','none'); h.Name = 'Selected Transmit Array Response Pattern'; wT = SteerVecTx(prm.CenterFreq,txBeamAng(:,txBeamID)); pattern(arrayTx,prm.CenterFreq,'PropagationSpeed',c,'Weights',wT); h = figure('Position',figposition([32 55 32 40]),'MenuBar','none'); h.Name = 'Selected Receive Array Response Pattern'; wR = SteerVecRx(prm.CenterFreq,rxBeamAng(:,rxBeamID)); pattern(arrayRx,prm.CenterFreq,'PropagationSpeed',c,'Weights',wR); % Plot MIMO scenario with tx, rx, scatterers, and determined beams prmScene = struct(); prmScene.txArraySize = prm.TxArraySize; prmScene.rxArraySize = prm.RxArraySize; prmScene.txElemPos = getElementPosition(arrayTx); % meters prmScene.rxElemPos = getElementPosition(arrayRx); % meters prmScene.txArrayPos = prm.posTx; prmScene.rxArrayPos = prm.posRx; prmScene.txAzAngles = -90:90; prmScene.rxAzAngles = [90:180 -179:-90]; prmScene.scatPos = prm.ScatPos; prmScene.lambda = lambda; prmScene.arrayScaling = 1; hPlotSpatialMIMOScene(prmScene,wT,wR); if ~prm.ElevationSweep view(2); end
Selected Beam pair with RSRP: 45.1516 dBm Transmit #8 (Azimuth: 60, Elevation: 0) Receive #6 (Azimuth: 90, Elevation: 0)
Эти графики подсвечивают шаблон направленности передачи, получают шаблон направленности и пространственную сцену, соответственно. Результаты зависят от отдельных направлений луча, используемых для разверток. Пространственная сцена открывает объединенный вид передачи, и получите массивы и соответствующие решительные лучи, наряду с рассеивателями.
Этот пример подсвечивает, что P-1 излучают процедуру управления при помощи блоков сигнала синхронизации для конца передачи и развертки луча получать-конца. Путем измерения мощности приемника опорного сигнала для SSBs можно идентифицировать лучшую ссылку пары луча для выбранной пространственной среды.
Пример позволяет изменение на частотном диапазоне, составе блока SSB, количестве SSBs, передачи, и получите размеры массивов, передайте и получите области значений развертки и измеряющийся режим. Чтобы видеть удар параметров на выборе луча, экспериментируйте с различными значениями. Получить обработка упрощена, чтобы подсветить beamforming аспекты для примера.
Для примера P-2 процедур развертки луча конца передачи с помощью сигналов CSI-RS для нисходящего канала смотрите, что Улучшение Луча Конца Передачи Нисходящего канала NR Использует CSI-RS (5G Toolbox). Можно использовать эти процедуры для улучшения луча и корректировки в связанном режиме, когда-то начальные ссылки пары луча устанавливаются [5], [6].
3GPP TR 38.802. "Исследование Новых Радио-технологических аспектов физического уровня доступа". Проект Партнерства третьего поколения; Сеть радиодоступа Technical Specification Group.
Джордани, M., М. Полезе, А. Рой, Д. Кэстор и М. Зорзи. "Пример на управлении лучом для 3GPP NR на mmWave частотах". Коммуникация IEEE Surveys & Tutorials, издание 21, № 1, 1 квартал 2019.
3GPP TS 38.211. "NR; Физические каналы и модуляция". Проект Партнерства третьего поколения; Сеть радиодоступа Technical Specification Group.
3GPP TS 38.215. "NR; измерения Физического уровня". Проект Партнерства третьего поколения; Сеть радиодоступа Technical Specification Group.
Джордани, M., М. Полезе, А. Рой, Д. Кэстор и М. Зорзи. "Автономное и неавтономное управление лучом для 3GPP NR в mmWaves". Коммуникация IEEE Мэг., апрель 2019, стр 123-129.
Onggosanusi, E. S. Md. Рахман, и др. "Модульный и канал с высоким разрешением утверждает информацию и управление лучом для 5G NR". Коммуникация IEEE Мэг., март 2018, стр 48-55.
function prm = validateParams(prm) % Validate user specified parameters and return updated parameters % % Only cross-dependent checks are made for parameter consistency. if strcmpi(prm.FreqRange,'FR1') if prm.CenterFreq > 7.125e9 || prm.CenterFreq < 410e6 error(['Specified center frequency is outside the FR1 ', ... 'frequency range (410 MHz - 7.125 GHz).']); end if strcmpi(prm.SSBlockPattern,'Case D') || ... strcmpi(prm.SSBlockPattern,'Case E') error(['Invalid SSBlockPattern for selected FR1 frequency ' ... 'range. SSBlockPattern must be one of ''Case A'' or ' ... '''Case B'' or ''Case C'' for FR1.']); end if ~((length(prm.SSBTransmitted)==4) || ... (length(prm.SSBTransmitted)==8)) error(['SSBTransmitted must be a vector of length 4 or 8', ... 'for FR1 frequency range.']); end if (prm.CenterFreq <= 3e9) && (length(prm.SSBTransmitted)~=4) error(['SSBTransmitted must be a vector of length 4 for ' ... 'center frequency less than or equal to 3GHz.']); end if (prm.CenterFreq > 3e9) && (length(prm.SSBTransmitted)~=8) error(['SSBTransmitted must be a vector of length 8 for ', ... 'center frequency greater than 3GHz and less than ', ... 'or equal to 7.125GHz.']); end else % 'FR2' if prm.CenterFreq > 52.6e9 || prm.CenterFreq < 24.25e9 error(['Specified center frequency is outside the FR2 ', ... 'frequency range (24.25 GHz - 52.6 GHz).']); end if ~(strcmpi(prm.SSBlockPattern,'Case D') || ... strcmpi(prm.SSBlockPattern,'Case E')) error(['Invalid SSBlockPattern for selected FR2 frequency ' ... 'range. SSBlockPattern must be either ''Case D'' or ' ... '''Case E'' for FR2.']); end if length(prm.SSBTransmitted)~=64 error(['SSBTransmitted must be a vector of length 64 for ', ... 'FR2 frequency range.']); end end prm.NumTx = prod(prm.TxArraySize); prm.NumRx = prod(prm.RxArraySize); if prm.NumTx==1 || prm.NumRx==1 error(['Number of transmit or receive antenna elements must be', ... ' greater than 1.']); end prm.IsTxURA = (prm.TxArraySize(1)>1) && (prm.TxArraySize(2)>1); prm.IsRxURA = (prm.RxArraySize(1)>1) && (prm.RxArraySize(2)>1); if ~( strcmpi(prm.RSRPMode,'SSSonly') || ... strcmpi(prm.RSRPMode,'SSSwDMRS') ) error(['Invalid RSRP measuring mode. Specify either ', ... '''SSSonly'' or ''SSSwDMRS'' as the mode.']); end % Select SCS based on SSBlockPattern switch lower(prm.SSBlockPattern) case 'case a' scs = 15; case {'case b', 'case c'} scs = 30; case 'case d' scs = 120; case 'case e' scs = 240; end prm.SCS = scs; end function rsrp = measureSSB(rxSSBGrid,mode,NCellID) % Compute the reference signal received power (RSRP) based on SSS, and if % selected, also PBCH DM-RS. sssInd = nrSSSIndices; % SSS indices numRx = size(rxSSBGrid,3); rsrpSSS = zeros(numRx,1); for rxIdx = 1:numRx % Extract signals per rx element rxSSBGridperRx = rxSSBGrid(:,:,rxIdx); rxSSS = rxSSBGridperRx(sssInd); % Average power contributions over all REs for RS rsrpSSS(rxIdx) = mean(rxSSS.*conj(rxSSS)); end if strcmpi(mode,'SSSwDMRS') pbchDMRSInd = nrPBCHDMRSIndices(NCellID); % PBCH DM-RS indices rsrpDMRS = zeros(numRx,1); for rxIdx = 1:numRx % Extract signals per rx element rxSSBGridperRx = rxSSBGrid(:,:,rxIdx); rxPBCHDMRS = rxSSBGridperRx(pbchDMRSInd); % Average power contributions over all REs for RS rsrpDMRS(rxIdx) = mean(rxPBCHDMRS.*conj(rxPBCHDMRS)); end end switch lower(mode) case 'sssonly' % Only SSS rsrp = max(rsrpSSS); % max over receive elements case 'ssswdmrs' % Both SSS and PBCH-DMRS, accounting for REs per RS rsrp = max((rsrpSSS*127+rsrpDMRS*144)/271); % max over receive elements end end