Получите импульсный ответ фильтра пути для уровня ссылки MIMO, исчезающий канал
pathFilters = getPathFilters(channel)pathFilters = getPathFilters(channel)channel. Используйте pathFilters вместе с выходным аргументом pathGains, возвращенным объектом канала восстановить совершенную оценку канала.
Восстановите импульсный ответ канала и выполните оценку смещения синхронизации с помощью фильтров пути модели канала Кластеризованной строки задержки (CDL) с профилем задержки CDL-D от Раздела TR 38.901 7.7.1.
Задайте конфигурационную структуру канала с помощью Системного объекта nrCDLChannel. Используйте профиль задержки CDL-D, распространение задержки 10 нс и скорость UT 15 км/ч:
v = 15.0; % UT velocity in km/h fc = 4e9; % carrier frequency in Hz c = physconst('lightspeed'); % speed of light in m/s fd = (v*1000/3600)/c*fc; % UT max Doppler frequency in Hz cdl = nrCDLChannel; cdl.DelayProfile = 'CDL-D'; cdl.DelaySpread = 10e-9; cdl.CarrierFrequency = fc; cdl.MaximumDopplerShift = fd;
Сконфигурируйте массив передачи как [M N P Mg Ын] = [2 2 2 1 1], представляя 1 панель (Mg=1, Ng=1) с антенной решеткой 2 на 2 (M=2, N=2) и углы поляризации P=2. Сконфигурируйте получить антенную решетку как [M N P Mg Ын] = [1 1 2 1 1], представляя одну пару перекрестных поляризованных соразмещенных антенн.
cdl.TransmitAntennaArray.Size = [2 2 2 1 1]; cdl.ReceiveAntennaArray.Size = [1 1 2 1 1];
Создайте случайную форму волны 1 длительности подкадра с 8 антеннами.
SR = 15.36e6; T = SR * 1e-3; cdl.SampleRate = SR; cdlinfo = info(cdl); Nt = cdlinfo.NumTransmitAntennas; txWaveform = complex(randn(T,Nt),randn(T,Nt));
Передайте входную форму волны через канал.
[rxWaveform,pathGains] = cdl(txWaveform);
Получите фильтры пути, используемые в фильтрации канала.
pathFilters = getPathFilters(cdl);
Выполните оценку смещения синхронизации с помощью nrPerfectTimingEstmate.
[offset,mag] = nrPerfectTimingEstimate(pathGains,pathFilters);
Постройте значение импульсного ответа канала.
[Nh,Nr] = size(mag); plot(0:(Nh-1),mag,'o:'); hold on; plot([offset offset],[0 max(mag(:))*1.25],'k:','LineWidth',2); axis([0 Nh-1 0 max(mag(:))*1.25]); legends = "|h|, antenna " + num2cell(1:Nr); legend([legends "Timing offset estimate"]); ylabel('|h|'); xlabel('Channel impulse response samples');
[1] 3GPP TR 38.901. “Исследование модели канала для частот от 0,5 до 100 ГГц”. Проект Партнерства третьего поколения; Сеть радиодоступа Technical Specification Group.