0

我正在做 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 中输入什么输入参数。

4

1 回答 1

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');
于 2014-01-13T09:10:00.217 回答