Q&A
ご丁寧な解説をありがとうございます。
配列など私の勉強不足で申し訳ありません…
入力信号の心電波形の画像を追加しました。
仰る通り、RR間隔の判別方法を変更した方が良いと思いました。プログラムを調べて考え直してみます。
前提・実現したいこと
入力された生体信号(心電波形)のRR間隔の時間を出力させるために
入力信号のデータを取得して、波形のグラフを作成したいです。
私の作成したコードではArduinoに入力する前にオシロスコープで取得している波形とは異なってしまいます。
コンパイル時のエラー等はありません。
入力波形は心電波形です。
解決したコード
int now=0 ;
long i=0 ; //配列番号の管理
byte j=0 ; //フラグ管理
unsigned long tNow,tPrev;
word peak;
unsigned long time ;
//int pin=10 ;
void setup() {
// put your setup code here, to run once:
Serial.begin(9600);
}
void loop() {
// put your main code here, to run repeatedly:
now=analogRead(A5) ;
time = millis() ;
Serial.print(time/1000) ;
Serial.print(",") ;
Serial.println(now*0.005);
//Serial.println(now);
//delay(10);
if(now>peak) {//ピーク値更新なら
peak=now;//保存して
tNow=millis();//時間を取得
j=1;//取得した
}
if(now<peak*7/10 && j==1) {//7割りに落ちた時にピーク確定。→集計
if(i){
//Serial.print("count:");
//Serial.print(i);
//Serial.print(" peak:");
//Serial.print(peak);
Serial.print(" RR:");
Serial.print(tNow-tPrev);
Serial.println("ms");
}
tPrev = tNow ;
peak=peak*8/10;//8割以上に上がったら再取得
i++;
j=0;
}
}
補足情報（FW/ツールのバージョンなど）
ArduinoUNO
回答4件
0
シリアルモニタのボーレートが57600になってますね
これはArduino側のコードで書かれたボーレートと一致させる必要があります
で、Serial.begin、と Serial.end、ってのは、例えば電話をかける、切る、と同じようなもんなんで、一つ送信するたびにbegin,endするんじゃなくて、最初にbeginしとけばそんでいいもんです
んで、endの方は実行しないでいいです
投稿2019/01/23 06:13
総合スコア87749
0
やりたいことは、直近の山の頂上から20%下がったときの時間間隔を測定したい、ということだと読み取りましたがよろしいですか?
まず。
コンパイル時のエラー等はありません。
とのことですが、
重大なバグがあります。
ArduinoのベースとなるC++言語では「コンパイルエラーにはならない(特別な使い方をするためにエラーにしない)けれど普通はやってはいけないこと」がゴマンとあって、その一つに「確保した配列の範囲外にアクセスすること」というのがあります。
Arduino
1unsigned long t[200] ; 2//<<<略>>> 3 t[i] = millis(); 4 5 //i=0の時は計算できないので分岐 6 7 if (i != 0) { 8 to = t[i] - t[i - 1]; //入れるに格納された秒数からRR間隔を計算 9 //ケーブルを通してPCに出力 10 Serial.println("RR"); 11 Serial.println(to); 12 } 13 i++; //次の配列へ
t[200]として宣言されていますので、配列の添字としては(原則として)0～199の範囲外で使ってはいけません。しかし、t[]の添字として使われているiは、何の制限もなくどんどん増えていきます。ダメです。
配列の範囲外にアクセス、特に書き込みを行ったときに何が起こるか、全く予想出来ない事態になります。他の変数の値を滅茶苦茶にしたり、あるいはプログラムの動作を不可能にしたり、いろいろ悪さをします。
この部分は「前回の時刻-今回の時刻」の計算をしているだけですから、過去の値を配列で全部保存している必要はないはずです。
Arduino
1unsigned long tNow, tPrev; 2//<<<略>>> 3 tNow = millis(); 4 5 //i=0の時は計算できないので分岐 6 7 if (i != 0) { 8 to = tNow - tPrev; //前回の時刻から引けばRR間隔が求まる 9 tPrev = tNow; //前回の時刻として保存 10 //ケーブルを通してPCに出力 11 Serial.println("RR"); 12 Serial.println(to); 13 } 14 i++;
で十分では。
つぎ。
Arduino
1 a = analogRead(A5); //アナログ入力の代入 2 //閾値を決める 3 if (a - x > 0) { 4 x = a; //波形の山が下り始めるところを見つける 5 if (x - threshold > 0) { 6 threshold = a; //R波のみ判定 7 } 8 }
この部分から読み取れることは、
if (a - x > 0) {
aの方が大きければxを書き換えている、つまりxは測定開始以来最大のシグナル値を保持している。
if (x - threshold > 0) {
xの方が大きければthresholdを書き換えている。つまり、thresholdも同じことで測定以来最大のシグナル値を保存している。
xの存在は無駄です。そしてなにより、波形の山が下がり始めるところをみつけているわけではありません。
結果として、
最初低い山があると(1)でthresholdが確定して、その0.8倍で最初のRR判定が行われます(2)。
次に高い山があると(3)thresholdが更新され、以降低い山があってもRR判定には引っかからない、ということになります。これは意図した動作でしょうか。Wikipediaの心電図の項の画像のようなのが典型パターンとすると、それでも破綻はしないのかも知れませんが、しかしRの高さが20%以上変動があるといろいろと支障が出そうな気がします。
この辺の判別方法は要検討ではないでしょうか。
なお、データをどう処理すればいいのかはプログラミングの問題ではありません。プログラムは「あなたがやらせたいことをコンピュータにわかるように書き下しただけの文書であって、プログラムがなにかを解決してくれるわけではなく、コンピュータに何をやらせればいいかはプログラム以前に決まっているはずの話です。
投稿2019/01/26 01:03
総合スコア7629
0
シリアルモニタのボーレートを9600にする。
Serial.begin(9600);をsetup()に移動。
y_waiwaiさんの回答見て存在に気づきましたけど
Serial.end()を抜く
投稿2019/01/23 05:58
編集2019/01/23 06:17
総合スコア13521
0
自己解決
こちらのコードにて解決することができました！
お力添えをしてくださった回答者様ありがとうございました。
int now=0 ;
long i=0 ; //配列番号の管理
byte j=0 ; //フラグ管理
unsigned long tNow,tPrev;
word peak;
unsigned long time ;
//int pin=10 ;
void setup() {
// put your setup code here, to run once:
Serial.begin(9600);
}
void loop() {
// put your main code here, to run repeatedly:
now=analogRead(A5) ;
time = millis() ;
Serial.print(time/1000) ;
Serial.print(",") ;
Serial.println(now*0.005);
//Serial.println(now);
//delay(10);
if(now>peak) {//ピーク値更新なら
peak=now;//保存して
tNow=millis();//時間を取得
j=1;//取得した
}
if(now<peak*7/10 && j==1) {//7割りに落ちた時にピーク確定。→集計
if(i){
//Serial.print("count:");
//Serial.print(i);
//Serial.print(" peak:");
//Serial.print(peak);
Serial.print(" RR:");
Serial.print(tNow-tPrev);
Serial.println("ms");
}
tPrev = tNow ;
peak=peak*8/10;//8割以上に上がったら再取得
i++;
j=0;
}
}
投稿2019/01/30 05:51
総合スコア12
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。
2019/01/23 06:27
2019/01/23 07:18