Matlab организует и отображает данные о температуре

Ответ 1 · рейтинг 1 · Rody Oldenhuis

Предполагая, что ваши местоположения не настроены так, чтобы быть ровно в 1 м друг от друга в прямоугольной сетке (это был бы странный офис ...), вы сталкиваетесь с проблемой интерполяции поверхности через разбросанные точки данных. Функция Matlab TriScatteredInterp будет тем, что вам нужно. Просто следуйте примеру по ссылке с некоторыми изменениями:

x = [x values of your locations]
y = [y values of your locations]
z = [all heat readings for all x,y for a single timestamp]

F = TriScatteredInterp(x,y,z);

и постройте как в примере. Вам нужно будет сделать это для всех временных меток, поэтому в псевдокоде:

x = [x values of your locations]
y = [y values of your locations] % assuming they don't change

F = cell(numel(data{1}{1}),1);
for t = 1:numel(data{1}{1}) % loop through all time stamps

    z = cellfun(@(p)p{1}(t), data);
    F{t} = TriScatteredInterp(x,y,z);

end

а затем вы можете построить первый F{1} и добавить ползунок к рисунку, чтобы выбрать другое время.

Обратите внимание, что это предполагает, что все узлы собирают данные с одинаковыми отметками времени. Если это не так (а я подозреваю, что это не так), вам нужно сделать еще один шаг: создать интерполяцию в измерении времени для каждой точки XY.

Это можно легко сделать с помощью spline. Например,

pp = spline(data{1}{1}, data{1}{2});

создает spline через все данные для первого местоположения, так что

z = ppval(pp, [any random time within the interval])

дает интерполированное значение тепла в любой момент времени в интервале. Вы можете сделать это за один раз, выпустив

z = spline(data{1}{1}, data{1}{2}, [any random vector of times] );

Итак, резюмируя:

% interpolate over time

% NOTE: use the maximum first time, and the minimum last time,
% to ensure these endpoints are included in all splines.
minTime = max( cellfun(@(p)p{1}(1), data) );
maxTime = min( cellfun(@(p)p{1}(1), data) );

trange = minTime : [some step] : maxTime;

npts = size(data,1);
z    = cell(npts,1);
for ii = 1:npts
    % creates interpolation for H(t) at equal times
    % trange for location ii
    z{ii} = spline(data{ii}{1}, data{ii}{2}, trange);
end

% interpolate spatially

x = [x values of your locations]
y = [y values of your locations] % assuming they don't change

nts = numel(trange)
F = cell(nts,1);
for t = 1:nts

    zed = cellfun(@(p)p(t),z);
    F{t} = TriScatteredInterp(x,y, zed);

end

% ... and further plotting commands

% get a list of all initial sampling times for all sensors
first_times = cellfun(@(x)   x{1}(1), data);
last_times  = cellfun(@(x) x{1}(end), data);

% find the max/min of these times.
% NOTE: the order of max/min might be counter-intuitive!
[m, im] = max(first_times);
[M, iM] = min(last_times);

% 'im' contains the number of the data set which has the latest initial
% sampling time. 'iM' contains the number of the data set which has the
% earliest final sampling time.

% Now find the indices to the actual sampling times in all data sets
% that correspond to these initial and final times:

proper_first = cellfun(@(x) find(x{1}==m,1), data);
proper_last  = cellfun(@(x) find(x{1}==m,1), data);

% index times and data in dataset k like so:
data{k}{1}( proper_first(k) : proper_last(k) ) % times
data{k}{2}( proper_first(k) : proper_last(k) ) % data