基于pburg和pwelch函数的Matlab HRV频域分析

Question

我正在做HRV分析，我遇到了一个问题。我有几组RR间隔(每组大约100-160个)，例如RR=[726 679 687 655 719 632 1071 978]等，每个RR大约有100-160个值。现在我需要使用函数pburg和pwelch进行频率分析，然后将这些函数的功率谱分解成频带(LF = 0.03 - 0.86 Hz；HF = 0.86 - 1.91 Hz)。我的采样频率是2000。

到目前为止，我的那部分代码看起来是这样的：

[Pxx1,w1]=pburg(RR,16,256,2000);
[Pxx2,w2]=pwelch(RR);

此函数返回y轴上的PSD和x轴上的频率(0-1 1kHz)。我只是想不出如何在这些频段评估PSD (我需要0.03 - 0.86 Hz，而x轴是0-1 1kHz，这似乎很极端)。之后，我可能需要以某种方式对该曲线下的空间进行itegrate？

第二件事是我不知道应该把什么输入参数输入到pwelch中。

Answer 1

如果我没记错的话，你想要在指定的频率范围内找到PSD。

如果我的假设是正确的，您可以执行以下操作。

%[Pxx,W] = PWELCH(double(data),window_length,overlap_size,[low_freq:high_freq], sampling_freq);

[Pxx,W] = PWELCH(double(RR), 512, 256, [30:860], 1000);

plot(W,(Pxx));

xlabel('Hz'); ylabel('Power');