特許 6298339 呼吸音解析装置及び呼吸音解析方法、並びにコンピュータプログラム及び記録媒体

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6298339

(24)【登録日】2018年3月2日

(45)【発行日】2018年3月20日

(54)【発明の名称】呼吸音解析装置及び呼吸音解析方法、並びにコンピュータプログラム及び記録媒体

(51)【国際特許分類】

   A61B 7/04 20060101AFI20180312BHJP

   A61B 5/08 20060101ALI20180312BHJP

【ＦＩ】

   A61B7/04 LZDM

   A61B7/04 U

   A61B5/08

【請求項の数】5

【全頁数】12

(21)【出願番号】特願2014-68856(P2014-68856)

(22)【出願日】2014年3月28日

(65)【公開番号】特開2015-188637(P2015-188637A)

(43)【公開日】2015年11月2日

【審査請求日】2017年2月7日

(73)【特許権者】

【識別番号】000005016

【氏名又は名称】パイオニア株式会社

(74)【代理人】

【識別番号】100104765

【弁理士】

【氏名又は名称】江上 達夫

(74)【代理人】

【識別番号】100107331

【弁理士】

【氏名又は名称】中村 聡延

(72)【発明者】

【氏名】大久保 英幸

(72)【発明者】

【氏名】亀谷 隆真

(72)【発明者】

【氏名】長谷部 剛

(72)【発明者】

【氏名】石戸谷 耕一

(72)【発明者】

【氏名】三浦 友博

【審査官】 樋熊 政一

(56)【参考文献】

【文献】 米国特許第６１３９５０５（ＵＳ，Ａ）

【文献】 米国特許出願公開第２００２／０１８３６４２（ＵＳ，Ａ１）

【文献】 特開２０１４−０５０６７２（ＪＰ，Ａ）

【文献】 特開２０１０−１８３９８０（ＪＰ，Ａ）

【文献】 特表２００１−５０５０８５（ＪＰ，Ａ）

【文献】 特開２０１３−１２３４９４（ＪＰ，Ａ）

【文献】 特表２００５−５０５３２９（ＪＰ，Ａ）

【文献】 酒井 智弥 Tomoya SAKAI，スパース表現による異常肺音の抽出 Extraction of Abnormal Lung Sounds via Sparse Representation，電子情報通信学会技術研究報告 Ｖｏｌ．１１１ Ｎｏ．２８ IEICE Technical Report，日本，社団法人電子情報通信学会 The Institute of Electro，第111巻

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ａ６１Ｂ ７／００−７／０４

Ａ６１Ｂ ５／０６−５／２２

(57)【特許請求の範囲】

【請求項1】

呼吸音から連続性ラ音を含む呼吸音成分を取得する成分取得手段と、
前記呼吸音成分に含まれる所定の特徴に対応する周波数を取得する周波数取得手段と、
前記所定の特徴に対応する周波数に基づいて、前記呼吸音成分に含まれる笛声音成分と類鼾音成分との割合を出力する出力手段と、
を備え、
前記出力手段は、前記周波数が高いほど、前記笛声音成分が含まれる割合が増加し、前記類鼾音成分が含まれる割合が減少するように前記割合を出力することを特徴とする呼吸音解析装置。

【請求項2】

前記所定の特徴は、極大値であることを特徴とする請求項１に記載の呼吸音解析装置。

【請求項3】

呼吸音から連続性ラ音を含む呼吸音成分を取得する成分取得工程と、
前記呼吸音成分に含まれる所定の特徴に対応する周波数を取得する周波数取得工程と、
前記所定の特徴に対応する周波数に基づいて、前記呼吸音成分に含まれる笛声音成分と類鼾音成分との割合を出力する出力工程と、
を備え、
前記出力工程では、前記周波数が高いほど、前記笛声音成分が含まれる割合が増加し、前記類鼾音成分が含まれる割合が減少するように前記割合を出力することを特徴とする呼吸音解析方法。

【請求項4】

呼吸音から連続性ラ音を含む呼吸音成分を取得する成分取得工程と、
前記呼吸音成分に含まれる所定の特徴に対応する周波数を取得する周波数取得工程と、
前記所定の特徴に対応する周波数に基づいて、前記呼吸音成分に含まれる笛声音成分と類鼾音成分との割合を出力する出力工程と、
をコンピュータに実行させ、
前記出力工程では、前記周波数が高いほど、前記笛声音成分が含まれる割合が増加し、前記類鼾音成分が含まれる割合が減少するように前記割合を出力することを特徴とするコンピュータプログラム。

【請求項5】

請求項４に記載のコンピュータプログラムが記録されていることを特徴とする記録媒体。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、例えば連続性ラ音を含む呼吸音を解析する呼吸音解析装置及び呼吸音解析方法、並びにコンピュータプログラム及び記録媒体の技術分野に関する。

【背景技術】

【0002】

この種の装置として、電子聴診器等によって検出される生体の呼吸音について、含まれている複数の音種（例えば、正常呼吸音と異常呼吸音）を夫々判別するものが知られている。例えば特許文献１では、スペクトル上の局所分散値に基づいて、正常呼吸音と連続性ラ音との判別を行うという手法が提案されている。

【先行技術文献】

【特許文献】

【0003】

【特許文献1】特開２００４−３５７７５８号公報

【発明の概要】

【発明が解決しようとする課題】

【0004】

しかしながら、上述した特許文献１に記載されているような技術では、連続性ラ音を笛声音と類鼾音とに分類することができない。即ち、呼吸音を正常呼吸音と異常呼吸音である連続性ラ音とに分離できたとしても、連続性ラ音について更なる分離を行えないという技術的問題点がある。

【0005】

本発明が解決しようとする課題には、上記のようなものが一例として挙げられる。本発明は、呼吸音に含まれる連続性ラ音を好適に解析可能な呼吸音解析装置及び呼吸音解析方法、並びにコンピュータプログラム及び記録媒体を提供することを課題とする。

【課題を解決するための手段】

【0006】

上記課題を解決するための呼吸音解析装置は、呼吸音から連続性ラ音を含む呼吸音成分を取得する成分取得手段と、前記呼吸音成分に含まれる所定の特徴に対応する周波数を取得する周波数取得手段と、前記所定の特徴に対応する周波数と、前記所定の特徴に対応する周波数に応じて変動する閾値との関係に基づいて、前記呼吸音成分に含まれる笛声音成分と類鼾音成分とを判定する判定手段とを備える。

【0007】

上記課題を解決するための呼吸音解析方法は、呼吸音から連続性ラ音を含む呼吸音成分を取得する成分取得工程と、前記呼吸音成分に含まれる所定の特徴に対応する周波数を取得する周波数取得工程と、前記所定の特徴に対応する周波数と、前記所定の特徴に対応する周波数に応じて変動する閾値との関係に基づいて、前記呼吸音成分に含まれる笛声音成分と類鼾音成分とを判定する判定工程とを備える。

【0008】

上記課題を解決するためのコンピュータプログラムは、呼吸音から連続性ラ音を含む呼吸音成分を取得する成分取得工程と、前記呼吸音成分に含まれる所定の特徴に対応する周波数を取得する周波数取得工程と、前記所定の特徴に対応する周波数と、前記所定の特徴に対応する周波数に応じて変動する閾値との関係に基づいて、前記呼吸音成分に含まれる笛声音成分と類鼾音成分とを判定する判定工程とをコンピュータに実行させる。

【0009】

上記課題を解決するための記録媒体は、上述したコンピュータプログラムが記録されている。

【図面の簡単な説明】

【0010】

【図1】本実施例に係る呼吸音解析装置の全体構成を示すブロック図である。

【図2】本実施例に係る呼吸音解析装置の動作を示すフローチャートである。

【図3】比較例に係る連続性ラ音判定における問題点を示す概念図である。

【図4】本実施例に係る閾値の変動を示すグラフである。

【図5】本実施例に係る閾値の具体的な値を示す表である。

【図6】本実施例に係る連続性ラ音判定結果の表示例を示すグラフ（その１）である。

【図7】本実施例に係る連続性ラ音判定結果の表示例を示すグラフ（その２）である。

【図8】本実施例及び比較例に係る呼吸音解析装置により笛声音を解析する例を示す概念図である。

【図9】本実施例及び比較例に係る呼吸音解析装置により類鼾音を解析する例を示す概念図である。

【発明を実施するための形態】

【0011】

＜１＞
本実施形態に係る呼吸音解析装置は、呼吸音から連続性ラ音を含む呼吸音成分を取得する成分取得手段と、前記呼吸音成分に含まれる所定の特徴に対応する周波数を取得する周波数取得手段と、前記所定の特徴に対応する周波数と、前記所定の特徴に対応する周波数に応じて変動する閾値との関係に基づいて、前記呼吸音成分に含まれる笛声音成分と類鼾音成分とを判定する判定手段とを備える。

【0012】

本実施形態に係る呼吸音解析装置によれば、その動作時には、先ず呼吸音から連続性ラ音を含む呼吸音成分が取得される。具体的には、正常呼吸音や異常呼吸音等の様々な音を含む呼吸音から、周知の技術等を利用して、連続性ラ音を含む呼吸音成分が抽出される。ただし、ここで取得される呼吸音成分は、連続性ラ音のみを含む呼吸音成分に限定されるものではない。即ち、取得される呼吸音成分は、連続性ラ音以外の音を含んでいても構わない。

【0013】

呼吸音成分が取得されると、呼吸音成分に含まれる所定の特徴に対応する周波数が取得される。なお、ここでの「所定の特徴」とは、呼吸音成分に含まれる音種に応じて特定の周波数に発生する特徴を意味しており、例えば周波数解析された信号に現れるピーク等である。

【0014】

所定の特徴に対応する周波数が取得されると、取得した周波数に基づいて、呼吸音成分に含まれる笛声音成分と類鼾音成分とが判定される。即ち、連続性ラ音を含む呼吸音成分に、連続性ラ音である笛声音成分及び類鼾音成分が含まれているか否かが判定される。なお、ここでの判定は、呼吸音成分中の笛声音成分及び類鼾音成分の存在に関する判定であれば特に限定されるものではない。例えば、呼吸音成分に、笛声音成分及び類鼾音成分のいずれが含まれているか、又は両方が含まれているのかを判定するものであってもよい。或いは、笛声音成分及び類鼾音成分がどのような割合で含まれているのかを判定するものであってもよい。また、笛声音成分及び類鼾音成分が含まれる可能性を判定するものであっても構わない。

【0015】

上述した判定は、所定の特徴に対応する周波数と、所定の特徴に対応する周波数に応じて変動する閾値との関係に基づいて行われる。例えば、取得された周波数と閾値とが互いに比較され、その大小関係に基づいて、笛声音成分及び類鼾音成分の判定が行われる。より具体的には、例えば取得された周波数が閾値より高い場合には笛声音成分であると判定され、低い場合には類鼾音成分と判定される。

【0016】

ここで特に、判定に用いられる閾値は、取得された所定の特徴に対応する周波数に応じて変動するものである。例えば閾値は、取得された周波数が高いほど、類鼾音よりも笛声音と判定する割合が大きくなるように変動する。このように、周波数に応じて閾値を変動させれば、固定値である閾値を用いる場合と比べて、より適切に笛声音及び類鼾音を判別することができる。なお、閾値をどのように変動させるかは、実験結果等に基づいて予め設定しておけばよい。また、閾値の変動特性を、被測定対象である生体の性別、年齢、身長、体重等に応じて決定するようにしてもよい。

【0017】

以上説明したように、本実施形態に係る呼吸音解析装置によれば、連続性ラ音である笛声音及び類鼾音を好適に判別することが可能である。

【0018】

＜２＞
本実施形態に係る呼吸音解析装置の一態様では、前記所定の特徴は、極大値である。

【0019】

この態様によれば、例えば呼吸音を示す信号に対して、高速フーリエ変換（ＦＦＴ：Fast Fourier Transform）等による周波数解析が実行され、解析結果の極大値（即ち、ピーク）に対応する周波数に関する情報が取得される。なお、周波数に関する情報は、極大値の位置に対応するものとして取得されるが、極大値の位置と完全に一致する周波数でなくとも、極大値の近傍位置に対応する周波数に関する情報として取得されてもよい。

【0020】

上述したように、呼吸音における所定の特徴として極大値を利用することで、より容易且つ的確に周波数を取得できる。

【0021】

＜３＞
本実施形態に係る呼吸音解析装置の他の態様では、前記判定手段は、前記呼吸音成分に含まれる前記笛声音成分と前記類鼾音成分との割合を判定する。

【0022】

この態様によれば、呼吸音成分に、笛声音成分と類鼾音成分とがどのような割合で含まれているかが判定される。よって、例えば呼吸音に基づく健康状態の診断等を適切に行うことが可能となる。なお、笛声音成分及び類鼾音成分以外の成分についても割合を判定するようにしてもよい。

【0023】

＜４＞
上述した笛声音成分と類鼾音成分との割合を判定する態様では、前記閾値は、前記所定の特徴に対応する周波数が高いほど、前記笛声音成分が含まれる割合が増加し、前記類鼾音成分が含まれる割合が減少するように判定させるものとして変動してもよい。

【0024】

この場合、類鼾音成分と比べて周波数が高い傾向にある笛声音成分（言い換えれば、笛声音成分と比べて周波数が低い傾向にある類鼾音成分）を精度よく判別することが可能である。なお、判定結果の割合の増減については、直線的な変動であってもよいし、所定の関数に応じた比較的複雑な変動であってもよい。或いは、段階的な変動であっても構わない。

【0025】

＜５＞
本実施形態に係る呼吸音解析方法は、呼吸音から連続性ラ音を含む呼吸音成分を取得する成分取得工程と、前記呼吸音成分に含まれる所定の特徴に対応する周波数を取得する周波数取得工程と、前記所定の特徴に対応する周波数と、前記所定の特徴に対応する周波数に応じて変動する閾値との関係に基づいて、前記呼吸音成分に含まれる笛声音成分と類鼾音成分とを判定する判定工程とを備える。

【0026】

本実施形態に係る呼吸音解析方法によれば、上述した本実施形態に係る呼吸音解析装置と同様に、連続性ラ音である笛声音及び類鼾音を好適に判別することができる。

【0027】

なお、本実施形態に係る呼吸音解析方法においても、上述した本実施形態に係る呼吸音解析装置における各種態様と同様の各種態様を採ることが可能である。

【0028】

＜６＞
本実施形態に係るコンピュータプログラムは、呼吸音から連続性ラ音を含む呼吸音成分を取得する成分取得工程と、前記呼吸音成分に含まれる所定の特徴に対応する周波数を取得する周波数取得工程と、前記所定の特徴に対応する周波数と、前記所定の特徴に対応する周波数に応じて変動する閾値との関係に基づいて、前記呼吸音成分に含まれる笛声音成分と類鼾音成分とを判定する判定工程とをコンピュータに実行させる。

【0029】

本実施形態に係るコンピュータプログラムによれば、上述した本実施形態に係る呼吸音解析方法と同様の処理をコンピュータに実行させることができるため、連続性ラ音である笛声音及び類鼾音を好適に判別することができる。

【0030】

なお、本実施形態に係るコンピュータプログラムにおいても、上述した本実施形態に係る呼吸音解析装置における各種態様と同様の各種態様を採ることが可能である。

【0031】

＜７＞
本実施形態に係る記録媒体は、上述したコンピュータプログラムが記録されている。

【0032】

本実施形態に係る記録媒体によれば、上述したコンピュータプログラムをコンピュータにより実行させることにより、連続性ラ音である笛声音及び類鼾音を好適に判別することができる。

【0033】

本実施形態に係る呼吸音解析装置及び呼吸音解析方法、並びにコンピュータプログラム及び記録媒体の作用及び他の利得については、以下に示す実施例において、より詳細に説明する。

【実施例】

【0034】

以下では、図面を参照して呼吸音解析装置及び呼吸音解析方法、並びにコンピュータプログラム及び記録媒体の実施例について詳細に説明する。

【0035】

＜全体構成＞
先ず、本実施例に係る呼吸音解析装置の全体構成について、図１を参照して説明する。ここに図１は、本実施例に係る呼吸音解析装置の全体構成を示すブロック図である。

【0036】

図１において、本実施例に係る呼吸音解析装置は、主な構成要素として、生体音センサ１１０と、信号記憶部１２０と、信号処理部１２５と、音声出力部１３０と、表示部１４０と、処理部２００とを備えて構成されている。

【0037】

生体音センサ１１０は、生体の呼吸音を検出可能に構成されたセンサである。生体音センサ１１０は、例えばＥＣＭ（Electret Condenser Microphone）やピエゾを利用したマイク、振動センサ等で構成されている。

【0038】

信号記憶部１２０は、例えばＲＡＭ（Random Access Memory）等のバッファとして構成されており、生体音センサ１１０で検出された呼吸音を示す信号（以下、適宜「呼吸音信号」と称する）を一時的に記憶する。信号記憶部１２０は、記憶した信号を、音声出力部１３０及び処理部２００に夫々出力可能に構成されている。

【0039】

信号処理部１２５は、生体音センサ１１０で取得した音を加工して音声出力部１３０に出力する。信号処理部１２５は、例えばイコライザーやフィルターとして機能し、取得した音を人が聴き易い状態に加工する。

【0040】

音声出力部１３０は、例えばスピーカやヘッドホンとして構成されており、生体音センサ１１０で検出され、信号処理部１２５で加工された呼吸音を出力する。

【0041】

表示部１４０は、例えば液晶モニタ等のディスプレイとして構成されており、処理部２００から出力される画像データを表示する。

【0042】

処理部２００は、複数の演算回路やメモリ等を含んで構成されている。処理部２００は、周波数解析部２１０、連続性ラ音検出部２２０、ピーク周波数検出部２３０、笛声音／類鼾音判定部２４０、及び画像生成部２５０を備えている。

【0043】

処理部２００の各部の動作については後に詳述する。

【0044】

＜動作説明＞
次に、本実施例に係る呼吸音解析装置の動作について、図２を参照して説明する。ここに図２は、本実施例に係る呼吸音解析装置の動作を示すフローチャートである。

【0045】

図２において、本実施例に係る呼吸音解析装置の動作時には、先ず生体音センサ１１０において呼吸音が検出され、処理部２００による呼吸音信号の取得が行われる（ステップＳ１０１）。

【0046】

呼吸音信号が取得されると、周波数解析部２１０において周波数解析（例えば、高速フーリエ変換）が実行される（ステップＳ１０２）。周波数解析が実行されると、その結果を用いて、連続性ラ音検出部２２０において、呼吸音信号に含まれる連続性ラ音を含む成分が検出される（ステップＳ１０３）。なお、連続性ラ音を検出する処理は、連続性ラ音以外の呼吸音（例えば、正常呼吸音や、他の異常呼吸音等）を検出する処理と並行して行われてもよい。ここまでの各処理については、周知の技術を利用して行うことが可能であるため、詳細な説明は省略する。

【0047】

連続性ラ音が検出されると（ステップＳ１０４：ＹＥＳ）、周波数ピーク検出部２３０において、連続性ラ音に対応する成分のピーク（極大値）の検出が実行され、ピーク位置に対応する周波数がピーク周波数として検出される（ステップＳ１０５）。ピーク検出する場合には、例えば所定の時間間隔（例えば、ＦＦＴをかける時間間隔等）で、周波数特性の領域において最大値をとる周波数を求めればよい。

【0048】

本実施例に係る呼吸音解析装置は、ここで検出されたピーク周波数に基づいて、連続性ラ音が笛声音であるか、又は類鼾音であるかを判定する。具体的には、ピーク周波数と、所定の閾値との大小を比較して、ピーク周波数が所定の周波数より高い場合には笛声音であると判定し、ピーク周波数が所定の周波数より低い場合には類鼾音であると判定する。ただし、本実施例では特に、笛声音及び類鼾音を判定するための閾値として、ピーク周波数に応じて変動する値を用いる。このため本実施例では、ピーク周波数が検出されると、笛声音及び類鼾音の判定に先立って、ピーク周波数に応じた閾値の決定が実行される（ステップＳ１０６）。

【0049】

ここで、上述した閾値の決定方法について、図３から図５を参照して具体的に説明する。ここに図３は、比較例に係る連続性ラ音判定における問題点を示す概念図である。また図４は、本実施例に係る閾値の変動を示すグラフであり、図５は、本実施例に係る閾値の具体的な値を示す表である。なお、図３に示されるスペクトルは、正常呼吸音に加えて笛声音を顕著に含む呼吸音のスペクトルである。

【0050】

図３において、笛声音は高音性連続性ラ音、類鼾音は低音性連続性ラ音と呼ばれるように、笛声音と類鼾音とは音の高さ（即ち、周波数）で判別することが可能である。しかしながら、笛声音及び類鼾音は、ピーク周波数が時間的に変化する。このため、図に示すような単一の閾値（即ち、値が変動しない一つの閾値）を利用して笛声音及び類鼾音を判定しようとする比較例では、時間の経過により、判定結果が変化してしまうことがある。例えば、図中に示す笛声音成分のように、ピーク周波数が判定閾値を跨ぐように変化してしまうと、それまでは正確に笛声音と判定されていたものが、誤って類鼾音として判定されることになってしまう。このため本実施例に係る呼吸音解析装置では、ピーク周波数に応じて閾値を変動させる。

【0051】

図４及び図５に示すように、本実施例に係る呼吸音解析装置では、笛声音と判定する割合及び類鼾音と判定する割合がピーク周波数に応じてなめらかに変化するように閾値が変動する。例えば、ピーク周波数が２００Ｈｚの場合には、笛声音が７％含まれ、類鼾音が９３％含まれると判定する。ピーク周波数が２５０Ｈｚの場合には、笛声音が５０％含まれ、類鼾音が５０％含まれると判定する。ピーク周波数が２８０Ｈｚの場合には、笛声音が７８％含まれ、類鼾音が２２％含まれると判定する。なお、ここでの具体的な数値はあくまで一例であり、異なる値を設定してもよい。また、測定対象である生体の性別、年齢、身長、体重等によって異なる変動特性を有するようにしてもよい。

【0052】

上述した変動する閾値を利用することで、図３で示したような場合に発生し得る誤判定を好適に防止することができる。即ち、本実施例に係る呼吸音解析装置では、笛声音及び類鼾音を判定するための閾値がピーク周波数に応じて適切な値にとなるよう変動するため、例えば変動しない単一の閾値を用いる場合と比較して、より正確な判定が行える。

【0053】

図２に戻り、笛声音／類鼾音判定部２４０では、上述したように変動する閾値を用いて、連続性ラ音が笛声音であるのか、又は類鼾音であるのかが判定される（ステップＳ１０７）。そして、画像生成部２５０では判定結果に基づく画像が生成され、表示部１４０において判定結果が表示される（ステップＳ１０８）。なお、解析結果を表示した後は、解析を継続するか否かが判定される（ステップＳ１０９）。解析を継続すると判定された場合（ステップＳ１０９：ＹＥＳ）、上述した処理がステップＳ１０１から再開される。一方、解析を継続しないと判定された場合（ステップＳ１０９：ＮＯ）、一連の処理は終了する。

【0054】

ここで、表示部１４０における解析結果の表示について、図６及び図７を参照して具体的に説明する。ここに図６及び図７は夫々、本実施例に係る連続性ラ音判定結果の表示例を示すグラフである。

【0055】

図６に示すように、解析結果である笛声音及び類鼾音の割合は、棒グラフとして表示部１５０に表示される。ただし、この表示方法は一例であり、他の表示態様で表示を行ってもよい。例えば、笛声音及び類鼾音の割合を円グラフとして表示してもよい。或いは、笛声音及び類鼾音の強度を数値化して表示してもよい。

【0056】

図７に示すように、笛声音及び類鼾音以外の音種（例えば、正常呼吸音、水泡音、捻髪音等）についても判定可能な場合には、それらの音種の割合も合わせて表示するようにしてもよい。

【0057】

なお、上述した画像としての出力に代えて或いは加えて、音声データによる出力も可能である。具体的には、笛声音と類鼾音とで別々に音声を出力することができる。或いは、笛声音及び類鼾音の一方だけを強調して音声を出力することができる。

【0058】

＜判定結果の具体例＞
最後に、図８及び図９において解析結果の具体例を挙げて、本実施例の利点について詳細に説明する。ここに図８は、本実施例及び比較例に係る呼吸音解析装置により笛声音を解析する例を示す概念図である。また図９は、本実施例及び比較例に係る呼吸音解析装置により類鼾音を解析する例を示す概念図である。なお、以下では、変動する閾値（図４及び図５を参照）を用いる本実施例に係る呼吸音解析装置、と、単一閾値（２５０Ｈｚ）を用いる比較例に係る呼吸音解析装置とで、笛声音及び類鼾音の各々スペクトルを解析する場合を考える。なお、笛声音及び類鼾音の判定は、所定の時間間隔ごと（時刻ｔ１〜ｔ７）で行われるものとする。

【0059】

図８に示すように、ピーク周波数が約３５０Ｈｚから約２００Ｈｚまで低下するような笛声音が解析対象であるとする。この場合、比較例に係る呼吸音解析装置では、時刻ｔ１からｔ６においては、解析対象である連続性ラ音が笛声音であると正確に判定されている。しかしながら、時刻ｔ７においては、笛声音のスペクトルのピーク周波数が低下している（即ち、閾値である２５０Ｈｚを下回っている）ため、類鼾音と誤って判定されている。

【0060】

一方、本実施例に係る呼吸音解析装置では、時刻ｔ１からｔ４においては、解析対象である連続性ラ音が１００％笛声音であると正確に判定されている。そして、時刻ｔ５からは類鼾音の割合が少しずつ増加するが、時刻ｔ７においても笛声音が含まれるものと判定され続ける。即ち、時刻ｔ７において完全に誤った判定をしてしまった比較例に対して、本実施例は笛声音も含まれる可能性があるという結果を提示することができる。

【0061】

図９に示すように、ピーク周波数が約２４０Ｈｚから約２６０Ｈｚまで上昇し、その後約１８０Ｈｚまで低下するような類鼾音が解析対象であるとする。この場合、比較例に係る呼吸音解析装置では、時刻ｔ１、及びｔ４からｔ７において、解析対象である連続性ラ音が類鼾音であると正確に判定されている。しかしながら、時刻ｔ２及びｔ３においては、類鼾音のスペクトルのピーク周波数が一時的に上昇している（即ち、閾値である２５０Ｈｚを上回っている）ため、笛声音と誤って判定されている。

【0062】

一方、本実施例に係る呼吸音解析装置では、ピーク周波数が上昇する時刻ｔ２及びｔ３においても類鼾音が含まれるものと判定され続ける。即ち、時刻ｔ２及びｔ３において完全に誤った判定をしてしまった比較例に対して、本実施例は類鼾音も含まれる可能性があるという結果を提示することができる。

【0063】

以上説明したように、本実施例に係る呼吸音解析装置によれば、ピーク周波数に応じて変動する閾値を利用することで、より好適に笛声音及び類鼾音を判定することができる。

【0064】

本発明は、上述した実施形態に限られるものではなく、特許請求の範囲及び明細書全体から読み取れる発明の要旨或いは思想に反しない範囲で適宜変更可能であり、そのような変更を伴う呼吸音解析装置及び呼吸音解析方法、並びにコンピュータプログラム及び記録媒体もまた本発明の技術的範囲に含まれるものである。

【符号の説明】

【0065】

１１０ 生体音センサ
１２０ 信号記憶部
１２５ 信号処理部
１３０ 音声出力部
１４０ 表示部
２００ 処理部
２１０ 周波数解析部
２２０ 連続性ラ音検出部
２３０ ピーク周波数検出部
２４０ 笛声音／類鼾音判定部
２５０ 画像生成部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】