JP4759727B2 - デジタル聴診解析システム - Google Patents
デジタル聴診解析システム Download PDFInfo
- Publication number
- JP4759727B2 JP4759727B2 JP2005080720A JP2005080720A JP4759727B2 JP 4759727 B2 JP4759727 B2 JP 4759727B2 JP 2005080720 A JP2005080720 A JP 2005080720A JP 2005080720 A JP2005080720 A JP 2005080720A JP 4759727 B2 JP4759727 B2 JP 4759727B2
- Authority
- JP
- Japan
- Prior art keywords
- analysis
- sound
- evaluation index
- waveform
- vibration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Description
本発明は、聴診音の異常を診断するデジタル聴診解析システムに係り、特に、心音データから特徴的なパターンを作成して解析するデジタル聴診解析システムに関する。
近年、日本では心臓病による死亡率が急速に増加しており、１９８５年からは脳卒中を抜いて第二位になっている。そして、心臓病の中では心不全や虚血性心疾患といわれる心筋梗塞による死亡が最も多く、また、原因不明の急性死が三割を占めている。このような心臓病を診断する代表的な方法として心音解析が挙げられるが、この心音解析に関する研究は心エコー図法や心臓カテーテルのような正確な検査方法の確立によって減少傾向にあった。しかしながら、近年のコンピュータの発達とデジタル信号処理技術の進歩によって、心音の録音や解析が容易にできるようになり、熟練を要する心音の聴診診断分野においても心音を自動診断する装置等が提案されている。
例えば、特許文献１には、患者の聴診音を入力する聴診マイクと、入力される聴診音の体外部雑音を濾過する濾過器と、聴診音のアナログ波形データをデジタルデータに変換するデジタル変換器と、聴診音のデジタルデータと予め入力された各種疾病の聴診音の標準データとを比較するコンピュータと、コンピュータの出力を記録する記録装置及び表示する複数のモニタから構成される自動判読記録診断装置が開示されている。
この特許文献１に開示された発明では、患者の聴診音がデジタル化され、このデジタルデータはコンピュータによって各種疾病の聴診音の標準データと比較されて病名が判読されるので、医師の個人別読力や主観的判断を介入せずに迅速で正確な病名を決定することができる。また、複数のモニタの設置によって、医師だけでなく、患者、看護師、医大生等の大勢の人が同時に見ることができるので対処診療の協同と迅速性を期することが可能である。
この特許文献１に開示された発明では、患者の聴診音がデジタル化され、このデジタルデータはコンピュータによって各種疾病の聴診音の標準データと比較されて病名が判読されるので、医師の個人別読力や主観的判断を介入せずに迅速で正確な病名を決定することができる。また、複数のモニタの設置によって、医師だけでなく、患者、看護師、医大生等の大勢の人が同時に見ることができるので対処診療の協同と迅速性を期することが可能である。
また、特許文献２には、聴診器に設置され受信した音声を被測定信号に変換して出力する第１の変換手段と、聴診器に設置され制御操作子の操作状態に基づいて単一周波数の制御信号を出力する制御信号発生手段と、制御信号によって操作され被測定信号を分析して記録する診断装置と、入力信号の周波数成分毎のレベルを検出する周波数成分検出手段と、周波数成分に基づいて入力信号と制御信号及び被測定信号を比較判定する判定手段とを具備する電子式聴診・診断装置が開示されている。
この特許文献２に開示された発明では、聴診器が取得した信号音はマイコン装置によって処理されて他の複数の聴診器においても聴診が可能となり、また、判読の結果が作成されるので診断の正確性を上げることができる。さらに、資料の記録計算や累積や伝送出力も可能になっている。
この特許文献２に開示された発明では、聴診器が取得した信号音はマイコン装置によって処理されて他の複数の聴診器においても聴診が可能となり、また、判読の結果が作成されるので診断の正確性を上げることができる。さらに、資料の記録計算や累積や伝送出力も可能になっている。
そして、聴診器について、特許文献３には、聴診音を検出するステート部と、ステート部に設置されスイッチング操作に応じて制御信号を送るスイッチと、聴診音を音信号に変換する変換手段と、制御信号に基づいて音信号の記録処理又は音信号の外部装置への伝送処理を行う信号処理装置を備えた聴診器が開示されている。
この特許文献３に開示された発明では、ステート部にスイッチが設置されており、このスイッチのスイッチング操作によって発生する制御信号に基づいて聴診音の音信号の録音や伝送を行うことができるので、看護師等の聴診を行う者は、聴診器の基本的な操作に加えて聴診音の音信号の録音又は伝送の操作を片手のみで容易に行うことができる。
この特許文献３に開示された発明では、ステート部にスイッチが設置されており、このスイッチのスイッチング操作によって発生する制御信号に基づいて聴診音の音信号の録音や伝送を行うことができるので、看護師等の聴診を行う者は、聴診器の基本的な操作に加えて聴診音の音信号の録音又は伝送の操作を片手のみで容易に行うことができる。
さらに、特許文献４には、聴診音を採取して聴診音信号を出力する聴診手段と、聴診音信号に圧縮符号化処理を行ってファイル化するデータ処理手段と、ファイル化された聴診音信号を通信手段を介して送信する送信手段を備えた聴診システムが開示されている。
この特許文献４に開示された発明では、聴診音信号を比較的安価に送信することができるので、訪問医療や在宅看護等の医療機関以外の看護や介護において使用すると医療の質を低下させることなく人件費の効率的な抑制が可能となる。
この特許文献４に開示された発明では、聴診音信号を比較的安価に送信することができるので、訪問医療や在宅看護等の医療機関以外の看護や介護において使用すると医療の質を低下させることなく人件費の効率的な抑制が可能となる。
しかしながら、特許文献１に記載された従来の技術では、患者から測定される聴診音に基づいてコンピュータが詳しい病名を判読するが、この判読は立会する医師の最終的な判断をより精確に行うためのものであり、一般家庭において個人が健康管理の目的で使用するには内容が高度であるという課題があった。
また、特許文献２に記載された従来の技術は、特許文献１の場合と同様に、医師等の専門家の診察を補助するためのものであって一般家庭で使用するには複雑な構成になっているという課題があった。また、病状の判定を行う病徴判定器についても記載されているが詳しい判定方法は開示されていない。
そして、特許文献３及び特許文献４に記載された従来の技術では、聴診器によって測定される聴診音を録音したり伝送したり送信したりすることはできるが、測定される聴診音を解析するものではない。
特開２００２−１６５７８９号公報
特公平４−３２６６１号公報
特開２００２−１５３４５９号公報
特開２００１−３３３８９９号公報
本発明はかかる従来の事情に対処してなされたものであり、一般家庭で簡単に心音診断が可能なデジタル聴診解析システムを提供することを目的とする。
上記目的を達成するため、請求項１記載の発明であるデジタル聴診解析システムは、解析入力情報格納部と、振動モデルを用いて特徴値波形を解析する心音波形解析部と、解析された前記特徴値波形を用いて評価指数を解析する評価指数解析部と、異常心音評価部と、解析結果情報格納部とを有するデジタル聴診解析システムであって、
解析入力情報格納部は、デジタル化された心音データと振動モデルの解析パラメータと評価指数の正常値範囲情報を読み出し可能に格納し、心音波形解析部は、心音データと解析パラメータを読み出して振動モデルを用いて心音の特徴値波形を解析する手段と、この特徴値波形を解析結果情報格納部に格納する手段とを備え、
振動モデルは、質量にばね系及び粘性減衰系を接続させた１の機械的振動系又はインダクタンスに電気容量と抵抗を接続させた電気的振動系からなる、あるいはこれらの振動系を複数連成させてなる鼓膜の振動応答に関する振動モデルであり、
解析入力情報格納部から読み出した解析パラメータ及び心音データを振動モデルに代入して鼓膜の振動応答を得て、これを特徴値波形とし、
評価指数は、特徴値波形から得られる心臓の僧帽弁、三尖弁、大動脈弁、肺動脈弁のいずれかの閉鎖音の持続時間、あるいはこれらの任意の弁に対して同一弁あるいは異なる２弁の閉鎖音間の持続時間、あるいはそれらの組合せとし、
評価指数解析部は、解析結果情報格納部から特徴値波形を読み出してこの特徴値波形から評価指数を演算する手段と、評価指数を解析結果情報格納部へ格納する手段とを備え、
異常心音評価部は、解析結果情報格納部から評価指数を読み出して解析入力情報格納部から読み出した評価指数の正常値範囲情報と比較する手段と、評価指数が正常値範囲内から外れる場合に警報情報信号を発する手段とを備えることを特徴とするデジタル聴診解析システムである。
上記構成のデジタル聴診解析システムでは、心音波形解析部において解析入力情報格納部に格納される心音データと解析パラメータが読み出され振動モデルに基づいて心音の特徴値波形が作成されるという作用を有する。
また、鼓膜の振動モデルを模擬した振動系を採用することで、より特徴的な心音の特徴値波形を作成するという作用を有する。
さらに、心音波形解析部において解析入力情報格納部に格納される心音データと解析パラメータが読み出され振動モデルに基づいて心音の特徴値波形が作成され、続いて評価指数解析部において、解析結果情報格納部から特徴値波形が読み出されて評価指数が演算される。そして、異常心音評価部において、解析結果情報格納部及び解析入力情報格納部から各々読み出される評価指数と評価指数の正常値範囲情報とが比較されて評価指数が正常値範囲内から外れる場合に警報情報信号が発せられるという作用を有する。
解析入力情報格納部は、デジタル化された心音データと振動モデルの解析パラメータと評価指数の正常値範囲情報を読み出し可能に格納し、心音波形解析部は、心音データと解析パラメータを読み出して振動モデルを用いて心音の特徴値波形を解析する手段と、この特徴値波形を解析結果情報格納部に格納する手段とを備え、
振動モデルは、質量にばね系及び粘性減衰系を接続させた１の機械的振動系又はインダクタンスに電気容量と抵抗を接続させた電気的振動系からなる、あるいはこれらの振動系を複数連成させてなる鼓膜の振動応答に関する振動モデルであり、
解析入力情報格納部から読み出した解析パラメータ及び心音データを振動モデルに代入して鼓膜の振動応答を得て、これを特徴値波形とし、
評価指数は、特徴値波形から得られる心臓の僧帽弁、三尖弁、大動脈弁、肺動脈弁のいずれかの閉鎖音の持続時間、あるいはこれらの任意の弁に対して同一弁あるいは異なる２弁の閉鎖音間の持続時間、あるいはそれらの組合せとし、
評価指数解析部は、解析結果情報格納部から特徴値波形を読み出してこの特徴値波形から評価指数を演算する手段と、評価指数を解析結果情報格納部へ格納する手段とを備え、
異常心音評価部は、解析結果情報格納部から評価指数を読み出して解析入力情報格納部から読み出した評価指数の正常値範囲情報と比較する手段と、評価指数が正常値範囲内から外れる場合に警報情報信号を発する手段とを備えることを特徴とするデジタル聴診解析システムである。
上記構成のデジタル聴診解析システムでは、心音波形解析部において解析入力情報格納部に格納される心音データと解析パラメータが読み出され振動モデルに基づいて心音の特徴値波形が作成されるという作用を有する。
また、鼓膜の振動モデルを模擬した振動系を採用することで、より特徴的な心音の特徴値波形を作成するという作用を有する。
さらに、心音波形解析部において解析入力情報格納部に格納される心音データと解析パラメータが読み出され振動モデルに基づいて心音の特徴値波形が作成され、続いて評価指数解析部において、解析結果情報格納部から特徴値波形が読み出されて評価指数が演算される。そして、異常心音評価部において、解析結果情報格納部及び解析入力情報格納部から各々読み出される評価指数と評価指数の正常値範囲情報とが比較されて評価指数が正常値範囲内から外れる場合に警報情報信号が発せられるという作用を有する。
請求項２の発明であるデジタル聴診解析システムは、請求項１に記載のデジタル聴診解析システムにおいて、特徴値波形、評価指数、警報情報の少なくとも１を出力及び／又は表示する表示部を有するものである。
本構成に係るデジタル聴診解析システムにおいては、特徴値波形、評価指数、警報情報を出力して、これらの情報に関して使用者に対して認知させたり、あるいはさらに情報を処理するために加工するという作用を有する。また、表示部は、使用者にこれらの情報を視認可能とする作用を有する。
本構成に係るデジタル聴診解析システムにおいては、特徴値波形、評価指数、警報情報を出力して、これらの情報に関して使用者に対して認知させたり、あるいはさらに情報を処理するために加工するという作用を有する。また、表示部は、使用者にこれらの情報を視認可能とする作用を有する。
本発明の請求項１及び請求項２に記載のデジタル聴診解析システムは、デジタル化した心音データから振動モデルを用いて特徴値波形を作成することができるので、心音データそのものよりも客観的な判断を行なうことが可能となる。
また、特徴値波形から演算される評価指数とこの評価指数の正常値範囲とを比較して異常を検知することができる。従って、専門的で詳細な病状診断はできないものの、評価指数を導入することで定量的で客観的な判断が可能であり、特に専門知識を有していない人でも扱い易く、日常の体調管理を手軽に行うことができる。
また、本発明に係るデジタル聴診解析システムは、簡単な構成であるので安価で提供できる。
また、特徴値波形から演算される評価指数とこの評価指数の正常値範囲とを比較して異常を検知することができる。従って、専門的で詳細な病状診断はできないものの、評価指数を導入することで定量的で客観的な判断が可能であり、特に専門知識を有していない人でも扱い易く、日常の体調管理を手軽に行うことができる。
また、本発明に係るデジタル聴診解析システムは、簡単な構成であるので安価で提供できる。
さらに、鼓膜の振動モデルを模擬した比較的単純な振動系を採用することで、容易に特徴的な心音成分を抽出することができる。また、心臓弁（僧帽弁、三尖弁、大動脈弁、肺動脈弁）の閉鎖音に着目して、これらの持続時間やこれらの閉鎖音間の持続時間等、特に心臓疾患に関連して変動し易い数値を評価指数として採用することで、より客観的で精度の高い診断を容易に実施することが可能である。
以下に、本発明に係るデジタル聴診解析システムの実施の形態を図１乃至図１５に基づき説明する。（請求項１及び請求項２に対応）
まず、最初に、本実施の形態に係るデジタル聴診解析システムの設置環境について図１を用いて説明する。
図１（ａ）は、本実施の形態に係るデジタル聴診解析システムを用いた無線聴診システムの送信部の構成図であり、（ｂ）は同じくデジタル聴診解析システムを用いた無線聴診システムの受信部の構成図である。
図１（ａ）において、無線聴診システム送信部１は、マイク入力２から採取された聴診音は、トリガー無線送信機３を介してプリアンプ４で音の補正が行われた後、ノイズフィルタ５によってノイズが除去され、オーディオ無線送信機６において無線信号として図１（ｂ）に示す無線聴診システム受信部１２に送信される構成になっている。
なお、トリガー無線送信機３では、聴診音に心音の変化等が検出されたときに無線送信を開始するように設定することが可能である。また、ノイズフィルタ５には、例えば、心拍（３０Ｈｚ〜４００Ｈｚ）や呼吸音（２０Ｈｚ〜１５００Ｈｚ）に対応するフィルタが設置される。
そして、ノイズフィルタ５を経た聴診音は、オーディオ無線送信機６以外にもパワーアンプ７とマイクロコントローラ１０に分岐されており、パワーアンプ７では音が増幅された後、イヤホン／コンピュータ８とＩＣレコーダ９に接続され、イヤホン／コンピュータ８では聴診音を聴いたり、保存したり、また、ＩＣレコーダ９では聴診音を記録したりすることができる。
また、マイクロコントローラ１０には外部メモリ１１が接続されており、マイクロコントローラ１０は外部メモリ１１を起動させて聴診音を外部メモリ１１に格納するために作動する。なお、外部メモリ１１には記録可能な記録媒体であれば特に限定されず、例えばメモリカード等が設置可能である。
まず、最初に、本実施の形態に係るデジタル聴診解析システムの設置環境について図１を用いて説明する。
図１（ａ）は、本実施の形態に係るデジタル聴診解析システムを用いた無線聴診システムの送信部の構成図であり、（ｂ）は同じくデジタル聴診解析システムを用いた無線聴診システムの受信部の構成図である。
図１（ａ）において、無線聴診システム送信部１は、マイク入力２から採取された聴診音は、トリガー無線送信機３を介してプリアンプ４で音の補正が行われた後、ノイズフィルタ５によってノイズが除去され、オーディオ無線送信機６において無線信号として図１（ｂ）に示す無線聴診システム受信部１２に送信される構成になっている。
なお、トリガー無線送信機３では、聴診音に心音の変化等が検出されたときに無線送信を開始するように設定することが可能である。また、ノイズフィルタ５には、例えば、心拍（３０Ｈｚ〜４００Ｈｚ）や呼吸音（２０Ｈｚ〜１５００Ｈｚ）に対応するフィルタが設置される。
そして、ノイズフィルタ５を経た聴診音は、オーディオ無線送信機６以外にもパワーアンプ７とマイクロコントローラ１０に分岐されており、パワーアンプ７では音が増幅された後、イヤホン／コンピュータ８とＩＣレコーダ９に接続され、イヤホン／コンピュータ８では聴診音を聴いたり、保存したり、また、ＩＣレコーダ９では聴診音を記録したりすることができる。
また、マイクロコントローラ１０には外部メモリ１１が接続されており、マイクロコントローラ１０は外部メモリ１１を起動させて聴診音を外部メモリ１１に格納するために作動する。なお、外部メモリ１１には記録可能な記録媒体であれば特に限定されず、例えばメモリカード等が設置可能である。
一方、図１（ｂ）において、無線聴診システム受信部１２では、無線聴診システム送信部１から送信される聴診音の無線信号が、オーディオ無線受信機１３とトリガー無線受信機１４で受信され、プリアンプ４で補正が行われて、ノイズフィルタ５でノイズが除去されると、オートチューニングアンプ１５とパワーアンプ７に分岐される。
オートチューニングアンプ１５は、送受信によって歪んだ無線信号を解析に適した状態まで増幅するものである。まず、受信される無線信号はオートチューニングアンプ１５を介してマイクロコントローラ１０に送られ、マイクロコントローラ１０において適切な解析信号として増幅度が達成されていないと判断した場合には、オートチューニングアンプ１５に増幅度をチューニングするように指令が出される。そして、この指令を受信したオートチューニングアンプ１５では指令に従って信号を解析に適切な増幅度に調節して、コンピュータ１８や外部メモリ１１に送信する。
なお、マイクロコントローラ１０は、外部メモリ１１が接続される場合には、外部メモリ１１を起動させて外部メモリ１１に信号を格納するために作動したり、さらには、聴診音のＡＤ変換を行ったりする。
また、コンピュータ１８には、本実施の形態に係るデジタル聴診解析システムが搭載されており、後述する聴診音による心音解析を行うことができるようになっている。
そして、パワーアンプ７では、音が増幅された後、イヤホン１６で音を聴いたり、オーディオ端子を有するコンピュータ１７で保存したり、ＩＣレコーダ９で記録したりすることができる。
なお、コンピュータ１７についてもデジタル聴診解析システムを搭載して心音解析を行うことが可能であるが、この場合、オートチューニングアンプ１５による無線信号の増幅度の調整はパワーアンプ７を用いて手動で設定する必要がある。また、マイクロコントローラ１０を介していないので、コンピュータ１７内で解析を行う前にＡＤ変換を施さなければならない。
オートチューニングアンプ１５は、送受信によって歪んだ無線信号を解析に適した状態まで増幅するものである。まず、受信される無線信号はオートチューニングアンプ１５を介してマイクロコントローラ１０に送られ、マイクロコントローラ１０において適切な解析信号として増幅度が達成されていないと判断した場合には、オートチューニングアンプ１５に増幅度をチューニングするように指令が出される。そして、この指令を受信したオートチューニングアンプ１５では指令に従って信号を解析に適切な増幅度に調節して、コンピュータ１８や外部メモリ１１に送信する。
なお、マイクロコントローラ１０は、外部メモリ１１が接続される場合には、外部メモリ１１を起動させて外部メモリ１１に信号を格納するために作動したり、さらには、聴診音のＡＤ変換を行ったりする。
また、コンピュータ１８には、本実施の形態に係るデジタル聴診解析システムが搭載されており、後述する聴診音による心音解析を行うことができるようになっている。
そして、パワーアンプ７では、音が増幅された後、イヤホン１６で音を聴いたり、オーディオ端子を有するコンピュータ１７で保存したり、ＩＣレコーダ９で記録したりすることができる。
なお、コンピュータ１７についてもデジタル聴診解析システムを搭載して心音解析を行うことが可能であるが、この場合、オートチューニングアンプ１５による無線信号の増幅度の調整はパワーアンプ７を用いて手動で設定する必要がある。また、マイクロコントローラ１０を介していないので、コンピュータ１７内で解析を行う前にＡＤ変換を施さなければならない。
このように構成された無線聴診システムにおいては、採取される聴診音は録音や保存が可能であるので、繰り返し再生することができ、また、患者の病状を管理する上で便利である。そして、聴診音を転送可能に設計されているので、聴診音を採取する場所と離れた場所で解析することができる。
なお、無線聴診システムの送信部及び受信部において複数の構成が並列に設けられているが、全ての構成を設置する必要はなく、いずれか１の構成又はこれらの組合せであってもよい。
また、この無線聴診システムにおいては、無線によって聴診音を転送しているが、転送方法は無線に限定されるものではなく、有線を利用することもできる。
なお、無線聴診システムの送信部及び受信部において複数の構成が並列に設けられているが、全ての構成を設置する必要はなく、いずれか１の構成又はこれらの組合せであってもよい。
また、この無線聴診システムにおいては、無線によって聴診音を転送しているが、転送方法は無線に限定されるものではなく、有線を利用することもできる。
以上のような無線聴診システムなどに搭載可能な本実施の形態に係るデジタル聴診解析システムの構成について図２を用いて説明する。
図２は、本実施の形態に係るデジタル聴診解析システムの構成図である。
図２において、デジタル聴診解析システム１９は、特徴値波形解析部２１、評価指数解析部２２及び異常心音評価部２３から構成される演算部２０と、表示部２４と、モデルパラメータ入力部２５と、心音データ２７、モデルパラメータ２８及び正常値範囲データ２９を格納する解析入力データベース２６と、特徴値波形３１、Ｔ１データ３２、Ｔ２データ３３、Ｔ１１データ３４及びＴ１２データ３５を格納する解析結果データベース３０とを有している。
特徴値波形解析部２１では、まず、解析入力データベース２６から心音データ２７とモデルパラメータ２８を読み出し、この読み出したモデルパラメータ２８と心音データ２７を、詳しくは後述するが、鼓膜の振動モデルに対応させて、心臓の僧帽弁と三尖弁の閉鎖音であるＩ音情報と大動脈弁と肺動脈弁の閉鎖音であるII音情報等、心音を構成する音情報を含む特徴的な心音波形（以下、特徴値波形という）を作成する。そして、作成した特徴値波形３１を解析結果データベース３０に格納する。
なお、解析入力データベース２６では、心音データ２７は被験者の氏名や病状等の個人情報や測定日等とともに格納されている。また、モデルパラメータ２８は、後述するが、固有振動数及び減衰比係数の各々についてデータテーブルが作成されて格納されている。
モデルパラメータ２８はモデルパラメータ入力部２５から解析入力データベース２６に格納され、特徴値波形解析部２１によって読み出されるほか、モデルパラメータ入力部２５から直接特徴値波形解析部２１に入力することも可能であり、その数値は所望の値に設定することができる。
また、心音データ２７は基本的にはデジタル化されたものであるが、特徴値波形解析部２１においてＡＤ変換を行うようにすればアナログの心音データを使用することも可能である。
そして、解析結果データベース３０では、作成された特徴値波形３１は使用した心音データ及びモデルパラメータに関連付けられて格納されている。
図２は、本実施の形態に係るデジタル聴診解析システムの構成図である。
図２において、デジタル聴診解析システム１９は、特徴値波形解析部２１、評価指数解析部２２及び異常心音評価部２３から構成される演算部２０と、表示部２４と、モデルパラメータ入力部２５と、心音データ２７、モデルパラメータ２８及び正常値範囲データ２９を格納する解析入力データベース２６と、特徴値波形３１、Ｔ１データ３２、Ｔ２データ３３、Ｔ１１データ３４及びＴ１２データ３５を格納する解析結果データベース３０とを有している。
特徴値波形解析部２１では、まず、解析入力データベース２６から心音データ２７とモデルパラメータ２８を読み出し、この読み出したモデルパラメータ２８と心音データ２７を、詳しくは後述するが、鼓膜の振動モデルに対応させて、心臓の僧帽弁と三尖弁の閉鎖音であるＩ音情報と大動脈弁と肺動脈弁の閉鎖音であるII音情報等、心音を構成する音情報を含む特徴的な心音波形（以下、特徴値波形という）を作成する。そして、作成した特徴値波形３１を解析結果データベース３０に格納する。
なお、解析入力データベース２６では、心音データ２７は被験者の氏名や病状等の個人情報や測定日等とともに格納されている。また、モデルパラメータ２８は、後述するが、固有振動数及び減衰比係数の各々についてデータテーブルが作成されて格納されている。
モデルパラメータ２８はモデルパラメータ入力部２５から解析入力データベース２６に格納され、特徴値波形解析部２１によって読み出されるほか、モデルパラメータ入力部２５から直接特徴値波形解析部２１に入力することも可能であり、その数値は所望の値に設定することができる。
また、心音データ２７は基本的にはデジタル化されたものであるが、特徴値波形解析部２１においてＡＤ変換を行うようにすればアナログの心音データを使用することも可能である。
そして、解析結果データベース３０では、作成された特徴値波形３１は使用した心音データ及びモデルパラメータに関連付けられて格納されている。
次に、評価指数解析部２２では、解析結果データベース３０に格納された特徴値波形３１を読み出して、特徴値波形３１に含まれる複数のＩ音情報及びII音情報を抽出する。そして、Ｉ音情報及びII音情報から各々Ｉ音ピーク及びII音ピークを選定し、これらのピークについて評価指数としてＩ音の持続時間Ｔ１、II音の持続時間Ｔ２、Ｉ音からII音までの持続時間Ｔ１２及びＩ音から次のＩ音までの持続時間Ｔ１１を演算し、演算した持続時間をデータとして各々解析結果データベース３０に格納する。
なお、このＩ音ピークやII音ピーク、あるいはＴ１，Ｔ２，Ｔ３，Ｔ４については図６を参照しながら後で説明する。
また、評価指数解析部２２では、演算した評価指数を用いて、例えば、Ｔ１とＴ２の関係やＴ１１とＴ１２の関係を示すグラフを作成したりすることもできる。
なお、解析結果データベース３０では、演算された評価指数は使用した特徴値波形と関連付けられて、各々がデータテーブルとして格納されている。
なお、このＩ音ピークやII音ピーク、あるいはＴ１，Ｔ２，Ｔ３，Ｔ４については図６を参照しながら後で説明する。
また、評価指数解析部２２では、演算した評価指数を用いて、例えば、Ｔ１とＴ２の関係やＴ１１とＴ１２の関係を示すグラフを作成したりすることもできる。
なお、解析結果データベース３０では、演算された評価指数は使用した特徴値波形と関連付けられて、各々がデータテーブルとして格納されている。
続いて、異常心音評価部２３では、まず、解析結果データベース３０に格納されている４個の評価指数であるＴ１データ３２、Ｔ２データ３３、Ｔ１１データ３４及びＴ１２データ３５を読み出す。続いて、解析入力情報データベース２６に格納されている正常値範囲データ２９から４個の評価指数各々に対する正常値範囲を読み出す。そして、４個の評価指数を順次対応する正常値範囲と比較して、正常値範囲を外れる評価指数がある場合は、警報情報信号を生成する。
なお、この警報情報信号とは警報情報を伝達する信号であり、警報情報には警報メッセージ、警報音及び警報ランプ等の視覚又は聴覚に訴える全ての情報を含むものである。
また、解析入力データベース２６では、４個の評価指数に対応する正常値範囲は各々データテーブルが作成されて正常値範囲データ２９に格納されている。
なお、この警報情報信号とは警報情報を伝達する信号であり、警報情報には警報メッセージ、警報音及び警報ランプ等の視覚又は聴覚に訴える全ての情報を含むものである。
また、解析入力データベース２６では、４個の評価指数に対応する正常値範囲は各々データテーブルが作成されて正常値範囲データ２９に格納されている。
そして、表示部２４では、演算部２０において演算される特徴値波形や、評価指数及び評価指数に基づくグラフ図を表示することができる。また、異常心音評価部において作成される警報メッセージや警報ランプを表示したり、警報音を鳴らしたりすることもできる。
ここで、心音について図３を参照しながら説明する。図３は、人の心臓の概念図である。
図３において、心臓３６は、左心房３７、左心室３８、右心房４１及び右心室４２の四つの部分に分かれており、左心房３７及び左心室３８では肺でガス交換した後の酸素をたくさん含んだ血液を全身に送る働きをし、右心房４１及び右心室４２では全身から戻る二酸化炭素を多く含む血液を肺に送る働きをしており、全体で収縮と弛緩を繰り返すことで血液を全身に循環させるポンプの役割を担っている。
また、左心房３７の入り口には僧帽弁３９、左心室３８の入り口には大動脈弁４０、右心室４２の入り口には三尖弁４３及び右心房４１の入り口には肺動脈弁４４と呼ばれる弁膜が各々ついており、これらの弁膜によって血液の逆流が防止されている。そして、心音とはこれらの弁膜が閉じる際に発する音である。
聴診による心音が正常であるか異常であるか、又は、心雑音があるのかないのかについては専門の知識を有する医師でないと厳密に聞き分けて診断するのは困難である。しかしながら、専門知識を有さない素人であっても、「ドック、ドック」と聞こえる正常な心音と、「グー、グー」と聞こえる僧帽弁３９の閉鎖不全による心雑音など比較的簡単に聞き分けることができるものもある。これは、鼓膜が聴診器からの音波で拍動された際に発生する低次モードの振動を人間の耳で感じ取りやすいためと考えられる。この発想に基づいて、本願発明の発明者は鼓膜を振動モデルに近似し、聴診器から採取した心音と振動モデルの振動応答との関係から心音の解析を行うことを試みた。
図３において、心臓３６は、左心房３７、左心室３８、右心房４１及び右心室４２の四つの部分に分かれており、左心房３７及び左心室３８では肺でガス交換した後の酸素をたくさん含んだ血液を全身に送る働きをし、右心房４１及び右心室４２では全身から戻る二酸化炭素を多く含む血液を肺に送る働きをしており、全体で収縮と弛緩を繰り返すことで血液を全身に循環させるポンプの役割を担っている。
また、左心房３７の入り口には僧帽弁３９、左心室３８の入り口には大動脈弁４０、右心室４２の入り口には三尖弁４３及び右心房４１の入り口には肺動脈弁４４と呼ばれる弁膜が各々ついており、これらの弁膜によって血液の逆流が防止されている。そして、心音とはこれらの弁膜が閉じる際に発する音である。
聴診による心音が正常であるか異常であるか、又は、心雑音があるのかないのかについては専門の知識を有する医師でないと厳密に聞き分けて診断するのは困難である。しかしながら、専門知識を有さない素人であっても、「ドック、ドック」と聞こえる正常な心音と、「グー、グー」と聞こえる僧帽弁３９の閉鎖不全による心雑音など比較的簡単に聞き分けることができるものもある。これは、鼓膜が聴診器からの音波で拍動された際に発生する低次モードの振動を人間の耳で感じ取りやすいためと考えられる。この発想に基づいて、本願発明の発明者は鼓膜を振動モデルに近似し、聴診器から採取した心音と振動モデルの振動応答との関係から心音の解析を行うことを試みた。
次に、鼓膜の振動モデルについて図４を用いて説明する。図４は、鼓膜の振動モデルを示す概念図である。
図４において、符号４５は鼓膜、符号４６はばね、そして、符号４７はダンパーであり、鼓膜の等価質量をｍ、ばね定数をｋ、粘性減衰係数をｃ及び聴診器からの心音をＳとすると、鼓膜の振動応答ｘは式（１）から算出される。
図４において、符号４５は鼓膜、符号４６はばね、そして、符号４７はダンパーであり、鼓膜の等価質量をｍ、ばね定数をｋ、粘性減衰係数をｃ及び聴診器からの心音をＳとすると、鼓膜の振動応答ｘは式（１）から算出される。
また、式（１）の両辺をｍで割り、さらに、心音信号入力値を式（２）で表し、
さらに、固有振動数をｐ及び減衰比係数をξとすると、式（１）は次式（３）で表される。
したがって、固有振動数ｐ及び減衰比係数ξを設定すると、式（３）から鼓膜の振動応答ｘを求めることができる。この固有振動数ｐと減衰比係数ξが先に説明したモデルパラメータである。
なお、本実施の形態においては図４に示される機械的振動系を振動モデルとして説明したが、等価質量ｍに代えてインダクタンスＬ、ばね（ばね定数）ｋに代えて電気容量Ｃ、ダンパー（粘性減衰係数）ｃに代えて抵抗Ｒを用いる電気的振動系を振動モデルとして採用してもよい。
なお、本実施の形態においては図４に示される機械的振動系を振動モデルとして説明したが、等価質量ｍに代えてインダクタンスＬ、ばね（ばね定数）ｋに代えて電気容量Ｃ、ダンパー（粘性減衰係数）ｃに代えて抵抗Ｒを用いる電気的振動系を振動モデルとして採用してもよい。
次に、鼓膜の振動モデルを用いて振動応答ｘを求めた結果について図５を参照しながら説明する。
図５（ａ）は、鼓膜の振動モデルに基づいて正常心音から振動応答ｘを求めた結果を示すチャート図であり、図５（ｂ）は同じく鼓膜の振動モデルに基づいて僧帽弁閉鎖不全心音から振動応答ｘを求めた結果を示すチャート図である。
なお、心音データは市販の心音聴診トレーニング教材に収録されたものを使用し、また、固有振動数ｐは１０Ｈｚ、減衰比係数ξは０．７０７とした。
図５（ａ）及び図５（ｂ）において、薄い灰色で示される波形が各々の心音Ｓの原波形であり、濃く太い実線で示される波形が振動応答ｘによる波形である。
図５（ａ）では正常心音の「ドック、ドック」と表現される波形が、また、図５（ｂ）では僧帽弁閉鎖不全心音の「グー、グー」と表現される波形が良好に抽出されていることがわかる。
このように、鼓膜の振動モデルに基づくと、測定される心音データから振動応答ｘによって特徴的な波形を作成することが可能であり、この振動応答ｘによる波形が特徴値波形となる。
なお、図５（ａ）及び図５（ｂ）においては、正負の振幅を持つ特徴値波形を示しているが、解析はいずれか一方の特徴値波形について行うとよい。
図５（ａ）は、鼓膜の振動モデルに基づいて正常心音から振動応答ｘを求めた結果を示すチャート図であり、図５（ｂ）は同じく鼓膜の振動モデルに基づいて僧帽弁閉鎖不全心音から振動応答ｘを求めた結果を示すチャート図である。
なお、心音データは市販の心音聴診トレーニング教材に収録されたものを使用し、また、固有振動数ｐは１０Ｈｚ、減衰比係数ξは０．７０７とした。
図５（ａ）及び図５（ｂ）において、薄い灰色で示される波形が各々の心音Ｓの原波形であり、濃く太い実線で示される波形が振動応答ｘによる波形である。
図５（ａ）では正常心音の「ドック、ドック」と表現される波形が、また、図５（ｂ）では僧帽弁閉鎖不全心音の「グー、グー」と表現される波形が良好に抽出されていることがわかる。
このように、鼓膜の振動モデルに基づくと、測定される心音データから振動応答ｘによって特徴的な波形を作成することが可能であり、この振動応答ｘによる波形が特徴値波形となる。
なお、図５（ａ）及び図５（ｂ）においては、正負の振幅を持つ特徴値波形を示しているが、解析はいずれか一方の特徴値波形について行うとよい。
続いて、特徴値波形の解析方法について図６を用いて説明する。
一般に、正常な心音の波形ではＩ音及びII音と呼ばれるピークが繰り返し現れることが知られている。そして、このＩ音は前述の心臓の弁膜である僧帽弁及び三尖弁の閉鎖によって生じ、また、II音は大動脈弁及び肺動脈弁の閉鎖並びに緊張によって生じるものであり、本実施の形態では、心音解析の基本であるＩ音とII音の持続時間を解析並びに評価することによって心音の異常を検知することを考案した。
図６（ａ）は、本実施の形態に係るデジタル聴診解析システムを用いて作成した特徴値波形からの評価指数の求め方を示す説明図であり、図６（ｂ）は評価指数の関係を示すグラフ図であり、図６（ｃ）は評価指数と頻度の関係を示すグラフ図である。
なお、図６（ａ）において、縦軸は振幅を、横軸は時間を示している。実線で示される特徴値波形は見やすくするために測定時間における２秒間分を拡大したものであり、Ｉ音ピーク４８及びII音ピーク４９が各々２本ずつ抽出されている。なお、薄い灰色で示される波形は１０秒間分の心音の波形であり、表記される横軸とは対応していない。
図６（ａ）において、Ｉ音ピーク４８及びII音ピーク４９の幅からＩ音及びII音の持続時間を求め、これらをそれぞれ評価指数Ｔ１及びＴ２と定義する。また、Ｉ音とII音を合わせた持続時間をみると、心臓の弁膜の閉鎖不全等によって現れるＩ音とII音の間の雑音に関する情報が得られることから、Ｉ音ピーク４８の開始からII音ピーク４９の終了までの持続時間を評価指数Ｔ１２と定義する。そして、不整脈や心拍リズムを評価するために、Ｉ音ピーク４８の開始から次のＩ音ピーク４８の開始までの持続時間を評価指数Ｔ１１と定義する。
一般に、正常な心音の波形ではＩ音及びII音と呼ばれるピークが繰り返し現れることが知られている。そして、このＩ音は前述の心臓の弁膜である僧帽弁及び三尖弁の閉鎖によって生じ、また、II音は大動脈弁及び肺動脈弁の閉鎖並びに緊張によって生じるものであり、本実施の形態では、心音解析の基本であるＩ音とII音の持続時間を解析並びに評価することによって心音の異常を検知することを考案した。
図６（ａ）は、本実施の形態に係るデジタル聴診解析システムを用いて作成した特徴値波形からの評価指数の求め方を示す説明図であり、図６（ｂ）は評価指数の関係を示すグラフ図であり、図６（ｃ）は評価指数と頻度の関係を示すグラフ図である。
なお、図６（ａ）において、縦軸は振幅を、横軸は時間を示している。実線で示される特徴値波形は見やすくするために測定時間における２秒間分を拡大したものであり、Ｉ音ピーク４８及びII音ピーク４９が各々２本ずつ抽出されている。なお、薄い灰色で示される波形は１０秒間分の心音の波形であり、表記される横軸とは対応していない。
図６（ａ）において、Ｉ音ピーク４８及びII音ピーク４９の幅からＩ音及びII音の持続時間を求め、これらをそれぞれ評価指数Ｔ１及びＴ２と定義する。また、Ｉ音とII音を合わせた持続時間をみると、心臓の弁膜の閉鎖不全等によって現れるＩ音とII音の間の雑音に関する情報が得られることから、Ｉ音ピーク４８の開始からII音ピーク４９の終了までの持続時間を評価指数Ｔ１２と定義する。そして、不整脈や心拍リズムを評価するために、Ｉ音ピーク４８の開始から次のＩ音ピーク４８の開始までの持続時間を評価指数Ｔ１１と定義する。
なお、これらの評価指数は、全てを用いなくてもよく、４個の持続時間のうち１の持続時間のみ用いてもよいし、また、４個の持続時間を如何様に組合せて用いることもできる。さらに、評価指数は本実施の形態においてはＩ音とII音に着目しているが、これら以外の心音の要素に着目してもよい。例えば、Ｉ音の僧帽弁及び三尖弁の閉鎖音、あるいはII音の大動脈弁及び肺動脈弁の閉鎖音が分離される場合には、これらの閉鎖音の持続時間や弁間の閉鎖音の持続時間やそれらを適宜組合わせて評価指数としてもよい。
また、Ｉ音やII音以外の、III音やIV音を用いてもよい。III音とは、II音にひき続いて起こる低調の音であり、急速充満期、すなわち拡張期のはじめに心房から心室内へ血液が急速に流入する際に生じるものであり、IV音とは、拡張期の終わり近く、すなわち前収縮期に左心房が収縮するときに心室壁が伸展されて生ずる音であり、心房音とも呼ばれるものである。これらの音に着目してそれぞれの持続時間等を評価指数として採用してもよい。
また、Ｉ音やII音以外の、III音やIV音を用いてもよい。III音とは、II音にひき続いて起こる低調の音であり、急速充満期、すなわち拡張期のはじめに心房から心室内へ血液が急速に流入する際に生じるものであり、IV音とは、拡張期の終わり近く、すなわち前収縮期に左心房が収縮するときに心室壁が伸展されて生ずる音であり、心房音とも呼ばれるものである。これらの音に着目してそれぞれの持続時間等を評価指数として採用してもよい。
次に、図６（ｂ）において、縦軸は評価指数であるＴ２及びＴ１２を示し、横軸は評価指数であるＴ１とＴ１１を示している。
図６（ｂ）では評価指数Ｔ１に対して評価指数Ｔ２をプロットし、また、評価指数Ｔ１１に対して評価指数Ｔ１２をプロットして分布を示しており、視覚的に心音の異常を判断することができる。すなわち、図６（ｂ）において点線で囲まれる正常値範囲の領域に評価指数の分布が入っている場合は、心音は正常であると判断でき、逆に、正常値範囲の領域に評価指数の分布が入っていない場合は、心音は異常であると判断できる。なお、心音データには個人差があるので、正常値範囲は、多数の健常者からデータを取得して統計的に決定することが望ましい。
そして、図６（ｃ）において、縦軸は評価指数を示し、縦軸は各々の評価指数の頻度を示している。このように各々の評価指数の値に対して、測定時間である１０秒間における各々の評価指数の出現頻度を棒グラフで示すと、各々の評価指数の出現頻度を簡単に比較することができる。例えば、評価指数Ｔ１の出現回数よりも評価指数Ｔ２の出現回数が明らかに少ない場合は、II音ピーク４９が欠損して観測されない不整脈の可能性を指摘することができる。また、図６（ｃ）における評価指数Ｔ１１は複数の棒で表現されているように評価指数の値にばらつきがあることが視覚的に簡単に判断できる。
このように、本実施の形態においては、特徴値波形を作成して算出される評価指数を評価することによって専門的な知識がなくても簡単に心音の異常を検知することができるのである。
図６（ｂ）では評価指数Ｔ１に対して評価指数Ｔ２をプロットし、また、評価指数Ｔ１１に対して評価指数Ｔ１２をプロットして分布を示しており、視覚的に心音の異常を判断することができる。すなわち、図６（ｂ）において点線で囲まれる正常値範囲の領域に評価指数の分布が入っている場合は、心音は正常であると判断でき、逆に、正常値範囲の領域に評価指数の分布が入っていない場合は、心音は異常であると判断できる。なお、心音データには個人差があるので、正常値範囲は、多数の健常者からデータを取得して統計的に決定することが望ましい。
そして、図６（ｃ）において、縦軸は評価指数を示し、縦軸は各々の評価指数の頻度を示している。このように各々の評価指数の値に対して、測定時間である１０秒間における各々の評価指数の出現頻度を棒グラフで示すと、各々の評価指数の出現頻度を簡単に比較することができる。例えば、評価指数Ｔ１の出現回数よりも評価指数Ｔ２の出現回数が明らかに少ない場合は、II音ピーク４９が欠損して観測されない不整脈の可能性を指摘することができる。また、図６（ｃ）における評価指数Ｔ１１は複数の棒で表現されているように評価指数の値にばらつきがあることが視覚的に簡単に判断できる。
このように、本実施の形態においては、特徴値波形を作成して算出される評価指数を評価することによって専門的な知識がなくても簡単に心音の異常を検知することができるのである。
続いて、本実施の形態に係るデジタル聴診解析システムを使用して心疾患症の心音を解析した結果を図７乃至図１１を用いて説明する。
なお、心疾患症の心音データには心音聴診トレーニング用ＣＤ教材及びインターネット上のオンライン臨床実習教材に収録されたものを使用し、また、固有振動数ｐは１０Ｈｚ、減衰比係数ξを０．７とした。
図７（ａ）は、本実施の形態に係るデジタル聴診解析システムを用いて作成した僧帽弁閉鎖不全症の特徴値波形であり、図７（ｂ）は図７（ａ）から求めた評価指数の関係を示すグラフ図である。
図７（ａ）において、僧帽弁閉鎖不全症の特徴値波形では、持続時間の長いＩ音ピークと音ピークが繰り返し観察される。僧帽弁閉鎖不全症では、僧帽弁の閉鎖に異常があるために、僧帽弁の閉鎖に関わるＩ音が止まりにくくII音まで続くことが知られており、図７（ａ）における特徴値波形は僧帽弁閉鎖不全症の症状を明確に示しているといえる。
また、図７（ｂ）において、評価指数の分布をみると、Ｔ１とＴ２の関係及びＴ１１とＴ１２の関係はいずれも正常値範囲内に入っておらず、心音に異常があることが示唆される。また、Ｔ１は０．１から０．５の近傍にまで分布しており、この値は正常な心音のＴ１２（約０．３〜０．５）と同じ程度まで観測されていることより、Ｉ音が長くII音まで続いていることが容易に推察される。
なお、心疾患症の心音データには心音聴診トレーニング用ＣＤ教材及びインターネット上のオンライン臨床実習教材に収録されたものを使用し、また、固有振動数ｐは１０Ｈｚ、減衰比係数ξを０．７とした。
図７（ａ）は、本実施の形態に係るデジタル聴診解析システムを用いて作成した僧帽弁閉鎖不全症の特徴値波形であり、図７（ｂ）は図７（ａ）から求めた評価指数の関係を示すグラフ図である。
図７（ａ）において、僧帽弁閉鎖不全症の特徴値波形では、持続時間の長いＩ音ピークと音ピークが繰り返し観察される。僧帽弁閉鎖不全症では、僧帽弁の閉鎖に異常があるために、僧帽弁の閉鎖に関わるＩ音が止まりにくくII音まで続くことが知られており、図７（ａ）における特徴値波形は僧帽弁閉鎖不全症の症状を明確に示しているといえる。
また、図７（ｂ）において、評価指数の分布をみると、Ｔ１とＴ２の関係及びＴ１１とＴ１２の関係はいずれも正常値範囲内に入っておらず、心音に異常があることが示唆される。また、Ｔ１は０．１から０．５の近傍にまで分布しており、この値は正常な心音のＴ１２（約０．３〜０．５）と同じ程度まで観測されていることより、Ｉ音が長くII音まで続いていることが容易に推察される。
次に、図８（ａ）は、本実施の形態に係るデジタル聴診解析システムを用いて作成した不整脈の特徴値波形であり、図８（ｂ）は図８（ａ）から求めた評価指数の関係を示すグラフ図であり、図８（ｃ）は図８（ａ）から求めた評価指数と頻度の関係を示すグラフ図である。これらの図の縦軸と横軸は図６を参照して説明したとおりである。
図８（ａ）において、不整脈の特徴値波形では、明瞭なＩ音ピークとII音ピークが繰り返し観察されるが、その間隔が不規則になっていることがわかる。また、図８（ｂ）においては、Ｔ１とＴ２の関係は正常値範囲内に入っているものの、Ｔ１１とＴ１２の関係は正常値範囲から外れ、特にＴ１１の値が大きく振れている。さらに、図８（ｃ）においても、Ｔ１１の値にばらつきが大きいことがよくわかる。すなわち、Ｔ１１の数値に表される心拍のリズムに異常があると判断することができ、不整脈の症状を的確に捉えているといえる。
図８（ａ）において、不整脈の特徴値波形では、明瞭なＩ音ピークとII音ピークが繰り返し観察されるが、その間隔が不規則になっていることがわかる。また、図８（ｂ）においては、Ｔ１とＴ２の関係は正常値範囲内に入っているものの、Ｔ１１とＴ１２の関係は正常値範囲から外れ、特にＴ１１の値が大きく振れている。さらに、図８（ｃ）においても、Ｔ１１の値にばらつきが大きいことがよくわかる。すなわち、Ｔ１１の数値に表される心拍のリズムに異常があると判断することができ、不整脈の症状を的確に捉えているといえる。
続いて、図９（ａ）は、本実施の形態に係るデジタル聴診解析システムを用いて作成した重症の大動脈弁閉鎖不全症の特徴値波形であり、図９（ｂ）は図９（ａ）から求めた評価指数の関係を示すグラフ図である。
図９（ａ）において、重症の大動脈弁閉鎖不全症の特徴値波形では、大動脈弁の閉鎖及び緊張に関わるII音ピークの後に雑音に起因するピークが認められており、大動脈弁閉鎖不全症の症状を示していることがわかる。また、図９（ｂ）では、Ｔ１１とＴ１２の関係が正常値範囲から大きく外れており、心音に異常があることが明確である。
図９（ａ）において、重症の大動脈弁閉鎖不全症の特徴値波形では、大動脈弁の閉鎖及び緊張に関わるII音ピークの後に雑音に起因するピークが認められており、大動脈弁閉鎖不全症の症状を示していることがわかる。また、図９（ｂ）では、Ｔ１１とＴ１２の関係が正常値範囲から大きく外れており、心音に異常があることが明確である。
また、図１０（ａ）は、本実施の形態に係るデジタル聴診解析システムを用いて作成した軽症の大動脈弁閉鎖不全症の特徴値波形であり、図１０（ｂ）は図１０（ａ）から求めた評価指数の関係を示すグラフ図である。
図１０（ａ）において、軽症の大動脈弁閉鎖不全症の特徴値波形では、比較的明瞭なＩ音ピークとII音ピークが繰り返し観察され、図９（ａ）に比べると小さいが、II音ピークの後の雑音がわずかに認められる。また、図１０（ｂ）において、評価指数の分布をみると、Ｔ１とＴ２及びＴ１１及びＴ１２のいずれの関係についてもほぼ正常値範囲内に入っており、心音の異常を察知するのが困難であると考えられる。そこで、減衰比係数ξを０．３と小さくして特徴値波形の解析を試みた。
図１０（ａ）において、軽症の大動脈弁閉鎖不全症の特徴値波形では、比較的明瞭なＩ音ピークとII音ピークが繰り返し観察され、図９（ａ）に比べると小さいが、II音ピークの後の雑音がわずかに認められる。また、図１０（ｂ）において、評価指数の分布をみると、Ｔ１とＴ２及びＴ１１及びＴ１２のいずれの関係についてもほぼ正常値範囲内に入っており、心音の異常を察知するのが困難であると考えられる。そこで、減衰比係数ξを０．３と小さくして特徴値波形の解析を試みた。
図１１（ａ）は、図１０において減衰比係数を小さくして求めた特徴値波形であり、図１１（ｂ）は図１１（ａ）中の符号Ａで囲まれた部分の拡大図である。
図１１（ａ）において、減衰比係数ξを０．３と小さくして求めた特徴値波形では、図１０（ａ）に比べると、II音ピークの後に雑音によるピークが強調されて示されている。さらに、図１１（ｂ）の拡大図によると、雑音によるピークが明瞭に現れていることがわかる。
このように本実施の形態では、モデルパラメータである減衰比係数ξを変更することによって、見逃しやすい軽症の心疾患症についても正確に解析して診断することが可能であることがわかった。
図１１（ａ）において、減衰比係数ξを０．３と小さくして求めた特徴値波形では、図１０（ａ）に比べると、II音ピークの後に雑音によるピークが強調されて示されている。さらに、図１１（ｂ）の拡大図によると、雑音によるピークが明瞭に現れていることがわかる。
このように本実施の形態では、モデルパラメータである減衰比係数ξを変更することによって、見逃しやすい軽症の心疾患症についても正確に解析して診断することが可能であることがわかった。
次に、本実施の形態に係るデジタル聴診解析システムにおける解析方法の妥当性の検証について図１２乃至図１４を参照しながら説明する。
まず、固有振動数の影響について説明する。
図１２（ａ）は、本実施の形態に係るデジタル聴診解析システムを用いて固有振動数を２０Ｈｚの条件で作成した正常心音の特徴値波形であり、図１２（ｂ）は同じく固有振動数を５Ｈｚの条件で作成した正常心音の特徴値波形であり、図１２（ｃ）は図１２（ｂ）から求めた評価指数の関係を示すグラフ図であり、図１２（ｄ）は図１２（ａ）から求めた評価指数の関係を示すグラフ図である。
なお、減衰比係数ξは０．７０７とした。
図１２（ａ）及び図１２（ｂ）において、固有振動数を変更すると、固有振動数の大きい図１２（ａ）では、固有振動数の小さい図１２（ｂ）に比べて特徴値波形の微細な振動を捉えていることがわかる。したがって、固有振動数は、微細な振動を必要とする場合は大きくするとよく、聴診器の特性、録音状態、心音の個人差及び心疾患の種類の応じて任意に変更することによって解析に適した特徴値波形を得ることができる。
一方、図１２（ｃ）及び図１２（ｄ）において、評価指数の分布をみると、いずれの固有振動数の場合においても評価指数は正常値範囲内に入っており、また、両者の評価指数はほぼ同値を示していることより、評価指数は固有振動数に依存しないといえる。なお、評価指数における固有振動数の影響を少なくするために、特徴値波形の解析においてノイズを除去する工夫を施している。
まず、固有振動数の影響について説明する。
図１２（ａ）は、本実施の形態に係るデジタル聴診解析システムを用いて固有振動数を２０Ｈｚの条件で作成した正常心音の特徴値波形であり、図１２（ｂ）は同じく固有振動数を５Ｈｚの条件で作成した正常心音の特徴値波形であり、図１２（ｃ）は図１２（ｂ）から求めた評価指数の関係を示すグラフ図であり、図１２（ｄ）は図１２（ａ）から求めた評価指数の関係を示すグラフ図である。
なお、減衰比係数ξは０．７０７とした。
図１２（ａ）及び図１２（ｂ）において、固有振動数を変更すると、固有振動数の大きい図１２（ａ）では、固有振動数の小さい図１２（ｂ）に比べて特徴値波形の微細な振動を捉えていることがわかる。したがって、固有振動数は、微細な振動を必要とする場合は大きくするとよく、聴診器の特性、録音状態、心音の個人差及び心疾患の種類の応じて任意に変更することによって解析に適した特徴値波形を得ることができる。
一方、図１２（ｃ）及び図１２（ｄ）において、評価指数の分布をみると、いずれの固有振動数の場合においても評価指数は正常値範囲内に入っており、また、両者の評価指数はほぼ同値を示していることより、評価指数は固有振動数に依存しないといえる。なお、評価指数における固有振動数の影響を少なくするために、特徴値波形の解析においてノイズを除去する工夫を施している。
次に、実測したデータをデジタル聴診解析システムを使用して解析した結果について図１３及び図１４を用いて説明する。
図１３（ａ）は、本実施の形態に係るデジタル聴診解析システムを用いて作成した正常心音の特徴値波形の一例であり、（ｂ）は（ａ）から求めた評価指数の関係を示すグラフ図であり、（ｃ）は（ａ）から求めた評価指数と頻度の関係を示すグラフ図である。
なお、固有振動数ｐは１０Ｈｚとし、減衰比係数ξは０．７０７として解析を行った。また、被験者は２０代後半の男性の健常者である。
図１３（ａ）において、正常心音の特徴値波形では、明瞭なＩ音ピークとII音ピークが繰り返して現れており、また、図１３（ｂ）において、特徴値波形から算出された評価指数の分布をみると、Ｔ１とＴ２及びＴ１１とＴ１２の関係はいずれも正常値範囲内に入っていることがわかる。さらに、図１３（ｃ）では、Ｔ１、Ｔ２、Ｔ１１及びＴ１２はいずれもばらつきが少なく、また、ほぼ同頻度であることより、特徴値波形のＩ音ピーク及びII音ピークは規則的に出現していることが推察される。したがって、この心音は正常であると判断できる。
図１３（ａ）は、本実施の形態に係るデジタル聴診解析システムを用いて作成した正常心音の特徴値波形の一例であり、（ｂ）は（ａ）から求めた評価指数の関係を示すグラフ図であり、（ｃ）は（ａ）から求めた評価指数と頻度の関係を示すグラフ図である。
なお、固有振動数ｐは１０Ｈｚとし、減衰比係数ξは０．７０７として解析を行った。また、被験者は２０代後半の男性の健常者である。
図１３（ａ）において、正常心音の特徴値波形では、明瞭なＩ音ピークとII音ピークが繰り返して現れており、また、図１３（ｂ）において、特徴値波形から算出された評価指数の分布をみると、Ｔ１とＴ２及びＴ１１とＴ１２の関係はいずれも正常値範囲内に入っていることがわかる。さらに、図１３（ｃ）では、Ｔ１、Ｔ２、Ｔ１１及びＴ１２はいずれもばらつきが少なく、また、ほぼ同頻度であることより、特徴値波形のＩ音ピーク及びII音ピークは規則的に出現していることが推察される。したがって、この心音は正常であると判断できる。
また、図１４（ａ）は、本実施の形態に係るデジタル聴診解析システムを用いて作成した正常心音の特徴値波形の一例であり、（ｂ）は（ａ）から求めた評価指数の関係を示すグラフ図であり、（ｃ）は（ａ）から求めた評価指数と頻度の関係を示すグラフ図である。
なお、図１３の場合と同様に、固有振動数ｐは１０Ｈｚとし、減衰比係数ξは０．７０７として解析を行った。また、被験者は２０代後半の男性の健常者である。
図１４（ａ）では、図１３（ａ）と同様に、明瞭なＩ音ピークとII音ピークが繰り返した特徴値波形が得られている。
一方、図１３（ｂ）においては、Ｔ１とＴ２の関係は正常値範囲内に入っているものの、Ｔ１１とＴ１２の関係は正常値範囲内からわずかに外れている。Ｔ１１とＴ１２の関係をよくみると、Ｔ１２は正常値範囲内であるが、Ｔ１１が正常値範囲内から外れており、これは、被験者が緊張しやすい性格であり、測定時に緊張して心拍が早くなったことが影響していると考えられる。
また、図１３（ｃ）では、いずれの評価指数もばらつきが少なく、ほぼ同頻度であることより、被験者は、特に規則正しい心拍リズムを有していると考えられる。
このように、正常な心音であっても、緊張等により心拍数が早くなることがあるので、正常値範囲は統計的な要素の他に個人の特性を考慮しながら設定することが好ましい。
なお、図１３の場合と同様に、固有振動数ｐは１０Ｈｚとし、減衰比係数ξは０．７０７として解析を行った。また、被験者は２０代後半の男性の健常者である。
図１４（ａ）では、図１３（ａ）と同様に、明瞭なＩ音ピークとII音ピークが繰り返した特徴値波形が得られている。
一方、図１３（ｂ）においては、Ｔ１とＴ２の関係は正常値範囲内に入っているものの、Ｔ１１とＴ１２の関係は正常値範囲内からわずかに外れている。Ｔ１１とＴ１２の関係をよくみると、Ｔ１２は正常値範囲内であるが、Ｔ１１が正常値範囲内から外れており、これは、被験者が緊張しやすい性格であり、測定時に緊張して心拍が早くなったことが影響していると考えられる。
また、図１３（ｃ）では、いずれの評価指数もばらつきが少なく、ほぼ同頻度であることより、被験者は、特に規則正しい心拍リズムを有していると考えられる。
このように、正常な心音であっても、緊張等により心拍数が早くなることがあるので、正常値範囲は統計的な要素の他に個人の特性を考慮しながら設定することが好ましい。
次に、図１５を参照しながら、鼓膜の振動モデルとして鼓膜の等価質量を２分割してそれぞれｍ１，ｍ２として、質量、ばね、ダンパーから構成される振動系を直列に連成させたものを採用した場合について説明する。図１５は、２分割された鼓膜の振動モデルを示す概念図である。
図１５において、符号４５ａ，４５ｂは鼓膜、符号４６ａ，４６ｂはばね、そして、符号４７ａ，４７ｂはダンパーであり、鼓膜の等価質量をｍ１，ｍ２、ばね定数をｋ１，ｋ２、粘性減衰係数をｃ１，ｃ２及び聴診器からの心音をＳとして、その振動応答をそれぞれｘ１，ｘ２とするものである。
なお、本図においても機械的振動系を用いているが、前述のとおり電気的振動系を用いてもよい。
図１５において、符号４５ａ，４５ｂは鼓膜、符号４６ａ，４６ｂはばね、そして、符号４７ａ，４７ｂはダンパーであり、鼓膜の等価質量をｍ１，ｍ２、ばね定数をｋ１，ｋ２、粘性減衰係数をｃ１，ｃ２及び聴診器からの心音をＳとして、その振動応答をそれぞれｘ１，ｘ２とするものである。
なお、本図においても機械的振動系を用いているが、前述のとおり電気的振動系を用いてもよい。
このような鼓膜の振動モデルを用いて振動応答ｘ１，ｘ２を求めた結果について図１６乃至図１９を参照しながら説明する。
図１６は、（ａ）は２分割された鼓膜の振動モデルを採用した本実施の形態に係るデジタル聴診解析システムを用いて作成した僧帽弁閉鎖不全の弁膜症の特徴値波形のうち、図１５の振動モデルの振動応答ｘ１による特徴値波形を示す例であり、（ｂ）は（ａ）から求めた評価指数の関係を示すグラフ図であり、（ｃ）は（ａ）から求めた評価指数と頻度の関係を示すグラフ図である。
一方、図１７の（ａ）乃至（ｃ）は、それぞれ図１６の（ａ）乃至（ｃ）に対して、振動モデルの振動応答ｘ１に代えてｘ２を採用した場合のグラフ図となっている。
これらの図を比較すると、それぞれ（ｂ），（ｃ）では大差がないものの、図１６（ａ）では、図７を参照しながら説明したように、僧帽弁閉鎖不全症が僧帽弁の閉鎖の異常によるものであるために、僧帽弁の閉鎖に関わるＩ音が止まりにくくII音まで続く状態が示されている。一方、図１７（ａ）では、僧帽弁の閉鎖不全による異常心音がピーク波形として現れている。従って、デジタル聴診解析システムでは抽出が容易であり、より精度の高い診断に供することが可能である。なお、図１６（ａ），図１７（ａ）に薄い灰色で示されるのは心音Ｓの原波形である。また、この後の図１８（ａ），図１９（ａ）についても同様である。
図１６は、（ａ）は２分割された鼓膜の振動モデルを採用した本実施の形態に係るデジタル聴診解析システムを用いて作成した僧帽弁閉鎖不全の弁膜症の特徴値波形のうち、図１５の振動モデルの振動応答ｘ１による特徴値波形を示す例であり、（ｂ）は（ａ）から求めた評価指数の関係を示すグラフ図であり、（ｃ）は（ａ）から求めた評価指数と頻度の関係を示すグラフ図である。
一方、図１７の（ａ）乃至（ｃ）は、それぞれ図１６の（ａ）乃至（ｃ）に対して、振動モデルの振動応答ｘ１に代えてｘ２を採用した場合のグラフ図となっている。
これらの図を比較すると、それぞれ（ｂ），（ｃ）では大差がないものの、図１６（ａ）では、図７を参照しながら説明したように、僧帽弁閉鎖不全症が僧帽弁の閉鎖の異常によるものであるために、僧帽弁の閉鎖に関わるＩ音が止まりにくくII音まで続く状態が示されている。一方、図１７（ａ）では、僧帽弁の閉鎖不全による異常心音がピーク波形として現れている。従って、デジタル聴診解析システムでは抽出が容易であり、より精度の高い診断に供することが可能である。なお、図１６（ａ），図１７（ａ）に薄い灰色で示されるのは心音Ｓの原波形である。また、この後の図１８（ａ），図１９（ａ）についても同様である。
さらに、図１８と図１９は、図１６と図１７の関係と同様に、弁膜症においてIII音が聞こえる場合において、それぞれ振動モデルの振動応答ｘ１とｘ２を採用した場合のグラフ図を示すものである。
この図１８と図１９の場合においては、それぞれの（ｂ）には大差が生じていないものの、（ａ）では、図１６と図１７と同様に、弁膜症に起因する異常心音がピーク波形で現れていることに加えて、III音も明確にピーク波形として現れていることがわかる。
従って、III音が現れていても、弁膜症に基づく異常心音との識別が容易である。
さらに、図１８及び図１９の（ｃ）を比較すると、図１９（ｃ）の方が、より評価指数のばらつきが少なくなっていることがわかるが、このことによってもより精度の高い解析を行なうことが可能であることがわかる。
また、図１６，１７を用いて説明したとおり、振動応答ｘ２を用いることで、異常心音やIII音の波形がピーク波形として現れるので、デジタル聴診解析システムでの抽出も容易であり、より精度の高い診断に供することが可能である。
なお、本実施の形態においては、振動系を直列に２つ連成させた振動モデルを用いたが、さらに数を増やしたり、並列に連成させてもよい。
この図１８と図１９の場合においては、それぞれの（ｂ）には大差が生じていないものの、（ａ）では、図１６と図１７と同様に、弁膜症に起因する異常心音がピーク波形で現れていることに加えて、III音も明確にピーク波形として現れていることがわかる。
従って、III音が現れていても、弁膜症に基づく異常心音との識別が容易である。
さらに、図１８及び図１９の（ｃ）を比較すると、図１９（ｃ）の方が、より評価指数のばらつきが少なくなっていることがわかるが、このことによってもより精度の高い解析を行なうことが可能であることがわかる。
また、図１６，１７を用いて説明したとおり、振動応答ｘ２を用いることで、異常心音やIII音の波形がピーク波形として現れるので、デジタル聴診解析システムでの抽出も容易であり、より精度の高い診断に供することが可能である。
なお、本実施の形態においては、振動系を直列に２つ連成させた振動モデルを用いたが、さらに数を増やしたり、並列に連成させてもよい。
以上より、本実施の形態では、質量とばね要素と減衰振動要素を組合わせた振動系を１乃至複数連成させた振動モデルに適当な固有振動数や減衰比係数を選択することにより、種々の心音データの解析を精確に行い、心音の異常を簡単に察知することが可能である。
また、実測された心音データについても的確に評価することができ、本解析方法の妥当性が見出せた。
また、実測された心音データについても的確に評価することができ、本解析方法の妥当性が見出せた。
最後に、本実施の形態に係るデジタル聴診解析システムの解析の処理方法について図２０を用いて説明する。
図２０は、本実施の形態に係るデジタル聴診解析システムの解析方法を示すフローチャート図である。
図２０において、ステップＳ１は心音データの読込み工程を、ステップＳ２はモデルパラメータの読込み工程を示している。このステップＳ１及びステップＳ２では、演算部の特徴値波形解析部によって、解析入力データベースに格納されている心音データ及びモデルパラメータが読込まれる。なお、モデルパラメータは解析入力データベース以外にもモデルパラメータ入力部から読込むこともできる。また、ステップＳ１とステップＳ２の順序はいずれが先であっても構わない。
図２０は、本実施の形態に係るデジタル聴診解析システムの解析方法を示すフローチャート図である。
図２０において、ステップＳ１は心音データの読込み工程を、ステップＳ２はモデルパラメータの読込み工程を示している。このステップＳ１及びステップＳ２では、演算部の特徴値波形解析部によって、解析入力データベースに格納されている心音データ及びモデルパラメータが読込まれる。なお、モデルパラメータは解析入力データベース以外にもモデルパラメータ入力部から読込むこともできる。また、ステップＳ１とステップＳ２の順序はいずれが先であっても構わない。
次に、ステップＳ３は特徴値波形の作成工程を示している。このステップＳ３では、特徴値波形解析部において、鼓膜の振動モデルに基づいてステップＳ２で読込まれたモデルパラメータを使用してステップＳ１で読込まれた心音データから特徴値波形が作成される。そして、図示していないが、作成された特徴値波形はいったん解析結果データベースに格納される。
そして、ステップＳ４はＴ１の算出工程を示している。このステップＳ４では、評価指数解析部は解析結果データベースに格納された特徴値波形を読込み、Ｉ音ピークからＩ音の持続時間Ｔ１を算出する。そして、算出されたＴ１を解析結果データベースに格納する。
続いて、ステップＳ５ではＴ１が正常値範囲であるか否かを判断する。このステップＳ５では、異常心音評価部によって、解析結果データベースに格納されているＴ１と、解析入力データベースに格納されているＴ１の正常値範囲データが読込まれ、Ｔ１とＴ１の正常値範囲を比較する。そして、Ｔ１が正常値範囲内である場合は、ステップＳ７に進む。また、Ｔ１が正常値範囲から外れる場合は、ステップＳ６に進み、僧帽弁閉鎖不全症又は大動脈弁狭窄症の可能性を忠告する警告メッセージを生成し、表示部又は出力部に表示又は出力する。なお、警告メッセージを警告音や警告ランプに替えることも可能である。
続いて、ステップＳ５ではＴ１が正常値範囲であるか否かを判断する。このステップＳ５では、異常心音評価部によって、解析結果データベースに格納されているＴ１と、解析入力データベースに格納されているＴ１の正常値範囲データが読込まれ、Ｔ１とＴ１の正常値範囲を比較する。そして、Ｔ１が正常値範囲内である場合は、ステップＳ７に進む。また、Ｔ１が正常値範囲から外れる場合は、ステップＳ６に進み、僧帽弁閉鎖不全症又は大動脈弁狭窄症の可能性を忠告する警告メッセージを生成し、表示部又は出力部に表示又は出力する。なお、警告メッセージを警告音や警告ランプに替えることも可能である。
次に、ステップＳ７はＴ２の算出工程を示している。このステップＳ７では、ステップＳ４と同様に、評価指数解析部によって解析結果データベースに格納された特徴値波形が読込まれて、II音ピークからII音の持続時間Ｔ２を算出する。そして、算出されたＴ２を解析結果データベースに格納する。
また、ステップＳ８ではＴ２が正常値範囲であるか否かを判断する。このステップＳ８では、ステップＳ５と同様に、解析結果データベースに格納されているＴ２と、解析入力データベースに格納されているＴ２の正常値範囲データが異常心音評価部によって読込まれてＴ２とＴ２の正常値範囲が比較される。そして、Ｔ２が正常値範囲内である場合は、ステップＳ１０に進み、一方、Ｔ２が正常値範囲から外れる場合は、ステップＳ９に進み、僧帽弁狭窄症又は大動脈弁閉鎖不全症の可能性を忠告する警告メッセージ等を生成し、表示部等に表示する。
また、ステップＳ８ではＴ２が正常値範囲であるか否かを判断する。このステップＳ８では、ステップＳ５と同様に、解析結果データベースに格納されているＴ２と、解析入力データベースに格納されているＴ２の正常値範囲データが異常心音評価部によって読込まれてＴ２とＴ２の正常値範囲が比較される。そして、Ｔ２が正常値範囲内である場合は、ステップＳ１０に進み、一方、Ｔ２が正常値範囲から外れる場合は、ステップＳ９に進み、僧帽弁狭窄症又は大動脈弁閉鎖不全症の可能性を忠告する警告メッセージ等を生成し、表示部等に表示する。
続いて、ステップＳ１０はＴ１１の算出工程を示し、ステップＳ１２ではＴ１１が正常値範囲であるか否かを判断する。まず、ステップＳ１０では、評価指数解析部において解析結果データベースに格納された特徴値波形が読込まれ、Ｉ音ピークと次のＩ音ピークまでの持続時間Ｔ１１を算出し、算出したＴ１１を解析結果データベースに格納する。
そして、ステップＳ１１では、異常心音評価部において解析結果データベースに格納されているＴ１１と、解析入力データベースに格納されているＴ１１の正常値範囲データが読込まれ、Ｔ１１とＴ１１の正常値範囲が比較される。Ｔ１１が正常値範囲内である場合は、ステップＳ１３に進み、Ｔ１１が正常値範囲から外れる場合は、ステップＳ１２に進む。そして、ステップＳ１２では、不整脈の可能性を忠告する警告メッセージ等を生成し、表示部等に表示するのである。
そして、ステップＳ１１では、異常心音評価部において解析結果データベースに格納されているＴ１１と、解析入力データベースに格納されているＴ１１の正常値範囲データが読込まれ、Ｔ１１とＴ１１の正常値範囲が比較される。Ｔ１１が正常値範囲内である場合は、ステップＳ１３に進み、Ｔ１１が正常値範囲から外れる場合は、ステップＳ１２に進む。そして、ステップＳ１２では、不整脈の可能性を忠告する警告メッセージ等を生成し、表示部等に表示するのである。
さらに、ステップＳ１３はＴ１２の算出工程を示し、ステップＳ１４ではＴ１２が正常値範囲であるか否かを判断する。同様に、ステップＳ１３では、評価指数解析部によって解析結果データベースに格納された特徴値波形が読込まれ、Ｉ音ピークとII音ピークからＩ音とII音の持続時間Ｔ１２が算出されて解析結果データベースに格納される。
そして、ステップＳ１４では、異常心音評価部によって解析結果データベースに格納されているＴ１２と、解析入力データベースに格納されているＴ１２の正常値範囲データが読込まれてＴ１２とＴ１２の正常値範囲が比較される。
Ｔ１２が正常値範囲から外れる場合は、ステップＳ１５に進み、僧帽弁閉鎖不全症の可能性を忠告する警告メッセージ等を生成し、表示部等に表示する。
一方、Ｔ１２が正常値範囲内である場合は、最後のステップＳ１６に進み、心音は異常ないことを知らせるメッセージを作成して、表示部に表示して、解析を終了する。
そして、ステップＳ１４では、異常心音評価部によって解析結果データベースに格納されているＴ１２と、解析入力データベースに格納されているＴ１２の正常値範囲データが読込まれてＴ１２とＴ１２の正常値範囲が比較される。
Ｔ１２が正常値範囲から外れる場合は、ステップＳ１５に進み、僧帽弁閉鎖不全症の可能性を忠告する警告メッセージ等を生成し、表示部等に表示する。
一方、Ｔ１２が正常値範囲内である場合は、最後のステップＳ１６に進み、心音は異常ないことを知らせるメッセージを作成して、表示部に表示して、解析を終了する。
なお、本実施の形態においては、算出される持続時間Ｔ１、Ｔ２、Ｔ１１及びＴ１２の全てについて正常値範囲と比較して異常を検知する方法を示したが、全ての持続時間を評価対象としなくても、いずれか１の持続時間又はその組合せについて評価する方法にしてもよい。また、本実施の形態においては、Ｔ１を算出してＴ１に関する正常値範囲と比較して、その範囲から外れる場合に警告メッセージ等を生成し、次にＴ２を算出するというように各評価指数を順次計算する工程を示したが、最初に４つの評価指数をすべて算出し、それぞれの正常値範囲と比較して、症状の判断を行ってもよい。このようなアルゴリズムでは、すべての症状の可能性について同時に判断することが可能である。
また、本実施の形態では、たとえば、Ｔ１において正常値範囲に入らなければステップＳ６で僧帽弁閉鎖不全症又は大動脈弁狭窄症の可能性と判断され、その後の評価指数の算出は実施されないようなフローとなっているが、もちろん、たとえステップＳ５でＮｏとなり、ステップＳ６に進んでも再度そこからステップＳ７に進むようなアルゴリズムであってもよい。これはステップＳ５，６の他、ステップＳ８，９やステップＳ１１，１２においても同様である。このようなアルゴリズムにおいても、すべての症状の可能性について判断することが可能である。
なお、Ｔ１，Ｔ２，Ｔ１１，Ｔ１２に代えて、前述のとおり、心臓弁のいずれかの弁の閉鎖音に関して持続時間や、同一弁あるいは異なる弁の閉鎖音間持続時間、あるいはこれらを適宜組合わせたものを評価するようにしてもよい。
また、本実施の形態では、たとえば、Ｔ１において正常値範囲に入らなければステップＳ６で僧帽弁閉鎖不全症又は大動脈弁狭窄症の可能性と判断され、その後の評価指数の算出は実施されないようなフローとなっているが、もちろん、たとえステップＳ５でＮｏとなり、ステップＳ６に進んでも再度そこからステップＳ７に進むようなアルゴリズムであってもよい。これはステップＳ５，６の他、ステップＳ８，９やステップＳ１１，１２においても同様である。このようなアルゴリズムにおいても、すべての症状の可能性について判断することが可能である。
なお、Ｔ１，Ｔ２，Ｔ１１，Ｔ１２に代えて、前述のとおり、心臓弁のいずれかの弁の閉鎖音に関して持続時間や、同一弁あるいは異なる弁の閉鎖音間持続時間、あるいはこれらを適宜組合わせたものを評価するようにしてもよい。
このように構成された本実施の形態においては、専門知識を有さないと診断しにくい聴診器からの心音を鼓膜の振動モデルに基づいて特徴値波形を作成し、この特徴値波形に含まれる心臓弁の動作に関わる評価指数から種々の診断を行うので、一般家庭等において簡単に心疾患等の病気の有無を察知することができる。
以上説明したように、本発明の請求項１及び請求項２に記載された発明は、心音の異常を簡単に診断するデジタル聴診解析システムを提供可能であり、病院や診療所以外にも一般家庭において使用可能である。
１…無線聴診システム送信部 ２…マイク入力 ３…トリガー無線送信機 ４…プリアンプ ５…ノイズフィルタ ６…オーディオ無線送信機 ７…パワーアンプ ８…イヤホン／コンピュータ ９…ＩＣレコーダ １０…マイクロコントローラ １１…外部メモリ １２…無線聴診システム受信部 １３…オーディオ無線受信機 １４…トリガー無線受信機 １５…オートチューニングアンプ １６…イヤホン １７…コンピュータ １８…コンピュータ １９…デジタル聴診解析システム ２０…演算部 ２１…特徴値波形解析部 ２２…評価指数解析部 ２３…異常心音評価部 ２４…表示部 ２５…モデルパラメータ入力部 ２６…解析入力データベース ２７…心音データ ２８…モデルパラメータ ２９…正常値範囲データ ３０…解析結果データベース ３１…特徴値波形 ３２…Ｔ１データ ３３…Ｔ２データ ３４…Ｔ１１データ ３５…Ｔ１２データ ３６…心臓 ３７…左心房 ３８…左心室 ３９…僧帽弁 ４０…大動脈弁 ４１…右心房 ４２…右心室 ４３…三尖弁 ４４…肺動脈弁 ４５，４５ａ，４５ｂ…鼓膜 ４６，４６ａ，４６ｂ…ばね ４７，４７ａ，４７ｂ…ダンパー ４８…Ｉ音ピーク ４９…II音ピーク ｍ，ｍ１，ｍ２…鼓膜の等価質量 ｋ，ｋ１，ｋ２…ばね定数 ｃ，ｃ１，ｃ２…粘性減衰係数 Ｓ…心音 ｘ，ｘ１，ｘ２…鼓膜の振動応答 ｐ…固有振動数 ξ…減衰比係数
Claims (2)
- 解析入力情報格納部と、振動モデルを用いて特徴値波形を解析する心音波形解析部と、解析された前記特徴値波形を用いて評価指数を解析する評価指数解析部と、異常心音評価部と、解析結果情報格納部とを有するデジタル聴診解析システムであって、
前記解析入力情報格納部は、デジタル化された心音データと前記振動モデルの解析パラメータと前記評価指数の正常値範囲情報を読み出し可能に格納し、
前記心音波形解析部は、前記心音データと解析パラメータを読み出して前記振動モデルを用いて心音の特徴値波形を解析する手段と、この特徴値波形を前記解析結果情報格納部に格納する手段とを備え、
前記振動モデルは、質量にばね系及び粘性減衰系を接続させた１の機械的振動系又はインダクタンスに電気容量と抵抗を接続させた電気的振動系からなる、あるいはこれらの振動系を複数連成させてなる鼓膜の振動応答に関する振動モデルであり、前記解析入力情報格納部から読み出した前記解析パラメータ及び前記心音データを前記振動モデルに代入して前記鼓膜の振動応答を得て、これを特徴値波形とし、
前記評価指数は、前記特徴値波形から得られる心臓の僧帽弁、三尖弁、大動脈弁、肺動脈弁のいずれかの閉鎖音の持続時間、あるいはこれらの任意の弁に対して同一弁あるいは異なる２弁の閉鎖音間の持続時間、あるいはそれらの組合せとし、
前記評価指数解析部は、前記解析結果情報格納部から前記特徴値波形を読み出してこの特徴値波形から前記評価指数を演算する手段と、前記評価指数を前記解析結果情報格納部へ格納する手段とを備え、
前記異常心音評価部は、前記解析結果情報格納部から前記評価指数を読み出して前記解析入力情報格納部から読み出した前記評価指数の正常値範囲情報と比較する手段と、前記評価指数が前記正常値範囲内から外れる場合に警報情報信号を発する手段とを備えることを特徴とするデジタル聴診解析システム。 - 前記特徴値波形、評価指数、警報情報の少なくとも１を出力及び／又は表示する表示部を有することを特徴とする請求項１に記載のデジタル聴診解析システム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2005080720A JP4759727B2 (ja) | 2004-03-19 | 2005-03-18 | デジタル聴診解析システム |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2004081506 | 2004-03-19 | ||
| JP2004081506 | 2004-03-19 | ||
| JP2005080720A JP4759727B2 (ja) | 2004-03-19 | 2005-03-18 | デジタル聴診解析システム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2005296643A JP2005296643A (ja) | 2005-10-27 |
| JP4759727B2 true JP4759727B2 (ja) | 2011-08-31 |
Family
ID=35328866
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2005080720A Expired - Lifetime JP4759727B2 (ja) | 2004-03-19 | 2005-03-18 | デジタル聴診解析システム |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4759727B2 (ja) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007046504A1 (ja) * | 2005-10-21 | 2007-04-26 | Yamaguchi University | 聴診心音信号の処理方法及び聴診装置 |
| JP4904488B2 (ja) * | 2006-01-17 | 2012-03-28 | 国立大学法人 長崎大学 | 肺音診断装置 |
| US8920343B2 (en) | 2006-03-23 | 2014-12-30 | Michael Edward Sabatino | Apparatus for acquiring and processing of physiological auditory signals |
| WO2008000254A1 (en) * | 2006-06-26 | 2008-01-03 | Coloplast A/S | Multi parametric classification of cardiovascular sounds |
| JP2008113936A (ja) * | 2006-11-07 | 2008-05-22 | Yasuaki Nakagawa | 生体音聴診装置 |
| JP5207174B2 (ja) * | 2008-03-31 | 2013-06-12 | 国立大学法人山口大学 | 心音周波数解析装置 |
| AU2013287260B2 (en) * | 2012-07-05 | 2015-04-30 | Pulmonary Apps, Llc | Wireless stethoscope and method of use thereof |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5010889A (en) * | 1988-02-04 | 1991-04-30 | Bloodline Technology | Intelligent stethoscope |
| US7300405B2 (en) * | 2003-10-22 | 2007-11-27 | 3M Innovative Properties Company | Analysis of auscultatory sounds using single value decomposition |
| US20070055151A1 (en) * | 2005-01-20 | 2007-03-08 | Shertukde Hemchandra M | Apparatus and methods for acoustic diagnosis |
-
2005
- 2005-03-18 JP JP2005080720A patent/JP4759727B2/ja not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JP2005296643A (ja) | 2005-10-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Thiyagaraja et al. | A novel heart-mobile interface for detection and classification of heart sounds | |
| US20210145306A1 (en) | Managing respiratory conditions based on sounds of the respiratory system | |
| CN1870941B (zh) | 使用奇异值分解的听诊声音分析 | |
| US20170258343A1 (en) | Intrinsic frequency hemodynamic waveform analysis | |
| US20040260188A1 (en) | Automated auscultation system | |
| US20060245597A1 (en) | Detection of coronary artery disease using an electronic stethoscope | |
| US20240324950A1 (en) | Systems and Methods for Remote Patient Screening and Triage | |
| US20050157887A1 (en) | System for outputting acoustic signal from a stethoscope | |
| JP4848524B2 (ja) | 聴診心音信号の処理方法及び聴診装置 | |
| JP2008528124A (ja) | 音声認識を使用した聴診音の分析 | |
| CN108492890B (zh) | 一种人体健康状态监测系统及方法 | |
| JP7639827B2 (ja) | 肺音分析システム | |
| JP4759727B2 (ja) | デジタル聴診解析システム | |
| US20220151582A1 (en) | System and method for assessing pulmonary health | |
| JP7552705B2 (ja) | 肺音分析システム | |
| Oktivasari et al. | Active filter analysis on designing electronic stethoscope | |
| WO2022044129A1 (ja) | 肺音分析システム | |
| US7998083B2 (en) | Method and device for automatically determining heart valve damage | |
| JP7533589B2 (ja) | 肺音分析システム | |
| Choi et al. | Development of wireless heart sound acquisition system for screening heart valvular disorder | |
| JP7662932B2 (ja) | 聴診システム、聴診器、及び、方法 | |
| JP2024143139A (ja) | 聴診装置およびこれを備えた聴診システム、聴診方法、聴診プログラム | |
| SK | SMART WEARABLE PHONOCARDIOGRAM FOR REAL TIME HEART SOUND ANALYSIS AND PREDICTIVE CARDIAC HEALTHCARE. | |
| US20180092606A1 (en) | Heart Sound Processing Method and System for Detecting Cardiopathy | |
| Mantji et al. | Electronic stethoscope design, prototyping and testing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20071205 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20110215 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110411 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20110506 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 4759727 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| EXPY | Cancellation because of completion of term |