(57)【特許請求の範囲】
【請求項１】
誘導方向の異なる第１の心電図信号および第２の心電図信号を取得する取得手段と、
前記第１の心電図信号から、前記第２の心電図信号の推定信号を生成する推定手段と、
前記第２の心電図信号と、前記推定手段が生成した推定信号との差信号を生成する減算手段と、
前記差信号の低周波成分を呼吸努力に関する信号として抽出する抽出手段と、
を有することを特徴とする心電図解析装置。
【請求項２】
前記第１の心電図信号と前記第２の心電図信号の一方がＸ軸方向の誘導であり、他方がＹ軸もしくはＺ軸方向の誘導であることを特徴とする請求項１に記載の心電図解析装置。
【請求項３】
前記第１の心電図信号と前記第２の心電図信号の一方がＣＣ５誘導であり、他方がＮＡＳＡ誘導であることを特徴とする請求項１または請求項２に記載の心電図解析装置。
【請求項４】
前記推定手段は、前記第１の心電図信号を入力信号とし、前記推定信号を出力信号とするシステムを多項式モデルによって同定することにより、前記推定信号を生成することを特徴とする請求項１から請求項３のいずれか１項に記載の心電図解析装置。
【請求項５】
前記多項式モデルがＡＲＸ(Auto Regressive with eXogenous)多項式モデル、ＯＥ(Output-Error)多項式モデル、ＢＪ(Box-Jenkins)多項式モデルのいずれかであることを特徴とする請求項４に記載の心電図解析装置。
【請求項６】
被検者に装着された電極を用いて、誘導方向の異なる第１の心電図信号および第２の心電図信号を測定する測定手段と、
測定された前記第１の心電図信号および前記第２の心電図信号を記憶する記憶手段と、
前記記憶手段に記憶された前記第１の心電図信号および前記第２の心電図信号を解析する請求項１から請求項５のいずれか１項に記載の心電図解析装置と、
を有することを特徴とする心電計。
【請求項７】
誘導方向の異なる第１の心電図信号および第２の心電図信号を取得する取得工程と、
前記第１の心電図信号から、前記第２の心電図信号の推定信号を生成する推定工程と、
前記第２の心電図信号と、前記推定工程で生成された推定信号との差信号を生成する減算工程と、
前記差信号の低周波成分を呼吸努力に関する信号として抽出する抽出工程と、
を有することを特徴とする心電図解析方法。
【請求項８】
コンピュータを、請求項１から請求項５のいずれか１項に記載の心電図解析装置が有する各手段として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は心電図解析装置、心電図解析方法、および心電計に関し、特には心電図信号から呼吸情報を抽出する技術に関する。
【背景技術】
【０００２】
睡眠時無呼吸（低呼吸）症候群(Sleep Apnea (Hypopnea) Syndrome: SASまたはSAHS)は、睡眠の分断による過度の日中傾眠を伴い、睡眠時に１０秒以上連続した無呼吸が５回／時間又は３０回／７時間以上発生する疾患である。以下では便宜上、ＳＡＳとＳＡＨＳを併せて単にＳＡＳと呼ぶ。
【０００３】
ＳＡＳの診断には、口鼻呼吸、いびき、動脈血酸素飽和度（ＳｐＯ２）、呼吸努力運動、脳波、眼球運動、頤筋筋電図、心電図、体位などの測定を終夜行う終夜睡眠ポリグラフ（Poly-Somno-Graphy:ＰＳＧ）検査が必要である。しかしながら、終夜睡眠ポリグラフ検査は、通常、病院での宿泊を要し、またカニューレや多数の電極、センサ類を装着する必要があるため、被検者の負担が比較的大きい。そのため、終夜睡眠ポリグラフ検査を行う前に、より簡便なスクリーニング検査（簡易ＰＳＧ検査）を行うことが多い。例えば特許文献１には、無呼吸及び低呼吸の有無、回数やＳｐＯ２値、いびき、呼吸努力運動といったスクリーニング検査に用いる睡眠時無呼吸検査装置が開示されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００６−３２０７３２号公報
【特許文献２】特開２０１０−５１３８７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
簡易ＰＳＧ検査によって被検者の負担は軽減されるが、鼻や口付近に呼吸センサを装着して就寝するのは煩わしく、より負担の少ない方法が望まれている。
【０００６】
一方、無呼吸や低呼吸によって心拍数が変化することに着目して、心電図から無呼吸や低呼吸を検出する方法も検討されている（特許文献２）。特許文献２では、無呼吸や低呼吸の周期と同期して、心拍数が徐脈と頻脈を繰り返す周期性心拍変動パターン（ＣＶＨＲ：cyclic variation of heart rate）を検出している。心拍数は例えば心電図（ＥＣＧ）のＲ波の間隔（Ｒ−Ｒ間隔）により得ることができるため、測定の負担が少ないホルタ心電図などを利用してＣＶＨＲを検出することができる。
【０００７】
しかしながら、心拍数の変化は必ずしも無呼吸や低呼吸によって生じるとは限らないため、心拍数が安定しない他の要因を有する被検者に対しては、ＣＶＨＲを精度良く測定することは困難である。例えば、心房細動の患者はもともと心拍数が不安定であるため、仮にＳＡＳに罹患していても、心拍数の変動が心房細動によるものか、ＳＡＳによるものかを心電図から判別することは困難である。
【０００８】
本発明はこのような従来技術の課題に鑑みてなされたものであり、被検者の呼吸情報を心電図信号から精度良く取得することが可能な心電図解析装置、心電図解析方法、および心電計の提供を目的とする。
【課題を解決するための手段】
【０００９】
上述の目的は、誘導方向の異なる第１の心電図信号および第２の心電図信号を取得する取得手段と、第１の心電図信号から、第２の心電図信号の推定信号を生成する推定手段と、第２の心電図信号と、推定手段が生成した推定信号との差信号を生成する減算手段と、差信号の低周波成分を呼吸努力に関する信号として抽出する抽出手段と、を有することを特徴とする心電図解析装置によって達成される。
【発明の効果】
【００１０】
このような構成により、本発明によれば、被検者の呼吸情報を心電図信号から精度良く取得することが可能な心電図解析装置、心電図解析方法、および心電計を提供できる。
【図面の簡単な説明】
【００１１】
【図１】本発明の実施形態に係る心電図解析装置の一例としてのホルタ心電計の機能構成例を示すブロック図である。
【図２】本発明の実施形態に係る解析処理を模式的に示したブロック図である。
【図３】図２における推定モデルの例を示す図である。
【図４】本発明の実施形態に係る心電図解析処理に関するフローチャートである。
【図５】本発明の実施例１を説明するための波形図である。
【図６】本発明の実施例１を説明するための波形図である。
【図７】本発明の実施例２を説明するための波形図である。
【図８】本発明の実施例２を説明するための波形図である。
【発明を実施するための形態】
【００１２】
以下、図面を参照して本発明をその例示的な実施形態に基づいて詳細に説明する。図１は、本発明の実施形態に係る心電図解析装置の一例としての、ホルタ心電計１００の機能構成例を示すブロック図である。なお、本実施形態の心電図解析装置は、心電図信号（誘導信号）を被検者から取得するための構成を有しているが、本発明において心電図信号を被検者から取得するための構成は必須でない。本発明に係る心電図解析装置は、予め測定された心電図データを、直接またはネットワークを介して接続された記憶装置や記憶媒体から取得するなど、任意の方法で取得可能であればよい。
【００１３】
図１において、ＣＰＵ１はＲＯＭ２に格納されている制御プログラムをＲＡＭ３に読み出して実行することにより、後述する心電図解析処理を始めとした、ホルタ心電計１００の各種機能を実現する。ＲＯＭ２はＣＰＵ１が実行するプログラムや、メニュー画面などを表示するためのＧＵＩデータ、ユーザ設定データ、初期設定データなど、処理に必要なパラメータ等を記憶する不揮発性メモリであり、少なくとも一部が書き換え可能であってよい。ＲＡＭ３はＣＰＵ１が実行するプログラムを展開する領域や、変数やデータ等の一時記憶領域として用いられる。メモリカード４は生体電極から入力される心電図信号や、別の生体電気信号やデータをデジタルデータの形式で記憶する記憶装置である。メモリカード４は、カードスロット５に対して着脱可能に装着される。
【００１４】
表示部６は例えば液晶表示装置である。操作部８は電源のオン、オフや測定の開始、停止、各種イベント入力などを行ったり、各種の設定を行なったりするためのスイッチ、ボタンなどからなる。操作部８は表示部６に設けられたタッチパネルを含んでもよい。
【００１５】
また、アナログ−デジタル変換器（Ａ／Ｄ変換器）９は心電図電極１２から入力されるアナログ信号をデジタル信号に変換する。センサ用Ｉ／Ｆ１０は、心電図電極１２からの心電図信号（誘導信号）を取得するためのインターフェースである。なお、センサ用Ｉ／Ｆ１０は心電図電極１２以外にも、ＳｐＯ２（動脈血酸素飽和度）センサや血圧・脈波測定用のカフなどを接続可能であってもよい。
通信インターフェース２０は例えばホストコンピュータやプリンタ等の外部機器４０と通信を行うための通信インターフェースであり、有線および／または無線通信規格に準拠した構成を有する。
【００１６】
このような構成を有するホルタ心電計１００を用いて心電図信号（誘導信号）の取得及び記録を行う場合、被検者の体表面の所定の位置に心電図電極１２の遠端に設けられた複数の電極パッドを取り付け、心電図電極１２の近端に設けられたコネクタをセンサ用Ｉ／Ｆ１０のコネクタに接続する。
【００１７】
本発明においては、心電図信号から呼吸信号を検出するため、呼吸努力に伴う胸部や腹部の筋肉の活動電位変化の影響の受けやすさの異なる、複数の誘導方向の心電図信号を用いる。典型的には、誘導方向が互いに直交する２方向を含むような誘導法を用いて測定する。
【００１８】
具体的には、左右（Ｘ軸）、上下（Ｙ軸）、および前後（Ｚ軸）方向の誘導方向のうち、少なくともＸ軸方向の誘導（Ｘ誘導）を含む、２つ以上の異なる誘導方向で測定することが好ましい。これは、呼吸努力による胸部や腹部の運動に起因する信号が、Ｘ軸方向の誘導に多く含まれると考えられるからである。
【００１９】
例えば、Ｘ誘導（ＣＭ５，ｍＣＭ５，ＣＣ５）、Ｙ誘導（ＮＡＳＡ，双極ａＶＦ，ＣＭｆ）、Ｚ誘導（ＮＡＳＡ，ｍＶ１，ＣＭ２，ＣＢ２，双極Ｖ３）の中から、Ｘ誘導と、他の誘導を選択することができる。また、標準１２誘導や、ベクトル心電図誘導法（例えばFrank誘導法）を用いることもできる。２チャンネルの心電図信号を測定する場合、Ｘ軸方向の誘導とＹ軸方向の誘導の組み合わせとして、例えばＣＣ５とＮＡＳＡ誘導の組み合わせを好適に用いることができる。
【００２０】
例えば操作部８の操作によって電源が投入されると、ＣＰＵ１は初期化処理など、心電図信号の取得を開始する前のタイミングで、心電図信号を取得するための動作モードの設定処理を行ことができる。また、ＣＰＵ１は必要に応じて、心電図信号の記録動作の開始前に、日時や、被検者を特定可能な情報（例えば患者ＩＤなど）を、ユーザに設定させるための入力画面などを表示してもよい。
【００２１】
そして、例えば操作部８からの記録開始指示の入力や、電源投入時からの所定時間経過などに応じて、ＣＰＵ１は心電図の記録動作を開始する。なお、心電図信号の取得自体は、記録開始指示より前に実行されてよい。例えば、表示部６に現在の動作モードに応じた波形表示を行うことで、動作モードの設定が正しく行われているかどうかを波形表示からユーザが確認することを可能にする。
【００２２】
心電図信号はＡ／Ｄ変換器９により所定のサンプリングレート（例えば１２５Ｈｚ〜数ＫＨｚ）でサンプリングされ、デジタル信号の形式でメモリカード４や外部機器４０に記録される。なお、心電図信号の測定ならびに記録動作は本発明と直接関係がなく、また従来と同様であってよいため、これ以上の説明は省略する。
【００２３】
（心電図信号の解析処理）
次に、誘導方向の異なる２つの心電図信号を解析して、呼吸努力に伴う胸部や腹部の動きを示す信号（以下、呼吸信号という）を呼吸情報として検出する方法に関して説明する。図２は、本実施形態に係る心電図解析処理を模式的に示したブロック図である。
【００２４】
解析部３００は、２チャンネルの心電図信号ＥＣ１，ＥＣ２から呼吸信号を抽出して出力する。解析部３００が有する各機能ブロック（後述）はソフトウェア、ハードウェア、あるいは両者の組み合わせによって実現することができる。ソフトウェアによって実現する部分は例えばＣＰＵ１がプログラムを実行することによって実施可能である。
【００２５】
解析部３００において、必要に応じて設けられるローパスフィルタ３０５および３０６は、心電図信号に含まれる不連続点（例えばＲ波）に存在する高周波成分を抑制することにより、後述する推定信号の低周波成分の推定精度を向上させる。ローパスフィルタ３０５および３０６のカットオフ周波数は例えば数Ｈｚ程度に設定することができる。
【００２６】
推定信号生成部３１０は、推定モデルを用いて、一方の心電図信号（ここではｃｈ１の心電図信号ＥＣ１）から、他方の心電図信号（ここではｃｈ２の心電図信号ＥＣ２）の推定信号（ここではｃｈ２の心電図信号ＥＣ２の推定信号＾ＥＣ２）を生成する。なお、ローパスフィルタ３０５および３０６が設けられる場合には、推定信号生成部３１０に入力するｃｈ１の心電図信号ＥＣ１はフィルタ処理後の信号であるため、ｃｈ２の心電図信号ＥＣ２の推定信号＾ＥＣ２もローパスフィルタ３０６が適用された後の心電図信号ＥＣ２に相当する。減算器３２０は、実測された心電図信号ＥＣ２と、推定信号＾ＥＣ２との差信号を生成する。呼吸信号抽出フィルタ３３０は、差信号からＲ波成分やバイアス誤差を除去するとともに、（呼吸信号成分に影響を与えない程度の）平滑化処理を行い、呼吸信号を抽出する。呼吸信号抽出フィルタ３３０は、例えばＲ波成分を除去するためのメジアンフィルタと、呼吸成分を通過させるバンドパスフィルタとの組み合わせによって実現できるが、他の組み合わせを用いてもよい。
【００２７】
呼吸信号と実測された心電図信号ＥＣ２との差信号が減算器３４０で生成され、推定誤差信号として推定信号生成部３１０にフィードバックされる。推定信号生成部３１０は、推定誤差信号が小さくなるように、推定モデルのパラメータを例えば最小二乗法（Least Mean Square：ＬＭＳ）、正規化最小二乗法（正規化ＬＭＳ：Normalized LMS）、予測誤差法などを用いて推定し、更新する。
【００２８】
ここで、推定信号生成部３１０について説明する。推定信号生成部３１０は、一方の心電図信号を入力とし、他方の心電図信号の推定信号を出力するシステムに推定モデルを適用し、推定モデルのパラメータを推定することによってシステムを同定し、推定信号＾ＥＣ２を生成する。本実施形態で推定信号生成部３１０は、推定モデルとして多項式モデルを用いる。
【００２９】
心電図電極１２を通じて測定される心電図信号ＥＣ１，ＥＣ２は、心臓の活動電位が体表面に伝達した信号である。そのため、心臓近傍における心電図信号ＥＣ０が、胸部モデルを通じて観測された信号であると考えることができる。ここでは、第１チャンネル（ｃｈ１）の心電図信号（第１の心電図信号）ＥＣ１が胸部モデル１ ２１０を通じて観測され、第２チャンネル（ｃｈ２）の心電図信号（第２の心電図信号）ＥＣ２が胸部モデル２ ２２０を通じて観測された信号であるとする。誘導方向が異なるため、チャンネルごとに異なる胸部モデルを想定している。
【００３０】
また、一方の心電図信号（ここではｃｈ２の心電図信号ＥＣ２）に、呼吸努力による胸部や腹部の運動に起因する信号（呼吸信号）が加算器２３０によって重畳されるものとする。なお、呼吸信号は両方の心電図信号に重畳しうるが、誘導方向が異なることにより、心電図信号に重畳する呼吸信号の強度には差が生じる。そのため、便宜上、呼吸信号の重畳が多い方の心電図信号についてのみ呼吸信号（呼吸による外乱）が重畳するものとしている。
【００３１】
このような前提に基づくと、測定される心電図信号ＥＣ１，ＥＣ２は以下の様に表すことができる。
ＥＣ１＝Ｓ１ＥＣ０            （１）
ＥＣ２＝Ｓ２ＥＣ０＋ｒｓ      （２）
ここで、ｒｓは呼吸信号である。さらに、式（２）からＥＣ０を除去すると、
ＥＣ２＝（Ｓ２／Ｓ１）ＥＣ１＋ｒｓ      （３）
となる。
従って、係数Ｓ２／Ｓ１を表す推定モデルのパラメータを推定することで、呼吸信号の含まれないＥＣ２の推定信号＾ＥＣ２を得ることができる。
【００３２】
この推定信号＾ＥＣ２と、実測されたＥＣ２との差信号は、ＥＣ２に含まれる心電図信号以外の成分、すなわち呼吸信号ｒｓの成分を多く含んだ信号である。そのため、差信号から、Ｒ波の成分やバイアス誤差などを除去することにより得られる信号を呼吸信号＾ｒｓとして抽出することができる。
【００３３】
なお、式（１）および（２）において、ＥＣ2からＥＣ１を推定した場合、呼吸        ＥＣ1＝（Ｓ１／Ｓ２）（ＥＣ2−ｒｓ）      （３）'
となる。この場合、係数Ｓ１／Ｓ２を表す推定モデルのパラメータを推定することにより、推定信号＾ＥＣ1と、実測されたＥＣ1との差信号は（Ｓ１／Ｓ２）ｒｓとなる。この場合、呼吸信号ｒｓを直接得ることはできないが、呼吸信号ｒｓを含んだ差信号が得られるため、呼吸信号を抽出できる。従って、計測された２つの心電図信号のどちらを推定信号生成部３１０に入力しても、呼吸波形を抽出することができる。
【００３４】
図３に、推定信号生成部３１０が推定モデルとして用いることのできる多項式モデルの例として、ＡＲＸ(Auto Regressive with eXogenous)多項式モデル、ＯＥ(Output-Error)多項式モデル、ＢＪ(Box-Jenkins)多項式モデルを示す。各モデルにおいて、ｕ（ｋ）が入力信号（ＥＣ１）、ｙ（ｋ）が出力信号（＾ＥＣ２）、ω（ｋ）がノイズ成分を示す。
【００３５】
なお、胸部モデルＳ１，Ｓ２がそれぞれＬＣＲ回路で表されるとすると、胸部モデルの伝達関数は双２次フィルタ(biquad filter)で表すことができる。従って、各モデルの次数は２を基本とし、心電図信号が含む高周波成分などに応じて調整することができる。
【００３６】
図４のフローチャートを参照して、本実施形態のホルタ心電計１００による心電図信号の解析処理動作について説明する。解析処理は心電図信号の記録時に実施してもよいし、記録済みの心電図信号に対して実施してもよい。また、上述したように、ホルタ心電計１００のような測定装置とは別の、測定機能を持たない装置で実施してもよい。以下では、予めＲＯＭ２に記憶された心電図信号の解析プログラムをＲＡＭ３に展開してＣＰＵ１が実行することによって解析処理を実行するものとするが、解析処理の少なくとも一部がＡＳＩＣ等の専用ハードウェアによって実現されてもよい。
【００３７】
Ｓ１０１でＣＰＵ１は、メモリカード４もしくは通信可能に接続された外部機器４０から、解析対象の心電図信号を取得する。取得する心電図信号は、上述の通り、少なくとも２つの、互いに異なる方向における心起電力の変化を表す心電図信号である。ここでは典型例として、Ｙ誘導を第１の心電図信号ＥＣ１として、またＸ誘導を第２の心電図信号ＥＣ２として取得するものとする。
【００３８】
なお、Ｓ１０１で取得する複数の心電図信号は同時に測定された信号であり、測定期間の一部の心電図信号を取得する場合も、同じ測定期間についての心電図信号を取得する。ＣＰＵ１は読み出した心電図信号をＲＡＭ３に記憶する。ＳＡＳのスクリーニングを行う場合のように解析区間が長時間にわたる場合、Ｓ１０１では所定時間分ずつ心電図信号を取得する。
【００３９】
次にＳ１０３でＣＰＵ１は、ＲＡＭ３に記憶されている第１の心電図信号ＥＣ１のサンプル列に対し、必要に応じてローパスフィルタ処理を適用してから推定信号生成部３１０に入力して推定モデルの同定を開始するとともに、第２の心電図信号ＥＣ２のサンプル列を減算器３２０に入力する。ここでは、推定モデルのパラメータ推定を、新たなサンプルが入力される都度行う（オンライン同定）するものとする。推定信号生成部３１０によって第２の心電図信号ＥＣ２の推定信号＾ＥＣ２のサンプルが順次生成され、減算器３２０に出力される。
【００４０】
Ｓ１０５で減算器３２０は実測された第２の心電図信号ＥＣ２から、推定信号生成部３１０が出力する推定信号＾ＥＣ２を減じて差信号を生成し、呼吸信号抽出フィルタ３３０に出力する。なお、ＣＰＵ１は、第１の心電図信号ＥＣ１のサンプル列が推定信号生成部３１０に入力される前にローパスフィルタ処理を適用されている場合には、第２の心電図信号ＥＣ２のサンプル列についても減算器３２０に入力する前に第１の心電図信号ＥＣ１に適用したのと同様のローパスフィルタ処理を適用する。
【００４１】
Ｓ１０７で呼吸信号抽出フィルタ３３０は、差信号から呼吸信号＾ｒｓを抽出して出力する。
【００４２】
以下、解析対象区間について、第１および第２の心電図信号ＥＣ１，ＥＣ２のサンプルについて順次Ｓ１０１〜Ｓ１０７の処理を繰り返し実行する。学習期間を設けずに推定モデルのパラメータを推定する場合、解析処理の開始直後は推定信号＾ＥＣ２と実測された第２の心電図信号ＥＣ２との整合性は高くないが、時間の経過と共に推定モデルのパラメータが最適化され、整合性が向上していく。
【００４３】
推定信号生成部３１０で推定モデルをオフライン同定する場合、予め定めた数の連続するサンプル群を用いて推定信号を生成する処理を、重複しない所定の一定期間（データ窓）に含まれるサンプル群ごとに順次実行する。例えば、１つ目のデータ窓に含まれるサンプル群を用いて対応区間の推定信号を生成したら、次のデータ窓に含まれるサンプル群を用いて次の区間の推定信号を生成する処理を繰り返し行う。データ窓の長さは例えば６０〜９０秒程度とすることができる。
【００４４】
ＣＰＵ１は、呼吸信号抽出フィルタ３３０が抽出した呼吸信号＾ｒｓをＲＡＭ３に保存した後、メモリカード４に記録したり、外部装置４０に送信したりする。あるいは、ＣＰＵ１は、呼吸信号＾ｒｓを用いて無呼吸や低呼吸の検出を並行して実行してもよい。なお、無呼吸や低呼吸の検出には、同じ解析対象区間に測定された他の信号（ＳｐＯ２）を必要に応じて用いてもよい。
【００４５】
（実施例）
予め記録された、ＣＣ５誘導とＮＡＳＡ誘導（サンプリングレート１２５Ｈｚ）の心電図信号を用い、上述の方法で呼吸信号＾ｒｓを得るまでの実例を、図５〜図８を用いて説明する。図５および図６は健常者（心拍数が変動する疾患のない被検者）についての実施例１を、図７および図８は心房細動の患者についての実施例２を示している。
【００４６】
図５〜図８において、ｅｃｇ１およびｅｃｇ２は実測された心電図信号を示し、ｅｃｇ１をＮＡＳＡ誘導、ｅｃｇ２をＣＣ５誘導とした。
【００４７】
また、ｅｃｇ１（ｆｉｌｔｅｒｅｄ）およびｅｃｇ２（ｆｉｌｔｅｒｅｄ）は、ｅｃｇ１およびｅｃｇ２にカットオフ周波数４Ｈｚのローパスフィルタ（図２のローパスフィルタ３０５、３０６に相当）を適用した心電図信号である。ｅｃｇ２（）はｅｃｇ１（ｆｉｌｔｅｒｅｄ）から推定したｅｃｇ２（ｆｉｌｔｅｒｅｄ）信号、ｅｒｒ（）は減算器３２０で得られる差信号、ｆｉｌｔｅｒｅｄ（）は差信号に呼吸信号抽出フィルタ３３０を適用して抽出した呼吸信号＾ｒｓである。
呼吸信号抽出フィルタ３３０は、残っているR波成分を除去するためのメジアンフィルタ （タップ数８０サンプル（０．６４秒））と、呼吸成分を残す程度の平滑化フィルタ（前後１５０サンプル（２．４秒）と、バイアス誤差を取り除くハイパスフィルタ（カットオフ周波数：０．０８Ｈｚ）との組み合わせとした。
【００４８】
また、ｒｅｓｐは、本実施形態による方法を評価するために、胸部に装着した努力呼吸センサで測定した呼吸信号を示している。なお、ｒｅｓｐとの比較を容易にするため、呼吸信号抽出フィルタ３３０の出力に−１を乗じたものを＾ｒｓとして図示している。
【００４９】
（）内のＡＲＸ、ＯＥ、ｒＯＥ、ＢＪはそれぞれ、推定信号生成部３１０で用いた推定モデルの種類を示している。なお、ｒＯＥ以外はデータ窓の長さを９０秒としたオフライン同定、ｒＯＥはオンライン同定を行った。同定に関する詳細は以下の通りである。
（次数）
ＡＲＸ（実施例１）  Ａｙ＝Ｂｕのとき  Ａ：７次，Ｂ：１０次
ＡＲＸ（実施例２）  Ａｙ＝Ｂｕのとき  Ａ：１０次，Ｂ：１０次
ＯＥ    ｙ＝Ｂ／Ｆｕ＋ｅのとき  Ｂ：２次，Ｆ：２次
ｒＯＥ                          Ｂ：２次，Ｆ：２次
ＢＪ    ｙ＝Ｂ／Ｆｕ＋Ｃ／Ｄｅのとき  Ｂ：２次，Ｆ：３次，Ｃ，Ｄ：５次
ＡＲＸモデル以外は実施例１，２とも同じ次数である。
（パラメータの同定法）
ＡＲＸ：最小２乗法
ＯＥ  ：予測誤差法
ｒＯＥ：正規化ＬＭＳ法
ＢＪ  ：予測誤差法
ステップサイズ：２．０ｅ−４
なお、ステップサイズの大きさは、パラメータの変化が拍動の影響を受けない程度の値として設定した一例である。
【００５０】
本実施形態の方法で得られた呼吸信号ｆｉｌｔｅｒｅｄ（）においても、努力呼吸センサで測定した呼吸信号ｒｅｓｐに類似した波形が得られていることが分かる。例えば健常者の例である図５〜図６における呼吸信号ｒｅｓｐの３０〜４０秒付近の区間に、無呼吸と思われる事象が現れているが、呼吸信号ｆｉｌｔｅｒｅｄ（）にも同様の事象が現れている。
【００５１】
なお、推定モデルにより呼吸信号ｆｉｌｔｅｒｅｄ（）とｒｅｓｐとの相関には差が見られ、複数の例による統計的な評価を行った結果、相関の高い順からＯＥ、ｒＯＥ、ＢＪ、ＡＲＸの順であった。ＯＥおよびｒＯＥモデルを用いて得られた推定信号の精度が良好であったのは、ＥＣ２＝（Ｓ２／Ｓ１）ＥＣ１＋ｒｓという関係と推定モデルとがうまく対応したことによるものと思われる。
【００５２】
さらに、図７〜図８に示したように、推定元となる心電図信号と、推定される心電図信号との関係を逆転させても、また、心房細動の患者の心電図信号であっても、良好な呼吸信号が得られた。
【００５３】
以上説明したように本実施形態によれば、誘導方向の異なる２つの心電図信号の一方から他方の推定信号を生成し、実測された心電図信号と推定信号との差信号の低周波成分を努力呼吸に関する信号（呼吸信号）として抽出する。そのため、心拍数の変動に依存せずに呼吸信号を得ることができ、心房細動の患者のように心拍数が変動する患者についても心電図信号からＳＡＳのスクリーニングを行うことが可能になる。
【００５４】
また、誘導方向の異なる心電図信号は、一般的な心電図の測定にも使用できる誘導法で得られる心電図信号から抽出することができるため、ホルタ心電計による心電図の測定結果などをＳＡＳのスクリーニング検査に利用することが可能である。従って、被検者の負担も従来の簡易ＰＳＧと比較して大幅に軽減できる。
【００５５】
なお、本発明に係る心電図解析装置は、一般的に入手可能な、パーソナルコンピュータのような汎用情報処理装置に、上述した動作を実行させるプログラム（アプリケーションソフトウェア）として実現することもできる。従って、このようなプログラムおよび、プログラムを格納した記憶媒体（ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ等の光学記録媒体や、磁気ディスクのような磁気記録媒体、半導体メモリカードなど）もまた本発明を構成する。
【図１】
0006624508.tif 000002
【図２】
0006624508.tif 000003
【図３】
0006624508.tif 000004
【図４】
0006624508.tif 000005
【図５】
0006624508.tif 000006
【図６】
0006624508.tif 000007
【図７】
0006624508.tif 000008
【図８】
0006624508.tif 000009