(57)【特許請求の範囲】
【請求項１】
腸が発する生体音を検知する生体音検知手段と、
上記生体音の周波数スペクトルを算出する周波数スペクトル算出手段と、
上記生体音の周波数スペクトルと、複数の蠕動音の標準周波数スペクトルのそれぞれとを個別にマッチングすることによって、複数のマッチング係数を算出するマッチング係数算出手段と、
上記複数のマッチング係数を演算処理することによって、上記生体音が蠕動音であるか否かを判定する蠕動音判定手段と、
を備えることを特徴とする蠕動音検出装置。
【請求項２】
上記複数の蠕動音の標準周波数スペクトルは、それぞれが特定の発生モードに起因して発生する蠕動音の周波数スペクトルであることを特徴とする請求項１に記載の蠕動音検出装置。
【請求項３】
上記蠕動音判定手段は、上記複数のマッチング係数のうち最大のマッチング係数が、規定の閾値よりも大きいときに、上記生体音が蠕動音であると判定することを特徴とする請求項１または２に記載の蠕動音検出装置。
【請求項４】
上記蠕動音判定手段は、上記複数のマッチング係数を全て合算した値が、規定の閾値よりも大きいときに、上記生体音が蠕動音であると判定することを特徴とする請求項１または２に記載の蠕動音検出装置。
【請求項５】
請求項１〜４のいずれか１項に記載の蠕動音検出装置を動作させるプログラムであって、コンピュータを上記の各手段として機能させるためのプログラム。
【請求項６】
請求項５に記載のプログラムを記録しているコンピュータ読取り可能な記録媒体。
【請求項７】
腸が発する生体音を検知する生体音検知工程と、
上記生体音の周波数スペクトルを算出する周波数スペクトル算出工程と、
上記生体音の周波数スペクトルと、複数の蠕動音の標準周波数スペクトルのそれぞれとを個別にマッチングすることによって、複数のマッチング係数を算出するマッチング係数算出工程と、
上記複数のマッチング係数を演算処理することによって、上記生体音が蠕動音であるか否かを判定する蠕動音判定工程と、
を備えることを特徴とする蠕動音検出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、腸が発する音が蠕動音であるか否かを判定する蠕動音検出装置、蠕動音検出方法、プログラム、および記録媒体に関する。
【背景技術】
【０００２】
これまで、消化器の活性度を評価する定量的な評価装置および評価方法は実用化されていない。そのため、消化器（例えば腸）が活性であるか否かの評価は、医師によって行われている。具体的には、医師は聴診器を用いて腹部の蠕動音を聞き、その蠕動音に基づいて腸が活性であるか否かを評価する。言い換えれば、腸の活性度は、評価者である医師の経験および主観に基づいて評価される。
【０００３】
このように、消化器の活性度を評価するためには、(i)当該消化器の発する音を検知する、(ii)検知した音を音響解析し評価する、という２つのステップを経る。
【０００４】
特許文献１には、消化器系が発する蠕動音をフーリエ変換することによって得られるフーリエ変換スペクトル（周波数スペクトル）に基づいて蠕動音を記録した時間区間を自動的に検出する技術が記載されている。この技術では、周波数スペクトルの１００〜１０００Ｈｚの範囲に有意な大きさのピークを持つ時間区間を、蠕動音を記録している区間だと考えている。
【０００５】
生体の発する音を検出および解析する試みとして、特許文献２には、入力された音声に対して、与えられた辞書の中から、その音声に対応する単語などを選択する音声認識手法が記載されている。当該技術において、対応する単語を選択する際に、ＤＰ(Dynamic Programming)マッチングを用いることができることが記載されている。ＤＰマッチングでは、音声認識の対象となる単語に対応した標準パターンをあらかじめ作成しておく。入力された音声を音響分析して得られる特徴量と、標準パターンとをマッチングさせ、例えば、音声の特徴量と最も類似する標準パターンに対応する単語が、音声の認識結果とされる。検出される音声パターンと、比較対象となる音声パターン（標準パターン）とをパターンマッチングする際、標準パターンは１つであることが一般的である。
【０００６】
特許文献３には、心音計で検出された弁閉鎖時の弁音波形の周波数分布として、Wigner分布を採用し、さらに、Wigner分布について各時刻毎に１次モーメントを求める技術が記載されている。
【０００７】
特許文献４に記載されている生体音検出処理装置は、生体音検出データ（呼吸音検出データ）をＦＦＴ処理し、振幅スペクトル、位相スペクトルおよびパワースペクトルを算出する。さらに、生体音検出処理装置はパワースペクトルから局所平均値および局所分散値を算出する。この局所分散値の大きさによって、振幅スペクトルは、正常呼吸音に対応するもの、または、連続性ラ音によるものに分類される。
【０００８】
特許文献５に記載されている咳嗽検出装置は、密着マイクロフォンを介して取得した音信号から、第１および第２の周波数帯域の音信号を抽出する。咳嗽検出装置は第１帯域信号から候補咳嗽を決定する。さらに、候補咳嗽と第２帯域信号との対応関係から、候補咳嗽が咳嗽であるか否かの判定を行う。
【先行技術文献】
【特許文献】
【０００９】
【特許文献１】特開２０００−２６２５２３号公報（２０００年９月２６日公開）
【特許文献２】特開平９−１６０５８５号公報（１９９７年６月２０日公開）
【特許文献３】特開平６−９０９１３号公報（１９９４年４月５日公開）
【特許文献４】特開２００４−３５７７５８号公報（２００４年１２月２４日公開）
【特許文献５】特開２００９−２３３１０３号公報（２００９年１０月１５日公開）
【発明の概要】
【発明が解決しようとする課題】
【００１０】
腸における蠕動運動の発生モードは複数ある。また、蠕動運動に伴い腸から発せられる生体音は個人差が大きい。これらに起因して、上記生体音の周波数スペクトルは大きなばらつきを有することがある。
【００１１】
したがって、周波数スペクトルの特定の範囲内に有意な大きさのピークがあるか否かによって蠕動音であるか否かを判定する特許文献１の技術では、蠕動音を精度良く検出することは困難である。また、特許文献２に記載される技術である検知される音声パターンと、１つの標準パターンとのパターンマッチングによっても、所望する音声を精度良く評価することは困難である。
【００１２】
一方、特許文献３〜５に記載されている技術は、音響解析のために複雑な演算処理を行う。そのため、音響解析に要する時間が長くなるという課題を有している。同時に、複雑な演算を行うために、高性能な処理装置が求められるという課題を有している。すなわち、特許文献３〜５に記載されている技術は装置が高価かつ大型になる、および、音響解析の解析結果を得るまでに時間がかかるという課題を有している。
【００１３】
本発明は、上記の課題を鑑みてなされたものである。本発明の目的は、腸が発する音を解析するために複雑な演算処理を行うことなく、その音が蠕動音であるか否かを精度良く判定できる蠕動音検出装置および蠕動音検出方法を提供することにある。
【課題を解決するための手段】
【００１４】
本発明の一態様に係る蠕動音検出装置は、上記の課題を解決するために、腸が発する生体音を検知する生体音検知手段と、上記生体音の周波数スペクトルを算出する周波数スペクトル算出手段と、上記生体音の周波数スペクトルと、複数の蠕動音の標準周波数スペクトルのそれぞれとを個別にマッチングすることによって、複数のマッチング係数を算出するマッチング係数算出手段と、上記複数のマッチング係数を演算処理することによって、上記生体音が蠕動音であるか否かを判定する蠕動音判定手段と、を備えることを特徴とする。
【００１５】
上記のように構成された蠕動音検出装置によれば、生体音検知手段を介して腸が発する生体音を検知する。周波数スペクトル算出手段は、上記生体音から生体音の周波数スペクトルを算出する。マッチング係数算出手段は、上記生体音の周波数スペクトルと標準周波数スペクトルとをマッチングすることによってマッチング係数を算出する。この際、上記標準周波数スペクトルは複数であり、上記生体音の周波数スペクトルは、それぞれの上記標準周波数スペクトルとマッチングされる。したがって、上記生体音の周波数スペクトルと、それぞれの上記標準周波数スペクトルとから算出される上記マッチング係数は複数である。
【００１６】
このように複数のマッチング係数を算出することによって、上記生体音に含まれる多様な蠕動音の要素を検知した上で、上記生体音が蠕動音であるか否かを判定することができる。また、蠕動音算出手段、マッチング係数算出、および、マッチング係数判定手段は、いずれも複雑な演算処理を必要としない。
【００１７】
したがって、複雑な演算処理を行うことなく、生体音が蠕動音であるか否かを精度良く判定できる。言い換えると、複雑な演算処理を行うことなく、蠕動音を精度良く検出することができる。
【００１８】
また、本発明の一態様に係る蠕動音検出装置では、さらに、上記複数の蠕動音の標準周波数スペクトルは、それぞれが特定の発生モードに起因して発生する蠕動音の周波数スペクトルであることが好ましい。
【００１９】
上記複数の標準周波数スペクトルのそれぞれは、複数の発生モードの蠕動運動に起因する蠕動音の周波数スペクトルである。したがって、本発明の一態様に係る蠕動音検出装置は、生体音が複数の発生モードに起因する蠕動音の要素を含んでいても、生体音からそれぞれの蠕動音の要素を検知することができる。
【００２０】
また、本発明の一態様に係る蠕動音検出装置では、さらに、上記蠕動音判定手段は、上記複数のマッチング係数のうち最大のマッチング係数が、規定の閾値よりも大きいときに、上記生体音が蠕動音であると判定する構成としてもよい。
【００２１】
上記の構成によれば、上記蠕動音判定手段は、複数の上記マッチング係数のうち、最大のマッチング係数が規定の閾値よりも大きいときに、上記生体音が蠕動音であると判定する。したがって、本発明の一態様に係る蠕動音検出装置は、生体音が複数の発生モードに起因する蠕動音の要素を含んでいても、生体音が蠕動音であるか否かを精度良く判定できる。
【００２２】
また、本発明の一態様に係る蠕動音検出装置では、さらに、上記蠕動音判定手段は、上記複数のマッチング係数を全て合算した値が、規定の閾値よりも大きいときに、上記生体音が蠕動音であると判定する構成としてもよい。
【００２３】
上記の構成によれば、上記蠕動音判定手段は、複数の上記マッチング係数を合算した値が、規定の閾値よりも大きいときに、上記生体音が蠕動音であると判定する。したがって、本発明の一態様に係る蠕動音検出装置は、複数の発生モードの要素を少しずつ含む生体音であっても、その生体音が蠕動音であるか否かを精度良く判定できる。
【００２４】
また、本発明の一態様に係る蠕動音検出装置が備えている各手段としてコンピュータを動作させるためのプログラム、および、それらのプログラムを記録したコンピュータ読み取り可能な記録媒体についても本発明の範疇に含まれる。
【００２５】
本発明の一態様に係る蠕動音検出方法は、上記の課題を解決するために、腸が発する生体音を検知する生体音検知工程と、上記生体音の周波数スペクトルを算出する周波数スペクトル算出工程と、上記生体音の周波数スペクトルと、複数の蠕動音の標準周波数スペクトルのそれぞれとを個別にマッチングすることによって、複数のマッチング係数を算出するマッチング係数算出工程と、上記複数のマッチング係数を演算処理することによって、上記生体音が蠕動音であるか否かを判定する蠕動音判定工程と、を備えることを特徴とする。
【００２６】
上記のように構成された蠕動音検出方法によれば、上記蠕動音検出装置と同様の効果を奏する。
【発明の効果】
【００２７】
本発明の一態様に係る蠕動音検出装置および蠕動音検出方法によれば、腸が発する音を解析するために複雑な演算処理を行うことなく、その音が蠕動音であるか否かを判定する際の精度を向上させることができる。
【図面の簡単な説明】
【００２８】
【図１】本発明の一実施形態に係る蠕動音検出装置の構成を示すブロック図である。
【図２】本発明の一実施形態に係る蠕動音検出装置において、生体音が蠕動音であるか否かを判定する流れを示すフローチャートである。
【図３】本発明の一実施形態に係る蠕動音検出装置がマッチング係数を算出する際に用いる標準周波数スペクトルを示す図である。（ａ）は第１の標準周波数スペクトルを示し、（ｂ）は第２の標準周波数スペクトルを示し、（ｃ）は第３の標準周波数スペクトルを示す。
【図４】本発明の一実施例に係る蠕動音検出装置、および、比較例に係る蠕動音検出装置による蠕動音の検出率を示す図である。
【図５】Midozolam、BuprenorphineおよびPropofolを投与されている被験者の腸蠕動音発生頻度を、本発明の一実施形態に係る蠕動音検出装置を用いて観察した結果を示す図である。
【図６】PropofolおよびDEXを投与されている被験者の腸蠕動音発生頻度を、本発明の一実施形態に係る蠕動音検出装置を用いて観察した結果を示す図である。
【発明を実施するための形態】
【００２９】
本発明の一実施形態に係る蠕動音検出装置１０について、図１〜３を参照しながら説明する。図１は、蠕動音検出装置１０の構成を示すブロック図である。図２は、蠕動音検出装置１０において、腸が発する生体音が蠕動音であるか否か判定する流れを示すフローチャートである。図３は、蠕動音検出装置１０がマッチング係数を算出する際に用いる第１〜第３の標準周波数スペクトルを示す図である。
【００３０】
図１に示すように蠕動音検出装置１０は、音響センサ１１、周波数スペクトル算出部１２、マッチング係数算出部１３、記憶部１４、および、蠕動音判定部１５を備えている。
【００３１】
蠕動音検出装置１０は、腸が発する生体音を検知および解析することによって、当該腸が発する音が蠕動音であるか否かを判定する装置である。言い換えると、蠕動音検出装置１０は腸が発する生体音を検出する装置である。以下において「腸が発する生体音」を単に「生体音」とも表記する。
【００３２】
以下に、図１および図３を参照しながら蠕動音検出装置１０が備える各部について説明する。蠕動音検出装置１０が、生体音が蠕動音であるか否かを判定する際の流れは、図２に示すフローチャートを参照しながら後述する。
【００３３】
（音響センサ１１）
生体音検知手段である音響センサ１１は、検知する音声信号を電気信号に変換するものである。音響センサ１１としては、たとえば密着型音響マイクを用いることができる。音響センサ１１は、蠕動音検出装置１０を操作する人によって、被験者の生体音を検知できる位置に固定されてもよい。生体音を検知する音響センサ１１の数は限定されず、１つでもよいし複数であってもよい。音響センサ１１は、検知する生体音を電気信号として周波数スペクトル算出部１２に出力する。
【００３４】
音響センサ１１は、電気信号を増幅するための増幅器（図示せず）を備えていてもよい。なお、上記増幅器は、音響センサ１１に設けられる構成に限定されず、周波数スペクトル算出部１２に設けられていてもよい。
【００３５】
（周波数スペクトル算出部１２）
周波数スペクトル算出手段である周波数スペクトル算出部１２は、アナログ／デジタル変換部（Ａ／Ｄ変換部）を備えていることが好ましい。Ａ／Ｄ変換部は、音響センサ１１より電気信号を受け取り、デジタルデータである生体音データに変換する。
【００３６】
次に、周波数スペクトル算出部１２は、生体音データに対して所定の時間ごとに高速フーリエ変換（Fast Fourier Transformation, ＦＦＴ）処理を行う。当該所定の時間を、以下においてＦＦＴ処理間隔とも呼ぶ。ＦＦＴ処理間隔は０．３秒から１．０秒の範囲内であることが好ましい。更には、ＦＦＴ処理間隔は０．３秒から０．５秒の範囲内であることがより好ましい。
【００３７】
ＦＦＴ処理間隔を長い時間に設定すると、複数の腸蠕動音がＦＦＴ処理間隔の時間内に含まれる虞がある。言い換えると、複数の腸蠕動音を１回の腸蠕動音と誤って数える可能性が生じる。この腸蠕動音の検知漏れを防ぐために、ＦＦＴ処理間隔は１．０秒以下であることが好ましく、更には０．５秒以下であることがより好ましい。
【００３８】
一方、ＦＦＴ処理時間を短い時間に設定すればするほど、ＦＦＴ処理を実行する頻度が増加する。すなわち、ＦＦＴ処理の回数増加に伴い、周波数スペクトル算出部１２における計算負荷が増大する。その結果として、後述するリアルタイム処理が難しくなる虞がある。したがって、ＦＦＴ処理間隔は０．３秒以上であることが好ましい。
【００３９】
なお、このＦＦＴ処理間隔の下限値は周波数スペクトル算出部１２の演算能力に依存する。周波数スペクトル算出部１２の演算能力が十分に高ければ、ＦＦＴ処理間隔を０．３秒より短い時間に設定してもリアルタイム処理を実現することができる。すなわち、ＦＦＴ処理時間の下限値は０．３秒に限定さるものではない。
【００４０】
本実施形態では、ＦＦＴ処理間隔を０．３２秒に設定している。すなわち、時間的に連続する生体音データを、周波数スペクトル算出部１２は０．３２秒ごとの生体音データに区切る。その上で、区切られた各生体音データを順次ＦＦＴ処理する。
【００４１】
ＦＦＴ処理の方法には、周知の方法からいずれかの方法を選択して用いることができる。本実施形態において周波数スペクトル算出部１２は、ＭａｔｈＷｏｒｋｓ社製ＭＡＴＬＡＢ（登録商標）に組み込まれている関数spectrogramを用いてＦＦＴ処理を実行する。
【００４２】
周波数スペクトル算出部１２は、生体音データをＦＦＴ処理することにより、生体音データの周波数スペクトル（以下において、生体音スペクトルとも記載する）を算出する。より詳しくは、周波数スペクトル算出部１２は、０．３２秒間の生体音データをＦＦＴ処理することによって、１つの生体音スペクトルを算出する。
【００４３】
周波数スペクトル算出部１２は、所定の時間ごとに算出される生体音スペクトルを、順次、マッチング係数算出部１３に出力する。なお、以下においては、各部の説明を分かりやすくするために、所定の時間ごとに算出される生体音スペクトルのうち１つの生体音スペクトルに着目して説明する。
【００４４】
（標準周波数スペクトル）
後述するマッチング係数算出部１３は、周波数スペクトル算出部１２から生体音スペクトルを受け取るとともに、生体音スペクトルとマッチングするための標準周波数スペクトルを記憶部１４から読み出す。この際、マッチング係数算出部１３は、記憶部１４から複数の標準周波数スペクトルを読み出す。本実施形態では、あらかじめ医師が蠕動音であると認識した生体音（以下において標準蠕動音と呼ぶ）を複数抽出し、それぞれの標準蠕動音をＦＦＴ処理したものを標準周波数スペクトルとして記憶部１４に格納している。
【００４５】
腸が蠕動運動する際、発生する蠕動音の発生モードは複数ある。そして、それぞれの発生モードに応じて、腸は異なる蠕動音を発する。したがって、マッチングに用いる標準周波数スペクトルは、複数の異なる発生モードに起因する蠕動音をＦＦＴ処理したものであることが好ましい。
【００４６】
蠕動音は個人によって大きくばらつくことがあるため、複数の人から抽出した複数の標準蠕動音をＦＦＴ処理することによって複数の標準周波数スペクトルとしてもよい。また、複数の標準蠕動音をＦＦＴ処理し、共通のスペクトルを特に抽出したものを標準周波数スペクトルとしてもよい。さらに、腸活動が良好な時に発生する生体音をＦＦＴ処理し、統計処理することで標準周波数スペクトルを算出してもよい。なお、複数の標準周波数スペクトルは、異なる発生モードに起因する標準蠕動音から得られる標準周波数スペクトルと、複数の人から抽出した標準蠕動音から得られる標準周波数スペクトルとを組み合わせることにより構成されてもよい。
【００４７】
本実施形態では、第１の標準周波数スペクトル（図３の（ａ））、第２の標準周波数スペクトル（図３の（ｂ））、および、第３の標準周波数スペクトル（図３の（ｃ））を比較参照用の標準周波数スペクトルとして用いる。
【００４８】
なお、マッチングに用いる標準周波数スペクトルは図３の（ａ）〜（ｃ）に示す３つに限られない。すなわち、２つ、または、４つ以上の標準周波数スペクトルをマッチングに用いてもよい。
【００４９】
（マッチング係数算出部１３）
マッチング係数算出手段であるマッチング係数算出部１３は、生体音スペクトルと、それぞれの標準周波数スペクトルとをマッチングする（相関を見る、比較する、合致点を抽出する、とも言う）ことによって、複数のマッチング係数（相関係数）を算出する。
【００５０】
より具体的には、生体音スペクトルと第１の標準周波数スペクトルとをマッチングすることによって、第１のマッチング係数Ｃ1を算出する。同様にマッチング係数算出部１３は、生体音スペクトルと第２の標準周波数スペクトルとをマッチングすることによって第２のマッチング係数Ｃ2を算出し、生体音スペクトルと第３の標準周波数スペクトルとをマッチングすることによって第３のマッチング係数Ｃ3を算出する。
【００５１】
マッチング係数とは、生体音スペクトルと、標準周波数スペクトルとにおける類似の程度を表す係数である。別の言い方をすれば、生体音スペクトルが含んでいる標準周波数スペクトルの要素の大きさを表す係数である。したがって、マッチング係数が大きいほど生体音スペクトルと、標準周波数スペクトルとが類似している、すなわち、生体音スペクトルが標準周波数スペクトルの要素を多く含んでいることを意味する。
【００５２】
生体音スペクトルおよび標準周波数スペクトルとのマッチングには、周知の手法を用いることができる。本実施形態においてマッチング係数算出部１３は、ＭａｔｈＷｏｒｋｓ社製ＭＡＴＬＡＢ（登録商標）に組み込まれている関数corrcoefを用いてマッチング係数を算出する。
【００５３】
マッチング係数算出部１３は、０〜１の範囲に含まれる実数となるマッチング係数を算出するように構成されていてもよい。マッチング係数算出部１３がこのように構成されている場合、マッチング係数が１であることは、蠕動音スペクトルと標準スペクトルとが一致していることを意味する。
【００５４】
マッチング係数算出部１３は、マッチングにより算出した複数のマッチング係数Ｃ1、Ｃ2およびＣ3を蠕動音判定部１５に出力する。
【００５５】
腸が発する生体音が、複数の発生モードのうち少なくとも１つの発生モードに起因する蠕動音の要素を含んでいると、蠕動音検出装置１０はその蠕動音の要素をＣ1、Ｃ2およびＣ3の少なくとも１つとして検知することができる。言い換えると、生体音が複数の発生モードに起因する蠕動音の組み合わせからなる場合であっても、蠕動音検出装置１０はそれぞれの蠕動音の要素を別個に検知することができる。したがって、蠕動音の検知漏れを低減し、生体音が蠕動音であるか否かを精度良く判定することができる。
【００５６】
なお、複数の人から抽出した蠕動音をＦＦＴ処理することによって複数の標準周波数スペクトルとすると、蠕動音の個人差に起因する蠕動運動の検出漏れを低減することができる。
【００５７】
（蠕動音判定部１５）
蠕動音判定手段である蠕動音判定部１５は、マッチング係数算出部１３より受け取る複数のマッチング係数Ｃ1、Ｃ2およびＣ3を演算処理することによって、判定用マッチング係数Ｃmを算出する。具体的には、Ｃ1、Ｃ2およびＣ3の大小関係を比較演算し、Ｃ1、Ｃ2およびＣ3のうち最大のマッチング係数を判定用マッチング係数Ｃmとする。
【００５８】
さらに、蠕動音判定部１５は、判定用マッチング係数Ｃm、および、規定の閾値Ｃthの大小関係を比較する。その結果、Ｃm＞Ｃthであれば生体音が蠕動音であると判定し、Ｃm≦Ｃthであれば生体音が蠕動音でないと判定する。
【００５９】
閾値Ｃthは正の実数であれば任意に設定することができ、経験的には０．５≦Ｃth≦０．９の範囲に設定されていることが好ましい。たとえば、閾値Ｃthを０．８とすることができる。閾値Ｃthを大きな値に設定すると、蠕動音であると判定する基準が厳しくなる。一方、閾値Ｃthを小さな値に設定すると、蠕動音であると判定する基準が甘くなる。言い換えると、閾値Ｃthが大きすぎると腸蠕動音を腸蠕動音ではないと誤検出する可能性が高くなり、閾値Ｃthが小さすぎると腸蠕動音以外の生体音を腸音として誤検出する可能性が高くなる。
【００６０】
なお、腸蠕動音を的確に検出できる閾値Ｃthの値は、マッチングに使用する標準周波数スペクトル、より具体的には複数の標準周波数スペクトルの組み合わせに依存している。したがって、マッチングに使用する複数の標準周波数スペクトルの組み合わせに応じて、実験的に最適な閾値Ｃthの値を定めることが好ましい。
【００６１】
閾値Ｃthは、記憶部１４にあらかじめ格納されており、必要に応じて蠕動音判定部１５が読み出す構成としてもよい。蠕動音検出装置１０は、外部からの操作によって閾値Ｃthを任意に変更できるように構成されていてもよい。また、蠕動音判定部１５が、図１には図示しないＲＯＭ（Read Only Memory）を備えており、閾値Ｃthは上記ＲＯＭにあらかじめ格納されている構成としてもよい。
【００６２】
（リアルタイム処理）
以上の説明では、各部の働きを分かりやすくするために、１つの蠕動音スペクトルに着目して説明してきた。実際には、周波数スペクトル算出部１２の項に記載してあるとおり、周波数スペクトル算出部１２は所定の時間ごとに継続して蠕動音スペクトルを算出している。
【００６３】
さらに、マッチング係数算出部１３および蠕動音判定部１５において行う処理は、複雑な解析処理を伴うものではないため上記所定の時間と比較して短時間のうちに処理できる。言い換えると、マッチング係数算出部１３および蠕動音判定部１５の処理は、蠕動音であるか否かの判定に遅延を生じさせる処理ではない。
【００６４】
したがって、蠕動音検出装置１０は、蠕動音の検知、蠕動音スペクトルの算出、複数のマッチング係数の算出、および、判定用マッチング係数Ｃmと閾値Ｃthとの比較による蠕動音であるか否かの判定という一連の処理を、リアルタイムに処理することができる。
【００６５】
（蠕動音回数の検出）
蠕動音検出装置１０は、単位時間当たりに腸が発する蠕動音の回数を検出することもできる。具体的には、所定の単位時間内に蠕動音判定部１５が蠕動音であると判定した回数をカウントしておけばよい。
【００６６】
蠕動音検出装置１０は、生体音が蠕動音であるか否かを精度良く判定することができるので、単位時間における蠕動音のカウントに際に重複してカウントすることを避けることができる。また、蠕動音検出装置１０は、検出した蠕動音の単位時間当たりの発生回数を算出することによって、腸活動状態の変化を調べることもできる。ここでいう腸活動とは、摂取物によって引き起こされた腸活動と、摂取物によらずに発生する腸活動との両方の意味を含んでいる。
【００６７】
単位時間当たりに腸が発する蠕動音の回数を検知し続けることによって、腸の活性状態を観察することに蠕動音検出装置１０を利用できる。
【００６８】
（蠕動音判定部１５の変形例）
ここでは、蠕動音判定部１５の変形例である蠕動音判定部１５’について説明する。
【００６９】
蠕動音判定部１５’は蠕動音判定部１５と比較して、判定用マッチング係数Ｃmの算出方法が異なっている。具体的には、蠕動音判定部１５’は、第１〜第３のマッチング係数Ｃ1、Ｃ2およびＣ3を全て合算した値を判定用マッチング係数Ｃmとして算出する。
【００７０】
蠕動音であるか否かを判定する演算処理に関しては、蠕動運動判定部１５と同様である。
【００７１】
腸の蠕動運動には、複数の発生モードが存在することは先に述べた。蠕動運動は、単一の発生モードのみによって発生するとは限らず、複数の発生モードの要素を少しずつ含む蠕動運動が発生することもある。蠕動音検出装置１０が蠕動音判定部１５’を備えることによって、生体音が複数の発生モードの要素を少しずつ含む場合であっても蠕動音として検知することができる。したがって、生体音が蠕動音であるか否かを判定する際に、その精度を向上させることができる。
【００７２】
蠕動音検出装置１０が蠕動音判定部１５’を備える際に用いる規定の閾値Ｃth’は、正の実数であれば任意に設定することができ、経験的には１≦Ｃth’≦１．８の範囲に設定されていることが好ましい。たとえば閾値Ｃth’を１．２５とすることができる。
【００７３】
閾値Ｃth’が大きすぎると腸蠕動音を腸蠕動音ではないと誤検出する可能性が高くなり、閾値Ｃth’が小さすぎると腸蠕動音以外の音を生体音として誤検出する可能性が高くなることは、蠕動音判定部１５における閾値Ｃthと同様である。また、マッチングに使用する複数の標準周波数スペクトルの組み合わせに応じて、実験的に最適な閾値Ｃth’の値を定めることが好ましいことも閾値Ｃthと同様である。
【００７４】
（蠕動音検出の流れ）
蠕動音検出装置１０において、生体音が蠕動音であるか否かを判定する方法を、図２に示すフローチャートを参照しながら説明する。
【００７５】
（Ｓ１０１）
蠕動音検出装置１０は、音響センサ１１を介して生体音を検知し、電気信号として周波数スペクトル算出部１２に出力する（生体音検知工程）。
【００７６】
（Ｓ１０２）
周波数スペクトル算出部１２は、蠕動音データを所定の時間ごとにＦＦＴ処理することによって、生体音の周波数スペクトル（生体音スペクトル）を算出する（周波数スペクトル算出工程）。
【００７７】
（Ｓ１０３）
記憶部１４には、複数の発生モードによって発生する標準蠕動音の周波数スペクトルである標準周波数スペクトルがあらかじめ格納されている。本実施形態においては、記憶部１４には３つの標準周波数スペクトルが格納されている。それぞれの発生モードに対応する標準周波数スペクトルを、第１の標準周波数スペクトル、第２の標準周波数スペクトル、および、第３の標準周波数スペクトルと呼ぶ。
【００７８】
（Ｓ１０４）
マッチング係数算出部１３は、周波数スペクトル算出部１２から生体音スペクトルを受け取り、記憶部１４から第１〜第３の標準周波数スペクトルを読み出す。マッチング係数算出部１３は、生体音スペクトルと第１の標準周波数スペクトルとをマッチングすることによって、第１のマッチング係数Ｃ1を算出する。同様にマッチング係数算出部１３は、生体音スペクトルと第２の標準周波数スペクトルとをマッチングすることによって第２のマッチング係数Ｃ2を算出し、生体音スペクトルと第３の標準周波数スペクトルとを比較演算することによって第３のマッチング係数Ｃ3を算出する（マッチング係数算出工程）。
【００７９】
（Ｓ１０５）
蠕動音判定部１５は、マッチング係数算出部１３より第１〜第３のマッチング係数Ｃ1、Ｃ2およびＣ3を受け取る。蠕動音判定部１５は、Ｃ1、Ｃ2およびＣ3を比較演算し、Ｃ1、Ｃ2およびＣ3のうち最大のマッチング係数を判定用マッチング係数Ｃmとして算出する。
【００８０】
さらに、蠕動音判定部１５は、判定用マッチング係数Ｃm、および、規定の閾値Ｃthとの大小関係を比較する（蠕動音判定工程）。
【００８１】
（Ｓ１０６）
判定用マッチング係数Ｃmが閾値Ｃthより大きい（Ｓ１０５がＹＥＳである）と、蠕動音判定部１５は、生体音が蠕動音であると判定する（蠕動音判定工程）。
【００８２】
（Ｓ１０７）
一方、判定用マッチング係数Ｃmが閾値Ｃth以下（Ｓ１０５がＮＯである）だと、蠕動音判定部１５は、生体音が蠕動音でないと判定する（マッチング係数判定工程）。
【００８３】
上述の各工程により生体音が蠕動音であるか否かを判定することによって、本発明の一実施形態に係る蠕動音検出方法は、複雑な演算処理を行うことなく判定精度を向上する効果を奏する。
【００８４】
〔実施例１〕
蠕動音検出装置１０を用いて、生体音が蠕動音であるか否かを判定した。なお、実施例１において用いる蠕動音検出装置１０は、蠕動音判定部１５を備えている。すなわち、Ｃ1、Ｃ2およびＣ3のうち最大のマッチング係数を判定用マッチング係数Ｃmとした。
【００８５】
腸の活性度を判定する際に用いた条件は以下の通りである。
・音響センサ１１として密着型音響マイクを用い、当該密着型音響マイクを被験者の腹部に密着固定した。
・周波数スペクトル算出部１２が生体音データをＦＦＴ処理する所定の時間間隔を０．３２秒とした。
・マッチング係数算出部１３にて用いる標準周波数スペクトルとして、第１の標準周波数スペクトル（図３の（ａ））、第２の標準周波数スペクトル（図３の（ｂ））、および、第３の標準周波数スペクトル（図３の（ｃ））を用いた。
・蠕動音判定部１５にて用いる規定の閾値Ｃthを０．８とした。
・３０分間、連続して蠕動音の検出を行った。
【００８６】
このような条件のもと、蠕動音検出装置１０を用いて生体音が蠕動音であるか否かを判定すると同時に、医師による判定を行った。その結果、医師が「蠕動音である」と判断した生体音の数に対して、蠕動音検出装置１０が「蠕動音である」と判定した数の割合は９８％であった。以下において、当該割合のことを検出率と呼ぶ。
【００８７】
実施例１と、後述する実施例２および比較例とにおける検出率をまとめた表を図４に示す。
【００８８】
〔実施例２〕
蠕動音検出装置１０を用いて、生体音が蠕動音であるか否かを判定した。なお、実施例１において用いる蠕動音検出装置１０は、蠕動音判定部１５’を備えている。すなわち、Ｃ1、Ｃ2およびＣ3を全て合算した値を判定用マッチング係数Ｃmとした。
【００８９】
腸の活性度を判定する際に用いた条件は、規定の閾値Ｃth’の値を除いて実施例１と同様である。実施例２において、閾値Ｃth’は１．２５とした。
【００９０】
その結果、実施例２において得られた検出率は９９％であった。
【００９１】
〔比較例〕
実施例１および２に対する比較例として、１つの標準周波数スペクトルを用いて蠕動音を検出した。実施例１および２において第１〜第３の標準周波数スペクトルを用いたのに対して、比較例では第１の標準周波数スペクトル（図３の（ａ））のみを用いた。それ以外の条件については、実施例１と同様である。
【００９２】
その結果、比較例において得られた検出率は８２％であった。
【００９３】
以上のように、１つの標準周波数スペクトルを用いて生体音が蠕動音であるか否かを判定した場合と比較して、第１〜第３の標準周波数スペクトルを用いた実施例１および２では大きく検出率が向上した。すなわち、複数の標準周波数スペクトルを用いて生体音が蠕動音であるか否かを判定することによって判定の精度が向上した。
【００９４】
〔実施例３〕
蠕動音検出装置１０を用いて、被験者に接種される薬剤が腸の消化活性に与える影響を調べた。具体的には、１分間のうちに蠕動音検出装置１０が「蠕動音である」と判定した回数を腸蠕動音発生回数として算出した。
【００９５】
本実施例において、蠕動音検出装置１０は蠕動音判定部１５’を備えている。すなわち、Ｃ1、Ｃ2およびＣ3を全て合算した値を判定用マッチング係数Ｃmとした。腸の活性度を判定する際に用いた条件は以下の通りである。
・音響センサ１１として４つの密着型音響マイクを用い、それぞれの密着型音響マイクを被験者の腹部に密着固定した。
・周波数スペクトル算出部１２が生体音データをＦＦＴ処理する所定の時間間隔を０．３２秒とした。
・マッチング係数算出部１３にて用いる標準周波数スペクトルとして、第１の標準周波数スペクトル（図３の（ａ））、第２の標準周波数スペクトル（図３の（ｂ））、および、第３の標準周波数スペクトル（図３の（ｃ））を用いた。
・蠕動音判定部１５にて用いる規定の閾値Ｃth’を１．５とした。
・２４時間連続して蠕動音の検出を行い、１分間当たりの腸蠕動音発生回数を算出した。
【００９６】
まず、被験者に鎮静・鎮痛剤としてMidazolamを４ｍｇ／ｈ、Buprenorphineを０．００８ｍｇ／ｈ、および、Propofolを１０ｍｇ／ｈ投与しながら、蠕動音検出装置１０を用いて腸蠕動音を検知した。得られた結果より、１分間当たりの腸蠕動音発生回数を算出した結果を図５に示す。
【００９７】
次の日に、同じ被験者に鎮静剤としてPropofolを３０ｍｇ／ｈ、および、DEXを０．６μｇ／ｋｇ／ｈ投与しながら、蠕動音検出装置１０を用いて腸蠕動音を検知した。検知した腸蠕動音の回数より、１分間当たりの腸蠕動音発生回数を算出した結果を図６に示す。
【００９８】
図５および図６を比較すると、図６の方が腸の蠕動音の発生回数が多いことがわかった。すなわち、「Midazolam+Buprenorphine+Propofol」よりも「Propofol+DEX」の方が腸の活性状態を高める作用が強いことを示す結果が得られた。
【００９９】
本実施例に加えて、蠕動音検出装置１０は、薬剤または栄養剤などの経口摂取に伴う腸活動の変化、血糖値と腸活動の相関、および、サイトカインと腸活動の相関の計測などに利用できる。
【０１００】
（付記事項）
上述した蠕動音検出装置１０の各ブロックは、集積回路（ＩＣチップ）上に形成された論理回路によってハードウェア的に実現してもよいし、ＣＰＵ（Central Processing Unit）を用いてソフトウェア的に実現してもよい。
【０１０１】
後者の場合、上記装置は、各機能を実現するプログラムの命令を実行するＣＰＵ、上記プログラムを格納したＲＯＭ（Read Only Memory）、上記プログラムを展開するＲＡＭ（Random Access Memory）、上記プログラム及び各種データを格納するメモリ等の記憶装置（記録媒体）などを備えている。そして、本発明の目的は、上述した機能を実現するソフトウェアである上記各装置の制御プログラムのプログラムコード（実行形式プログラム、中間コードプログラム、ソースプログラム）をコンピュータ読み取り可能に記録した記録媒体を、上記装置に供給し、そのコンピュータ（またはＣＰＵやＭＰＵ）が記録媒体に記録されているプログラムコードを読み出し実行することによっても、達成可能である。
【０１０２】
上記記録媒体としては、例えば、磁気テープやカセットテープ等のテープ類、フロッピー（登録商標）ディスク／ハードディスク等の磁気ディスクやＣＤ−ＲＯＭ／ＭＯ／ＭＤ／ＤＶＤ／ＣＤ−Ｒ等の光ディスクを含むディスク類、ＩＣカード（メモリカードを含む）／光カード等のカード類、マスクＲＯＭ／ＥＰＲＯＭ／ＥＥＰＲＯＭ／フラッシュＲＯＭ等の半導体メモリ類、あるいはＰＬＤ（Programmable logic device）やＦＰＧＡ（Field Programmable Gate Array）等の論理回路類などを用いることができる。
【０１０３】
また、上記各装置を通信ネットワークと接続可能に構成し、上記プログラムコードを通信ネットワークを介して供給してもよい。この通信ネットワークは、プログラムコードを伝送可能であればよく、特に限定されない。例えば、インターネット、イントラネット、エキストラネット、ＬＡＮ、ＩＳＤＮ、ＶＡＮ、ＣＡＴＶ通信網、仮想専用網（Virtual Private Network）、電話回線網、移動体通信網、衛星通信網等が利用可能である。また、この通信ネットワークを構成する伝送媒体も、プログラムコードを伝送可能な媒体であればよく、特定の構成または種類のものに限定されない。例えば、ＩＥＥＥ１３９４、ＵＳＢ、電力線搬送、ケーブルＴＶ回線、電話線、ＡＤＳＬ（Asymmetric Digital Subscriber Line）回線等の有線でも、ＩｒＤＡやリモコンのような赤外線、Ｂｌｕｅｔｏｏｔｈ（登録商標）、ＩＥＥＥ８０２．１１無線、ＨＤＲ（High Data Rate）、ＮＦＣ（Near Field Communication）、ＤＬＮＡ（Digital Living Network Alliance）、携帯電話網、衛星回線、地上波デジタル網等の無線でも利用可能である。
【０１０４】
本発明は上述した各実施形態に限定されるものではなく、請求項に示した範囲で種々の変更が可能であり、異なる実施形態にそれぞれ開示された技術的手段を適宜組み合わせて得られる実施形態についても本発明の技術的範囲に含まれる。
【産業上の利用可能性】
【０１０５】
本発明は、腸が発する音が蠕動音であるか否かを判定する蠕動音検出装置および蠕動音検出方法として利用できる。
【符号の説明】
【０１０６】
１０  蠕動音検出装置
１１  音響センサ（生体音検知手段）
１２  周波数スペクトル算出部（周波数スペクトル算出手段）
１３  マッチング係数算出部（マッチング係数算出手段）
１４  記憶部
１５  蠕動音判定部（蠕動音判定手段）
【図１】
0005877511.tif 000002
【図２】
0005877511.tif 000003
【図３】
0005877511.tif 000004
【図４】
0005877511.tif 000005
【図５】
0005877511.tif 000006
【図６】
0005877511.tif 000007