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 の 2 番目のことは、どの入力引数を pwelch に入力すればよいかわからないことです。