MATLAB: How to automatic segmentation signal

segmentationsignal

hi guys
i have a pcg signal
i want to segment the signal into 4 segment, s1 & s2 = 1 segment
how to do it?
thankss
shannon.jpg

Best Answer

The findchangepts (link) function (introduced in R2016a) is perfect for this:
D = load('matlab.mat');
data = D.data;
Fs = D.fs;
t = linspace(0, 1, numel(data))/Fs;
cp = findchangepts(data, 'Statistic','std', 'MinDistance',2E+3, 'MaxNumChanges',7);
figure
plot(t, data)
hold on
YL = ylim;
plot([t(cp); t(cp)], YL(:)*ones(1,numel(cp)), '-r', 'LineWidth',1)
hold off
producing:
how to automatic segmentation signal - 2019 01 22.png
The ‘cp’ vector are the change points corresponding to the vertical red lines in the plot. The findchangepts function takes about 15 seconds to complete this (on my desktop, using tic and toc),
Experiment to get the result you want.