Where could I find a MATLAB code for estimating the PRV-high frequency (0.04 - 0.15 Hz) power from a PPG signal? Thanks in advance!
Hi Erick,
a package that comprises both, beat detection from PPG and HRV analysis function, should be this one: https://physionet.org/physiotools/physionet-cardiovascular-signal-toolbox/
Greetings, Sebastian
Hi
The methodology (code) is not different from HRV analyses having RR-interval from ECG.
fs = 4;
time = (Rtimes(1):1/fs:Rtimes(end));
IBIint = interp1(Rtimes,IBIs,time','spline')';
y = detrend(IBIint, 'linear'); % you can apply different detrending methods
window = min(300*fs, length(y));
noverlap = window/2;
nfft = max(256,2^nextpow2(window));
if rem(nfft, 2) == 0
DFT = (nfft/2)+1;
else
DFT = (nfft+1)/2;
end
% you can decide which DFT you want for example DFT = 0:0.001:2;
% for AR metod
order = 20; % you can use different order
[Pyy, f] = pburg(y,order,DFT,fs, 'onesided');
% for FFT method
[Pyy, f] = pwelch(y,window,noverlap,DFT,fs,'PSD','onesided');
VLF_ab = sum(Pyy(f >= 0.003 & f < 0.04));
LF_ab = sum(Pyy(f >= 0.04 & f < 0.15));
HF_ab = sum(Pyy(f >= 0.15 & f < 0.4));
- Badges
- Science topic
- Similar topics
- Biological Science
- Physiology
- Biosignals
- Biomedical Signal Processing