Интерполяция географических данных
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")
