(57)【要約】加圧空気供給システム内の吸入ノイズを特性化する方法。当該方法は、吸入ノイズに基づいて吸入ノイズモデル（９１２，１０１２）を生成すること、少なくとも１つの吸入ノイズバーストを有する吸入ノイズを含む入力信号（８０２）を受信すること、前記入力信号を前記吸入ノイズモデルと比較（８１０）して類似性尺度を得ること、前記類似性尺度を少なくとも１つの閾値（８３２，８３４）と比較して前記少なくとも１つの吸入ノイズバーストを検出すること、および前記検出された少なくとも１つの吸入ノイズバーストを特性化すること、を備える。
2008504587.tif 000002
【特許請求の範囲】
【請求項１】
加圧空気供給システム内の吸入ノイズを特性化する方法であって、当該方法は、
加圧空気供給システムにより生成された吸入ノイズに基づき吸入ノイズモデルを生成することであって、吸入ノイズの少なくとも１つのデジタルサンプルから導出された少なくとも１つの自己相関シーケンスに基づき決定性デジタルフィルタとして表される前記吸入ノイズモデルを生成すること、
少なくとも１つの吸入ノイズバーストを有する吸入ノイズを含む入力信号を受け取ること、
前記入力信号を前記吸入ノイズモデルと比較して類似性尺度を得ること、
前記類似性尺度を少なくとも１つの閾値と比較して前記少なくとも１つの吸入ノイズバーストを検出すること、
前記検出された少なくとも１つの吸入ノイズバーストを特性化すること、
を備える、方法。
【請求項２】
前記特性化に基づく少なくとも１つのパラメータを判定することをさらに備える、請求項１記載の方法。
【請求項３】
前記少なくとも１つのパラメータは、使用者の呼吸数と吸入空気流量のうちの少なくとも１つを含む、請求項２記載の方法。
【請求項４】
前記検出された吸入ノイズバーストは、二値信号として特性化される、請求項１記載の方法。
【請求項５】
前記検出された吸入ノイズバーストは、
前記検出された少なくとも１つの吸入ノイズバーストの開始時間と、
前記検出された少なくとも１つの吸入ノイズバーストの終了時間と、
前記検出された少なくとも１つの吸入ノイズバーストの継続時間と、
を含む少なくとも１つの要因により特性化され、
前記少なくとも１つのパラメータは、前記特性化の要因のうちの少なくとも１つの関数である、請求項２記載の方法。
【請求項６】
前記決定性デジタルフィルタは、自己相関係数のセットから生成される線形予測符号化（ＬＰＣ）係数のセットに基づくＬＰＣフィルタであり、
前記吸入ノイズモデルを生成することは、
吸入ノイズをサンプリングして、前記吸入ノイズの前記少なくとも１つのデジタルサンプルを生成すること、
前記少なくとも１つのデジタルサンプルをウィンドウ化すること、
前記ウィンドウ化された少なくとも１つのデジタルサンプルから自己相関係数のセットを判定すること、
前記自己相関係数のセットに基づき前記ＬＰＣ係数のセットを生成すること、
ＬＰＣ係数のセットから前記ＬＰＣフィルタを生成すること
を備える、請求項１記載の方法。
【請求項７】
前記検出された少なくとも１つの吸入ノイズバーストに基づき前記吸入ノイズモデルを更新することをさらに備える、請求項１記載の方法。
【請求項８】
前記決定性デジタルフィルタは、第１の自己相関係数のセットに基づき生成される線形予測符号化（ＬＰＣ）係数のセットに基づくＬＰＣフィルタであり、
前記吸入ノイズモデルを更新することがさらに、
前記検出された少なくとも１つの吸入ノイズバーストをサンプリングして、該検出された少なくとも１つの吸入ノイズバーストの少なくとも１つのデジタルサンプルを生成すること、
前記検出された少なくとも１つの吸入ノイズバーストの前記少なくとも１つのデジタルサンプルから第２の自己相関係数のセットを判定すること、
前記第１及び前記第２の自己相関係数のセットの関数として前記第１の自己相関係数のセットを更新すること、
前記更新された自己相関係数のセットに基づいて前記ＬＰＣ係数のセットを更新すること、
前記更新されたＬＰＣ係数のセットに基づいて前記ＬＰＣフィルタを更新すること、
を含む、請求項７記載の方法。
【請求項９】
加圧空気供給システム内の吸入ノイズを特性化する装置であって、
処理素子と、
前記処理素子に接続された記憶素子であって、以下のステップ：
加圧空気供給システムにより生成された吸入ノイズに基づき吸入ノイズモデルを生成するステップであって、吸入ノイズの少なくとも１つのデジタルサンプルから導出された少なくとも１つの自己相関シーケンスに基づき決定性デジタルフィルタとして表される前記吸入ノイズモデルを生成するステップと、
少なくとも１つの吸入ノイズバーストを有する吸入ノイズを含む入力信号を受け取るステップと、
前記入力信号を前記吸入ノイズモデルと比較して類似性尺度を得るステップと、
前記類似性尺度を少なくとも１つの閾値と比較して前記少なくとも１つの吸入ノイズバーストを検出するステップと、
前記検出された少なくとも１つの吸入ノイズバーストを特性化するステップと、
を実行するように処理装置に指示するためのコンピュータプログラムを記憶する、前記記憶素子と、
を備える装置。
【請求項１０】
吸入ノイズを特性化するシステムであって、
加圧空気供給システムと、
前記加圧空気供給システムに接続された通信システムであって、当該通信システムが、
前記加圧空気供給システムにより生成された吸入ノイズに基づき吸入ノイズモデルを生成することであって、吸入ノイズの少なくとも１つのデジタルサンプルから導出された少なくとも１つの自己相関シーケンスに基づき決定性デジタルフィルタとして表される前記吸入ノイズモデルを生成すること、
少なくとも１つの吸入ノイズバーストを有する吸入ノイズを含む入力信号を受け取ること、
前記入力信号を前記吸入ノイズモデルと比較して類似性尺度を得ること、
前記類似性尺度を少なくとも１つの閾値と比較して前記少なくとも１つの吸入ノイズバーストを検出すること、
前記検出された少なくとも１つの吸入ノイズバーストを特性化すること、
を含む、通信システムと、
を備えることを特徴とするシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、一般的に通信システムに連結される加圧空気供給システムに関する。
【背景技術】
【０００２】
消火のような危険な環境下での活動に従事する人員間の適切かつ確実な通信は、彼ら自身の健康と安全を維持しつつ職務を達成するために不可欠である。作業状況によっては、たとえば自給式呼吸器（ＳＣＢＡ）マスクと空気供給システム、自給式潜水呼吸器（ＳＣＵＢＡ）マスクと空気供給システム、または航空機酸素マスクシステムなどの加圧空気供給システムの使用を要する場合がある。しかしながら、人員が上記加圧空気供給システムを使用している間でも、適切かつ確実な通信が維持され、人員の健康と安全が有効に監視されることが望ましい。
【０００３】
図１は、通信システム１３０に連結される加圧空気供給システム１１０を含む従来技術システム１００の概略ブロック図である。加圧空気供給システムは従来、ＳＣＢＡマスクなどの呼吸マスク１１２、エアシリンダ（図示せず）、レギュレータ１１８、レギュレータ１１８をエアシリンダに接続する高圧ホース１２０を有する。使用される空気供給システム１１０の種類に応じて、システム１１０はたとえば、清潔な呼吸空気を使用者に提供したり、使用者の肺に有害毒素が届かないようにしたり、燃焼する構造物内の過熱空気によって使用者の肺が焼けるのを防いだり、使用者の肺を水から守ったり、顔や呼吸器官の火傷から使用者を守ったりすることによって、使用者を保護する。さらに、概して、空気はマスクの着用者が吸入するときにだけ通常供給されるため、マスクはプレッシャーデマンド型呼吸システムとみなされる。
【０００４】
通信システム１３０は通常、マスク着用者の音声を記録するように設計され、マスクの内部に搭載されるか、マスクの外部に装着されるか、あるいはマスク１１２のボイスミッターポート上で手で保持し得る従来のマイクロフォン１３２を有する。さらに、通信システム１３０は、マスク着用者が音声をたとえば他の通信部に伝えるために使用可能な双方向無線機などの通信部１３４も有する。マスクマイクロフォン装置１３２は、無線機１３４に直接、あるいは仲介電子処理装置１３８を介して接続される。この接続は従来のワイヤケーブル（たとえば１３６）を通じて行うことも、従来のＲＦ、赤外線、または超音波短距離送受信システムを用いて無線で行うこともできる。仲介電子処理装置１３８はたとえば、デジタル信号処理装置として実装され、インタフェース電子装置、オーディオアンプ、および仲介電子処理装置用ならびにマスクマイクロフォン用のバッテリー電源を含んでもよい。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
システム１００などのシステムの使用に付随する欠点がいくつかある。説明の簡易化のために、図１に示されるシステム１００のマスク−無線機音声経路を示す図２のブロック図を参照し、これらの制限を説明する。唇からの音声入力２１０（たとえば、Ｓｉ（ｆ））が、音響共鳴およびヌルを特徴とする音響伝達関数２２０（たとえば、ＭＳＫ（ｆ））を有するマスク（たとえば、ＳＣＢＡマスク）に入る。これらの共鳴とヌルは、マスクの空洞体積とマスク内表面からの音の反響による。伝達関数ＭＳＫ（ｆ）を特徴とするこれらの効果は、入力された音声波形Ｓｉ（ｆ）を歪ませ、そのスペクトル成分を変化させる。もう１つの音源は呼吸器具から生成されたノイズ２３０（たとえば、レギュレータ吸入ノイズ）で、マスクにも入り、ＭＳＫ（ｆ）の影響を受ける。もう１つの伝達関数２４０（たとえば、ＮＰｋ（ｆ））は、ノイズが音声の位置とわずかに異なるマスク内の位置から生成されるという事実の説明となる。音声とノイズＳ（ｆ）は、固有の伝達関数２５０（たとえば、ＭＩＣ（ｆ））を有するマイクロフォンにより音響エネルギーから電気信号に変換される。その後、マイクロフォン信号は通常、オーディオアンプ、および伝達関数２６０（たとえば、ＭＡＡ（ｆ））を有するその他の回路を通過する。次に、ＭＡＡ（ｆ）からの出力信号２７０（たとえば、Ｓｏ（ｆ））は、さらなる処理と伝送のために無線機に入力される。
【０００６】
システム１００のようなシステムの欠点に戻ると、このような欠点の例は、動作の１部として大きな音響ノイズがこうしたシステムにより生成されることに関連する。より具体的には、これらのノイズは、特に無線機などの電子システムとともに使用されるとき、通信の質を大幅に低下させる可能性がある。ＳＣＢＡシステムのような加圧空気供給システムによって生じる顕著な音響的アーティファクトであるノイズは、図２で枠２３０として示されるレギュレータ吸入ノイズである。
【０００７】
レギュレータ吸入ノイズは、マスク着用者が吸入する度に生じる広帯域ノイズバーストとして発生する。マスク内が負圧になると、エアレギュレータバルブが開放され、高圧空気がマスクに入り、大きなヒス音を引き起こす。このノイズがその後に続く音声と共にマスク通信システムのマイクロフォンに補捉され、音声とはほぼ同じエネルギーを有する。吸入ノイズは通常吸入後にのみ発生するので、ふつうは音声にはかぶらない。しかしながら、この吸入ノイズが問題を引き起こすこともあり、その例を以下に示す。たとえば、吸入ノイズはＶＯＸ（音声動作スイッチ）回路を始動させることにより、無線チャンネルを開放し占有して、同じ無線チャンネル上の他の話し手を妨害する可能性がある。さらに、デジタル無線機を使用する通信システムでは、吸入ノイズがＶＡＤ（音声活動検出器）アルゴリズムを始動させて、無線信号処理チェーンのさらに下流でノイズ抑制アルゴリズムにおいてノイズ推定に混乱をきたす場合がある。また、概して吸入ノイズは聞き手を困らせる。
【０００８】
システム１００などのシステムの第２の欠点を以下に述べる。これらのシステムは通常、鼻、口、または顔全体を覆うマスクを使用する。空気システムのマスクは、マスクの体積と内反射面の形状寸法の関数である特定の音響共鳴および反共振（ヌル）のセットを呈示し、マスク内で生成される音声のスペクトル特性を変える固定形状寸法の閉鎖された空気腔を形成する。より具体的には、空気マスク音声経路（図２）を特性化する際、システムの最も難しい部分は、話し手の唇からマスクマイクロフォンへの音響伝達関数（２２０）である。これらのスペクトル歪は、装着される音声通信システム、特に劣化した音声を処理するように最適化されていないパラメータデジタル符復号器を用いるシステムの性能を大幅に低下させかねない。音響マスク歪は、特にパラメータデジタル符復号器が含まれる場合、通信システムの質と了解度に影響を及ぼすことが証明されている。一般的には、吸入ノイズだけでなく、空気システムはマスクの不十分な音響により、発生する音声の質を最大限悪化させる。
【０００９】
図３は、マスク内（３２０）およびマスクマイクロフォン出力（３１０）で測定されたスペクトル振幅反応と、マスク、マイクロフォン、およびマイクロフォン増幅器に関して算出された総合伝達関数（３３０）の例を示す。これらの具体的なデータは、頭と胴のシミュレータに搭載されたＳＣＢＡマスクを用いて得られた。音響の励起は人工マウスシミュレータを駆動する３Ｈｚ〜１０ＫＨｚの掃引正弦波からなる。図３に示されるように、スペクトルは主にマイクロフォン内のプリアンプ帯域通過フィルタのために、５００Ｈｚ未満および４．０ＫＨｚ超の周波数で大幅に減衰され、５０〜４．０ＫＨｚの有意の音声通過帯域領域に多数の顕著なスペクトルの山と谷を含む。これらのスペクトルの山と谷は通常、くし型フィルタ（ｃｏｍｂ  ｆｉｌｔｅｒｉｎｇ）を生じるマスク内部の反射と空洞共振状態とによって引き起こされる。有意なスペクトルの山と谷は、通過帯域を前後に移動しつつ、音声ピッチの成分とフォルマントを変調し、結果的に質の悪化と音声の歪みを招く。このようなシステムを特性化する１つまたは複数の伝達関数を決定することが望ましい。その伝達関数は音声歪を低減する均等化システムを定義するために使用される。
【００１０】
システムの伝達関数を適切に決定し伝送チャンネルを均等化するためには、実績ある手法が多数存在する。システムの伝達関数を決定する有効な方法の１つが、システムを励起し、システムパラメータを決定する広帯域基準信号の使用である。多くの音声伝達環境の伝達関数を推定する際の問題は、適切な広帯域励起信号が容易に利用可能でないことである。一般的なアプローチの１つは、基準として長期平均音声スペクトルを使用することである。しかしながら、特に音声入力がまばらであるとき、この基準を使用する順応時間が長くなる可能性がある。また、長期音声スペクトルは、音声スペクトルをかなり変えるおそれのある叫び声や感情的ストレスを含むことの多い公共サービス活動においては、個人間で相当な変動がある。
【００１１】
システム１００などのシステムに付随する別の欠点は、たとえば、生体計測パラメータを含むマスク着用者の特定のパラメータを測定するより有効な方法と装置がないことである。危険な環境で働き、システム１００のようなシステムを使用しているかもしれない個人の上記パラメータの測定は、彼らの安全と動作を監視するために重要である。たとえば、個人の呼吸数と空気消費量の測定は、その個人の作業量、生理的適合性、ストレスレベル、および保管された供給空気の消費量（すなわち、有効作業時間）を特徴付ける重要なパラメータである。呼吸を測定する従来の方法には、胸部インピーダンスプレチスモグラフ法の使用または熱感度器を使用する気流温度測定などがある。しかしながら、消火活動などの肉体的に厳しい環境で作業する個人から、これらの従来の方法を用いて確実な測定値を得ることは、通常測定のために使用される身体装着センサおよびアーティファクトのずれを引き起こす可能性のある激しい肉体運動のためにより困難である。
【課題を解決するための手段】
【００１２】
よって、吸入ノイズを有効に検出および減衰し、音声を均等化し（すなわち、歪み効果を除去し）、通信システムに連結される加圧空気供給システムを含むシステムの使用者に関連するパラメータを測定するための方法と装置が必要とされる。
【発明を実施するための最良の形態】
【００１３】
以下、添付の図面を参照して、本発明の好適な実施形態を説明する。本発明は多くの様々な形を取る実施形態の代表的なものであり、図面で示し、本明細書で特定の実施形態を詳細に説明する。本記載は本発明の原則の例としてみなされ、図示され説明される特定の実施形態に本発明を制限することを意図していないことを理解し得る。さらに、本明細書で使用される用語および単語は限定的ではなく、単に説明的なものとみなされるべきである。説明の簡潔性と明瞭性のために、図面に示される構成要素は必ずしも一律の縮尺に従って描かれているとは限らない。たとえば、一部の構成要素の寸法は互いに誇大に描かれている。また、適切と考えられる場合、対応する要素を示すために、各図面間で参照符号が繰り返される。
【００１４】
本発明の様々な側面を詳細に説明する前に、先に略述したエアレギュレータ吸入ノイズをより詳細に説明することが本発明を理解する上で有益である。吸入ノイズは、人が息を吸い込み、レギュレータバルブが開放されるときにＳＣＢＡまたはその他の加圧空気供給システムマスクに入る高圧空気の結果である。バルブでの乱流は、ＳＣＢＡマスクに直接結合される非常に大きな広帯域ヒス雑音を引き起こし、それは音声信号を有するマイクロフォンでの振幅に匹敵する。ＳＣＢＡマスクおよびその広帯域スペクトグラム５００内に記録される典型的な吸入ノイズ４００の例をそれぞれ図４および５に示す。
【００１５】
図５に示されるように、ノイズスペクトルは、およそ５００、１７００、２７００、および６０００Ｈｚで発生するスペクトルの顕著な山を有する広帯域である。山はマスク内の共振と、マスク内部の反射によるくし型フィルタから起こり、各種マスクモデル、サイズ、および構造に応じて周波数や大きさが変動する場合がある。いったんマスクが顔に装着されると、内部の全体的な幾何学的配置はほぼ一定であるため、ノイズスペクトルの配色は、マスク／着用者の特定の組み合わせに関して通常固定される。このことは、所与のＳＣＢＡマスクを着用する同一の話し手に関して、ＳＣＢＡマスクのマイクロフォンから生じる３つの別々の吸入ノイズ（６１０、６２０および６３０）が重畳して示される図６に実証される。この一致は、異なる話し手および異なるメーカー製のマスクに関しても観察された。さらに、同一のマスクを着用する、異なる話し手から生じるエアレギュレータノイズのスペクトルの高い類似性も観察された。
【００１６】
最後に、図７はＳＣＢＡシステムから記録される音声７１０の例を示す。図７が実証するように、人は通常息を吸い込む間、話そうとはしないため、吸入ノイズ７２０の影響は音声自体には及ぼされない。しかしながら、ノイズは、音声検出器と無線機におけるノイズ抑制回路とに関わる問題を引き起こし、聴音しづらくするのに十分なエネルギーとスペクトルとを有する。
【００１７】
本発明の第１の態様は、図１に示されるシステム１００のような通信システムに連結される加圧空気供給システムにおける吸入ノイズを検出し除去する方法と装置である。本発明のこの実施形態による方法は、本明細書ではＡＲＩＮＡ（エアレギュレータ吸入ノイズ減衰器）とも称される。エアレギュレータの吸入ノイズを特定し除去するＡＲＩＮＡ法の根拠は、音声や、たとえば様々な環境ノイズなどのその他の種類のノイズと比べた場合の相対的定常性である。ＡＲＩＮＡ法８００のブロック図は図８に示され、４つの部分、ノイズモデルマッチング部８１０、ノイズ検出部８３０、ノイズ減衰部８５０、およびノイズモデル更新部８７０に分割することができる。
【００１８】
ＡＲＩＮＡ法８００の基本的な方法論は以下のように要約される。方法８００は、好ましくはデジタルフィルタ（たとえば、全極型線形予測符号化（ＬＰＣ）デジタルフィルタ）を用いて吸入ノイズをシミュレーションする。次に、方法８００は、ノイズモデルフィルタの逆関数を用いる音声入力信号（すなわち、マスクマイクロフォンにより補捉される音声およびノイズ）を濾波し、逆ノイズモデルフィルタの出力エネルギーと、入力信号またはその他のエネルギー基準のエネルギーとを比較する。接近したスペクトル一致が入力信号とモデルとの間で生じる信号期間中、入力信号を備えるレギュレータ吸入ノイズは所望レベルに減衰させることができる。
【００１９】
図８に示されるＡＲＩＮＡ法８００の特性に移ると、処理の第１ステップは、方法８００のノイズモデルマッチング部８１０を介して入力信号８０２を基準ノイズモデルと継続的に比較することにより吸入ノイズの発生を検出することである。これは、許容可能な実行の複雑さに応じて図９または図１０に従い好適な実施形態において実行することができる。しかしながら、当業者であれば、代替のスペクトルマッチング方法を利用可能なことを認識し得る。図９および図１０に示される２つの好適なマッチング方法は、本明細書では正規化モデル誤差（Ｎｏｒｍａｌｉｚｅｄ  Ｍｏｄｅｌ  Ｅｒｒｏｒ、すなわちＮＭＥ）法およびイタクラ−サイトウ（Ｉｔａｋｕｒａ−Ｓａｉｔｏ、すなわちＩ−Ｓ）歪み法と称される。両方法では、基準ノイズモデルが、吸入ノイズのスペクトル特徴を概算するデジタルフィルタ（９１２、１０１２）により表される。好適な実施形態では、このモデルは、ＬＰＣ係数のセットにより特定される全極型（自己回帰）フィルタとして表される。しかしながら、当業者であれば、たとえば、既知のＡＲＭＡ（自己回帰移動平均）モデルなどの全極型モデルの代わりに代替のフィルタモデルが利用可能なことを認識し得る。
【００２０】
基準ノイズモデルフィルタ係数は、吸入ノイズの少なくとも１つのデジタルサンプル（ｄｉｇｉｔｉｚｅｄ  ｓａｍｐｌｅ）から導いた自己相関係数のセットから得られる。最初のノイズサンプルとそれに対応する最初の自己相関係数（８７２）は、任意数の前もって取られたノイズ録音からオフラインで取得することができ、本発明の実施に必須ではない。さらに、実験が示すように、あるＳＣＢＡマスクからの最初のノイズサンプルは、たとえば同じ設計の他のマスクに関しても、また場合によっては異なる設計のマスクに関しても適切に働く。自己相関係数は、未処理のサンプルノイズデータから直接算出することもできるし、あるいは当業者にとって十分既知な共通の方法を用いてＬＰＣまたはその他の反射係数などのその他一般的に使用されるスペクトルパラメータ表示から導き出すこともできる。
【００２１】
好適な実施形態では、ノイズモデルの自己相関係数は、以下の標準式に従い算出される。
【００２２】
【数１】
2008504587.tif 000003
ただし、Ｒｉは最大ｐ個の自己相関係数のｉ次の係数、ｘｎは最大Ｎ個のサンプルがある標準的な吸入ノイズ信号サンプルセグメントのｎ次のサンプル、Ｒ０はセグメント全体のエネルギーを表す。自己相関係数ｐの次数は通常１０〜２０で、好適な実施形態での値は１４である。さらに、理想的には、Ｎ個の信号サンプルは、スペクトル推定を円滑化するために自己相関が実行される前に、ハミング窓（Ｈａｍｍｉｎｇ  ｗｉｎｄｏｗ）を用いてウィンドウ化される。ハミング窓は下式で表される。
【００２３】
【数２】
2008504587.tif 000004
当業者であれば、その他のウィンドウ化方法も利用できることを認識し得る。
【００２４】
次に、ノイズモデル自己相関係数を使用して、以下のｚ−領域表示伝達関数を有する全極型線形予測モデルフィルタを表す、１０次のノイズモデルＬＰＣ係数ａ１、ａ２、．．．、ａｐのセットが決定される。
【００２５】
【数３】
2008504587.tif 000005
ただし、ｚ＝ｅ−ｊｎωＴはｚ変換変数である。この例では、１０次のＬＰＣ係数が決定された。しかしながら、特定の実施に基づき、別の次数のＬＰＣ係数を選択することもできる。自己相関−ＬＰＣパラメータ変換（ステップ９１２、１０１２）は、当業者にとって既知な任意数のパラメータ変換手法を用いて実行することができる。好適な実施形態では、ＬＰＣパラメータは、当業者にとって十分既知なダービン（Ｄｕｒｂｉｎ）方法を用いて自己相関パラメータから導かれる。
【００２６】
次に、図９に示されるＮＭＥスペクトルマッチング方法の詳細に移ると、導出された全極型ＬＰＣノイズフィルタモデルは、逆ＬＰＣフィルタ（ステップ９１４）を形成するように逆転される。
【００２７】
【数４】
2008504587.tif 000006
理想的には、マスクマイクロフォンから得られ且つ音声および吸入ノイズを含む、低域濾波されサンプリングされた音声入力信号８０２、Ｓ（ｚ）が逆フィルタＨ（ｚ）を通過することにより（ステップ９１４）、出力信号が得られる。
【００２８】
【数５】
2008504587.tif 000007
次に、逆フィルタの入力および出力信号のエネルギーＥｉｎ、Ｅｏｕｔが（それぞれステップ９１８および９１６で）算出され、歪み尺度Ｄはステップ９２０で算出される。歪み尺度Ｄは、ノイズモデルと入力信号との間の類似性尺度として機能する。Ｄの理論上の下限は無限次元ではゼロだが、実質上、下限は、入力信号とそれが有限次元のノイズモデルにどのくらい適切にマッチするかとによって決定される。この実施では、歪み尺度は正規化モデル誤差（ＮＭＥ)と称され、Ｅｉｎに対するＥｏｕｔの比により定義され、ステップ９２０で次式のように算出される。
【００２９】
【数６】
2008504587.tif 000008
次に入力信号のエネルギーは、それがノイズモデルにどのくらい適切にマッチするかに従い除去される。好適な実施形態では、上述の信号フィルタリングは、先行する式で示されるように周波数領域で行うこともできるが、時間領域における畳み込みを介して行われる。
【００３０】
ＡＲＩＮＡ法８００の信号処理は通常、セグメント化フレーム毎に行われる。好適な実施形態では、入力信号８０２は低域濾波され、８．０ＫＨｚでサンプリングされ、８０サンプルのブロックにバッファされ（１０ｍｓｅｃ）、逆ノイズモデルフィルタを通過する（式５）。よって、すべてのフィルタリングは理想的には、入力信号８０２の連続８０サンプルセグメントで行われる。次に、逆ノイズモデルフィルタの正規化モデル誤差（ＮＭＥ）が、フィルタ出力フレームエネルギーを入力信号フレームエネルギーで割ることにより算出される（式６）。しかしながら、この計算は、時間分解能を向上させるため理想的には、サブクレーム毎に行われる。したがって、各８０ポイントフレームがたとえば４つの２０ポイントサブフレームに分割されるが、代替のサブフレームへの分割は必要とされる精度に応じて利用することができる。次に、最後の１６のサブフレームの出力フィルタエネルギーＥｏｕｔの平均を取り、その量を対応するタイムアラインされた１６のサブフレームの入力フィルタエネルギーＥｉｎの平均で割ることにより、正規化モデル誤差信号（ＮＭＥ）全体を均等化することができる。これは分析を遅延させず、レギュレータノイズスペクトルを変質させる場合のあるその他の大きな背景ノイズの一時的なドロップアウトや影響を除去するのに役立つ。それによって、本発明の本実施では、平均ＮＭＥ値は、入力信号スペクトルの類似性に対するノイズモデルの尺度として使用される。
【００３１】
好適な実施形態では、図１０に示される第２の、より複雑だがより正確なノイズモデルマッチング方法８１０は、イタクラ−サイトウ歪み法の変形である。２つの信号間のスペクトル類似性を判定するＩ−Ｓ法は、当業者にとって周知である。この方法では、残りのノイズモデル逆フィルタエネルギーが、前述のＮＭＥ法のように入力信号エネルギーとではなく、「最適な」信号フィルタの残留エネルギーと比較される。フィルタは、現在の信号セグメントのスペクトルと最もよくマッチするという意味において「最適」である。
【００３２】
最適に濾波された信号に対応する残留エネルギーは、ステップ１０１８〜１０２４を用いて算出される。Ｉ−Ｓ法では、ステップ１０１８で、理想的には２つの連続する入力信号８０２の８０のサンプルバッファが、単独の１６０のサンプルセグメントに結合される。１６０のサンプルセグメントは、好ましくは次式により与えられる１６０ポイントのハミング窓を用いてウィンドウ化される。
【００３３】
【数７】
2008504587.tif 000009
ウィンドウ化された信号データは次に、式１に記載される方法を用いて自己相関させられる。ステップ１０１８で生成される自己相関係数は、Ｒｉ、ｉ＝０、１、２、．．．、ｐとして指定される。ＬＰＣ係数の対応セットは、好ましくはステップ１０１２で基準ノイズモデルパラメータを生成するために使用されたのと同じように、ステップ１０２０でダービンアルゴリズムを用いて自己相関係数から導かれる。ステップ１０２０で生成される信号モデルＬＰＣ係数は、ａｉｉ＝１、２、．．．、ｐとして指定される。ステップ１０２２では、これらのＬＰＣ係数（ステップ１０２０）が、ｂｉを導き出す下記式９に従い自己相関させられる。これらのパラメータを用いて、このフィルタを通過する信号Ｅｓの残留エネルギーが次式によりステップ１０２４で算出される。
【００３４】
【数８】
2008504587.tif 000010
【００３５】
【数９】
2008504587.tif 000011
ノイズモデルを通過する入力信号のエネルギーは、ステップ１０１２〜１０１６を用いて算出される。ステップ１０１２では、ノイズモデルＬＰＣ係数が上述するようにノイズモデル自己相関係数（８７４）から算出される。ステップ１０１２で生成されるこれらのＬＰＣ係数は、ａｉｉ＝１、２、．．．、ｐとして指定される。ステップ１０１４では、ＬＰＣ係数（ステップ１０１２から）が、ｂｉを導き出す下記式１１により自己相関させられる。これらのパラメータとステップ１０１８で算出される自己相関シーケンス、Ｒｉを用いて、基準ノイズモデルを通過する信号のエネルギーは式１０により与えられるようにステップ１０１６で算出される。
【００３６】
【数１０】
2008504587.tif 000012
【００３７】
【数１１】
2008504587.tif 000013
基準ノイズモデルに対する「最適」信号モデルのスペクトル歪の尺度Ｄは、次式で定義されるようにステップ１０２８で算出される。
【００３８】
【数１２】
2008504587.tif 000014
信号モデルが基準ノイズモデルに類似するほど、歪み尺度は下限の１．０に近づいていく。この歪み尺度は、吸入ノイズの存在を判定するために、ＡＲＩＮＡ法８００のノイズ検出部８３０によって使用される。Ｉ−Ｓの歪み尺度は、好適な実施形態では１６０のサンプルを用いて算出される。Ｉ−Ｓの歪み尺度によって判定されるような吸入ノイズの分類は、１６０サンプルセグメントの各８０サンプルフレームに関連付けられる。さらに、ステップ１０１２および１０１４は、最初のノイズモデル（たとえば、最初の自己相関係数８７２に基づき）を生成する、あるいは上述され、以下詳述されるノイズモデル更新部８７０に従いノイズモデルを更新するために行うだけでよい。
【００３９】
ＡＲＩＮＡ法８００のノイズ検出部８３０では、スペクトルマッチ８１０から導かれる値（すなわち、入力信号とノイズモデル間の類似性尺度を表すＮＭＥまたはＩ−Ｓの歪み尺度）は次に、実験的に導かれる閾値（たとえば、Ｄｍｉｎ１)と比較される。この検出閾値は、音声や他の種類のノイズを誤って吸入ノイズと分類せずに吸入ノイズの存在を検出するために選択される。
【００４０】
さらに、ノイズフィルタモデルの特定性、吸入ノイズのスペクトル変動、および音声とノイズモデルの類似性に応じて、たとえば、誤検出が起こる可能性がある。したがって、実際のエアレギュレータ吸入ノイズの継続時間が音声アーティファクトと比べてかなり長いため、ノイズ継続時間閾値テストも理想的には適用される（ステップ８３４)。よって、検出閾値は、検出が有効と確認される前に、所定数の連続フレーム「Ｋ１」（たとえば、４フレーム）に対して満たされねばならない。相対的な信号エネルギー、波形のゼロ交差、およびその他の特徴的なパラメータ情報を、音声／吸入ノイズの区別を向上させるために検出スキームに含めてもよい。よって、両方の閾値基準が満たされる場合（ステップ８３２および８３４)、スペクトルマッチは許容可能に近似しているとみなされ、吸入ノイズが現在存在すると推定される。
【００４１】
ＡＲＩＮＡ法８００のノイズ減衰部８５０では、ノイズ検出部８３０の出力が、出力信号乗算器（８５２）をゲートするために使用され、この出力信号乗算器（８５２）を入力信号８０２が通過する。吸入ノイズが検出されれば、乗算器ゲインＧがステップ８５４で所望の減衰値「Ｇｍｉｎ」に設定される。この減衰ゲイン値は完全にノイズを除去するため０．０であってもよいし、ノイズを完全には除去しないが抑制するためにより高い値に設定してもよい。エアレギュレータが機能していることを聞き手に納得させるために、完全な抑制は望ましくないかもしれない。好適な実施形態では、Ｇｍｉｎは０.０５の値を有する。そうでなく、吸入ノイズが検出されない場合、ゲインＧは音声信号を減衰しないように理想的には１．０に設定される。このゲーティング／乗算スキームは種々に変形することができる。たとえば、ゲーティングのアタック減衰が突然に起こらないように、吸入ノイズの前後に直接生じる場合がある音声を減衰する可能性を減少させることによって、感知する音声の質を向上させるような変形を採用することができる。さらに、方法８００から容易に分かるように、本発明の重要な利点は、レギュレータノイズが検出されるときを除き、従来の継続的なノイズフィルタリング方法と異なり、元の信号が変質させられないことである。
【００４２】
ＡＲＩＮＡ法８００の重要な要素は、検出のためにノイズモデルを定期的に更新可能なことである。たとえば、時間が経過すると、顔の上の空気マスクの移動によって、空気マスクが音響伝達関数に及ぼす影響に変化が生じる場合がある。また、様々な人々による空気マスクの着用または異なるマスクの使用は、最初の基準ノイズモデルのスペクトルが実際の吸入ノイズスペクトルから逸脱する場合があることを意味する。元の基準ノイズモデルを定期的に更新することにより、正確な現在の基準ノイズモデルを維持できる。したがって、ＡＲＩＮＡ法８００のノイズモデル更新部８７０はノイズモデルの更新に使用される。
【００４３】
ノイズモデル更新部８７０は、レギュレータ吸入ノイズの基準ＬＰＣフィルタモデルを更新すべき時を判定するためにノイズ検出部８３０の出力を使用する。たとえば、ノイズ検出部８３０からの出力は、ノイズモデルを更新するかどうか判定するステップ８７６で、第２の実験的に決定される閾値（たとえば、Ｄｍｉｎ２）と比べることができる。閾値が満たされると、吸入ノイズとして検出される連続サブフレームの数がカウントされて（ステップ８７８)、各サブフレーム内の信号サンプルがバッファに記憶される。連続するノイズサブフレームの数が閾値「Ｋ２」（たとえば、好適な実施形態では、８サブフレーム、１６０サンプル）を超えると、ステップ８８０でノイズモデルを更新する決定がなされる。非ノイズサブフレームが検出されると（たとえば、ステップ８３２、８３４および８７６のいずれかで)、ノイズフレームのカウント値がステップ８８４でゼロにリセットされ、ステップ８７８で更新される。次に、現在検出されている吸入ノイズを表す「Ｋ２」個の連続信号サブフレームの自己相関係数は、前述の式１および式２を用いてステップ８８２で算出することができる。
【００４４】
これらの新たな自己相関係数は、ステップ８７４でノイズモデル自己相関係数を更新するために使用される。理想的には、ステップ８８２で算出される自己相関係数は、たとえば次式のような単純な重み付け式を用いてステップ８７４で、先のノイズモデル自己相関係数で平均を取る。
【００４５】
【数１３】
2008504587.tif 000015
ただし、ＲｉＲＥＦは現在の基準ノイズモデルの自己相関係数で、ＲｉＮＥＷは現在検出されている吸入ノイズサンプルの自己相関係数で、αは最初の基準モデルがどのくらい速く更新されるかを決定する１．０〜０．０の重み係数である。この重み係数は、吸入ノイズのスペクトル特性の変動の速さに応じて調節することができ、前述したようにその速さは通常遅い。次に、ノイズモデル逆フィルタのＬＰＣ係数の新しいセットは、ステップ９１２および１０１２で更新されたモデル自己相関から算出することができる。ノイズモデルからの大きなずれが誤検出により生じないように、ノイズモデルの調節に制限を課すことができる。また、必要に応じてシステムが初期モデル状態にリセット可能なように最初の基準ノイズモデル係数（８７２）が記憶される。ノイズモデル更新部８７０を参照することにより、上述の方法８００の調整性能で、システムは特定のマスクおよびレギュレータの特性を調整し、最適な検出性能を実現することができる。
【００４６】
ＡＲＩＮＡ法８００の利点として、従来の継続的フィルタリングを採用するアルゴリズムの場合のように、アルゴリズムを処理することによって音声信号自体が不可逆的に影響を及ぼされることがない。別の利点としては、ここで使用されるＬＰＣモデリングは、簡単で、リアルタイムで容易に調整可能であり、単純でコンピュータでの計算が効率的である。当業者であれば、上記の利点が本発明のＡＲＩＮＡの実施形態に関連するすべての利点を包含することを意味するのではなく、その代表的なものとしての機能を果たすことを意味しているにすぎないことを認識し得る。
【００４７】
本発明の第２の態様は、図１に示されるシステム１００のような通信システムに連結される加圧空気供給システムで音声信号を均等化する方法と装置である。本発明の本実施形態による方法は、本明細書ではＡＭＳＥ（空気マスク音声イコライザ）法とも称する。均等化のＡＭＳＥ法の根拠は、音声やたとえば様々な環境ノイズなどの他の種類のノイズと比べた場合のノイズの相対的定常性である。同じマスクの共鳴状況がレギュレータノイズと音声信号の両方に影響を及ぼすため、音響反射によるピークとヌルは音声とノイズとの間のソース位置の違いによりノイズと音声との間でわずかに異なるが、ノイズを均等化することは、音声信号の均等化に適したイコライザをもたらすことにもなるはずである。
【００４８】
ＡＭＳＥ法は、マスクの空洞と構造により生成される音響共鳴スペクトルのピークとヌル（すなわち、スペクトル規模音響伝達関数）を推定するため、すべてのマスク式加圧空気呼吸システム（たとえば、ＳＣＢＡ）に存在する広帯域エアレギュレータ吸入ノイズを利用する。次に、このスペクトル情報が、リアルタイムでデジタル逆フィルタの補正を構成するために使用される。デジタル逆フィルタはスペクトル面で歪んだ音声信号を均等化し、マスクなしに生成される非歪み音声を概算する出力信号を生成するのに適用される。この動作によって、マスクマイクロフォンから得られる音声の質が向上し、通信の了解度を向上させることができる。
【００４９】
ＡＭＳＥ法の詳細に移ると、方法１１００のブロック図を図１１に示す。その方法はノイズモデルマッチング部１１１０、ノイズ検出部１１３０、マスク音声均等化部１１５０、およびノイズモデル更新部１１７０の４つの部分に分割することができる。ＡＭＳＥ法のノイズモデルマッチング部、ノイズ検出部、およびノイズモデル更新部は理想的には、先に詳述したＡＲＩＮＡ法の対応部分と同一である。したがって、簡潔さのため、これら３つの部分の詳細な説明はここでは繰り返さない。しかしながら、以下、ＡＭＳＥ法１１００のマスク音声均等化部１１５０（点線領域）を詳述する。
【００５０】
ＡＭＳＥ法１１００のマスク音声均等化部１１５０を用いて、上述の式３を用いるステップ１１５２でノイズのｎ次ＬＰＣモデルを生成するため、吸入ノイズ基準自己相関係数が使用される。ステップ１１５２で生成されるＬＰＣモデルはマスクの伝達関数、たとえば、図２のＭＳＫ（ｆ）を特性化し、吸入ノイズに関してはノイズ経路伝達関数ＮＰ（ｆ）も含む。好ましくは、１４次モデルが適切だが、任意の次数を使用することができる。当業者であれば、たとえば、既知のＡＲＭＡ（自己回帰移動平均）モデルなどの全極型モデルの代わりに代替のフィルタモデルを使用してもよいことを認識し得る。さらに、本発明の好適な実施形態に関して上述した時間領域のフィルタリング作業とは全く異なり、周波数領域でフィルタリング作業を実行してもよい。
【００５１】
次に、ＬＰＣモデル係数は好ましくは（式４に従い）逆フィルタ内で使用され、音声信号はステップ１１５６で逆フィルタを通過する。音声信号が逆フィルタを通過することで入力信号が有効に均等化され、それによって図２のマスク伝達関数ＭＳＫ（ｆ）で生じたスペクトル歪（山と谷）が除去される。適切な固定化後フィルタを用いる、ステップ１１５８でのポストフィルタリングは理想的には、吸入ノイズの非白色度を補正し、以下の特定の符復号器または無線機の要件に最適に一致する規定の音調の質を音声信号に与えるために、均等化された信号上で実行される。このポストフィルタリングは、図２のノイズ経路伝達関数ＮＰ（ｆ）を補正するためにも使用することができる。
【００５２】
ＡＭＳＥ法８００のイコライザがエアレギュレータノイズに及ぼす影響を、２つの異なる次数の均等化フィルタに関して図１２に示す。具体的には、図１２は、均等化前の吸入ノイズバーストのスペクトル表示１２１０を示す。さらに、１４次の均等化フィルタ（１２２０）および２０次の均等化フィルタ（１２３０）を用いて、均等化後の吸入ノイズのスペクトルも示される。図示されるように、スペクトルの山は、２０次の均等化フィルタにより極めて十分に、且つ１４次均等化フィルタによりもある程度は十分に平坦化される。さらに、これらのフィルタにより均等化されるマスク音声に関する調音テストが示すように、音質は均等化されていない音声と比べ、均等化フィルタの使用により大幅に向上される。加えて、知覚された音質は２つのフィルタの次数間ではほとんど差が見られない。
【００５３】
ＡＭＳＥアルゴリズムアプローチの利点は、１）マスク音響共鳴特性を判定するための励起源として空気−マスクシステムに固有の定期的な、スペクトル上安定した、広帯域レギュレータノイズを使用すること、２）システム伝達関数モデリングが、簡単で定着した有効な手法を用いてリアルタイムで達成されること、３）均等化が同じ有効な手法を用いてリアルタイムで達成されること、および４）システム伝達関数モデルが変動し続ける状況に対してリアルタイムで継続的に順応可能なことである。当業者であれば、上記の利点が本発明のＡＭＳＥの実施形態に関連するすべての利点を包含することを意味するのではなく、その代表的なものとしての機能を果たすことを意味しているにすぎないことを認識し得る。
【００５４】
本発明の第３の態様は、図１に示されるシステム１００のような通信システムに連結される加圧空気供給システムにおいて、継続時間と吸入ノイズの周波数とを判定し、呼吸数と空気使用量を判定する方法と装置である。本発明の本実施形態による方法は、ここではＩＮＲＲＡ（吸入ノイズ呼吸数分析器）法とも称する。ＩＮＲＲＡ法は、人からの呼吸音を測定する代わりに、エアレギュレータにより生成される音を監視することにより呼吸を測定する実質上非間接的な方法である。ＩＮＲＲＡ法の根拠は、ＳＣＢＡなどの加圧空気呼吸システムが一方向の気流を有することである。空気は空気源とレギュレータを通じてのみシステム内に入ることができ、排気バルブを通じてのみシステム外に出ることができる。吸込および排気バルブは同時に開放不能である。よって、レギュレータ吸込バルブの動作は、使用者の呼吸サイクルに直接関連する。
【００５５】
レギュレータ吸込バルブの開放の１つの指標はレギュレータ吸入ノイズである。吸入ノイズは、ＳＣＢＡまたはその他の加圧空気供給システムマスクには入る高圧空気の産物である。人が息を吸い込むとき、レギュレータバルブを開放させ、圧縮タンク空気を導入させる、マスク内のわずかな負圧を生成するように、マスクは気密にされる。バルブを横断する乱気流は、ＳＣＢＡマスクに直接連結され、マイクロフォンに捕捉され、吸入毎に生じる大きな広帯域ヒス雑音を引き起こす。上述したように、ノイズは突然生じ、吸入の継続時間中に渡って一定の振幅を有し、非常に優れた開始および終了時間分解能を提供する。所与のマスクの種類と着用者に関しては、吸入ノイズのスペクトル特性は、開口サイズ、声道状況、および肺の気流などの要因に応じて相当変動する直接的な人間の呼吸音とは対照的に極めて安定している。ＩＮＲＲＡは、呼吸数の尺度としてエアレギュレータ吸入ノイズの安定性を利用する。
【００５６】
ＩＮＲＲＡは、全体のスペクトル特性により吸入ノイズの存在を特定する適合フィルタリングスキームを使用する。加えて、ＩＮＲＲＡは、万一生じた場合の、ノイズのスペクトル特性の変動に対して順応することにより、吸入ノイズとその他の音声との間の最適な識別を提供することができる。各吸入の開始を算出することによって、瞬間的な呼吸数とその時間平均を吸入ノイズの発生から容易に算出することができる。さらに、各吸入ノイズの終点を測定し、継続時間を算出し、予測可能なマスクレギュレータ流速についての情報を提供することにより、システムは空気流量を推定することができる。これは、吸入ノイズを録音するマイクロフォンからの信号のみを用いて達成可能である。
【００５７】
ＩＮＲＲＡ法１３００のブロック図を図１３に示す。ＩＮＲＲＡ法は、ノイズモデルマッチング部１３１０、ノイズ検出部１３３０、吸入息定義部１３５０、パラメータ推定部１３７０、およびノイズモデル更新部１３９０の５つの部分に分割可能である。ＩＮＲＲＡ法のノイズモデルマッチング部、ノイズ検出部、およびノイズモデル更新部は理想的には、先に詳述したＡＲＩＮＡ法の対応部分と同一である。したがって、簡潔にするため、これらの３つの部分の詳細な説明はここでは繰り返さない。しかしながら、ＩＮＲＲＡ法１３００の吸入息定義部１３５０およびパラメータ推定部１３７０を以下に詳細に説明する。
【００５８】
まず、吸入息定義部１３５０を説明する。ＩＮＲＲＡ法１３００の吸入息定義部１３５０の目的は、少なくとも１つの要因、たとえば、この場合、吸入息に対応する１つかまたはそれ以上の吸入ノイズバーストの終点と継続時間のセットに基づき吸入ノイズを特性化することである。吸入ノイズ検出部１３３０の決定は、１および０の値を用いる時係数ｍの関数として吸入ノイズの有無を表すステップ１３５２で、好ましくは二値信号、ＩＮＭｍ、ｍ＝０、１、２、．．．、Ｍ−１を生成するために使用される。この二値信号は、Ｍ個のサンプルの長さの回転バッファに記憶され、Ｍは少なくとも２つの吸入ノイズバースト、または最も遅い予測呼吸速度での呼吸の期間を包含するのに十分な二値信号のサンプルを記憶できるほど大きい。好適な実施形態では、これは結局およそ１５秒になる。この二値信号の時間分解能とＭの値とは、上述したように吸入ノイズモデルマッチング部に依存する、吸入ノイズ検出部１３３０で使用される最も短いサブフレーム時間によって決定され、どのスペクトルマッチング方法をステップ１３１０で使用するかに応じて、２０サンプル（２．５ｍｓｅｃ）または８０サンプル（１０ｍｓｅｃ）のいずれかである。
【００５９】
１３３０からの吸入ノイズ検出器の出力は必ずしも完全であるとは限らないため、吸入ノイズの検出中に誤検出が発生して、ノイズの実際の開始と継続時間に関して曖昧さを招く可能性がある。よって、ステップ１３５２により生成される二値吸入ノイズ信号は、ステップ１３５４で周知の移動平均型またはその他の適切なフィルタを用いて統合される。このフィルタによって、短い継続時間の誤検出が除去され、呼吸に対応する完全な吸入ノイズバーストを定義する、より正確な信号が生成される。ステップ１３５４で生成されるこの信号から、各ノイズバーストの正確な開始時間Ｓｉ、終了時間Ｅｉ、および呼吸継続時間Ｄｉのうち少なくとも１つの要因を、ステップ１３５６において処理フレーム継続時間精度内で判定することができる。二値信号ＩＮＭｍによって表されるような吸入ノイズバーストの開始時間および終了時間は、信号バッファ内の相対指数に着目することにより得られる。継続時間Ｄｉは次式のように単独の吸入ノイズバーストに対して定義される。
【００６０】
【数１４】
2008504587.tif 000016
ただし、ｉは、長さＭおよび期間Ｔ秒の二値信号バッファに存在するＩＴ吸入ノイズバーストのｉ番目を指定する。これらの吸入ノイズバースト係数値は理想的には、回転する有限長バッファ、ノイズバースト／呼吸ごとの１セットのパラメータに記憶される。ＩＮＲＲＡアルゴリズム部１３１０、１３３０、１３５２、および１３５４により処理されるＳＣＢＡマスクマイクロフォン音声の結果は図１４〜１７に示され、ＳＣＢＡを着用する男性の話し手の、静かな部屋で録音された音声に基づく。図１４は、ノイズバースト１４１０の混在する入力音声１４２０を示す。図１５は、吸入ノイズモデルマッチング部１３１０の出力のスペクトル歪尺度Ｄの時間−振幅表示１５００を示す。図１６は、吸入ノイズ検出器１３３０の二値出力時間−振幅表示１６００を示す。図１７は、未処理の検出器出力を組み込み、各吸入の継続時間を正確に定義する呼吸定義アルゴリズム１３５０の移動平均フィルタ成分１３５４の出力の時間−振幅表示１７００を示す。
【００６１】
パラメータ推定部１３７０は、吸入息定義部１３５０による吸入ノイズの特性係数に基づき推定可能なパラメータの例を示す。判定することのできる２つの上記パラメータの例は、使用者の呼吸数と概算吸入空気流量である。呼吸数は、吸入息定義部で判定可能な連続吸入ノイズバーストの順次的な開始時間Ｓｉを用いて容易に判定することができる。たとえば、分毎の「瞬時的」呼吸数は次式に従い算出できる。
【００６２】
【数１５】
2008504587.tif 000017
ただし、Ｓｉは秒単位の２つの連続するノイズバースト（吸入息）開始時間である。従って、平均呼吸数は次式より算出できる。
【００６３】
【数１６】
2008504587.tif 000018
ただし、ＩＴは、規定時間Ｔにおける、検出された連続する呼吸（吸入ノイズバースト）数である。
【００６４】
吸入中の概算空気流量は、吸入息定義部によって判定可能な呼吸の継続時間、および、たとえば、オフラインで判定可能な初期空気タンク充填圧力とレギュレータの平均流速とに関する追加情報から推定することができる。吸込バルブが開放されるとき、空気供給タンク圧がエアレギュレータの最低入力圧レベルを超えている限り、エアレギュレータはマスクにほぼ一定の圧力で（周囲空気／水圧の関数）空気量を送る。さらに、マスクに入る空気流速は、マスクレギュレータ吸込バルブが開放されている間、ほぼ一定である。よって、タンクから出されて、呼吸する人に運ばれる空気量は、吸込バルブが開放されている時間に比例する。バルブが開放されている時間は、各吸入ノイズの継続時間により測定可能である。
【００６５】
充填される際の供給タンク内の初期空気量は、タンク体積Ｖ０、充填圧力Ｐ０、ガス温度Ｔ０、一般ガス定数Ｒ、ガスのモル質量Ｎｍの関数であり、周知の理想的なガス式ＰＶ＝ＮｍＲＴにより算出することができる。初期充填圧力とタンクシリンダー体積は知ることができるので、タンクガスとマスクガスの温度が同一だと仮定すると、マスク圧力での呼吸に利用可能な空気量は次式により与えられる。
【００６６】
【数１７】
2008504587.tif 000019
吸入事象ｉ中に使用者に供給される空気の概算量は次式のようになる。
【００６７】
【数１８】
2008504587.tif 000020
ただし、ＩＶｉは空気量、Ｄｉは吸入ノイズから判定される吸入事象の継続時間、ＫＲは特定のエアレギュレータの気流速度に関連する校正係数である。ＫＲは個々のシステムについて実験的に導き出すこともできるし、あるいはメーカーのデータから判定することもできる。個々の吸入量ＩＶｉより、時間Ｔまでに使用される空気の概算総量ＶＴは次式に従い定義することができる。
【００６８】
【数１９】
2008504587.tif 000021
ただし、ＩＴは時間Ｔまでの吸入総数である。残りのタンク供給空気はしたがって次式のようになる。
【００６９】
【数２０】
2008504587.tif 000022
ＩＮＲＲＡ法の利点は、最小限の音声帯域で呼吸ノイズを捕捉するマイクロフォン信号を使用できることと、特別なセンサが不要なことである。別の利点は、呼吸検出器が、特徴が可変である人間の呼吸ノイズではなく、安定したスペクトル特性を有するエアレギュレータにより生成されるノイズを検出することに基づく。さらに別の利点は、他の種類の音響呼吸分析器のように特定周波数を調べることに固定されないことである。さらに、システムは、環境や異なる使用者および加圧空気呼吸マスクシステムの変動に自動的に順応する。よって、ＩＮＲＲＡ法は、システム１００の外部に、たとえば、監視するための無線データチャネルを介して自動的に送ることのできる貴重な情報である瞬時的または平均呼吸数および概算空気使用量のデータを継続的に提供することができる。当業者であれば、上記の利点が本発明のＩＮＲＲＡの実施形態に関連するすべての利点を包含することを意味するのではなく、その代表的なものとしての機能を果たすことを意味しているにすぎないことを認識し得る。
【００７０】
本発明による３つの方法（ＡＲＩＮＡ、ＡＭＳＥ、およびＩＮＲＲＡ）はすべて、好ましくは記憶装置に記憶されるソフトウェアアルゴリズム（上述のシステム１００によるシステムに含まれる）として実行され、そのステップはシステム１００のＤＳＰ１３８のような適切な処置装置において実行される。本発明の自己相関とＬＰＣフィルタリング法に相当するアルゴリズムは処理時間の半数を占め得る。しかしながら、ＡＲＩＮＡ、ＡＭＳＥ、およびＩＮＲＲＡ法に相当する、これらのアルゴリズムまたはアルゴリズムの全体は、小さなハードウェア占有場所で有効に実行することができる。さらに、本発明の別の実施形態では、ＡＭＳＥ法はＡＲＩＮＡ法のように多くの方法を利用するため、有効に組み合わせることができる。
【００７１】
本発明を特定の実施形態と組み合わせて説明してきたが、当業者は追加の利点と変形を容易に想到し得る。したがって、本発明はより広い側面では、図示され記載されるような特定の詳細、代表的な装置、および図示される例には限定されない。様々な変更、変形、および置換は、上述の説明に鑑み当業者にとっては自明である。たとえば、吸入ノイズを特定し減衰する方法を上述したが、本発明に関して説明した方法は、呼気ノイズ、または上述の方法を用いる有効な検出に役立つ擬似固定的スペクトル特性を有するその他の種類のノイズに適用させることもできる。よって、本発明は上述の説明に限定されず、添付の特許請求の範囲の思想と範囲に従うすべての変更、変形、および変種を包含するものと理解されるべきである。
【図面の簡単な説明】
【００７２】
【図１】通信システムに連結される呼吸用の加圧空気供給システムを含む従来技術システムの概略ブロック図である。
【図２】図１に示されるシステムのマスク−無線機音声経路を示す。
【図３】マスク内でおよびマスクマイクロフォン出力で測定されたスペクトル振幅反応と、マスク、マイクロフォン、およびマイクロフォン増幅器に関して算出された総合伝達関数の例を示す。
【図４】ＳＣＢＡエアレギュレータによって生じる吸入ノイズの例を示す。
【図５】図４に示される吸入ノイズの長期のスペクトルの大きさを示す。
【図６】所与のＳＣＢＡマスクを着用する単独の話し手から生じる吸入ノイズの４つの重複スペクトルを示す。
【図７】音声と混在する吸入ノイズバーストを示す、ＳＣＢＡマイクロフォンからの音声出力を示す。
【図８】本発明の一実施形態による吸入ノイズを検出し除去する方法の概略ブロック図である。
【図９】図８の方法で使用されるスペクトル照合回路の一実施形態の概略ブロック図である。
【図１０】図８の方法で使用されるスペクトル照合回路の別の実施形態の概略ブロック図である。
【図１１】本発明の別の実施形態による音声信号を均等化する方法の概略ブロック図である。
【図１２】本発明による１４次および２０次のＬＰＣ逆フィルタ均等化後のスペクトルと比べた場合の均等化前の吸入ノイズスペクトルを示す。
【図１３】生体計測パラメータの測定に使用される本発明の別の実施形態による、吸入ノイズの周波数の継続時間を判定し、呼吸数および空気使用量を判定するための方法の概略ブロック図である。
【図１４】音声およびエアレギュレーションの吸入ノイズを含む、マイクロフォン入力からの信号を示す。
【図１５】図１３に示される方法により判定される図１４に示される信号の平均正規化モデル誤差を示す。
【図１６】図１３に示される方法により生成される吸入ノイズ検出器の出力信号を示す。
【図１７】図１３に示される方法により生成される総合吸入検出器出力を示す。
【図１】
2008504587.tif 000023
【図２】
2008504587.tif 000024
【図３】
2008504587.tif 000025
【図４】
2008504587.tif 000026
【図５】
2008504587.tif 000027
【図６】
2008504587.tif 000028
【図７】
2008504587.tif 000029
【図８】
2008504587.tif 000030
【図９】
2008504587.tif 000031
【図１０】
2008504587.tif 000032
【図１１】
2008504587.tif 000033
【図１２】
2008504587.tif 000034
【図１３】
2008504587.tif 000035
【図１４】
2008504587.tif 000036
【図１５】
2008504587.tif 000037
【図１６】
2008504587.tif 000038
【図１７】
2008504587.tif 000039
【国際調査報告】
2008504587.tif 000040
2008504587.tif 000041