我正在进行HRV分析，我遇到了一个问题。我有一组RR间隔（每个约100-160），例如RR=[726 679 687 655 719 632 1071 978]等，每个RR中有大约100-160个值。现在我需要使用函数pburg和pwelch进行频率分析，然后将这些函数的PSD分成频带（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-1kHz）。我只是想弄清楚如何在这些频段中评估PSD（我需要0.03 - 0.86 Hz，x轴是0-1kHz，这似乎是极端的）。之后，我需要以某种方式迭代该曲线下的空间？
我想知道pburg第二件事是我不知道我应该输入pwelch的输入参数。
答案 0 :(得分：0)
如果我正确找到你，你想在指定的频率范围内找到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');