JPH0751240A

JPH0751240A - 脳波の多次元相関分析方法

Info

Publication number: JPH0751240A
Application number: JP31756892A
Authority: JP
Inventors: Toshimitsu Musha; 利光武者
Original assignee: Individual
Current assignee: Individual
Priority date: 1992-11-26
Filing date: 1992-11-26
Publication date: 1995-02-28

Abstract

(57)【要約】【目的】脳波の相関分析を時間的、空間的に、且つ明
確に行える。【構成】脳内における活動ニューロンの電気的状態
と、人間の意識状態とを対応づける相関を分析する方法
であって、頭皮面において計測した電位の波形から、極
値を形成するスペクトルを抽出する段階を先ず設け、前
記のスペクトルが位置的に分布する相関を、時間的に分
析する段階を続けて設け構成しており、前記の相関を時
間的な分布図に表示し、前記の分布図を格納し得る様に
なっていることを特徴とする。

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は、ＭＥ（メディカル・エ
レクトロニクス：医用電子機器）分野において、人間の
脳波を計測し脳内での活動の相関を分析する脳波の相関
分析方法、特に、この活動の部位及びその時間的、空間
的（即ち、多次元）な分布をディジタル信号処理により
分析する脳波の相関分析方法であり、更に、これらの相
関の在り方は思考・認識・記憶の想起・快不快及び・精
神的な疲労緊張などの内容と現象論的に対応づけられ
る。

【０００２】

【従来の技術】頭皮面に現れる電位は、脳内にある多数
のニューロン（神経細胞）の電気的な活動を反映してお
り、思考過程、認識、記憶の想起、またある種の感情に
ともなって、これらニューロンの活動に連携が発生する
ことは医学者らにより既に広く知られていた。

【０００３】図１１は従来の脳波の相関分析方法を示す
図である。従来の脳波の相関分析方法は、図１１に示す
様に脳波を頭皮面の測定器Ｓで電気信号に変換し単一の
チャネルを有する計測装置ＤＥＴ１を介して得られたデ
ータを表示装置１に導入する。表示装置１にはこのデー
タを準備１する段階２を先ず設け、そのまま画像に変換
し表示する段階３を続けて設けて、この画像を医学者が
観察して整理１する段階４を最後に設け構成しており、
医学者は自らの知識に従いこの画像から特徴を見出し自
らの経験により傾向を分類する様になっていた。

【０００４】この様な方法により、頭皮面の電位は意識
状態及び心理的な状態に従い局所的に強まることが、研
究者らによりこの電気信号の波形において観測され学会
等に発表されていた。

【０００５】しかし、従来の脳波の相関分析方法には次
の様な問題点があった。（１）単一のチャネルを介して散発的に計測を繰り返す
ので、多数のグラフを人手により整理しなければならず
時間がかかる。

【０００６】そこで、本発明人は既に前記の事実と問題
点に基づき、これらニューロンの活動による電位の波形
を頭皮面において計測する測定装置、及び脳内の活動部
位の位置・強度を推定する方法を開発していた。この測
定装置及び方法は頭皮面に複数個の電極を付け頭皮，頭
蓋，脳内の影響を除く処理を施して、直接またはテレメ
ータなどによって間接的に多チャンネルの増幅器を経由
してデータレコーダに時間的に変化する電位の波形を記
録するものである。

【０００７】これらの特許出願済の装置及び方法の詳細
については、次のそれぞれの広報に記載されている。（ａ）計測方法：特開平２−３１７３９号公報（ｂ）計測装置：特開平２−３１７３６号公報（ｃ）計測方法：特開平３−９９６３０号公報

【０００８】ここで、本発明人は次の仮定１を提起し
た。仮定１：脳内にあるニューロン（神経細胞）の電気的な
活動を反映するいくつかの電気的な双極子がシナプスに
発生する。即ち、これらの装置及び方法を使用し、人間
の脳内にある多数の電気的な双極子から発生する脳波を
頭皮面から求めるのである。また、この脳波による電位
の変化は自発的なものＥＥＧと外部の要因による誘発電
位とが有り、本発明人によるこれらの装置及び方法はこ
のいづれにも適用する様になっている。

【０００９】

【発明が解決しようとする課題】しかし、本発明人のこ
の測定装置及び方法をもってしても、前記した従来の脳
波の相関分析方法に使用したのでは依然として次の様な
問題点があった。（２）個人個人の知識と経験により分析を行うので、医
学者同志で意見が異なり結論がまとまらない。（３）また、分析した根拠や過程の保管を人間の記憶に
頼る事になるので、再現しにくく科学的とは言えない。

【００１０】本発明は以上の問題点に鑑み、この分布を
多次元（即ち、時間及び空間）的にディジタル信号処理
し、時間的、空間的な相関の状態について多様性に富む
明確で再現性のある分析を行うことを目的とする。

【００１１】

【課題を解決する手段】

（１）脳内における活動ニューロンの電気的状態と、
人間の意識状態とを対応づける相関を分析する方法であ
って、頭皮面において計測した電位の波形から、特徴と
なる要素を抽出する段階を先ず設け、前記の要素が位置
的に分布する相関を、時間的に分析する段階を続けて設
け構成しており、前記の相関を時間的な分布図に表示
し、前記の分布図を格納し得る様になっていることを特
徴とする脳波の相関分析方法。（２）前記の特徴となる要素は極値を形成するスペク
トルである前項（１）に記載の脳波の相関分析方法。（３）前記の特徴となる要素は前記のスペクトルを有
する信号のポテンシャル及び位相である前項（１）及び
（２）に記載の脳波の相関分析方法。（４）前記の分析する段階は位置的に結合するコヒー
レンスの係数を時間を追って分析する手段である前項
（１）ないし（３）に記載の脳波の相関分析方法。（５）前記の分析する段階は前記のスペクトルを時間
を追ってバリアンス分析する手段である前項（１）ない
し（３）に記載の脳波の相関分析方法。

【００１２】

【作用】頭皮面において計測した電位の波形から極値を
形成するスペクトルを有する信号をディジタル信号処理
により抽出し、この特徴的な信号が分布する相関を統計
学的な演算処理により時間を追い分析して、分析した相
関を画像処理により分布図に表示することができる。

【００１３】図１１は従来の脳波の相関分析方法を示す
図である。従来の脳波の相関分析方法は、図１１に示す
様に脳波を頭皮面の測定器Ｓで電気信号に変換し単一の
チャネルを有する計測装置ＤＥＴ１を介して得られたデ
ータを表示装置１に導入する。表示装置１にはこのデー
タを準備１する段階２を先ず設け、そのまま画像に変換
し表示する段階３を続けて設けて、この画像を医学者が
観察して整理１する段階４を最後に設け構成しており、
医学者は自らの知識に従いこの画像から特徴を見出し自
らの経験により傾向を分類する様になっていた。

【００１４】

【実施例】以下、本発明の方法の実施例を図面を参照し
て説明する。図１は本発明の方法による実施例を示す原
理図で、この実施例におけるシステム５の主要部６は図
１１の従来例における表示する段階３を削除し、代わり
に、分析に必要な要素をディジタル信号処理により抽出
する段階７を先ず設け、この要素を画像処理により分析
する段階８を続けて設け構成しており、本発明人が開発
した前記の測定装置ＤＥＴ２を介して導入したデータを
準備２する段階９をこの主要部６の前に有し、この要素
について分析した結果を整理２する段階１０をこの主要
部６の後に続けて有する様になっている。

【００１５】ここで、以降の説明における理解を容易に
する為に、一般的な信号の波形に含まれる成分を分類す
る手順について述べる。図２は脳波を測定して得られた
電位の波形Ｅが変化する図である。横軸は時間ｔの推移
を示し、縦軸は電位ｅの変化を示して、両軸の交点Ｏは
時間ｔと電位ｅとが共に零である原点を示しており、電
位の波形Ｅが時間ｔと共に変化する様子を現す様になっ
ている。

【００１６】この電位の波形Ｅは、計測するそれぞれの
チャネルについて仮に次に示す数１で定義され、

【００１７】

【数１】

【００１８】即ち、時間と共に変化する関数である振幅
Ｅ_tと位相Φ_tと周波数ｆとのそれぞれの成分を情報の
源であるニューロンの活動が変調する様になっている。

【００１９】この様な電位の波形Ｅは、脳の精神活動に
よる特異な変化の存在しない平穏なコントロール状態
と、この特異な変化が生ずる活性化された状態とではこ
れら振幅Ｅ_t位相Φ_t周波数ｆの成分の値を大きく変動
し、それぞれの特有な分布に従って極値化する変化を構
成しており、脳の精神活動が変遷するに連れ、頭皮面の
いずれかの位置で時間と共にそれぞれの成分を極大化す
る様になっている。

【００２０】従って、この変遷を把握する為には、電位
の波形Ｅを構成するこれらの成分を所定の時点において
スペクトルに解析し、このスペクトルを極大化する周波
数ｆを含む帯域の信号を選別し、この信号がそれぞれの
成分を頭皮面において変動する分布を求める必要があ
る。

【００２１】図３はこの実施例の具体的なプログラムを
概略的に示す流れ図である。そこで、抽出する段階７は
導入した電位の波形Ｅを所定の時間においてそれぞれの
チャネルについて周波数ｆのスペクトルに解析する手段
１１を前置して設け、周波数ｆの帯域を全てのチャネル
について限定する手段１２を次に設けて構成しており、
それぞれのチャネルについてスペクトルを極大化する周
波数ｆの帯域を抽出する様になっている。

【００２２】続いて、分析する段階８は、限定された帯
域においてそれぞれのチャネル相互間におけるポテンシ
ャルを求めるポテンシャル分析の手段１３と、同様にし
て相関の大きさから相関を求めるコヒーレンス分析の手
段１４と、同じくバリアンスを求めるバリアンス分析の
手段１５と、同じくこの極大化する周波数ｆの近傍にお
ける周波数ｆの微小な変移の分布を求めると共に、更に
その位相Φ_tが変化する分布を求める位相分析の手段１
６とを並列に設けて構成しており、それぞれの手段を分
析者が指定して選択し得る様になっている。尚、１７は
システムを起動する段階で、１８はファイルを取り出す
段階で、１９は分析者が診断する段階で、２０はハード
・コピーをする段階で、２１はファイルを格納する段階
で、２２はシステムを終了する段階である。

【００２３】次に、本発明を実施する為の装置について
詳しく説明する。図４は実施例の装置を概略的に示すブ
ロック図である。この実施例の装置は電位の波形Ｅを導
入する入力部ＩＮと、導入された電位の波形Ｅを所定の
手順に従って処理するシステム２３と、処理して得られ
たデータをファイル化し格納するファイル・メモリＦＭ
と、このメモリから画面データを随時に取り出して表示
する表示部２４と、この画面データを印刷するプリンタ
２５と、処理する手順を設定するパラメータを指示する
キーボード２６とを設け、相互にバス接続して構成され
ている。

【００２４】システム２３は周波数ｆのスペクトルを解
析する手段１１と、周波数ｆの帯域を限定する手段１２
と、解析されたそれぞれの成分に従う分布を求めるそれ
ぞれの分析する手段１３，１４，１５，１６と、入力部
ＩＮと、表示部２４とプリンタ２５とキーボード２６と
を制御する手段Ｉ／Ｏとをプログラム・モジュール化し
て設け構成しており、指示されるコマンドを複号化しシ
ステム全体を管理するシーケンサ２７により所定の手順
に従ってそれぞれを制御する様になっている。

【００２５】ファイル・メモリＦＭは数値データをそれ
ぞれのテーブルに分類して格納する数値ファイル２８
と、イメージデータ２９をそれぞれの画面に分類して格
納する処理データ・ファイルとに構成している。

【００２６】この数値ファイルは所定の時間におけるそ
れぞれのチャネルについて解析された周波数のスペクト
ルと、限定された周波数の帯域について極大値を生ずる
周波数ｆとこの極大値に有するチャネル番号とを配列す
る極大値テーブルと、キーボードから指定するパラメー
タを配列するパラメータ・テーブルとを有して構成して
おり、処理部におけるそれぞれのプログラム・モジュー
ルによりアドレス制御してそれぞれのファイルをアクセ
スする様になっている。

【００２７】処理データ・ファイルはコヒーレンス分析
により相関の強さを求めるコヒーレンス係数の分布と、
バリアンス分析により求めるバリアンスの分布と、周波
数分析により周波数ｆの微小な変位を求める分布と、位
相Φ_t分析により位相Φ_tの変化を求める分布とを有し
て構成している。

【００２８】続いて、この実施例の方法により相関を分
析するプログラムについて図１と３及び図５を参照して
詳しく説明する。このプログラムは図１における準備２
と抽出と分析と整理２とのそれぞれの段階９，７，８，
１０をモジュールに細分化して構成しており、それぞれ
のモジュールは観察者が指定するコマンドを判定して実
行する様になっている。

【００２９】先ず、図１における準備２の段階９は本発
明のシステムを起動する手段１７を先ず設け、続いてこ
のシステム上にデータを準備するファイル・ロードの手
段１８を次に設けて構成している。起動する手段１７は
前記して述べた図４の装置を始動し、所定のメニュー画
面等を使用して本発明のシステムを起動する様になって
いる。

【００３０】ファイル・ロードの手段１８はファイル・
メモリをアクセスする処理と、電位の波形Ｅを数値化し
た新しいデータ・ファイルを本発明のシステムに導入す
る処理とにより構成しており、観察者は電位の波形Ｅを
期待する画面の様式で表示する為に、ファイルのデータ
を標本化する個数とチャネルの本数と画面上の拡大率と
をパラメータの数値に現して指定し、ファイルを導入す
るコマンドＦＬをキーボード等から指示し得る様になっ
ている。

【００３１】図６は表示する画面を示す図であり、イメ
ージを表示する図形の領域４１と、この画面で指示し得
るそれぞれのコマンドを表示するコマンド用の窓４２
と、この画面で設定するパラメータの入力を促すメッセ
ージと入力されたパラメータの数値とを表示するパラメ
ータ用の窓４３とを設けて表示するフォーマットを構成
しており、以下それぞれのモジュールの手段に使用する
画面においても同様のフォーマットを使用する様になっ
ている。

【００３２】図７は電位の波形Ｅを表示する領域を示す
図であり、図６における図形の領域４１に使用する様に
なっている。図７においてそれぞれのチャネルにおける
電位の波形Ｅは測定器で計測したチャネル番号Ｃｈ．を
縦方向に、計測した時間を横方向に配置して構成してお
り、新しいデータ・ファイルに数値化されている電位の
波形Ｅを、これらのチャネル番号Ｃｈ．に従ってそれぞ
れに配列し、変化する電位の値を波形化し時系列に表示
する様になっている。

【００３３】続いて、本発明の主要部である抽出と分析
の段階について図１と４とを参照し詳しく説明する。図
５はこの実施例のモジュールにおけるプログラムを概略
的に示す流れ図である。

【００３４】先ず、分析者が期待する分析の方式を指令
する処理３１を始めに設け、次に指令を判定する処理３
２を続いて設け、更に指令された方式に従いデータを演
算する処理３３を肯定側に続けて設け、否定側には次の
モジュールに分岐する処理３４を設け、また演算した結
果を相互に連続して補間する処理３５をこの演算する処
理に続いて設け、更に補間した数値を図形に階調する処
理３６を続けて設け、最後に階調した図形を表示する画
面に展開する処理３７を続いて設け構成している。

【００３５】従って、分析者は指令する処理においてこ
の手段の前に続く前画面３８を使用し、必要なパラメー
タを予め設定し指令するコマンドを予め指定する様にな
っており、この手段の後に続く後画面３９を次の手段に
おける指令する処理に使用する様になっている。

【００３６】図３において使用されるコマンド群を表１
により説明し、それぞれのコマンドを区別する分類と符
号とに従ってそれぞれに付された名称と指示する機能と
を説明している。

【００３７】

【表１】

【００３８】さて、前記して述べた様に抽出する段階
７は解析する手段１１と限定する手段１２とで構成して
おり、以下、抽出と分析とのそれぞれの段階７，８にお
けるそれぞれの手段において相互に異なる処理について
のみ詳しく説明する。

【００３９】先ず、解析する手段１１は指令するコマン
ドをＦＴに指定し、演算する処理３３をフーリエ変換Ｆ
Ｔに設定し、時間を追って表示するスペクトルに後画面
３９を設定してこの解析する手段１１を実行する。

【００４０】このフーリエ変換ＦＴは電位の波形Ｅを周
波数ｆの成分について解析し、スペクトルの数値を演算
して求め、スペクトルがそれぞれのチャネルにおいて分
布する画面を表示しており、それぞれのチャネルについ
て所定の時間における周波数ｆの成分が有するポテンシ
ャルＰの強度を分析者が識別する様になっている。

【００４１】ここで、ｊ番目のチャネルから導入した電
位の波形Ｅを時間の関数Ｕj(t)とすると、ポテンシャル
Ｐの大きさは次の数２の様に定義される。

【００４２】

【数２】

【００４３】ポテンシャルＰの強度の関数Ｐj(f)は電位
の波形Ｅの関数Ｕj(t)を所定の時間間隔Ｔについてフー
リエ変換ＦＴし、次にこの結果の２乗平均値を求めてお
り、それぞれのチャネルにおけるポテンシャルＰの強度
が正負の符号に拘らず周波数ｆの成分に従う変化をこの
時間間隔Ｔで分割するスペクトルになる様になってい
る。

【００４４】図６におけるこのスペクトルを表示する図
形４４は、全てのチャネルに対応するそれぞれのスペク
トルを重ねて合成し表示している。このスペクトルを表
示する図形４４において、分析者が周波数ｆの上限と加
減とを指定して限定する帯域を設定し、フィルタリング
ＦＬＴのコマンドをキーボードから指示し得る様になっ
ている。

【００４５】続いて、限定する手段は前記した図５にお
いてコマンドをＦＬＴに指定し、演算する処理をフィル
タリング演算に設定し、周波数ｆの変化を分布する画面
を後画面に設定してこの手段を実行する様になってい
る。

【００４６】このフィルタリング演算は合成する図形に
おいて設定した周波数ｆの帯域について極大値を形成す
る周波数ｆをそれぞれのチャネル毎に演算して求め、周
波数ｆを分布する図形を表示しており、周波数ｆが極大
値を有する任意の帯域を分析者がキーボード等から設定
し得る様になっている。

【００４７】この際に、周波数ｆの変化を分布する図形
も選択して表示し得る様にしてあり、極大値を有する周
波数ｆが実際には頭皮面で微妙に変化する様子を、前記
した配置する図形を使用し分析者が識別する様になって
いる。

【００４８】続いて、前記して述べた分析の段階はポテ
ンシャル分析する手段とコヒーレンス分析する手段とバ
リアンス分析する手段と位相分析する手段とを相互に並
列に設けて構成しており、分析者が指定するコマンドに
従っていずれかひとつの手段を選択しこの分析する段階
を実行する様になっている。ポテンシャル分析する手段
はコマンドにＰＳＤを指定し、演算する処理をポテンシ
ャルＰの強度演算に設定し、このポテンシャルＰの強度
を変化する画面を後画面に設定してこの手段を実行する
様になっている。このポテンシャルＰの強度演算は前記
した周波数ｆの帯域において極大値を有するそれぞれの
チャネルのポテンシャルＰの強度を演算して求め、この
ポテンシャルＰの強度の変化を分布する図を表示してい
る。図８はポテンシャルＰの強度と、結合の係数と、バ
リアンスとのそれぞれの変化を分布するフォーマットを
共通に示す図である。

【００４９】コヒーレンス分析する手段はコマンドにＣ
ＣＤを指定し、演算する処理を結合の係数演算に設定
し、この結合の係数を変化する画面を後画面に設定して
この手段を実行する様になっている。この結合の係数演
算は前記した周波数ｆの帯域において極大値を有するそ
れぞれのチャネル相互における相関の強さをそれぞれの
チャネル毎に係数として演算して求め、この係数の変化
を分布する図を前記した図８のフォーマットで表示して
おり、設定された帯域で極大値を有する位置が頭皮面で
相互に相関する様子を、前記した配置する図形を使用し
分析者が識別する様になっている。

【００５０】バリアンス分析する手段はコマンドにＭＶ
Ｄを指定し、演算する処理をバリアンス演算に設定し、
バリアンスの変化を分布する画面を後画面に設定してこ
の手段を実行する様になっている。このバリアンス演算
は前記した周波数ｆの帯域において極大値を包含する波
形の面積をそれぞれのチャネル毎に演算して求め、結果
の値を分布する図を前記した図８のフォーマットで表示
しており、頭皮面で変化する様子を前記した配置する図
形を使用し分析者が識別する様になっている。

【００５１】位相分析する手段はコマンドにＰＦＤを指
定し、演算する処理を位相演算に設定し、位相Φ_tの変
化を分布する画面を後画面に設定してこの手段を実行す
る様になっている。この位相演算は前記した合成する図
形において極大値を形成する特定の周波数ｆを指定し、
フーリエ変換で派生した位相の成分をそれぞれのチャネ
ル毎に基準のチャネルの位相Φ_tからの差として演算し
て求め、位相Φ_tの変化を分布する図形を表示してい
る。

【００５２】図９は位相Φ_tの変化を分布する図形であ
り、極大値を有する周波数ｆの位相Φ_tが頭皮面で変化
する様子を、前記した配置する図形を使用し分析者が識
別する様になっている。尚、図１０はそれぞれのチャネ
ルが担当する頭皮面の位置を配置する図形に縮小したス
ペクトルを分散して表示する分布の図形である。

【００５３】最後に、図１において述べた様に整理２の
段階１０は前記した分析の段階８に続いて、分析した結
果を分析者が診断する手段１９と、このシステム上のデ
ータをハード・コピーする手段２０と、同じく保管する
ファイル・セーブの手段２１とを次に設け、このシステ
ムを停止する終了の手段２２を最後に設け構成してい
る。終了の手段は所定のメニュー画面等を使用してこの
システムを終了し前記した図４の装置を停止する様にな
っている。

【００５４】ファイル・セーブの手段はファイル・メモ
リをアクセスする処理と、分析して作成したデータ・フ
ァイルをこのメモリに格納する処理とにより構成してお
り、観察者は作成したデータ・ファイルを期待するファ
イルの様式で格納する為に、ファイルの名称と格納先の
番地等をパラメータの数値に現して指定し、ファイルを
格納するコマンドＳＶをキーボード等から指示し得る様
になっている。

【００５５】尚、本発明の方法は前述にのみ限定される
ものではなく、本発明の要旨を逸脱しない範囲で種々の
変更を加えうることは勿論である。

【００５６】

【発明の効果】頭皮面における計測した電位の波形から
極値を形成する周波数を有する信号をディジタル信号処
理により抽出するので、多数のチャネルを介して連続的
に収集する多量のデータを処理することができる。ま
た、抽出したこの特徴的な信号が分析する相関を統計学
的な演算処理により時間を追って分析するので、脳波が
変遷する傾向を精確に把握し時間的に識別することがで
きる。更に、相関分析した結果を画像処理により分布図
に表示するので、視覚により総合的に診断することがで
き、再現性もあり信頼することができる。

【図面の簡単な説明】

【図１】本発明の方法による実施例を示す図である。

【図２】脳波による電位の波形を示す図である。

【図３】図１における具体的なプログラムを概略的に示
す流れ図である。

【図４】図１における具体的な装置を概略的に示すブロ
ック図である。

【図５】図３におけるそれぞれの分析の手段を概略的に
示す共通の流れ図である。

【図６】図５における前及び後画面の表示する様式を示
す説明図である。

【図７】図６におけるイメージの領域に脳波による電位
の波形を多チャネルで表示する説明図である。

【図８】図６におけるイメージの領域にスペクトルの相
関を表示する図形を展開する分布図である。

【図９】図６における頭皮面の配置図に図５において位
相を表示する図形を展開する分布図である。

【図１０】図６における頭皮面の配置図に図５において
スペクトルを表示する図形を展開する分布図である。

【図１１】従来の脳波の分析方法を説明する図である。

【符号の説明】

Ｓ１測定器ＤＥＴ２計測装置９準備２する段階７抽出する段階８分析する段階１０整理２する段階１３ポテンシャル分析する手段１４コヒーレンス分析する手段１５バリアンス分析する手段１６位相分析する手段

Claims

【特許請求の範囲】

【請求項１】脳内における活動ニューロンの電気的状
態と、人間の意識状態とを対応づける相関を分析する方
法であって、頭皮面において計測した電位の波形から、特徴となる要
素を抽出する段階を先ず設け、前記の要素が位置的に分布する相関を、時間的に分析す
る段階を続けて設け構成しており、前記の相関を時間的な分布図に表示し、前記の分布図を
格納し得る様になっていることを特徴とする脳波の相関
分析方法。
【請求項２】前記の特徴となる要素はスペクトルにお
いて極値を形成する周波数の近傍である請求項１に記載
の脳波の相関分析方法。
【請求項３】前記の特徴となる要素は前記のスペクト
ルを有する信号のポテンシャル及び位相である請求項１
及び２に記載の脳波の相関分析方法。
【請求項４】前記の分析する段階は位置的に結合する
コヒーレンスの係数を時間を追って分析する手段である
請求項１ないし３に記載の脳波の相関分析方法。
【請求項５】前記の分析する段階は前記のスペクトル
を時間を追ってバリアンス分析する手段である請求項１
ないし３に記載の脳波の相関分析方法。