rr =[2289.1 742.19 738.28 750 722.66 761.72 800.78 828.13 800.78 824.22 832.03 808.59 746.09 757.81 781.25 781.25 746.09 750 785.16 785.16 812.5 804.69 828.13 1648.4 ]; % RR interval
fs = 256;
NFFT = 1024*4;
[PSD,F] = pwelch(rr,length(rr),[],NFFT,fs); %uses a hamming window
VLF = [0 0.04];
LF = [0.04 0.15];
HF = [0.15 0.4];
% find the indexes corresponding to the VLF, LF, and HF bands
iVLF= (F>=VLF(1)) & (F<=VLF(2));
iLF = (F>=LF(1)) & (F<=LF(2));
iHF = (F>=HF(1)) & (F<=HF(2));
% calculate areas, within the freq bands (ms^2)
aVLF=trapz(PSD(iVLF));
aLF=trapz(PSD(iLF));
aHF=trapz(PSD(iHF));
aTotal=aVLF+aLF+aHF;
%calculate LF/HF ratio
lfhf = aLF/aHF;
disp(['ratio = ',num2str(lfhf)])
LFnu = (aLF / (aLF + aHF)) * 100;
HFnu = (aHF / (aLF + aHF)) * 100;
Что-то не так с этим кодом, я не могу найти пик vlf или lf или hf на графике plot(F,PSD); Пожалуйста, помогите мне найти vlf, lf, hf из интервала RR.