Создание контейнера данных распространения RF
Используйте propagationData объект для импорта и визуализации геолокационных данных распространения. Данные измерения могут быть данными о потерях в тракте, измерениями уровня сигнала, данными отношения сигнал/шум (SNR) или сотовой информацией.
plot | Печать данных распространения на карте |
contour | Отображение карты контуров |
location | Координаты местоположения данных |
getDataVariable | Получение значений переменных данных точек данных в объекте данных распространения |
interp | Интерполяция географических данных |
Запустите Site Viewer с базовыми картами и файлами зданий для Манхэттена. Для получения дополнительной информации о файле osm см. [1].
viewer = siteviewer("Basemap","streets_dark",... "Buildings","manhattan.osm");
Показать площадку датчика на здании.
tx = txsite("Latitude",40.7107,... "Longitude",-74.0114,... "AntennaHeight",80); show(tx)
Создание приемных площадок вдоль близлежащих улиц.
latitude = [linspace(40.7088, 40.71416, 50), ... linspace(40.71416, 40.715505, 25), ... linspace(40.715505, 40.7133, 25), ... linspace(40.7133, 40.7143, 25)]'; longitude = [linspace(-74.0108, -74.00627, 50), ... linspace(-74.00627 ,-74.0092, 25), ... linspace(-74.0092, -74.0110, 25), ... linspace(-74.0110, -74.0132, 25)]'; rxs = rxsite("Latitude", latitude, "Longitude", longitude);
Вычислить уровень сигнала в каждом местоположении приемника.
signalStrength = sigstrength(rxs, tx)';
Создать propagationData объект для хранения данных вычисленного уровня сигнала.
tbl = table(latitude, longitude, signalStrength); pd = propagationData(tbl);
Постройте график данных уровня сигнала на карте в виде цветных точек.
legendTitle = "Signal" + newline + "Strength" + newline + "(dB)"; plot(pd, "LegendTitle", legendTitle, "Colormap", parula);
Приложение
[1] OSM-файл загружается из https://www.openstreetmap.org, что обеспечивает доступ к данным карты, полученным от толпы, по всему миру. Данные лицензированы по лицензии Open Data Commons Open Database License (ODbL), https://opendatacommons.org/licenses/odbl/.
Определите названия и местоположения сайтов вокруг Бостона.
names = ["Fenway Park","Faneuil Hall","Bunker Hill Monument"]; lats = [42.3467,42.3598,42.3763]; lons = [-71.0972,-71.0545,-71.0611];
Создайте массив узлов передатчика.
txs = txsite("Name",names,... "Latitude",lats,... "Longitude",lons, ... "TransmitterFrequency",2.5e9); show(txs)
Создайте карту отношения сигнал-помеха-плюс-шум (SINR), где источник сигнала для каждого местоположения выбирается в качестве узла передатчика с самым сильным сигналом.
sv1 = siteviewer("Name","SINR map"); sinr(txs,"MaxRange",5000)
Возвращает данные распространения SINR.
pd = sinr(txs,"MaxRange",5000); [sinrDb,lats,lons] = getDataVariable(pd,"SINR");
Вычислить ёмкость с помощью теоремы Шеннона - Хартли.
bw = 1e6; % Bandwidth is 1 MHz sinrRatio = 10.^(sinrDb./10); % Convert from dB to power ratio capacity = bw*log2(1+sinrRatio)/1e6; % Unit: Mbps
Создайте новые данные распространения для карты мощностей и просмотрите контурный график.
pdCapacity = propagationData(lats,lons,"Capacity",capacity); sv2 = siteviewer("Name","Capacity map"); legendTitle = "Capacity" + newline + "(Mbps)"; contour(pdCapacity,"LegendTitle",legendTitle);
readtable | rxsite | siteviewer | txsite