Scribd
Upload
277 views2 pages

Ecg Matlab Code

This document contains code to analyze an ECG signal. It loads ECG data, sets the sampling rate and time vector. It plots the original signal and extracts a single period and PQRST cycle to plot. It then analyzes heart rate by counting dominant peaks above a threshold in the signal and calculates beats per minute (BPM) based on the peak count and duration of the recording.

Uploaded by

mak
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
277 views2 pages

Ecg Matlab Code

This document contains code to analyze an ECG signal. It loads ECG data, sets the sampling rate and time vector. It plots the original signal and extracts a single period and PQRST cycle to plot. It then analyzes heart rate by counting dominant peaks above a threshold in the signal and calculates beats per minute (BPM) based on the peak count and duration of the recording.

Uploaded by

mak
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
You are on page 1/ 2

close all;

clear all;
clc;
%
%S=load('I:\BIOM_Signal_processing\Hw5\ECGsignal_1.mat')
y1=xlsread('J:\BIOM_Signal_processing\Hw5\ECGsignal_1.xls');
fs = 250 % find the sampling rate or frequency
T = 1/fs;% sampling rate or frequency
% find the length of the data per second
N = length(y1);
ls = size(y1);
t = (0 : N-1)/fs;% sampling period

figure; %subplot(1,2,1);
%subplot(1,2,1);
plot(t,y1);
title ('plot of the orignal ECG signal');
xlabel ('time (sec)');
ylabel ('ECG Amplitute (mv)');
grid on;
figure;
%%% Create a period
y1new = y1(10:350);
ax = axis; axis([ax(1:2) -3.2 3.2])
t2 = (0 : length (y1new)-1)/fs;% sampling period
%subplot(1,2,2);
plot (y1new);
title ('plot one typical period of the signal amplitude spectrume of the ECG
signal');
xlabel ('time (msec)');
ylabel ('ECG Amplitute (mv)');

figure;
%%% Create a period
y1new1 = y1(40:200);
ax = axis; axis([ax(1:2) -3.2 3.2])
t2 = (0 : length (y1new)-1)/fs;% sampling period
%subplot(1,2,2);
plot (y1new1);
title ('plot PQRST cycle of the ECG signal');
xlabel ('time (msec)');
ylabel ('ECG Amplitute (mv)');

%% heart rate analysis


% count the dominat peak
beat_count =0;
for k = 2 : length(y1)-1
%the peak has to be greater than 1 and greater than the value before it and
greater then the value after it.
if(y1(k)> y1(k-1) && y1(k) > y1(k+1) && y1(k)> 1)
beat_count = beat_count +1;
end
end
display (k);
disp('dominant peaks');
%% divide the peak count by the duration in minute
duration_in_sec = N/fs;
duration_in_minute = duration_in_sec/60;
BPM = beat_count/duration_in_minute;

You might also like