Как найти тренд в последовательности в MATLAB?

Question

У меня есть вектор доступности с 1 и 0. 1 означает, что погода хорошая; и 0 означает, что погода не очень хорошая, и место недоступно.

У меня есть step_duration (например, 10 часов). Учитывая step_index (начало шага) как 101, мне нужно найти окно с прямыми 10 часами хорошей погоды.

Ожидаемое решение: с 10-часовой ожидаемой погодой, скажем, вектор доступности [0 0 1 1 1 1 0 0 1 1 1 1 1 1 1 1 1 1]. Итак, индексы, которые мы находим, находятся от 109-118. А соответствующая погодная задержка (учитывая, что мы не могли найти часы подряд) - от индекса 101-108 (т. Е. 9 часов). Мне нужно написать код для такого алгоритма.

Некоторые примеры кода, которые я мог думать о том, как следует (хотя это не совсем то, что я хочу):

window = 0;  % Initialize the counter for finding the weather window. 
     for tt = step_index 
      if Accessibility(tt) == 0 
       % bad weather, move to next index 
       % reset window. 
       window = 0; 
       k = k + 1; 
       possible_window(k) = window; 
       continue 
      else 
       % Increase window 
       window = window + 1; 
       % sote window history 
       k = k + 1; 
       possible_window (k) = window; 
      end 
     end 
     tt = tt + 1; 
end 

Answer 1

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

w = [0 0 1 1 1 1 0 0 1 1 1 1 1 1 1 1 1 1] 
n = 10; 
find(conv(w, ones(1,n)) == n) - n + 1 

это дает вам показатели, где группы 10 начала

Answer 2

Я думаю, что вы ищете пер йод с хорошей погодой, где это будет определено как имеющее достаточное количество в участке 10.

В этом случае вы можете использовать фильтр, чтобы проверить, сколько из них происходит на каждом участке 10. Это может быть что-то вроде этого :

f= filter(ones(10,1),1,x) 

Не уверен в синтаксисе, но проверяйте его. Фильтр обязательно должен вас туда достать.

Вы могли бы следить за нахождение хороших периодов, как так:

find(f==10) 
find(f>=8) 
find(f==max(f)) 

Answer 3

Решение на базе strfind найти начиная индексы N или N + последовательные из них -

%// sample accessibility vector (different from the one in question) 
access_vec = [1 1 1 1 1 1 1 1 1 1 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 0] 

N = 10; %// N or N+ consecutive ones are to be detected 
indices = strfind([0 access_vec],[0 ones(1,N)]) 

Выход -

indices = 
    1 13 

---

Другой пример -

access_vec = [ 0 1 zeros(1,10) 1] 

Выход -

indices = 
    [] 

Answer 4

Пусть данные будут

w = [0 0 1 1 1 1 0 0 1 1 1 1 1 1 1 1 1 1]; 
step_index = 101; 
step_duration = 10; 

Тогда:

d = diff([0 w 0]); 
starts = find(d==1); %// indices of starts of runs of 1's 
ends = find(d==-1); %// indices of ends of runs of 1's (actually starts of 0's) 
keep = find(ends-starts>=step_duration); %// detect runs that are long enough 
starts = starts(keep)+step_index-1; %// keep only those runs and adjust with "-1" 
ends = ends(keep)+step_index-2; %// "-2" because of how ends is defined 

В результате в векторах starts в nd ends.В вашем примере есть только один такое окно, так как они содержат только один элемент:

starts = 
    109 
ends = 
    118