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

Сконфигурируйте канал для профиля задержки CDL-B из раздела 7.7.1 TR 38.901 с 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