JPH0654813A - 複数箇所の体表面から採取した脈波を用いる診断方法並びに装置 - Google Patents

Abstract

(57)【要約】 【目的】 体表面の複数箇所から採取した脈波データか
らカオスアトラクターと１次、２次のリアプノフ数をそ
れぞれ導出し、これらを採取箇所ごとに比較することに
より、心身状態の診断の資料とする。 【構成】 上記を実行するために、診断装置を、被験者
の体表面に同時に装着可能の複数個の脈波センサと、各
センサで採取した上記データからカオスアトラクターと
１次、２次のリアプノフ数をそれぞれ導出する演算手段
と、上記データ及び数値を記憶する記憶手段と、上記ア
トラクター、各数値及び左右比較演算の結果を表示する
表示手段とで構成した。

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は、複数箇所の体表面から
採取した脈波を用いる診断方法並びに装置に関するもの
である。

【０００２】

【従来の技術】従来、エレクトロニクス技術の進展によ
り、脳波図や心電図等の電気的な計測結果に基づいて心
身の異常等を診断することが行われているが、これは、
医師が脳波図や心電図等を観察して診断を下すことが行
われているだけであり、そのため、個別の医師によって
診断が異なったり、診断内容が曖昧になったりするとい
う問題がある。

【０００３】そこで、これらの問題を解決するために、
本願出願人らにより「体表面から採取した脈波及び／又
は心拍を用いる診断方法並びに装置」が発明され、特願
平２−３４０６３５号として出願されている。

【０００４】

【発明が解決しようとする課題】ところが、上記発明を
実施・再検討しているうちに、体表面の対称箇所で採取
した脈波データを比較することによって、上記データに
新たな秩序が含まれていることを発見し、これの解析に
着手したが、前記したように上記発明は体表面の一箇所
から脈波を採取しており、対称箇所からデータを採取す
るには少なくとも２回の測定を要するので各回の測定に
時間差が生じ、そのため、対称箇所のデータではあるが
互いに異なった条件下でのデータになってしまい、特
に、脈波は僅かな刺激でも変動するので、データ解析が
困難になるという問題がある。

【０００５】

【課題を解決するための手段】本発明では、被験者の体
表面の複数箇所から同時に脈波データを採取し、上記デ
ータを採取箇所ごとにそれぞれ数空間に埋めこんで得た
写像と、各採取箇所ごとの上記データがカオスの定義条
件に適合する程度を示す数値とを算出し、各採取箇所間
の上記写像と数値とを比較することにより被験者の心身
の状態を把握することを特徴とする複数箇所の体表面か
ら採取した脈波データを用いる診断方法と、被験者の体
表面に同時に装着可能の複数個の脈波センサと、上記各
センサで採取したデータをそれぞれ数空間に埋めこんで
写像を算出すると共に、各採取箇所ごとの上記データが
カオスの定義条件に適合する程度を示す数値を算出する
演算手段と、各採取箇所の上記写像と数値とを比較演算
する演算手段と、脈波データ及び上記数値を記憶する記
憶手段と、上記各写像、各数値及び比較演算の結果とを
表示する表示手段とを具備することを特徴とする体表面
の複数箇所から同時に採取した脈波データを用いる診断
装置とを提供せんとするものである。

【０００６】

【実施例】本発明の実施例を図面を参照して説明する。

【０００７】図１は、診断装置Ａを示しており、同診断
装置Ａは、左右の脈波センサ 1L,1R、信号増幅器 3L,3
R、アナログ／ディジタル信号変換器 4L,4Rを上記の順
にそれぞれ並列状態に接続し、これらに演算手段及び記
憶手段としての小型コンピュータ５、表示手段としての
ＣＲＴディスプレイ６及びプリンタ７と、これらに電力
を供給する電源８とを接続している。

【０００８】脈波センサ 1L,1Rは、被験者の指尖部に装
着して、同指尖部から脈波を採取するためのもので、図
２で示すように、柔軟かつ遮光性を有する素材を有底筒
状に形成して、被験者の指尖部11を挿入できるようにし
たカバー部12と、同カバー部12の内面に設けたフォトセ
ンサ14とで構成されている。

【０００９】フォトセンサ14は、カバー部12の内面に、
被験者の指尖部11の指腹部13に、それぞれ当接可能とし
た赤外線発光ダイオード15とフォトトランジスタ16とで
構成されており、赤外線発光ダイオード15とフォトトラ
ンジスタ16との光軸15a,16aを指尖部11の内部で、２０
°〜３０°の角度で交差させており、赤外線発光ダイオ
ード15から投射した波長 940nmの赤外線を上記指尖部11
の内部で反射させ、この反射光をフォトトランジスタ16
に入射させることにより、被験者の指尖部11における脈
波を計測し、これを電圧として出力することができる。

【００１０】なお、赤外線発光ダイオード15とフォトト
ランジスタ16との先端部には、それぞれ、凸レンズが装
着されている。

【００１１】図３はフォトセンサ14の回路図である。

【００１２】脈波センサ 1L,1Rは上記のように構成され
ており、被験者の指尖部11に装着するカバー部12を、柔
軟かつ遮光性を有する素材で有底筒状に形成したこと
で、指尖部11への着脱が容易であり、脈波の計測に外光
の影響がなく、フォトセンサ14の装着位置を安定させ、
正確に脈波の計測を行うことができる。

【００１３】なお、脈波センサ 1L,1Rには、体表面に伝
達された拍動を敏感な圧力センサで検出するものなど、
要は体表面から脈波を検出できるものであればよい。

【００１４】信号増幅器 3L,3Rは、上記脈波センサ 1L,
1Rの出力電圧をそれぞれ正確に一定の増幅率で増幅し
て、次のアナログ／ディジタル信号変換器 4L,4Rに出力
するものである。

【００１５】アナログ／ディジタル信号変換器 4L,4R
は、信号増幅器 3L,3Rの出力電圧を１２ビットのデジタ
ル信号に変換し、脈波のデータとして、次の小型コンピ
ュータ５に出力するものであり、同アナログ／ディジタ
ル信号変換器 4L,4Rに小型コンピュータ５からの出力要
求信号が入力するたびに、上記デジタル信号を出力する
ようにしている。

【００１６】小型コンピュータ５は、次に説明する各種
演算機能と、演算の結果等をＣＲＴディスプレイ６及び
プリンタ７に出力する機能と、アナログ／ディジタル信
号変換器 4L,4Rから入力した脈波データと、上記演算の
結果等を記憶する機能とを有している。

【００１７】まず、小型コンピュータ５の全体的な処理
手順を図４を参照して説明する。

【００１８】小型コンピュータ５をスタート(100) させ
ると、まず初期設定(101) が行われ、後述する４次元数
空間の視線方向を設定(102) し、脈波のデータ処理の各
種演算に頻繁に用いられる定数を予め算出して記憶させ
(103) 、ＣＲＴディスプレイ６をオープン(104) して、
メニューを表示させる(105) 。

【００１９】メニューには、図５で示すように、アナロ
グ／ディジタル信号変換器 4L,4Rからのパワーオンレス
ポンスを受信して診断装置Ａを使用可能な状態にするパ
ワーオンシーケンス(107) 、アトラクター表示ウインド
ゥのクリア(108) 、アトラクター表示サイズを大きくす
るズームイン(109) 、アトラクター表示サイズを小さく
するズームアウト(110) 、脈波データをウインドゥ内に
収まるようにして表示させるパルスウエーブ(111) 、デ
ータの保存(112) 、保存したデータの読み込みと、その
データのアトラクターの表示(116) 、４次元数空間の視
線方向変更によるアトラクターの回転(113) （なお、現
在までの試行で発見した脈波データの構造が把握しやす
い角度を３パターン登録し、これらを容易に指定できる
ようにしている）、２００Ｈｚのサンプリング周期で指
定された点数の脈波データを収集して、アトラクターと
してウインドゥ上に描く次データ(114) 、診断装置Ａの
作動を終了させるＥＮＤ(115) 等がある。

【００２０】次に、アトラクター表示のための演算につ
いて説明する。この演算は左右脈波センサ 1L,1Rからの
データについて各自独立して行われ、その一方について
説明する。

【００２１】図６は、上記演算の概要を示しており、脈
波データを２００Ｈｚのサンプリング周期で採取し(12
0) 、この脈波データを４次元数空間に埋めこんでアト
ラクターを作成する(121) 。

【００２２】４次元数空間のアトラクターを３次元数空
間に投影する(122) 。

【００２３】３次元数空間に投影されたアトラクターを
２次元数空間に写影し、画面に出力する(123) 。

【００２４】という順序で行われる。

【００２５】なお、上記演算中に、メニューを呼出し
て、アトラクターを回転させて任意の方向から見ること
ができるようにすることと、アトラクターの拡大、縮小
及びデータの保存と、保存したデータの読み込みと、そ
のデータのアトラクターの表示とが可能である。

【００２６】また、演算速度を高めるために、アナログ
／ディジタル信号変換器 4L,4Rからの脈波データを整数
型とし、上記演算に頻繁に用いられる定数を予め算出し
て記憶させている。

【００２７】なお、上記演算において、脈波データを４
次元数空間に埋めこむというのは、小型コンピュータ５
に次々と入力する脈波データのある時点での脈波の数値
を第１の軸の数値Ｘとし、この時点から（例えば埋め込
み間隔τ＝１０とすれば）１０個目の数値を第２の軸の
数値Ｙ、２０個目の数値を第３の軸の数値Ｚ、３０個目
の数値を第４の軸の数値Ｗとして、これらの数値で４次
元のベクトルを形成させ、このようにして次のベクトル
を１個目、１１個目、２１個目、３１個目の数値でＸ，
Ｙ，Ｚ，Ｗを形成し、こうして作成した多数のベクトル
で脈波データのアトラクターを４次元数空間内に形成す
る。

【００２８】そして、４次元数空間の視方向の単位ベク
トル n1 ＝（n1，n2，n3，n4）を第４の軸に一致させる
ために、次の行列計算を行って、各ベクトル（Ｘ，Ｙ，
Ｚ，Ｗ），(n1 ，n2，n3，n4）を（Ｘ',Ｙ',Ｚ',Ｗ')，
(n1',n2'，n3',n4')に変換する。

【００２９】

【数１】

【００３０】ところで、４次元数空間に形成されたアト
ラクターは，ＣＲＴディスプレイ６では図形として表示
することができないので、次の行列計算を行い３次元数
空間への投影点の座標Ｘ",Ｙ",Ｚ" を算出する。

【００３１】 そして、次式の計算により、上記３次元座標を２次元数
空間に写影した座標Ｘ"'，Ｙ"'を算出する。

【００３２】Ｘ"'＝Ｙ"cosβ−Ｘ"cosα Ｙ"'＝Ｚ" −Ｘ"'tan β ただし、α，βは３次元空間のＸ，Ｙ軸が２次元空間の
Ｘ，Ｙ軸となす角度である。

【００３３】このようにして得た２次元の座標Ｘ"'，
Ｙ"'を表示手段としてのＣＲＴディスプレイ６、又はプ
リンタ７にアトラクターとして出力させる。

【００３４】なお、上記アトラクターは前記したよう
に、左右脈波センサ 1L,1Rからのデータを個別に演算処
理して得られたものであり、ＣＲＴディスプレイ６の画
面に設定した左右ウインドゥに同時に表示され診断の資
料に供される。

【００３５】次に、脈波データがカオスに適合する程度
を示す数値であるリアプノフ数を算出する演算について
説明する。この演算も左右個別に行われる。

【００３６】なお、現実に被験者の体表面から採取した
データと、既に発表されたリアプノフ数計算法のモデル
として用いられるデータとは、次のように異なるため、
上記の計算法を実施することができない。

【００３７】すなわち、既に発表されたリアプノフ数計
算法のモデルとして用いられるデータは、公式（数学モ
デル）によってカオスのデータを生成させており、連続
的な無限のデータをとることができるので、次にのべる
各条件を満たしているのに対し、現実に被験者の体表面
から採取したデータでは、必ずしも数空間の望む位置に
データをとれるとは限らないし、データ数も有限であ
る。

【００３８】 現実に被験者の体表面から採取したデ
ータでは、例えば、前述した脈波から生成した４次元ベ
クトル軌道のある点ＡにおけるベクトルＡ´（単位ベク
トル）が、軌道の直交位置に必ずしもあるとは限らな
い。

【００３９】また、点Ｂにおける次のベクトルＢ”（単
位ベクトル）が点Ｂにおける軌道の直交位置に見つかる
とは限らない。

【００４０】また、ベクトルＢ´とＢ”のなす角が必ず
しも小さいとは限らない。

【００４１】 そこで、各点における最も近似したベ
クトルを探すのであるが、例えば、前記ベクトルＢ”を
探すとき、単にＢ´に近いという条件だけで検索する
と、点Ｂの次の点におけるベクトルをとる可能性があ
る。

【００４２】 データ数が有限であるので、点Ａがデ
ータの終端付近にあると、次の点（点Ａにτを加えた
点）のデータがとれない。

【００４３】 アナログデータがアナログ／ディジタ
ル変換によってディジタル化されるが、これも有限の数
であるので、別々のデータが同じ値になっている可能性
がある。

【００４４】そこで本発明では、リアプノフ数計算に用
いるデータの採用条件を次のように定めている。

【００４５】すなわち、４次元空間において、上記軌道
上に頂点を置き、かつベクトルＢ´を中心線とした頂角
が小さい円錐を設定して、この円錐内にあるベクトルを
採用するようにした。

【００４６】上記によって、前記ベクトルＢ”が非常に
小さくて、軌道方向と略同一方向でないかぎり、隣の点
のベクトルを取る可能性を回避できる。

【００４７】また、各ベクトルの軌道に対する角度θ
は、各データのベクトルを（Ｘ，Ｙ，Ｚ，Ｗ）、軌道方
向のベクトルを（Ｘ',Ｙ',Ｚ',Ｗ')とすると、

【００４８】 で算出することができ、これによってベクトルと軌道の
直交性をチェックすることができる。

【００４９】また、各ベクトルの絶対値の上限と下限と
を設定して、この範囲内にあるベクトルを採用するよう
にすることによって、前記ベクトルＢ”が非常に小さく
て、軌道方向と略同一方向であっても、隣の点のベクト
ルを取る可能性を回避することができる。

【００５０】なお、上記範囲は前記カオスアトラクター
を観察して、最も収束がよい値を設定することができ
る。

【００５１】このようにしてとった各ベクトルの伸び率
を２を底とした対数に変換し、これの算術平均値を１次
のリアプノフ数λ1 とする。

【００５２】上記の計算を実行するために、小型コンピ
ュータ５内に、第７a 図と第７b 図で示すプログラム(5
0)を格納している。

【００５３】すなわち、リアプノフ数λ1 の演算をスタ
ートさせると(51)、まず、前記のベクトル軌道の始端近
傍に、データの採用条件を判断するための基準となる点
Ａを設定する(52)。なおこの点Ａは計算の進行に伴って
先送りされる。

【００５４】次に、次の点Ｂ（τだけ後の点）をとる余
裕があるか否かを判断し(53)、余裕がある場合は(53)Y
この点Ｂを次の点として採用し(54)、この点Ｂで仮のデ
ータを探し(55)、見つかれば(55)Y 、次の点から仮のデ
ータを探し(56)、見付からなければ先頭のデータから検
索する(57)。

【００５５】次に、このようにして見つけたデータが前
述のデータの採用条件に適合していれば(58)Y 、このデ
ータを採用し(59)、なければ(58)N 、データ採用条件中
のベクトルの大きさの範囲を更新し(60)、この範囲の上
限を逸脱していなければ(61)N 、ステップ(55)に戻り、
逸脱していると(61)Y 、基準点Ａを次の点Ｂに移して(6
2)、ステップ(53)に戻る。

【００５６】このようにして採用したデータは、前述し
たベクトルＡ´がベクトルＢ´に発展する余裕があるか
否かを判断され(63)、余裕があれば(63)Y 、ベクトルＢ
´の軌道との直交性をチェックする(64)。

【００５７】なお、ステップ(63)で余裕がない場合(63)
N 及びステップ(64)で直交していない場合(64)N は、ス
テップ(55)に戻る。

【００５８】そして、後述のリトライ中であるか否かを
判断して(65)、リトライ中であれば(65)Y 、前回のベク
トルとのなす角を算出して(66)、この角が大きい場合(6
6)Lは、ステップ(55)に戻り、小さい場合(66)S 及び前
記ステップ(65)でリトライ中でない場合(65)N は、ベク
トルＢ´をベクトルＡ´から発展したベクトルとして確
定する(67)。

【００５９】そして、リトライのために点Ａとベクトル
Ａ´とを記憶させておく(68)。

【００６０】なお、演算の終末に際し、ステップ(53)で
次の点をとる余裕がなくなれば(53)N 、演算を終了させ
る(E) 。

【００６１】次に、点Ｂにおける仮のベクトルＢ”をさ
がし(70)、このベクトルＢ”がデータ採用条件に適合し
ていなければ(71)N 、上記条件のベクトルの大きさの範
囲を更新し(72)、それでも上限を逸脱していると(74)Y
、リトライフラグをセットして (74)、ステップ(55)
に戻る。

【００６２】また、ステップ(73)で上限を逸脱していな
ければ(74)N 、ステップ(70)に戻る。

【００６３】そして、ステップ(71)でデータ採用条件に
適合したベクトルＢ”があると(71)Y 、このベクトル
Ｂ”が軌道と直交しているか否かを判断し(75)、直交し
ていなければ(75)N 、ステップ(70)に戻り、直交してい
ると(75)Y 、点Ｂにおける各ベクトルＢ´，Ｂ”のなす
角度を算出し(76)、この角度が充分に小さいか否か、す
なわち、前述の円錐内にＢ”が入っているか否かを判断
し(77)、角度が小さくない場合(77)N は、ステップ(70)
に戻り、小さい場合(77)Y は、点Ａから点Ｂに移動した
ことによるベクトルＡ´からベクトルＢ´への伸び率
を、２を底とした対数に変換し(78)、この数値を算術平
均して１次のリアプノフ数λ1 とする(79)。そして、現
在計算している点Ｂがデータの終端であるか否かを判断
し(80)、終端でなければ(80)N 、上記の点Ｂを次回計算
の基準点（前記の点Ａに相当）に代入し(81)、ベクトル
Ｂ”を単位ベクトルに変換して次回計算の基準ベクトル
（前記のベクトルＢ´に相当）に代入し(82)、リトライ
のために点ＢとベクトルＢ´とを記憶させ(83)、点Ｂの
次の点を算出し(84)、算出された次の点の基準ベクトル
から発展したベクトルを算出し(85)、次回計算において
用いるベクトルの直交チェック用の座標を更新して(8
6)、ステップ(70)に戻り、上記計算を繰返す。

【００６４】なお、ステップ(80)で点Ｂがデータの終端
であると判断されると(80)Y 、演算を終了させる(E) 。

【００６５】つまり、ステップ(52)〜(68)では、主とし
て、次の点におけるベクトルＢ´を参照してデータ採用
条件に適合する基準点を検索し、ステップ(70)〜(77)で
は同条件に適合する発展したベクトルを検索するように
しており、上記のように、採用するデータの採用条件を
設定し、この条件に適合したデータを検索することによ
って、不適格なデータが計算に算入されることを防止
し、しかも、適格なデータが見つからない場合は、その
点におけるデータ検索を放棄するのではなく、採用条件
のベクトルの大きさの設定範囲を広げて検索して、次善
のデータを採用するようにしたことで、生体からの採取
データから高い精度で１次のリアプノフ数λ1 を算出す
ることができる。

【００６６】次に、２次のリアプノフ数λ2 の計算につ
いて説明する。

【００６７】２次のリアプノフ数λ2 の計算も、基本的
には前述した１次のリアプノフ数λ1 の計算と同じであ
り、点Ａ1 においてデータ採用条件に適合したベクトル
Ａ2とＡ3 をとって４次元数空間内に三角形Ａ1-Ａ2-Ａ3
を形成し、次に、ベクトルＡ2 とＡ3 とが次の点Ｂ1
までに発展したベクトルＢ2 とＢ3 とで三角形Ｂ1-Ｂ2-
Ｂ3 を形成し、更に、点Ｂ1 で新たにとった適格なベク
トルＢ´とＢ”とで三角形Ｂ1-Ｂ'-Ｂ”を形成して、三
角形Ｂ1-Ｂ2-Ｂ3 と三角形Ｂ1-Ｂ'-Ｂ”のなす角度が充
分に小さいとき、三角形Ａ1-Ａ2-Ａ3 から三角形Ｂ1-Ｂ
2-Ｂ3 への面積の伸び率を算出し、これを２を底とした
対数に変換し、この数値を算術平均することによって２
次のリアプノフ数λ2 を算出する。

【００６８】なお、各三角形のなす角度は、ベクトルＢ
2 とＢ3 の合成ベクトルと、ベクトルＢ´とＢ”の合成
ベクトルとがなす角度をもって、上記三角形がなす角度
とした。

【００６９】次に、本診断装置Ａを用いた診断例につい
て図９〜図１１を参照して説明する。

【００７０】図９は、健康者のアトラクターを示してお
り、明らかに左右の形状が相似している。

【００７１】図１０は左側脳梗塞患者のカオスアトラク
ター、図１１は右側脳出血患者のカオスアトラクター、
図１２は左側脳梗塞患者のカオスアトラクター、図１３
は脳梗塞の病歴を有する人のカオスアトラクターであ
り、いずれも左右のアトラクターの形状が著しく異なっ
ている。

【００７２】上記のように、健康者では左右のアトラク
ターが相似しているのに対し、脳の一方に障害を有する
人のアトラクターが異なっており、この違いを読み取る
ことによって、脳の障害を診断することができる。

【００７３】本実施例では上記のように、左右指尖部に
装着した脈波センサ 1L,1Rで採取した被験者の脈波デー
タから、左右の脈波のカオスアトラクターと、１次、２
次のリアプノフ数λ1,λ2 とを個別に算出してそれぞれ
ＣＲＴディスプレイの左右ウインドゥ上に表示すること
ことができ、かかるカオスアトラクターと、１次、２次
のリアプノフ数λ1,λ2 とを、それぞれ左右比較するこ
とにより、被験者の心身の状態を診断することが可能で
あり、特に、一見ランダムに見える上記データから、論
理的に秩序を抽出し、同秩序に基づいて、心身の状態を
示す図形及び数値を表現し、これを左右比較することに
よって、診断者の違いによる差異のない、極めて客観的
な診断を下すことができる。

【００７４】

【発明の効果】本発明によれば、体表面の複数箇所から
同時に採取した脈波データからカオスアトラクターとリ
アプノフ数とをそれぞれ導出し、これらを採取箇所ごと
に比較することによって、一箇所から採取したデータで
は検出できない症状や心身の状態を検知することができ
る。

【００７５】特に、脈波センサを２個用いて、体表面の
対称箇所、例えば、左右手の同一指尖部から同時に脈波
データを採取することにより、採取の時間差による状態
の変動のない同一条件下における対称箇所の脈波データ
採取することができ、かかる脈波データを解析すること
によって被験者の状態を正確に把握することが可能にな
るという効果がある。

【図面の簡単な説明】

【図１】本発明に係る診断装置の構成を示す説明図。

【図２】脈波センサの構造を示す断面説明図。

【図３】フォトセンサの回路図。

【図４】脈波データの全体的な処理手順を示す説明図。

【図５】メニューの説明図。

【図６】アトラクター表示のための演算手順の説明図。

【図７】リアプノフ数を求める演算処理のフローチャー
ト。

【図８】リアプノフ数を求める演算処理のフローチャー
ト。

【図９】健康者のカオスアトラクター。

【図１０】左側脳梗塞患者のカオスアトラクター。

【図１１】右側脳出血患者のカオスアトラクター。

【図１２】左側梗塞患者のカオスアトラクター。

【符号の説明】

Ａ：診断装置 １：脈波センサ ５：小型コンピュータ（演算手段、記憶手段） ６：ＣＲＴディスプレイ（表示手段） 11：指尖部 14：フォトセンサ 15：赤外線発光ダイオード 16：フォトトランジスタ

Claims (2)

【特許請求の範囲】
1. 【請求項１】 被験者の体表面の複数箇所から同時に脈
   波データを採取し、上記データを採取箇所ごとにそれぞ
   れ数空間に埋めこんで得た写像と、各採取箇所ごとの上
   記データがカオスの定義条件に適合する程度を示す数値
   とを算出し、各採取箇所間の上記写像同志並びに数値同
   志を比較することにより被験者の心身の状態を把握する
   ことを特徴とする複数箇所の体表面から採取した脈波を
   用いる診断方法。
2. 【請求項２】 被験者の体表面に同時に装着可能の複数
   個の脈波センサと、上記各センサで採取したデータをそ
   れぞれ数空間に埋めこんで写像を算出すると共に、各採
   取箇所ごとの上記データがカオスの定義条件に適合する
   程度を示す数値を算出する演算手段と、各採取箇所間の
   上記写像同志並びに数値同志を比較演算する演算手段
   と、脈波データ及び上記数値を記憶する記憶手段と、上
   記各写像、各数値及び比較演算の結果とを表示する表示
   手段とを具備することを特徴とする体表面から採取した
   脈波を用いる診断装置。