心拍解析にかんして

現在、心拍を計測して取得したRR間隔のデータをローレンツプロットをもちいて解析したいと考えているのですが、解析に役立つアプリ等はないでしょうか。

プログラミングをしてみているのですが、うまく結果が出せないので、何かあれば教えてほしいです。

[Shunichi Kusano]

どのような解析をしたいのでしょうか。うまく結果が出せない、とありますがどのような結果がほしいのでしょうか。

今どんなアルゴリズムを組んでいて、こういう結果までは出せているけど、もっとこういう解析をしたい、といった情報を提供いただけると回答がつきやすくなるかと思います。パッと見た限りそのものずばりのアプリはないと思いますが、もっと具体的にやりたいことがわかれば何かアドバイスがつくかもしれないなと思ったので、コメントしました～。

[司 小川]

ローレンツプロットを用いてストレスの判定をしようと考えております。現状としては以下にコードと出力結果を上げますので見ていただければと思うのですが、値がかなり大きな値となってしまっております。計測データには不整脈は検出されておりませんでしたが、外れ値に引っ張られて結果がうまく出ません。

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]

RRがうまくとれていないことが問題と理解いたしました（違ったらすみません）。RRはどのようなアルゴリズムで計測しているのでしょうか。