Интерполяция географических данных
interpvalue = interp(pd,lat,lon)NaN.
interpvalue = interp(pd,Name,Value)
Определите имена и местоположения сайтов вокруг Бостона.
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);
Вычислите мощность приемника данные для каждой площадки передатчика.
maxr = 20000; pd1 = coverage(txs(1),"MaxRange",maxr); pd2 = coverage(txs(2),"MaxRange",maxr); pd3 = coverage(txs(3),"MaxRange",maxr);
Вычислите прямоугольник, содержащий местоположения всех данных.
locs = [location(pd1); location(pd2); location(pd3)]; [minlatlon, maxlatlon] = bounds(locs);
Создайте сетку расположений по прямоугольнику.
gridlength = 300; latv = linspace(minlatlon(1),maxlatlon(1),gridlength); lonv = linspace(minlatlon(2),maxlatlon(2),gridlength); [lons,lats] = meshgrid(lonv,latv); lats = lats(:); lons = lons(:);
Получите данные для каждого передатчика в местоположениях сетки с помощью интерполяции.
v1 = interp(pd1,lats,lons); v2 = interp(pd2,lats,lons); v3 = interp(pd3,lats,lons);
Создайте данные распространения, содержащие минимальные мощности приемника значения.
minReceivedPower = min([v1 v2 v3],[],2,"includenan"); pd = propagationData(lats,lons,"MinReceivedPower",minReceivedPower);
Постройте график минимальной мощности приемника, которая показывает самый слабый сигнал, полученный от любого сайта передатчика. Показанная область может соответствовать области обслуживания триангуляции с использованием трех участков передатчика.
sensitivity = -110; contour(pd,"Levels",sensitivity:-5,"Type","power")
