MATLAB Answers

心拍解析にかんして

14 views (last 30 days)
司 小川
司 小川 on 20 Jan 2021
Commented: Shunichi Kusano on 22 Jan 2021
現在、心拍を計測して取得したRR間隔のデータをローレンツプロットをもちいて解析したいと考えているのですが、解析に役立つアプリ等はないでしょうか。
プログラミングをしてみているのですが、うまく結果が出せないので、何かあれば教えてほしいです。

  3 Comments

Shunichi Kusano
Shunichi Kusano on 21 Jan 2021
どのような解析をしたいのでしょうか。うまく結果が出せない、とありますがどのような結果がほしいのでしょうか。
今どんなアルゴリズムを組んでいて、こういう結果までは出せているけど、もっとこういう解析をしたい、といった情報を提供いただけると回答がつきやすくなるかと思います。パッと見た限りそのものずばりのアプリはないと思いますが、もっと具体的にやりたいことがわかれば何かアドバイスがつくかもしれないなと思ったので、コメントしました~。
司 小川
司 小川 on 21 Jan 2021
ローレンツプロットを用いてストレスの判定をしようと考えております。現状としては以下にコードと出力結果を上げますので見ていただければと思うのですが、値がかなり大きな値となってしまっております。計測データには不整脈は検出されておりませんでしたが、外れ値に引っ張られて結果がうまく出ません。
RawRRには1列目に時間[㎳],2列目にRR間隔のデータが入っております。
dt = 60;
nbin = floor(RawRR(end,1)/dt); % 解析時間dtごとに解析を何回行うか計算
tindLP = zeros(nbin,1);
LP = zeros(nbin,1);
for ii = 1:nbin
trange = [dt*(ii-1)+1, dt*ii]; % 15 [s]ごとの解析時間を指定
seg_index = find(RawRR(:,1) >= trange(1) & RawRR(:,1) < trange(2)); % 解析時間中に存在するRR間隔データのインデックスを抽出
tindLP(ii,1) = mean(RawRR(seg_index,1)); % 結果表示のための時間インデックス(時間帯の平均値)
segRR = RawRR(seg_index,2); % 指定した時間帯のRR間隔のデータの作成
Pvec = [segRR(1:end-1,1), segRR(2:end,1)];
SD1 = std(Pvec*[1,-1]'/sqrt(2));
SD2 = std(Pvec*[1,1]'/sqrt(2));
LP(ii,1) = SD1*SD2*pi();
end
figure(1);
plot(tindLP, LP, 'o-');
Shunichi Kusano
Shunichi Kusano on 22 Jan 2021
RRがうまくとれていないことが問題と理解いたしました(違ったらすみません)。RRはどのようなアルゴリズムで計測しているのでしょうか。

Sign in to comment.

Answers (0)

Categories

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!