Импорт и визуализация данных созданий

Импорт файла OpenStreetMap (.osm), соответствующего Canary Wharf в Лондоне, Великобритания. Файл был загружен из https://www.openstreetmap.org, что обеспечивает доступ к картографическим данным, полученным по всему миру. Данные лицензированы под лицензией Open Data Commons Open Database License (ODbL), https://opendatacommons.org/licenses/odbl/. Информация о созданиях, содержащаяся в файле OpenStreetMap, импортируется и визуализируется в Средство просмотра сайта.

viewer = siteviewer("Buildings","canarywharf.osm","Basemap","topographic");

Определение площадки передатчика

Задайте сайт передатчика, чтобы смоделировать сценарий небольшой камеры в плотном городском окружении. Узел передатчика представляет базовую станцию, которая размещена на шесте, обслуживающем окружающую зону, которая включает в себя соседний парк. Передатчик использует изотропную антенну по умолчанию и работает на несущей частоте 28 ГГц с уровнем степени 5 Вт.

tx = txsite("Name","Small cell transmitter", ...
    "Latitude",51.50375, ...
    "Longitude",-0.01843, ...
    "AntennaHeight",10, ...
    "TransmitterPower",5, ...
    "TransmitterFrequency",28e9);
show(tx)

Просмотр карты покрытия для распространения линии видимости

Создайте модель распространения трассировки лучей с помощью метода съёмки и прыгающего луча (SBR). Модель распространения SBR использует анализ трассировки лучей для вычисления путей распространения и их соответствующих потерь пути. Потери пути вычисляются из потерь свободного пространства, потерь отражения из-за материала и потерь поляризации антенны.

Установите максимальное количество отражений равным 0 в порядок, чтобы ограничить начальный анализ только путями распространения линии визирования. Установите типы материалов создания и местности, чтобы смоделировать идеальное отражение.

rtpm = propagationModel("raytracing", ...
    "Method","sbr", ...
    "MaxNumReflections",0, ...
    "BuildingsMaterial","perfect-reflector", ...
    "TerrainMaterial","perfect-reflector");

Просмотрите соответствующую карту покрытия для максимальной области значений 250 метров от базовой станции. Карта покрытия показывает мощность приемника для приемника в каждом месте земли, но не вычисляется для верхних частей или сторон создания.

coverage(tx,rtpm, ...
    "SignalStrengths",-120:-5, ...
    "MaxRange",250, ...
    "Resolution",3, ...
    "Transparency",0.6)

Определите месторасположение приемника в местоположении, не связанном с линией видимости

Карта покрытия для распространения линии визирования показывает затенение из-за препятствий. Определите сайт приемника, чтобы смоделировать мобильный приемник в затрудненном месте. Постройте график линии видимости, чтобы показать препятствующий путь от передатчика к приемнику.

rx = rxsite("Name","Small cell receiver", ...
    "Latitude",51.50216, ...
    "Longitude",-0.01769, ...
    "AntennaHeight",1);

los(tx,rx)

Постройте график пути распространения с помощью трассировки лучей

Скорректируйте модель распространения трассировки лучей, включив в нее пути с одним отражением, и постройте график лучей. Результат показывает распространение сигнала вдоль пути с одним отражением. Путь не заканчивается точно в узле приемника, потому что метод трассировки лучей SBR вычисляет приблизительные пути. Выберите нанесенный на график путь, чтобы просмотреть соответствующие характеристики распространения, которые включают мощность приемника, изменение фазы, расстояние и углы отправления и прибытия.

rtpm.MaxNumReflections = 1;
clearMap(viewer)
raytrace(tx,rx,rtpm)

Анализируйте силу сигнала и эффект материалов

Вычислите мощность приемника с помощью модели распространения, которая ранее была сконфигурирована, чтобы смоделировать идеальное отражение. Затем присвойте более реалистичный вид материала и повторно вычислите мощность приемника. Обновление лучей, показанных в Средство Просмотра. Использование реалистичного отражения материала приводит к потере степени примерно на 8 дБ по сравнению с идеальным отражением.

ss = sigstrength(rx,tx,rtpm);
disp("Received power using perfect reflection: " + ss + " dBm")

Received power using perfect reflection: -70.3924 dBm

rtpm.BuildingsMaterial = "concrete";
rtpm.TerrainMaterial = "concrete";

raytrace(tx,rx,rtpm)
ss = sigstrength(rx,tx,rtpm);
disp("Received power using concrete materials: " + ss + " dBm")

Received power using concrete materials: -78.9591 dBm

Включите потерю погоды

Добавление ухудшений погоды к модели распространения и повторное вычисление мощности приемника результатов в еще 1,5 дБ потерь.

rtPlusWeather = rtpm + propagationModel("gas") + propagationModel("rain");
raytrace(tx,rx,rtPlusWeather)

ss = sigstrength(rx,tx,rtPlusWeather);
disp("Received power including weather loss: " + ss + " dBm")

Received power including weather loss: -80.4766 dBm

Просмотр карты покрытия с одноотражающими путями

Используйте сконфигурированную модель распространения и повторно сгенерируйте карту покрытия, включая пути с одним отражением и ухудшения погоды. Код для повторной генерации результатов покрытия включается, но комментируется. Результаты, получаемые при запуске кода, загружаются из файла, чтобы сэкономить несколько минут времени расчета в презентации примера. Получившаяся карта покрытия показывает мощность приемника в области вокруг приемника без линии участка, проанализированного выше.

clearMap(viewer)

% Load coverage results and plot. Coverage results were generated using
% commented coverage call below, which takes a couple of minutes to complete.
show(tx)

coverageResults = load("coverageResults.mat");
contour(coverageResults.propDataSingleRef, ...
    "Type","power", ...
    "Transparency",0.6)

% coverage(tx,rtPlusWeather, ...
%     "SignalStrengths",-120:-5, ...
%     "MaxRange", 250, ...
%     "Resolution",2, ...
%     "Transparency",0.6);

Графические пути распространения, включающие два отражения

Разверните анализ «точка-точка», чтобы включить пути с двумя отражениями и выбрать меньшее угловое разделение между запущенными лучами для метода SBR. Визуализация показывает два кластера путей распространения, и общая мощность приемника увеличивается приблизительно на 3 дБ по сравнению с путями с одним отражением.

rtPlusWeather.PropagationModels(1).MaxNumReflections = 2;
rtPlusWeather.PropagationModels(1).AngularSeparation = "low";

ss = sigstrength(rx, tx, rtPlusWeather);
disp("Received power with two-reflection paths: " + ss + " dBm")

Received power with two-reflection paths: -77.1445 dBm

clearMap(viewer)
raytrace(tx,rx,rtPlusWeather)

Просмотр карты покрытия с двумя-отраженными Путями

Визуализируйте предварительно вычисленную карту покрытия, которая включает два пути отражения и ухудшения погоды, как определено комбинированной моделью распространения выше. Результирующая карта покрытия показывает мощность приемника в большей области, чем карта покрытия с одноотражающими путями выше. Заметьте, что некоторые области остаются там, где не появляется сигнал.

clearMap(viewer)

% Use pre-loaded coverage results to plot. Coverage results were generated 
% using commented coverage call below,  which may take half hour or so to
% complete depending on the computer hardware. 
show(tx)

contour(coverageResults.propDataDoubleRef, ...
    "Type","power", ...
    "Transparency",0.6)

% coverage(tx,rtPlusWeather, ...
%     "SignalStrengths",-120:-5, ...
%     "MaxRange", 250, ...
%     "Resolution",2, ...
%     "Transparency",0.6);

Просмотр карты покрытия с четырьмя путями отражения

Учитывайте больше путей распространения и генерируйте более точную карту покрытия, увеличивая максимальное количество отражений для анализа трассировки лучей до 4. Снова визуализируйте предварительно вычисленную карту покрытия, которая показывает почти полное покрытие для области вокруг площадки передатчика.

rtPlusWeather.PropagationModels(1).MaxNumReflections = 4;

clearMap(viewer)

% Use pre-loaded coverage results to plot. Coverage results were generated using
% commented coverage call below, which may take one hour or so to
% complete depending on the computer hardware. 
show(tx)

contour(coverageResults.propDataFourRef, ...
    "Type","power", ...
    "Transparency",0.6)

% coverage(tx,rtPlusWeather, ...
%     "SignalStrengths",-120:-5, ...
%     "MaxRange", 250, ...
%     "Resolution",2, ...
%     "Transparency",0.6);

Используйте рулевое управление лучом, чтобы увеличить мощность приемника

Многие современные коммуникационные системы используют методы, чтобы управлять антенной передатчика, чтобы достичь оптимального качества ссылки. Этот раздел использует Toolbox™ Phased Array System, чтобы оптимально управлять лучом, чтобы максимизировать мощность приемника для нелинейной ссылки зрения.

Определите пользовательскую антенну из Report ITU-R M.2412 [1] для оценки 5G радиотехнологий. Создайте равномерный прямоугольный массив 8 на 8 из шаблона элемента, заданного в разделе 8.5 отчета, укажите ее на юг и просмотрите диаграмму направленности излучения.

tx.Antenna = helperM2412PhasedArray(tx.TransmitterFrequency);
tx.AntennaAngle = -90;

clearMap(viewer)
show(rx)
pattern(tx,"Transparency",0.6)
hide(tx)

Функции raytrace с выходом для доступа к вычисленным лучам. Возвращенный comm.Ray объекты включают как геометрические, так и связанные с распространением характеристики каждого луча.

rtPlusWeather.PropagationModels(1).MaxNumReflections = 1;
ray = raytrace(tx,rx,rtPlusWeather);
disp(ray{1})

  Ray with properties:

      PathSpecification: 'Locations'
       CoordinateSystem: 'Geographic'
    TransmitterLocation: [3×1 double]
       ReceiverLocation: [3×1 double]
            LineOfSight: 0
    ReflectionLocations: [3×1 double]
              Frequency: 2.8000e+10
         PathLossSource: 'Custom'
               PathLoss: 117.4546
             PhaseShift: 3.8184

   Read-only properties:
       PropagationDelay: 6.6489e-07
    PropagationDistance: 199.3293
       AngleOfDeparture: [2×1 double]
         AngleOfArrival: [2×1 double]
         NumReflections: 1

Получите угол отхода для пути с одним отражением и приложите этот угол, чтобы направить антенну в оптимальном направлении, чтобы достичь более высокой мощности приемника. Азимут угла отхода смещен азимутом угла физического антенны, чтобы преобразовать его в азимут вектора управления, заданный в локальной системе координат антенны фазированной решетки.

aod = ray{1}.AngleOfDeparture;
steeringaz = wrapTo180(aod(1)-tx.AntennaAngle(1));
steeringVector = phased.SteeringVector("SensorArray",tx.Antenna);
sv = steeringVector(tx.TransmitterFrequency,[steeringaz;aod(2)]);
tx.Antenna.Taper = conj(sv);

Постройте график диаграммы направленности излучения, чтобы показать энергию антенны, направленную вдоль пути распространения. Новая мощность приемника увеличивается примерно на 20 дБ. Увеличенная мощность приемника соответствует пиковому усилению антенны.

pattern(tx,"Transparency",0.6)
raytrace(tx,rx,rtPlusWeather);
hide(tx)
    
ss = sigstrength(rx, tx, rtPlusWeather);
disp("Received power with beam steering: " + ss + " dBm")

Received power with beam steering: -57.5126 dBm

Заключение

Этот пример использовал трассировку лучей для анализа ссылки и покрытия в городском окружении. Анализ показывает:

Этот пример проанализировал мощность приемника и пути для ссылок и покрытия. Чтобы увидеть, как использовать трассировку лучей для конфигурации модели канала для симуляции уровня канала, смотрите пример Indoor MIMO-OFDM Communication Link с использованием Ray Tracing (Communications Toolbox).

Ссылки

[1] Отчет по M.2412 МСЭ-Р «Руководство по оценке технологий радиоинтерфейсов для IMT-2020,» 2017. https://www.itu.int/pub/R-REP-M.2412