Обнаружение изменений в массиве данных

Question

У меня есть файл test.xlsx, для которого ссылка может быть найдена ниже. Это диаграмма

Согласно диаграмме, после 62-й точки диаграмма имеет сдвиг вверх. Есть ли способ, которым мы можем распознать этот момент посредством обработки изображений?

Я использовал medfilt1(list,5) (link) в MATLAB с окнами разного размера, но результат был не очень хорошим.

Не могли бы вы дать мне какие-либо советы или комментарии?

https://drive.google.com/file/d/0B5QzBOKhSsdgY3Juc3BlaDVuaTA/view?usp=sharing

Answer 1

Смотрите ниже код. Код проводит статистическое сравнение, используя ztest, между скользящим средним данных (индекс-k: индекс) и данными (индексом). Ссылка на ztest: here. ztest тесты для нулевой гипотезы, что data(index) = mean(data(index-k:index)). Если ztest возвращает h=1, нулевая гипотеза отклоняется, в результате чего data(index) <> mean(data(index-k:index)). В случае с вашими данными, ztest будет обнаруживать переключения на более высокие или пониженные значения w.r.t. заданный уровень alpha.

Попробуйте этот код с вашими данными. Код alpha=0.05, код может обнаружить изменения на index = 48 и index = 66. Это места, где происходит переключение (или понижение). Обратите внимание, что описанная выше перестройка происходит в index=66 и НЕ на index=62, как вы заявляли.

Также попробуйте использовать разные значения alpha и k.

data=xlsread('test.xlsx'); 
m = length(data); 

start_index = 11; 
alpha = 0.05; 
k = 10; 

h_array = zeros(length(data),1); 

for index = start_index:m 
    avg = mean(data(index-k:index)); 
    stdev = std(data(index-k:index)); 
    if index > start_index   
     val = data(index); 
     h = ztest(val,avg,stdev,'alpha',alpha); 
     h_array(index) = h; 
    end 
end 

Answer 2

Просто используйте хранимую функцию для этой цели

detrend(list);