Cum se detectează vârfurile în MATLAB

MATLAB este un pachet software tehnic folosit pentru prelucrarea și analiza semnalelor. O procedură comună în analiza semnalelor este detectarea vârfurilor, adică identificarea maximelor locale – valori mai mari decât punctele de date adiacente – în semnale zgomotoase. De obicei, este necesar să restricționăm detectarea vârfurilor la maximă locală cu lățime (durată în domeniul timp) și înălțime sau magnitudine specifice.

Pasul 1 – Importarea datelor

Definim o sursă de date importând semnalul în MATLAB. De exemplu, creăm un semnal sinus cu zgomot aleator:

my_signal = sin(0:0.1:10) + rand(1,101);

Pasul 2 – Detectarea vârfurilor

Utilizăm funcția findpeaks() cu interpolare quadratică:

[peak_value, peak_location] = findpeaks(my_signal);

Pasul 3 – Înălțime minimă a vârfurilor

Se pot specifica vârfuri cu o înălțime minimă:

[peak_value, peak_location] = findpeaks(my_signal,'minpeakheight',2.5);

Pasul 4 – Distanță minimă între vârfuri

Se caută vârfuri separate de cel puțin un număr de indici:

[peak_value, peak_location] = findpeaks(my_signal,'minpeakdistance',5);

Pasul 5 – Pragul diferenței

Se pot seta vârfuri cu diferență minimă față de punctele adiacente:

[peak_value, peak_location] = findpeaks(my_signal,'threshold',0.5);

Pasul 6 – Numărul maxim de vârfuri

Se limitează numărul de vârfuri detectate:

[peak_value, peak_location] = findpeaks(my_signal,'npeaks',5);

Pasul 7 – Sortarea rezultatelor

Se poate sorta lista de vârfuri:

[peak_value, peak_location] = findpeaks(my_signal,'sortstr','ascend');

Valorile permise pentru sortstr sunt 'ascend', 'descend' și 'none'.