Разбиение кусочного сигнала на разные части

Question

Я пытаюсь проанализировать и применить фильтр к сигналу (измерениям). Проблема заключается во время записи измерений, между двумя временными метками, значение не было записано, поэтому программа, используемая для записи этих измерений, только что связала эти 2 свободных конца, это проблема для моего алгоритма. Я хочу разбить сигнал на несколько частей на основе этой аномалии. Мне нужно сделать это автоматически.

Ниже приводится иллюстрация такой аномалии:

Ниже у меня есть ПИК plot(diff(t)):

я хотел бы некоторую помощь с обнаружением этого ненормальности и разрезая сигнал в точке. Поэтому в моем примере сигнал будет разрезан на 3 части.

Ps - есть 2 нарушения в моем 1-ПОС, я только увеличенный на 1, чтобы показать его правильно

Вот код, который я придумал до сих пор:

a = diff(t); 
[value, ind] = findpeaks(a,'thresold',0.5); 
for i = 1:length(ind) 
    T(i) = %not sure how to use here 
end 

редактировать с пиковыми изображениями :

1st peak of 2

Значение Ind: [3557; 7550]

Answer 1

Как насчет использования функции findpeaks?

Ниже приведен пример.

a=[0 0 0 0 0 1 0 0 0 0 1 0 0 0]; 
findpeaks(a,'threshold',0.5); 

Результат будет ниже.

Чтобы найти индекс,

[value, ind]=findpeaks(a,'threshold',0.5); 

Здесь ind будет индекс пика.

В вашем случае,

a=diff(t); 
[value, ind]=findpeaks(a,'threshold',0.5); 
firstT=t(1:ind(1)-1); 
secondT=t(ind(1):ind(2)-1); 
ThirdT=t(ind(2):end);