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

Сконфигурируйте канал для профиля задержки CDL-B от Раздела TR 38.901 7.7.1 с распространением задержки на 1 000 нс, 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 (10 МГц, 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