(57)【要約】【課題】ハウリングの抑制とエントレインメント耐性の向上を両立可能なフィードバックキャンセラ及び補聴器を提供する。【解決手段】フィードバックキャンセラは、レシーバ１１と、マイクロホン１２と、ＦＩＲフィルタ１３及び係数更新部１４からなる適応フィルタ２０と、減算部１７と、補聴処理部１０と、白色化フィルタ１５、１６を備えており、係数更新部１４は、白色化フィルタ１５、１６のそれぞれの出力信号に基づき、適応フィルタ２０の係数を更新するハイパーＨ∞フィルタを含んで構成される。【選択図】図１
2019087815.tif 000002
【特許請求の範囲】
【請求項１】
電気信号を音に変換する第１の変換手段と、
音を電気信号に変換する第２の変換手段と、
ＦＩＲフィルタ及び係数更新部により構成され、前記第１の変換手段から前記第２の変換手段へのフィードバック伝達関数を適応的に推定する適応フィルタと、
前記第２の変換手段の出力信号から前記適応フィルタの出力信号を減算して第１の信号を生成する減算部と、
前記第１の信号に所定の信号処理を施して、前記第１の変換手段に入力される第２の信号を生成する信号処理部と、
前記第１の信号を白色化する第１の白色化フィルタと、
前記第１の白色化フィルタと共通の白色化フィルタ係数により、前記第２の信号を白色化する第２の白色化フィルタと、
を備え、
前記係数更新部は、前記第１及び第２の白色化フィルタのそれぞれの出力信号に基づき、下記の（１）、（２）、（３）、（４）式で表されるフィルタリングアルゴリズムに従って、前記適応フィルタの係数を更新するハイパーＨ∞フィルタを含んで構成されており、
【数１２】
2019087815.tif 000015
ここで、
2019087815.tif 000016
前記ハイパーＨ∞フィルタの評価基準は、下記の（６）、（７）、（８）式で表される状態空間モデルに対し、フィルタ誤差に対応する項Ｔ１と、初期状態の誤差に対応する項Ｔ２と、システム雑音に対応する項Ｔ３と、観測雑音に対応する項Ｔ４と、評価関数の重みρとを用いた次の（１０）式、
（Ｔ１／ρ）／（Ｔ２＋Ｔ３＋Ｔ４／ρ）                  （１０）
の最大値が、予め与えられた上限値（γｆ２）より小さく抑えるように設定されることを特徴とするフィードバックキャンセラ。
【数１３】
2019087815.tif 000017
ただし、
ｘｋ：状態ベクトル（または状態）
ｘ０：初期状態
ｗｋ：システム雑音
ｖｋ：観測雑音
ｙｋ：観測信号
ｚｋ：出力信号
Ｈｋ：観測行列
ｘ＾ｋ｜ｋ：観測信号ｙ０〜ｙｋまでを用いた時刻ｋの状態ｘｋの状態推定値
ｘ＾０｜０：状態の初期推定値
Ｋｓ，ｋ＋１：フィルタゲイン
Σｗ ｋ：システム雑音の共分散行列に対応
Ｐ＾ｋ｜ｋ−１：ｘ＾ｋ｜ｋ−１の誤差の共分散行列に対応
Ｐ＾１｜０：初期状態の誤差の共分散行列に対応
「フィルタ誤差」：（１）式のｚｖｋ｜ｋと（８）式のｚｋ＝Ｈｋｘｋの差（（５）式のｅｆ，ｉに対応）の重み付きノルムの２乗の和
「初期状態の誤差」：時刻ｋ＝０のときの（６）式の状態ｘkの初期値ｘ０とその推定値の差の重み付きノルム
「システム雑音に対応する項」：（６）式のシステム雑音ｗｋの重み付きノルムの２乗の和
「観測雑音に対応する項」：（７）式の観測雑音ｖｋの重み付きノルムの２乗の和
【請求項２】
前記係数更新部は、前記ハイパーＨ∞フィルタを用いた第１の係数更新部と、ＮＬＭＳ(Normalized Least Mean Square)を用いた第２の係数更新部とを並列に接続して構成されることを特徴とする請求項１に記載のフィードバックキャンセラ。
【請求項３】
前記第１及び第２の白色化フィルタの各々は、ＮＬＭＳ又はＬＭＳによる適応アルゴリズムに従って動作する線形予測フィルタであることを特徴とする請求項１又は２に記載のフィードバックキャンセラ。
【請求項４】
前記ハイパーＨ∞フィルタは、下記の（１１）乃至（１９）式で表される演算を実行し、前記適応フィルタの係数を更新する高速Ｈ∞フィルタ（ＦＨＦ）であることを特徴とする請求項１から３のいずれか１項に記載のフィードバックキャンセラ。
【数１４】
2019087815.tif 000018
ここで、ｍｋ：Ｎ×２行列、μｋ：１×２行列であり、
ｃｋ：２×Ｎ行列Ｃｋ＝［ｃｋ，・・・，ｃｋ−Ｎ＋１］の第１列ベクトルである。
ただし、ｃｋ−ｉ＝０、ｋ−ｉ＜０と仮定する。
また、０＜ρ＝１−χ（γｆ）≦１、γｆ＞１とし、χ（γｆ）は、
χ（１）＝１、χ（∞）＝０を満たす任意の単調減少関数である。
また初期値は、
Ｋ−１＝０、ａ−１＝０、Ｓ−１＝１／ε０、ε０＞＞０、ｂ−１＝０、ｘ＾−１｜−１＝０である。
【請求項５】
前記ハイパーＨ∞フィルタは、下記の（２２）乃至（２８）式で表される演算を実行し、前記適応フィルタの係数を更新するＪ−高速Ｈ∞フィルタ（Ｊ−ＦＨＦ）あることを特徴とする請求項１から３のいずれか１項に記載のフィードバックキャンセラ。
【数１５】
2019087815.tif 000019
ただし、
【数１６】
2019087815.tif 000020
であり、ｄｉａｇ{・}は対角行列、Ｒｅ，ｋ（１，１）は行列Ｒｅ，ｋの（１，１）成分を表す。
【請求項６】
前記ハイパーＨ∞フィルタは、入力信号としてＭ次（Ｍ：正の整数）の自己回帰信号を用いることにより計算量を削減したＰ−高速Ｈ∞フィルタ（Ｐ−ＦＨＦ）であることを特徴とする請求項４に記載のフィードバックキャンセラ。
【請求項７】
前記Ｐ−ＦＨＦは、前記自己回帰信号を使用しない（Ｍ＝０）設定とすることにより白色化した白色化高速Ｈ∞フィルタ（ＷＦＨＦ）であることを特徴とする請求項６に記載のフィードバックキャンセラ。
【請求項８】
前記第１及び第２の白色化フィルタの収束速度に関わる動作パラメータは、前記ハイパーＨ∞フィルタの収束速度に応じた所定の範囲内で適宜に予め調節可能であることを特徴とする請求項１から７のいずれか１項に記載のフィードバックキャンセラ。
【請求項９】
前記動作パラメータは、前記第１及び第２の白色化フィルタにおける前記適応アルゴリズムで用いるステップサイズパラメータであることを特徴とする請求項８に記載のフィードバックキャンセラ。
【請求項１０】
請求項１から９のいずれか１項に記載のフィードバックキャンセラを備え、
前記第１の変換手段はレシーバであり、
前記第２の変換手段はマイクロホンであり、
前記信号処理部は、前記第１の信号に所定の補聴処理を施す補聴処理部である、
ことを特徴とする補聴器。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ハウリングの発生を抑制可能な構成を備えるフィードバックキャンセラ及び補聴器に関する。
【背景技術】
【０００２】
一般的な補聴器は、外部空間から伝わる音を収集するマイクロホンと、使用者の外耳道に音を出力するレシーバとを具備する。補聴器の使用時には、レシーバから出力された音が外耳道内から外部空間に漏れてマイクロホンにフィードバックすることにより、ハウリングが発生する場合がある。このようなハウリングの発生を抑制するための手段として、フィードバック伝達関数を適応的に推定する適応フィルタを用いたフィードバックキャンセラが広く知られている。この種のフィードバックキャンセラは、通常のハウリングの発生を抑圧するには有効であるが、周期的な信号が入力されたときに適応動作の不具合を招く場合がある。すなわち、適応フィルタへの入力信号が周期性を有する（自己相関が高い）場合、いわゆるエントレインメントを生じ、その入力信号が歪んで、異音を発生する現象が知られている。このようなエントレインメントに対する対策としては、例えば、適応フィルタの入力側に、白色化フィルタ（周波数等化ユニット）を挿入し、入力信号を白色化する構成が提案されている（例えば、特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特表２００７−５２５９１７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
前述したように、白色化フィルタを用いて適応フィルタへの入力信号を白色化することにより、フィードバックキャンセラにおけるエントレインメント（異音）の発生を抑制することができる。しかし、フィードバックキャンセラにおいてハウリングが発生した場合を想定すると、通常は適応フィルタの収束速度に比べて白色化フィルタの収束速度が速いために、まだ適応フィルタが収束しない時点でハウリング成分が白色化される状態になり、ハウリングの抑制が不十分となってハウリング抑制時間が長くなるという問題がある。一方、この問題を回避するため、白色化フィルタの収束速度が十分に遅くなるように制御する対策は可能であるが、この場合、ハウリングを抑制できたとしても、前述のエントレインメント耐性が大きく劣化するため、異音の発生が避けらない。
【０００５】
本発明はこれらの問題を解決するためになされたものであり、白色化フィルタと後述のハイパーＨ∞フィルタを用いたシンプルな構成により、複雑な制御を行うことなく、ハウリングの抑制とエントレインメント耐性の向上を両立し得るフィードバックキャンセラ及び補聴器を提供するものである。
【課題を解決するための手段】
【０００６】
上記課題を解決するために、本発明のフィードバックキャンセラは、電気信号を音に変換する第１の変換手段（１１）と、音を電気信号に変換する第２の変換手段（１２）と、ＦＩＲフィルタ（１３）及び係数更新部（１４）により構成され、前記第１の変換手段から前記第２の変換手段へのフィードバック伝達関数を適応的に推定する適応フィルタ（２０）と、前記第２の変換手段の出力信号から前記適応フィルタの出力信号を減算して第１の信号を生成する減算部（１７）と、前記第１の信号に所定の信号処理を施して、前記第１の変換手段に入力される第２の信号を生成する信号処理部と、前記第１の信号を白色化する第１の白色化フィルタ（１５）と、前記第１の白色化フィルタと共通の白色化フィルタ係数により、前記第２の信号を白色化する第２の白色化フィルタ（１６）とを備え、前記係数更新部は、前記第１及び第２の白色化フィルタのそれぞれの出力信号に基づき、後述の（１）、（２）、（３）、（４）式で表されるフィルタリングアルゴリズムに従って、前記適応フィルタの係数を更新するハイパーＨ∞フィルタを含んで構成されており、前記ハイパーＨ∞フィルタの評価基準は、後述の（６）、（７）、（８）式で表される状態空間モデルに対し、フィルタ誤差に対応する項Ｔ１と、初期状態の誤差に対応する項Ｔ２と、システム雑音に対応する項Ｔ３と、観測雑音に対応する項Ｔ４と、評価関数の重みρとを用いた後述の（１０）式の最大値が、予め与えられた上限値（γｆ２）より小さく抑えるように設定されることを特徴としている。
【０００７】
本発明のフィードバックキャンセラによれば、第１及び第２の白色化フィルタの動作によりエントレインメント（異音）の発生を抑制するとともに、ハウリングが発生した場合にはハイパーＨ∞フィルタのフィルタリングアルゴリズムにより短時間で適応フィルタが収束するので、ハウリング成分が完全に白色化される前に適応フィルタにより迅速かつ確実にハウリングを抑制することができる。よって、白色化フィルタの複雑な動作制御を行うことなく、ハウリングの抑制性能とエントレインメント耐性の両方を向上させることが可能となる。
【０００８】
本発明の係数更新部は、ハイパーＨ∞フィルタを用いた第１の係数更新部（１４）と、ＮＬＭＳ(Normalized Least Mean Square)を用いた第２の係数更新部（４１）とを並列に接続して構成することができる。このような構成により、従来のＮＬＭＳアルゴリズムと、本発明に特徴的なハイパーＨ∞フィルタとの両方の特性を反映させつつ適応フィルタ（４０）を動作させることができる。
【０００９】
本発明の第１及び第２の白色化フィルタとして、ＮＬＭＳ又はＬＭＳによる適応アルゴリズムに従って動作する線形予測フィルタを採用することができる。例えば、遅延部（３０）と、適応フィルタ（３１）と、減算部（３２）とにより構成される線形予測フィルタを適応白色化フィルタとして用いることができる。
【００１０】
本発明の係数更新部に用いるハイパーＨ∞フィルタとして、多様な下位概念のフィルタを採用することができる。例えば、後述の（１１）〜（１９）式で表される演算を実行し、適応フィルタの係数を更新する高速Ｈ∞フィルタ（ＦＨＦ）を用いることができる。また例えば、後述の（２２）〜（２８）式で表される演算を実行し、適応フィルタの係数を更新するＪ−高速Ｈ∞フィルタ（Ｊ−ＦＨＦ）を用いることができる。また例えば、ＦＨＦへの入力信号としてＭ次（Ｍ：正の整数）の自己回帰信号を用いることにより計算量を削減したＰ−高速Ｈ∞フィルタ(Ｐ−ＦＨＦ)を用いることができる。また例えば、Ｐ−ＦＨＦにおける自己回帰信号を使用しない（Ｍ＝０）設定とすることにより白色化した白色化高速Ｈ∞フィルタ（ＷＦＨＦ）を用いることができる。
【００１１】
本発明の第１及び第２の白色化フィルタの収束速度に関わる動作パラメータは、ハイパーＨ∞フィルタの収束速度に応じた所定の範囲内で適宜に予め調節しておくことが望ましい。これにより、第１及び第２の白色化フィルタは、そのエントレインメント抑制を維持する範囲内で、ハイパーＨ∞フィルタの収束速度に適合する収束速度で動作させることができる。この場合、第１及び第２の白色化フィルタにおける適応アルゴリズムで用いるステップサイズパラメータを典型的な動作パラメータの１つとして挙げることができる。
【００１２】
また、上記課題を解決するために、本発明の補聴器は、いずれかの前記フィードバックキャンセラを備え、前記第１の変換手段はレシーバであり、前記第２の変換手段はマイクロホンであり、前記信号処理部は、前記第１の信号に所定の補聴処理を施す補聴処理部であることを特徴としている。本発明の補聴器によれば、前述のフィードバックキャンセラの作用効果に基づき、エントレインメントによる異音の発生を防止し、かつ十分なハウリング抑制を確保することで、使用者にとって快適な補聴器を実現することができる。
【発明の効果】
【００１３】
以上述べたように、本発明によれば、エントレインメント対策として白色化フィルタを設けることに加え、白色化フィルタと適応フィルタの収束速度の違いに起因するハウリング成分の白色化を防ぐべく、係数更新部にハイパーＨ∞フィルタを導入したため、迅速かつ確実なハウリングの抑制とエントレインメント耐性の向上の両方を実現することが可能となる。
【図面の簡単な説明】
【００１４】
【図１】本発明を適用した第１実施形態の補聴器に関し、ディジタル信号処理に関連する具体的な構成例を示すブロック図である。
【図２】図１の白色化フィルタの具体例として、適応白色化フィルタの構成を示す図である。
【図３】本発明を適用した第２実施形態の補聴器に関し、ディジタル信号処理に関連する具体的な構成例を示すブロック図である。
【図４】第１実施形態の補聴器において、図１の構成例に関しての一変形例を示すブロック図である。
【図５】第２実施形態の補聴器において、図３の構成例に関しての一変形例を示すブロック図である。
【発明を実施するための形態】
【００１５】
以下、本発明を適用した２つの実施形態について添付図面を参照しながら説明する。以下の各実施形態では、フィードバックキャンセラを備えた補聴器に対して本発明を適用する例について説明する。
【００１６】
［第１実施形態］
図１は、本発明を適用した第１実施形態の補聴器に関し、ディジタル信号処理に関連する具体的な構成例を示すブロック図である。図１の構成例には、補聴処理部１０と、レシーバ１１と、マイクロホン１２と、ＦＩＲ(Finite Impulse Response)フィルタ１３及び係数更新部１４からなる適応フィルタ２０と、２つの白色化フィルタ１５、１６と、減算部１７とが示されている。このうち、レシーバ１１及びマイクロホン１２以外の構成要素は、例えば、ディジタル信号処理を実行可能なＤＳＰ(Digital Signal Processor)による信号処理によって実現することができる。図１の各構成要素は、補聴器の内部に搭載された電池（不図示）から供給される電源により動作する。なお、図示は省略しているが、レシーバ１１の入力側には、ディジタル信号をアナログ信号に変換するＤＡ変換器が設けられ、マイクロホン１２の出力側には、アナログ信号をディジタル信号に変換するＡＤ変換器が設けられる。
【００１７】
以上の構成において、補聴処理部１０は、減算部１７から出力される誤差信号ｅ（ｎ）を増幅するとともに、各々の使用者に適合して個別に設定された所定の補聴処理を施し、補聴処理後の信号ｓ（ｎ）を出力する。補聴処理部１０による補聴処理は、図１に示す伝達関数Ｇ（ｚ）で表す。補聴処理部１０によって適用可能な補聴処理としては、誤差信号ｅ（ｎ）に対する所定のゲインの付与に加えて、例えば、誤差信号ｅ（ｎ）に対するマルチバンドコンプレッション、ノイズリダクション、トーンコントロール、出力制限処理など、補聴器の使用者の聴力特性や使用環境に合わせた多様な処理を挙げることができる。
【００１８】
なお、補聴処理部１０に入力される誤差信号ｅ（ｎ）は、本発明の第１の信号に相当し、補聴処理部１０から出力される信号ｓ（ｎ）は、本発明の第２の信号に相当する。
【００１９】
レシーバ１１（本発明の第１の変換手段）は、例えば、使用者の外耳道内に設置され、補聴処理部１０から出力される前述の信号ｓ（ｎ）を音に変換して外耳道内の空間に出力する。レシーバ１１としては、例えば、電磁型のレシーバを用いることができる。また、マイクロホン１２（本発明の第２の変換手段）は、補聴器の外部空間から伝わる音を収集して電気信号に変換し、それを所望信号ｄ（ｎ）として出力する。マイクロホン１２としては、ＭＥＭＳ（Micro Electro Mechanical Systems）やコンデンサ型のマイクロホンを用いることができる。
【００２０】
図１において、マイクロホン１２には外部の環境音のみが入力されるのが理想的であるが、実際にはレシーバ１１から出力される音が外耳道から外部空間を経由してマイクロホン１２に回り込み、フィードバック音となる。このとき、レシーバ１１の入力からマイクロホン１２の出力に至るフィードバック経路におけるフィードバック伝達関数Ｆ（ｚ）を想定することができる。なお、レシーバ１１及びマイクロホン１２もそれぞれ固有の伝達関数を有するが、いずれもフィードバック伝達関数Ｆ（ｚ）に含めて考えることができる。フィードバック伝達関数Ｆ（ｚ）は、例えば、補聴器の構造、使用者の挙動（例えば、使用者が補聴器に手を近づける場合）、周囲の環境（例えば、自動車内）などに応じて変化し、補聴器のハウリングの要因となる。第１実施形態では、補聴器におけるハウリングの発生を抑制するために、後述するように、フィードバック伝達関数Ｆ（ｚ）を推定することでフィードバック成分をキャンセルする構成を有するが、詳細は後述する。
【００２１】
適応フィルタ２０のうちのＦＩＲフィルタ１３は、補聴処理後の信号ｓ（ｎ）に対し、係数更新部１４から供給される係数Ｗ（ｎ）を用いて、前述のフィードバック伝達関数Ｆ（ｚ）に対応する伝達関数Ｗ（ｚ）を適応的に推定するフィルタ演算を行い、出力信号ｙ（ｎ）を生成する。ＦＩＲフィルタ１３は、所定のタップ数（例えば、３２タップ）を有して構成される。なお、減算部１７は、所望信号ｄ（ｎ）からＦＩＲフィルタ１３の出力信号ｙ（ｎ）を減算し、それを前述の誤差信号ｅ（ｎ）として出力する。一方、ＦＩＲフィルタ１３と一体的に適応フィルタ２０を構成する係数更新部１４は、ＦＩＲフィルタ１３のフィルタ演算に用いる前述の係数Ｗ（ｎ）を順次更新する。
【００２２】
本発明の適用に際しては、係数更新部１４として、一般的なＮＬＭＳ(Normalized Least Mean Square)又はＬＭＳ(Least Mean Square)に基づくアルゴリズムを用いた構成に比べ、より収束速度の高速化が可能なハイパーＨ∞フィルタを採用する。なお、ハイパーＨ∞フィルタの具体的なアルゴリズム及びバリエーションについては後述する。
【００２３】
また、１対の白色化フィルタ１５、１６の主な役割は、無相関化された信号に応じて係数更新部１４を動作させることにより、エントレインメントに起因する異音の発生を抑制することにある。一方の白色化フィルタ１５は、誤差信号ｅ（ｎ）を入力し、誤差信号ｅ（ｎ）を白色化（無相関化）した出力信号を生成し、他方の白色化フィルタ１６は、信号ｓ（ｎ）を入力し、信号ｓ（ｎ）を白色化（無相関化）した信号を生成し、１対の白色化フィルタ１５、１６のそれぞれの出力信号が係数更新部１４に供給される。一方の白色化フィルタ１５に設定された白色化フィルタ係数は、他方の白色化フィルタ１６に対しても同様に設定されるので、２つの白色化フィルタ１５、１６は同一の特性を有している。白色化フィルタ１５、１６の具体的な構成例については後述する。
【００２４】
以下、図１の係数更新部１４としてのハイパーＨ∞フィルタのアルゴリズムについて具体的に説明する。まず、以下の説明で用いる記号を、次のように定義する。
ｘｋ：状態ベクトル（または状態）
ｘ０：初期状態
ｗｋ：システム雑音
ｖｋ：観測雑音
ｙｋ：観測信号
ｚｋ：出力信号
Ｈｋ：観測行列
ｘ＾ｋ｜ｋ：観測信号ｙ０〜ｙｋまでを用いた時刻ｋの状態ｘｋの状態推定値
ｘ＾０｜０：状態の初期推定値
Ｋｓ，ｋ＋１：フィルタゲイン
Σｗ ｋ：システム雑音の共分散行列に対応
Ｐ＾ｋ｜ｋ−１：ｘ＾ｋ｜ｋ−１の誤差の共分散行列に対応
Ｐ＾１｜０：初期状態の誤差の共分散行列に対応
【００２５】
上記記号のうち、ｘｋ、ｘ０、ｗｋ、ｖｋ、ｚｋはいずれも未知であり、ｙｋ、Ｈｋはいずれも既知である。なお、上記各記号に対し、上部に付される“＾”は、推定値であることを意味し、入力の都合上、記号の右上に記載しているが、実際には以下の数式に示すように記号の真上に記載される。また、Ｌ、Ｈ、Ｐ、Ｋ等は行列を意味し、原則として大文字で記載している。
【００２６】
第１実施形態のハイパーＨ∞フィルタの基本的なフィルタリングアルゴリズムは、下記の（１）、（２）、（３）、（４）式で表すことができる。なお、（２）式はフィルタ方程式であり、（３）式はフィルタゲインの式であり、（４）式はリカッチ方程式である。
【数１】
2019087815.tif 000003
ここで、（１）〜（４）式において、下記の（５）式を前提とする。
2019087815.tif 000004
【００２７】
また、次の（６）、（７）、（８）式で表される状態空間モデルを考える。
【数２】
2019087815.tif 000005
【００２８】
上記の状態空間モデルに対して、次の（９）式で表すＨ∞評価基準を想定する。
【数３】
2019087815.tif 000006
【００２９】
（９）式のＨ∞評価基準を簡潔化すると、フィルタ誤差に対応する項Ｔ１と、初期状態の誤差に対応する項Ｔ２と、システム雑音に対応する項Ｔ３と、観測雑音に対応する項Ｔ４と、評価関数の重みρとを用いた次の（１０）式の最大値が、予め与えられた上限値（γｆ２）より小さく抑えるように設定することに対応する。
（Ｔ１／ρ）／（Ｔ２＋Ｔ３＋Ｔ４／ρ）                  （１０）
ただし、「フィルタ誤差」とは、（１）式のｚｖｋ｜ｋと（８）式のｚｋ＝Ｈｋｘｋの差（（５）式のｅｆ，ｉに対応）の重み付きノルムの２乗の和であり、「初期状態の誤差」とは、時刻ｋ＝０のときの（６）式の状態ｘkの初期値ｘ０とその推定値の差の重み付きノルムであり、「システム雑音に対応する項」とは、（６）式のシステム雑音ｗｋの重み付きノルムの２乗の和であり、「観測雑音に対応する項」とは、（７）式の観測雑音ｖｋの重み付きノルムの２乗の和である。
【００３０】
前述のハイパーＨ∞フィルタは、その下位概念として、基本的な原理を共通としつつ異なる特徴を有する多様なバリエーションが存在する。以下では、ハイパーＨ∞フィルタ（Hyper H∞ Filter）のバリエーションとして、高速Ｈ∞フィルタ（Fast H∞ Filter）、Ｊ−高速Ｈ∞フィルタ（J-Fast H∞ Filter）、Ｐ−高速Ｈ∞フィルタ（Predictor-Based Fast H∞ Filter）、白色化高速Ｈ∞フィルタ（Whiten Fast H∞ Filter）の４つについて説明する。なお、基本概念であるハイパーＨ∞フィルタについての詳細は、例えば、特許第４０６７２６９号に記載されている。
【００３１】
まず、高速Ｈ∞フィルタ（以下、「ＦＨＦ」という）は、前述のハイパーＨ∞フィルタと厳密に等価であって、ハイパーＨ∞フィルタのゲイン行列Ｋｋの更新により、単位時間当たり１１Ｎの乗算回数で実行可能とすることで計算量を削減したアルゴリズムである。具体的には、ＦＨＦのアルゴリズムにおいては、下記の（１１）〜（１９）式で表される演算を実行する。なお、ＦＨＦについての詳細は、例えば、特許第４０６７２６９号に記載されている。
【数４】
2019087815.tif 000007
【００３２】
ここで、ｍｋ：Ｎ×２行列、μｋ：１×２行列であり、ｃｋ：２×Ｎ行列Ｃｋ＝［ｃｋ，・・・，ｃｋ−Ｎ＋１］の第１列ベクトルである。ただし、ｃｋ−ｉ＝０、ｋ−ｉ＜０と仮定する。また、０＜ρ＝１−χ（γｆ）≦１、γｆ＞１とし、χ（γｆ）は、χ（１）＝１、χ（∞）＝０を満たす任意の単調減少関数である。例えば、χ（γｆ）＝γｆ−２である。初期値は、Ｋ−１＝０、ａ−１＝０、Ｓ−１＝１／ε０、ε０＞＞０、ｂ−１＝０、ｘ＾−１｜−１＝０である。
【００３３】
また、ＦＨＦの存在条件は、次の（２０）式に示すスカラーの不等式と等価である。
【数５】
2019087815.tif 000008
ここで、（２０）式中のρ及びΞ＾ｉは、次の（２１）式で定義される。
【数６】
2019087815.tif 000009
【００３４】
次に、Ｊ−高速Ｈ∞フィルタ（以下、「Ｊ−ＦＨＦ」という）は、前述のハイパーＨ∞フィルタと厳密に等価であって、前述のＦＨＦよりも更に計算量を削減したアルゴリズムである。具体的には、Ｊ−ＦＨＦのアルゴリズムにおいては、下記の（２２）〜（２８）式で表される演算を実行する。なお、Ｊ−ＦＨＦについての詳細は、例えば、特許第４４４４９１９号に記載されている。
【数７】
2019087815.tif 000010
ただし、
【数８】
2019087815.tif 000011
であり、ｄｉａｇ{・}は対角行列、Ｒｅ，ｋ（１，１）は行列Ｒｅ，ｋの（１，１）成分を表す。
【００３５】
次に、Ｐ−高速Ｈ∞フィルタ（以下、「Ｐ−ＦＨＦ」という）は、前述のＦＨＦへの入力信号をＭ次（Ｍ：正の整数）の自己回帰信号を用いることにより、ＦＨＦよりも更に計算量を削減したアルゴリズムである。なお、Ｐ−ＦＨＦについての詳細は、例えば、特許第５１９６６５３号に記載されている。
【００３６】
次に、白色化高速Ｈ∞フィルタ（以下、「ＷＦＨＦ」という）は、前述のＰ−ＦＨＦにおいて、自己回帰信号を用いない０次（Ｍ＝０）に設定した場合のアルゴリズムである。また、Ｍ＝０のときのＰ−ＦＨＦに対応する次の（３４）式のＷＦＨＦも提案されている。
【数９】
2019087815.tif 000012
ここで
【数１０】
2019087815.tif 000013
ただし、Ｓ−１＝ε０-1である。
【００３７】
以下、図１の白色化フィルタ１５、１６の具体例について説明する。例えば、白色化フィルタ１５、１６としては、図２に示す構成を有する適応白色化フィルタを採用することができる。図２に示す適応白色化フィルタは線形予測フィルタとして機能し、遅延部３０と、適応フィルタ３１と、減算部３２とにより構成され、入力信号ｓｉ（ｎ）を白色化した出力信号ｓｏ（ｎ）を生成する。遅延部３０は、入力信号ｓｉ（ｎ）を１サンプルだけ遅延させた信号ｓｉ（ｎ−１）を出力する。適応フィルタ３１は、遅延部３０からの信号ｓｉ（ｎ）に対し、例えば、ＮＬＭＳ(Normalized Least Mean Square)アルゴリズムを用いてＷａ（ｚ）を適応的に推定し、信号ｓａ（ｎ）を生成する。減算部３２は、入力信号ｓｉ（ｎ）から、適応フィルタ３１の信号ｓａ（ｎ）を減算し、それを前述の出力信号ｓｏ（ｎ）として出力する。
【００３８】
ここで、図２の適応フィルタ３１におけるフィルタ係数Ｗａ（ｎ）は、例えば、次の（３５）式に従って順次更新される。
【数１１】
2019087815.tif 000014
なお、（３５）式に含まれるステップサイズμａは、図１の適応フィルタ２０の収束速度とのバランスを勘案して適切な値に設定される。一例として、μ＝０．００５に設定することができる。
【００３９】
従来のＮＬＭＳアルゴリズムを用いて係数更新を行う場合に対して、前述のＦＨＦ、Ｊ−ＦＨＦ及びＷＦＨＦを用いて係数更新を行う場合では、パスチェンジ性能及び耐エントレインメント性能をより向上することができる。
【００４０】
耐エントレインメント性能は、主に白色化フィルタ１５、１６の動作に基づくものである。一方、ハウリング抑制性能（パスチェンジ性能）は、従来のＮＬＭＳと比べて収束速度が格段に早いＦＨＦ、Ｊ−ＦＨＦ及びＷＦＨＦを用いて適応フィルタ２０を構成したことで向上したものである。例えば、ＮＬＭＳを用いて適応フィルタ２０を構成した場合には、仮にエントレインメント対策として白色化フィルタ１５、１６を用いたとしても、適応フィルタ２０の収束速度が遅くなるため、ハウリング成分が白色化されてしまい、ハウリング抑制が不十分となってハウリング抑制時間が長くなる弊害がある。このように、本発明の構成を採用すれば、ハウリング状態の検出等の複雑な制御を行うことなく、ハウリングの抑制性能とエントレインメント耐性の両方を向上可能であることが確認された。
【００４１】
［第２実施形態］
次に図３は、本発明を適用した第２実施形態の補聴器に関し、ディジタル信号処理に関連する具体的な構成例を示すブロック図である。図３の構成例において、補聴処理部１０、レシーバ１１、マイクロホン１２、ＦＩＲフィルタ１３、係数更新部１４、２つの白色化フィルタ１５、１６、減算部１７に関しては、第１実施形態と共通であり、それぞれの説明を省略する。図３のうち、図１と異なるのは、ハイパーＨ∞フィルタを用いた係数更新部１４に対し、ＮＬＭＳを用いた係数更新部４１を並列させた点である。
【００４２】
具体的には、図３に示すように、ＦＩＲフィルタ１３及び２つの係数更新部１４、４１により、適応フィルタ４０が構成される。すなわち、第２実施形態では、ハイパーＨ∞フィルタを用いた係数更新部１４（本発明の第１の係数更新部）とＮＬＭＳアルゴリズムを用いた係数更新部４１（本発明の第２の係数更新部）を組み合わせて、ＦＩＲフィルタ１３のフィルタ演算に用いる係数Ｗ（ｎ）を更新する。また、１対の白色化フィルタ１５、１６の各々の出力信号は、並列する係数更新部１４、４１の両方に入力される。そして、ＦＩＲフィルタ１３に供給される係数Ｗ（ｎ）は、２つの係数更新部１４、４１から供給されるそれぞれの係数の和で表すことができる。
【００４３】
第２実施形態の構成によれば、適応フィルタ４０のフィルタ演算に用いる係数Ｗ（ｎ）が、ハイパーＨ∞フィルタとＮＬＭＳアルゴリズムの両方に依存して決定されるので、適応フィルタ４０の動作には、第１実施形態のハイパーＨ∞フィルタの特性と従来のＮＬＭＳアルゴリズムの両方の特性が反映されることになる。ハイパーＨ∞フィルタを用いた係数更新部１４としては、例えば、収束速度の高速化と計算量の削減の両面で最も優位性がある前述のＷＦＨＦを採用することができる。第２実施形態の構成の基本的な効果は第１実施形態と概ね共通であるが、フィードバックパス変化によるハウリング抑制性能と、エントレインメント発生の抑制性能とのそれぞれの達成目標に応じて、第１実施形態と第２実施形態とを適切に使い分けることができる。
【００４４】
［変形例］
本発明の第１実施形態及び第２実施形態の補聴器を構成する場合、図１及び図３の構成例に限られることなく、多様な変形が可能である。図４は、第１実施形態の補聴器において、図１の構成例に関しての一変形例を示すブロック図である。図４の変形例においては、図１とは接続形態が変更されており、補聴処理部１０と、レシーバ１１と、マイクロホン１２と、２つのＦＩＲフィルタ１３、５０と、係数更新部１４と、２つの白色化フィルタ１５ａ、１６と、２つの減算部１７、５１とが示されている。本変形例では、図１とは異なり、マイクロホン１２の出力側の構成が２系統に分岐している。
【００４５】
すなわち、所望信号ｄ（ｎ）が白色化フィルタ１５ａを介して所望信号ｄ’（ｎ）に変換され、減算部５１は、所望信号ｄ’（ｎ）からＦＩＲフィルタ５０の出力信号ｙ’（ｎ）を減算し、それを誤差信号ｅ’（ｎ）として出力する。係数更新部１４は、誤差信号ｅ’（ｎ）と白色化フィルタ１６から出力される信号ｓ’（ｎ）とに基づいて、ＦＩＲフィルタ５０におけるフィルタ演算に用いる係数Ｗ（ｎ）を更新する。図４において、一方の白色化フィルタ１５ａに設定された白色化フィルタ係数は、他方の白色化フィルタ１６に対しても同様に設定されるとともに、一方のＦＩＲフィルタ５０に対して設定された係数Ｗ（ｎ）は、他方のＦＩＲフィルタ１３に対しても同様に設定される。なお、図４において、ＦＩＲフィルタ１３及び係数更新部１４が図１の適応フィルタ２０に対応すると同時に、ＦＩＲフィルタ５０及び係数更新部１４も図１の適応フィルタ２０に対応する。
【００４６】
次に図５は、第２実施形態の補聴器において、図３の構成例に関しての一変形例を示すブロック図である。図５の変形例においても、図３とは接続形態が変更されているが、図４の変形例と比べた場合、ハイパーＨ∞フィルタを用いた係数更新部１４に対し、ＮＬＭＳを用いた係数更新部４１を並列させた点のみが異なっている。従って、２つの係数更新部１４、４１を並列させる構成は、図３で説明した通りである。他の構成要素に関しても、図４と共通であり、説明を省略する。なお、図５において、ＦＩＲフィルタ１３及び２つの係数更新部１４、４１が図３の適応フィルタ４０に対応すると同時に、ＦＩＲフィルタ５０及び２つの係数更新部１４、４１も図３の適応フィルタ４０に対応する。
【００４７】
以上のように、図１の構成に代えて、図４の変形例を採用する場合であっても、第１実施形態と同様の作用、効果を得ることができる。また、図３の構成に代えて、図５の変形例を採用する場合であっても、第２実施形態と同様の作用、効果を得ることができる。
【００４８】
以上、本実施形態に基づき本発明の内容を具体的に説明したが、本発明は前述の実施形態に限定されるものではなく、その要旨を逸脱しない範囲で多様な変更を施すことができる。例えば、本実施形態では、フィードバックキャンセラを備えた補聴器に対して本発明を適用する場合を説明したが、本発明のフィードバックキャンセラを備えていれば、多様な機器やシステムに対しても広く本発明を適用することができる。例えば、本発明のフィードバックキャンセラをエコーキャンセラ（会議システム）に対しても適用することができる。また、レシーバは、音を発するレシーバに限らず、頭蓋骨などの振動を通じて知覚するための骨伝導レシーバであってもよい。その他の点についても上記実施形態により本発明の内容が限定されるものではなく、本発明の作用効果を得られる限り、上記実施形態に開示した内容には限定されることなく適宜に変更可能である。
【符号の説明】
【００４９】
１０…補聴処理部
１１…レシーバ
１２…マイクロホン
１３、５０…ＦＩＲフィルタ
１４、４１…係数更新部
１５、１５ａ、１６…白色化フィルタ
１７、５１…減算部
２０、４０…適応フィルタ
３０…遅延部
３１…適応フィルタ
３２…減算部
【図１】
2019087815.tif 000021
【図２】
2019087815.tif 000022
【図３】
2019087815.tif 000023
【図４】
2019087815.tif 000024
【図５】
2019087815.tif 000025