Практическая оценка канала
[ выполняет практическую оценку канала на полученной ресурсной сетке h,nVar,info] = nrChannelEstimate(rxGrid,refInd,refSym)rxGrid при помощи ссылочной ресурсной сетки, содержащей ссылочные символы refSym в местах нахождения refInd. Функция возвращает оценку канала h, оценка отклонения шума nVarи дополнительную информацию info .
[ задает опции с помощью одного или нескольких аргументов пары "имя-значение" в дополнение к входным параметрам в любом из предыдущих синтаксисов.h,nVar,info] = nrChannelEstimate(___,Name,Value)
Сгенерировать символы опорного сигнала демодуляции (DM-RS) физического широковещательного канала (PBCH) для идентификационного номера 42 камеры физического слоя. Зависящая от времени часть инициализации скремблирования DM-RS составляет 0.
ncellid = 42; ibar_SSB = 0; dmrsSym = nrPBCHDMRS(ncellid,ibar_SSB);
Получите индексы ресурсного элемента для PBCH DM-RS.
dmrsInd = nrPBCHDMRSIndices(ncellid);
Создайте ресурсную сетку, содержащую сгенерированные символы DM-RS.
nrb = 20; scs = 15; carrier = nrCarrierConfig('NSizeGrid',nrb,'SubcarrierSpacing',scs); nTxAnts = 1; txGrid = nrResourceGrid(carrier,nTxAnts); txGrid(dmrsInd) = dmrsSym;
Модулируйте ресурсную сетку с помощью заданных длины БПФ и длины циклического префикса.
ofdmInfo = nrOFDMInfo(carrier); nulls = [1:136 377:512].'; txWaveform = nrOFDMModulate(carrier,txGrid);
Создайте модель канала TDL-C с заданными свойствами.
channel = nrTDLChannel;
channel.NumReceiveAntennas = 1;
channel.SampleRate = ofdmInfo.SampleRate;
channel.DelayProfile = 'TDL-C';
channel.DelaySpread = 100e-9;
channel.MaximumDopplerShift = 20;Получите максимальное количество задержанных выборок из пути канала с помощью наибольшей задержки и задержки реализации фильтра канала.
chInfo = info(channel); maxChDelay = ceil(max(chInfo.PathDelays*channel.SampleRate)) + chInfo.ChannelFilterDelay;
Чтобы смыть задержанные выборки из канала, добавьте нули в конце переданной формы волны, соответствующей максимальному количеству задержанных выборок и количеству передающих антенн. Передайте заполненную форму волны через модель канала TDL-C.
[rxWaveform,pathGains] = channel([txWaveform; zeros(maxChDelay,nTxAnts)]);
Оцените смещение синхронизации для передачи, используя символы DM-RS в качестве опорных символов. Модуляция OFDM опорных символов использует начальное число пазов 0.
initialSlot = 0; offset = nrTimingEstimate(carrier,rxWaveform,txGrid);
Синхронизируйте принятую форму волны в соответствии с оцененным временным смещением.
rxWaveform = rxWaveform(1+offset:end,:);
Создайте полученную ресурсную сетку, содержащую демодулированную и синхронизированную полученную форму волны.
cpFraction = 0.55;
rxGrid = nrOFDMDemodulate(carrier,rxWaveform,'CyclicPrefixFraction',cpFraction);Получите практическую оценку канала.
H = nrChannelEstimate(rxGrid,dmrsInd,dmrsSym);
Получите идеальную оценку канала.
pathFilters = getPathFilters(channel); H_ideal = nrPerfectChannelEstimate(carrier,pathGains,pathFilters,offset);
Сравните практические и совершенные оценки канала.
figure; subplot(1,2,1); imagesc(abs(H)); xlabel('OFDM Symbol'); ylabel('Subcarrier'); title('Practical Estimate Magnitude'); subplot(1,2,2); imagesc(abs(H_ideal)); xlabel('OFDM Symbol'); ylabel('Subcarrier'); title('Perfect Estimate Magnitude');
