Реконструировать импульсную характеристику канала и выполнять оценку временного смещения с использованием фильтров тракта lte3DChannel Системный объект.

Сконфигурируйте канал для профиля задержки CDL-B из TR 38.901, раздел 7.7.1, с 1000 нс расширением задержки, 2 передающими антеннами и 1 приемной антенной.

lte3d = lte3DChannel.makeCDL('CDL-B',1000e-9);
lte3d.Seed = 11;
lte3d.TransmitAntennaArray.Size = [1 1 2];
lte3d.ReceiveAntennaArray.Size = [1 1 1];

Создайте сигнал LTE для опорного канала измерения (RMC) R.10 (10MHz, QPSK, R = 1/3, 2 порта CRS).

rmc = lteRMCDL('R.10');
rmc.TotSubframes = 1;
data = [1; 0; 0; 1];
[txWaveform,~,txInfo] = lteRMCDLTool(rmc,data);
lte3d.SampleRate = txInfo.SamplingRate;

Пропустите сигнал через канал.

[rxWaveform,pathGains] = lte3d(txWaveform);

Оценка смещения синхронизации выполняется с использованием опорных сигналов, специфичных для соты.

corrcfg.CellRS = 'On';
offset = lteDLFrameOffset(rmc,rxWaveform,corrcfg)

offset = 15

Получите фильтры трактов, используемые при фильтрации каналов.

pathFilters = getPathFilters(lte3d);

Восстановите импульсную характеристику канала, используя фильтры тракта и коэффициенты усиления тракта. Возьмите среднее значение коэффициента усиления тракта для всех отсчетов времени (первое измерение). Создайте импульсную характеристику h для каждой передающей и приемной антенны. Суммирование откликов для каждой передающей антенны.

[~,Np,P,R] = size(pathGains);
Nh = size(pathFilters,1);
h = zeros([Nh P R]);
pathGains = permute(mean(pathGains,1),[2 3 4 1]);
for np = 1:Np
    h = h + pathFilters(:,np) .* pathGains(np,:,:);
end
h = permute(sum(h,2),[1 3 2]);
mag = abs(h);

Постройте график величины импульсной характеристики канала.

plot(mag,'o:')
title('Magnitude of Channel Impulse Response')
xlabel('Samples')
ylabel('Magnitude')

Оцените смещение по времени, найдя пик величины импульсной характеристики.

offset_ref = find(mag==max(mag)) - 1

offset_ref = 15