(57)【要約】【課題】睡眠時無呼吸予測方法を提供する。【解決手段】本発明は、ニューラルネットワークを用いた機械学習を利用していびきから睡眠時の無呼吸を予測する方法である。これは、原信号をキャプチャするステップａと、原信号中から少なくとも１つのいびき信号を検索し、且ついびき信号を分割した後、一次元特徴ベクトルを有するいびき信号ベクトルを出力するいびき信号分割アルゴリズム実行ステップｂと、いびき信号の特徴からいびき信号ベクトルを抽出し、いびき信号ベクトルを、二次元特徴ベクトルを有するいびき機能マトリックスに変換する特徴抽出アルゴリズム実行ステップｃと、いびき機能マトリックスを分類し、いびき機能マトリックスが分類学習を経た後、いびき信号の意味を解釈して呼吸停止回数または呼吸低下回数を予測するニューラルネットワークアルゴリズム実行ステップｄと、を含む。【選択図】図１
2020185390.tif 000002
【特許請求の範囲】
【請求項１】
原信号をキャプチャするステップと、
前記原信号中から少なくとも１つのいびき信号を検索し、且つ前記いびき信号を分割した後、一次元特徴ベクトルを有するいびき信号ベクトルを出力するいびき信号分割アルゴリズム実行ステップと、
前記いびき信号の特徴から前記いびき信号ベクトルを抽出し、前記いびき信号ベクトルを、二次元特徴ベクトルを有するいびき機能マトリックスに変換する特徴抽出アルゴリズム実行ステップと、
前記いびき機能マトリックスを分類し、前記いびき機能マトリックスが分類学習を経た後、前記いびき信号の意味を解釈して呼吸停止回数または呼吸低下回数を予測するニューラルネットワークアルゴリズム実行ステップと、を含むことを特徴とする睡眠時無呼吸予測方法。
【請求項２】
前記いびき信号分割アルゴリズム実行ステップのいびき信号分割アルゴリズムは、第一閾値及び第二閾値を分割基準とし、且つスライディングウィンドウを設定して前記原信号のリニア検査を行い、前記スライディングウィンドウ内の前記原信号の最大値を計算するステップを更に含み、
前記原信号の最大値が前記第二閾値より大きい場合、前記いびき信号が発生したと認定し、前記スライディングウィンドウの停留位置を前記いびき信号の位置として設定し、
次いで、前記スライディングウィンドウが右に向けて継続的に前記いびき信号のリニア検査を行い、同時に前記スライディングウィンドウ内で前記いびき信号の絶対値の合計を計算し、算出された前記いびき信号の絶対値の合計が前記第一閾値未満である場合、前記スライディングウィンドウの停留位置を前記いびき信号の右側停止位置として設定し、
その後、前記スライディングウィンドウが左に向けて継続的に前記いびき信号のリニア検査を行い、同時に前記スライディングウィンドウ内で前記いびき信号の絶対値の合計を計算し、算出された前記いびき信号の絶対値の合計が前記第一閾値未満である場合、前記スライディングウィンドウの停留位置を前記いびき信号の左側停止位置として設定し、且つ前記いびき信号の前記右側停止位置及び前記いびき信号の前記左側停止位置から前記いびき信号ベクトルを分割することを特徴とする請求項１に記載の睡眠時無呼吸予測方法。
【請求項３】
前記第一閾値の式は次式で表され、
【数１】
2020185390.tif 000012
Ｍは前記第一閾値を示し、meanは平均値を示し、ｆ()はダウンサンプリング関数を示し、Ｙｉは前記原信号のベクトルを示し、
また、前記第二閾値の式は次式で表され、
【数２】
2020185390.tif 000013
【数３】
2020185390.tif 000014
Ｘは前記第二閾値を示し、meanは平均値を示し、stdは標準偏差を示し、Ｎは自然数を示し、sortは小から大までの順序を示し、式（３）は前記Ｙｉが分割されたｎ個のベクトルを示すことを特徴とする請求項２に記載の睡眠時無呼吸予測方法。
【請求項４】
前記いびき信号ベクトルの長さは25000に設定されることを特徴とする請求項２に記載の睡眠時無呼吸予測方法。
【請求項５】
前記スライディングウィンドウの長さは1000に設定されることを特徴とする請求項２に記載の睡眠時無呼吸予測方法。
【請求項６】
前記特徴抽出アルゴリズムはメル周波数ケプストラム係数アルゴリズムであり、
前記メル周波数ケプストラム係数アルゴリズムはプリエンファシスプロセスと、フレーミングとウィンドウプロセスと、高速フーリエ変換プロセスと、メルフィルターバンクと、非線形変換と、離散コサイン変換プロセスと、を含むことを特徴とする請求項１に記載の睡眠時無呼吸予測方法。
【請求項７】
前記ニューラルネットワークアルゴリズムは畳み込みニューラルネットワークアルゴリズムであり、
前記畳み込みニューラルネットワークアルゴリズムではDenseブロックネットワークモデルが識別モデルとして選択されることを特徴とする請求項１に記載の睡眠時無呼吸予測方法。
【請求項８】
前記Denseブロックネットワークモデルは複数のDenseブロックと、複数のTransitionレイヤーと、分類レイヤーとを含むことを特徴とする請求項７に記載の睡眠時無呼吸予測方法。
【請求項９】
前記Transitionレイヤーは畳み込みプロセス及びプーリングプロセスを含み、
前記分類レイヤーはSoftmaxレイヤーであることを特徴とする請求項８に記載の睡眠時無呼吸予測方法。
【請求項１０】
前記DenseブロックはDenseレイヤーと、BR BN-ReLU-Convレイヤーと、成長率とを含むことを特徴とする請求項８に記載の睡眠時無呼吸予測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ニューラルネットワークを用いた機械学習を利用していびきから睡眠時の無呼吸を予測する方法（A method for predicting sleep apnea from neural networks）に関し、更に詳しくはニューラルネットワークを用いた機械学習を利用していびき信号の意味を解釈して呼吸停止回数及び呼吸低下回数を予測する方法に関する。
【背景技術】
【０００２】
いびき及び閉塞性睡眠時無呼吸症候群（obstructive sleep apnea, OSA）は睡眠時に上気道の一部が閉塞されることで起こり、共に同じ疾病に属するが、重症度は異なり、少なくとも人口の４％以上が罹患している。罹患すると日中に疲労や眠気を感じたり、記憶力が減退したり、うつ病に罹ったりする。さらには交通事故の発生率を上昇させたり、精神面の疾患、動脈性高血圧、心血管の疾病、脳卒中、及びメタボリックシンドロームなどを引き起こしたりすることもある。
【０００３】
現在の臨床上の診断基準は、睡眠ポリグラフ検査に基づいている。この検査は、病院に一泊して睡眠をとり、睡眠効率及び睡眠の深さを観測するものであり、睡眠時の呼吸停止回数及び血中酸素濃度等の検査項目が含まれる。しかしながら、睡眠ポリグラフ検査は病院で行う必要があるが、これは日数がかかる上、費用も人手もかかり、且つ検査費用も高かった。このため、医師や一般市民にとっては、迅速且つ簡単に診断を行うことができ、睡眠時の呼吸停止の重症度を毎日測れる道具があれば、非常に有益である。
【０００４】
いびきは閉塞性睡眠時無呼吸症候群の１つの主要な症状であり、患者のうち約９５％にはいびきの症状がある。いびきをかく患者達にとっては、自宅で自らいびき検知ができれば、便利で有用である。単純ないびき及び睡眠時無呼吸症候群によるいびきについて、頻度及び音の違い等を基に、過去に分析がなされており、睡眠時の閉塞部位の評価に応用されている。しかしながら、過去の医師や研究では、いびきから閉塞性睡眠時無呼吸症候群による睡眠時の呼吸停止（Sleep apnea: SA）回数や、睡眠時の呼吸低下（Sleep hypopnea: SH）回数を測ることできなかった。なお、１時間あたりの平均呼吸停止回数及び呼吸低下回数は「無呼吸低呼吸指数」（Apnea-hypopnea Index: AHI）と呼ばれている。
【発明の概要】
【発明が解決しようとする課題】
【０００５】
そこで、本発明者は上記従来技術の欠点が改善可能と考えた。ニューラルネットワークを用いた機械学習を利用して睡眠時の呼吸停止回数または睡眠時の呼吸低下回数を分類学習させることが、本発明の課題である。
【０００６】
本発明は、上述の従来技術の問題に鑑みてなされたものである。上記の課題解決のため、本発明は、ニューラルネットワークを用いた機械学習を利用していびきから睡眠時の無呼吸を予測する方法を提供することを目的とする。
換言すれば、本発明はニューラルネットワークを用いた機械学習により睡眠時の呼吸停止回数または睡眠時の呼吸低下回数を分類学習し、過去の医師や研究がいびきから睡眠時の呼吸停止回数や呼吸低下回数を予測できなかったという問題点を解決し、いびき信号の意味を解釈して呼吸停止回数または呼吸低下回数により人間の睡眠時の呼吸が正常であるか病気を患っているかを予測するという効果を提供する。
【課題を解決するための手段】
【０００７】
上記課題を解決するために、本発明のある態様の睡眠時無呼吸予測方法は、
原信号をキャプチャするステップと、
前記原信号中から少なくとも１つのいびき信号をキャプチャし、且つ前記いびき信号を分割した後、一次元特徴ベクトルを有するいびき信号ベクトルを出力するいびき信号分割アルゴリズム実行ステップと、
前記いびき信号の特徴から前記いびき信号ベクトルを抽出し、前記いびき信号ベクトルを、二次元特徴ベクトルを有するいびき機能マトリックスに変換する特徴抽出アルゴリズム実行ステップと、
前記いびき機能マトリックスを分類し、前記いびき機能トリックスが分類学習を経た後、前記いびき信号の意味を解釈して呼吸停止回数または呼吸低下回数を予測するニューラルネットワークアルゴリズム実行ステップと、を含む。
【０００８】
本発明の好適例において、前記いびき信号分割アルゴリズムは、第一閾値及び第二閾値を分割基準とし、且つスライディングウィンドウを設定して前記原信号のリニア検査を行い、前記スライディングウィンドウ内の前記原信号の最大値を計算するステップを更に含む。
前記原信号の最大値が前記第二閾値より大きい場合、前記いびき信号が発生したと認定し、前記スライディングウィンドウの停留位置を前記いびき信号の位置として設定する。
次いで、前記スライディングウィンドウが右に向けて継続的に前記いびき信号のリニア検査を行い、同時に前記スライディングウィンドウ内で前記いびき信号の絶対値の合計を計算し、算出された前記いびき信号の絶対値の合計が前記第一閾値未満である場合、前記スライディングウィンドウの停留位置を前記いびき信号の右側停止位置として設定する。
その後、前記スライディングウィンドウが左に向けて継続的に前記いびき信号のリニア検査を行い、同時に前記スライディングウィンドウ内で前記いびき信号の絶対値の合計を計算し、算出された前記いびき信号の絶対値の合計が前記第一閾値未満である場合、前記スライディングウィンドウの停留位置を前記いびき信号の左側停止位置として設定し、且つ前記いびき信号の前記右側停止位置及び前記いびき信号の前記左側停止位置から前記いびき信号ベクトルを分割する。
【０００９】
本発明の好適例において、前記第一閾値の式は式（１）で表され、
【数１】
2020185390.tif 000003
Ｍは前記第一閾値を示し、meanは平均値を示し、ｆ()はダウンサンプリング関数を示し、Ｙｉは前記原信号のベクトルを示す。
また、前記第二閾値の式は式（２）で表され、
【数２】
2020185390.tif 000004
【数３】
2020185390.tif 000005
Ｘは前記第二閾値を示し、meanは平均値を示し、stdは標準偏差を示し、Ｎは自然数を示し、sortは小から大までの順序を示し、式（３）はＹｉが分割されたn個のベクトル示す。
【００１０】
本発明の好適例において、前記いびき信号ベクトルの長さは25000に設定されている。
【００１１】
本発明の好適例において、前記スライディングウィンドウの長さは1000に設定されている。
【００１２】
本発明の好適例において、前記特徴抽出アルゴリズムはメル周波数ケプストラム係数アルゴリズムであり、
前記メル周波数ケプストラム係数アルゴリズムはプリエンファシスプロセスと、フレーミングとウィンドウプロセスと、高速フーリエ変換プロセスと、メルフィルターバンクと、非線形変換と、離散コサイン変換プロセスと、を含む。
【００１３】
本発明の好適例において、前記ニューラルネットワークアルゴリズムは畳み込みニューラルネットワークアルゴリズムであり、
前記畳み込みニューラルネットワークアルゴリズムはDenseブロックネットワークモデルが識別モデルとして選択されている。
【００１４】
本発明の好適例において、前記Denseブロックネットワークモデルは複数のDenseブロックと、複数のTransitionレイヤーと、分類レイヤーと、を含む。
【００１５】
本発明の好適例において、前記Transitionレイヤーは畳み込みプロセス及びプーリングプロセスを含む。前記分類レイヤーはSoftmaxレイヤーである。
【００１６】
本発明の好適例において、前記DenseブロックはDenseレイヤーと、BR BN-ReLU-Convレイヤーと、成長率と、を含む。
【発明の効果】
【００１７】
上記技術手段により、前記いびき信号分割アルゴリズム、前記特徴抽出アルゴリズム、及び前記ニューラルネットワークアルゴリズムを統合することで、前記原信号を有効に処理し、前記いびき信号の呼吸停止回数または呼吸低下回数を予測する。前記いびき信号の呼吸停止回数または呼吸低下回数から人間の睡眠時の呼吸が正常であるか病気にかかっているかを予測し、過去の医師や研究がいびきから睡眠時の呼吸停止回数や呼吸低下回数を予測できなかったという問題点を解決する。
【図面の簡単な説明】
【００１８】
【図１】本発明の一実施形態を説明するためのフローチャートである。
【図２】本発明の一実施形態を説明するためのブロック図である。
【図３Ａ】本発明の一実施形態に係る原信号のグラフである。
【図３Ｂ】本発明の一実施形態に係る正規化後の原信号のグラフである。
【図３Ｃ】本発明の一実施形態の一部の原信号のグラフである。
【図３Ｄ】本発明の一実施形態に係る正規化後の原信号といびき信号を検索中のグラフである。
【図３Ｅ】本発明の一実施形態に係るいびき信号分割を示すグラフである
【図３Ｆ】本発明の第一単一いびき信号ベクトルが分割されたグラフである。
【図３Ｇ】本発明の第二単一いびき信号ベクトルが分割されたグラフである。
【図３Ｈ】本発明の第三単一いびき信号ベクトルが分割されたグラフである。
【図３Ｉ】本発明の第四単一いびき信号ベクトルが分割されたグラフである。
【図３Ｊ】本発明の第五単一いびき信号ベクトルが分割されたグラフである。
【図３Ｋ】本発明の第六単一いびき信号ベクトルが分割されたグラフである。
【図３Ｌ】本発明の第七単一いびき信号ベクトルが分割されたグラフである。
【図４】本発明に係るメル周波数ケプストラム係数アルゴリズムを示すブロック図である。
【図５】本発明に係る緻密ネットワークモジュールを示すブロック図である。
【図６】本発明のDenseブロックを示す模式図である。
【図７Ａ】本発明の一実施形態に係る睡眠時無呼吸モデル識別のグラフである。
【図７Ｂ】本発明の一実施形態に係る睡眠時無呼吸モデル識別により正常信号または病気信号を予測するグラフである。
【発明を実施するための形態】
【００１９】
以下、図面を参照しつつ、本発明に係る睡眠時無呼吸予測方法について具体的に説明する。なお、本発明は、以下に説明する実施形態に限定されるものではない。
【００２０】
（実施形態）
以下、本発明を図１〜７Ｂに基づいて説明する。本発明の睡眠時無呼吸予測方法は、以下のステップから構成されている。
【００２１】
＜ステップａ＞
原信号（Ｙ）をキャプチャする。本実施例では、各患者に対し睡眠ポリグラフ検査（PSG）を行う。
睡眠研究変量は無呼吸低呼吸指数（Apnea-hypopnea index, AHI）評価と、いびき指数と、最小酸素飽和度（MOS）と、を含む。
前記無呼吸低呼吸指数評価とは、睡眠１時間あたりの閉塞性呼吸停止（Apnea）及び呼吸低下（hypopnea）発作の総数を指す。
呼吸停止（Apnea）とは、少なくとも１０秒間呼吸が停止した状態を指す。
呼吸低下（hypopnea）とは、基線通気値が５０％以上低下した状態が１０秒以上続き、且つ酸素飽和度が４％以上低下した状態を指す。
PSG実施中は、胸骨上切痕上方に取り付けた小型マイクのデジタル音量計を使用していびきを記録するが、この方法に限定するわけではない。
【００２２】
＜ステップｂ＞
いびき信号分割アルゴリズム（Ｇ1）により前記原信号（Ｙ）中から少なくとも１つのいびき信号（Ｂ）をキャプチャし、且つ前記いびき信号（Ｂ）を分割した後、一次元特徴ベクトルを有するいびき信号ベクトル（Ｓ）を出力する。
前記原信号（Ｙ）は終夜録音され、学習前にはデータに対する前処理を施す必要がある。また、目標はいびきであるため、音声周波数中の各いびきを自動的に分割する必要があり、全自動分割アルゴリズムを設計し、最初の原信号を一段一段いびきの段に分割するが、この方法に限定するわけではない。
【００２３】
上述したように、本実施例では、前記いびき信号分割アルゴリズム（Ｇ1）は、第一閾値（Ｍ）及び第二閾値（Ｘ）を分割基準とし、スライディングウィンドウ（Ｗ）を設定して前記原信号（Ｙ）のリニア検査を行うステップを更に含む。
図３Ｄに示されるように、前記スライディングウィンドウ（Ｗ）内の前記原信号（Ｙ）の最大値（Ｘｉ）の計算を行う。
前記原信号（Ｙ）の最大値（Ｘｉ）が前記第二閾値（Ｘ）より大きい場合、前記いびき信号（Ｂ）が発生したと認定し、前記スライディングウィンドウ（Ｗ）の停留位置を前記いびき信号（Ｂ）の位置と設定する。
次いで、前記スライディングウィンドウ（Ｗ）が右に向けて継続的に前記いびき信号（Ｂ）のリニア検査を行い（図３Ｅを参照）、同時に前記スライディングウィンドウ（Ｗ）内で前記いびき信号の絶対値の合計（Ｍｉ）を計算して算出された絶対値の合計（Ｍｉ）が前記第一閾値（Ｍ）未満である場合、前記スライディングウィンドウ（Ｗ）の停留位置を前記いびき信号（Ｂ）の右側停止位置（Ｒ）と設定する。
そして、前記スライディングウィンドウ（Ｗ）が左に向けて継続的に前記いびき信号（Ｂ）のリニア検査を行い、同時に前記スライディングウィンドウ（Ｗ）内で前記いびき信号の絶対値の合計（Ｍｉ）の計算を行い、算出された絶対値の合計（Ｍｉ）が前記第一閾値（Ｍ）未満である場合、前記スライディングウィンドウ（Ｗ）の停留位置を前記いびき信号（Ｂ）の左側停止位置（Ｌ）と設定し、且つ前記いびき信号（Ｂ）の右側停止位置（Ｒ）及び前記いびき信号（Ｂ）の左側停止位置（Ｌ）から前記いびき信号ベクトル（Ｓ）を分割する。
図３Ｆ〜図３Ｌに示されるように、分割後に第一単一いびき信号ベクトル（Ｓ1）、第二単一いびき信号ベクトル（Ｓ2）、第三単一いびき信号ベクトル（Ｓ3）、
第四単一いびき信号ベクトル（Ｓ4）、第五単一いびき信号ベクトル（Ｓ5）、
第六単一いびき信号ベクトル（Ｓ6）、及び第七単一のいびき信号ベクトル（Ｓ7）が形成される。
また、各前記単一のいびき信号ベクトル（Ｓ1〜Ｓ7）は同じ長さである必要があるため、前記いびき信号ベクトル（Ｓ）をキャプチャする長さを調整し、前記いびき信号ベクトル（Ｓ）の長さを25000に設定するが、この方法に限定するわけではない。
【００２４】
上述したように、本実施例では、前記第一閾値（Ｍ）の式（１）は下記に示す。
【数１】
2020185390.tif 000006
ここでは、Ｍは前記第一閾値（Ｍ）を示し、meanは平均値を示し、ｆ()はダウンサンプリング関数を示し、Ｙｉは原信号（Ｙ）を２分間毎の音声周波数に分割したフレームベクトル（Ｙｉ）を示す。前記フレームベクトル（Ｙｉ）をサンプリングが必要な次元まで等分割する方式で４００次元までダウンサンプリングを行い、その後各段全ての最大値のみ取得すると、元の前記フレームベクトル（Ｙｉ）が1*400ベクトルまで縮小される。ダウンサンプリングにより更に信頼性の高い前記第一閾値（Ｍ）を獲得できるが、この方法に限定するわけではない。
前記第二閾値（Ｘ）の式（２）及び（３）は下記に示す。
【数２】
2020185390.tif 000007
【数３】
2020185390.tif 000008
ここでは、Ｘは前記第二閾値（Ｘ）を示し、meanは平均値を示し、stdは標準偏差を示し、Ｎは自然数を示し、sortは小から大までの順序を示し、式（３）はＹｉを等分割したｎ個のベクトルを示す。換言すれば、前記ｎ個のベクトルは前記フレームベクトル（Ｙｉ）の長さから前記スライディングウィンドウ（Ｗ）の大きさを除算したものであるから、前記スライディングウィンドウ（Ｗ）の長さが1000に設定され、且つ式（３）から前記自然数（Ｎ）を獲得した後、式（２）に代入されれば前記第二閾値（Ｘ）が算出されるが、この方法に限定するわけではない。
【００２５】
＜ステップｃ＞
特徴抽出アルゴリズム（Ｇ2）により前記いびき信号ベクトル（Ｓ）の特徴を抽出し、前記いびき信号ベクトル（Ｓ）を、二次元特徴ベクトルを有するいびき機能マトリックス（Ａ）に変換する。これにより、全て分割が完了した後、前記原信号（Ｙ）が複数の各前記単一のいびき信号ベクトル（Ｓ1〜Ｓ7）に分割され、その後、各前記単一のいびき信号ベクトル（Ｓ1〜Ｓ7）に前記特徴抽出アルゴリズムを使用する。
前記特徴抽出アルゴリズム（Ｇ2）はメル周波数ケプストラム係数（Mel-Frequency Cepstral Coefficients, MFCC）アルゴリズムにより特徴の抽出を行う。これは図４に示されるように、プリエンファシスプロセス（Ｇ21）と、フレーミングとウィンドウプロセス（Ｇ22、窓関数プロセス）と、高速フーリエ変換プロセス（Ｇ23）と、メルフィルターバンク（Ｇ24）と、非線形変換プロセス（Ｇ25）と、離散コサイン変換プロセス（Ｇ26）と、を含むが、この方法に限定するわけではない。
【００２６】
＜Ｇ21ステップ＞
前記プリエンファシスプロセス（Ｇ21）では、音声中の高周波部分を補償するためにプリエンファシスが用いられ、信号の高周波が強調される。式（４）にプリエンファシスを示す。
【数４】
2020185390.tif 000009
ここでは、Hpreemはプリエンファシス処理後の出力であり、αpreemは入力音声信号である。
【００２７】
＜Ｇ22ステップ＞
前記フレーミングとウィンドウプロセス（Ｇ22、窓関数プロセス）では、まず信号をフレーム毎に分割する。すなわち、信号が多くの段に分割され、分割後の信号の長さが約２０ｍｓ〜４０ｍｓとなる。各段はサウンドフレームと呼ばれる。これは規定の数値ではなく、使用者の需要に応じて微調整可能である。
隣接する２つのサウンドフレームの変化が過大になることを回避するため、隣接する２つのサウンドフレーム間に１０ｍｓの重複領域が設けられ、且つサウンドフレームの左右両端の連続性を強化するために各サウンドフレームに漢明ウィンドウ（ハミング窓関数）が乗算されることにより、サウンドフレーム内の信号の２つの境界辺がゆっくりと減衰し、境界上に明らかな不連続現象が生じなくなり、ノイズのエネルギースペクトルの強度が弱くなり、正弦波のピークが突出していることを示す。
フーリエ変換を行う際に、サウンドフレームの境界の左右に不連続の変化があった場合、フーリエ変換により原信号ではないエネルギー分布が生じ、分析判断に誤差が生じるため、信号に漢明ウィンドウ（ハミング窓関数）を乗算することでこの現象を低減する。
【００２８】
＜Ｇ23ステップ＞
前記高速フーリエ変換プロセス（Ｇ23）では、信号を時間領域から周波数領域に変換するために高速フーリエ変換（FFT）が用いられている。FFTは離散フーリエ変換（DFT）の高速アルゴリズムである。
【００２９】
＜Ｇ24ステップ＞
前記メルフィルターバンク（Ｇ24）では、フィルターバンクは互いに重複するバンドパスフィルターである。メル尺度に基づいて、１kHz未満では線形を呈し、それ以上では対数を呈する。メルスケーリングプロセスは式（５）を満たす。
【数５】
2020185390.tif 000010
ここでは、melはメルフィルターバンクの出力であり、ｆはフィルターの入力であり、2595及び700は固定値である。エネルギースペクトルは１組２０個の三角バンドパスフィルターが乗算されることで成り立ち、メル尺度をこの２０個のフィルターのスペクトルとする。
【００３０】
＜Ｇ25ステップ＞
前記離散コサイン変換プロセス（Ｇ25）では、MFCCの各フレームに対する計算に離散コサイン変換（DCT）が用いられている。DCTプロセスは式（６）を満たす。
【数６】
2020185390.tif 000011
前記メル周波数ケプストラム係数アルゴリズムを使用して前記いびき信号（Ｂ）を前記いびき機能マトリックス（Ａ）に変換する。
【００３１】
＜ステップｄ＞
ニューラルネットワークアルゴリズム（Ｇ3）により前記いびき機能マトリックス（Ａ）を分類し、前記いびき機能マトリックス（Ａ）が分類学習を経た後、前記いびき信号（Ｂ）の呼吸停止回数または呼吸低下回数を取得する。これにより、前記特徴抽出アルゴリズム（Ｇ2）により各前記単一のいびき信号ベクトル（Ｓ1〜Ｓ7）の特徴を抽出した後、前記いびき信号（Ｂ）から二次元の特徴ベクトルを取得する。
ディープラーニング関連技術はイメージ分類分野で大きな成功を収めており、多くのイメージ分類問題は全て前記ニューラルネットワークアルゴリズム（Ｇ3）を使用してイメージの特徴を抽出した後に分類を行っている。
前記ニューラルネットワークアルゴリズム（Ｇ3）によりイメージの特徴を抽出できるが、前記特徴抽出アルゴリズム（Ｇ2）により特徴の処理がなされた前記いびき機能マトリックス（Ａ）を検索できると考え、前記ニューラルネットワークアルゴリズム（Ｇ3）を使用して前記いびき機能マトリックス（Ａ）の分類を行ったが、この方法に限定するわけではない。
【００３２】
上述したように、本実施例では、前記ニューラルネットワークアルゴリズム（Ｇ3）は畳み込みニューラルネットワークアルゴリズムであり、前記畳み込みニューラルネットワークアルゴリズムはDenseブロックネットワーク（DenseNet）モデル（ＤＮ）が識別モデルとして選択されている。
図５に示されるように、前記Denseブロックネットワークモデル（ＤＮ）は複数のDenseブロック（Dense Block）（Ｄ）と、複数のTransitionレイヤー（Transition Layer）（Ｔ）と、分類レイヤー（Ｅ）と、を含む。前記Transitionレイヤー（Ｔ）は畳み込み（convolution）プロセス（Ｔ1)及びプーリング（pooling）プロセス（Ｔ2）を含み、且つ前記分類レイヤー（Ｅ）はSoftmaxレイヤーである。
図６に示されるように、前記Denseブロック（Ｄ）はDenseレイヤー（Ｉ）と、BR BN-ReLU-Convレイヤー（Batch Normalization-rectified linear units-Convolution）（ＢＲ）と、成長率（ｋ）と、を含み、前記成長率（ｋ）はレイヤーから出力された特徴写像の数量を示す。これにより、前記Denseブロックネットワークモデル（ＤＮ）は主に接続される各前記Denseブロック（Ｄ）及び各前記Transitionレイヤー（Ｔ）で構成され、最後に前記分類レイヤー（Ｅ）に接続される。ゆえに、前記Denseブロック（Ｄ）は密集して接続される畳み込みニューラルネットワークであり、前記Denseブロックネットワークモデル（ＤＮ）が構築されている。
その後、まず各いびき信号（Ｂ）を分割し、且ついびきが呼吸停止または呼吸低下であるかラベルを付ける。例えば、前記いびき機能マトリックス（Ａ）が呼吸停止または呼吸低下ではない場合、正常ラベル（Ａ1）を付け、前記いびき機能マトリックス（Ａ）が呼吸停止または呼吸低下である場合、病気ラベル（Ａ2）を付ける。
ラベル付けが完了した後、モデルにフィードして分類学習を行う。
最後に、識別睡眠時無呼吸モデル（Ｆ）が取得されるが、この方法に限定するわけではない。
【００３３】
さらに、各前記Denseブロック（Ｄ）中の任意の２つのレイヤー間を直接接続し、即ちネットワーク中の各レイヤーの入力を全ての先行レイヤーの全ての出力とし、且つ各レイヤーが学習した特徴写像が背後の全てのレイヤーに直接伝送されたものも全てこれらの入力となる。この方法により前記Denseブロックネットワークモデル（ＤＮ）が全ての特徴を更に有効的に利用し、各レイヤー間の伝送を強化する。
各前記Denseブロック（Ｄ）の間に挟まれた前記Transitionレイヤー（Ｔ）は特徴マトリックスのサイズを縮小するためのものであり、前記Denseブロック（Ｄ）の最後のレイヤーが累積した全ての先行レイヤーの情報を出力するため、モデルが非常に巨大となる。よって、前記Transitionレイヤー（Ｔ）が前記Denseブロック（Ｄ）及び前記Denseブロック（Ｄ）の間のサイズを縮小させることによってパラメータの数量が大幅に減少する。
上述の特性により、前記Denseブロックネットワークモデル（ＤＮ）が、ネットワークアーキテクチャが深すぎる場合に深いニューラルネットワークで直面する勾配の消失問題を解決する。モデルパラメータの数量が大幅に減少する以外、過剰適合に対する優れた抵抗力を更に有するが、但しこの方法に限定するわけではない。
【００３４】
このような構成により、図２に戻って、前記識別睡眠時無呼吸モデル（Ｆ）により正常信号（Ｆ1）または病気信号（Ｆ2）を予測する。
また、図７Ａに示されるように、前記識別睡眠時無呼吸モデル（Ｆ）のグラウンドトゥルース（Ground Truth）は、前記原信号（Ｙ）が青色曲線を呈することで構築された前記正常信号（Ｆ1）は正常ないびき（Normal snoring）を示すと共に緑色曲線を呈し、前記病気信号（Ｆ２）が閉塞性睡眠時無呼吸症候群（obstructive sleep apnea：OSA）を示すと共にピンク色の曲線を呈する。
また、図７Ｂ図に示されるように、前記識別睡眠時無呼吸モデル（Ｆ）に前記いびき信号（Ｂ）が入力されて赤色曲線を呈する場合に、前記正常信号（Ｆ1）及び前記病気信号（Ｆ2）のうちのどちらであるかを予測（Prediction）し、前記ニューラルネットワークアルゴリズム（Ｇ3）を利用していびきから睡眠時の呼吸停止を予測するが、この方法に限定するわけではない。
【００３５】
上記の実施形態は、本発明の理解を容易にするためのものであり、本発明を限定して解釈するためのものではない。本発明は、その趣旨を逸脱することなく、変更、改良され得ると共に、本発明にはその等価物が含まれることは言うまでもない。
【符号の説明】
【００３６】
ａ，ｂ，ｃ，ｄ      ステップ
Ｙ      原信号
Ｙｉ      フレームベクトル
Ｂ      いびき信号
Ｓ      いびき信号ベクトル
Ｍ      第一閾値
Ｘ      第二閾値
Ｘｉ      原信号の最大値
Ｍｉ      いびき信号の絶対値の合計
Ｒ      右側停止位置
Ｌ      左側停止位置
Ｗ      スライディングウィンドウ
Ｓ1      第一単一いびき信号ベクトル
Ｓ2      第二単一いびき信号ベクトル
Ｓ3      第三単一いびき信号ベクトル
Ｓ4      第四単一いびき信号ベクトル
Ｓ5      第五単一いびき信号ベクトル
Ｓ6      第六単一いびき信号ベクトル
Ｓ7      第七単一いびき信号ベクトル
Ａ      いびき機能マトリックス
Ａ1      正常ラベル
Ａ2      病気ラベル
Ｇ1      いびき信号分割アルゴリズム
Ｇ2      特徴抽出アルゴリズム
Ｇ21      プリエンファシスプロセス
Ｇ22      フレーミングとウィンドウプロセス
Ｇ23      高速フーリエ変換プロセス
Ｇ24      メルフィルターバンク
Ｇ25      非線形変換プロセス
Ｇ26      離散コサイン変換プロセス
Ｇ3      ニューラルネットワークアルゴリズム
ＤＮ      緻密ネットワークモジュール
Ｄ      Denseブロック
Ｔ      Transitionレイヤー
Ｔ1      畳み込みプロセス
Ｔ2      プーリングプロセス
Ｅ      分類レイヤー
Ｉ      Denseレイヤー
ＢＲ      BR BN-ReLU-Convレイヤー
ｋ      成長率
Ｆ      睡眠時無呼吸モデル
Ｆ1      正常信号
Ｆ2      病気信号
【図１】
2020185390.tif 000015
【図２】
2020185390.tif 000016
【図３Ａ】
2020185390.tif 000017
【図３Ｂ】
2020185390.tif 000018
【図３Ｃ】
2020185390.tif 000019
【図３Ｄ】
2020185390.tif 000020
【図３Ｅ】
2020185390.tif 000021
【図３Ｆ】
2020185390.tif 000022
【図３Ｇ】
2020185390.tif 000023
【図３Ｈ】
2020185390.tif 000024
【図３Ｉ】
2020185390.tif 000025
【図３Ｊ】
2020185390.tif 000026
【図３Ｋ】
2020185390.tif 000027
【図３Ｌ】
2020185390.tif 000028
【図４】
2020185390.tif 000029
【図５】
2020185390.tif 000030
【図６】
2020185390.tif 000031
【図７Ａ】
2020185390.tif 000032
【図７Ｂ】
2020185390.tif 000033