MATLAB: I’m working on a PCG signal I want to calculate and plot the heart rate variability HRV any help please

hearthrvpcgphonocardiogramratevariability

I'm working on a PCG signal I want to calculate and plot the heart rate variability HRV any help please ? is it possible to Work with the same method used with hrv of ECG

Best Answer

‘I'm working on a PCG signal I want to calculate and plot the heart rate variability HRV any help please ?’

Not possible from a phonocardiogram.

‘Is it possible to Work with the same method used with hrv of ECG’

That is the only possible way to do it, since HRV requires a normal QRS complex preceded by a normal P-wave to identify a ventricular depolarization. It is only possible to estimate HRV from the EKG.

Related Solutions

MATLAB: How do i get the y value for a given x value

You have to create a frequency vector as well to plot and interpret the fft correctly as a function of frequency, and for that you need the sampling frequency, (Fs) of ‘recording’. Then use interp1 to get the value at any specific frequency.

For example:

Fs = 1000;                                  % Sampling Frequency
Ts = 1/Fs;
Tv = linspace(0,10*pi,505);
recording = sin(290*Tv).*cos(310*Tv)+rand(size(Tv));
spect = fft(recording);
specta = abs(spect);
Fn = Fs/2;                                  % Nyquist Frequency
Fv = linspace(0, 1, length(spect)/2+1)*Fn;
Iv = 1:length(Fv);                          % Index vector
a200 = interp1(Fv, specta(Iv), 200);        % Find Amplitude At 200 Hz
figure(1)
plot(Fv, specta(Iv), 200,a200,'pr')         % Plot fft &  Amplitude at 200 Hz
grid

MATLAB: Preprocessing ECG signals from 3 channels

You have to have a fourth (reference or ‘ground’ electrode (ideally) to make any sense of those data. First, look up a reference on ‘vectorcardiography’. It’s too complicated to explain in detail here, so I’ll just provide an outline. You did not describe in sufficient detail your experimental setup, so my ability to reply to it is limited.

The cardiac electrophysiologic cycle traces a series of three principal 3D loops, corresponding to the P-wave, QRS complex, and T-wave. The scalar EKG (that you are recording) are projections of this loop on the axes between any two electrodes. Your ‘EKG1-EKG2’ correspond approximately to Lead I, ‘EKG1-EKG3’ to Lead II, and ‘EKG2-EKG3’ to Lead III of the standard scalar electrocardiogram.

Each lead is ‘correct’ and each will give you the ‘heartbeat’. The only way to merge them correctly is to combine them to form a vectorcardiogram. I went into this in considerable detail in Detecting QRS Complex in ECG Signal. I refer you to it.

See the appropriate sections of Braunwald's Heart Disease: A Textbook of Cardiovascular Medicine, 9th Edition for a full discussion.

Best Answer

Related Solutions

MATLAB: How do i get the y value for a given x value

MATLAB: Preprocessing ECG signals from 3 channels

Related Question