(57)【特許請求の範囲】
【請求項１】
コンピュータによって実行される、生理学的現象を示す非定常信号を、基本成分としてＱＧＫ（準ガウシアンカーネル）を用いて機能要素のセットへと分解する方法であって、各々の機能要素を特徴づけるパラメータを推定するために、以下の工程を含む方法：
信号の時間経過における動的トレンド及び基線トレンドを推定するため、マルチパス移動平均フィルタを適用し、トレンドを用いて、機能要素が時間的に配列されている信号セグメントを示すガイド関数を抽出する工程（当該工程において、ゼロ交差点並びにガイド関数の絶対値の大域的な最小値及び局所的な最小値の各々を分割点の配列の１つとして定義するように、前記ガイド関数は、分割点の配列の定義によって信号断片へと分割される）；
上記非定常信号をガイド関数の分割に基づいて断片に分け、上記非定常信号を、上記ガイド関数の分割点に対応する点において分割する工程；
先行する機能要素に適合するＱＧＫからの重複を補正し、ＥＣＫ（経験的成分カーネル）を形成する工程であって、再帰的なサブ工程は以下を含む工程：
相似基底関数（ＳＢＦ）アルゴリズムを用いて、周波数領域へＥＣＫを変換する工程；
ＥＣＫの正規化振幅スペクトルとＴ（ｖ）＝ｅｘｐ（−ｖ２）で解析的に定義されるガウシアンプロファイル周波数領域のテンプレートとの間の最善の適合を見出し、評価することによって、κ（加重）及びσ（分散）パラメータを推定し、かつ、Ｅ、Ｇ、Ｓ、Ｗ及びＵの分類階級を用いたパラメータ信頼性を検証する工程；
位相スペクトルの定義された値に対応する特性点を見出すことによって、ＥＣＫの複素スペクトルの虚部から、特定のκ及びσパラメータと共に、ＱＧＫの開始時間τを推定する工程；
新たに特定したＱＧＫを加えることによって、再帰サイクル中の先行するサイクルにおいて特定されたＱＧＫの和から構成された、非定常信号モデルを拡張する工程；及び
特定結果の評価のため、分類パラメータξを用いて、不都合なＱＧＫを除去した後、特定された要素を再構成する工程；
ＱＧＫの集合の選択によって、非定常信号の部分モデルを作成する工程であって、上記ＱＧＫの集合のパラメータは予め定義した階級に属する工程。
【請求項２】
以下に従って分類パラメータξは、１から５の整数値を取る、請求項１に記載の方法：
【数１】
0006483107.tif 000209
式中、Ｅ−優れている、Ｇ−よい、Ｓ−満足できる、Ｗ−不十分、及びＵ−不成功である。
【請求項３】
非定常信号にローパス無限インパルス応答フィルタを適用することによって、動的トレンドを推定する、請求項１に記載の方法。
【請求項４】
上記ガイド関数の信号セグメントへの分割は、ゼロ交差点のみに基づいた分割点により定義される、請求項１又は３に記載の方法。
【請求項５】
各々のＥＣＫに適合する、各々のＱＧＫの初期値は、当該各々のＥＣＫの初期値における分割点と一致している、請求項１に記載の方法。
【請求項６】
推定工程は、ＥＣＫの複素スペクトルの位相及びガウシアン複素スペクトルの位相の間において、最善の適合を見出すことによって、各ＥＣＫに適合する各ＱＧＫの修正された初期値を推定する、サブ工程をさらに含む、請求項５に記載の方法。
【発明の詳細な説明】
【発明の詳細な説明】
【０００１】
〔技術分野〕
本発明は、全体的に、信号処理、及び、特に非定常生物医学的信号の機能要素への分解、及び、上記機能要素の各々を特徴づけるパラメータの推定に関する。
【０００２】
〔背景技術〕
非定常信号は、時間に依存する特性値を伴う系によって発生する。理論的立場、実験的立場及び診断的立場から上記非定常信号を理解し、分析する必要性が、生物医学的研究のほぼすべての分野において、生じている。しかし、非定常信号は、一般信号理論に従わず、また、非定常信号を扱う好ましい数学的手順は、かなり限られている。
【０００３】
現在の信号処理手順は、一般的に、線形定常系の理論に基づいている。従って、実データの非定常特性は、しばしば無視されるか、又は上記効果は無視できるものとみなされる。時間的不変性及び線形性に関する当該仮説によって、いくつかの応用例に関して許容できる結果が提供されているが、多くの生理現象は、極めて可変的及び非線形的であり、従って、理にかなった正確な線形近似が認められていない。一方で、膨大な数の科学的及び診断的証拠により、心電図（ＥＣＧ）及び脳波図（ＥＥＧ）のような、主要な非定常生理信号に関する、より総合的な解析は、異なる正常条件及び病的条件下における、様々な生理学的システムのメカニズム及び状態に関する、より深い理解を導き得ることが示唆される。このような背景から、ヒトの脳及び体のイメージング、コンピュータを導入した健康モニタリング、並びにＥＥＧベースの脳−コンピュータ間インターフェースといった、ますます挑戦的な多数の応用例への実現可能なアプローチとして、非定常信号処理に関する新たなアルゴリズムの開発が浮上している。同時に、コンピュータを用いた生物医学的システムの継続的な進歩によって、速さ、サイズ及び費用といった点で、洗練されたアルゴリズムの開発が、実用的かつ優れたコストパフォーマンスを示すようになってきた。
【０００４】
非定常信号解析は、線形定常系理論の場合のような適切な理論及び計算の枠組みによって裏付けられていない。特定のタイプの生理信号の解析においては、数学的に厳格な方法よりもむしろ発見的方法が、一般的に受け入れられている。例えば、ＥＥＧ信号解析の際、上記信号は、局所定常性概念（Ｂａｒｌｏｗ，  Ｊ．Ｓ．  “Ｍｅｔｈｏｄｓ  ｏｆ  ａｎａｌｙｓｉｓ  ｏｆ  ｎｏｎ−ｓｔａｔｉｏｎａｒｙ  ＥＥＧｓ，  ｗｉｔｈ  ｅｍｐｈａｓｉｓ  ｏｎ  ｓｅｇｍｅｎｔａｔｉｏｎ  ｔｅｃｈｎｉｑｕｅｓ：  ａ  ｃｏｍｐａｒａｔｉｖｅ  ｒｅｖｉｅｗ”．  Ｊ．  Ｃｌｉｎ．  Ｎｅｕｒｏｐｈｙｓｉｏｌ．  １９８５，  ｖｏｌ．  ２，  ｐｐ．  ２６７−３０４  参照）を用いて解析され得る。上記概念は、信号は本質的には非定常であるが、極少時間内において、信号パターンがごくわずかに定常でなくなることを推測している（Ｐｒｉｅｓｔｌｅｙ，  Ｍ．Ｂ．：  “Ｎｏｎ−Ｌｉｎｅａｒ  ａｎｄ  Ｎｏｎ−Ｓｔａｔｉｏｎａｒｙ  Ｔｉｍｅ  Ｓｅｒｉｅｓ  Ａｎａｌｙｓｉｓ”．  Ａｃａｄｅｍｉｃ  Ｐｒｅｓｓ，  Ｌｏｎｄｏｎ，  １９８８  参照）。このような準定常成分は、明確な非線形系によって発生しうる信号の種々の要素と結びついている。
【０００５】
これまで、非定常信号の成分の検出に関する問題点には、データ駆動型信号処理方法論及びモデル駆動型信号処理方法論を用いて対処してきた。
【０００６】
データ駆動型信号処理方法の利点は、複雑かつ実験的な波形の動きを詳細に記述できる点である。上記データ駆動型方法の一例は、Ｈｕａｎｇらによって構築されてきた、経験的モード分解（ＥＭＤ）であり（Ｈｕａｎｇ  ｅｔ  ａｌ．：  “Ｔｈｅ  ｅｍｐｉｒｉｃａｌ  ｍｏｄｅ  ｄｅｃｏｍｐｏｓｉｔｉｏｎ  ａｎｄ  ｔｈｅ  Ｈｉｌｂｅｒｔ  ｓｐｅｃｔｒｕｍ  ｆｏｒ  ｎｏｎｌｉｎｅａｒ  ａｎｄ  ｎｏｎ−ｓｔａｔｉｏｎａｒｙ  ｔｉｍｅ  ｓｅｒｉｅｓ  ａｎａｌｙｓｉｓ”．  Ｐｒｏｃ．  Ｒ．  Ｓｏｃ．  Ｌｏｎｄｏｎ，  Ｓｅｒ．  Ａ，  １９９８，  ｖｏｌ．  ４５４，  ｐｐ．  ９０３−９９５  参照）、適用原理の下、内因性モード関数と呼ばれる、「振る舞いがよい」振動関数の和に非定常信号を分解する。以下の開発において、各々内因性モード関数に対する瞬時周波数を提供するＨｉｌｂｅｒｔ変換によって、ＥＭＤは、裏付けられている。これらの測定を用いて、本技術は、瞬時周波数の変化という点に関して本システムを記述する。しかし、本解決策は、多くの実験的パラメータに依存し、厳格な基準によってその選択が立証されていないため、一意的に定義されていない。ＥＭＤに関する強固な数学的基盤なしでは、この純粋に経験的及び直観的なアプローチは、多くの難解な方法論的課題を生み出す。例えば、内因性モード関数は、解析的定義なしで、反復アルゴリズムの出力として単に数字的に定義されてしまう。これらの条件下で、方法論の異なる局面についての判断は、経験論に基づく状況に応じた比較によって、行われている。
【０００７】
生物医学的信号処理に対するモデル駆動型方法には、関与する系に対する正確なモデルが必要である。課題の具体的な局面は、生理信号が、複数の微視的スケールの情報源によって導かれた、大域的スケールの変数である点である。物理的モデルの作成には、関連する微視的スケールの成分の厳格な特定及び合成、これら根拠に基づいた大域的モデルの作成が必要となる。しかし、関与する細胞現象及び分子現象の独特なパターンを正確に記載することが不可能であるため、上記の厳格な決定論的概念を用いた大域的スケールから微視的スケールへの変換は、一意的な解を有さない。このような不確定性を考えた場合、大域的スケールの信号に着目した細胞機構のモデル作成に関する膨大な試みは、モデルに関する生理学的詳細及び解剖学的詳細だけでなく、基礎数学的手順という点で、大きく異なる顕著に多様な発見的手順によって、裏付けられている。本モデルの矛盾の程度は知られていない。
【０００８】
従って、従来の方法論では成功していない、非定常時系列信号の成分を効果的に特定する必要性は未だ存在しており、それゆえ、当該モデルは、信号の特性値を分類づけるのに効果的に利用することができるだろう。
【０００９】
〔概要〕
本発明の目的は、既存の手順の一つ又は複数の欠点を、実質的に克服、少なくとも、改善する点にある。
【００１０】
本発明の第一の態様において、基本成分として準ガウシアンカーネル（ＱＧＫ）を用いて自己相似機能要素のセットへと非定常信号を分解する方法であって、個々の機能要素を特徴づけるパラメータを推定するために、以下の工程を含む方法を提供する：
非定常信号の時間経過における動的トレンド及び基線のトレンドを推定し、ガイド関数及びその適合した分割を抽出する工程；
非定常信号を要素の発達（development）の領域に分ける工程；
先行する機能要素に適合するＱＧＫからの重複を再帰的に補正し、再帰サイクルに含まれる各サイクルにおいて経験的成分カーネル（ＥＣＫ）を形成する工程、そのサブ工程は以下を含む：
周波数領域へＥＣＫを変換する工程；
ＱＧＫのパラメータの推定及びそれらの検証を行う工程；及び、
非定常信号の大域的モデル及び部分的モデルを作成する工程。
【００１１】
本発明の他の態様においては、前記方法を実施するための装置を提供する。
【００１２】
本発明に係る他の態様も開示される。
【００１３】
〔図面の簡単な説明〕
本発明の一つ又は複数の実施形態は、以下の図面を参照して記載されるであろう：
図１Ａは、頭皮の表面に取り付けた二つのセンサーの間で測定される電位の典型的な波形を示す；
図１Ｂは、ＥＥＧ信号が示す極めて非定常的な特徴を例示する；
図２Ａは、聴覚オッドボールパラダイムにおいて検出される典型的なＥＲＰ信号を示す；
図２Ｂは、図２Ａに示した各セグメントを、そのセグメントに適合するパラメータを有する準ガウシアンカーネルを用いて置換した典型的な結果を示す；
図２Ｃは、準ガウシアンカーネル成分の和として得られたモデルを、比較のための図２Ａに示した元のＥＲＰ信号とともに示す；
図３は、生理信号を機能要素に適応的に分解し、その後、各機能要素を特徴づけるパラメータを決定する、波形分析方法の概略フロー図を示す；
図４Ａは、図３に記載の方法中で応用した、動的トレンド及び基線のトレンドの推定の概略フロー図を示す；
図４Ｂは、典型的な生理信号、マルチパスの移動窓による平均化を離散的な生理信号に適用して得られた基線のトレンド、及びトレンドを除去した基線の信号を示す；
図５Ａ及び図５Ｂは、ゼロ交差点から結果的に得られた分割点、並びに大域的及び局所的な最小値での分割点を各々、示す；
図５Ｃは、ゼロ交差点（実線の縦線）と最小値（点線の縦線）とにおいて、それぞれ、ＥＥＧ信号を分割した例を示す；
図６Ａは、高い時間分解能を持つ分割点を提供する初期信号及びガイド関数を示す；
図６Ｂ及び図６Ｃは、図６Ａに記載のガイド関数及び初期信号を、縦線で示した分割点とともに再度示す；
図７Ａは、図６Ａに記載の初期信号、及び低い時間分解能による信号分割を提供するガイド関数を示す；
図７Ｂ及び図７Ｃは、図７Ａに記載のガイド関数及び初期関数を、縦線で示す分割点とともに再度示す；
図８Ａは、準ガウシアンカーネルの構成に用いたガウシアン成分を示す；
図８Ｂは、準ガウシアンカーネルの一例を示す；
図９Ａは、有限成分の和への挿入箇所の分解を例示する；
図９Ｂは、三つの三角形基底関数の和へのハット関数の分解を示す；
図９Ｃは、同様の基底関数アルゴリズム中で使用した、三角形基底関数を示す；
図９Ｄは、三角形基底関数の複素スペクトルの実部と虚部とを示す；
図１０Ａは、パラメータσの推定に用いた振幅テンプレートを示す；
図１０Ｂは、パラメータσの推定過程における、正規化した実験振幅スペクトル、並びに異なる位置の振幅テンプレート、Ｔ１、Ｔ２及びＴ３を例示する；
図１１Ａは、階級Ｅに属する正規化振幅スペクトルに、テンプレートを重ね合わせた結果を示す；
図１１Ｂは、階級Ｇに属する正規化振幅スペクトルに、テンプレートを重ね合わせた結果を示す；
図１２は、図中央のパネルに記載の、ＥＥＧ信号の異なるセグメントの正規化振幅スペクトルに重ね合わせた振幅テンプレートのいくつかの結果を例示する；
図１３Ａは、低い時間分解能を用いたＥＥＧ分割の例を示す；
図１３Ｂは、図１３Ａに記載の１０セグメントの振幅スペクトルを示す；
図１３Ｃは、振幅と周波数との両方を正規化した後、図１３Ｂの振幅スペクトルを重ね合わせたものを示す；
図１４Ａは、パラメータτの推定値が得られることを基準にした場合の横座標周波数スケール上の点ａ及びｂを示す；
図１４Ｂは、横座標スケール上の点ｂが、複素スペクトルの虚部の最小値によって定義された場合を例示する；
図１５Ａは、話している人の呼吸の波形例を示す；
図１５Ｂは、図１５Ａの信号において特定した５４個のＱＧＫの波形を示す；
図１５Ｃは、特定した５４個のＱＧＫの和として合成されたモデルによって、図１５Ａの信号を再現した場合の、正確性を示す。
【００１４】
〔発明の詳細な説明〕
工程及び／又は特徴に係る添付の図面の一つ又は複数において、参照がなされ、同一の参照番号が含まれるが、逆の意図が示されていない限り、これらの工程及び／又は特徴は、本説明の目的のため、同一の機能又は作用を有している。
【００１５】
非定常信号を記載するため、重要な工程がモデルを構築している。この課題に対する従来アプローチは、「物理学に基づいた」モデルによって裏付けられている：その目的は、発生源の明白な物理モデルという点で、信号を記述することである（Ｃａｎｄｙ  Ｊ．Ｖ．  “Ｍｏｄｅｌ−ｂａｓｅｄ  ｓｉｇｎａｌ  ｐｒｏｃｅｓｓｉｎｇ”．  Ｊｏｈｎ  Ｗｉｌｅｙ  ＆  Ｓｏｎｓ，  Ｎｅｗ  Ｊｅｒｓｅｙ，  ２００６  参照）。しかし、非定常信号は、通常、複数の非線形的発生源から発生し、これらの相関性は未だ明確でない。
【００１６】
これらの物理学に基づいたアプローチとは異なり、本明細書に開示の方法は、カオスに基づいたモデリング（ＣＢＭ）アプローチによって裏付けられる。これらは、発生源の振る舞いを定義する決定論的要素及び確率論的要素の両方を組み合わせている。このアプローチの一般的なフレームワークは、生物医学的信号の発生の基礎を成す一般的なメカニズムに関する文献によって紹介されている。
【００１７】
生物医学的信号解析に関する一般的な側面は、生理学的プロセスは、巨視的スケールで振る舞う複数の初期発生源によって生成される大域的スケール変数である、という点である。この観点において、生理学的プロセスは、分子的要素及び細胞的要素の非常に複雑な系の外部に対するパフォーマンス、つまり、決定論的要素及び確率論的要素で定義されるパフォーマンスを反映している。必要とされる機能要素数は非常に多く、細胞機構の詳細は非常に複雑であるため、該当する細胞全体のモデルは実現不可能である、という事実により、大きな困難が生じている。しかし、複数の微視的スケール変数によって生成される集団効果に関する、数多くの理論的研究及び計算機科学的研究が、大域的スケールにおいて重要である細胞機構に関する一般的側面に、そうでないものから、光を投じた。よって、局所的な微視的スケールのボリュームにおける、多くの初期発生源の同期的活性化は、大域的スケールにおいて測定される信号における検出可能な変化を生成するための必要条件である、というのが一般的な物理学的概念である。これらの処理に関する主要な側面は、脳波図（ＥＥＧ）、心電図（ＥＣＧ）及び筋電図（ＥＭＧ）のような、臨床的に重要な生物医学的信号に関して、幅広く研究がなされてきた。明確な生理学的要素の指標として、ピークを持つ大域的スケールの波形を示唆する一般的なフレームワークは、ＥＥＧの生成の基礎を成すメカニズムの参照によって示され得る。
【００１８】
図１Ａ中に示すように、ＥＥＧは頭皮の表面に取り付けたセンサーを用いて測定され得る。電位は、頭皮下の脳の皮質領域内のニューロンの活性によって引き起こされる。単一ニューロンの活性は、十分高くないため、頭皮の表面に取り付けたセンサー間で測定可能な電位を引き起こすことができない。むしろ、接近して存在する皮質ニューロンの大きな集団が同期して作用すると、活性の和によって、細胞外双極子が形成され、その結果、測定電位が得られる。
【００１９】
頭皮の表面で記録されるこれらの集合電位の変動及びダイナミクスは、タイムラグなしで皮質の神経処理装置の活性に従う。この観点において、ＥＥＧ信号は、ヒトの脳における情報処理の基本メカニズムをリアルタイムに研究するための、最も重要な客観的指標の一つとしての機能を果たしている。図１Ｂ中に示すＥＥＧ信号の一例は、脳の電気的活性は、典型的に極めて不規則かつ非定常的であり、様々な活性パターンを生成することを示す。
【００２０】
ＥＥＧ信号という観点で情報処理を理解する鍵は、ＥＥＧ信号中の変化及び認知的課題中の特定刺激との時間関係を検出することによって提供される。ＥＥＧ信号のこのような特定の断片は、事象関連電位（ＥＲＰｓ）と呼ばれている。図２Ａは、聴覚オッドボールパラダイムにおいて誘発される典型的なＥＲＰを示す。オッドボールパラダイムは、２つのタイプの音程、すなわち、バックグラウンド刺激及びターゲット刺激を用いてＥＥＧ信号の特徴に影響を与える試験である。図２Ａ中に示すＥＥＧ信号において、ｔ＝０において、ターゲット刺激を与える。この状況の後に、連続した正負の波形の揺れが生じる。重要な実験的及び臨床的証拠を元に、刺激を与えてから特定の時間間隔において観察される波形のピークは、様々な皮質ニューロンネットワークに関連する機能エンティティとして見なされている。聴覚オッドボールパラダイムにおいて、このようなピーク電位は、各々Ｎ１００、Ｐ２００、Ｎ２００及びＰ３００と示されている後期潜時(late latency)ＥＲＰとして知られている。数字は、刺激後の負のピーク（Ｎ）又は正のピーク（Ｐ）の出現のおおよその時間（単位はミリ秒）を表している。Ｐ３００の特定の側面は、反復ＥＲＰの平均を検討する際、後期潜時要素は、巨大な一つのピークとして現れる点である。しかし、単一のＥＲＰ波形、すなわち、元のＥＥＧ信号を検討する際、後期潜時ＥＲＰであるＰ３００は、通常、Ｐ３ａ及びＰ３ｂと示されている２つの補要素からなる複雑な電位を示し得る。
【００２１】
上記様式で名づけられる異なる正負の変動は、重要な生理学的システムに密接に関連する一連の活性化処理の客観的指標として現れる。特定の構造機能機構を有する機能系によって生成される要素の指標としてのピークを持つ波形の概念は、科学的及び臨床的に非常に重要なほぼすべての生理学的プロセスへと従来応用されている。広く受け入れられたアプローチを考えると、種々の異なる手動又はコンピュータを導入した測定技術（ピーク感知手段）は、生理信号をピークの振幅及びピークの潜時へと変換する。これら広く実験的な技術の致命的な限界は、離れた時間点におけるピークの振幅及び潜時への信号の変換によって、重要な機能的情報及び臨床的情報を含んだ生物学的波形のダイナミクスを特定できない点である。計算的観点において、これらすべての技術の欠点は、要素の時間的重複を分解できない点である。重複を修正しなければ、振幅及び潜時のパラメータ測定は、非常に不確定なエラーの影響を受け得る。
【００２２】
図３は、単なる波形のピークではなく、特定の波形を持つ（正又は負の）振れとして要素を捉える波形解析の方法１００の概念フロー図を示す。より具体的には、方法１００は、非定常信号ｖ（ｔ）を機能要素に適応的に分解し、その後、時間−周波数解析及びパラメータ推定手順を適用し、動的モデル及び関連パラメータという形で、各々の機能要素を特定する。非定常信号ｖ（ｔ）は、ＥＥＧ信号、心電図（ＥＣＧ）信号、筋電図（ＥＭＧ）信号、呼吸由来の信号等であり得る。
【００２３】
コンピュータシステム内で実行される、一つ又は複数のアプリケーションプログラムのようなソフトウェアを用いて、方法１００を実施し得る。代わりに、一つ又は複数の集積回路のような方法１００の機能又は副次機能を果たす、専用ハードウェアにおいて、方法１００を実行してもよい。このような専用ハードウェアは、デジタル信号処理装置、又は一つもしくは複数のマイクロプロセッサ及び関連メモリを含み得る。
【００２４】
方法１００は、工程１０１から始まり、アナログ信号である非定常信号ｖ（ｔ）を調整し、標本値ｖｍ＝ｖ（ｔｍ）からなる離散入力数列とする。より具体的には、非定常信号ｖ（ｔ）を標準電圧の範囲内で増幅し、一定のサンプリング間隔ｔｍ＝ｍΔ（ｍは整数値であり、Δはサンプリング間隔である）において抽出する。この手順により、標本セット
【００２５】
【数１】
0006483107.tif 000002
【００２６】
が形成される。非定常信号ｖ（ｔ）は、帯域制限信号として見なすことができ、サンプリング間隔Δは、抽出された非定常信号ｖ（ｔ）の正確な再現に必要な最小のサンプリングレートである、ナイキストレートを超えるように選択する。
【００２７】
詳細は以下に記載するが、その後、方法１００は、非定常信号ｖ（ｔ）の離散的表現である、離散入力数列
【００２８】
【数２】
0006483107.tif 000003
【００２９】
上で作用し、入力数列
【００３０】
【数３】
0006483107.tif 000004
【００３１】
をセグメントに分割し、各々のセグメントに時間−周波数解析を適用し、セグメントに対応するパラメータと動的モデルとを、各々のセグメントにおいて特定する。パラメータ数列は以上のように形成され、パラメータ数列は、使用者が定義した基準に応じた各々のセグメントを分類し、非定常信号ｖ（ｔ）又は定義されたセグメントを合成するために使用される。
【００３２】
よって、非定常信号ｖ（ｔ）、及び離散対応値
【００３３】
【数４】
0006483107.tif 000005
【００３４】
は、典型的に高い周波数のノイズ及び低い周波数のバックグラウンド活性を含んでいる。これらの不適切な活性は、通常、ローパス及びハイパスフィルタを用いて、更なる解析から取り除かれる。しかし、重要な機能要素及び不適切な活性の両方の、非常に不規則及び非定常的な性質が原因で、時間不変的周波数の特徴を持つフィルタの適用は、機能要素のダイナミクスを、著しくかつ予期せず歪ませ、機能的に重要な情報が失われ得る。
【００３５】
この課題に対する一般的な解決策は、適応フィルタの使用である。このような適応フィルタは、信号のダイナミクスにおいてトレンドを変えるために動作を適応させなくてはならない。最も一般的にトレンド推定に用いられるデジタルフィルタは、移動平均フィルタ又はボックスカーフィルタであり、標本各々を、問題になっている標本の平均及び任意の数の連続する前後の点に、置換する有限インパルス応答を用いたデジタルフィルタである。しかし、ボックスカーフィルタの周波数応答は、最適ではなく、単純な移動平均は、トレンド軌道において著しい歪みの原因となる。このような歪みを防ぐため、工程１０２では、入力数列ｖから、高い周波数のノイズ及び低い周波数のバックグラウンド活性を省くため、マルチパス移動平均（ＭＰＭＡ）フィルタを用いる。より具体的には、２つのＭＰＭＡフィルタは、入力数列ｖにおける動的及び基線トレンド推定のため、入力数列ｖに適用される。
【００３６】
ＭＰＭＡは、反復処理であり、入力数列に移動窓平均を適用し、一回目の反復の入力数列は入力数列ｖであり、それに続く次の反復の入力数列は先行の反復の結果となる。従って、ｋ回目の反復において、各々、入力数列は
【００３７】
【数５】
0006483107.tif 000006
【００３８】
であり、出力数列は
【００３９】
【数６】
0006483107.tif 000007
【００４０】
であり、式中、標本数はＭ＋１、ｋは１からＫの値を取る反復数であり、
【００４１】
【数７】
0006483107.tif 000008
【００４２】
である。
【００４３】
ｍ個目の標本に対する出力数列
【００４４】
【数８】
0006483107.tif 000009
【００４５】
は、以下のように与えられる：
【００４６】
【数９】
0006483107.tif 000010
【００４７】
式中、Ｐは、窓の広さを定義し、当該窓の広さはＮＷ＝２Ｐ＋１である。
【００４８】
Ｐの連続したその前の点及びＰの連続したその後の点によって
【００４９】
【数１０】
0006483107.tif 000011
【００５０】
が定義されることは、方程式（１）より明らかである。しかし、数列ｖの定義によると、問題となっている処理の第一の標本は、
【００５１】
【数１１】
0006483107.tif 000012
【００５２】
である。結果として、ｍ＋ｊ＜０である値
【００５３】
【数１２】
0006483107.tif 000013
【００５４】
は、計算から省かれる。これは、ｎが−１から−Ｐの値を取るｖｎ＝０のゼロ値、及び拡張数列ｖからなるゼロパディングとしてみなされ得る。
【００５５】
従来の移動窓平均によって、結果的に得られた数列のスペクトル中の高周波数の歪みが生成される一方で、工程１０２において適用されるＭＰＭＡは、スペクトルから徐々にこのようなずれを除去する。ｍがＰからＭ−Ｐの値を取る残差
【００５６】
【数１３】
0006483107.tif 000014
【００５７】
は、各反復に従って小さくなる。ｋ回反復後のトレンド推定における改善は、反復間の標本値変化を測定することで評価する。正確性は、以下の平均二乗誤差で表す：
【００５８】
【数１４】
0006483107.tif 000015
【００５９】
反復は、所望の改善が達成されると終了する。ＥＥＧ、ＥＲＰ、ＥＭＧ、及びＥＣＧ等の、典型的な電気生理信号に適用する場合は、以下の条件を満たす場合、ｋ回反復後に、反復は終了する：
【００６０】
【数１５】
0006483107.tif 000016
【００６１】
式中、Ｍｖは出力信号ｖｋの平均値であり、Ｔｒは、無次元の閾値である。これら適用において、０．０５が妥当な閾値Ｔｒである場合、反復は、典型的に５反復後に終了する。
【００６２】
図４Ａは、工程１０２中に適用する動的及び基線のトレンド推定の概要フロー図を示す。異なる２つのＭＰＭＡは、フィルタリングされた数列を形成するため、ＭＰＭＡ１ユニット及びＭＰＭＡ２ユニット中で入力数列ｖに適用される：動的トレンドは
【００６３】
【数１６】
0006483107.tif 000017
【００６４】
であり、基線のトレンドは
【００６５】
【数１７】
0006483107.tif 000018
【００６６】
である。ＭＰＭＡ１及びＭＰＭＡ２は、異なる窓パラメータＮＷ１及びＮＷ２を有し、
【００６７】
【数１８】
0006483107.tif 000019
【００６８】
の条件を満たさなくてはならない。入力数列ｖ、並びに変換数列ｄ及びｂは、結合され、２つの出力数列ｇ及びｕが生成される。
【００６９】
数列
【００７０】
【数１９】
0006483107.tif 000020
【００７１】
は、
【００７２】
【数２０】
0006483107.tif 000021
【００７３】
と共に、基線の再トレンド化信号、すなわち、基線のトレンドｂを差し引いた後の初期数列ｖを意味する。
【００７４】
図４Ｂは、基線のトレンドが除去された生理学的信号ｕの計算を例示する。実線は、Δ＝０．００４ｓで抽出した典型的なＥＥＧ信号ｖの２０１点での一部分である。点線は、３１点の窓にＭＰＭＡ２を適用し、５回反復を用いた結果得られた基線ｂである。ＥＥＧ信号ｖから基線ｂを差し引いた後の、トレンドを除去した数列ｕは破線で示されている。
【００７５】
【数２１】
0006483107.tif 000022
【００７６】
を満たす数列
【００７７】
【数２２】
0006483107.tif 000023
【００７８】
は、ガイド関数と呼ばれている。ガイド関数は、信号の時間経過から、不適切な高周波数及び低周波数の要素を取り除いたものを用いて、得られる。これらの純化処理は、重要な機能要素の探索のガイドとしての数列ｇの役割を裏付ける。
【００７９】
図３を再度見ると、工程１０２中のトレンド除去に続いて、工程１０３が続き、ガイド関数ｇに基づいた分割が実施される。工程１０３中に実施される分割の目的は、機能要素に関連した信号のセグメントを定義する分割点を推定することである。分割点は、図５Ａに図示されたゼロ交差（ＺＣ）、並びに図５Ｂに図示されたガイド関数の絶対値の大域的な最小値及び局所的な最小値として定義される。図５Ｃは、縦の実線で示されたゼロ交差点を有するガイド関数ｇ及び縦の破線で示された最小値の例を示す。
【００８０】
より具体的には、点ｔｍは、以下の条件では、ゼロ交差点として限定される：
【００８１】
【数２３】
0006483107.tif 000024
【００８２】
点ｔｍは、以下の条件では、
【００８３】
【数２４】
0006483107.tif 000025
【００８４】
の大域的な最小値又は局所的な最小値である：
【００８５】
【数２５】
0006483107.tif 000026
【００８６】
分割点の数列
【００８７】
【数２６】
0006483107.tif 000027
【００８８】
はこのように形成され、セグメントの長さ
【００８９】
【数２７】
0006483107.tif 000028
【００９０】
は、以下を満たす分割点
【００９１】
【数２８】
0006483107.tif 000029
【００９２】
によって定義される：
【００９３】
【数２９】
0006483107.tif 000030
【００９４】
一般的に、ガイド関数ｇは、トレンドを除去した信号ｕの近似として考えられ得る。近似の正確性は、ＭＰＭＡ１の異なるパラメータを用いて修正することができる。従って、分割手順の時間分解能は、ＭＰＭＡ１のパラメータの適切な選択によって規定される。
【００９５】
図６Ａ及び図７Ａは、ＭＰＭＡ１の異なるパラメータを用いて、トレンドを除去したあるＥＥＧ信号又はトレンドを除去した同一のＥＥＧ信号ｕ（灰色線）に適合した２つのガイド関数（破線）を示す。より具体的には、図６Ｂ及び図７Ｂ中に別々に示されたガイド関数ｇ１及びガイド関数ｇ２は、各々三点の窓、及び七点の窓を用いて計算されている。両ケースとも、２反復を用いて計算されている。図６Ａに記載の結果は、トレンドを除去した数列ｕ及びガイド関数ｇ１の間のずれが少ないことを示す。よって、ガイド関数ｇ１は、トレンドを除去した信号を詳細に再現し、結果として、トレンドを除去した数列ｕの実質的に厳密な再現として見なすことができる。縦線で示す分割点は１１のゼロ交差点及び８の最小値を含む。
【００９６】
逆に、図７Ａに記載の結果、及び図７Ｂとの図６Ｂの比較は、より広い窓を用いた同じＭＰＭＡ手順によって、分割点から最小値を取り除き、結果として、複雑なピークを持つ波形を巨大な一つの波形によって置換することを示す。いくつかの場合、処理の重要な点は、機能的に不適切な活性を取り除き、ＥＥＧ、ＥＭＧ、ＥＣＧ等の、大域的スケールの生理学的信号の時間領域において、有力な動的トレンドを強調する点である。しかし、不規則的な高周波数の信号の変動は、物理学的に重要な発生源から生成され得る。例えば、ＥＥＧの場合、３０Ｈｚを超える周波数を示す活性は、ガンマリズムとして知られている。よって、装置の特徴に依存して、高分解能の分割及び低分解能の分割の両方が、これらの異なる組み合わせと同様に、複素非定常信号の要素構成の検討のための強力な計算道具を提供する。
【００９７】
工程１０３中で適用される信号変換は、時間領域において実施されるが、周波数領域の概念とも関連し得る。図６Ａ及び図７Ａから明らかなように、ＭＰＭＡ１中で使用される狭い窓は、より高い周波数の帯域幅においてローパスフィルタリングを実施し、一方で、ＭＰＭＡ２中で使用される広い窓は、比較的低い周波数の帯域幅においてローパスフィルタリングを実施する。しかし、当初の信号は、非定常プロセスであるため、固定された「カットオフ」周波数とともに従来の周波数帯に用いることは不明瞭である。
【００９８】
再度、図３を参照すると、工程１０３中の分割に続いて、方法１００は、セグメント再帰として、工程１０４から１０８を実行し、当該再帰の連続したサイクルは、定義された信号セグメントの連続を扱う。各々のセグメントに対して、要素モデルが特定され、当該要素モデルのパラメータは、対応する単一信号断片に適合している。このセグメント中、信号に適合した要素モデル各々は、以下のセグメントへと拡張されるため、工程１０４から１０８の反復工程の目的は、時間的要素の重複を解消することでもある。要素の時間的重複の解消は、一般的な方法では裏付けられていない信号処理の多くの分野において緊急課題である。この課題への如何なるアプローチも、要素の発達において、予想できる構造及び予想パターンを記述するため、及び過去の軌道から要素の今後の値を予測するため、要素モデル作成に関与しなくてはならない点に、主要な難しさがある。
【００９９】
「物理学に基づいた」アプローチ（Ｃａｎｄｙ  Ｊ．Ｖ．  “Ｍｏｄｅｌ−ｂａｓｅｄ  ｓｉｇｎａｌ  ｐｒｏｃｅｓｓｉｎｇ”．  Ｊｏｈｎ  Ｗｉｌｅｙ  ＆  Ｓｏｎｓ，  Ｎｅｗ  Ｊｅｒｓｅｙ，  ２００６  参照）の欠点は、発生源の明確な物理学モデルという点において、通常、信号に係る十分な情報を得ることができないという事実に生じている。物理学に基づいたアプローチに対して、本明細書に開示の方法は、決定論的カオスという形で、非定常信号の主要要素を定義するカオスに基づいたモデリング（ＣＢＭ）を用いた一意的な方法論を紹介している。このカオスに基づいたモデリングアプローチは、正規化効果の実験論的結果（Ｍｅｌｋｏｎｉａｎ  Ｄ，  Ｂｌｕｍｅｎｔｈａｌ  Ｔ，  Ｇｏｒｄｏｎ  Ｅ．  “Ｎｕｍｅｒｉｃａｌ  Ｆｏｕｒｉｅｒ  ｔｒａｎｓｆｏｒｍ  ｓｐｅｃｔｒｏｓｃｏｐｙ  ｏｆ  ＥＭＧ  ｈａｌｆ−ｗａｖｅｓ：  ｆｒａｇｍｅｎｔａｒｙ−ｄｅｃｏｍｐｏｓｉｔｉｏｎ−ｂａｓｅｄ  ａｐｐｒｏａｃｈ  ｔｏ  ｎｏｎｓｔａｔｉｏｎａｒｙ  ｓｉｇｎａｌ  ａｎａｌｙｓｉｓ”．  Ｂｉｏｌ  Ｃｙｂｅｒｎ．，  １９９９，  ｖｏｌ．  ８１，  ｐｐ．  ４５７−４６７  参照）によって、奨励されている。この研究結果は、決定論的かつ不規則な要素の干渉が原因となり、微視的スケールからの発生源の活性が、初期発生源の物理学的性質よりむしろ、大域的スケールにおいて確率論的不変性を反映した統計学的手順へと、変換されることを示す。よって、大域的スケールの生理学的要素は、複数の初期発生源の累積統計的集合として、発達する。
【０１００】
この発見的基準の普遍的かつ大域的スケールの要素モデルを数式化するため、本方法は、ランダム変数の集合の漸近挙動を定義する手段として、中心極限定理を引用する（Ｇｎｅｄｅｎｋｏ  Ｂ，  Ｋｏｌｍｏｇｏｒｏｖ  Ａ．  “Ｌｉｍｉｔ  ｄｉｓｔｒｉｂｕｔｉｏｎｓ  ｆｏｒ  ｓｕｍｓ  ｏｆ  ｉｎｄｅｐｅｎｄｅｎｔ  ｒａｎｄｏｍ  ｖａｒｉａｂｌｅｓ”．  Ａｄｄｉｓｏｎ−Ｗｅｓｌｅｙ  Ｐｕｂｌｉｓｈｉｎｇ  Ｃｏｍｐａｎｙ，  Ｒｅａｄｉｎｇ，  Ｍａｓｓ．，  １９５４  参照）。平均
【０１０１】
【数３０】
0006483107.tif 000031
【０１０２】
及び分散
【０１０３】
【数３１】
0006483107.tif 000032
【０１０４】
とともに、一定数のランダム変数
【０１０５】
【数３２】
0006483107.tif 000033
【０１０６】
の和と仮定すると、指数フーリエ変換
【０１０７】
【数３３】
0006483107.tif 000034
【０１０８】
の漸近挙動は、
【０１０９】
【数３４】
0006483107.tif 000035
【０１１０】
のようなガウス関数に対して極めて一般的な条件下において、収束し（ｓｅｅ  Ｐａｐｏｕｌｉｓ  Ａ．  “Ｔｈｅ  Ｆｏｕｒｉｅｒ  ｉｎｔｅｇｒａｌ  ａｎｄ  ｉｔｓ  ａｐｐｌｉｃａｔｉｏｎｓ”．  ＭｃＧｒａｗ−Ｈｉｌｌ  Ｂｏｏｋ  Ｃｏ，  ＮｅｗＹｏｒｋ，  １９６２  参照）
【０１１１】
【数３５】
0006483107.tif 000036
【０１１２】
式中、Ｆ（ｉω）は、指数フーリエ変換（複素変数ｉωの関数、
【０１１３】
【数３６】
0006483107.tif 000037
【０１１４】
）、
【０１１５】
【数３７】
0006483107.tif 000038
【０１１６】
及び
【０１１７】
【数３８】
0006483107.tif 000039
【０１１８】
である。
【０１１９】
よって、単一発生源の影響は、大域的スケールの変数として意味を持つ概略パラメータσ及びηへと変換される。
【０１２０】
物理学用語では、
【０１２１】
【数３９】
0006483107.tif 000040
【０１２２】
は、多数の微視的スケールの発生源から生成される大域的スケールの時間プロセスの振幅スペクトルとして、見なされ得る。
【０１２３】
この理論的予測の実験的立証は、聴覚の事象関連電位の研究によって、提供されてきた（Ｍｅｌｋｏｎｉａｎ  Ｄ，  Ｇｏｒｄｏｎ  Ｅ，  Ｒｅｎｎｉｅ  Ｃ，  Ｂａｈｒａｍａｌｉ  Ｈ．  Ｄｙｎａｍｉｃ  ｓｐｅｃｔｒａｌ  ａｎａｌｙｓｉｓ  ｏｆ  ｅｖｅｎｔ−ｒｅｌａｔｅｄ  ｐｏｔｅｎｔｉａｌｓ．  Ｅｌｅｃｔｒｏｅｎｃｅｐｈａｌｏｇｒａｐｈｙ  ａｎｄ  Ｃｌｉｎｉｃａｌ  Ｎｅｕｒｏｐｈｙｓｉｏｌｏｇｙ，  １９９８，  ｖｏｌ．１０８，  ｐｐ．  ２５１−２５９  参照）。動的スペクトル解析は、ガウス関数が、聴覚の事象関連電位の主要な後期潜時要素の振幅スペクトルのプロファイルに関する普遍的な記述を提供することを示した。
【０１２４】
方程式（７）の周波数領域の複素スペクトルの時間領域対応値は、無限時間スケールの間隔に関連した該当する文献において定義される。方程式（８）と同様に、時間領域の解決策は以下の形のガウス関数である傾向がある：
【０１２５】
【数４０】
0006483107.tif 000041
【０１２６】
逆に、開示の方法の信号処理のルーティーンは、生理学的に実現可能なシステム中において、一過性と見なされ得る因果時間関数に対応する。
【０１２７】
無限時間スケールから半無限時間スケールへと時間領域解析を移動させるため、開示の方法は、ｔ≧０において、２つのガウス関数の和として定義される、原因関数ｗ（ｔ）を導入し、
【０１２８】
【数４１】
0006483107.tif 000042
【０１２９】
式中、
【０１３０】
【数４２】
0006483107.tif 000043
【０１３１】
【数４３】
0006483107.tif 000044
【０１３２】
【数４４】
0006483107.tif 000045
【０１３３】
であり、Ｕ（ｔ）は単位ステップ関数であり、εは無次元定数である。
【０１３４】
方程式（１１）及び（１２）と方程式（７）を比較すると、方程式（７）からのパラメータη（平均）は、定数εを用いて、パラメータσ（分散）の比として、方程式（１１）及び（１２）中に記載されることがわかる。定数εの値は、非定常信号のタイプに依存する。従って、開示の方法を用いた、ヒト及び動物のＥＥＧに関する解析及びコンピュータシミュレーションは、これらの応用例に対する普遍定数として、ε＝２であることを示唆する。この推定値は、定数εのデフォルト値だと考えられる。
【０１３５】
パラメータσ＝１及び定数ε＝２を満たす関数ｇＡ（ｔ）及び−ｇＢ（ｔ）を、図８Ａ中に各々、破線及び点線で示す。これら数式の両端は、
【０１３６】
【数４５】
0006483107.tif 000046
【０１３７】
である。
【０１３８】
図８Ａ中、実線で示す、関数ｇ（ｔ）は、無限範囲の奇関数である。関数ｇ（ｔ）のフーリエ変換は以下であり：
【０１３９】
【数４６】
0006483107.tif 000047
【０１４０】
式中、
【０１４１】
【数４７】
0006483107.tif 000048
【０１４２】
【数４８】
0006483107.tif 000049
【０１４３】
である。
【０１４４】
ｇ（ｔ）は奇関数であるため、
【０１４５】
【数４９】
0006483107.tif 000050
【０１４６】
及び
【０１４７】
【数５０】
0006483107.tif 000051
【０１４８】
である。
【０１４９】
方程式（１７）からの最後の三角積分は、以下の表記が用いられるｗ（ｔ）のフーリエ正弦変換である。
【０１５０】
【数５１】
0006483107.tif 000052
【０１５１】
方程式（１０）によってｗ（ｔ）を置換した後、この積分は以下の形式のように、解析的に記載される（Ｅｒｄｅｌｙｉ，  Ａ．  （Ｅｄ）．  “Ｔａｂｌｅｓ  ｏｆ  ｉｎｔｅｇｒａｌ  ｔｒａｎｓｆｏｒｍｓ“，  Ｖｏｌ．  １．  ＭｃＧｒａｗ−Ｈｉｌｌ，  Ｎｅｗ  Ｙｏｒｋ，  １９５４，  ｐ．  ７８  参照）：
【０１５２】
【数５２】
0006483107.tif 000053
【０１５３】
この式は、複素スペクトル
【０１５４】
【数５３】
0006483107.tif 000054
【０１５５】
の虚数部分として見なすことができ、式中、
【０１５６】
【数５４】
0006483107.tif 000055
【０１５７】
は、振幅スペクトルであり、かつ、
【０１５８】
【数５５】
0006483107.tif 000056
【０１５９】
は、位相スペクトルである。
【０１６０】
従って、方程式（２１）及び（２２）の周波数関数の組み合わせは、ｔ≧０において、ｗ（ｔ）の時間領域信号を補足的に定義する周波数領域に対応するものと見なすことができる。この関数を図８Ｂに示す。ｗ（ｔ）が最大に達する近似値は、
【０１６１】
【数５６】
0006483107.tif 000057
【０１６２】
である。最大値は、
【０１６３】
【数５７】
0006483107.tif 000058
【０１６４】
である。
【０１６５】
物理的には、ｗ（ｔ）は、複数の微視的スケールの発生源の集合の活性化によって、ｔ＝０において誘導される一過的なプロセスと見なすことができる。別々の発生源の振る舞いは、確率論的要素及び決定論的要素の両方に規定され、決定論的カオスとして特徴づけることができる。このような振る舞いの形態の主要な特徴は、決定論的カオスを生成している系の非線形性である。この観点において、方程式（１１）で記載されるプロセスは、以下の非線形微分方程式システムで規定される動的システムのインパルス応答と見なすことができる：
【０１６６】
【数５８】
0006483107.tif 000059
【０１６７】
式中、
【０１６８】
【数５９】
0006483107.tif 000060
【０１６９】
は、
【０１７０】
【数６０】
0006483107.tif 000061
【０１７１】
の導関数を意味する。
【０１７２】
合成非定常プロセスのモデルは、以下の形式で定義される：
【０１７３】
【数６１】
0006483107.tif 000062
【０１７４】
式中、κｉ及びτｉは、加重及び開始時間を意味する自由なパラメータである。
【０１７５】
準ガウシアンカーネル（ＱＧＫ）と呼ばれる、要素モデルは以下のように定義される：
【０１７６】
【数６２】
0006483107.tif 000063
【０１７７】
上記式は、３つのパラメータのセット、すなわち、
【０１７８】
【数６３】
0006483107.tif 000064
【０１７９】
によって特徴づけられている。
【０１８０】
一般的な観点では、ｖＮ（ｔ）の推定は、パラメータ
【０１８１】
【数６４】
0006483107.tif 000065
【０１８２】
がランダム変数であるという事実に従っている。非定常性の特定の側面は、新たな発生源が、基礎を成すシステムに結合した際、時刻τｉにおいてｖ（ｔ）の変化を生み出すシステムの構造的要素であるということである。決定論的カオスの重要な側面は、方程式（１０）〜（１２）のガウス依存性(Gaussian dependences)は、決定論的プロセスであるが、これら時間発展は、基礎を成す初期ランダム事象の漸近挙動として発展する、という事実によって、明らかにされている。微視的スケールにおける該当する変換の確率論的側面は、マルコフ過程、特に不均一な出生死滅過程という形で計算されてきた（Ｍｅｌｋｏｎｉａｎ  Ｄ，  Ｂｌｕｍｅｎｔｈａｌ  Ｔ．  Ｄｙｎａｍｉｃｓ  ｏｆ  ｔｈｅ  ｅｙｅｂｌｉｎｋ  ＥＭＧ  ａｔ  ｇｌｏｂａｌ  ａｎｄ  ｍｉｃｒｏｓｃｏｐｉｃ  ｓｃａｌｅｓ．  Ｔｈｅ  Ｏｐｅｎ  Ｃｙｂｅｒｎｅｔｉｃｓ  ａｎｄ  Ｓｙｓｔｅｍｉｃｓ  Ｊｏｕｒｎａｌ，２０１３，  ｖｏｌ．  ７，  ｐｐ．  １１−２２  参照）。
【０１８３】
方法１００は、工程１０４から１０８を実行し、方程式（２４）の各々ＱＧＫ
【０１８４】
【数６５】
0006483107.tif 000066
【０１８５】
に対するパラメータ
【０１８６】
【数６６】
0006483107.tif 000067
【０１８７】
を特定する。特定アルゴリズムは再帰として構成され、当該再帰のガイド原則は方程式（２４）中の「ｎ」と「ｎ−１」との間の以下のような関係から導かれる：
【０１８８】
【数６７】
0006483107.tif 000068
【０１８９】
関数ｖｎ（ｔ）は、概念上はモデルであるため、この関数ｖｎ（ｔ）は、トレンドを除去した数列ｕ（ｔ）の近似、すなわち、ｔ＞０において
【０１９０】
【数６８】
0006483107.tif 000069
【０１９１】
であると見なすことができる。ＱＧＫ
【０１９２】
【数６９】
0006483107.tif 000070
【０１９３】
の近似への寄与は、ｔ＝τｎにおける
【０１９４】
【数７０】
0006483107.tif 000071
【０１９５】
の初期値から始まる。つまり、以下のようである。
【０１９６】
【数７１】
0006483107.tif 000072
【０１９７】
τ１＝０とすると、
【０１９８】
【数７２】
0006483107.tif 000073
【０１９９】
ｎ＞１のとき、
【０２００】
【数７３】
0006483107.tif 000074
【０２０１】
である。
【０２０２】
方程式（２８）は、前項の影響がない場合において、再帰の最初のサイクルに対応している。方程式（２９）は、関数
【０２０３】
【数７４】
0006483107.tif 000075
【０２０４】
によって前項の影響を受ける再帰ｎ回目に対応している。この関数は、再帰の前のサイクルにおいて得られる、主要なパラメータ
【０２０５】
【数７５】
0006483107.tif 000076
【０２０６】
（ｉ＝１，・・，ｎ−１）によって解析的に定義される。
【０２０７】
【数７６】
0006483107.tif 000077
【０２０８】
の役割は、時間要素の修正である。この観点において、方程式（２９）の右項は、解析的ＱＧＫ  γｎ（ｔ）の重複を修正した準実験的対応値を表す。
【０２０９】
ＱＧＫ  γｎ（ｔ）のパラメータは、複素非線形関係によって、時間領域において相互に関連がある。周波数領域の利点は、周波数領域関係が、連続オーダーで実施される別々の手順によって、これらの関係を分離し、パラメータκｎ、σｎ及びτｎの連続を推定する方法を提供する点である。対応する周波数領域の操作の重要な局面は、κｎ及びσｎの推定値が、解析下で、時間関数の時間シフトに対して不変である振幅スペクトルのプロファイルによって、定義される、という点である。これは、パラメータτｎの暫定的な推定値（κｎ及びσｎの推定段階中）として、分割点
【０２１０】
【数７７】
0006483107.tif 000078
【０２１１】
の値を用いる方法を提供する。従って、方程式（２８）及び（２９）は、以下の形式中、時間から周波数への変換の目的のために、提示される：
【０２１２】
【数７８】
0006483107.tif 000079
【０２１３】
【数７９】
0006483107.tif 000080
【０２１４】
標準的な手順によって、時間から周波数への領域変換を裏付けるため、方程式（３１）の波形要素を
【０２１５】
【数８０】
0006483107.tif 000081
【０２１６】
だけ移動させ、単一の方程式として記載する：
【０２１７】
【数８１】
0006483107.tif 000082
【０２１８】
間隔は、
【０２１９】
【数８２】
0006483107.tif 000083
【０２２０】
であり、
【０２２１】
【数８３】
0006483107.tif 000084
【０２２２】
を満たす。この方程式は、経験的成分カーネル（ＥＣＫ）と名付けられる。
【０２２３】
工程１０５中、ＥＣＫ  ｙｎ（ｔ）は、有限余弦フーリエ変換及び有限正弦フーリエ変換の数値推定を用いて、時間領域から周波数領域へと変換される。より具体的には、Ｍｅｌｋｏｎｉａｎ  （２０１０）中（Ｍｅｌｋｏｎｉａｎ  Ｄ．  “Ｓｉｍｉｌａｒ  ｂａｓｉｓ  ｆｕｎｃｔｉｏｎ  ａｌｇｏｒｉｔｈｍ  ｆｏｒ  ｎｕｍｅｒｉｃａｌ  ｅｓｔｉｍａｔｉｏｎ  ｏｆ  Ｆｏｕｒｉｅｒ  ｉｎｔｅｇｒａｌｓ”．  Ｎｕｍｅｒ  Ａｌｇｏｒ．，  ２０１０，  ｖｏｌ．  ５４，  ｐｐ．  ７３−１００  参照）に開示の、相似基底関数(similar basis function)（ＳＢＦ）アルゴリズムを用いており、その内容は参照によって組み込まれる。
【０２２４】
ＳＢＦアルゴリズムは、多項式挿入を用いた三角積分の数値推定における、最大精度を提供する、ファイロン様式の原型である。高速フーリエ変換（ＦＦＴ）及びその改良によって実施される、従来のデジタルスペクトル解析に反して、ＳＢＦアルゴリズムは、時間領域及び周波数領域の両方において、計算点の厳密なグリッドを扱う必要がない。同じことは、工程１０５中、時間から周波数への変換に従う信号セグメントの長さが可変である本ケースにおいても、特に有利である。周波数領域において、計算点の柔軟な選択の主な利点は、対数周波数スケールを用いた周波数領域依存性を表すことができる点である。
【０２２５】
ＳＢＦアルゴリズムは、普遍的な手順で実現されているため、添え字「ｎ」（再帰における回数）は、アルゴリズムの以下の記載において省略する。
【０２２６】
ＳＢＦアルゴリズムは、以下の有限余弦及び正弦フーリエ変換の推定に対応している：
【０２２７】
【数８４】
0006483107.tif 000085
【０２２８】
【数８５】
0006483107.tif 000086
【０２２９】
式中、ω＝２πｆは角速度、及びｆは周波数である。これら方程式中に使用される表記
【０２３０】
【数８６】
0006483107.tif 000087
【０２３１】
は、ｔ≧λにおいて、ｙ（ｔ）＝０と仮定する。これらは、ｙ（ｔ）の複素スペクトルの実数部分及び複素部分であり、
【０２３２】
【数８７】
0006483107.tif 000088
【０２３３】
を満たし、式中、
【０２３４】
【数８８】
0006483107.tif 000089
【０２３５】
である。
【０２３６】
対応する振幅スペクトルは以下である：
【０２３７】
【数８９】
0006483107.tif 000090
【０２３８】
デジタル形式へのｙ（ｔ）の変換は、
【０２３９】
【数９０】
0006483107.tif 000091
【０２４０】
を満たす有限数の節点
【０２４１】
【数９１】
0006483107.tif 000092
【０２４２】
中の、標本値
【０２４３】
【数９２】
0006483107.tif 000093
【０２４４】
によってｙ（ｔ）を特定する。節点は均等に分けられる必要はない。よって、計算は、データ点のセット
【０２４５】
【数９３】
0006483107.tif 000094
【０２４６】
に対応している。これらデータ点を用いて、ｙ（ｔ）は、区分線形多項式ｈ（ｔ）によって標本に挿入される。挿入は、以下の挿入条件を満たさなければいけない：
【０２４７】
【数９４】
0006483107.tif 000095
【０２４８】
式中、
【０２４９】
【数９５】
0006483107.tif 000096
【０２５０】
である。
【０２５１】
【数９６】
0006483107.tif 000097
【０２５２】
の条件下で、ＳＢＦアルゴリズムは、相似基底関数：
【０２５３】
【数９７】
0006483107.tif 000098
【０２５４】
（式中、ａｉは加重係数であり、
【０２５５】
【数９８】
0006483107.tif 000099
【０２５６】
は、相似基底関数（ＳＢＦ）である）の和として、間隔［０，λ］間に形成される。ＳＢＦは、相似関係によって定義される：
【０２５７】
【数９９】
0006483107.tif 000100
【０２５８】
この単純にパラメータ化された時間スケーリングは、「三角形基底関数」（ＴＢＦ）と呼ばれる基本的な有限要素、
【０２５９】
【数１００】
0006483107.tif 000101
【０２６０】
より、ＳＢＦのファミリーを生成する。ＴＢＦは、図９Ｃ中に図示されている単位直角三角形である。
【０２６１】
有限要素を用いて、挿入式ｈ（ｔ）は、以下の形式で表される：
【０２６２】
【数１０１】
0006483107.tif 000102
【０２６３】
式中、
【０２６４】
【数１０２】
0006483107.tif 000103
【０２６５】
は、ハット関数である。ハット関数
【０２６６】
【数１０３】
0006483107.tif 000104
【０２６７】
は、以下のように定義される：
【０２６８】
【数１０４】
0006483107.tif 000105
【０２６９】
上記式は、メッシュ
【０２７０】
【数１０５】
0006483107.tif 000106
【０２７１】
における式である。図９Ａは、方程式（４１）中に記載の様式で、ハット関数
【０２７２】
【数１０６】
0006483107.tif 000107
【０２７３】
の連続（標本１、２、３、４）に続く、「０」標本へと取り付けられたＴＢＦの和として、挿入式ｈ（ｔ）の分解を例示する。
【０２７４】
方程式（３７）及び（４１）を比較すると、方程式（４１）のハット関数は、ＴＢＦによって方程式（３７）中で置換されていることがわかる。この結合において、有用だとわかる幾何学法則は、ハット関数（三角形ａｂｃ）を３つのＴＢＦ（三角形ｏｈｃ、三角形ｆｏｄ及び三角形ｇｏａ）の和に分解することを示す、図９Ｂによって、図示されている。方程式（３７）中のすべてのハット関数にこの分解を適用すると、挿入係数の推定のため、以下の方程式となる：
【０２７５】
【数１０７】
0006483107.tif 000108
【０２７６】
式中、
【０２７７】
【数１０８】
0006483107.tif 000109
【０２７８】
である。
【０２７９】
数値解という観点において、方程式（４２）の構築の主な利点は、ＴＢＦからの余弦フーリエ積分及び正弦フーリエ積分が、以下のような解析的形式で表される点である：
【０２８０】
【数１０９】
0006483107.tif 000110
【０２８１】
【数１１０】
0006483107.tif 000111
【０２８２】
図９Ｄは、Ｒｅ（方程式（４５））及びＩｍ（方程式（４６））と図示した、これら方程式を示す。
【０２８３】
フーリエ変換の理論（相似定理）によると、時間領域における横座標の圧縮は、横座標の拡張に加えて周波数領域における縦座標の縮小に対応する。これらの処理は、ＳＢＦ及びＴＢＦのフーリエ積分間の以下の相似関係を明らかにする：
【０２８４】
【数１１１】
0006483107.tif 000112
【０２８５】
【数１１２】
0006483107.tif 000113
【０２８６】
これらの方程式を用いて、方程式（３８）の補間式(interpolant)の余弦積分及び正弦積分は、以下の解析方程式で表される：
【０２８７】
【数１１３】
0006483107.tif 000114
【０２８８】
【数１１４】
0006483107.tif 000115
【０２８９】
ｈ（ｔ）はｙ（ｔ）への補間式であるため、これらのフーリエ変換は、方程式（３３）及び（３４）の有限フーリエ積分の推定値として見なされる。結果として、
【０２９０】
【数１１５】
0006483107.tif 000116
【０２９１】
この観点において、方程式（４９）及び（５０）中の関数Ｒｃ及びＲｓの、その解析式である（４５）及び（４６）による置換は、以下の計算方程式を導く：
【０２９２】
【数１１６】
0006483107.tif 000117
【０２９３】
【数１１７】
0006483107.tif 000118
【０２９４】
方程式（４３）からの荷重係数推定後、方程式（５２）及び（５３）は、データ点へと直接適用される。
【０２９５】
ＳＢＦアルゴリズムの利点は、角周波数の任意のセットにおいて、周波数特性を計算できる点である。以下のパラメータ特定手順は、対数周波数スケールを用いる。従って、周波数特性の標本
【０２９６】
【数１１８】
0006483107.tif 000119
【０２９７】
は、複数点
【０２９８】
【数１１９】
0006483107.tif 000120
【０２９９】
において推定され、
式中、１を超えるＣは、対数横座標スケールにおいてサンプリングレートを定義する定数であり、ωＡは、初期周波数である。
【０３００】
１０に対して、１０、２０、５０及び１００標本に対応するＣの典型値は、各々１．２５８９、１．１２２、１．０４７１、及び１．０２３３である。
【０３０１】
振幅スペクトルは、同じ点において以下のように計算される：
【０３０２】
【数１２０】
0006483107.tif 000121
【０３０３】
ｎ番目の再帰サイクルと仮定すると、対応する振幅スペクトル
【０３０４】
【数１２１】
0006483107.tif 000122
【０３０５】
は、κｎ（加重）パラメータ及びσｎ（分散）パラメータの特定の役割を果たす。この目的のため、方法１００（図３）は、次に工程１０６に進む。
【０３０６】
角周波数ωＡが、
【０３０７】
【数１２２】
0006483107.tif 000123
【０３０８】
のように選択される場合、加重は、
【０３０９】
【数１２３】
0006483107.tif 000124
【０３１０】
のように定義される。
【０３１１】
この観点で決定的に重要なものは、実験的振幅スペクトルのモデルとしての方程式（２１）のガウシアン振幅スペクトルの妥当性試験である。当該試験は、方程式（５５）から計算された正規化振幅スペクトルに適用され、方程式（５６）からパラメータκｎの推定値を用いる：
【０３１２】
【数１２４】
0006483107.tif 000125
【０３１３】
方程式（２１）のモデルと、このスペクトルとの比較は、テンプレートとして、ガウシアンスペクトルを用いて、対数周波数スケールにおいて、実施される。後者は、以下のように解析的に定義され、
【０３１４】
【数１２５】
0006483107.tif 000126
【０３１５】
式中、νは、実数の無次元変数である。この関数を図１０Ａに示す。数値計算のため、テンプレートは、複数点の標本値によって定義され、
【０３１６】
【数１２６】
0006483107.tif 000127
【０３１７】
式中、１を超えるＣは、方程式（５４）中と同じ定数であり、Ｍ＋１は、テンプレートの標本数である。
【０３１８】
よって、テンプレートは、
【０３１９】
【数１２７】
0006483107.tif 000128
【０３２０】
である、セグメント
【０３２１】
【数１２８】
0006483107.tif 000129
【０３２２】
に拡張する。境界点
【０３２３】
【数１２９】
0006483107.tif 000130
【０３２４】
は、矢印Ａ及びＢによって図１０Ａに図示されている。
【０３２５】
テンプレートと比較される実験的な正規化振幅スペクトルの標本は、対数角周波数スケール上の以下の点に対応している：
【０３２６】
【数１３０】
0006483107.tif 000131
【０３２７】
周波数範囲（
【０３２８】
【数１３１】
0006483107.tif 000132
【０３２９】
）及びサンプリングレートは、Ｉパラメータ及びＣパラメータによって定義される。これらパラメータの適切な選択は、デジタル形式において、振幅スペクトルの十分正確な表示方法を提供する。振幅スペクトルの周波数の範囲は、テンプレートが定義した周波数の範囲を超えなくてはいけない。従って、Ｉ＞Ｍである。
【０３３０】
テンプレート適合は、横座標スケールを通して、ガウシアン振幅テンプレートを動かすことによって、実行される。これは、テンプレート標本が、変化しない間、横座標上の位置が、定数Ｃによる、周波数の乗法によって、連動して移動する段階的処理である。よって、第一の移動は、点
【０３３１】
【数１３２】
0006483107.tif 000133
【０３３２】
の数列への方程式（５７）の複数点の数列の変化によって実行される。ｋ番目の移動の後の周波数点の数列は、
【０３３３】
【数１３３】
0006483107.tif 000134
【０３３４】
である。図１０Ｂは、実験的な正規化振幅スペクトル
【０３３５】
【数１３４】
0006483107.tif 000135
【０３３６】
（黒い実線）、並びに点Ｔ１、Ｔ２及びＴ３の異なる移動におけるガウシアン振幅テンプレートのセットを例示する。
【０３３７】
テンプレートスペクトル及び実験的な振幅スペクトルの間のずれは、平均二乗誤差（ＭＳＥ）を用いて、各々のテンプレートの位置において評価され：
【０３３８】
【数１３５】
0006483107.tif 000136
【０３３９】
式中、ｋは位置番号である（ｋ＝１、・・・、Ｉ−Ｍ＋１）。
【０３４０】
角周波数の点では、テンプレートの初期位置は、実験的スペクトルの周波数範囲
【０３４１】
【数１３６】
0006483107.tif 000137
【０３４２】
を含んでおり、式中、
【０３４３】
【数１３７】
0006483107.tif 000138
【０３４４】
である。ｋ回の移動後、テンプレートに含まれる周波数範囲は、
【０３４５】
【数１３８】
0006483107.tif 000139
【０３４６】
であり、式中
【０３４７】
【数１３９】
0006483107.tif 000140
【０３４８】
である。
【０３４９】
最良のテンプレート適合は、ＭＳＥの最小値によって示される。対応する移動回数をＫとすると、パラメータ
【０３５０】
【数１４０】
0006483107.tif 000141
【０３５１】
は、対数角周波数スケールにおいて、ＭＳＥが最小であるテンプレートの位置を定義する。この適合が、方程式（５５）の実験的な振幅スペクトルの適切なモデルとして、方程式（２１）の振幅スペクトルを実証するかどうかを決定するため、これら関数間の残差を、標本によって定義される無次元誤差関数
【０３５２】
【数１４１】
0006483107.tif 000142
【０３５３】
として見なし、
【０３５４】
【数１４２】
0006483107.tif 000143
【０３５５】
式中、ｊは、０からＭの値を取り、
【０３５６】
【数１４３】
0006483107.tif 000144
【０３５７】
である。
【０３５８】
ここでは、
【０３５９】
【数１４４】
0006483107.tif 000145
【０３６０】
の両方は、等しい初期値
【０３６１】
【数１４５】
0006483107.tif 000146
【０３６２】
を持つ無次元関数である。従って、
【０３６３】
【数１４６】
0006483107.tif 000147
【０３６４】
と仮定する。周波数の増加に伴う誤差関数の変化は、周波数の範囲
【０３６５】
【数１４７】
0006483107.tif 000148
【０３６６】
を推定するため、閾値εＴを用いて評価され、上記範囲においては、
【０３６７】
【数１４８】
0006483107.tif 000149
【０３６８】
である。
【０３６９】
境界角周波数ωＪにおける
【０３７０】
【数１４９】
0006483107.tif 000150
【０３７１】
の最終値は、
【０３７２】
【数１５０】
0006483107.tif 000151
【０３７３】
として示される。この推定値は、以下の図式に従って１から５の整数値を取るカテゴリー化パラメータξを用いて、テンプレート適合の質を確認するために用いられる：
【０３７４】
【数１５１】
0006483107.tif 000152
【０３７５】
階級、すなわち、ξの値によって、モデリングの結果は以下のように評価される：Ｅ−優れている(excellent)、Ｇ−よい(good)、Ｓ−満足できる(satisfactory)、Ｗ−不十分(weak)、Ｕ−不成功(unsuccessful)。
【０３７６】
指令に依存して、階級Ｗ又は階級Ｕの場合、加重は無視される、すなわち、κｎ＝０である。
【０３７７】
図１１Ａ及び図１１Ｂは、閾値
【０３７８】
【数１５２】
0006483107.tif 000153
【０３７９】
を用いてテンプレート適合の２つの実施例の誤差関数を図示する。対応する境界周波数は、図１１Ａ及び図１１Ｂ中、縦の点線によって図示されている。これらの線は、図１１Ａ及び図１１Ｂ中の振幅スペクトルと、各々、０．０６５及び０．３０１の値において交差する。その結果、対応するモデルは、各々、階級Ｅ及びＧに属する。
【０３８０】
最良な適合は、テンプレート標本が、物理エンティティ、つまり、対数横座標スケール上の
【０３８１】
【数１５３】
0006483107.tif 000154
【０３８２】
の対応サンプルにおける角周波数の値、と関連し得るテンプレート標本における、テンプレートのいくつかの固定位置を定義する。この観点において、テンプレート標本は、以下の形式の角周波数の関数として、表され、
【０３８３】
【数１５４】
0006483107.tif 000155
【０３８４】
式中、下付き文字「Ｆ」は、テンプレートの最終位置を示し、σｎは、ｎ番目のサイクルの再帰における方程式（２１）より得られたσパラメータを示す。
【０３８５】
任意のωｍにおいて、所要のパラメータσｎは、以下のように計算される
【０３８６】
【数１５５】
0006483107.tif 000156
【０３８７】
計算の安定した正確性のため、角周波数ωｍは、
【０３８８】
【数１５６】
0006483107.tif 000157
【０３８９】
が０．４から０．８の値を取る周波数の範囲から選択される。
【０３９０】
図１２は、テンプレート適合の典型的な結果を例示する。中央のパネルは、縦の点線で示された分割点を有するＥＥＧの一部を示す。分割点間で抽出された６個の連続したＥＥＧ断片の振幅スペクトルは、上側のパネル及び下側のパネルにおいて、実線で表される。対応するテンプレートは、灰色の線で記載される。テンプレート適合の例は各々、以下の情報によって、補足される：モデルの階級（Ｅ又はＢ）、角周波数
【０３９１】
【数１５７】
0006483107.tif 000158
【０３９２】
及び値
【０３９３】
【数１５８】
0006483107.tif 000159
【０３９４】
である。
【０３９５】
図１３Ａ中に同じ信号を再描写するため、分割点の間の信号断片の振幅スペクトルが、図１３Ｂ中に示される。図１３Ｃは、振幅と周波数との両方を正規化した後のこれらスペクトル特性の重ね合わせを示す。これらの点は、ガウシアン振幅テンプレート（黒い線）との適合が、低い周波数の範囲又は中間の周波数の範囲において優れているという点において、典型的である。コンピュータシミュレーションは、高い周波数の偏差は、機能要素と混同した不適切な活性によって生成されることを示してきた。
【０３９６】
方法１００（図３）に再び目を向けると、工程１０６に続いて、方法１００は、開始時間τｎが推定される工程１０７へと続く。この手順は、方程式（５５）を用いた振幅スペクトルの正規化と同様の方法において正規化した後、正規化振幅スペクトル
【０３９７】
【数１５９】
0006483107.tif 000160
【０３９８】
及び複素スペクトルの対応する虚部を用いている：
【０３９９】
【数１６０】
0006483107.tif 000161
【０４００】
式中、下付き文字ｎは、再帰サイクル数である。
【０４０１】
図１４Ａ及び図１４Ｂは、ＡＳ（正規化振幅スペクトル
【０４０２】
【数１６１】
0006483107.tif 000162
【０４０３】
）及びＩｍ（複素スペクトル
【０４０４】
【数１６２】
0006483107.tif 000163
【０４０５】
の正規化虚部）と示すこれら関数の典型的な形を例示する。これら関数の一般的関係は、
【０４０６】
【数１６３】
0006483107.tif 000164
【０４０７】
の形式中、位相スペクトル
【０４０８】
【数１６４】
0006483107.tif 000165
【０４０９】
を含む。
【０４１０】
方程式（２２）の位相スペクトルと同様、
【０４１１】
【数１６５】
0006483107.tif 000166
【０４１２】
は、以下のように、角周波数への線形従属として表される。
【０４１３】
【数１６６】
0006483107.tif 000167
【０４１４】
方程式（２２）及び方程式（６９）の位相スペクトル間の差は、シフト係数
【０４１５】
【数１６７】
0006483107.tif 000168
【０４１６】
によって、説明される。
【０４１７】
物理的に、この移動は、要素開始時間τｎ及び分割点
【０４１８】
【数１６８】
0006483107.tif 000169
【０４１９】
の値の間の差を示す。この観点で、
【０４２０】
【数１６９】
0006483107.tif 000170
【０４２１】
である。
【０４２２】
先行の再帰サイクルにおける、ＳＢＦアルゴリズムの手順は、解析下において、ＥＣＫの未知の開始時間の置換として、分割点を取る。従って、分割点は、数値積分が開始する時刻としての機能を果たす。κｎ及びσｎパラメータの元である振幅スペクトルは、対応する時間関数の時間移動に対して不変であるため、この条件は、当初定義したパラメータ値において、ほとんど効果を有さない。
【０４２３】
τｎ推定は、
【０４２４】
【数１７０】
0006483107.tif 000171
【０４２５】
に従った方程式（７１）を参照する。
【０４２６】
この方程式を用いるため、未知のパラメータβｎを推定すべきである。方程式（６８）及び方程式（６９）を参照すると、パラメータβｎは、関数
【０４２７】
【数１７１】
0006483107.tif 000172
【０４２８】
の間の以下の関係によって、完全に定義される：
【０４２９】
【数１７２】
0006483107.tif 000173
【０４３０】
この方程式からパラメータβｎを見出すため、角周波数の単一点に対して
【０４３１】
【数１７３】
0006483107.tif 000174
【０４３２】
の推定値を有することは、十分である。しかし、計算誤差を減らすため、２つの推定値は、異なる角周波数
【０４３３】
【数１７４】
0006483107.tif 000175
【０４３４】
において、計算される。
【０４３５】
【数１７５】
0006483107.tif 000176
【０４３６】
と仮定すると、これら周波数の選択は、ω１ｎからω２ｎの間隔において、
【０４３７】
【数１７６】
0006483107.tif 000177
【０４３８】
は、ωの増加と共に減少すると推定される。
【０４３９】
図１４Ａ及び図１４Ｂ中、横座標スケール上、文字ａで示される縦の点線は、
【０４４０】
【数１７７】
0006483107.tif 000178
【０４４１】
を満たす点として、角周波数ω１ｎを定義する。方程式（７３）への参照は、この点において
【０４４２】
【数１７８】
0006483107.tif 000179
【０４４３】
であることを示す。方程式（７２）へのこの値の挿入は、開始時間の第一の推定値を提供する：
【０４４４】
【数１７９】
0006483107.tif 000180
【０４４５】
方程式（２１）及び方程式（２２）の形式中、周波数領域モデルに従うと、
【０４４６】
【数１８０】
0006483107.tif 000181
【０４４７】
の期待される形式は、図１４Ａ中に図示されるように、極小値に続くゼロ交差によって特徴づけられる。このような場合を考えると、角周波数ω２ｎは、
【０４４８】
【数１８１】
0006483107.tif 000182
【０４４９】
が横座標スケールと交差する、すなわち、
【０４５０】
【数１８２】
0006483107.tif 000183
【０４５１】
を満たす点として定義される。
【０４５２】
図１４Ａ中、ω２ｎは、横座標スケール上のｂから始まる縦の点線によって示される。位相スペクトル
【０４５３】
【数１８３】
0006483107.tif 000184
【０４５４】
の値は、
【０４５５】
【数１８４】
0006483107.tif 000185
【０４５６】
を定義する。
【０４５７】
方程式（７２）は、以下のように、位相移動の第二の推定値を提供する。
【０４５８】
【数１８５】
0006483107.tif 000186
【０４５９】
いくつかの場合、ＥＣＫの不適切な要素は、
【０４６０】
【数１８６】
0006483107.tif 000187
【０４６１】
の期待される周波数領域のプロファイルを歪める。図１４Ｂ中のプロットＩｍは、横座標スケールを交差する代わりに、最小に到達する
【０４６２】
【数１８７】
0006483107.tif 000188
【０４６３】
の特徴的な例である。対応する角周波数は、横座標スケール上のｂから始まる縦の点線によって示される。開始時間の第二の推定値を得るため、この周波数は、方程式（７５）中、ω２ｎを置換する。
【０４６４】
τｎ［１］及びτｎ［２］を仮定すると、開始時間の最終的な選択は、σ＝σｎを満たす方程式（１９）から
【０４６５】
【数１８８】
0006483107.tif 000189
【０４６６】
の間のずれの、異なるτｎの値に対する評価が元になっている。当該ずれは、以下の平均二乗誤差によって計算され、
【０４６７】
【数１８９】
0006483107.tif 000190
【０４６８】
式中、標本
【０４６９】
【数１９０】
0006483107.tif 000191
【０４７０】
は、τｎの推定値であるτｎ［１］又はτｎ［２］を用いて計算され、Ｎは、
【０４７１】
【数１９１】
0006483107.tif 000192
【０４７２】
のように、選択される整数である。
【０４７３】
ＭＳＥ［１］及びＭＳＥ［２］が、τｎ［１］及びτｎ［２］と関連があると仮定すると、ＭＳＥ［１］≦ＭＳＥ［２］の場合、τｎ＝τｎ［１］である。それ以外は、τｎ＝τｎ［２］である。
【０４７４】
よって、工程１０６及び工程１０７は、検討中のセグメントに対するパラメータ
【０４７５】
【数１９２】
0006483107.tif 000193
【０４７６】
を推定する。工程１０８は、方法１００がその後、工程１０４から工程１０８が、次の分割のパラメータ
【０４７７】
【数１９３】
0006483107.tif 000194
【０４７８】
を推定する工程１０４へと戻る前に、方程式（２６）を用いて、以下の部分モデルを推定する。
【０４７９】
【数１９４】
0006483107.tif 000195
【０４８０】
図２Ｂは、対応するＥＣＫに適合するパラメータ
【０４８１】
【数１９５】
0006483107.tif 000196
【０４８２】
を有するＱＧＫと共に、図２Ａ（Ｎ１００、ｐ２００、Ｎ２００、Ｐ３ａ及びＰ３ｂ）中に示す各々機能要素の置換の典型的な結果を図示する。セグメントは、図２Ａ中の実験的な対応結果として、同じ記号Ｎ１００、Ｐ２００、Ｎ２００、Ｐ３ａ及びＰ３ｂによって、明示されている。図２Ｃは、比較として図２Ａ中に示す当初のＥＲＰ信号と共に、要素ＱＧＫの和、すなわち、大域的モデルを図示する。信号分解及び特定要素のパラメータは、発明の方法の一般的な手順から完全に導かれたものであり、以降で適用される現象に適合させるための如何なる修正も行っていないことから、これは、満足できる程度の適合である。
【０４８３】
この結論の指示における更なる証拠は、異なる発生源によって生成される非定常信号を用いて、得られている。この点を例示するため、図１５Ａは、Δ＝０．０５６ｓにおいて抽出された、話している人の呼吸の波形の６０秒セグメントを示す。プロットは、３点又は４７点の窓フィルタに対して５つのパスを用いて、呼吸センサーからトレンドを除去した信号を表す。ｙ軸ユニットは、生の呼吸センサーの離散的な出力である。
【０４８４】
図１５Ｂは、特定した５４個のＱＧＫの重ね合わせを描く。すべてのモデル要素の和として合成された呼吸波形の大域的モデルは、図１５Ｃ中、実線で示される。同じ図において、点線は、図１５Ａから再現された当初の呼吸波形である。
【０４８５】
ＥＥＧと同様に、呼吸の波形は、高度な非定常信号である。しかし、電気の発生源（皮質ニューロン）によって生成されるＥＥＧと異なり、呼吸器の動きは、胸壁の動きを規定する数多くの生体機械的要素に依存する。よって、呼吸器の動きの基礎となる物理学的プロセスは、脳波の電気力学とは、著しく異なっている。それにもかかわらず、本明細書に開示の方法は、両方の適用に対し、等しく効果的な要素の特定ツールを提供する。実験的波形への顕著に正確なＱＧＫの適合は、異なる生物種及び異なる実験条件下において記録される、様々な電気生理学的信号（ＥＥＧ、ＥＭＧ、ＥＣＧ）及び呼吸の波形にとって典型的である。これらの結果は、発明の方法論のカオスに基づくモデリング（ＣＢＭ）のアプローチを強く裏付ける。
【０４８６】
すべてのセグメントが処理され、パラメータξ、κ、σ、τが工程１０３において形成された各々のセグメントに対して推定された後、方法１００は、モデリングの結果が、指標としてパラメータξを用いて検討される、工程１０９へと続く。各々のセグメントに対するパラメータξの値は、定義された閾値ξＴと比較される。ξｎ≧ξＴの場合、対応する要素のモデルは、更なる処理段階から除かれる。ξＴの典型値は、４及び５である。
【０４８７】
ＮＩの不適当な要素を除いた後、残っている要素は、ＮＳ−１が分割点の数である、ｎ＝１からＭ＝ＮＳ−ＮＩの連続的かつ順番に、番号が付け替えられる。対応するパラメータのセットは、以下のように記述され、
【０４８８】
【数１９６】
0006483107.tif 000197
【０４８９】
式中、
【０４９０】
【数１９７】
0006483107.tif 000198
【０４９１】
は、対応するパラメータのＭ値を各々含むベクトルである。
【０４９２】
ベクトルτの要素が、定義されたモデル要素の発達が開始する時間点の有限集合を定義する一方で、ベクトルκ及びベクトルσの要素は、モデル要素の波形を定義する。この観点において、ベクトル
【０４９３】
【数１９８】
0006483107.tif 000199
【０４９４】
は、データ点のセット
【０４９５】
【数１９９】
0006483107.tif 000200
【０４９６】
を提供し、これは時系列として、パラメータκの値の数列を表す。さらに、ベクトル
【０４９７】
【数２００】
0006483107.tif 000201
【０４９８】
は、パラメータσの値の時系列を提供する。
【０４９９】
時系列は、信号処理の多くの方法が適用される数値データの標準形式である。よって、信号を標準パラメータのセットに分解する開示の方法１００は、データ正規化、データ変換、主成分分析、自己相関、自己回帰モデル、データクラスタリング、予測等、といった、非定常信号処理の、従来の時系列解析の方法へと適用する方法を提供する。
【０５００】
要素モデル（信号解析）のセットへの非定常信号の分解の結果を用いて、方法１００は、異なる問題を解決するための方法を提供し、これが信号合成である。これは、あらかじめ定義された階級に属する要素の複合体の形成に関与する。
【０５０１】
合成のための、要素の収集は、工程１１０によって実行されるパラメータ階級分け手順から開始する。開始時間は、変化しないままで、階級分けは、各々特定された要素モデルに対して、パラメータξ、κ及びσへと適用される。より具体的には、パラメータの部分集合は、以下の条件に従って、モデルパラメータＰＭの全体集合から定義され：
【０５０２】
【数２０１】
0006483107.tif 000202
【０５０３】
式中、
【０５０４】
【数２０２】
0006483107.tif 000203
【０５０５】
は、対応するパラメータの境界（最初と最後との）値を定義する定数である。
【０５０６】
方程式（７８）より、パラメータの部分集合は、以下のように示され、
【０５０７】
【数２０３】
0006483107.tif 000204
【０５０８】
式中、オーバーハット（＾）は、選択された要素を意味し、Ｌは、要素番号である。明らかに、Ｌ≦Ｍを満たす。
【０５０９】
最後に、方法１００は、工程１１０からのパラメータの部分集合
【０５１０】
【数２０４】
0006483107.tif 000205
【０５１１】
が、以下の形式において非定常信号の部分モデルの合成のため使用される、工程１１０において終了する。
【０５１２】
【数２０５】
0006483107.tif 000206
【０５１３】
階級分け及び合成は、複数の目的を有する手順であり、異なる目的を追求できる。例えば、極端に小さい振幅、又は極端に大きい振幅を持つ要素の除外は、信号の不適当な要素を取り除くことができる。κ及びσの特定値の選択によって、振動要素等を強調できる。
【０５１４】
モデルパラメータ
【０５１５】
【数２０６】
0006483107.tif 000207
【０５１６】
のフルセット及び／又は、選択されたパラメータのセット
【０５１７】
【数２０７】
0006483107.tif 000208
【０５１８】
は、更なる処理の前に、コンピュータシステムの記憶媒体に蓄積されるか、外部コンピュータシステムへと伝達され得る。更なる処理は、最適デコンボリューション、発見的クラスタリングアルゴリズム、再帰的なベイズ推定等、といった異なったモデルに基づいた信号処理方法を含み得る。
【０５１９】
前記の事項は、本発明のいくつかの実施形態のみを記述しており、発明の目的及び方針から外れない限りは、改良及び／又は変更を行ってもよく、実施形態は、例示されるにすぎず、これに制限されるものではない。
【図面の簡単な説明】
【０５２０】
【図１Ａ】図１Ａは、頭皮の表面に取り付けた二つのセンサーの間で測定される電位の典型的な波形を示す。
【図１Ｂ】図１Ｂは、ＥＥＧ信号が示す極めて非定常的な特徴を例示する。
【図２Ａ】図２Ａは、聴覚オッドボールパラダイムにおいて検出される典型的なＥＲＰ信号を示す。
【図２Ｂ】図２Ｂは、図２Ａに示した各セグメントを、そのセグメントに適合するパラメータを有する準ガウシアンカーネルを用いて置換した典型的な結果を示す。
【図２Ｃ】図２Ｃは、準ガウシアンカーネル成分の和として得られたモデルを、比較のための図２Ａに示した元のＥＲＰ信号とともに示す。
【図３】図３は、生理信号を機能要素に適応的に分解し、その後、各機能要素を特徴づけるパラメータを決定する、波形分析方法の概略フロー図を示す。
【図４Ａ】図４Ａは、図３に記載の方法中で応用した、動的トレンド及び基線のトレンドの推定の概略フロー図を示す。
【図４Ｂ】図４Ｂは、典型的な生理信号、マルチパスの移動窓による平均化を離散的な生理信号に適用して得られた基線のトレンド、及びトレンドを除去した基線の信号を示す。
【図５Ａ】図５Ａ及び図５Ｂは、ゼロ交差点から結果的に得られた分割点、並びに大域的及び局所的な最小値での分割点を各々、示す。
【図５Ｂ】図５Ａ及び図５Ｂは、ゼロ交差点から結果的に得られた分割点、並びに大域的及び局所的な最小値での分割点を各々、示す。
【図５Ｃ】図５Ｃは、ゼロ交差点（実線の縦線）と最小値（点線の縦線）とにおいて、それぞれ、ＥＥＧ信号を分割した例を示す。
【図６】図６Ａは、高い時間分解能を持つ分割点を提供する初期信号及びガイド関数を示す；図６Ｂ及び図６Ｃは、図６Ａに記載のガイド関数及び初期信号を、縦線で示した分割点とともに再度示す。
【図７】図７Ａは、図６Ａに記載の初期信号、及び低い時間分解能による信号分割を提供するガイド関数を示す；図７Ｂ及び図７Ｃは、図７Ａに記載のガイド関数及び初期関数を、縦線で示す分割点とともに再度示す。
【図８Ａ】図８Ａは、準ガウシアンカーネルの構成に用いたガウシアン成分を示す。
【図８Ｂ】図８Ｂは、準ガウシアンカーネルの一例を示す。
【図９】図９Ａは、有限成分の和への挿入箇所の分解を例示する；図９Ｂは、三つの三角形基底関数の和へのハット関数の分解を示す；図９Ｃは、同様の基底関数アルゴリズム中で使用した、三角形基底関数を示す；図９Ｄは、三角形基底関数の複素スペクトルの実部と虚部とを示す。
【図１０Ａ】図１０Ａは、パラメータσの推定に用いた振幅テンプレートを示す。
【図１０Ｂ】図１０Ｂは、パラメータσの推定過程における、正規化した実験振幅スペクトル、並びに異なる位置の振幅テンプレート、Ｔ１、Ｔ２及びＴ３を例示する。
【図１１Ａ】図１１Ａは、階級Ｅに属する正規化振幅スペクトルに、テンプレートを重ね合わせた結果を示す。
【図１１Ｂ】図１１Ｂは、階級Ｇに属する正規化振幅スペクトルに、テンプレートを重ね合わせた結果を示す。
【図１２】図１２は、図中央のパネルに記載の、ＥＥＧ信号の異なるセグメントの正規化振幅スペクトルに重ね合わせた振幅テンプレートのいくつかの結果を例示する。
【図１３Ａ】図１３Ａは、低い時間分解能を用いたＥＥＧ分割の例を示す。
【図１３Ｂ】図１３Ｂは、図１３Ａに記載の１０セグメントの振幅スペクトルを示す。
【図１３Ｃ】図１３Ｃは、振幅と周波数との両方を正規化した後、図１３Ｂの振幅スペクトルを重ね合わせたものを示す。
【図１４Ａ】図１４Ａは、パラメータτの推定値が得られることを基準にした場合の横座標周波数スケール上の点ａ及びｂを示す。
【図１４Ｂ】図１４Ｂは、横座標スケール上の点ｂが、複素スペクトルの虚部の最小値によって定義された場合を例示する。
【図１５Ａ】図１５Ａは、話している人の呼吸の波形例を示す。
【図１５Ｂ】図１５Ｂは、図１５Ａの信号において特定した５４個のＱＧＫの波形を示す。
【図１５Ｃ】図１５Ｃは、特定した５４個のＱＧＫの和として合成されたモデルによって、図１５Ａの信号を再現した場合の、正確性を示す。
【図１Ａ】
0006483107.tif 000210
【図１Ｂ】
0006483107.tif 000211
【図２Ａ】
0006483107.tif 000212
【図２Ｂ】
0006483107.tif 000213
【図２Ｃ】
0006483107.tif 000214
【図３】
0006483107.tif 000215
【図４Ａ】
0006483107.tif 000216
【図４Ｂ】
0006483107.tif 000217
【図５Ａ】
0006483107.tif 000218
【図５Ｂ】
0006483107.tif 000219
【図５Ｃ】
0006483107.tif 000220
【図６】
0006483107.tif 000221
【図７】
0006483107.tif 000222
【図８Ａ】
0006483107.tif 000223
【図８Ｂ】
0006483107.tif 000224
【図９】
0006483107.tif 000225
【図１０（Ａ）】
0006483107.tif 000226
【図１０（Ｂ）】
0006483107.tif 000227
【図１１Ａ】
0006483107.tif 000228
【図１１Ｂ】
0006483107.tif 000229
【図１２】
0006483107.tif 000230
【図１３Ａ】
0006483107.tif 000231
【図１３Ｂ】
0006483107.tif 000232
【図１３Ｃ】
0006483107.tif 000233
【図１４（Ａ）】
0006483107.tif 000234
【図１４（Ｂ）】
0006483107.tif 000235
【図１５Ａ】
0006483107.tif 000236
【図１５Ｂ】
0006483107.tif 000237
【図１５Ｃ】
0006483107.tif 000238