(57)【要約】【課題】ノイズ除去を行って極めて高精度で呼吸インピーダンスの連続測定を可能とする。【解決手段】ラウドスピーカ２１により口腔内に複数の異なる周波数から間引きして残った周波数成分のみを有するように周波数間引きされたオシレーション波であって呼気時と吸気時に応じてパルスの正負を異ならせて駆動するパルス信号によるオシレーション波による空気振動圧を加え、口腔内の圧力を検出し、呼吸の流量を検出し、この得られた信号をフーリエ変換手段３２にてフーリエ変換してスペクトルを得て、このフーリエ変換結果について間引きした周波数成分対応のスペクトルにより雑音として寄与する呼吸高周波成分を抽出手段３３にて求め、間引きにより残った周波数成分対応のスペクトルからこの呼吸高周波成分を減算してオシレーション波成分を抽出し、この抽出結果について周波数毎に圧力成分を流量成分で除算する演算を演算手段３４において行って呼吸インピーダンスを得る。【選択図】図１
2010264235.tif 000002
【特許請求の範囲】
【請求項１】
口腔内に空気振動圧を加えるための加圧手段と、
口腔内の圧力を検出する圧力検出手段と、
呼吸による流量を検出する流量検出手段と、
前記加圧手段を呼気時と吸気時に応じてパルスの正負を異ならせて駆動するパルス信号であって、複数の異なる周波数から間引きして残った周波数成分のみを有するように周波数間引きされたパルス信号であるオシレーション波による空気振動圧を生じさせる制御手段と、
前記加圧手段による加圧状態下において前記圧力検出手段と前記流量検出手段により得られる信号を得て、この得られた信号をフーリエ変換してスペクトルを得るフーリエ変換手段と、
このフーリエ変換手段による変換結果について間引きした周波数成分対応のスペクトルにより呼吸高周波成分を求め、間引きにより残った周波数成分対応のスペクトルから減算してオシレーション波成分を取り出す抽出手段と、
この抽出手段による抽出結果について周波数毎に圧力成分を流量成分で除算する演算手段と
を具備することを特徴とする呼吸インピーダンス測定装置。
【請求項２】
制御手段は周波数間引きとして、周期Ｔ1のパルス波を与えることにより、ｎ／Ｔ1（ｎ：整数、Ｔ1：実数）の周波数成分のみを有するオシレーション波による空気振動圧を生じさせることを特徴とする請求項１に記載の呼吸インピーダンス測定装置。
【請求項３】
制御手段は、加圧手段の入力信号と出力信号及び前記加圧手段の伝達関数を用いた逆演算に基づき所望の圧波形のオシレーション波が出力信号となるように前記加圧手段へ入力信号を与える信号入力手段を具備することを特徴とする請求項１または２に記載の呼吸インピーダンス測定装置。
【請求項４】
信号入力手段は、逆演算により得られた信号の各周波数成分に一定の値を加えるか、前記出力信号のオンセット部分にインパルスを加えた信号を逆演算して、得られた信号を前記加圧手段へ入力信号を与えることを特徴とする請求項３に記載の呼吸インピーダンス測定装置。
【請求項５】
信号処理区間Ｔ2（Ｔ2＝ｍＴ1（ｍは１以上の整数））において、信号処理を行うことを特徴とする請求項１乃至４のいずれか１項に記載の呼吸インピーダンス測定装置。
【請求項６】
呼吸インピーダンス測定装置により測定された呼吸インピーダンスに基づき表示装置に表示を行う呼吸インピーダンス表示方法において、
インピーダンス軸と周波数軸と時間軸とにより三次元に呼吸インピーダンスの値をとって三次元の表示を行い、
間引かれた周波数について補間処理を行って得られる呼吸インピーダンスを前記三次元の表示に含めた画像を作成して表示を行い、
呼気フローと吸気フローを検出するフロー検出手段により得られるフロー波形を前記インピーダンスと共に表示することを特徴とする呼吸インピーダンス表示方法。
【請求項７】
呼気の期間と吸気の期間を、画面の背景において縦方向に帯状に色分けして表示を行うことを特徴とする請求項５に記載の呼吸インピーダンス表示方法。
【請求項８】
複数回の呼吸における呼吸インピーダンスを経過時間毎に平均して平均値を求め、この平均値を三次元画像化して表示し、
所定数の周波数における前記平均値を、それぞれ線分として表示すると共に、前記所定数の周波数における最大値、最小値、最大最小差を求めて、文字により表示することを特徴とする請求項５または請求項６に記載の呼吸インピーダンス表示方法。
【請求項９】
時間軸方向の長さを、呼気と吸気との組を少なくとも二組繰り返す長さとして、画像を作成して表示を行うことを特徴とする請求項５乃至７のいずれか１項に記載の呼吸インピーダンス表示方法。
【請求項１０】
インピーダンス値の大小を、色の変化及び／または濃淡変化により表現した画像を作成して表示を行うことを特徴とする請求項５乃至８のいずれか１項に記載の呼吸インピーダンス表示方法。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、人などの呼吸インピーダンスを連続測定することが可能な呼吸インピーダンス測定装置及び呼吸インピーダンス表示方法に関するものである。
【背景技術】
【０００２】
従来、この種の装置は、呼吸器系に正弦波空気振動圧を負荷するための正弦波加圧装置と、呼吸器系の気流速度を検出するための気流速度検出器と、呼吸器系の気圧を検出するための気圧検出器と上記気流速度検出器および気圧検出器で検出した気流速度および気圧から呼吸抵抗を算出する抵抗演算部とを有したものが知られている。
【０００３】
上記従来の装置は、正弦波加圧装置が負荷する正弦波空気振動圧の信号を基準信号に変換するための基準信号変換器と、この基準信号変換器からの正弦波空気振動圧の基準信号により気流速度の信号を処理して上記基準信号と同じ周波数の成分のみを取り出すベクトル演算器とが設けられ、このベクトル演算器で得られた気流速度の信号と上記気圧検出器で検出した気圧の信号とから抵抗演算部で呼吸抵抗を算出するように構成したものである。
【０００４】
この装置は上記の通り、ベクトル演算器で得られた気流速度の信号と上記気圧検出器で検出した気圧の信号とから抵抗演算部で呼吸抵抗を測定するようにしているので、換気量が少なくまた換気数が多い呼吸であってもノイズを除去することができ、精度の高い呼吸抵抗の測定を行うことができる利点を有するものである（特許文献１参照）。
【０００５】
しかしながら、上記従来の装置によってもノイズ除去は十分とは言えず、更に高性能な呼吸インピーダンス測定装置の実現が求められている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開平０３−３９１４０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
本発明はこのような呼吸インピーダンス測定における現状に鑑みてなされたもので、その目的は、極めて高精度で呼吸インピーダンス測定が可能であり、呼吸系の疾患の判定に広く応用可能な呼吸インピーダンス測定装置及び呼吸インピーダンス表示方法を提供することである。
【課題を解決するための手段】
【０００８】
本発明に係る呼吸インピーダンス測定装置は、口腔内に空気振動圧を加えるための加圧手段と、口腔内の圧力を検出する圧力検出手段と、呼吸による流量を検出する流量検出手段と、前記加圧手段を呼気時と吸気時に応じてパルスの正負を異ならせて駆動するパルス信号であって、複数の異なる周波数から間引きして残った周波数成分のみを有するように周波数間引きされたパルス信号であるオシレーション波による空気振動圧を生じさせる制御手段と、前記加圧手段による加圧状態下において前記圧力検出手段と前記流量検出手段により得られる信号を得て、この得られた信号をフーリエ変換してスペクトルを得るフーリエ変換手段と、このフーリエ変換手段による変換結果について間引きした周波数成分対応のスペクトルにより呼吸高周波成分を求め、間引きにより残った周波数成分対応のスペクトルから減算してオシレーション波成分を取り出す抽出手段と、この抽出手段による抽出結果について周波数毎に圧力成分を流量成分で除算する演算手段とを具備することを特徴とする。
【０００９】
本発明に係る呼吸インピーダンス測定装置では、制御手段は周波数間引きとして、周期Ｔ1のパルス波を与えることにより、ｎ／Ｔ1（ｎ：整数、Ｔ1：実数）の周波数成分のみを有するオシレーション波による空気振動圧を生じさせることを特徴とする。
【００１０】
本発明に係る呼吸インピーダンス測定装置では、制御手段は、加圧手段の入力信号と出力信号及び前記加圧手段の伝達関数を用いた逆演算に基づき所望の圧波形のオシレーション波が出力信号となるように前記加圧手段へ入力信号を与える信号入力手段を具備することを特徴とする。
【００１１】
本発明に係る呼吸インピーダンス測定装置では、信号入力手段は、逆演算により得られた信号の各周波数成分に一定の値を加えるか、前記出力信号のオンセット部分にインパルスを加えた信号を逆演算して、得られた信号を前記加圧手段へ入力信号を与えることを特徴とする。
【００１２】
本発明に係る呼吸インピーダンス測定装置では、信号処理区間Ｔ2（Ｔ2＝ｍＴ1（ｍは１以上の整数））において、信号処理を行うことを特徴とする。
【００１３】
本発明に係る呼吸インピーダンス表示方法は、呼吸インピーダンス測定装置により測定された呼吸インピーダンスに基づき表示装置に表示を行う呼吸インピーダンス表示方法において、インピーダンス軸と周波数軸と時間軸とにより三次元に呼吸インピーダンスの値をとって三次元の表示を行い、間引かれた周波数について補間処理を行って得られる呼吸インピーダンスを前記三次元の表示に含めた画像を作成して表示を行い、呼気フローと吸気フローを検出するフロー検出手段により得られるフロー波形を前記インピーダンスと共に表示することを特徴とする。
【００１４】
本発明に係る呼吸インピーダンス表示方法では、呼気の期間と吸気の期間を、画面の背景において縦方向に帯状に色分けして表示を行うことを特徴とする。
【００１５】
本発明に係る呼吸インピーダンス表示方法では、複数回の呼吸における呼吸インピーダンスを経過時間毎に平均して平均値を求め、この平均値を三次元画像化して表示し、所定数の周波数における前記平均値を、それぞれ線分として表示すると共に、前記所定数の周波数における最大値、最小値、最大最小差を求めて、文字により表示することを特徴とする。
【００１６】
本発明に係る呼吸インピーダンス表示方法では、時間軸方向の長さを、呼気と吸気との組を少なくとも二組繰り返す長さとして、画像を作成して表示を行うことを特徴とする。
【００１７】
本発明に係る呼吸インピーダンス表示方法では、インピーダンス値の大小を、色の変化及び／または濃淡変化により表現した画像を作成して表示を行うことを特徴とする。
【発明の効果】
【００１８】
本発明によれば、口腔内に周波数間引きされたオシレーション波による空気振動圧を加え、口腔内の圧力を検出し、呼吸の流量を検出し、この得られた信号をフーリエ変換してスペクトルを得て、このフーリエ変換結果について間引きした周波数成分対応のスペクトルにより雑音として寄与する呼吸高周波成分を求め、間引きにより残った周波数成分対応のスペクトルからこの呼吸高周波成分を減算してオシレーション波成分を抽出し、この抽出結果について周波数毎に圧力成分を流量成分で除算する演算を行って呼吸インピーダンスを得るので、呼吸高周波成分を確実に除去したオシレーション波成分を用いて呼吸インピーダンスを得ることができる上に、オシレーション波による空気振動圧にパルス信号を採用し、該パルスの正負を呼気時と吸気時に応じて異ならせて駆動するので、吸気と呼気で生理学的に同一条件において呼吸抵抗を測定することができ、適切な測定により極めて高精度で呼吸インピーダンス測定が可能となる。
【００１９】
本発明によれば、周期Ｔ1のパルスを与えることにより、ｎ／Ｔ1（ｎ：整数Ｔ1：実数）の周波数成分のみを有するオシレーション波による空気振動圧を生じさせるので、間引いた周波数成分対応のスペクトルにより呼吸高周波成分を求め、間引きにより残った周波数成分対応（ｎ／Ｔ1周波数以外）のスペクトルからこの呼吸高周波成分を減算して呼吸高周波成分を確実に除去し、極めて高精度で呼吸インピーダンス測定が可能となる。
【００２０】
本発明によれば、異なる複数周波数の正弦波を複数合成することにより、周波数成分間引きされたオシレーション波による空気振動圧を生じさせるので、間引きした周波数成分対応のスペクトルは呼吸高周波成分のみが含まれる。従って、オシレーション信号に含まれる呼吸高周波成分を、この呼吸高周波成分信号から推定できる。この推定した信号を用いて呼吸高周波成分を確実に除去し、極めて高精度で呼吸インピーダンス測定が可能となる。
【００２１】
本発明によれば、加圧の入力信号と出力信号及び加圧実行部分の伝達関数を用いた逆演算に基づき所望波形のオシレーション波が出力信号となるように加圧実行部分へ入力信号を与えるので、所望の圧波形のオシレーション波を用いて測定を行うことが可能となり、極めて高精度で呼吸インピーダンス測定が可能となる。
【００２２】
本発明によれば、逆演算により得られた信号の各周波数成分に一定の値を加えるか、前記出力信号のオンセット部分にインパルスを加えた信号を逆演算して、得られた信号を入力信号とするので、逆演算の結果の信号波形を安定させ、これによって所望波形のオシレーション波を用いて測定を行うことが可能となり、極めて高精度で呼吸インピーダンス測定が可能となる。
【００２３】
本発明によれば、信号処理区間Ｔ2（Ｔ2＝ｍＴ1（ｍは１以上の整数））において、信号処理を行うので、間引きにより残った周波数成分対応のスペクトルと、間引きされた周波数成分対応のスペクトルとを適切に得て所要の処理を行うことができる。
【００２４】
本発明に係る呼吸インピーダンス表示方法によれば、呼吸インピーダンス測定装置により測定された呼吸インピーダンスに基づき表示装置に表示を行う呼吸インピーダンス表示方法において、インピーダンス軸と周波数軸と時間軸とにより三次元に値をとって三次元の表示を行い、間引かれた周波数について補間処理を行って得られる呼吸インピーダンスを前記三次元に値をとって表示する場合に含めた画像を作成して表示を行い、呼気フローと吸気フローを検出するフロー検出手段により得られるフロー波形を前記インピーダンスと共に表示するので、補間処理の結果についても画像化して表示するため、インピーダンス値の変化がきめ細かく滑らかに表示され、周波数全体についてインピーダンスの把握を適切に行うことが可能となり、呼気フローと吸気フローによるフロー波形が共に表示されており、呼吸系の疾患の判定に広く応用可能である。
【００２５】
本発明に係る呼吸インピーダンス表示方法では、呼気の期間と吸気の期間を、画面の背景において縦方向に帯状に色分けして表示を行うので、呼吸インピーダンスの表示において呼気の期間と吸気の期間が一目瞭然となり、呼吸系の疾患の判定などに便利である。
【００２６】
本発明に係る呼吸インピーダンス表示方法では、複数回の呼吸における呼吸インピーダンスを経過時間毎に平均して平均値を求め、この平均値を三次元画像化して表示し、所定数の周波数における前記平均値を、それぞれ線分として表示すると共に、前記所定数の周波数における最大値、最小値、最大最小差を求めて、文字により表示するので、各呼吸毎に異なる呼吸インピーダンスを平均化して文字により取得でき、比較などに便利である。
【００２７】
本発明に係る呼吸インピーダンス表示方法は、時間軸方向の長さを、呼気と吸気との組を少なくとも二組繰り返す長さとして、画像を作成して表示を行うので、突発的な変化ではなくある程度のスパンをもった観測が可能となり、適切な観測が確保される。
【００２８】
本発明に係る呼吸インピーダンス表示方法は、インピーダンス値の大小を、色の変化または濃淡変化により表現した画像を作成して表示を行うので、インピーダンス値の大小を一目瞭然に識別することが容易に可能となり、呼吸インピーダンスによる各種の研究、検査、更には患者への視角的な説明（アピール）に極めて役立つことが期待される。
【図面の簡単な説明】
【００２９】
【図１】本発明の実施例に係る呼吸インピーダンス測定装置の構成図。
【図２】本発明の実施例に係る呼吸インピーダンス測定装置に用いるオシレーション波である三角パルス波の一例を示す図。
【図３】本発明の実施例に係る呼吸インピーダンス測定装置に用いるオシレーション波であるハニングパルス波の一例を示す図。
【図４】本発明の実施例に係る呼吸インピーダンス測定装置において、パルスの正負を呼気時と吸気時に応じて異ならせて駆動する例を示す図。
【図５】本発明の実施例に係る呼吸インピーダンス測定装置に用いるオシレーション波を逆演算により生成する過程を説明するための図。
【図６】本発明の実施例に係る呼吸インピーダンス測定装置によりオシレーション波である三角パルス波を用いて呼吸インピーダンスを求める過程を示す図。
【図７】本発明の実施例に係る呼吸インピーダンス測定装置に採用されているフィルタの周波数特性を示す図。
【図８】本発明の実施例に係る呼吸インピーダンス測定装置によりオシレーション波である三角パルス波を用いて呼吸インピーダンスを求める過程を示す図。
【図９】本発明の実施例に係る呼吸インピーダンス測定装置によりオシレーション波である三角パルス波を用いて呼吸インピーダンスを求める過程を示す図。
【図１０】本発明の実施例に係る呼吸インピーダンス測定装置によりオシレーション波である三角パルス波を用いて呼吸インピーダンスを求める過程を示す図。
【図１１】本発明の実施例に係る呼吸インピーダンス測定装置により得られた呼吸インピーダンスを示す図。
【図１２】本発明の手法を用いない呼吸インピーダンス測定装置により得られた呼吸インピーダンスを示す図。
【図１３】本発明の実施例に係る呼吸インピーダンス測定装置を用いてＣＯＰＤ患者（stage１）について呼吸インピーダンスを表示した一例を示す図。
【図１４】本発明の実施例に係る呼吸インピーダンス測定装置を用いて正常者について呼吸インピーダンスを表示した一例を示す図。
【図１５】本発明の実施例に係る呼吸インピーダンス測定装置を用いてＣＯＰＤ患者（stage２）について呼吸インピーダンスを表示した一例を示す図。
【図１６】本発明の実施例に係る呼吸インピーダンス測定装置を用いて喘息患者について呼吸インピーダンスを表示した一例を示す図。
【図１７】本発明の実施例に係る呼吸インピーダンス測定装置を用いて呼吸同期毎に加算平均して平均値を求め、この平均値を三次元画像化して表示した一例を示す図。
【図１８】図１７の三次元画像の表示について、所定数の周波数における平均値を、それぞれ線分として表示した一例を示す図。
【図１９】図１８の表示について、所定数の周波数における最大値、最小値、最大最小差を求めて、文字により表示した一例を示す図。
【発明を実施するための形態】
【００３０】
以下、添付図面を参照して本発明に係る呼吸インピーダンス測定装置及びその測定方法の実施例を説明する。図１に、本発明に係る呼吸インピーダンス測定装置の実施例構成図を示す。この呼吸インピーダンス測定装置は、人の口腔に先端が取り付けられ、呼吸流が流れるチューブ１１と、チューブ１１に取り付けられ口腔内の圧力を検出する圧力検出手段を構成する圧力センサ１２と、圧力センサ１２と同位置において呼吸による流量を検出する流量検出手段を構成する流量センサ１３と、口腔内に空気振動圧を加えるための加圧手段を構成するラウドスピーカ２１と、コンピュータ３０とを主な構成要素とする。
【００３１】
圧力センサ１２の出力信号はアンプ１４により増幅され、Ａ／Ｄ変換器１５によりディジタル化されてコンピュータ３０に取り込まれる。また、流量センサ１３の出力信号はアンプ１６により増幅され、Ａ／Ｄ変換器１７によりディジタル化されてコンピュータ３０に取り込まれる。
【００３２】
コンピュータ３０には、制御手段３１、フーリエ変換手段３２、抽出手段３３、演算手段３４が備えられている。また、制御手段３１は、信号入力手段３５を備えている。制御手段３１は、加圧手段であるラウドスピーカ２１を駆動する信号を出力し、奇数周波数成分または偶数周波数成分のみを有するオシレーション波による空気振動圧を生じさせるものである。制御手段３１の出力はＤ／Ａ変換器２２によりディジタル化されてドライバ２３へ送られ、ドライバ２３がラウドスピーカ２１を駆動して、口腔内に空気振動圧が加えられる。
【００３３】
上記において、制御手段３１は、Ｔ1秒周期のパルス波を与えることにより、ｎ／Ｔ1（ｎ：整数、Ｔ1：実数）の周波数成分を有するオシレーション波による空気振動圧を生じさせるものである（周波数間引き）。パルス波としては様々な波形を用いることができる。例えば三角パルスは、図２（ａ）に示すようにベースレベルの時間幅が25ｍｓ程度のものである。この三角パルスを例えばＴ1=0.5 秒周期で出力すると2,4,6,8Hz,・・・のスペクトルを有する三角パルス波を与えることができる（図２（ｂ））。また、上記の三角パルスを例えばＴ1=0.333 秒周期で出力すると3,6,9,12Hz, ・・・のスペクトルを有する三角パルス波を与えることができる。勿論、Ｔ1=0.25秒周期で出力するようにしても良く、この場合には4,8,12,16Hz, ・・・のスペクトルを有する三角パルス波を与えることができる。
【００３４】
また、他のパルス波の例としてのハニングパルスは、図３に示すようにベースレベルの時間幅が25ms程度のものである。三角パルス波の場合と同様にして作成され、出力される。
【００３５】
制御手段３１に備えられている信号入力手段は、ラウドスピーカ２１の入力信号と出力信号及びラウドスピーカ２１の伝達関数を用いた逆演算に基づき所望波形のオシレーション波が出力信号となるようにラウドスピーカ２１へ入力信号を与えるものである。
【００３６】
上記において、パルス波の正負とフロー（呼吸流）の方向との関係は生理学的に重要である。今、フロー曲線において負の領域にある呼気フローの方向（下に凸）とパルス波が同一方向で負のパルス（下に凸）とする。この場合、フローは加速され、その速度は大きくなる。このため、気管支内圧がパルス波によるオシレーション波の影響で減衰し、気管支の径は短縮する。結果として得られる呼吸抵抗は増大する。一方、この負のパルスを吸気においても使用した場合、吸気フローはフロー曲線の正の領域にあり、方向は上に凸であり、呼気と逆向きであるから上記負のパルス波はフローを減速する方向に働き、気管支内圧は高くなり気管支径は拡張する。結果として得られる呼吸抵抗は小さくなる。
【００３７】
斯して、呼気と吸気に拘りなく同一方向に凸のパルスを用いた場合には、測定される呼吸抵抗はオシレーション波として加えられるパルスの呼吸フローに対する方向性に依存してしまい、吸気と呼気で生理学的に異なる条件において呼吸抵抗の測定をしたことになり、正しい測定とならない。これを避けるために、制御手段３１は、呼気と吸気とに応じてパルスの正負を異ならせてパルスを出力する。
【００３８】
その第１の手法は、パルス波の方向（上に凸または下に凸）を吸気、呼気のいずれにおいてもフローと同一方向とする。即ち、フロー曲線（図４（ｃ））が負の領域となる呼気の場合には、負のパルス波を用い、フロー曲線が正の領域となる吸気の場合には、正のパルス波を用いるようにする。図４（ａ）に、このようにしてパルス波の正負を切り換えた場合を示す。
【００３９】
その第２の手法は、パルス波の方向（上に凸または下に凸）を吸気、呼気のいずれにおいてもフローと逆方向とする。即ち、フロー曲線（図４（ｃ））が負の領域となる呼気の場合には、正のパルス波を用い、フロー曲線が正の領域となる吸気の場合には、負のパルス波を用いるようにする。図４（ｂ）に、このようにしてパルス波の正負を切り換えた場合を示す。
【００４０】
制御手段３１は、流量センサ１３の出力に基づき、呼気と吸気についてその切り替わり時点を検出する。すなわち、流量センサ１３の出力は、呼吸に対応して基本的に図４（ｃ）に示されるように変化するものであるから、上記流量センサ１３の出力信号のゼロクロス点を境にしてパルス波の正負を切り換える。
【００４１】
具体的に、例えば三角パルスを用いて測定を行う動作について詳細に説明する。以下では、正のパルス波についての説明であるが、負のパルス波についても正のパルス波と同様の処理が行われる。図５（ａ）に示されるような三角パルスを入力してラウドスピーカ２１を駆動した場合には、ラウドスピーカ２１の出力信号は、図５（ｂ）に示されるようなゼロレベルの上下に極大点を有する信号となる。そこで、図５（ｃ）に示されるようなモデルを考える。ラウドスピーカ２１の伝達関数をＨ（ω）、入力信号をＸ（ω）、出力信号をＹ（ω）とすると、以下のようになるので、逆変換によりｘ’（ｔ）を求めて駆動信号とする。
【００４２】
【数１】
2010264235.tif 000003
【００４３】
求めたＹ（ω）は、高い周波数までの成分を持たないので、（式１）から得られるｘ’（ｔ）は不安定となる。このため（式２）のように分母に定数Ａ0 を加えたものをフーリエ逆変換してｘ’（ｔ）を求めて駆動信号とする。これはまた、図５（ｂ）に示されるようなラウドスピーカ２１の出力信号に対し、図５（ｄ）に示されるようにオンセット部分にインパルスを加えた信号を逆演算して図５（ｅ）に示される信号ｘ’（ｔ）を得ることも出来る。
【００４４】
以上の説明では、三角パルスの場合を説明したが、ハニングパルスについても同様にして逆演算により信号を得て、この信号によりラウドスピーカ２１を駆動することができる。
【００４５】
パルス波のいずれを用いるかに関しては、図示しないキーボードなどによりコンピュータ３０へ指示を与えることができ、これに応じて制御手段３１が選択した信号波形を出力する。
【００４６】
コンピュータ３０に備えられているフーリエ変換手段３２、抽出手段３３、演算手段３４について説明する。フーリエ変換手段３２は、ラウドスピーカ２１が上記のようにして駆動されることによって生じる口腔内の加圧状態下において、圧力センサ１２と流量センサ１３により信号を得て、この得られた信号をフーリエ変換してスペクトルを得るものである。フーリエ変換手段３２の前段には、ＣＩＣフィルタ３６が設けられ、圧力センサ１２と流量センサ１３により得られる呼吸信号とオシレーション成分との分離が行われる。また、フーリエ変換手段３２は処理の前に必要であれば、ハニング窓により信号の取り出しを行う。
【００４７】
抽出手段３３は、フーリエ変換手段３２による変換結果について間引きした周波数成分対応のスペクトルにより呼吸高周波成分を求め、間引きにより残った周波数成分対応のスペクトルから減算してオシレーション波成分を取り出すものである。周波数間引きに対応して、フーリエ変換手段３２により得られたスペクトルについてｎ／Ｔ1（ｎ：整数）の周波数成分を除く他の周波数対応のスペクトルにより呼吸高周波成分を求め、間引きにより残った周波数成分（ｎ／Ｔ1の周波数成分）対応のスペクトルから減算してオシレーション波成分を取り出すものである。
【００４８】
演算手段３４は、抽出手段３３による抽出結果について周波数毎に圧力成分を流量成分で除算することにより呼吸インピーダンスを算出するものである。即ち、呼吸インピーダンスをＺ（ω）、口腔内圧力のオシレーション波成分をＰ（ω）、流量のオシレーション波成分をＦ（ω）とし、呼吸インピーダンスＺ（ω）が抵抗成分Ｒ（ω）とリアクタンス成分Ｘ（ω）からなるものとすると、次の式により呼吸インピーダンスをＺ（ω）が求められる。
【００４９】
【数２】
2010264235.tif 000004
【００５０】
演算手段３４により求められた呼吸インピーダンスＺ（ω）は、コンピュータ３０に接続されたＬＣＤなどの表示部４０用の表示信号とされて表示部４０へ出力され、表示がなされる。
【００５１】
以上の通りに構成された呼吸インピーダンス測定装置による動作を説明する。この例では、三角パルス波が選択されて測定動作が開始される。制御手段３１及び信号入力手段３５により、逆演算された波形によりスピーカ２１がＴ1秒周期で（例えば0.5秒間隔で）駆動される。
【００５２】
このとき、圧力センサ１２と流量センサ１３により得られる信号の波形は、いずれも図６（ａ）または図６（ｂ）に示されるように、呼吸信号に三角パルス波が重畳した波形となっており、これがＣＩＣフィルタ３６を通過させられて呼吸波とオシレーション波（三角パルス波）の分離が行われる。図７にＣＩＣフィルタ３６の周波数特性を示す。ＣＩＣフィルタ３６により位相の変化なく、分離を行うことができる。但し、呼吸信号には高周波成分（オシレーション信号と同一の周波数帯）が含まれるために完全には分離できない。
【００５３】
ＣＩＣフィルタ３６による分離の後には図６（ｃ）に示されるように、負のパルスが反転されてすべて正のパルスとされたオシレーション波について、二つの三角パルスの中間点から１秒間（信号処理区間Ｔ2）の間を取り出し、信号処理に用いる。次に図８に示されるように、Ｔ2秒の区間を取り出した各パルスにハニング窓による処理を行い、パルスの取出しを行う。
【００５４】
ハニング窓による処理に続き、フーリエ変換手段３２によるフーリエ変換が行われてスペクトルが得られる。このとき、得られるスペクトルについては、例えば、Ｔ1＝0.5秒、Ｔ2＝1.0秒の場合(一般的には、Ｔ2＝ｍＴ1（ｍは１以上の整数））の場合、図９に示されるように、間引きした周波数成分対応の１、３、５、・・・の奇数周波数のスペクトルにはオシレーション波成分が含まれない呼吸信号スペクトルとなっている。また、間引きにより残った周波数成分対応の２、４、６、・・・の偶数周波数のスペクトルにはオシレーション波成分と呼吸信号成分が含まれている。
【００５５】
そこで、抽出手段３３では、図１０に示すように上記偶数周波数のスペクトルから、奇数周波数のスペクトルより推定された雑音成分を減算し、オシレーション波成分を取り出す。
【００５６】
上記抽出手段３３による処理によって、従来、呼吸信号に含まれていないとされていた３Ｈｚ以上の呼吸高周波信号が除去され、高精度な呼吸インピーダンス測定が可能となっている。次に演算手段３４は、抽出手段３３による抽出結果について周波数毎に式（３）により示した通り、圧力成分を流量成分で除算することにより呼吸インピーダンスを算出する。算出された呼吸インピーダンスの表示信号が作成され、表示部４０へ出力され、表示がなされる。
【００５７】
このようにして測定され表示された健常者の呼吸インピーダンスを図１１に示す。また、図１２に呼吸高周波信号の除去を行わない場合の呼吸インピーダンスを示す。これらの図では、横軸は１目盛りが１Ｈｚの周波数軸であり、縦軸がインピーダンスである。斜めの軸が時間軸であり、図の上側に純抵抗分を表示し、図の下側にリアクタンス分を表示している。ここでは０．５秒毎に三角パスル波を与え続けることにより、新たなインピーダンスの表示が次々に現れ更新されることによりインピーダンスの連続測定が行われる。この図１１と図１２から明らかな通り、本実施例の装置によればノイズが除去され高精度な呼吸インピーダンス測定が可能となっていることが分かる。なお、上記抽出手段３３による減算処理から明らかな通り、減算により残る成分は、間引きにより残った周波数成分対応の２、４、６、・・・の偶数周波数の成分であり、間引きした周波数成分対応の１、３、５、・・・の奇数周波数の成分は存在しない。そこで、演算手段３４が補間処理を行い、存在しない成分に関しても呼吸インピーダンス表示が可能となっている。
【００５８】
本発明の実施例では、演算手段３４が表示装置に表示を行うための画像を作成して、表示を行うことにより、呼吸インピーダンス表示方法が実現される。即ち、演算手段３４は、自らが前述の通りに算出した呼吸インピーダンスについて、例えば、各周波数を画面の奥側から手前側に値をとるように座標を定め（Ｙ軸）、各周波数毎に抵抗成分Ｒrsを取り出し、これを表示装置の画面における高さ方向にプロットし（Ｚ軸）、測定時間を画面の右方向として（Ｘ軸）、図１３に示すような三次元画像を作成して表示装置に表示する。即ち、インピーダンス軸と周波数軸と時間軸とにより三次元に値をとって三次元の表示を行うものである。
【００５９】
上記の画像作成に際しては、間引かれた周波数について補間処理を行って得られる呼吸インピーダンスを上記三次元に値をとって表示する場合に含めた画像を作成して表示を行う。例えば、奇数の周波数を間引いた場合には、間引いた奇数に隣接する偶数の周波数に対応する２つの呼吸インピーダンス値が求まっているから、この２つの呼吸インピーダンス値の平均を求め、これを間引いた周波数対応の呼吸インピーダンス値とする。このようにして、補間処理の結果についても画像化して表示するため、呼吸インピーダンス値の変化がきめ細かく滑らかに表示され、周波数全体について呼吸インピーダンスの把握を適切に行うことが可能となる。
【００６０】
信号処理の時間分解能は０．５秒であり、図１３〜図１６に示すように、時間軸方向の長さを、呼気と吸気との組を少なくとも二組繰り返す長さとして、画像を作成して表示を行う。図１３〜図１６の例では、呼気と吸気との組を六組繰り返す長さとしている。また、この場合において、呼気の期間と吸気の期間を、画面の背景において縦方向に帯状に色分けして表示を行う。更に、呼気フローと吸気フローを検出するフロー検出手段である流量センサ１３や圧力センサ１２により得られるフロー波形（呼吸信号波形）ＢＳを上記呼吸インピーダンス値と共に表示する。
【００６１】
更に、呼吸インピーダンス値の大小を、色の変化及び（または）濃淡変化により表現した画像を作成して表示を行う。図１３〜図１６では、抵抗値Ｒrsについて、図１３〜図１６の下側に示すカラースケールＣＬによる色付けを行い、画像を作成して表示している。
【００６２】
以上の各処理により得られる画像を表示するので、被検者は安静呼吸を繰り返しているだけで、自動的に図１３〜図１６に示されるような画像が時系列に作成されて表示され、しかも呼吸インピーダンスの変化を間引いた周波数対応の部分を含めて、色の変化及び（または）濃淡変化により表現した画像として目視観察することができる。
【００６３】
呼吸抵抗は計測される呼吸インピーダンスの実部であり、その虚部はリアクタンスである。本実施例に係る装置ではこの両者を三次元カラー表示し、さらにフロー波形ＢＳと対比して見ることができるようにしている（図１３〜図１６）。図１４は正常者について測定を行った結果を示し、図１５はＣＯＰＤ患者（慢性閉塞性肺疾患）について測定を行った結果を示し、図１６は喘息患者について測定を行った結果の表示例を示している。
【００６４】
図１３は、ＣＯＰＤ患者について測定を行った結果の表示例であるが吸気区間における呼吸抵抗が緑色を呈して表示されており、これは正常者と同じ状態の色により表示されている。よって、該患者は吸気期間において呼吸困難を感じないと思料される。これによりＣＯＰＤのstage1 の患者と判断できる。図１５もＣＯＰＤ患者について測定を行った結果の表示例であるが、吸気期間において赤色に近い色によって抵抗が高いことが表示されている。よって、重症でありＣＯＰＤにおけるstage2の患者であると判断できる。このように本表示法（３Ｄカラーグラフィック表示）により呼吸インピーダンスの周波数依存性と呼吸周期依存性が同時に表現され、呼吸系に関する正常、異常の判定はもとより、異常の程度の判定も容易に行い得ることが理解できる。
【００６５】
本実施例により表示がなされた画面を示す図１３は、前述の通りＣＯＰＤ患者の呼吸インピーダンスを三次元表示した例である。図１３の表示例では、吸気の呼吸抵抗は緑色であり健常者と同一色で表示されているのであるが、呼気の部分においては赤色となって表示されており呼吸抵抗が非常に高い値となっていることが分る。よって、該図１３により表示された呼吸インピーダンスを呈する患者は、呼気の場合に困難を伴うことがこの実測データ（表示）から確認することができる。一方、斯様な呼気の場合に困難を伴う患者について気管支のＣＴスキャンを撮り、気管支の径を調べて見ると吸気時において気管支が拡張し、呼気時においては狭窄していることが確認されるという知見が、例えば、文献（Kurosawa, et al. N Engl J Med.350:1036, 2004)などに紹介されている。よって、この呼吸抵抗の呼吸周期依存性は気管支の径の変化や狭窄などをかなり反映していると解釈でき、これが本実施例の装置により確認できるようになり、臨床上非常に有益である。
【００６６】
ところで、臨床ではこれらの変化をできる限り少数のパラメータで表現し数値化することの要求がなされている。そこで、本実施例の装置においては、演算手段３４が複数回の呼吸における呼吸インピーダンスを経過時間毎に平均して平均値を求め、この平均値を三次元画像化して表示装置に表示を行う。即ち、呼吸インピーダンスは１呼吸毎に多少変化するため、呼吸毎の加算平均を行って三次元画像化して表示装置に図１７に示すように表示を行う。これにより、一呼吸に対する平均的な呼吸インピーダンスを視覚的に把握することができ、比較検討などの場合に好適である。
【００６７】
更に、演算手段３４は、所定数の周波数における上記平均値を、それぞれ線分として表示すると共に、上記所定数の周波数における最大値、最小値、最大最小差を求めて、文字により表示装置に表示する。この表示例を図１８、図１９に示す。具体的には、上記において得られた加算平均パターンを用いて、その代表値として5[Hz]、20[Hz]、共振周波数（resと表現）の抵抗値の時間変化を算出したものである。なお、ここで共振周波数とはリアクタンスが零（式（２）でX(ω)=0）となる周波数の値である。
【００６８】
実際の表示にあっては、緑色のR5が5[Hz]を示し、赤色のR20が20[Hz]を示し、青色のRresが共振周波数における抵抗値を示している。図１９のテーブルはそれらの抵抗値の呼吸周期における最小値(min)、最大値（max）、さらに両者の差（sub）を現すものである。これらの数値は生体の生理学的意味合いを表現しており、臨床上極めて有益である。図１９のテーブルにおいて、Fresは共振周波数を表すものである。
【００６９】
演算手段３４は、例えば、図１７に示される一呼吸に対する平均的な呼吸インピーダンスを一画面の左半分側に表示し、図１８のグラフを一画面の右半分側上部に表示し、図１９のテーブルを一画面の右半分側下部に表示するように機能する。このような表示により、呼吸インピーダンスの変動傾向を平均表示やグラフより把握し、実際の数値についてはテーブルを参照して確認できるので、臨床の場などにおいて有効に利用され得るものである。
【符号の説明】
【００７０】
１１  チューブ
１２  圧力センサ
１３  流量センサ
２１  ラウドスピーカ
３０  コンピュータ
３１  制御手段
３２  フーリエ変換手段
３３  抽出手段
３４  演算手段
３５  信号入力手段
３６  ＣＩＣフィルタ
４０  表示部
【図１】
2010264235.tif 000005
【図２】
2010264235.tif 000006
【図３】
2010264235.tif 000007
【図４】
2010264235.tif 000008
【図５】
2010264235.tif 000009
【図６】
2010264235.tif 000010
【図７】
2010264235.tif 000011
【図８】
2010264235.tif 000012
【図９】
2010264235.tif 000013
【図１０】
2010264235.tif 000014
【図１１】
2010264235.tif 000015
【図１２】
2010264235.tif 000016
【図１３】
2010264235.tif 000017
【図１４】
2010264235.tif 000018
【図１５】
2010264235.tif 000019
【図１６】
2010264235.tif 000020
【図１７】
2010264235.tif 000021
【図１８】
2010264235.tif 000022
【図１９】
2010264235.tif 000023