Создайте контейнер данных распространения RF
Используйте propagationData объект для импорта и визуализации геолокированных данных распространения. Данные измерения могут быть данными потерь пути, измерениями уровня сигнала, данными отношения сигнал/шум (ОСШ) или сотовой информацией.
plot | Постройте данные распространения на карте |
contour | Отобразите контурную карту |
location | Координаты местоположения данных |
getDataVariable | Получите значения переменных данных точек данных в объекте данных распространения |
interp | Интерполяция географических данных |
Запустите Средство Просмотра с основными картами и файлами создания для Manhattan. Для получения дополнительной информации о файле 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
