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