Findpeaks function won't plot all peaks.
14 ビュー (過去 30 日間)
古いコメントを表示
I am currently trying to plot all the peaks of my data. Unfortunatley the last peak fails to be plotted using this function. It also happens to be the penultimate value in my data set so I'm not sure if this is where the problem is. It would be apprciated if someone could take a look and see if they can figure out where I am going wrong.
Thanks very much.
ECG = load('ECG.csv')
Frequency = 350; %Frequency of ECG [Hz]
Time = (0:length(ECG)-1)/Frequency; %Number of samples divided by frquency
Amp = ECG(:,1); %ECG Amplitude
figure
plot(Time,Amp)
findpeaks(Amp,Time,'MinPeakProminence',100);
[pks,locs] = findpeaks(Amp,Time,'MinPeakProminence',100);
meanCycle = mean(diff(locs));
BPM = 60/meanCycle;
fprintf('%.2f', BPM)
0 件のコメント
採用された回答
Star Strider
2020 年 1 月 6 日
Thje last R-wave is not actually a peak as findpeaks defines it.
It is necessary to get creative in order to isolate and define it:
ECG = load('ECG.csv')
Frequency = 350; %Frequency of ECG [Hz]
Time = (0:length(ECG)-1)/Frequency; %Number of samples divided by frquency
Amp = ECG(:,1); %ECG Amplitude
figure
% plot(Time,Amp)
findpeaks(Amp,Time,'MinPeakProminence',100);
[pks,locs] = findpeaks(Amp,'MinPeakProminence',100); % Return Indices Instead of Times
ofst = locs(end)+1; % Last Peak ‘loc’ + 1
[Rwaves,idx] = max(Amp(ofst:end)); % Maximum Of Remaining Segment Of ‘Amp’
meanCycle = mean(diff(Time(locs)));
BPM = 60/meanCycle;
fprintf('%.2f\n', BPM)
figure
plot(Time,Amp)
hold on
plot(Time(locs), Amp(locs), '^r', 'MarkerFaceColor','r') % Plot Peaks That ‘findpeaks’ Returns
plot(Time(idx+ofst), Rwaves, 'pr', 'MarkerFaceColor','r') % Plot Last Incomplete Peak
hold off
Producing:
With the last ‘peak’ plotted with a pentagram.
4 件のコメント
その他の回答 (0 件)
参考
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!