Matlab Заполните предыдущее значение, если отсутствует значение (или ноль) - PullRequest
3 голосов
/ 08 сентября 2011

У меня есть вектор, содержащий временной ряд с различными значениями и некоторые пропущенные значения между ними, которые установлены на ноль:

X=[0,0,2,0,5,0,0,0,4,0];

Я хочу создать новый вектор, в котором будут заполнены отсутствующие значения (нули)по предыдущему значению, если оно существует, так что я получаю новый вектор, похожий на:

Z=[0,0,2,2,5,5,5,5,4,4];

Я просматривал справку и форумы Matlab, как это, чтобы найти аккуратную и подходящую функцию, которая бы решала этуя с решением одной линии или аналогичным, но я не смог сделать это.Я могу решить эту проблему, выполнив несколько различных шагов в соответствии с нижеприведенным описанием, но я предполагаю, что должно быть доступно лучшее и более простое решение?

Текущее решение:

X=[0,0,2,0,5,0,0,0,4,0];
ix=logical(X);
Y = X(ix);
ixc=cumsum(ix);
Z=[zeros(1,sum(~logical(ixc))) Y(ixc(logical(ixc)))];

Это помогает., но это кажется слишком сложным решением простой проблемы, так может кто-нибудь помочь мне с лучшей?Спасибо.

Ответы [ 3 ]

5 голосов
/ 08 сентября 2011

Вот несколько более простая версия с использованием cumsum:

X=[0,0,2,0,5,0,0,0,4,0];

%# find the entries where X is different from zero
id = find(X); 

%# If we want to run cumsum on X directly, we'd 
%# have the problem that the non-zero entry to the left
%# be added to subsequent non-zero entries. Thus, 
%# subtract the non-zero entries from their neighbor 
%# to the right 
X(id(2:end)) = X(id(2:end)) - X(id(1:end-1));

%# run cumsum to fill in values from the left
Y = cumsum(X)

Y =
     0     0     2     2     5     5     5     5     4     4
2 голосов
/ 08 сентября 2011

Вот кое-что, что я написал. Это делает трюк?

% INPUT: the array you would like to populate
% OUTPUT: the populated array

function popArray = populate(array)

popArray = array;

% Loops through all the array elements and if it equals zero, replaces it
% with the previous element
%
% Since there is no element before the first to potentially populate it, this
% starts with the second element.
for ii = 2:length(popArray)
    if array(ii) == 0;
        popArray(ii)= popArray(ii-1);
    end
end

disp(popArray);
1 голос
/ 08 сентября 2011

Позвольте мне предложить другое векторизованное решение (хотя мне оно больше нравится @ Jonas ):

X = [0 0 2 0 5 0 0 0 4 0]

id = find(X);
X(id(1):end) = cell2mat( arrayfun(@(a,b)a(ones(1,b)), ...
    X(id), [diff(id) numel(X)-id(end)+1], 'UniformOutput',false) )
...