D = load('t1.mat');
t1 = double(D.t1);
[cp,res] = findchangepts(t1(:,3), 'MaxNumChanges',2, 'Statistic','std');
figure
plot(t1(:,1), t1(:,3))
hold on
plot(repmat(t1(cp,1),1,2)', repmat(ylim',1,numel(cp)), 'LineWidth',2)
hold off
grid
See the documentation for the function for detailks on the various options.
Best Answer