(57)【特許請求の範囲】
【請求項１】
睡眠中の生体が発する音声に起因する対象信号を取得する信号取得部と、
前記対象信号に基づいて、いびきの有無を含む前記音声の種類を推定する推定部と、
前記推定部の推定結果に基づいて、前記音声の種類を判定する種類判定部と、
前記種類判定部の判定結果に基づいて、睡眠中の前記生体の呼吸状態を判定する呼吸状態判定部と、
を備え、
前記呼吸状態判定部は、前記生体がいびきを発生していないと推定された第１区間の前後の第２区間および第３区間の前記音声の種類の判定結果に基づき、前記第１区間の前記生体の呼吸状態を判定し、
前記推定部は、前記音声の種類としていびきの種類を推定する呼吸状態判定装置。
【請求項２】
前記種類判定部は、前記生体が不安定いびきを発生しているか否かを判定する請求項１に記載の呼吸状態判定装置。
【請求項３】
前記呼吸状態判定部は、前記第１区間の前記生体の呼吸状態が閉塞性無呼吸状態か否かを判定する請求項２に記載の呼吸状態判定装置。
【請求項４】
前記呼吸状態判定部は、前記第１区間の前記生体の呼吸状態が無呼吸または低呼吸の状態であるか否かを判定する請求項２に記載の呼吸状態判定装置。
【請求項５】
前記呼吸状態判定部は、前記種類判定部の第２区間および第３区間の判定結果のうち、１以上の不安定いびきの発生に応じて、前記第１区間の呼吸状態を判定する請求項３または４に記載の呼吸状態判定装置。
【請求項６】
前記種類判定部の判定結果を蓄積するバッファ部を備え、
前記呼吸状態判定部は、前記バッファ部で蓄積した複数の前記判定結果に基づいて、複数の区間の前記生体の呼吸状態を判定する請求項５に記載の呼吸状態判定装置。
【請求項７】
前記呼吸状態判定部は、複数の前記判定結果が、予め定められた時間の間に予め定められた回数以上の前記不安定いびきが発生していることに応じて、警告信号を出力する請求項６に記載の呼吸状態判定装置。
【請求項８】
前記呼吸状態判定部の判定結果に基づいて、睡眠時無呼吸症候群の診断パラメータを算出する診断パラメータ算出部を更に備える請求項６または７に記載の呼吸状態判定装置。
【請求項９】
前記診断パラメータ算出部は、無呼吸低呼吸指数および睡眠障害指数の少なくとも１つを前記診断パラメータとして算出する請求項８に記載の呼吸状態判定装置。
【請求項１０】
前記診断パラメータ算出部は、前記種類判定部の第２区間および第３区間の判定結果に基づく予備パラメータを算出する請求項８または９に記載の呼吸状態判定装置。
【請求項１１】
前記対象信号の特徴量を算出する特徴量算出部と、
前記特徴量の非定常性を示す非定常性特徴量を算出する非定常性特徴量算出部と、
を備え、
前記推定部は、前記非定常性特徴量に基づき、前記音声の種類を推定する請求項１から１０のいずれか一項に記載の呼吸状態判定装置。
【請求項１２】
前記推定部は、前記非定常性特徴量に基づき、前記音声が複数の種類のいずれかに該当するかを示すスコアを生成するスコア生成部を有する請求項１１に記載の呼吸状態判定装置。
【請求項１３】
前記スコア生成部は、前記非定常性特徴量と予め生成した複数のモデルのそれぞれとの照合結果に基づいて、前記音声の種類を推定するための前記スコアを生成する請求項１２に記載の呼吸状態判定装置。
【請求項１４】
前記スコア生成部は、前記非定常性特徴量が予め定められた閾値を超えたか否かに応じて異なるスコア値を与える請求項１２に記載の呼吸状態判定装置。
【請求項１５】
前記特徴量算出部は、前記特徴量として前記音声の音高を示す第１特徴量、音の強さを示す第２特徴量、およびスペクトル波形を示す第３特徴量を算出する請求項１１から１４のいずれか一項に記載の呼吸状態判定装置。
【請求項１６】
前記非定常性特徴量算出部は、前記非定常性特徴量として前記第１特徴量の分散を示す第４特徴量、前記第２特徴量の分散を示す第５特徴量、および前記第３特徴量のエントロピーを示す第６特徴量のうち少なくとも１つを算出する請求項１５に記載の呼吸状態判定装置。
【請求項１７】
機器の作動方法であって、
睡眠中の生体が発する音声に起因する対象信号を取得する信号取得段階と、
前記対象信号に基づいて、いびきの有無を含む前記音声の種類を推定する推定段階と、
前記音声の種類の推定結果に基づいて、前記音声の種類を判定する種類判定段階と、
前記音声の種類の判定結果に基づいて、睡眠中の前記生体の呼吸状態を判定する呼吸状態判定段階と、
を備え、
前記呼吸状態判定段階は、前記生体がいびきを発生していないと推定された第１区間の前後の第２区間および第３区間の前記音声の種類の判定結果に基づき、前記第１区間の前記生体の呼吸状態を判定し、
前記推定段階は、前記音声の種類としていびきの種類を推定する作動方法。
【請求項１８】
コンピュータに、請求項１から１６のいずれか一項に記載の呼吸状態判定装置として機能させるプログラム。
【請求項１９】
請求項１８に記載のプログラムを記憶するコンピュータ読み取り可能な媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、呼吸状態判定装置、機器の作動方法、およびプログラムに関する。
【背景技術】
【０００２】
従来、睡眠時無呼吸症候群（ＳＡＳ：Ｓｌｅｅｐ  Ａｐｎｅａ  Ｓｙｎｄｒｏｍｅ）は、口および鼻の気流、血中酸素飽和濃度、心拍、胸および腹部の体動（呼吸パターン）等の生体信号に基づき、専門家が診断していた。また、睡眠中の生体が発するいびきに基づき、異常を検出する試みが知られていた（関連技術として、例えば、特許文献１および２を参照）。
特許文献１  特表２０１１−５２４７６９号公報
特許文献２  特開２００５−３０４９４１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００３】
しかしながら、専門家が実行するＳＡＳの診断は、特殊な装置を生体に装着して生体情報を取得するので、容易に診断することは困難であった。また、睡眠中の生体が発する音声に基づく診断は、睡眠中の雑音等をいびきとして検出してしまい、正確かつ簡便に診断することは困難であった。
【課題を解決するための手段】
【０００４】
本発明の第１の態様においては、睡眠中の生体が発する音声に起因する対象信号を取得する信号取得部と、対象信号に基づいて、いびきの有無を含む音声の種類を推定する推定部と、推定部の推定結果に基づいて、音声の種類を判定する種類判定部と、種類判定部の判定結果に基づいて、睡眠中の生体の呼吸状態を判定する呼吸状態判定部と、を備え、呼吸状態判定部は、生体がいびきを発生していないと推定された第１区間の前後の第２区間および第３区間の音声の種類の判定結果に基づき、第１区間の生体の呼吸状態を判定する呼吸状態判定装置、呼吸状態判定方法、およびプログラムを提供する。
【０００５】
本発明の第２の態様においては、睡眠中の生体が発する音声に起因する対象信号を取得する信号取得部と、対象信号に基づいて、睡眠時無呼吸症候群の診断パラメータを出力する診断パラメータ出力部と、を備える呼吸状態判定装置を提供する。
【０００６】
なお、上記の発明の概要は、本発明の必要な特徴の全てを列挙したものではない。また、これらの特徴群のサブコンビネーションもまた、発明となりうる。
【図面の簡単な説明】
【０００７】
【図１】本実施形態に係る呼吸状態判定装置１０の構成例を示す。
【図２】本実施形態に係る呼吸状態判定装置１０の動作フローを示す。
【図３】本実施形態に係る信号取得部１１０が取得する対象信号の一例を示す。
【図４】本実施形態に係る特徴量算出部１２０が算出した第１特徴量の一例を示す。
【図５】本実施形態に係る特徴量算出部１２０が算出した第２特徴量の一例を示す。
【図６】本実施形態に係る特徴量算出部１２０が算出した第４特徴量の一例を示す。
【図７】本実施形態に係る非定常性特徴量算出部１３０が算出した第５特徴量の一例を示す。
【図８】本実施形態に係る非定常性特徴量算出部１３０が算出した、安定いびきの領域の第６特徴量の一例を示す。
【図９】本実施形態に係る非定常性特徴量算出部１３０が算出した、不安定いびきの領域の第６特徴量の一例を示す。
【図１０】本実施形態に係る推定部１４０がいびきの有無を判断した結果の一例を示す。
【図１１】本実施形態に係るスコア生成部１４２が等間隔に並ぶスコアの値を各データに与えた結果を示す。
【図１２】本実施形態に係る種類判定部１５０が図１１に示す対象信号のスコアを平均化した結果の一例を示す。
【図１３】本実施形態に係る種類判定部１５０が図３に示す対象信号に対していびきの種類を判定した結果の一例を示す。
【図１４】本実施形態に係る呼吸状態判定装置１０の変形例を示す。
【図１５】本実施形態に係る呼吸状態判定装置１０として機能するコンピュータ１９００のハードウェア構成の一例を示す。
【発明を実施するための形態】
【０００８】
以下、発明の実施の形態を通じて本発明を説明するが、以下の実施形態は特許請求の範囲にかかる発明を限定するものではない。また、実施形態の中で説明されている特徴の組み合わせの全てが発明の解決手段に必須であるとは限らない。
【０００９】
図１は、本実施形態に係る呼吸状態判定装置１０の構成例を示す。呼吸状態判定装置１０は、睡眠中の生体が発する音声に起因する対象信号に基づき、口および鼻の気流、血中酸素飽和濃度、心拍、胸および腹部の動き等の生体信号を用いずに、音声の種類を判定する。また、呼吸状態判定装置１０は、音声の種類の判定結果に基づき、生体の呼吸状態を判定する。なお、本実施形態において、睡眠中の人間が発する音声を、いびき、寝息等の睡眠音とする。また、呼吸状態判定装置１０は、当該音声に起因する対象信号に基づき、いびきの種類を推定する例を説明する。
【００１０】
呼吸状態判定装置１０は、推定装置１００と、種類判定部１５０と、バッファ部１６０と、呼吸状態判定部１７０とを備える。推定装置１００は、対象信号に基づき、生体の音声の種類を推定する。ここで、対象信号は、睡眠中の生体が発する音声に応じた信号であればよく、例えば、音声信号、生体を伝わる振動信号、予め定められた体の部位の動きの検出信号、鼻孔または口からの気流の検出信号等である。推定装置１００は、信号取得部１１０と、特徴量算出部１２０と、非定常性特徴量算出部１３０と、推定部１４０とを有する。
【００１１】
信号取得部１１０は、睡眠中の生体が発する音声に起因する対象信号を取得する。信号取得部１１０は、データベース等の外部の記憶装置に接続され、予め定められた形式で記憶された対象信号を読み出して取得してよい。また、信号取得部１１０は、ネットワーク等に接続され、当該ネットワークを介して対象信号を取得してもよい。これに代えて、信号取得部１１０は、マイク、ＩＣレコーダ、携帯電話、タブレット型ＰＣ、パソコン、動き検出装置等の外部の装置に接続され、対象信号を取得してもよい。
【００１２】
また、信号取得部１１０は、外部の装置等に内蔵され、対象信号を取得してもよい。この場合、推定装置１００または呼吸状態判定装置１０全体が外部の装置等に内蔵されてもよい。これに代えて、信号取得部１１０は、マイク、動き検出装置等の検出装置を含み、当該検出装置の出力をデジタル変換して対象信号を取得してもよい。この場合、検出装置は、生体に対して非接触で音声、生体の動き、振動等を検出するセンサであることが望ましい。本実施形態において、信号取得部１１０が、睡眠中の生体の音声信号を取得する例を説明する。
【００１３】
特徴量算出部１２０は、信号取得部１１０に接続され、取得した対象信号の特徴量を算出する。特徴量算出部１２０は、例えば、対象信号の周波数領域におけるピーク周波数を特徴量とする。また、特徴量算出部１２０は、対象信号が周波数領域で複数のピークを有し、当該複数のピークが基本周波数成分と基本周波数の高調波成分を含む場合（即ち、調波構造を有する場合）、当該基本周波数の周波数を特徴量としてよい。特徴量算出部１２０は、一例として、対象信号が音声信号の場合、特徴量としてこのようなピーク周波数および／または基本周波数の周波数を音声の音高とし、当該音高を第１特徴量として算出する。
【００１４】
また、特徴量算出部１２０は、ピーク周波数および／または基本周波数の振幅値を特徴量としてよい。特徴量算出部１２０は、一例として、対象信号が音声信号の場合、このような振幅値を音の強さとし、当該音の強さを第２特徴量として算出する。また、特徴量算出部１２０は、対象信号を周波数に変換した周波数スペクトルを特徴量としてもよく、一例として、このようなスペクトル波形を第３特徴量として算出する。
【００１５】
非定常性特徴量算出部１３０は、特徴量算出部１２０に接続され、特徴量の非定常性を示す非定常性特徴量を算出する。非定常性特徴量算出部１３０は、例えば、特徴量の分散、標準偏差、エントロピー、最大値および最小値の範囲、四分位数等の、データのバラツキの度合いを非定常性特徴量として算出する。非定常性特徴量算出部１３０は、一例として、非定常性特徴量として第１特徴量の分散を示す第４特徴量、第２特徴量の分散を示す第５特徴量、および第３特徴量のエントロピーを示す第６特徴量のうち少なくとも１つを算出する。
【００１６】
推定部１４０は、対象信号に基づいて、いびきの有無を含む音声の種類を推定する。即ち、推定部１４０は、非定常性特徴量算出部１３０に接続され、非定常性特徴量に基づき、音声の種類を推定する。本実施形態において、推定部１４０は、音声の種類としていびきの種類を推定する。推定部１４０は、一例として、いびきの種類を「いびき無し（いびきをしていない状態）」、「安定いびき」、および「不安定いびき」として推定する。
【００１７】
ここで、安定いびきは、健康上、問題の無いいびきであり、不安定いびきは、生体が無呼吸または低呼吸の状態において発するいびきである。より具体的に説明すると、安定いびきは、気道が比較的軽く締まり、音の強さが安定し、呼吸停止が無く、寝息が混じり、基本周波数および高調波成分との分離が容易な周波数特性となる等の特徴を有する。なお、「基本周波数および高調波成分との分離が容易な周波数特性」とは、周波数スペクトルに特定の周波数帯域の振幅強度が強くなって有色となるといった（例えば、「調波構造を有する」）特性である。
【００１８】
これに対し、不安定いびきは、気道が強く締まり、不規則な音の強さで、呼吸停止を伴い、基本周波数および高調波成分との分離が比較的困難となる周波数特性となる等の特徴を有する。なお、「基本周波数および高調波成分との分離が比較的困難となる周波数特性」とは、周波数スペクトルの振幅強度が略一様に広がって白色性が強くなった特性である。
【００１９】
推定部１４０は、推定結果をスコア値として出力してよい。推定部１４０は、スコア生成部１４２および記憶部１４４を含む。
【００２０】
スコア生成部１４２は、非定常性特徴量に基づき、音声が複数の種類のいずれかに該当するかを示すスコアを生成する。スコア生成部１４２は、例えば、「いびき無し」と推定する場合、より低いスコア値を出力し、「安定いびき」と推定する場合、「いびき無し」と推定するスコア値よりも高いスコア値を出力し、「不安定いびき」と推定する場合、「安定いびき」と推定するスコア値よりも高いスコア値を出力する。スコア生成部１４２のスコア生成については後述する。
【００２１】
記憶部１４４は、スコア生成部１４２に接続され、スコア生成部１４２がスコア生成に用いるデータ、閾値、モデル等を記憶する。記憶部１４４は、スコア生成部１４２の要求に応じて、当該データ等を供給する。
【００２２】
また、記憶部１４４は、推定装置１００または呼吸状態判定装置１０の一時記憶領域として機能してもよい。記憶部１４４は、例えば、信号取得部１１０に接続され、当該信号取得部１１０が取得した対象信号を記憶する。また、記憶部１４４は、特徴量算出部１２０が算出する特徴量および非定常性特徴量算出部１３０が算出する非定常性特徴量を記憶してもよい。また、記憶部１４４は、推定装置１００または呼吸状態判定装置１０の動作の過程で生成する（または利用する）中間データ、算出結果、およびパラメータ等をそれぞれ記憶してもよい。また、記憶部１４４は、呼吸状態判定装置１０内の各部の要求に応じて、記憶したデータを要求元に供給してよい。
【００２３】
種類判定部１５０は、推定部１４０に接続され、推定部１４０の推定結果に基づき音声の種類を判定する。種類判定部１５０は、一例として、推定部１４０から受けとるスコア値に基づき、いびきの種類を判定する。なお、推定装置１００および種類判定部１５０は、対象信号に基づき、睡眠中の生体の音声の種類を判定する種類判定装置を構成してよい。
【００２４】
バッファ部１６０は、種類判定部１５０に接続され、種類判定部１５０の判定結果を蓄積する。一例として、種類判定部１５０は、時系列に判定した判定結果をバッファ部１６０に供給し、バッファ部１６０は受けとった複数の判定結果を順次記憶する。また、バッファ部１６０は受けとった複数の判定結果のうち、先に受けとった判定結果をさきに出力するＦＩＦＯ（Ｆｉｒｓｔ  Ｉｎ  Ｆｉｒｓｔ  ＯＵＴ）の機能を有してよい。
【００２５】
呼吸状態判定部１７０は、バッファ部１６０に接続され、種類判定部１５０の判定結果に基づいて、睡眠中の生体の呼吸状態を判定する。呼吸状態判定部１７０は、バッファ部１６０で蓄積した複数の区間の判定結果に基づいて、複数の区間の生体の呼吸状態を判定してよい。呼吸状態判定部１７０は、例えば、閉塞性無呼吸状態か否かを判定する。また、呼吸状態判定部１７０は、生体が無呼吸または低呼吸の状態であるか否かを判定してもよい。呼吸状態判定部１７０は、判定結果を、呼吸状態判定装置１０の出力とする。
【００２６】
以上の本実施形態に係る呼吸状態判定装置１０は、対象信号に基づき、睡眠中の生体のいびきの種類を判定し、当該判定結果に基づいて睡眠中の生体の呼吸状態を判定する。呼吸状態判定装置１０によるいびきの種類の判定および呼吸状態の判定について、図２を用いて説明する。図２は、本実施形態に係る呼吸状態判定装置１０の動作フローを示す。
【００２７】
まず、信号取得部１１０は、対象信号を取得する（Ｓ２００）。信号取得部１１０は、対象信号として、音声信号の時系列データを取得してよい。図３は、本実施形態に係る信号取得部１１０が取得する対象信号の一例を示す。図３の横軸は時間、縦軸は振幅を示す。図３は、０から略２００秒程度までの範囲で不安定いびきが発生し、略２５０秒から略６００秒までの範囲で安定いびきが発生した例を示す。
【００２８】
次に、特徴量算出部１２０は、特徴量を算出する（Ｓ２１０）。特徴量算出部１２０は、例えば、第１特徴量として基本周波数を算出する。特徴量算出部１２０は、対象信号にフーリエ変換等の周波数変換を施し、変換した周波数領域から複数のピーク値を取得して基本周波数を選択してよい。図４は、本実施形態に係る特徴量算出部１２０が算出した第１特徴量の一例を示す。図４の横軸は時間、縦軸は音高（基本周波数の周波数）を示す。安定いびきの領域の音高は、略一定の周波数値を示すが、不安定いびきの領域の音高は、安定いびきの音高に比べてバラツキが大きいことがわかる。
【００２９】
また、特徴量算出部１２０は、例えば、第２特徴量として基本周波数の振幅を算出する。特徴量算出部１２０は、対象信号を周波数変換したデータから、当該振幅値を取得してよい。図５は、本実施形態に係る特徴量算出部１２０が算出した第２特徴量の一例を示す。図５の横軸は時間、縦軸は音の強さ（基本周波数の振幅値）を示す。不安定いびきの領域の音の強さは、安定いびきの領域の音の強さと比較して、バラツキが大きいことがわかる。
【００３０】
また、特徴量算出部１２０は、例えば、第３特徴量として対象信号を周波数に変換した周波数スペクトルを算出する。特徴量算出部１２０は、第１特徴量および第２特徴量を算出する場合に算出した周波数変換の結果を、第３特徴量としてよい。
【００３１】
次に、非定常性特徴量算出部１３０は、非定常性特徴量を算出する（Ｓ２２０）。非定常性特徴量算出部１３０は、予め定められた時間間隔における非定常性特徴量を算出してよい。非定常性特徴量算出部１３０は、例えば、第１特徴量の分散を第４特徴量として算出する。図６は、本実施形態に係る特徴量算出部１２０が算出した第４特徴量の一例を示す。図６の横軸は時間、縦軸は音高の分散を示す。安定いびきの領域の分散は、略一定の周波数値を反映して分散値が小さくなるが、不安定いびきの領域の音高は、安定いびきの音高に比べてバラツキが大きいことがわかる。
【００３２】
また、非定常性特徴量算出部１３０は、例えば、第２特徴量の分散を第５特徴量として算出する。図７は、本実施形態に係る非定常性特徴量算出部１３０が算出した第５特徴量の一例を示す。図７の横軸は時間、縦軸は音の強さの分散を示す。不安定いびきの領域の音の強さの分散は、安定いびきの領域の分散と比較して、バラツキが大きいことがわかる。
【００３３】
また、非定常性特徴量算出部１３０は、例えば、第３特徴量のスペクトラルエントロピーを第６特徴量として算出する。ここで、エントロピーは、調波構造を有する（即ち、なんらかの秩序がある）スペクトルのエントロピーの値は小さくなり、背景雑音、突発的な音等はエントロピーの値が大きくなるので、バラツキ（乱雑さ）の度合いを示す特徴量として用いることが望ましい。非定常性特徴量算出部１３０は、一例として、次式で示されるスペクトラルエントロピーを算出して、第６特徴量としてよい。
【数１】
0006535186.tif 000002
【００３４】
ここで、ｆは周波数を示し、Ｓｆは周波数ｆにおける振幅強度を示す。また、（数１）式は、エントロピーを算出する式に（−１）を掛けて符合を反転させている。即ち、本実施形態のスペクトラルエントロピーは、調波構造を有する音の場合に値が大きくなり、雑音等の場合に値が小さくなる特徴量である。
【００３５】
図８は、本実施形態に係る非定常性特徴量算出部１３０が算出した、安定いびきの領域の第６特徴量の一例を示す。また、図９は、本実施形態に係る非定常性特徴量算出部１３０が算出した、不安定いびきの領域の第６特徴量の一例を示す。図８および図９の横軸は時間、縦軸は（数１）式のスペクトラルエントロピーを示す。安定いびきの領域のスペクトラルエントロピーは、略一定の大きさのピークが略等間隔に出現しているのに対して、不安定いびきの領域のスペクトラルエントロピーは、突発的なピークが生じていることがわかる。また、安定いびきの領域のスペクトラルエントロピーは、不安定いびきの領域のスペクトラルエントロピーと比較して、全体的に値が大きくなっていることがわかる。
【００３６】
次に、推定部１４０は、非定常特徴量に基づき、いびきの種類を推定する。本実施例において、推定部１４０のスコア生成部１４２は、いびきの種類を推定するスコアを生成する。まず、推定部１４０は、スコア生成部１４２がスコアを生成すべく、第６特徴量に基づき、生体がいびきを発生しているか否かを判別する（Ｓ２３０）。推定部１４０は、例えば、図８および図９に示したようなスペクトラルエントロピーに対して、予め定められた閾値を用いて、いびきの有無を判別する。推定部１４０は、一例として、当該閾値とスペクトラルエントロピーとを比較し、当該閾値を超えるデータが出現する領域を、生体がいびきを発生している期間と判断する。
【００３７】
生体が安定いびきを発生している場合、対象信号の音高および音の強さは安定しているものの、振幅値の値そのものは小さいので、背景雑音、寝息、寝返り等のノイズ成分との分離が困難である。しかしながら、本実施形態の推定部１４０は、対象信号をスペクトラルエントロピーに変換してからいびきの有無を判断するので、当該ノイズ成分と容易に分離することができ、いびきの有無の判断を正確に実行することができる。
【００３８】
また、推定部１４０は、図３に示すような対象信号の時間波形に対して、予め定められた閾値を用いて、いびきの有無を判別してもよい。生体が不安定いびきを発生している場合、対象信号の音高および音の強さは不安定ではあるが、振幅値の値そのものは大きく変動するので、不安定いびきの領域は時間波形の段階においても閾値を用いることで容易に判断することができる。推定部１４０は、スペクトラルエントロピーによるいびきの有無の判断に加えて、時間波形に対するいびきの有無の判断を実行してもよい。この場合、推定部１４０は、時間波形の振幅強度の大きさに応じていずれかを選択的に用いてよく、これに代えて、両者の判断を用いてもよい（重みを付けて判断してもよい）。
【００３９】
図１０は、本実施形態に係る推定部１４０がいびきの有無を判断した結果の一例を示す。図１０は、図３に示した対象信号の時間波形に対して、いびきの有無を判断した結果であり、Ｘ軸と実線とで囲む領域がいびきと判断したいびき領域である。例えば、図１０における区間Ａ、Ｃで示す領域はいびき領域である。また、区間Ｂ、Ｄで示す領域はいびき無しの領域である。また、区間Ｅで示す領域は、いびき領域といびき無しの領域が混在していることを示す。
【００４０】
図１０に示すように、推定部１４０は、区間Ａ、Ｃで示されるような振幅値の小さい安定いびきの領域においても、精度よくいびき領域を判断することができる。また、推定部１４０は、区間Ｅで示されるような複数の領域が混在する複雑な領域も精度よく判断できていることがわかる。
【００４１】
次に、スコア生成部１４２は、推定部１４０がいびきの有無を判断した結果と、非定常性特徴量に基づき、対象信号に含まれるデータのそれぞれにスコアを与える（Ｓ２４０からＳ２８０）。スコア生成部１４２は、「いびき無し」と判断したデータに対して、（Ｓ２４０：Ｎｏ）、対応する第１スコアを与える（Ｓ２４２）。スコア生成部１４２は、例えば、区間Ｂ、Ｄ、およびＥの一部で示す領域のデータに対して、第１スコアを与える。スコア生成部１４２は、一例として、第１スコアの値として０を与える。
【００４２】
また、スコア生成部１４２は、「いびきあり」と判断した領域（例えば、区間Ａ、Ｃ、およびＥの残りの一部）のデータに対しては（Ｓ２４０：Ｙｅｓ）、音高の分散の値が大きいか否かを更に判断する（Ｓ２５０）。スコア生成部１４２は、例えば、第４特徴量が予め定められた閾値を超えたか否かに応じて、音高の分散が大きいか否かを判断する。そして、スコア生成部１４２は、第４特徴量が閾値以下のデータに対して（Ｓ２５０：Ｎｏ）、対応する第２スコアを与える（Ｓ２５２）。スコア生成部１４２は、一例として、第１スコア値よりも大きいスコア値を第２スコア値として与える。
【００４３】
また、スコア生成部１４２は、第４特徴量が閾値を超えたデータに対して（Ｓ２５０：Ｙｅｓ）、音の強さの分散の値が大きいか否かを更に判断する（Ｓ２６０）。スコア生成部１４２は、例えば、第５特徴量が予め定められた閾値を超えたか否かに応じて、音の強さの分散が大きいか否かを判断する。そして、スコア生成部１４２は、第５特徴量が閾値以下のデータに対して（Ｓ２６０：Ｎｏ）、対応する第３スコアを与える（Ｓ２６２）。スコア生成部１４２は、一例として、第２スコア値よりも大きいスコア値を第３スコア値として与える。
【００４４】
また、スコア生成部１４２は、第５特徴量が閾値を超えたデータに対して（Ｓ２６０：Ｙｅｓ）、対応する第４スコアを与える（Ｓ２７０）。スコア生成部１４２は、一例として、第３スコア値よりも大きいスコア値を第４スコア値として与える。
【００４５】
このように、スコア生成部１４２は、「いびきあり」と判断した領域のデータに対して、非定常性特徴量が予め定められた閾値を超えたか否かに応じて異なるスコア値を与える。スコア生成部１４２は、対象信号の時間波形データに対して、対応するスコアをそれぞれ与えるまで、Ｓ２４０からＳ２８０のステップを繰り返してよい（Ｓ２８０：Ｎｏ）。この場合、スコア生成部１４２は、「いびきあり」と判断した領域のデータに対して、時系列に並ぶ順にそれぞれスコアを与えてよい。以上のように、スコア生成部１４２は、音高の分散および音の強さの分散が大きいと判断したデータに対してより大きいスコア値を与えるので、スコアの高いデータほど不安定いびきと推定することができる。
【００４６】
なお、各スコア値は、記憶部１４４に予め記憶され、スコア生成部１４２がそれぞれ読みだして用いてよい。また、記憶部１４４は、スコア生成部１４２が与えて生成した対象信号のスコアを記憶してもよい。ここで、第１スコア値から第４スコア値は、０、１、２、３と等間隔に並ぶ値としてよく、これに代えて、スコアの間隔に重みをつけて不当間隔に並ぶ値としてもよい。
【００４７】
図１１は、本実施形態に係るスコア生成部１４２が生成した対象信号のスコアの一例を示す。図１１は、スコア生成部１４２が図３および１０に示す対象信号に対して生成したスコアの一例であり、横軸が時間、縦軸がスコア値である。図１１は、スコア生成部１４２が等間隔に並ぶ第１スコア値から第４スコア値のうちいずれか一つの値を各データにそれぞれ与えた結果を示す。
【００４８】
スコア生成部１４２は、Ｓ２４０からＳ２８０の動作を繰り返して、時系列に並ぶいびきの有無が判断されたデータに対して、データ毎に順次スコアを与えてよい。スコア生成部１４２のスコア生成が終了すると（Ｓ２８０：Ｙｅｓ）、種類判定部１５０は、対象信号のスコアに応じて、いびきの種類を判定する（Ｓ２９０）。種類判定部１５０は、例えば、対象信号のスコアを平均化して、いびきの種類を判定する。種類判定部１５０は、時間軸において０ではない（非ゼロの）スコア値が連続する区間の平均を算出してよい。例えば、種類判定部１５０は、図１０に示す区間ＡおよびＣに示す領域に対応するスコアの区間について、それぞれ平均を算出する。
【００４９】
また、種類判定部１５０は、２つの非ゼロのスコア値の領域に挟まれるスコア値ゼロの領域の時間間隔が、予め定められた時間間隔以下の場合、２つの非ゼロのスコア値の領域およびゼロのスコア値の領域を同一区間とし、当該同一区間の平均を計算してよい。即ち、種類判定部１５０は、短期間で非ゼロのスコア値とゼロのスコア値が交互に繰り返される領域を、一つの区間とみなして平均化する。これにより、種類判定部１５０は、図１０に示す区間Ｅ示す領域に対応するスコアの区間について、一つの区間とみなして平均を算出する。これに代えて、種類判定部１５０は、移動平均を算出してもよい。図１２は、本実施形態に係る種類判定部１５０が図１１に示す対象信号のスコアを平均化した結果の一例を示す。なお、図１２において、予め定められた閾値の一例を点線で示す。
【００５０】
そして、種類判定部１５０は、対象信号のスコアを平均化した後に、スコアの平均値が予め定められた閾値を超えている区間を、生体が不安定いびきをしている不安定いびき区間と判定する。また、種類判定部１５０は、スコアの平均値が予め定められた閾値以下の区間を、生体が安定いびきをしている安定いびき区間と判定する。このように、種類判定部１５０は、生体が不安定いびきを発生しているか否かを判定し、これにより、対象信号は、いびき無し、安定いびき、および不安定いびきのいずれかに判定されることになる。
【００５１】
図１３は、本実施形態に係る種類判定部１５０が図３に示す対象信号に対していびきの種類を判定した結果の一例を示す。図１２に示す閾値を用いることで、種類判定部１５０は、区間ＡおよびＣを安定いびき、区間ＢおよびＤをいびき無し、区間Ｅを不安定いびきと判定することができる。
【００５２】
以上のように、本実施形態に係る推定装置１００は、非定常性特徴量を用いることにより、対象信号に対していびきの種類を推定することができる。また、本実施形態に係る呼吸状態判定装置１０（または種類判定装置）は、対象信号に対していびきの種類を判定することができる。図１３に示すような呼吸状態判定装置１０が判定したいびきの種類の判定結果は、専門家による判定結果と比較することで、９０％に至る一致率が得られた例もあり、専門家の判断程度の正確な結果を十数倍から数十倍の高速化させて得られることわかった。
【００５３】
次に、呼吸状態判定部１７０は、いびきの種類の判定結果に基づき、生体の呼吸状態を判定する（Ｓ３００）。呼吸状態判定部１７０は、時間軸上において、生体がいびき無しと判定された区間を挟む前後の２つの区間のいびきの種類に応じて、いびき無しの区間の呼吸状態を判定する。呼吸状態判定部１７０は、生体がいびきを発生していないと推定された第１区間の前後の第２区間および第３区間の音声の種類の判定結果に基づき、第１区間の生体の呼吸状態を判定する。呼吸状態判定部１７０は、複数の区間を、順次判定してよい。
【００５４】
呼吸状態判定部１７０は、第１区間の生体の呼吸状態が閉塞性無呼吸状態か否かを判定してよい。また、呼吸状態判定部１７０は、第１区間の生体の呼吸状態が無呼吸または低呼吸の状態であるか否かを判定してもよい。呼吸状態判定部１７０は、種類判定部１５０の第２区間および第３区間の判定結果のうち、１以上の不安定いびきの発生に応じて、第１区間の呼吸状態を判定する。例えば、呼吸状態判定部１７０は、種類判定部１５０が第２区間および第３区間のいずれか一方を不安定いびきと判定したことに応じて、第１区間の呼吸状態が閉塞性無呼吸状態と判定する。また、呼吸状態判定部１７０は、第２区間および第３区間が不安定いびきと判定されたことに応じて、第１区間の呼吸状態が閉塞性無呼吸状態と判定する。
【００５５】
また、種類判定部１５０は、第１区間の期間が予め定められた期間を超えることを更に条件として、第１区間の呼吸状態が閉塞性無呼吸状態と判定してもよい。種類判定部１５０は、一例として、第１区間が略１０秒を超える期間の場合に、当該第１区間の呼吸状態を閉塞性無呼吸状態と判定する。このように、呼吸状態判定部１７０は、不安定いびきの有無と、無音または寝息の状態の継続時間に応じて閉塞性無呼吸状態を判断するので、より正確に呼吸状態を判断することができる。
【００５６】
例えば、図１３において、種類判定部１５０は、区間ＡおよびＣを安定いびき区間に、区間Ｅを不安定いびき区間に判定したので、呼吸状態判定部１７０は、区間Ａ、Ｃ、およびＥによって挟まれる区間ＢおよびＤの呼吸状態を判定する。この場合、種類判定部１５０は、区間ＡおよびＣが安定いびき区間と判定されたことに応じて、区間Ｂの呼吸状態を有呼吸状態（寝息区間）と判定する。また、種類判定部１５０は、区間Ｃが安定いびき区間、および区間Ｅが不安定いびき区間と判定されたことに応じて、区間Ｄの呼吸状態を閉塞性無呼吸状態と判定する。
【００５７】
以上のように、本実施形態に係る呼吸状態判定装置１０は、いびきが検出されない区間の生体の呼吸状態を、いびきが検出される区間のいびきの種類に応じて判定するので、音がほとんど検出されない区間の呼吸状態を判定することができる。また、呼吸状態判定装置１０は、いびきが検出される区間のいびきの種類を、自動で精度よく判定することができるので、いびきが検出されない区間の呼吸状態も自動で精度よく判定することができる。
【００５８】
なお、専門家は、８時間等の予め定められた期間に、不安定いびきと判断したいびきの回数をカウントし、予め定められた回数以上カウントしたことに応じて、生体がＳＡＳであると診断する場合がある。そこで、呼吸状態判定部１７０は、複数の判定結果が、予め定められた時間の間に予め定められた回数以上の不安定いびき区間が発生していることに応じて、警告信号を出力してもよい。これによって、呼吸状態判定装置１０は、ＳＡＳセンサとして機能することができる。
【００５９】
また、この場合、呼吸状態判定装置１０は、ＳＡＳをリアルタイムでセンシングすることができる。これにより、呼吸状態判定装置１０は、例えば、自動圧調整式持続陽圧呼吸療法装置（Ａｕｔｏ−ＣＰＡＰ）、可動式枕、可動式ベッド等の外部の装置に直接警告信号を供給することで、生体の危険状態を直ちに回避させることもできる。
【００６０】
以上の本実施形態に係るスコア生成部１４２は、いびきの有無、音高の分散の大小、および音の強さの分散の大小を、予め定められた対応する閾値と比較してスコアを生成することを説明した。これに代えて、スコア生成部１４２は、非定常性特徴量と予め生成した複数のモデルのそれぞれとの照合結果に基づいて、音声の種類を推定するためのスコアを生成してもよい。
【００６１】
スコア生成部１４２は、例えば、隠れマルコフモデルとして知られる統計モデルを用いて、対象信号に対応するスコアを与える。また、スコア生成部１４２は、パターンマッチング、サポートベクターマシン等のパターン認識技術を用いて対象信号に対応するスコアを与えてもよい。スコア生成部１４２は、非定常性特徴量を用いるので、これらの技術を用いても、いびきの種類を推定するためのスコアを生成することができ、当該スコアを用いることで、呼吸状態判定装置１０は、精度の高いいびき種類の判定を実行することができる。
【００６２】
図１４は、本実施形態に係る呼吸状態判定装置１０の変形例を示す。本変形例の呼吸状態判定装置１０は、無呼吸低呼吸指数（ＡＨＩ：Ａｐｎｅａ  Ｈｙｐｏｐｎｅａ  Ｉｎｄｅｘ）、睡眠障害指数（ＲＤＩ：Ｒｅｓｐｉｒａｔｏｒｙ  Ｄｉｓｔｕｒｂａｎｃｅ  Ｉｎｄｅｘ）等のＳＡＳの診断パラメータを出力する。本変形例の呼吸状態判定装置１０は、診断パラメータ算出部２００を更に備える。
【００６３】
本変形例の呼吸状態判定部１７０は、生体がいびきを発生していないと推定された第１区間の前後の第２区間および第３区間のうち、少なくとも一方の区間を種類判定部１５０が不安定いびきと判定したことに応じて、当該生体は第１区間において無呼吸または低呼吸の状態であると判定する。診断パラメータ算出部２００は、呼吸状態判定部１７０に接続され、呼吸状態判定部１７０の判定結果に基づいて、睡眠時無呼吸症候群の診断パラメータを算出する。診断パラメータ算出部２００は、無呼吸低呼吸指数および睡眠障害指数の少なくとも１つを診断パラメータとして算出してよい。
【００６４】
診断パラメータ算出部２００は、例えば、時間的に連続する無呼吸または低呼吸と判断された区間を一つの領域として計数し、１時間当たりの計数値の平均値をＡＨＩとして算出する。診断パラメータ算出部２００は、バッファ部１６０に蓄積された判定結果に基づき、予め定められた時間以上の呼吸状態判定部１７０の判定結果を用いて診断パラメータを算出する。診断パラメータ算出部２００は、生体の全睡眠時間の呼吸状態判定部１７０による判定結果を用いることが好ましい。
【００６５】
また、診断パラメータ算出部２００は、種類判定部１５０の第２区間および第３区間の判定結果に基づく予備パラメータを算出してもよい。診断パラメータ算出部２００は、無呼吸低呼吸指数および睡眠障害指数の他に、症状の診断の参考となる情報等を予備パラメータとして算出してよい。診断パラメータ算出部２００は、例えば、不安定いびきと判定された期間の音高（最高値および／または平均値）、寝入り時間および寝入り時間からの経過時間等の情報を、予備パラメータとして算出して出力する。このように、本変形例の呼吸状態判定装置１０は、対象信号に基づいて、睡眠時無呼吸症候群の診断パラメータおよび予備パラメータを出力することができる。
【００６６】
以上で説明した呼吸状態判定装置１０は、生体の音声信号または音声に起因する信号であれば、当該音声がいびきであるか否か、そして不安定いびきであるか否かを判断することができる。したがって、呼吸状態判定装置１０を生体に対して非接触の装置とすることもできる。また、呼吸状態判定装置１０は、生体の音声信号または音声に起因する信号を取得する、例えば、マイク、ＩＣレコーダ等に搭載することもできる。また、呼吸状態判定装置１０は、マイク、ＩＣレコーダ等に予め記憶された音声信号等に対して判定することもできる。
【００６７】
呼吸状態判定装置１０は、このような生体の音声に起因する信号に基づき、睡眠中の雑音等といびきを分別してから呼吸状態を判定するので、正確かつ簡便に診断することができる。また、呼吸状態判定装置１０は、呼吸状態を判定する区間の前後のいびき種類に応じて、呼吸状態を判定するので、より正確な判定を実行することができる。
【００６８】
図１５は、本実施形態に係る呼吸状態判定装置１０として機能するコンピュータ１９００のハードウェア構成の一例を示す。本実施形態に係るコンピュータ１９００は、ホスト・コントローラ２０８２により相互に接続されるＣＰＵ２０００、ＲＡＭ２０２０、グラフィック・コントローラ２０７５、および表示装置２０８０を有するＣＰＵ周辺部と、入出力コントローラ２０８４によりホスト・コントローラ２０８２に接続される通信インターフェイス２０３０、ハードディスクドライブ２０４０、およびＤＶＤドライブ２０６０を有する入出力部と、入出力コントローラ２０８４に接続されるＲＯＭ２０１０、フレキシブルディスク・ドライブ２０５０、および入出力チップ２０７０を有するレガシー入出力部と、を備える。
【００６９】
ホスト・コントローラ２０８２は、ＲＡＭ２０２０と、高い転送レートでＲＡＭ２０２０をアクセスするＣＰＵ２０００およびグラフィック・コントローラ２０７５とを接続する。ＣＰＵ２０００は、ＲＯＭ２０１０およびＲＡＭ２０２０に格納されたプログラムに基づいて動作し、各部の制御を行う。グラフィック・コントローラ２０７５は、ＣＰＵ２０００等がＲＡＭ２０２０内に設けたフレーム・バッファ上に生成する画像データを取得し、表示装置２０８０上に表示させる。これに代えて、グラフィック・コントローラ２０７５は、ＣＰＵ２０００等が生成する画像データを格納するフレーム・バッファを、内部に含んでもよい。
【００７０】
入出力コントローラ２０８４は、ホスト・コントローラ２０８２と、比較的高速な入出力装置である通信インターフェイス２０３０、ハードディスクドライブ２０４０、ＤＶＤドライブ２０６０を接続する。通信インターフェイス２０３０は、ネットワークを介して他の装置と通信する。ハードディスクドライブ２０４０は、コンピュータ１９００内のＣＰＵ２０００が使用するプログラムおよびデータを格納する。ＤＶＤドライブ２０６０は、ＤＶＤ−ＲＯＭ２０９５からプログラムまたはデータを読み取り、ＲＡＭ２０２０を介してハードディスクドライブ２０４０に提供する。
【００７１】
また、入出力コントローラ２０８４には、ＲＯＭ２０１０と、フレキシブルディスク・ドライブ２０５０、および入出力チップ２０７０の比較的低速な入出力装置とが接続される。ＲＯＭ２０１０は、コンピュータ１９００が起動時に実行するブート・プログラム、および／または、コンピュータ１９００のハードウェアに依存するプログラム等を格納する。フレキシブルディスク・ドライブ２０５０は、フレキシブルディスク２０９０からプログラムまたはデータを読み取り、ＲＡＭ２０２０を介してハードディスクドライブ２０４０に提供する。入出力チップ２０７０は、フレキシブルディスク・ドライブ２０５０を入出力コントローラ２０８４へと接続すると共に、例えばパラレル・ポート、シリアル・ポート、キーボード・ポート、マウス・ポート等を介して各種の入出力装置を入出力コントローラ２０８４へと接続する。
【００７２】
ＲＡＭ２０２０を介してハードディスクドライブ２０４０に提供されるプログラムは、フレキシブルディスク２０９０、ＤＶＤ−ＲＯＭ２０９５、またはＩＣカード等の記録媒体に格納されて利用者によって提供される。プログラムは、記録媒体から読み出され、ＲＡＭ２０２０を介してコンピュータ１９００内のハードディスクドライブ２０４０にインストールされ、ＣＰＵ２０００において実行される。
【００７３】
プログラムは、コンピュータ１９００にインストールされ、コンピュータ１９００を信号取得部１１０、特徴量算出部１２０、非定常性特徴量算出部１３０、推定部１４０、スコア生成部１４２、記憶部１４４、種類判定部１５０、バッファ部１６０、呼吸状態判定部１７０、および診断パラメータ算出部２００として機能させる。
【００７４】
プログラムに記述された情報処理は、コンピュータ１９００に読込まれることにより、ソフトウェアと上述した各種のハードウェア資源とが協働した具体的手段である信号取得部１１０、特徴量算出部１２０、非定常性特徴量算出部１３０、推定部１４０、スコア生成部１４２、記憶部１４４、種類判定部１５０、バッファ部１６０、呼吸状態判定部１７０、および診断パラメータ算出部２００として機能する。そして、この具体的手段によって、本実施形態におけるコンピュータ１９００の使用目的に応じた情報の演算または加工を実現することにより、使用目的に応じた特有の呼吸状態判定装置１０が構築される。
【００７５】
一例として、コンピュータ１９００と外部の装置等との間で通信を行う場合には、ＣＰＵ２０００は、ＲＡＭ２０２０上にロードされた通信プログラムを実行し、通信プログラムに記述された処理内容に基づいて、通信インターフェイス２０３０に対して通信処理を指示する。通信インターフェイス２０３０は、ＣＰＵ２０００の制御を受けて、ＲＡＭ２０２０、ハードディスクドライブ２０４０、フレキシブルディスク２０９０、またはＤＶＤ−ＲＯＭ２０９５等の記憶装置上に設けた送信バッファ領域等に記憶された送信データを読み出してネットワークへと送信し、もしくは、ネットワークから受信した受信データを記憶装置上に設けた受信バッファ領域等へと書き込む。このように、通信インターフェイス２０３０は、ＤＭＡ（ダイレクト・メモリ・アクセス）方式により記憶装置との間で送受信データを転送してもよく、これに代えて、ＣＰＵ２０００が転送元の記憶装置または通信インターフェイス２０３０からデータを読み出し、転送先の通信インターフェイス２０３０または記憶装置へとデータを書き込むことにより送受信データを転送してもよい。
【００７６】
また、ＣＰＵ２０００は、ハードディスクドライブ２０４０、ＤＶＤドライブ２０６０（ＤＶＤ−ＲＯＭ２０９５）、フレキシブルディスク・ドライブ２０５０（フレキシブルディスク２０９０）等の外部記憶装置に格納されたファイルまたはデータベース等の中から、全部または必要な部分をＤＭＡ転送等によりＲＡＭ２０２０へと読み込ませ、ＲＡＭ２０２０上のデータに対して各種の処理を行う。そして、ＣＰＵ２０００は、処理を終えたデータを、ＤＭＡ転送等により外部記憶装置へと書き戻す。このような処理において、ＲＡＭ２０２０は、外部記憶装置の内容を一時的に保持するものとみなせるから、本実施形態においてはＲＡＭ２０２０および外部記憶装置等をメモリ、記憶部、または記憶装置等と総称する。本実施形態における各種のプログラム、データ、テーブル、データベース等の各種の情報は、このような記憶装置上に格納されて、情報処理の対象となる。なお、ＣＰＵ２０００は、ＲＡＭ２０２０の一部をキャッシュメモリに保持し、キャッシュメモリ上で読み書きを行うこともできる。このような形態においても、キャッシュメモリはＲＡＭ２０２０の機能の一部を担うから、本実施形態においては、区別して示す場合を除き、キャッシュメモリもＲＡＭ２０２０、メモリ、および／または記憶装置に含まれるものとする。
【００７７】
また、ＣＰＵ２０００は、ＲＡＭ２０２０から読み出したデータに対して、プログラムの命令列により指定された、本実施形態中に記載した各種の演算、情報の加工、条件判断、情報の検索・置換等を含む各種の処理を行い、ＲＡＭ２０２０へと書き戻す。例えば、ＣＰＵ２０００は、条件判断を行う場合においては、本実施形態において示した各種の変数が、他の変数または定数と比較して、大きい、小さい、以上、以下、等しい等の条件を満たすかどうかを判断し、条件が成立した場合（または不成立であった場合）に、異なる命令列へと分岐し、またはサブルーチンを呼び出す。
【００７８】
また、ＣＰＵ２０００は、記憶装置内のファイルまたはデータベース等に格納された情報を検索することができる。例えば、第１属性の属性値に対し第２属性の属性値がそれぞれ対応付けられた複数のエントリが記憶装置に格納されている場合において、ＣＰＵ２０００は、記憶装置に格納されている複数のエントリの中から第１属性の属性値が指定された条件と一致するエントリを検索し、そのエントリに格納されている第２属性の属性値を読み出すことにより、所定の条件を満たす第１属性に対応付けられた第２属性の属性値を得ることができる。
【００７９】
以上に示したプログラムまたはモジュールは、外部の記録媒体に格納されてもよい。記録媒体としては、フレキシブルディスク２０９０、ＤＶＤ−ＲＯＭ２０９５の他に、ＤＶＤ、Ｂｌｕ−ｒａｙ（登録商標）、またはＣＤ等の光学記録媒体、ＭＯ等の光磁気記録媒体、テープ媒体、ＩＣカード等の半導体メモリ等を用いることができる。また、専用通信ネットワークまたはインターネットに接続されたサーバシステムに設けたハードディスクまたはＲＡＭ等の記憶装置を記録媒体として使用し、ネットワークを介してプログラムをコンピュータ１９００に提供してもよい。
【００８０】
以上、本発明を実施の形態を用いて説明したが、本発明の技術的範囲は上記実施の形態に記載の範囲には限定されない。上記実施の形態に、多様な変更または改良を加えることが可能であることが当業者に明らかである。その様な変更または改良を加えた形態も本発明の技術的範囲に含まれ得ることが、特許請求の範囲の記載から明らかである。
【００８１】
特許請求の範囲、明細書、および図面中において示した装置、システム、プログラム、および方法における動作、手順、ステップ、および段階等の各処理の実行順序は、特段「より前に」、「先立って」等と明示しておらず、また、前の処理の出力を後の処理で用いるのでない限り、任意の順序で実現しうることに留意すべきである。特許請求の範囲、明細書、および図面中の動作フローに関して、便宜上「まず、」、「次に、」等を用いて説明したとしても、この順で実施することが必須であることを意味するものではない。
【符号の説明】
【００８２】
１０  呼吸状態判定装置、１００  推定装置、１１０  信号取得部、１２０  特徴量算出部、１３０  非定常性特徴量算出部、１４０  推定部、１４２  スコア生成部、１４４  記憶部、１５０  種類判定部、１６０  バッファ部、１７０  呼吸状態判定部、２００  診断パラメータ算出部、１９００  コンピュータ、２０００  ＣＰＵ、２０１０  ＲＯＭ、２０２０  ＲＡＭ、２０３０  通信インターフェイス、２０４０  ハードディスクドライブ、２０５０  フレキシブルディスク・ドライブ、２０６０  ＤＶＤドライブ、２０７０  入出力チップ、２０７５  グラフィック・コントローラ、２０８０  表示装置、２０８２  ホスト・コントローラ、２０８４  入出力コントローラ、２０９０  フレキシブルディスク、２０９５  ＤＶＤ−ＲＯＭ
【図１】
0006535186.tif 000003
【図２】
0006535186.tif 000004
【図３】
0006535186.tif 000005
【図４】
0006535186.tif 000006
【図５】
0006535186.tif 000007
【図６】
0006535186.tif 000008
【図７】
0006535186.tif 000009
【図８】
0006535186.tif 000010
【図９】
0006535186.tif 000011
【図１０】
0006535186.tif 000012
【図１１】
0006535186.tif 000013
【図１２】
0006535186.tif 000014
【図１３】
0006535186.tif 000015
【図１４】
0006535186.tif 000016
【図１５】
0006535186.tif 000017