効率的なサブバンド適応ＦＩＲフィルタリング
【課題】音声チャネルを備えた音声システムの提供。
【解決手段】デジタルフィルタが各チャネルに実装されるデジタル音声プロセッサを備える音声システムであって、各チャネルのデジタルフィルタが広帯域入力音声信号を受信し、受信音声信号を複数サブバンドに分割し、等帯域幅のサブバンド信号を提供する様に構成され、かつ、サブバンド信号スペクトルが入力音声信号のスペクトルを構成する解析フィルタバンク２２、各サブバンド用で、各フィルタリングサブバンド信号を提供するサブバンドＦＩＲフィルタ、フィルタリングサブバンド信号を受信し、広帯域出力音声信号を提供する為に結合する様に構成される合成フィルタバンク２２‘を備え、遅延が各サブバンド信号に関連付けられ、サブバンド信号の内の１つの遅延が解析フィルタバンクの上流の広帯域入力音声信号に適用され、残留遅延が解析フィルタバンクの下流の残サブバンド信号に適用される、音声システム。
【解決手段】デジタルフィルタが各チャネルに実装されるデジタル音声プロセッサを備える音声システムであって、各チャネルのデジタルフィルタが広帯域入力音声信号を受信し、受信音声信号を複数サブバンドに分割し、等帯域幅のサブバンド信号を提供する様に構成され、かつ、サブバンド信号スペクトルが入力音声信号のスペクトルを構成する解析フィルタバンク２２、各サブバンド用で、各フィルタリングサブバンド信号を提供するサブバンドＦＩＲフィルタ、フィルタリングサブバンド信号を受信し、広帯域出力音声信号を提供する為に結合する様に構成される合成フィルタバンク２２‘を備え、遅延が各サブバンド信号に関連付けられ、サブバンド信号の内の１つの遅延が解析フィルタバンクの上流の広帯域入力音声信号に適用され、残留遅延が解析フィルタバンクの下流の残サブバンド信号に適用される、音声システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、音声アプリケーションでのサブバンド信号処理に関し、特に、所望されるターゲット転送関数の計算上効率的な実装を提供するサブバンドＦＩＲフィルタのセットを設計する方法に関する。
【背景技術】
【０００２】
ＩＩＲフィルタ（無限インパルス応答フィルタ）と対照的に、ＦＩＲフィルタ（融点インパルス応答フィルタ）は、所望される転送関数（つまり、大きさ及び位相応答）を有するデジタルフィルタを実現する可能性を提供する。したがって、所望される転送関数は、最小位相となるように設計できるだけではなく、最大位相、又は混合相となるように設計することもできる。さらに、しばしば音声信号処理で所望される、線形位相転送関数も実装できる。
【発明の概要】
【発明が解決しようとする課題】
【０００３】
しかしながら、特に音声信号処理では、動作中のメモリ要件だけではなく計算上の作業も増加させる、ＦＩＲフィルタを使用するときには、必要とされるフィルタ長は、しばしばかなり長くなる。これは、おもに２つの物理的な要因による。第１に、大部分の音声用途に存在する室内インパルス応答の減衰期間はかなり長く、その結果、ＦＩＲフィルタを使用するときには、相応して長いフィルタ長につながる。第２に、人間の聴覚系は、異なる周波数バンドで不均一な周波数分解能を提供するように適応されている。人間の聴覚系は、低周波はきわめてよく分解する。つまり、周波数差は低い絶対周波数でよく認識されるのに対し、高周波はそれほど容易に識別されない。例えば、両方の場合とも周波数差は１００Ｈｚであるが、１００Ｈｚのトーンが２００Ｈｚのトーンから容易に区別され得るのに対し、人間の耳では５１００Ｈｚトーンから５０００Ｈｚのトーンを区別するのが困難である。すなわち、人間の聴覚系の周波数分解能は周波数が高くなるにつれて減少する。この現象は周知であり、バーク尺度、メル尺度、及びＥＲＢ（同等矩形帯域幅）尺度等の人間の聴覚系に適応された心理音響周波数尺度を形成する。
【０００４】
研究は、内装（カーペット、布張りをした家具等）の結果として、エネルギーの分解が低速であるため、試聴室の室内インパルス応答が、特に低周波で著しく長いことを示している。音声再生システムによって生じる音圧がバス周波数範囲（つまり、２００Ｈｚ未満）で最大であるのに対し、人間の聴覚系は低周波音声信号により敏感でなくなるという事実によって、この影響はさらに強くなっている。
【０００５】
これら全ての要因を統合した見解は、典型的な試聴室の特性だけではなく、人間の聴覚系の特性によって、音声信号処理システムで十分な音声品質を提供するためには、ＦＩＲフィルタの長さが特定の最小長さよりも短くないことが要求される。例えば、バス周波数範囲で約１０Ｈｚの必須周波数分解能を提供するには、４４１００Ｈｚのサンプリング周波数での音声チャネルごとに、４４１０フィルタ係数のＦＩＲフィルタが必要となり、それによって現代的な音声システムは最大１６のチャネルを有する。
【０００６】
長いＦＩＲフィルタが動作中に多大な計算作業及び／又は高メモリ要件を伴うので、標準的なデジタル音声信号処理の使用を可能にするためには、音声アプリケーションでのＦＩＲフィルタの効率的な実装に対するニーズがある。
【課題を解決するための手段】
【０００７】
少なくとも１つの音声チャネルのある音声システムが開示される。音声システムは、チャネルごとに少なくとも１つのデジタルフィルタが実装されるデジタル音声プロセッサを含む。各チャネルのデジタルフィルタは、広帯域入力音声信号を受信し、その入力音声信号を複数のサブバンドに分割し、このようにして等しい帯域幅を有するサブバンド信号を提供し、サブバンド信号のスペクトルが入力音声信号のスペクトルを作るように構成される解析フィルタバンクと、各サブバンド用であり、このようにしてそれぞれフィルタリングされたサブバンド信号を提供するサブバンドＦＩＲフィルタと、フィルタリングされたサブバンド信号を受信し、１つの広帯域出力音声信号を提供するためにそれらを結合するように構成される合成フィルタバンクとを含み、遅延は各サブバンド信号に関連付けられ、サブバンド信号の内の１つの遅延が解析フィルタバンクの上流にある広帯域入力音声信号に適用され、残留遅延が解析フィルタバンクの下流の残りのサブバンド信号に適用される。対応するフィルタリング方法が開示される。
【０００８】
本発明は、以下の図面及び説明を参照するとよりよく理解できる。図中、類似する参照番号は、対応する部分を示す。
【０００９】
例えば、本願発明は以下の項目を提供する。
（項目１）
少なくとも１つのデジタルフィルタが各チャネルに実装されるデジタル音声プロセッサを備える少なくとも１つの音声出力チャネルのある音声システムであって、各チャネルの前記デジタルフィルタが、
広帯域入力音声信号を受信し、前記受信した音声信号を複数のサブバンドに分割し、このようにして等しい帯域幅を有するサブバンド信号を提供するように構成される解析フィルタバンクであって、前記サブバンド信号の前記スペクトルが前記入力音声信号の前記スペクトルを構成する解析フィルタバンクと、
各サブバンド用のサブバンドＦＩＲフィルタであって、このようにしてそれぞれのフィルタリングされたサブバンド信号を提供するサブバンドＦＩＲフィルタと、
前記フィルタリングされたサブバンド信号を受信し、広帯域出力音声信号を提供するためにそれらを結合するように構成される合成フィルタバンクと、
を備え、
遅延が各サブバンド信号に関連付けられ、前記サブバンド信号の内の１つの前記遅延が前記解析フィルタバンクの上流の前記広帯域入力音声信号に適用され、前記残留遅延が前記解析フィルタバンクの下流の前記残りのサブバンド信号に適用される、
音声システム。
（項目２）
前記遅延されていない入力音声信号を受信し、追加サブバンド信号として、前記解析フィルタバンクでタッピングされていない最低周波数サブバンドの前記サブバンド信号を提供する低域通過フィルタと、
前記それぞれのサブバンドＦＩＲフィルタの上流又は加硫で前記最低周波数サブバンド信号に適用される遅延と、
をさらに備え、
前記広帯域入力音声信号に適用される前記遅延が、前記解析フィルタバンクでタッピングされる前記サブバンドと関連付けられる最短の遅延に一致する、
上記の項目に記載のシステム。
（項目３）
前記サブバンドＦＩＲフィルタのフィルタ係数の前記数が、銭サブバンドフィルタの前記周波数分解能が、少なくともほぼ前記人間の耳の前記周波数分解能に相当するように設定される、
上記項目のいずれか１つの項目に記載のシステム。
（項目４）
少なくとも第１の音声チャネルと、第２の音声チャネルとを備え、
各広帯域入力音声信号が、前記解析フィルタバンクでタッピングされる全てのチャネルの全てのサブバンド信号の最短の遅延によって与えられる等しい遅延で遅延される、
上記項目のいずれか１つの項目に記載のシステム。
（項目５）
前記解析フィルタバンク（複数の場合がある）及び前記合成フィルタバンク（複数の場合がある）が、前記入力音声信号をブロック単位で処理する偶数で積み重ねられるＧＤＴＦフィルタであり、前記サブバンド信号が、前記ブロックサイズに対応する係数で、前記広帯域信号に比較してダウンサンプリングされ、
前記処理されたサブバンド信号の前記最低周波数及び前記最高周波数のサブバンド信号が実数値であり、前記他のサブバンド信号が複素数値である、
上記項目のいずれか１つの項目に記載のシステム。
（項目６）
前記低域通過フィルタが、前記解析フィルタバンク（複数の場合がある）によって使用されるプロトタイプフィルタに相当する、上記項目のいずれか１つの項目に記載のシステム。
（項目７）
音声システムの少なくとも１つの音声出力チャネル内に提供される少なくとも１つの音声信号をフィルタリングするための方法であって、各チャネルのために、
所与の周波数バンド内にあるスペクトルを有する、フィルタリングされる広帯域入力音声信号を受信することと、
前記入力音声信号の前記周波数バンドを、複数のサブバンドに分割し、このようにしてサブバンド信号を提供することであって、前記サブバンド信号の前記スペクトルが前記入力音声信号の前記スペクトルを構成する、提供することと、
各サブバンド信号をフィルタリングし、このようにしてサブバンドごとにそれぞれのＦＩＲ又はＩＩＲフィルタを使用してそれぞれのフィルタリングされたサブバンド信号を提供することと、
広帯域出力音声信号を提供するために、前記フィルタリングされたサブバンド信号を結合することと、
を含み、
遅延が各サブバンドに関連付けられ、前記サブバンドの内の１つの前記遅延が、前記周波数バンドを分割する前に前記広帯域入力音声信号に適用され、前記残留遅延が前記周波数バンドを分割後に前記残りのサブバンド信号に適用される、
方法。
（項目８）
前記周波数バンドを分割することが、前記入力信号をブロック単位で処理する第１の偶数で積み重ねられたＧＤＴＦ解析フィルタバンクを使用して実行され、
前記周波数バンドを分割することが、前記広帯域信号と比較して、前記ブロックサイズに相当する係数分、前記サブバンド信号を前記ダウンサンプリングすることをさらに含み、
前記フィルタリングされたサブバンド信号を結合することが、前記第１のフィルタバンクに対応し、前記広帯域出力信号をブロック単位で提供する第２の偶数で積み重ねられたＧＤＴＦ合成フィルタバンクを使用して実行される、
上記の項目に記載の方法。
（項目９）
前記遅延していない入力音声信号を低域通過フィルタリングし、追加サブバンド信号として、前記最低周波数サブバンドの前記サブバンド信号を提供することと、
それぞれのＦＩＲ又はＩＩＲフィルタを使用して前記最低周波数サブバンド信号をフィルタリングし、このようにして対応するフィルタリングされたサブバンド信号を提供することと、
ＦＩＲフィルタリングの前又は後に前記最低周波数サブバンド信号を表示することと、
をさらに含む、上記項目のいずれか１つの項目に記載の方法。
（項目１０）
前記サブバンドＦＩＲフィルタのフィルタ係数の前記数が、前記サブバンドフィルタの前記周波数分解能が前記人間の耳の前記周波数分解能に少なくともほぼ一致するように設定される、上記項目のいずれか１つの項目に記載の方法。
（項目１１）
前記音声システムが、第１の音声出力チャネル及び第２の音声出力チャネルを備え、
各広帯域入力音声信号が、前記周波数バンドを除算することによって得られる全てのチャネルの全てのサブバンド信号の最短遅延によって与えられる等しい遅延で遅延される、
上記項目のいずれか１つの項目に記載の方法。
（項目１２）
前記低域通過フィルタリングが、低域通過フィルタとして、前記周波数バンドを分割するために使用される前記ＧＤＦＴフィルタバンク（複数の場合がある）によって使用されるプロトタイプフィルタに相当するフィルタを使用することを含む、
上記項目のいずれか１つの項目に記載の方法。
【００１０】
（摘要）
少なくとも１つの音声チャネルのある音声システムが開示される。音声システムは、チャネルごとに少なくとも１つのデジタルフィルタが実装されるデジタル音声プロセッサを含む。各チャネルのデジタルフィルタは、広帯域入力音声信号を受信し、複数のサブバンドに入力音声信号を分割し、このようにして等しい帯域幅を有するサブバンドを提供するように構成される解析フィルタバンクであって、サブバンドの信号のスペクトルが入力音声信号のスペクトルを構成する解析フィルタバンクと、各サブバンド用であり、このようにしてそれぞれのフィルタリングされたサブバンド信号を提供するサブバンドＦＩＲフィルタと、フィルタリングされたサブバンド信号を受信し、それらを結合し、広帯域出力音声信号を提供するように構成される合成フィルタバンクとを含み、遅延が各サブバンド信号と関連付けられ、サブバンド信号の内の１つの遅延が解析フィルタバンクの上流で広帯域入力音声信号に適用され、残留遅延が解析フィルタバンクの下流で残りのサブバンド信号に適用される。
【図面の簡単な説明】
【００１１】
【図１】ＦＩＲフィルタのフィルタ係数の適応計算のための基本信号処理構造を示し、ＦＩＲフィルタは所定のターゲット関数Ｐ（ｚ）にほぼ一致する転送関数Ｇ（ｚ）を表す。
【図２】図１の構造の変型を示し、それによって図１の全帯域ＦＩＲフィルタＧ（ｚ）がサブバンドＦＩＲフィルタＧｍ（ｚ）のセットで置き換えられる。
【図３】全帯域フィルタＧ（ｚ）に代わることができるサブバンドＦＩＲフィルタＧｍ（ｚ）のセットを含む信号処理構造を示す図である。
【図４】それぞれがサブバンドＦＩＲフィルタ及び遅延回線を含む、多くのサブバンドを含む適応ＦＩＲフィルタ設計のための信号処理構造を示す図である。
【図５】サブバンドＦＩＲフィルタの長さを適応するときに低周波サブバンドを強調するためのバーク尺度に基づいた重み付け係数を示す図である。
【図６】図偶数で積み重ねられたフィルタバンク及び奇数で積み重ねられたフィルタバンクの規模応答を概略で示す図である。
【図７ａ】複数のサブバンドフィルタ内で分割されるＦＩＲフィルタの構造を概略で示す図である。
【図７ｂ】図７ａの例の簡略化された代替策を概略で示す図である。７ａ
【図８】メモリ要件が削減された図６のフィルタバンクの代替策を示す図である。
【図９】図８の例を改善するためのオプションを示す図である。
【図１０】２つの別々の音声チャネル内のフィルタを含む例を示す図である。
【発明を実施するための形態】
【００１２】
図１は、ＦＩＲフィルタ２０のフィルタ係数ｇｋ （ｋ ＝ ０， １， ．．．， Ｋ−１）の適応計算のための基本信号処理構造を示し、それによって下付き文字ｋはフィルタ係数の指数を示し、Ｋはフィルタ長を示す。ＦＩＲフィルタ２０は、フィルタ係数ｇｋの適応が成功した後に、基準システム１０の所定のターゲット関数Ｐ（ｚ）にほぼ一致する（離散）転送関数Ｇ（ｚ）を有する。適応フィルタ設計手順を実行するために、基準フィルタ１０及びＦＩＲフィルタ２０には、例えば、ホワイトノイズ又はＦＩＲターゲット転送関数Ｐ（ｚ）の通過帯域を含む帯域幅を有する任意の他の信号である、試験信号（入力信号ｘ［ｎ］）が供給される。ＦＩＲフィルタ２０の出力信号ｙ［ｎ］は、基準システム１０の出力信号から、つまり所望される信号ｄ［ｎ］から減算される（減算器３０）。差ｄ［ｎ］−ｙ［ｎ］は、エラー信号ｅ［ｎ］として使用され、適応部２１に供給される。適応部２１は、各サンプル時間間隔中に、エラー信号及び入力信号ｘ［ｎ］（本文脈では基準信号としても示される）からＦＩＲフィルタ係数ｇｋの更新されたセットを計算するように構成される。最小二乗平均（ＬＭＳ）アルゴリズム又は正規化最小二乗平均（ＮＬＭＳ）アルゴリズムが、フィルタ係数の適応のために利用され得る。ただし、この目的のためには、さまざまな適応アルゴリズムが活用されてもよい。適応アルゴリズムの収束後、ＦＩＲフィルタ係数ｇｋは、ターゲット転送関数Ｐ（ｚ）の最適近似である転送関数Ｇ（ｚ）を表す。
【００１３】
ＦＩＲフィルタを使用するときの計算作業を削減するための１つのオプションは、フィルタリングされる信号のスペクトルを多くの狭帯域信号（サブバンド信号）に分割し、各狭帯域信号を別々にフィルタリングすることである。全帯域信号のいくつかのサブバンド信号への分割は、いわゆる解析フィルタバンク（ＡＦＢ）によって実現され得る。同様に、サブバンド信号は、対応する合成フィルタバンク（ＳＦＢ）を用いて単一の全帯域信号に（再）結合され得る。以下では、全帯域信号は、例えば所望される信号ｄ［ｎ］のように下付き文字なしで示され、ｎは時間インデックスである。さらに、例えばｄｍ［ｎ］のような下付き文字を有する信号は、対応する全帯域信号ｄ［ｎ］の分解であるサブバンド信号のセットを示す。それによって、下付き文字ｍは、サブバンドの数（ｍ ＝ １， ２， ．．．， Ｍ）を示す。同様に、離散全帯域転送関数Ｇ（ｚ）は、多くのサブバンド転送関数Ｇｍ（ｚ）に分解され得る。
【００１４】
図２は、図１の基本信号処理構造に類似した、それによって適応ＦＩＲフィルタ２０が、サブバンドＦＩＲフィルタのセット２０’で置き換えられる基本信号処理構造を示す。この目的のため、全帯域入力信号ｘ［ｎ］は、ＡＦＢ２２を使用することによって数Ｍのサブバンド入力信号ｘｍ［ｎ］（ｍ ＝ １， ２， ．．．， Ｍである）に分割される。同様に、所望される全帯域信号ｄ［ｎ］は、ＡＦＢ１１を使用して（再び、ｍ ＝ １， ２， ．．．， Ｍ）数Ｍのサブバンド信号ｄｍ［ｎ］に分けられる。各サブバンドＦＩＲフィルタは、狭帯域転送関数Ｇｍ（ｚ）を実現し、下付き文字ｍは再びサブバンドの数を示す。各サブバンドフィルタＧｍ（ｚ）は、そのフィルタ係数ｇｍｋによってあらわされてもよく、それによってｋは、ｋ＝０からｋ＝Ｋｍ−１（Ｋｍは、ｍ番目のサブバンド内のフィルタＧｍ（ｚ）のフィルタ長さである）に及ぶフィルタ係数のインデックスを示す。各ＦＩＲフィルタＧｍ（ｚ）は、対応するエラー信号ｅｍ［ｎ］ ＝ ｄｍ［ｎ］ − ｙｍ［ｎ］を受信し、それぞれのサブバンドｍごとに更新されたフィルタ係数ｇｍｋ （ｋ ＝ １， ２， ．．．， Ｋｍ−１） のそれぞれのセットを計算する、適応部と関連付けられている（適応部のセットは図２で数字２１’で示される）。
【００１５】
図１の例と同様に、Ｍ個のサブバンドＦＩＲフィルタＧｍ（ｚ）の各１つのフィルタ係数ｇｍｋは、適応アルゴリズムの収束後に、全てのサブバンド転送関数Ｇｍ（ｚ）の組み合わせから生じる総転送特性が所定のターゲット関数Ｐ（ｚ）に一致するように適応される。
【００１６】
適切なフィルタ係数ｇｍｋの計算後、ＦＩＲフィルタＧｍ（ｚ）のセット２０’は、音声信号をフィルタリングするための図３に示されるように、解析フィルタバンク（ＡＦＢ２２）と対応する合成フィルタバンク（ＳＦＢ２２’）との間で操作されてよい。この場合、ＡＦＢ２２、ＦＩＲフィルタバンク２０’、及びＳＦＢ２２’は、例えば音声システムの等化フィルタを表すことがある、ターゲット関数Ｐ（ｚ）にほぼ一致する転送関数Ｇ（Ｚ）をともに実装する。現代の音声システムでのように、リスナーにとって所望されるサウンドインプレッションを生じさせるために、大きさだけではなく位相が等化にさらされる。したがって、ターゲット関数Ｐ（ｚ）は、一般的に、非線形位相特性のある非最小位相フィルタを表す。
【００１７】
解析フィルタバンクと、サブバンドフィルタのセットと、合成フィルタバンクとを含む、図３の信号処理構造が、（計算作業及びメモリ要件という点で）同じ転送関数の「通常の」ＦＩＲフィルタよりもより効率的であるのかという質問に対する回答は、とりわけ解析フィルタバンク及び合成フィルタバンクの効率的な実装形態の可用性に依存している。人間の聴覚系の不均一な周波数分解能を説明するために、低い中心周波数があるサブバンドの帯域幅が、より高い中心周波数のあるサブバンドの帯域幅よりも狭いフィルタバンクを使用することが望ましいだろう。かかる心理音響的に動機付けされる全帯域信号の、帯域幅が可聴周波数範囲内のそれぞれのサブバンドの位置に依存するサブバンド信号のセットへの分割を実現するいくつかの手法が存在する。しかしながら、異なる帯域幅のサブバンドのセットに１つの入力スペクトルを不均一に分割することを可能にする効率的なフィルタバンクは既知ではない。それにも関わらず、技術では、全帯域信号の、等しい帯域幅のサブバンド信号のセットへの分割を可能にする他の方法が既知である。一例は、単一のプロトタイプフィルタ、及びほぼあらゆる信号処理環境で利用できるＦＦＴアルゴリズムを利用する、Ｓ．Ｗｅｉｓｓらによって説明されるオーバーサンプリングされたＧＤＦＴフィルタバンクの高速実装形態である（ＩＥＥ Ｅｌｅｃｔｒｏｎｉｃｓ Ｌｅｔｔｅｒｓの第３６巻、１５０２−１５０３ページ、２０００年のＳ． Ｗｅｉｓｓ、「Ｆａｓｔ Ｉｍｐｌｅｍｅｎｔａｔｉｏｎ ｏｆ Ｏｖｅｒｓａｍｐｌｅｄ Ｍｏｄｕｌａｔｅｄ Ｆｉｌｔｅｒ Ｂａｎｋｓ」を参照すること）。
【００１８】
上記を要約すると、効率的な実装形態は不均一な帯域幅のサブバンドを処理するために利用できないので、等しい帯域幅のサブバンドを用いて動作するフィルタバンクを使用する必要がある。ただし、等しく幅広いサブバンドの機能不全を緩和するために、それぞれのサブバンドに割り当てられるＦＩＲフィルタの異なるフィルタ長が選ばれてよい。すなわち、ＦＩＲフィルタは、高周波分解能が要求されるサブバンドにおいてよりも、低周波分解能が要求されるサブバンド内でより少ないフィルタ係数を含む。高周波分解能が要求されるサブバンドは、通常、、可聴周波数範囲の下部にあるサブバンドである。したがって、人間の聴覚系の周波数分解能に相当する周波数分解能は、等しく幅広いサブバンドを用いて動作する効率的なフィルタバンクを使用することによって達成され得る。
【００１９】
上述されたように、ターゲット関数Ｐ（ｚ）は、一般に、周波数上での非線形グループ遅延特性を有する非最小位相フィルタである。異なるサブバンドでの異なるグループ遅延から生じる異なる信号伝搬遅延を補償するために、遅延回線は、遅延回線の上流又は下流にある各サブバンドＦＩＲフィルタに接続されてよい。したがって、遅延等化は、計算上、非効率的である追加のＦＩＲフィルタ係数を使用して達成される必要はない。フィルタ係数の数だけではなく、遅延値も実現されるターゲット転送関数Ｐ（ｚ）（つまり、大きさ及び位相応答）に依存するので、フィルタ係数の数（つまり、フィルタ「タップ」の数、及び遅延値は、「適応タップ割り当て」アルゴリズム及び「適応遅延割り当て」アルゴリズムを使用して、本明細書中でいかに説明されるように、サブバンドごとに適応的に決定され得る。その結果、Ｍ個のサブバンドＦＩＲフィルタＧｍ（ｚ）を設計するときには、フィルタ係数（図１の係数ｇｍｋを参照すること）だけではなく、係数の数ＫＭ及び追加遅延△ｍも適応的に決定される。図４は、サブバンドＦＩＲフィルタ設計に使用される信号処理構造を示す。図４の例は、図２の構造の改良版であり、各サブバンド信号経路に追加の遅延、及びグローバル「適応割り当て及び遅延割り当て部」を備えている。４は、図２の構造の改良版を示し、４の例は、図２の構造の改良版であり、各サブバンド信号経路に追加の遅延、及びグローバル「適応割り当て及び遅延割り当て部」を備えている。
【００２０】
図２の例でのように、全体入力信号ｘ［ｎ］（例えば、帯域が制限されたホワイトノイズ）が、ターゲット転送関数うＰ（ｚ）のあるシステム１０に供給され、このようにして所望される信号ｄ［ｎ］を生成する。入力信号ｘ［ｎ］だけではなく、所望される信号ｄ［ｎ］も、それぞれ数Ｍのサブバンド信号ｄｍ［ｎ］ 及び ｘｍ［ｎ］に分割される。図４の例は、最初のサブバンド及び最後のサブバンド（ｍ ＝ １及びｍ ＝ Ｍ）と関連付けられた構成要素及び信号だけを示す。サブバンド入力信号ｘｍ［ｎ］は、転送関数Ｇｍ（ｚ）がある適応ＦＩＲフィルタに供給され、このようにしてフィルタリングされたサブバンド信号ｙｍ［ｎ］を生成する。フィルタリングされた各サブバンド信号ｙｍ［ｎ］は、対応する所望される信号ｄｍ［ｎ］から差し引かれ、サブバンドごとにエラー信号ｅｍ［ｎ］を生じさせる。適応部は、それぞれのＦＩＲフィルタＧｍ（ｚ）のフィルタ係数ｇｍｋ（つまり、フィルタのインパルス応答｛ｇｍ０， ｇｍ１， ．．．， ｇｍ（Ｋ−１）｝）を最適化するための各ＦＩＲフィルタＧｍ（ｚ）に割り当てられ、それによってフィルタ係数ｇｍｋの最適セットがそれぞれのエラー信号ｅｍ［ｎ］のノルム（例えば、累乗）を最小限に抑える。
【００２１】
遅延Dｍを提供する遅延回線は、各サブバンドＦＩＲフィルタＧｍ（ｚ）に上流で又は下流で接続される。さらに、以下にさらに詳細に説明される適応タップ割り当て及び適応遅延割り当てアルゴリズムに従って、遅延回線Dｍの対応する遅延値だけではなく、ＦＩＲフィルタＧｍ（ｋ）のフィルタ長Ｋｍも動的に適応するように構成される、「適応タップ割り当て及び適応遅延割り当て部」４０が設けられる。
【００２２】
上述された適応タップ割り当て（つまり、ＦＩＲフィルタ長の適応）にはさまざまな手法が考えられ得る。１つの手法は、総エラー信号ｅ［ｎ］（それによって、ｅ［ｎ］ ＝ ｅ１［ｎ］ ＋ ｅ２［ｎ］ ＋ ．．．＋ ｅＭ［ｎ］）が最小値に達するまで、サブバンドＦＩＲフィルタＧｍ（ｚ）のフィルタ長Ｋｍを変えることである。実際には、かかる方法は優れた結果を出すが、、フィルタ係数の数の各変更後に、適応フィルタが再び収束する時間を要するので、きわめて時間がかかる。上述の手法と同様に優れた結果を出すが、要する時間がはるかに少ない別の手法は、Ｓ個の一番端のフィルタ係数ｇｍ（Ｋｍ−１）， ｇｍ（Ｋｍ−２）， ．．．， ｇｍ（Ｋｍ−Ｓ）のエネルギーを検討する。サブバンドのフィルタＧｍ（ｚ）のフィルタ長Ｋｍは、上述のＳ個の一番端のフィルタ係数のエネルギーがほぼ等しくなるまで変えられる。この手法は、サブバンドフィルタが経時的に急激に減衰することを必要とし、それは現実世界のシステムではつねに当てはまらなくてはならない。各サブバンドフィルタのＳ個の一番端のフィルタタップのエネルギーを比較することによって、サブバンドフィルタＧｍ（ｚ）がどれほどよくターゲット関数Ｐ（ｚ）を近似し、信号減衰挙動がターゲット関数Ｐ（ｚ）のインパルス応答の信号減衰挙動に似ているサブバンドフィルタインパルス応答を達成するように、サブバンドフィルタ全体でフィルタ係数を分散し直すためのガイドラインとなり、最小限に抑えることができるエラーに関して最適とみなすことができる。
【００２３】
以下に、適応タップ割り当てアルゴリズムの例がさらに詳しく説明される。他のＭ／２サブバンド信号は、第１のＭ／２サブバンドの信号の共役複素数のコピーであるため、実数値の全帯域入力信号ｘ［ｎ］の場合（図２を参照すること）、Ｍ／２サブバンドだけを処理する必要があることに留意するべきである。したがって、それぞれのサブバンドＦＩＲフィルタ転送関数は、以下の関係性に準拠し、
【００２４】
【数１】
それによって、アスタリスクは共役複素数演算子を示す。
【００２５】
したがってサブバンドフィルタＧｍ（ｚ）のフィルタ長Ｋｍ（それによって、ｍ ＝ １， ２， ．．．， Ｍ／２）は、Ｑ個のサンプル（つまり、サブバンドシステム内のサンプル）の期間で修正される。しかしながら、全てのサブバンドフィルタＧｍ（ｚ）のフィルタ係数の総数ｇｍｋ［ｎ］は一定のままである。すなわち、１つ又は複数のサブバンドフィルタのフィルタ長が増加する場合、別のサブバンドフィルタのフィルタ長は、フィルタ係数の総数を一定に保つように削減されなければならない。したがって、Ｑ個のサンプルの期間では、Ｍ／２サブバンドＦＩＲフィルタのそれぞれの長さは、ΔＫ 個の係数分削減される。したがって、さらに以下に説明される特定の基準に従って、Ｍ個のサブバンドフィルタ全体で分散し直されるΔＫ・Ｍ／２個の「自由な」係数がある。
【００２６】
上記の「分散し直し」は、以下の方程式で表されてよい。
【００２７】
【数２】
上式ではｍ ＝ １， ２， ．．．であり、Ｍ／２はサブバンドの数を示す。式 ｃｍ［ｎ／Ｑ］は、フィルタタップ（つまり、フィルタ係数）の分散の上述された基準を表す。上述されたように、１つの有用な基準は、サブバンドエラー信号ｅｍ［ｎ］のエネルギーである。この場合、ｃｍ［ｎ／Ｑ］は以下のように表すことができる。
【００２８】
【数３】
上式では、ｍ ＝ １， ２， ．．．であり、Ｍ／２及びＲはエラー信号が平均化されるサンプルの数である。ただし、適応サブバンドＦＩＲフィルタは、方程式（３）が評価される前に集束しなければならず、結果的に、ＲはＱよりもはるかに小さくなる必要がある（つまり、Ｒ ＜＜ Ｑ）。別の基準は、それぞれのサブバンドＦＩＲフィルタの一番端のＳ個のフィルタ係数のエネルギーを検討する。この場合、ｃｍ［ｎ／Ｑ］ は以下として表すことができ、
【００２９】
【数４】
上式ではｍ ＝ １， ２， ．．．であり、Ｍ／２及びＫｍはそれぞれのサブバンドでの現在のフィルタ長Ｋｍ［ｎ］である。代わりに、サブバンド入力信号ｘｍ［ｎ］のｎエネルギーは、一番端のＳ個のフィルタ係数とともに検討され得る（方程式５を参照すること）。この場合、ｃｍ［ｎ／Ｑ］は、以下として表すことができる。
【００３０】
【数５】
上述されたように、方程式（３）に従った基準は最善の結果を生じさせるが、評価するのに時間がかかる。方程式（５）に従った基準は、ＡＥＣ（音響エコーキャンセル）システムの場合でのように、ターゲットシステムが時間的に変化し、入力信号が任意に着色されるときに使用されるべきである。本ケースでは、つまり入力信号ｘ［ｎ］がホワイトノイズとなるように設計者によって選択されることがあるサブバンドＦＩＲフィルタ設計の場合、方程式（４）は優れた品質結果を出すと同時に、高速適応を可能にする。
【００３１】
心理音響的な態様を説明するために、方程式（３）、（４）、又は（５）に定義される式ｃｍは、対応する重み付け係数ｗｍで重み付けされてよい。つまり、式 ｃｍ［ｎ／Ｑ］は、方程式（３）、（４）、又は（５）の ｗｍ・ｃｍ［ｎ／Ｑ］で置き換えられる。重み付け係数ｗｍは、人間の聴覚系の周波数分解能が検討されるように選ばれる必要がある。バーク尺度を使用すると、係数ｗｍ は、以下のように計算でき、
【００３２】
【数６】
上式では、ｆｃ，ｍ は、ｆｃ，ｍ ＝ （２ｍ−１）・ｆＳ／（２・Ｍ）として計算されてよく、ｆＳ がＨｚ単位のサンプリング周波数であるｍ番目のサブバンドの（Ｈｚ単位の）中心周波数を示す。方程式（６）に従って計算された正規化重み付け係数ｗｍの例のセットは、図５に示され、それによって多くの１６のサブバンドが検討された（Ｍ／２＝１６）。
【００３３】
上述されたように、方程式（２）に従った係数の分散し直しは、サブバンドエラー信号ｅｍ［ｎ］（方程式（３）と比較すること）のエネルギー又はサブバンドフィルタ係数ｇｍｋ（ｋ ＝ ｛Ｋｍ−Ｓ， ．．．， Ｋｍ−１｝であり、方程式（４）と比較すること）の部分（例えば、一番端のＳ個のフィルタ係数）のエネルギー、又はサブバンドフィルタ係数ｇｍｋ及びサブバンド入力信号ｘｍ［ｎ］ （方程式（５）と比較すること）のエネルギーの組み合わせを「等化する」。したがって、異なるエネルギーを元にした基準が考えられる。アルゴリズムの「目標」は、サブバンド全体で信号エネルギーを等しく分散することである。この目標は最適化プロセスの所望される結果であるので、上述されたアルゴリズムも最小化（又は最大化）タスクとして見ることができる。
【００３４】
理想的な場合、上述された基準ｃｍ［ｎ／Ｑ］（方程式（３）から（５）を参照すること）は、各サブバンドＭで等しい値を有し、それによって以下となる。
【００３５】
【数７】
すなわち、サブバンドと見なされるＭ／２の内のそれぞれ１つで、基準ｃｍ［ｎ］は、全てのサブセット上で総計された ｃｍ［ｎ］という合計値よりも小さい係数２／Ｍである。上述された最小化タスクに戻ると、最小限に抑えられる品質関数は、以下のように書くことができる。
【００３６】
【数８】
上述されたように、ＦＩＲフィルタＧｍ（ｚ）は、一般的に、異なる信号伝搬遅延を生じさせる異なるサブバンドでの異なるグループ遅延を有する。かかる異なる遅延を補償するために、各ＦＩＲフィルタは各サブバンドでフィルタの上流で又は加硫で調整可能な遅延回線に接続される。再び、サブバンドの遅延値は、「適応遅延割り当て」として以下に説明されるように反復して最適化できる。
【００３７】
「適応タップ割り当て」として上述されたサブバンド内のＦＩＲフィルタに適応的に割り当てられるＦＩＲフィルタ係数は、有限インパルス応答を生じさせる矩形ウィンドウ関数での乗算により切り捨てられた無限インパルス応答と見なすことができる。各サブバンドＦＩＲフィルタのために、この矩形ウィンドウ関数が（時間軸に沿って）シフトされてよく、それによってそれぞれの時間シフトは対応するＦＩＲフィルタに接続される遅延回線の有効遅延を表す。
【００３８】
本明細書に説明される適応遅延割り当てアルゴリズムは、各サブバンドに、それぞれのＦＩＲフィルタの最大エネルギー有限インパルス応答を生じさせる、遅延値（つまり、上述された矩形ウィンドウの時間シフト）を見つけることをターゲットにしている。すなわち、それぞれの遅延回線に遅延値を割り当てるとき、エネルギー又は対応するＦＩＲフィルタ係数だけが検討される。しかし、必ずしもこれが唯一の成功する手法であるわけはない。代替策として、遅延値は、総エネルギー信号ｅ［ｎ］ ＝ ｅ１［ｎ］ ＋ ｅ２［ｎ］ ＋ ．．．＋ ｅＭ／２［ｎ］のノルムが最小限に抑えられるように選ばれてよい。すでに適応タップ割り当てアルゴリズムに関して上述されたように、かかる基準を使用すると、十分に優れた結果が生じるが、収束は遅くなる。代わりに、遅延値は、それぞれのＦＩＲフィルタ係数の総エネルギーが最大となるように選ばれてよい。
【００３９】
しかしながら、サブバンドフィルタのＦＩＲインパルス応答（つまり、ＦＩＲフィルタ係数）は、一般的に、サブバンドへの変換に起因する原因因子と因果因子を含むので、最適遅延値を見つけるための他の基準はそれぞれのインパルス応答の原因部分及び因果部分のエネルギーを等化するようにターゲットされ得る。つまり、それぞれのサブバンドＦＩＲインパルス応答の第１の半分のエネルギーと第２の半分のエネルギーの差異は最小限に抑えられなければならない。かかる最適化戦略は、以下のように簡潔に要約されてよい。
【００４０】
サブバンドｍごとに、それぞれのサブバン路ＦＩＲインパルス応答ｇｍｋの第１の部分及び第２の部分を定義し、ＦＩＲインパルス応答の第１の部分及び第２の部分のエネルギーを計算する。Ｖ個のサンプルの期間で以下を繰り返す。
（Ａ） ＦＩＲインパルス応答 ｇｍｋのエネルギーが、インパルス応答の第２の部分のエネルギーよりも大きい場合には、それぞれの遅延回線Dｍの遅延をＰタップ（例えば、Ｐ＝１）分増加する。
（Ｂ） ＦＩＲインパルス応答ｇｍｋの第１の部分のエネルギーが、インパルス応答の第２の部分のエネルギーよりも低い場合には、それぞれの遅延回線DｍをＰタップ分（例えば、Ｐ＝１）減少する。ただし、最小遅延はゼロである。
（Ｃ）両方のエネルギーが等しい場合は、遅延を未変更のままにする。
【００４１】
ＦＩＲインパルス応答を等しい長さの２つの部分に分割したくなることがある。しかし、これは最大値がインパルス応答のほぼ真中に位置するインパルス応答に対してだけ良好な結果を生じさせる。左側にその最大値を有する原因インパルス応答を検討すると、上記の最適化戦略は、エネルギーが等しく分散されるまで、上述の矩形ウィンドウを右に時間シフトし、このようにしてインパルス応答の左側にゼロを挿入し、これは明らかに望ましくない。したがって、ＦＩＲインパルス応答の上述された第１の部分及び第２の部分がそれぞれのインパルス応答の全長をカバーするのではなく、変数Ｕによって定められ、０＜Ｕ＜１である、その上部マージン及び下部マージンだけをカバーする場合に、実践では、よりすぐれた結果が達成される。例を示すと、Ｕ＝０．６である場合は、エネルギーを計算するときにフィルタ係数の６０パーセントだけが検討され、それによって下部３０パーセントは上述の第１の部分を形成し、上部３０パーセントは対応する第２の部分を形成する。したがって、ｍ番目のサブバンドでは、それぞれ第１の部分及び第２の部分のエネルギーＥｍＬ ａｎｄ ＥｍＵは以下の通りである。
【００４２】
【数９】
【００４３】
【数１０】
したがって、ＥｍＵ ＜ ＥｍＬの場合、遅延値Dｍは増加され、ＥｍＵ ＞ ＥｍＬの場合減少され、ＥｍＵ ＝ ＥｍＬの場合未変更のままとされる、又は以下の数学用語で表される。
【００４４】
【数１１】
一般に、適応遅延割り当てアルゴリズムは、最小化タスクとして見ることができ、それによって、各サブバンドで、それぞれの遅延値Dｍが、品質基準 ｃ’ｍ［ｎ／Ｖ］ ＝ ｜ＥｍＵ−ＥｍＬ｜が最小値に達するまで変えられ、それによってインデックスｎ／Ｖは、基準がＶ番目のサンプル時間ごとしか評価されないことを示す。
【００４５】
タップの数分１つのサブバンドの遅延回線の遅延値Dｍを調整する（例えば、増加する）ときには、対応するフィルタインパルス応答（ＦＩＲフィルタ係数）が左にシフトされ、このようにして（遅延回線及びＦＩＲフィルタの遅延を含む）総遅延を一定に保つことに留意する必要がある。これによって、遅延値Dｍの変動に起因するＦＩＲフィルタ係数の長い再適応化に対するニーズが回避される。ただし、適応タップ割り当てだけではなく適応遅延割り当てが実行される適応期間は、ＬＭＳ／ＮＬＭＳアルゴリズムの適応ステップ幅よりもかなり長くなければならず、したがってＦＩＲフィルタ係数ｇｍｋの適応は、適応タップ割り当て及び適応遅延割り当てが開始されるときには定常状態にある。
【００４６】
各サブバンドの追加遅延の適応調整だけではなく、フィルタ長の適応調整も有用に組み合わせることができるが、それぞれの適応方法（適応タップ割り当て及び適応遅延割り当て）は、互いとは無関係に単独で使用され得る。例えば、フィルタ係数の数は（例えば、バーク尺度に従って）サブバンドごとに事前に定義されてよく、遅延 Dｍは、上述されたように適応的に調整されてよい。さらに、フィルタ長の適応決定は、別々の遅延部を使用しなくても、又は一定の所定の遅延を用いて、上述されたように適応タップ割り当て方法に従って形成されてよい。これは、特に、所望されるターゲット関数Ｐ（ｚ）が線形位相又は最小位相転送関数である場合に適切であり得る。
【００４７】
効率的なサブバンドＦＩＲフィルタ構造を設計するためのフィルタ設計方法を説明してきたが、かかるサブバンドＦＩＲフィルタ構造を効率的に実装するための例が、さらに詳しく説明される。第１に、解析フィルタバンク及び合成フィルタバンクが詳しく説明される。サブバンドＦＩＲフィルタを実装するための効率的な信号処理構造は、次いで、図７及び図８に関して説明される。
【００４８】
上述されたように、解析フィルタバンク及び対応する合成フィルタバンク（図３のフィルタバンク２２及び２２‘を比較すること）は、実数値の全帯域入力信号ｘ［ｎ］を複素数値サブバンド信号のセットｘｍ［ｎ］に分割し、複素数値のフィルタリングされたサブバンド信号ｙｍ［ｎ］を実数値の出力信号ｙ［ｎ］に再結合する。しかし、実数値サブバンド信号だけと動作する他のフィルタバンクも使用できる。この場合、サブバンドＦＩＲフィルタはより効率的に実装され得る。しかし、実数値フィルタバンクは、対応する複素数値フィルタバンクの２倍の多さのＭＡＣ（乗算―蓄積動作）を必要とする。さらに、サブバンドは複雑なサブバンド信号を使用するとき係数Ｍによってアンダーサンプリングされることがあるのに対し、対称的に、実数値フィルタバンクを使用するときは、最大アンダーサンプリング係数はＭ／２に過ぎない。最後に、サブバンドＦＩＲフィルタの必須フィルタ長は実数値サブバンド信号を使用するときに必要となるフィルタ長の半分に過ぎない。これは、複素数値ＡＦＢ及びＳＦＢの使用には、実数値フィルタバンクよりも少ない計算作業が必要になることを示している。しかし、両方のタイプのフィルタバンクとも、本発明での使用に適用できる。
【００４９】
複素数値フィルタバンクを使用するとき、偶数で積み重ねられたフィルタバンク又は奇数で積み重ねられたフィルタバンクのどちらかを選ぶことができる。両方のタイプとも、汎用化離散フーリエ変換アルゴリズム（ＧＤＦＴアルゴリズム）によって効率的に実装され得る。図６は、奇数に積み重ねられた（図６ｂ）フィルタバンクだけではなく偶数で積み重ねられた（図６ａ）のサブバンドも示す。一見して、Ｍ個のサブバンドの半分は、Ｍこのサブバンドの他方の半分の複素共役コピーであり、結果的に数Ｍ／２の複素数値サブバンド信号だけが処理されなければならないので、奇数で積み重ねられた（図６ｂ）のフィルタバンクの使用が好ましく見える。封数で積み重ねられたＧＤＦＴフィルタバンクを使用するとき、最低のサブバンド及び最高のサブバンドは実数値であるが、他のサブバンドの半分の帯域幅しか有していない。その結果、数Ｍ／２＋１のサブバンド信号が処理されなければならず、それによって最低のサブバンドの信号は実数値である。結果として、複素数値サブバンドと比較すると、メモリ要件は半分にすぎず、必須ＭＡＣは最低のサブバンドの４分の１にすぎない。最低のサブバンドのＦＩＲフィルタ （Ｇ１（ｚ））が、（例えば、バーク又は任意の他の心理音響的に適応された周波数尺度に従って周波数分解能と比較されるとき）最も多くのフィルタ係数ｇ１ｋ を有するという事実を考えると、処理されるサブバンドの総数はＭ／２からＭ／２＋１に増加するが、総効率は必要とされるメモリに関して２５パーセント向上し、必要とされるＭＡＣに関して５０パーセント向上する。
【００５０】
以下では、複数のサブバンドフィルタＧｍ（ｚ）及びそれぞれの遅延Dｍに分解された（所望される）デジタルＦＩＲフィルタの実装形態が、説明例として音声システムのアプリケーションを使用して説明される。図７ａに示される構造は、音声チャネルごとに実装されてよい。すなわち、各音声チャネルに解析フィルタバンク２２及び対応する合成フィルタバンク２２‘が提供される。この文脈では、「音声チャネル」は増幅器によって提供される出力チャネルを指し、入力信号の数は指さない（増幅器が２つの入力チャネル（ステレオ入力）しか有さないが、例えば５．１サラウンドサウンド出力には６つの出力チャネルを有してよいことに留意すること）。デジタルＦＩＲフィルタは、上述されたフィルタ設計方法に従って設計され得る。しかし、任意のＦＩＲフィルタは、本明細書に説明されるフィルタ構造及び方法を使用して実装されてよい。解析フィルタバンク２２（ｍ＝１，２，３、．．．、Ｍ）によって提供されるＭ個のサブバンド信号ｘｍ［ｎ］ のそれぞれは、（対応するサブバンドフィルタＧｍ（ｚ）を適用することによって）別々にフィルタリングされ、（対応する遅延Dｍを適用することによって）遅延される。上述されたように、ＦＩＲフィルタＧｍ（ｚ）のフィルタ長が、人間の聴覚系の周波数特性を考慮に入れるように設定されるのに対し、等しく幅広いサブバンドが使用される場合、フィルタバンクは計算上効率的にしか実現できないので、個々のサブバンドの帯域幅は等しい。それぞれの遅延Dｍは、ＦＩＲフィルタＧｍ（ｚ）の必要とされるフィルタタップを削減するために各サブバンドに適用される。一般的に、フィルタ長及び遅延は、任意の設定することができ、上述される適応タップ割り当て方法及び適応遅延割り当て方法に従って調整されてよい。遅延を除き、図７ａのフィルタ構造は、図３に関してすでに説明された基本構造に相当する。
【００５１】
フィルタバンク２２及び２２‘は、さらに詳細に以下に説明される。図６に関してすでに説明されたように、偶数で積み重ねられたＧＤＦＴフィルタバンクの適用は、人間の耳の周波数分解能に一致するサブバンドフィルタ長を使用するときに特に望ましい。すなわち、最低サブバンドに適用されるＦＩＲフィルタは、最高の周波数分解能、したがって最高の数のフィルタ係数を有さなければならない。偶数で積み重ねられたＧＤＦＴフィルタバンクを使用するとき、弧の最低のサブバンドは実数値であり、したがってメモリ要件は４分の１削減され、必要される計算時間は、奇数で積み重ねられたＧＤＦＴフィルタバンクが使用される代替策に比較して半減される。２つ以上の音声チャネルには、通常、サブバンドで実装されるＦＩＲフィルタが提供され、したがってフィルタバンクのこの選択はかなりの改善を伴う。
【００５２】
メモリ要件は、サブバンド数と必要とされるメモリとの間の依存性を検討するときにさらに削減されてよい。一方、サブバンドの増加する数は、個々のサブバンドの帯域幅の減少を伴う。その結果、プロトタイプフィルタは、必ずさらに多くのフィルタ係数を生じさせるより険しい通過帯域／ストップバンド遷移を有さなければならず、したがって特に遅延回線を実装するために必要とされるメモリ量が増加する。他方、多数のサブバンドは、（ＦＦＴアルゴリズムで使用される）より大きなブロックサイズを伴い、その結果として、内部メモリを外部メモリにスワップするためにより多くの計算時間を使用することができ、このようにして必要とされる内部メモリを削減する。サブバンド数を決定するときには妥協策を見つけ出さなければならないことが上記から明らかとなるはずである。分析は、音声用途の場合、Ｍ＝１６、Ｍ＝３２、及びＭ＝６４のサブバンドのために良好な結果を得ることができ、それによってＮ＝Ｍ／２のアンダーサンプリングが適用される。すなわち、Ｍ＝３２サブバンドを使用するとき、サブバンド信号のサンプリングレートは、Ｎ＝１６という係数で削減されてよく、ＦＦＴブロックサイズも１６に設定される。後者の例（Ｍ＝３２）では、１７のサブバンド（Ｍ／２＋１）だけが実際に処理されなければならず、それによって（最高のサブバンドだけではなく）最低のサブバンドは上述されたように実数値である。
【００５３】
大量の計算時間を消費するだろう長い初期化を生じさせずに、可能な限り多くの内部メモリをより安価な外部メモリにスワップすることも望ましい。サブバンド遅延Dｍを実装するために必要とされるメモリは、外部メモリにスワップされるには適切であると思われる。しかし、これは多数のメモリ転送初期化を生じさせるだろう。例えば、１０の音声チャネルの内の１７のサブバンドがブロックごとに（読み取り及び書き込みのために）１７回のメモリ転送初期化を必要とする。これは初期化毎に３ミップスを仮定するとき、１０２０ミップスの計算力が必要とされることを意味する。この問題は、図７ｂに示される完全帯域（広帯域）ドメインの中に最も短いサブバンド遅延を「移動させる」ことによって、図７ａに示される構造を再配置するときに緩和できる。図７ｂの例では、最も短いサブバンド遅延DＭＩＮが第１の遅延D１である。広帯域ドメインの対応する遅延は、「バルク遅延」として示され、DＢＵＬＫ＝Ｎ・DＭＩＮに等しい。残りのサブバンド遅延は、図７ｂに示されるようにDＭＩＮ分削減される。
【００５４】
上述されたように、１７回の初期化の代わりに（読み取りアクセス用及び書き込みアクセス用に）１つの転送初期化だけしか必要でないので、広帯域バルク遅延は、外部メモリにスワップされてよい。残りの―ただし、今はさらに短い―サブバンド遅延が内部メモリに残る。ただし、図７ｂに示されている解決策はさらに改善され得る。最低サブバンドのＦＩＲフィルタが最長インパルス応答（つまり、フィルタタップの最大数）、したがって残りの遅延よりもはるかに短い最短遅延を有するので、他のサブバンド遅延のために必要とされる（内部）メモリは依然として多い。この状況は、サブバンドフィルタ構造から最低のサブバンドを「抽出し」、図８に示されるように別々にそれを処理することによって、図７ｂの構造を修正するときにさらに改善されてよい。
【００５５】
図７ｂの例と対照的に、バルク遅延DＢＵＬＫは、（通常は最低のサブバンドの遅延D１である）最短のサブバンド遅延に基づいて形成されるのではなく、むしろ例えばD２のような第２の最短の遅延に基づいて形成される。遅延D３−D２、D４−D２、 D５−D２等の間の差は図７ｂの例においてよりはるかに小さいので、必要とされるメモリの達成可能な削減はより大きくなる。したがって、広帯域入力信号ｘ［ｎ］は、フィルタバンクプロトタイプフィルタ２２‘を称してフィルタリングされ、次いで係数Ｎ（デシメータ２５）でダウンサンプリングされ 、係数２Ｍ（利得２６）で乗算され、最後に元のサブバンド遅延D１で遅延され、第１のサブバンドの出力信号ｙ１［ｎ］を得るために対応する（実数値の）サブバンドＦＩＲフィルタでフィルタリングされる。残りのサブバンド出力信号ｙ２［ｎ］， ｙ３［ｎ］， ．．．， ｙＭ／２＋１［ｎ］は図７ｂの例においてと同様に取得される。唯一の相違点は、図７ｂの例においてよりも大きく、相応してより短いサブバンド遅延D３−D２， D４−D２， D５−D２を生じさせるバルク遅延DＢＵＬＫ，のサイズである。全てのサブバンド出力信号ｙｍ［ｎ］は、前の例と同じように広帯域出力信号ｙ［ｎ］に結合される。
【００５６】
ＦＩＲフィルタ２２‘（つまり、ＡＦＢプロトタイプフィルタＰ（ｚ））を実装するための計算作業、及び以後のダウンサンプリングを削減するために、ＦＩＲフィルタ２２’及びデシメータ２５は、計算上より効率的である対応する多相フィルタによって置き換えられてよい（図９を参照すること）。図９の図では、フィルタ転送関数Ｐ０（ｚ）、Ｐ１（ｚ）等はフィルタ係数 ［ｐ０， ｐＮ， ｐ２Ｎ， ．．．］、［ｐ１， ｐＮ＋１， ｐ２Ｎ＋１， ．．．］等を含むフィルタを表し、それによって係数ｐ０、ｐ１、 ｐ２，等がＡＦＢプロトタイプフィルタＰ（ｚ）のフィルタ係数を表している。
【００５７】
図８又は図９の例を使用すると、本例（１７のサブバンド、１０のチャネル）では３４０回の初期化から２０回の初期化へ、メモリ転送初期化の大幅な削減が可能になる。これらは、全てのチャネルのバルク遅延をスワップし、一度にスワップし直すことができるように、バルク遅延を等化することによってさらに２に削減されてよい。この手順は図１０に示される。
【００５８】
本発明及びその優位点は詳細に説明されてきたが、添付特許請求の範囲によって定められる本発明の精神及び範囲から逸脱することなく、多様な変更、置換、及び改変を本明細書に加えることができることが理解されるべきである。
【００５９】
さらに、本発明の範囲は、本明細書に説明されるプロセス、機械、製造、物質の組成、手段、方法、及びステップの特定の実施形態に制限されることが意図されていない。統合者は、本発明の開示から、本明細書に説明される対応する実施形態と実質的に同じ機能を実行する、又は実質的に同じ結果を達成する、現在存在している、又は後に開発されるプロセス、機械、製造、物質の組成、手段、方法、及びステップが本発明に従って活用され得ることを容易に理解する。したがって、添付特許請求の範囲は、それらの範囲内で、かかるプロセス、機械、製造、物質の組成、手段、方法、及びステップを含むことが意図される。
【技術分野】
【０００１】
本発明は、音声アプリケーションでのサブバンド信号処理に関し、特に、所望されるターゲット転送関数の計算上効率的な実装を提供するサブバンドＦＩＲフィルタのセットを設計する方法に関する。
【背景技術】
【０００２】
ＩＩＲフィルタ（無限インパルス応答フィルタ）と対照的に、ＦＩＲフィルタ（融点インパルス応答フィルタ）は、所望される転送関数（つまり、大きさ及び位相応答）を有するデジタルフィルタを実現する可能性を提供する。したがって、所望される転送関数は、最小位相となるように設計できるだけではなく、最大位相、又は混合相となるように設計することもできる。さらに、しばしば音声信号処理で所望される、線形位相転送関数も実装できる。
【発明の概要】
【発明が解決しようとする課題】
【０００３】
しかしながら、特に音声信号処理では、動作中のメモリ要件だけではなく計算上の作業も増加させる、ＦＩＲフィルタを使用するときには、必要とされるフィルタ長は、しばしばかなり長くなる。これは、おもに２つの物理的な要因による。第１に、大部分の音声用途に存在する室内インパルス応答の減衰期間はかなり長く、その結果、ＦＩＲフィルタを使用するときには、相応して長いフィルタ長につながる。第２に、人間の聴覚系は、異なる周波数バンドで不均一な周波数分解能を提供するように適応されている。人間の聴覚系は、低周波はきわめてよく分解する。つまり、周波数差は低い絶対周波数でよく認識されるのに対し、高周波はそれほど容易に識別されない。例えば、両方の場合とも周波数差は１００Ｈｚであるが、１００Ｈｚのトーンが２００Ｈｚのトーンから容易に区別され得るのに対し、人間の耳では５１００Ｈｚトーンから５０００Ｈｚのトーンを区別するのが困難である。すなわち、人間の聴覚系の周波数分解能は周波数が高くなるにつれて減少する。この現象は周知であり、バーク尺度、メル尺度、及びＥＲＢ（同等矩形帯域幅）尺度等の人間の聴覚系に適応された心理音響周波数尺度を形成する。
【０００４】
研究は、内装（カーペット、布張りをした家具等）の結果として、エネルギーの分解が低速であるため、試聴室の室内インパルス応答が、特に低周波で著しく長いことを示している。音声再生システムによって生じる音圧がバス周波数範囲（つまり、２００Ｈｚ未満）で最大であるのに対し、人間の聴覚系は低周波音声信号により敏感でなくなるという事実によって、この影響はさらに強くなっている。
【０００５】
これら全ての要因を統合した見解は、典型的な試聴室の特性だけではなく、人間の聴覚系の特性によって、音声信号処理システムで十分な音声品質を提供するためには、ＦＩＲフィルタの長さが特定の最小長さよりも短くないことが要求される。例えば、バス周波数範囲で約１０Ｈｚの必須周波数分解能を提供するには、４４１００Ｈｚのサンプリング周波数での音声チャネルごとに、４４１０フィルタ係数のＦＩＲフィルタが必要となり、それによって現代的な音声システムは最大１６のチャネルを有する。
【０００６】
長いＦＩＲフィルタが動作中に多大な計算作業及び／又は高メモリ要件を伴うので、標準的なデジタル音声信号処理の使用を可能にするためには、音声アプリケーションでのＦＩＲフィルタの効率的な実装に対するニーズがある。
【課題を解決するための手段】
【０００７】
少なくとも１つの音声チャネルのある音声システムが開示される。音声システムは、チャネルごとに少なくとも１つのデジタルフィルタが実装されるデジタル音声プロセッサを含む。各チャネルのデジタルフィルタは、広帯域入力音声信号を受信し、その入力音声信号を複数のサブバンドに分割し、このようにして等しい帯域幅を有するサブバンド信号を提供し、サブバンド信号のスペクトルが入力音声信号のスペクトルを作るように構成される解析フィルタバンクと、各サブバンド用であり、このようにしてそれぞれフィルタリングされたサブバンド信号を提供するサブバンドＦＩＲフィルタと、フィルタリングされたサブバンド信号を受信し、１つの広帯域出力音声信号を提供するためにそれらを結合するように構成される合成フィルタバンクとを含み、遅延は各サブバンド信号に関連付けられ、サブバンド信号の内の１つの遅延が解析フィルタバンクの上流にある広帯域入力音声信号に適用され、残留遅延が解析フィルタバンクの下流の残りのサブバンド信号に適用される。対応するフィルタリング方法が開示される。
【０００８】
本発明は、以下の図面及び説明を参照するとよりよく理解できる。図中、類似する参照番号は、対応する部分を示す。
【０００９】
例えば、本願発明は以下の項目を提供する。
（項目１）
少なくとも１つのデジタルフィルタが各チャネルに実装されるデジタル音声プロセッサを備える少なくとも１つの音声出力チャネルのある音声システムであって、各チャネルの前記デジタルフィルタが、
広帯域入力音声信号を受信し、前記受信した音声信号を複数のサブバンドに分割し、このようにして等しい帯域幅を有するサブバンド信号を提供するように構成される解析フィルタバンクであって、前記サブバンド信号の前記スペクトルが前記入力音声信号の前記スペクトルを構成する解析フィルタバンクと、
各サブバンド用のサブバンドＦＩＲフィルタであって、このようにしてそれぞれのフィルタリングされたサブバンド信号を提供するサブバンドＦＩＲフィルタと、
前記フィルタリングされたサブバンド信号を受信し、広帯域出力音声信号を提供するためにそれらを結合するように構成される合成フィルタバンクと、
を備え、
遅延が各サブバンド信号に関連付けられ、前記サブバンド信号の内の１つの前記遅延が前記解析フィルタバンクの上流の前記広帯域入力音声信号に適用され、前記残留遅延が前記解析フィルタバンクの下流の前記残りのサブバンド信号に適用される、
音声システム。
（項目２）
前記遅延されていない入力音声信号を受信し、追加サブバンド信号として、前記解析フィルタバンクでタッピングされていない最低周波数サブバンドの前記サブバンド信号を提供する低域通過フィルタと、
前記それぞれのサブバンドＦＩＲフィルタの上流又は加硫で前記最低周波数サブバンド信号に適用される遅延と、
をさらに備え、
前記広帯域入力音声信号に適用される前記遅延が、前記解析フィルタバンクでタッピングされる前記サブバンドと関連付けられる最短の遅延に一致する、
上記の項目に記載のシステム。
（項目３）
前記サブバンドＦＩＲフィルタのフィルタ係数の前記数が、銭サブバンドフィルタの前記周波数分解能が、少なくともほぼ前記人間の耳の前記周波数分解能に相当するように設定される、
上記項目のいずれか１つの項目に記載のシステム。
（項目４）
少なくとも第１の音声チャネルと、第２の音声チャネルとを備え、
各広帯域入力音声信号が、前記解析フィルタバンクでタッピングされる全てのチャネルの全てのサブバンド信号の最短の遅延によって与えられる等しい遅延で遅延される、
上記項目のいずれか１つの項目に記載のシステム。
（項目５）
前記解析フィルタバンク（複数の場合がある）及び前記合成フィルタバンク（複数の場合がある）が、前記入力音声信号をブロック単位で処理する偶数で積み重ねられるＧＤＴＦフィルタであり、前記サブバンド信号が、前記ブロックサイズに対応する係数で、前記広帯域信号に比較してダウンサンプリングされ、
前記処理されたサブバンド信号の前記最低周波数及び前記最高周波数のサブバンド信号が実数値であり、前記他のサブバンド信号が複素数値である、
上記項目のいずれか１つの項目に記載のシステム。
（項目６）
前記低域通過フィルタが、前記解析フィルタバンク（複数の場合がある）によって使用されるプロトタイプフィルタに相当する、上記項目のいずれか１つの項目に記載のシステム。
（項目７）
音声システムの少なくとも１つの音声出力チャネル内に提供される少なくとも１つの音声信号をフィルタリングするための方法であって、各チャネルのために、
所与の周波数バンド内にあるスペクトルを有する、フィルタリングされる広帯域入力音声信号を受信することと、
前記入力音声信号の前記周波数バンドを、複数のサブバンドに分割し、このようにしてサブバンド信号を提供することであって、前記サブバンド信号の前記スペクトルが前記入力音声信号の前記スペクトルを構成する、提供することと、
各サブバンド信号をフィルタリングし、このようにしてサブバンドごとにそれぞれのＦＩＲ又はＩＩＲフィルタを使用してそれぞれのフィルタリングされたサブバンド信号を提供することと、
広帯域出力音声信号を提供するために、前記フィルタリングされたサブバンド信号を結合することと、
を含み、
遅延が各サブバンドに関連付けられ、前記サブバンドの内の１つの前記遅延が、前記周波数バンドを分割する前に前記広帯域入力音声信号に適用され、前記残留遅延が前記周波数バンドを分割後に前記残りのサブバンド信号に適用される、
方法。
（項目８）
前記周波数バンドを分割することが、前記入力信号をブロック単位で処理する第１の偶数で積み重ねられたＧＤＴＦ解析フィルタバンクを使用して実行され、
前記周波数バンドを分割することが、前記広帯域信号と比較して、前記ブロックサイズに相当する係数分、前記サブバンド信号を前記ダウンサンプリングすることをさらに含み、
前記フィルタリングされたサブバンド信号を結合することが、前記第１のフィルタバンクに対応し、前記広帯域出力信号をブロック単位で提供する第２の偶数で積み重ねられたＧＤＴＦ合成フィルタバンクを使用して実行される、
上記の項目に記載の方法。
（項目９）
前記遅延していない入力音声信号を低域通過フィルタリングし、追加サブバンド信号として、前記最低周波数サブバンドの前記サブバンド信号を提供することと、
それぞれのＦＩＲ又はＩＩＲフィルタを使用して前記最低周波数サブバンド信号をフィルタリングし、このようにして対応するフィルタリングされたサブバンド信号を提供することと、
ＦＩＲフィルタリングの前又は後に前記最低周波数サブバンド信号を表示することと、
をさらに含む、上記項目のいずれか１つの項目に記載の方法。
（項目１０）
前記サブバンドＦＩＲフィルタのフィルタ係数の前記数が、前記サブバンドフィルタの前記周波数分解能が前記人間の耳の前記周波数分解能に少なくともほぼ一致するように設定される、上記項目のいずれか１つの項目に記載の方法。
（項目１１）
前記音声システムが、第１の音声出力チャネル及び第２の音声出力チャネルを備え、
各広帯域入力音声信号が、前記周波数バンドを除算することによって得られる全てのチャネルの全てのサブバンド信号の最短遅延によって与えられる等しい遅延で遅延される、
上記項目のいずれか１つの項目に記載の方法。
（項目１２）
前記低域通過フィルタリングが、低域通過フィルタとして、前記周波数バンドを分割するために使用される前記ＧＤＦＴフィルタバンク（複数の場合がある）によって使用されるプロトタイプフィルタに相当するフィルタを使用することを含む、
上記項目のいずれか１つの項目に記載の方法。
【００１０】
（摘要）
少なくとも１つの音声チャネルのある音声システムが開示される。音声システムは、チャネルごとに少なくとも１つのデジタルフィルタが実装されるデジタル音声プロセッサを含む。各チャネルのデジタルフィルタは、広帯域入力音声信号を受信し、複数のサブバンドに入力音声信号を分割し、このようにして等しい帯域幅を有するサブバンドを提供するように構成される解析フィルタバンクであって、サブバンドの信号のスペクトルが入力音声信号のスペクトルを構成する解析フィルタバンクと、各サブバンド用であり、このようにしてそれぞれのフィルタリングされたサブバンド信号を提供するサブバンドＦＩＲフィルタと、フィルタリングされたサブバンド信号を受信し、それらを結合し、広帯域出力音声信号を提供するように構成される合成フィルタバンクとを含み、遅延が各サブバンド信号と関連付けられ、サブバンド信号の内の１つの遅延が解析フィルタバンクの上流で広帯域入力音声信号に適用され、残留遅延が解析フィルタバンクの下流で残りのサブバンド信号に適用される。
【図面の簡単な説明】
【００１１】
【図１】ＦＩＲフィルタのフィルタ係数の適応計算のための基本信号処理構造を示し、ＦＩＲフィルタは所定のターゲット関数Ｐ（ｚ）にほぼ一致する転送関数Ｇ（ｚ）を表す。
【図２】図１の構造の変型を示し、それによって図１の全帯域ＦＩＲフィルタＧ（ｚ）がサブバンドＦＩＲフィルタＧｍ（ｚ）のセットで置き換えられる。
【図３】全帯域フィルタＧ（ｚ）に代わることができるサブバンドＦＩＲフィルタＧｍ（ｚ）のセットを含む信号処理構造を示す図である。
【図４】それぞれがサブバンドＦＩＲフィルタ及び遅延回線を含む、多くのサブバンドを含む適応ＦＩＲフィルタ設計のための信号処理構造を示す図である。
【図５】サブバンドＦＩＲフィルタの長さを適応するときに低周波サブバンドを強調するためのバーク尺度に基づいた重み付け係数を示す図である。
【図６】図偶数で積み重ねられたフィルタバンク及び奇数で積み重ねられたフィルタバンクの規模応答を概略で示す図である。
【図７ａ】複数のサブバンドフィルタ内で分割されるＦＩＲフィルタの構造を概略で示す図である。
【図７ｂ】図７ａの例の簡略化された代替策を概略で示す図である。７ａ
【図８】メモリ要件が削減された図６のフィルタバンクの代替策を示す図である。
【図９】図８の例を改善するためのオプションを示す図である。
【図１０】２つの別々の音声チャネル内のフィルタを含む例を示す図である。
【発明を実施するための形態】
【００１２】
図１は、ＦＩＲフィルタ２０のフィルタ係数ｇｋ （ｋ ＝ ０， １， ．．．， Ｋ−１）の適応計算のための基本信号処理構造を示し、それによって下付き文字ｋはフィルタ係数の指数を示し、Ｋはフィルタ長を示す。ＦＩＲフィルタ２０は、フィルタ係数ｇｋの適応が成功した後に、基準システム１０の所定のターゲット関数Ｐ（ｚ）にほぼ一致する（離散）転送関数Ｇ（ｚ）を有する。適応フィルタ設計手順を実行するために、基準フィルタ１０及びＦＩＲフィルタ２０には、例えば、ホワイトノイズ又はＦＩＲターゲット転送関数Ｐ（ｚ）の通過帯域を含む帯域幅を有する任意の他の信号である、試験信号（入力信号ｘ［ｎ］）が供給される。ＦＩＲフィルタ２０の出力信号ｙ［ｎ］は、基準システム１０の出力信号から、つまり所望される信号ｄ［ｎ］から減算される（減算器３０）。差ｄ［ｎ］−ｙ［ｎ］は、エラー信号ｅ［ｎ］として使用され、適応部２１に供給される。適応部２１は、各サンプル時間間隔中に、エラー信号及び入力信号ｘ［ｎ］（本文脈では基準信号としても示される）からＦＩＲフィルタ係数ｇｋの更新されたセットを計算するように構成される。最小二乗平均（ＬＭＳ）アルゴリズム又は正規化最小二乗平均（ＮＬＭＳ）アルゴリズムが、フィルタ係数の適応のために利用され得る。ただし、この目的のためには、さまざまな適応アルゴリズムが活用されてもよい。適応アルゴリズムの収束後、ＦＩＲフィルタ係数ｇｋは、ターゲット転送関数Ｐ（ｚ）の最適近似である転送関数Ｇ（ｚ）を表す。
【００１３】
ＦＩＲフィルタを使用するときの計算作業を削減するための１つのオプションは、フィルタリングされる信号のスペクトルを多くの狭帯域信号（サブバンド信号）に分割し、各狭帯域信号を別々にフィルタリングすることである。全帯域信号のいくつかのサブバンド信号への分割は、いわゆる解析フィルタバンク（ＡＦＢ）によって実現され得る。同様に、サブバンド信号は、対応する合成フィルタバンク（ＳＦＢ）を用いて単一の全帯域信号に（再）結合され得る。以下では、全帯域信号は、例えば所望される信号ｄ［ｎ］のように下付き文字なしで示され、ｎは時間インデックスである。さらに、例えばｄｍ［ｎ］のような下付き文字を有する信号は、対応する全帯域信号ｄ［ｎ］の分解であるサブバンド信号のセットを示す。それによって、下付き文字ｍは、サブバンドの数（ｍ ＝ １， ２， ．．．， Ｍ）を示す。同様に、離散全帯域転送関数Ｇ（ｚ）は、多くのサブバンド転送関数Ｇｍ（ｚ）に分解され得る。
【００１４】
図２は、図１の基本信号処理構造に類似した、それによって適応ＦＩＲフィルタ２０が、サブバンドＦＩＲフィルタのセット２０’で置き換えられる基本信号処理構造を示す。この目的のため、全帯域入力信号ｘ［ｎ］は、ＡＦＢ２２を使用することによって数Ｍのサブバンド入力信号ｘｍ［ｎ］（ｍ ＝ １， ２， ．．．， Ｍである）に分割される。同様に、所望される全帯域信号ｄ［ｎ］は、ＡＦＢ１１を使用して（再び、ｍ ＝ １， ２， ．．．， Ｍ）数Ｍのサブバンド信号ｄｍ［ｎ］に分けられる。各サブバンドＦＩＲフィルタは、狭帯域転送関数Ｇｍ（ｚ）を実現し、下付き文字ｍは再びサブバンドの数を示す。各サブバンドフィルタＧｍ（ｚ）は、そのフィルタ係数ｇｍｋによってあらわされてもよく、それによってｋは、ｋ＝０からｋ＝Ｋｍ−１（Ｋｍは、ｍ番目のサブバンド内のフィルタＧｍ（ｚ）のフィルタ長さである）に及ぶフィルタ係数のインデックスを示す。各ＦＩＲフィルタＧｍ（ｚ）は、対応するエラー信号ｅｍ［ｎ］ ＝ ｄｍ［ｎ］ − ｙｍ［ｎ］を受信し、それぞれのサブバンドｍごとに更新されたフィルタ係数ｇｍｋ （ｋ ＝ １， ２， ．．．， Ｋｍ−１） のそれぞれのセットを計算する、適応部と関連付けられている（適応部のセットは図２で数字２１’で示される）。
【００１５】
図１の例と同様に、Ｍ個のサブバンドＦＩＲフィルタＧｍ（ｚ）の各１つのフィルタ係数ｇｍｋは、適応アルゴリズムの収束後に、全てのサブバンド転送関数Ｇｍ（ｚ）の組み合わせから生じる総転送特性が所定のターゲット関数Ｐ（ｚ）に一致するように適応される。
【００１６】
適切なフィルタ係数ｇｍｋの計算後、ＦＩＲフィルタＧｍ（ｚ）のセット２０’は、音声信号をフィルタリングするための図３に示されるように、解析フィルタバンク（ＡＦＢ２２）と対応する合成フィルタバンク（ＳＦＢ２２’）との間で操作されてよい。この場合、ＡＦＢ２２、ＦＩＲフィルタバンク２０’、及びＳＦＢ２２’は、例えば音声システムの等化フィルタを表すことがある、ターゲット関数Ｐ（ｚ）にほぼ一致する転送関数Ｇ（Ｚ）をともに実装する。現代の音声システムでのように、リスナーにとって所望されるサウンドインプレッションを生じさせるために、大きさだけではなく位相が等化にさらされる。したがって、ターゲット関数Ｐ（ｚ）は、一般的に、非線形位相特性のある非最小位相フィルタを表す。
【００１７】
解析フィルタバンクと、サブバンドフィルタのセットと、合成フィルタバンクとを含む、図３の信号処理構造が、（計算作業及びメモリ要件という点で）同じ転送関数の「通常の」ＦＩＲフィルタよりもより効率的であるのかという質問に対する回答は、とりわけ解析フィルタバンク及び合成フィルタバンクの効率的な実装形態の可用性に依存している。人間の聴覚系の不均一な周波数分解能を説明するために、低い中心周波数があるサブバンドの帯域幅が、より高い中心周波数のあるサブバンドの帯域幅よりも狭いフィルタバンクを使用することが望ましいだろう。かかる心理音響的に動機付けされる全帯域信号の、帯域幅が可聴周波数範囲内のそれぞれのサブバンドの位置に依存するサブバンド信号のセットへの分割を実現するいくつかの手法が存在する。しかしながら、異なる帯域幅のサブバンドのセットに１つの入力スペクトルを不均一に分割することを可能にする効率的なフィルタバンクは既知ではない。それにも関わらず、技術では、全帯域信号の、等しい帯域幅のサブバンド信号のセットへの分割を可能にする他の方法が既知である。一例は、単一のプロトタイプフィルタ、及びほぼあらゆる信号処理環境で利用できるＦＦＴアルゴリズムを利用する、Ｓ．Ｗｅｉｓｓらによって説明されるオーバーサンプリングされたＧＤＦＴフィルタバンクの高速実装形態である（ＩＥＥ Ｅｌｅｃｔｒｏｎｉｃｓ Ｌｅｔｔｅｒｓの第３６巻、１５０２−１５０３ページ、２０００年のＳ． Ｗｅｉｓｓ、「Ｆａｓｔ Ｉｍｐｌｅｍｅｎｔａｔｉｏｎ ｏｆ Ｏｖｅｒｓａｍｐｌｅｄ Ｍｏｄｕｌａｔｅｄ Ｆｉｌｔｅｒ Ｂａｎｋｓ」を参照すること）。
【００１８】
上記を要約すると、効率的な実装形態は不均一な帯域幅のサブバンドを処理するために利用できないので、等しい帯域幅のサブバンドを用いて動作するフィルタバンクを使用する必要がある。ただし、等しく幅広いサブバンドの機能不全を緩和するために、それぞれのサブバンドに割り当てられるＦＩＲフィルタの異なるフィルタ長が選ばれてよい。すなわち、ＦＩＲフィルタは、高周波分解能が要求されるサブバンドにおいてよりも、低周波分解能が要求されるサブバンド内でより少ないフィルタ係数を含む。高周波分解能が要求されるサブバンドは、通常、、可聴周波数範囲の下部にあるサブバンドである。したがって、人間の聴覚系の周波数分解能に相当する周波数分解能は、等しく幅広いサブバンドを用いて動作する効率的なフィルタバンクを使用することによって達成され得る。
【００１９】
上述されたように、ターゲット関数Ｐ（ｚ）は、一般に、周波数上での非線形グループ遅延特性を有する非最小位相フィルタである。異なるサブバンドでの異なるグループ遅延から生じる異なる信号伝搬遅延を補償するために、遅延回線は、遅延回線の上流又は下流にある各サブバンドＦＩＲフィルタに接続されてよい。したがって、遅延等化は、計算上、非効率的である追加のＦＩＲフィルタ係数を使用して達成される必要はない。フィルタ係数の数だけではなく、遅延値も実現されるターゲット転送関数Ｐ（ｚ）（つまり、大きさ及び位相応答）に依存するので、フィルタ係数の数（つまり、フィルタ「タップ」の数、及び遅延値は、「適応タップ割り当て」アルゴリズム及び「適応遅延割り当て」アルゴリズムを使用して、本明細書中でいかに説明されるように、サブバンドごとに適応的に決定され得る。その結果、Ｍ個のサブバンドＦＩＲフィルタＧｍ（ｚ）を設計するときには、フィルタ係数（図１の係数ｇｍｋを参照すること）だけではなく、係数の数ＫＭ及び追加遅延△ｍも適応的に決定される。図４は、サブバンドＦＩＲフィルタ設計に使用される信号処理構造を示す。図４の例は、図２の構造の改良版であり、各サブバンド信号経路に追加の遅延、及びグローバル「適応割り当て及び遅延割り当て部」を備えている。４は、図２の構造の改良版を示し、４の例は、図２の構造の改良版であり、各サブバンド信号経路に追加の遅延、及びグローバル「適応割り当て及び遅延割り当て部」を備えている。
【００２０】
図２の例でのように、全体入力信号ｘ［ｎ］（例えば、帯域が制限されたホワイトノイズ）が、ターゲット転送関数うＰ（ｚ）のあるシステム１０に供給され、このようにして所望される信号ｄ［ｎ］を生成する。入力信号ｘ［ｎ］だけではなく、所望される信号ｄ［ｎ］も、それぞれ数Ｍのサブバンド信号ｄｍ［ｎ］ 及び ｘｍ［ｎ］に分割される。図４の例は、最初のサブバンド及び最後のサブバンド（ｍ ＝ １及びｍ ＝ Ｍ）と関連付けられた構成要素及び信号だけを示す。サブバンド入力信号ｘｍ［ｎ］は、転送関数Ｇｍ（ｚ）がある適応ＦＩＲフィルタに供給され、このようにしてフィルタリングされたサブバンド信号ｙｍ［ｎ］を生成する。フィルタリングされた各サブバンド信号ｙｍ［ｎ］は、対応する所望される信号ｄｍ［ｎ］から差し引かれ、サブバンドごとにエラー信号ｅｍ［ｎ］を生じさせる。適応部は、それぞれのＦＩＲフィルタＧｍ（ｚ）のフィルタ係数ｇｍｋ（つまり、フィルタのインパルス応答｛ｇｍ０， ｇｍ１， ．．．， ｇｍ（Ｋ−１）｝）を最適化するための各ＦＩＲフィルタＧｍ（ｚ）に割り当てられ、それによってフィルタ係数ｇｍｋの最適セットがそれぞれのエラー信号ｅｍ［ｎ］のノルム（例えば、累乗）を最小限に抑える。
【００２１】
遅延Dｍを提供する遅延回線は、各サブバンドＦＩＲフィルタＧｍ（ｚ）に上流で又は下流で接続される。さらに、以下にさらに詳細に説明される適応タップ割り当て及び適応遅延割り当てアルゴリズムに従って、遅延回線Dｍの対応する遅延値だけではなく、ＦＩＲフィルタＧｍ（ｋ）のフィルタ長Ｋｍも動的に適応するように構成される、「適応タップ割り当て及び適応遅延割り当て部」４０が設けられる。
【００２２】
上述された適応タップ割り当て（つまり、ＦＩＲフィルタ長の適応）にはさまざまな手法が考えられ得る。１つの手法は、総エラー信号ｅ［ｎ］（それによって、ｅ［ｎ］ ＝ ｅ１［ｎ］ ＋ ｅ２［ｎ］ ＋ ．．．＋ ｅＭ［ｎ］）が最小値に達するまで、サブバンドＦＩＲフィルタＧｍ（ｚ）のフィルタ長Ｋｍを変えることである。実際には、かかる方法は優れた結果を出すが、、フィルタ係数の数の各変更後に、適応フィルタが再び収束する時間を要するので、きわめて時間がかかる。上述の手法と同様に優れた結果を出すが、要する時間がはるかに少ない別の手法は、Ｓ個の一番端のフィルタ係数ｇｍ（Ｋｍ−１）， ｇｍ（Ｋｍ−２）， ．．．， ｇｍ（Ｋｍ−Ｓ）のエネルギーを検討する。サブバンドのフィルタＧｍ（ｚ）のフィルタ長Ｋｍは、上述のＳ個の一番端のフィルタ係数のエネルギーがほぼ等しくなるまで変えられる。この手法は、サブバンドフィルタが経時的に急激に減衰することを必要とし、それは現実世界のシステムではつねに当てはまらなくてはならない。各サブバンドフィルタのＳ個の一番端のフィルタタップのエネルギーを比較することによって、サブバンドフィルタＧｍ（ｚ）がどれほどよくターゲット関数Ｐ（ｚ）を近似し、信号減衰挙動がターゲット関数Ｐ（ｚ）のインパルス応答の信号減衰挙動に似ているサブバンドフィルタインパルス応答を達成するように、サブバンドフィルタ全体でフィルタ係数を分散し直すためのガイドラインとなり、最小限に抑えることができるエラーに関して最適とみなすことができる。
【００２３】
以下に、適応タップ割り当てアルゴリズムの例がさらに詳しく説明される。他のＭ／２サブバンド信号は、第１のＭ／２サブバンドの信号の共役複素数のコピーであるため、実数値の全帯域入力信号ｘ［ｎ］の場合（図２を参照すること）、Ｍ／２サブバンドだけを処理する必要があることに留意するべきである。したがって、それぞれのサブバンドＦＩＲフィルタ転送関数は、以下の関係性に準拠し、
【００２４】
【数１】
それによって、アスタリスクは共役複素数演算子を示す。
【００２５】
したがってサブバンドフィルタＧｍ（ｚ）のフィルタ長Ｋｍ（それによって、ｍ ＝ １， ２， ．．．， Ｍ／２）は、Ｑ個のサンプル（つまり、サブバンドシステム内のサンプル）の期間で修正される。しかしながら、全てのサブバンドフィルタＧｍ（ｚ）のフィルタ係数の総数ｇｍｋ［ｎ］は一定のままである。すなわち、１つ又は複数のサブバンドフィルタのフィルタ長が増加する場合、別のサブバンドフィルタのフィルタ長は、フィルタ係数の総数を一定に保つように削減されなければならない。したがって、Ｑ個のサンプルの期間では、Ｍ／２サブバンドＦＩＲフィルタのそれぞれの長さは、ΔＫ 個の係数分削減される。したがって、さらに以下に説明される特定の基準に従って、Ｍ個のサブバンドフィルタ全体で分散し直されるΔＫ・Ｍ／２個の「自由な」係数がある。
【００２６】
上記の「分散し直し」は、以下の方程式で表されてよい。
【００２７】
【数２】
上式ではｍ ＝ １， ２， ．．．であり、Ｍ／２はサブバンドの数を示す。式 ｃｍ［ｎ／Ｑ］は、フィルタタップ（つまり、フィルタ係数）の分散の上述された基準を表す。上述されたように、１つの有用な基準は、サブバンドエラー信号ｅｍ［ｎ］のエネルギーである。この場合、ｃｍ［ｎ／Ｑ］は以下のように表すことができる。
【００２８】
【数３】
上式では、ｍ ＝ １， ２， ．．．であり、Ｍ／２及びＲはエラー信号が平均化されるサンプルの数である。ただし、適応サブバンドＦＩＲフィルタは、方程式（３）が評価される前に集束しなければならず、結果的に、ＲはＱよりもはるかに小さくなる必要がある（つまり、Ｒ ＜＜ Ｑ）。別の基準は、それぞれのサブバンドＦＩＲフィルタの一番端のＳ個のフィルタ係数のエネルギーを検討する。この場合、ｃｍ［ｎ／Ｑ］ は以下として表すことができ、
【００２９】
【数４】
上式ではｍ ＝ １， ２， ．．．であり、Ｍ／２及びＫｍはそれぞれのサブバンドでの現在のフィルタ長Ｋｍ［ｎ］である。代わりに、サブバンド入力信号ｘｍ［ｎ］のｎエネルギーは、一番端のＳ個のフィルタ係数とともに検討され得る（方程式５を参照すること）。この場合、ｃｍ［ｎ／Ｑ］は、以下として表すことができる。
【００３０】
【数５】
上述されたように、方程式（３）に従った基準は最善の結果を生じさせるが、評価するのに時間がかかる。方程式（５）に従った基準は、ＡＥＣ（音響エコーキャンセル）システムの場合でのように、ターゲットシステムが時間的に変化し、入力信号が任意に着色されるときに使用されるべきである。本ケースでは、つまり入力信号ｘ［ｎ］がホワイトノイズとなるように設計者によって選択されることがあるサブバンドＦＩＲフィルタ設計の場合、方程式（４）は優れた品質結果を出すと同時に、高速適応を可能にする。
【００３１】
心理音響的な態様を説明するために、方程式（３）、（４）、又は（５）に定義される式ｃｍは、対応する重み付け係数ｗｍで重み付けされてよい。つまり、式 ｃｍ［ｎ／Ｑ］は、方程式（３）、（４）、又は（５）の ｗｍ・ｃｍ［ｎ／Ｑ］で置き換えられる。重み付け係数ｗｍは、人間の聴覚系の周波数分解能が検討されるように選ばれる必要がある。バーク尺度を使用すると、係数ｗｍ は、以下のように計算でき、
【００３２】
【数６】
上式では、ｆｃ，ｍ は、ｆｃ，ｍ ＝ （２ｍ−１）・ｆＳ／（２・Ｍ）として計算されてよく、ｆＳ がＨｚ単位のサンプリング周波数であるｍ番目のサブバンドの（Ｈｚ単位の）中心周波数を示す。方程式（６）に従って計算された正規化重み付け係数ｗｍの例のセットは、図５に示され、それによって多くの１６のサブバンドが検討された（Ｍ／２＝１６）。
【００３３】
上述されたように、方程式（２）に従った係数の分散し直しは、サブバンドエラー信号ｅｍ［ｎ］（方程式（３）と比較すること）のエネルギー又はサブバンドフィルタ係数ｇｍｋ（ｋ ＝ ｛Ｋｍ−Ｓ， ．．．， Ｋｍ−１｝であり、方程式（４）と比較すること）の部分（例えば、一番端のＳ個のフィルタ係数）のエネルギー、又はサブバンドフィルタ係数ｇｍｋ及びサブバンド入力信号ｘｍ［ｎ］ （方程式（５）と比較すること）のエネルギーの組み合わせを「等化する」。したがって、異なるエネルギーを元にした基準が考えられる。アルゴリズムの「目標」は、サブバンド全体で信号エネルギーを等しく分散することである。この目標は最適化プロセスの所望される結果であるので、上述されたアルゴリズムも最小化（又は最大化）タスクとして見ることができる。
【００３４】
理想的な場合、上述された基準ｃｍ［ｎ／Ｑ］（方程式（３）から（５）を参照すること）は、各サブバンドＭで等しい値を有し、それによって以下となる。
【００３５】
【数７】
すなわち、サブバンドと見なされるＭ／２の内のそれぞれ１つで、基準ｃｍ［ｎ］は、全てのサブセット上で総計された ｃｍ［ｎ］という合計値よりも小さい係数２／Ｍである。上述された最小化タスクに戻ると、最小限に抑えられる品質関数は、以下のように書くことができる。
【００３６】
【数８】
上述されたように、ＦＩＲフィルタＧｍ（ｚ）は、一般的に、異なる信号伝搬遅延を生じさせる異なるサブバンドでの異なるグループ遅延を有する。かかる異なる遅延を補償するために、各ＦＩＲフィルタは各サブバンドでフィルタの上流で又は加硫で調整可能な遅延回線に接続される。再び、サブバンドの遅延値は、「適応遅延割り当て」として以下に説明されるように反復して最適化できる。
【００３７】
「適応タップ割り当て」として上述されたサブバンド内のＦＩＲフィルタに適応的に割り当てられるＦＩＲフィルタ係数は、有限インパルス応答を生じさせる矩形ウィンドウ関数での乗算により切り捨てられた無限インパルス応答と見なすことができる。各サブバンドＦＩＲフィルタのために、この矩形ウィンドウ関数が（時間軸に沿って）シフトされてよく、それによってそれぞれの時間シフトは対応するＦＩＲフィルタに接続される遅延回線の有効遅延を表す。
【００３８】
本明細書に説明される適応遅延割り当てアルゴリズムは、各サブバンドに、それぞれのＦＩＲフィルタの最大エネルギー有限インパルス応答を生じさせる、遅延値（つまり、上述された矩形ウィンドウの時間シフト）を見つけることをターゲットにしている。すなわち、それぞれの遅延回線に遅延値を割り当てるとき、エネルギー又は対応するＦＩＲフィルタ係数だけが検討される。しかし、必ずしもこれが唯一の成功する手法であるわけはない。代替策として、遅延値は、総エネルギー信号ｅ［ｎ］ ＝ ｅ１［ｎ］ ＋ ｅ２［ｎ］ ＋ ．．．＋ ｅＭ／２［ｎ］のノルムが最小限に抑えられるように選ばれてよい。すでに適応タップ割り当てアルゴリズムに関して上述されたように、かかる基準を使用すると、十分に優れた結果が生じるが、収束は遅くなる。代わりに、遅延値は、それぞれのＦＩＲフィルタ係数の総エネルギーが最大となるように選ばれてよい。
【００３９】
しかしながら、サブバンドフィルタのＦＩＲインパルス応答（つまり、ＦＩＲフィルタ係数）は、一般的に、サブバンドへの変換に起因する原因因子と因果因子を含むので、最適遅延値を見つけるための他の基準はそれぞれのインパルス応答の原因部分及び因果部分のエネルギーを等化するようにターゲットされ得る。つまり、それぞれのサブバンドＦＩＲインパルス応答の第１の半分のエネルギーと第２の半分のエネルギーの差異は最小限に抑えられなければならない。かかる最適化戦略は、以下のように簡潔に要約されてよい。
【００４０】
サブバンドｍごとに、それぞれのサブバン路ＦＩＲインパルス応答ｇｍｋの第１の部分及び第２の部分を定義し、ＦＩＲインパルス応答の第１の部分及び第２の部分のエネルギーを計算する。Ｖ個のサンプルの期間で以下を繰り返す。
（Ａ） ＦＩＲインパルス応答 ｇｍｋのエネルギーが、インパルス応答の第２の部分のエネルギーよりも大きい場合には、それぞれの遅延回線Dｍの遅延をＰタップ（例えば、Ｐ＝１）分増加する。
（Ｂ） ＦＩＲインパルス応答ｇｍｋの第１の部分のエネルギーが、インパルス応答の第２の部分のエネルギーよりも低い場合には、それぞれの遅延回線DｍをＰタップ分（例えば、Ｐ＝１）減少する。ただし、最小遅延はゼロである。
（Ｃ）両方のエネルギーが等しい場合は、遅延を未変更のままにする。
【００４１】
ＦＩＲインパルス応答を等しい長さの２つの部分に分割したくなることがある。しかし、これは最大値がインパルス応答のほぼ真中に位置するインパルス応答に対してだけ良好な結果を生じさせる。左側にその最大値を有する原因インパルス応答を検討すると、上記の最適化戦略は、エネルギーが等しく分散されるまで、上述の矩形ウィンドウを右に時間シフトし、このようにしてインパルス応答の左側にゼロを挿入し、これは明らかに望ましくない。したがって、ＦＩＲインパルス応答の上述された第１の部分及び第２の部分がそれぞれのインパルス応答の全長をカバーするのではなく、変数Ｕによって定められ、０＜Ｕ＜１である、その上部マージン及び下部マージンだけをカバーする場合に、実践では、よりすぐれた結果が達成される。例を示すと、Ｕ＝０．６である場合は、エネルギーを計算するときにフィルタ係数の６０パーセントだけが検討され、それによって下部３０パーセントは上述の第１の部分を形成し、上部３０パーセントは対応する第２の部分を形成する。したがって、ｍ番目のサブバンドでは、それぞれ第１の部分及び第２の部分のエネルギーＥｍＬ ａｎｄ ＥｍＵは以下の通りである。
【００４２】
【数９】
【００４３】
【数１０】
したがって、ＥｍＵ ＜ ＥｍＬの場合、遅延値Dｍは増加され、ＥｍＵ ＞ ＥｍＬの場合減少され、ＥｍＵ ＝ ＥｍＬの場合未変更のままとされる、又は以下の数学用語で表される。
【００４４】
【数１１】
一般に、適応遅延割り当てアルゴリズムは、最小化タスクとして見ることができ、それによって、各サブバンドで、それぞれの遅延値Dｍが、品質基準 ｃ’ｍ［ｎ／Ｖ］ ＝ ｜ＥｍＵ−ＥｍＬ｜が最小値に達するまで変えられ、それによってインデックスｎ／Ｖは、基準がＶ番目のサンプル時間ごとしか評価されないことを示す。
【００４５】
タップの数分１つのサブバンドの遅延回線の遅延値Dｍを調整する（例えば、増加する）ときには、対応するフィルタインパルス応答（ＦＩＲフィルタ係数）が左にシフトされ、このようにして（遅延回線及びＦＩＲフィルタの遅延を含む）総遅延を一定に保つことに留意する必要がある。これによって、遅延値Dｍの変動に起因するＦＩＲフィルタ係数の長い再適応化に対するニーズが回避される。ただし、適応タップ割り当てだけではなく適応遅延割り当てが実行される適応期間は、ＬＭＳ／ＮＬＭＳアルゴリズムの適応ステップ幅よりもかなり長くなければならず、したがってＦＩＲフィルタ係数ｇｍｋの適応は、適応タップ割り当て及び適応遅延割り当てが開始されるときには定常状態にある。
【００４６】
各サブバンドの追加遅延の適応調整だけではなく、フィルタ長の適応調整も有用に組み合わせることができるが、それぞれの適応方法（適応タップ割り当て及び適応遅延割り当て）は、互いとは無関係に単独で使用され得る。例えば、フィルタ係数の数は（例えば、バーク尺度に従って）サブバンドごとに事前に定義されてよく、遅延 Dｍは、上述されたように適応的に調整されてよい。さらに、フィルタ長の適応決定は、別々の遅延部を使用しなくても、又は一定の所定の遅延を用いて、上述されたように適応タップ割り当て方法に従って形成されてよい。これは、特に、所望されるターゲット関数Ｐ（ｚ）が線形位相又は最小位相転送関数である場合に適切であり得る。
【００４７】
効率的なサブバンドＦＩＲフィルタ構造を設計するためのフィルタ設計方法を説明してきたが、かかるサブバンドＦＩＲフィルタ構造を効率的に実装するための例が、さらに詳しく説明される。第１に、解析フィルタバンク及び合成フィルタバンクが詳しく説明される。サブバンドＦＩＲフィルタを実装するための効率的な信号処理構造は、次いで、図７及び図８に関して説明される。
【００４８】
上述されたように、解析フィルタバンク及び対応する合成フィルタバンク（図３のフィルタバンク２２及び２２‘を比較すること）は、実数値の全帯域入力信号ｘ［ｎ］を複素数値サブバンド信号のセットｘｍ［ｎ］に分割し、複素数値のフィルタリングされたサブバンド信号ｙｍ［ｎ］を実数値の出力信号ｙ［ｎ］に再結合する。しかし、実数値サブバンド信号だけと動作する他のフィルタバンクも使用できる。この場合、サブバンドＦＩＲフィルタはより効率的に実装され得る。しかし、実数値フィルタバンクは、対応する複素数値フィルタバンクの２倍の多さのＭＡＣ（乗算―蓄積動作）を必要とする。さらに、サブバンドは複雑なサブバンド信号を使用するとき係数Ｍによってアンダーサンプリングされることがあるのに対し、対称的に、実数値フィルタバンクを使用するときは、最大アンダーサンプリング係数はＭ／２に過ぎない。最後に、サブバンドＦＩＲフィルタの必須フィルタ長は実数値サブバンド信号を使用するときに必要となるフィルタ長の半分に過ぎない。これは、複素数値ＡＦＢ及びＳＦＢの使用には、実数値フィルタバンクよりも少ない計算作業が必要になることを示している。しかし、両方のタイプのフィルタバンクとも、本発明での使用に適用できる。
【００４９】
複素数値フィルタバンクを使用するとき、偶数で積み重ねられたフィルタバンク又は奇数で積み重ねられたフィルタバンクのどちらかを選ぶことができる。両方のタイプとも、汎用化離散フーリエ変換アルゴリズム（ＧＤＦＴアルゴリズム）によって効率的に実装され得る。図６は、奇数に積み重ねられた（図６ｂ）フィルタバンクだけではなく偶数で積み重ねられた（図６ａ）のサブバンドも示す。一見して、Ｍ個のサブバンドの半分は、Ｍこのサブバンドの他方の半分の複素共役コピーであり、結果的に数Ｍ／２の複素数値サブバンド信号だけが処理されなければならないので、奇数で積み重ねられた（図６ｂ）のフィルタバンクの使用が好ましく見える。封数で積み重ねられたＧＤＦＴフィルタバンクを使用するとき、最低のサブバンド及び最高のサブバンドは実数値であるが、他のサブバンドの半分の帯域幅しか有していない。その結果、数Ｍ／２＋１のサブバンド信号が処理されなければならず、それによって最低のサブバンドの信号は実数値である。結果として、複素数値サブバンドと比較すると、メモリ要件は半分にすぎず、必須ＭＡＣは最低のサブバンドの４分の１にすぎない。最低のサブバンドのＦＩＲフィルタ （Ｇ１（ｚ））が、（例えば、バーク又は任意の他の心理音響的に適応された周波数尺度に従って周波数分解能と比較されるとき）最も多くのフィルタ係数ｇ１ｋ を有するという事実を考えると、処理されるサブバンドの総数はＭ／２からＭ／２＋１に増加するが、総効率は必要とされるメモリに関して２５パーセント向上し、必要とされるＭＡＣに関して５０パーセント向上する。
【００５０】
以下では、複数のサブバンドフィルタＧｍ（ｚ）及びそれぞれの遅延Dｍに分解された（所望される）デジタルＦＩＲフィルタの実装形態が、説明例として音声システムのアプリケーションを使用して説明される。図７ａに示される構造は、音声チャネルごとに実装されてよい。すなわち、各音声チャネルに解析フィルタバンク２２及び対応する合成フィルタバンク２２‘が提供される。この文脈では、「音声チャネル」は増幅器によって提供される出力チャネルを指し、入力信号の数は指さない（増幅器が２つの入力チャネル（ステレオ入力）しか有さないが、例えば５．１サラウンドサウンド出力には６つの出力チャネルを有してよいことに留意すること）。デジタルＦＩＲフィルタは、上述されたフィルタ設計方法に従って設計され得る。しかし、任意のＦＩＲフィルタは、本明細書に説明されるフィルタ構造及び方法を使用して実装されてよい。解析フィルタバンク２２（ｍ＝１，２，３、．．．、Ｍ）によって提供されるＭ個のサブバンド信号ｘｍ［ｎ］ のそれぞれは、（対応するサブバンドフィルタＧｍ（ｚ）を適用することによって）別々にフィルタリングされ、（対応する遅延Dｍを適用することによって）遅延される。上述されたように、ＦＩＲフィルタＧｍ（ｚ）のフィルタ長が、人間の聴覚系の周波数特性を考慮に入れるように設定されるのに対し、等しく幅広いサブバンドが使用される場合、フィルタバンクは計算上効率的にしか実現できないので、個々のサブバンドの帯域幅は等しい。それぞれの遅延Dｍは、ＦＩＲフィルタＧｍ（ｚ）の必要とされるフィルタタップを削減するために各サブバンドに適用される。一般的に、フィルタ長及び遅延は、任意の設定することができ、上述される適応タップ割り当て方法及び適応遅延割り当て方法に従って調整されてよい。遅延を除き、図７ａのフィルタ構造は、図３に関してすでに説明された基本構造に相当する。
【００５１】
フィルタバンク２２及び２２‘は、さらに詳細に以下に説明される。図６に関してすでに説明されたように、偶数で積み重ねられたＧＤＦＴフィルタバンクの適用は、人間の耳の周波数分解能に一致するサブバンドフィルタ長を使用するときに特に望ましい。すなわち、最低サブバンドに適用されるＦＩＲフィルタは、最高の周波数分解能、したがって最高の数のフィルタ係数を有さなければならない。偶数で積み重ねられたＧＤＦＴフィルタバンクを使用するとき、弧の最低のサブバンドは実数値であり、したがってメモリ要件は４分の１削減され、必要される計算時間は、奇数で積み重ねられたＧＤＦＴフィルタバンクが使用される代替策に比較して半減される。２つ以上の音声チャネルには、通常、サブバンドで実装されるＦＩＲフィルタが提供され、したがってフィルタバンクのこの選択はかなりの改善を伴う。
【００５２】
メモリ要件は、サブバンド数と必要とされるメモリとの間の依存性を検討するときにさらに削減されてよい。一方、サブバンドの増加する数は、個々のサブバンドの帯域幅の減少を伴う。その結果、プロトタイプフィルタは、必ずさらに多くのフィルタ係数を生じさせるより険しい通過帯域／ストップバンド遷移を有さなければならず、したがって特に遅延回線を実装するために必要とされるメモリ量が増加する。他方、多数のサブバンドは、（ＦＦＴアルゴリズムで使用される）より大きなブロックサイズを伴い、その結果として、内部メモリを外部メモリにスワップするためにより多くの計算時間を使用することができ、このようにして必要とされる内部メモリを削減する。サブバンド数を決定するときには妥協策を見つけ出さなければならないことが上記から明らかとなるはずである。分析は、音声用途の場合、Ｍ＝１６、Ｍ＝３２、及びＭ＝６４のサブバンドのために良好な結果を得ることができ、それによってＮ＝Ｍ／２のアンダーサンプリングが適用される。すなわち、Ｍ＝３２サブバンドを使用するとき、サブバンド信号のサンプリングレートは、Ｎ＝１６という係数で削減されてよく、ＦＦＴブロックサイズも１６に設定される。後者の例（Ｍ＝３２）では、１７のサブバンド（Ｍ／２＋１）だけが実際に処理されなければならず、それによって（最高のサブバンドだけではなく）最低のサブバンドは上述されたように実数値である。
【００５３】
大量の計算時間を消費するだろう長い初期化を生じさせずに、可能な限り多くの内部メモリをより安価な外部メモリにスワップすることも望ましい。サブバンド遅延Dｍを実装するために必要とされるメモリは、外部メモリにスワップされるには適切であると思われる。しかし、これは多数のメモリ転送初期化を生じさせるだろう。例えば、１０の音声チャネルの内の１７のサブバンドがブロックごとに（読み取り及び書き込みのために）１７回のメモリ転送初期化を必要とする。これは初期化毎に３ミップスを仮定するとき、１０２０ミップスの計算力が必要とされることを意味する。この問題は、図７ｂに示される完全帯域（広帯域）ドメインの中に最も短いサブバンド遅延を「移動させる」ことによって、図７ａに示される構造を再配置するときに緩和できる。図７ｂの例では、最も短いサブバンド遅延DＭＩＮが第１の遅延D１である。広帯域ドメインの対応する遅延は、「バルク遅延」として示され、DＢＵＬＫ＝Ｎ・DＭＩＮに等しい。残りのサブバンド遅延は、図７ｂに示されるようにDＭＩＮ分削減される。
【００５４】
上述されたように、１７回の初期化の代わりに（読み取りアクセス用及び書き込みアクセス用に）１つの転送初期化だけしか必要でないので、広帯域バルク遅延は、外部メモリにスワップされてよい。残りの―ただし、今はさらに短い―サブバンド遅延が内部メモリに残る。ただし、図７ｂに示されている解決策はさらに改善され得る。最低サブバンドのＦＩＲフィルタが最長インパルス応答（つまり、フィルタタップの最大数）、したがって残りの遅延よりもはるかに短い最短遅延を有するので、他のサブバンド遅延のために必要とされる（内部）メモリは依然として多い。この状況は、サブバンドフィルタ構造から最低のサブバンドを「抽出し」、図８に示されるように別々にそれを処理することによって、図７ｂの構造を修正するときにさらに改善されてよい。
【００５５】
図７ｂの例と対照的に、バルク遅延DＢＵＬＫは、（通常は最低のサブバンドの遅延D１である）最短のサブバンド遅延に基づいて形成されるのではなく、むしろ例えばD２のような第２の最短の遅延に基づいて形成される。遅延D３−D２、D４−D２、 D５−D２等の間の差は図７ｂの例においてよりはるかに小さいので、必要とされるメモリの達成可能な削減はより大きくなる。したがって、広帯域入力信号ｘ［ｎ］は、フィルタバンクプロトタイプフィルタ２２‘を称してフィルタリングされ、次いで係数Ｎ（デシメータ２５）でダウンサンプリングされ 、係数２Ｍ（利得２６）で乗算され、最後に元のサブバンド遅延D１で遅延され、第１のサブバンドの出力信号ｙ１［ｎ］を得るために対応する（実数値の）サブバンドＦＩＲフィルタでフィルタリングされる。残りのサブバンド出力信号ｙ２［ｎ］， ｙ３［ｎ］， ．．．， ｙＭ／２＋１［ｎ］は図７ｂの例においてと同様に取得される。唯一の相違点は、図７ｂの例においてよりも大きく、相応してより短いサブバンド遅延D３−D２， D４−D２， D５−D２を生じさせるバルク遅延DＢＵＬＫ，のサイズである。全てのサブバンド出力信号ｙｍ［ｎ］は、前の例と同じように広帯域出力信号ｙ［ｎ］に結合される。
【００５６】
ＦＩＲフィルタ２２‘（つまり、ＡＦＢプロトタイプフィルタＰ（ｚ））を実装するための計算作業、及び以後のダウンサンプリングを削減するために、ＦＩＲフィルタ２２’及びデシメータ２５は、計算上より効率的である対応する多相フィルタによって置き換えられてよい（図９を参照すること）。図９の図では、フィルタ転送関数Ｐ０（ｚ）、Ｐ１（ｚ）等はフィルタ係数 ［ｐ０， ｐＮ， ｐ２Ｎ， ．．．］、［ｐ１， ｐＮ＋１， ｐ２Ｎ＋１， ．．．］等を含むフィルタを表し、それによって係数ｐ０、ｐ１、 ｐ２，等がＡＦＢプロトタイプフィルタＰ（ｚ）のフィルタ係数を表している。
【００５７】
図８又は図９の例を使用すると、本例（１７のサブバンド、１０のチャネル）では３４０回の初期化から２０回の初期化へ、メモリ転送初期化の大幅な削減が可能になる。これらは、全てのチャネルのバルク遅延をスワップし、一度にスワップし直すことができるように、バルク遅延を等化することによってさらに２に削減されてよい。この手順は図１０に示される。
【００５８】
本発明及びその優位点は詳細に説明されてきたが、添付特許請求の範囲によって定められる本発明の精神及び範囲から逸脱することなく、多様な変更、置換、及び改変を本明細書に加えることができることが理解されるべきである。
【００５９】
さらに、本発明の範囲は、本明細書に説明されるプロセス、機械、製造、物質の組成、手段、方法、及びステップの特定の実施形態に制限されることが意図されていない。統合者は、本発明の開示から、本明細書に説明される対応する実施形態と実質的に同じ機能を実行する、又は実質的に同じ結果を達成する、現在存在している、又は後に開発されるプロセス、機械、製造、物質の組成、手段、方法、及びステップが本発明に従って活用され得ることを容易に理解する。したがって、添付特許請求の範囲は、それらの範囲内で、かかるプロセス、機械、製造、物質の組成、手段、方法、及びステップを含むことが意図される。
【特許請求の範囲】
【請求項１】
少なくとも１つのデジタルフィルタが各チャネルに実装されるデジタル音声プロセッサを備える少なくとも１つの音声出力チャネルのある音声システムであって、各チャネルの前記デジタルフィルタが、
広帯域入力音声信号を受信し、前記受信した音声信号を複数のサブバンドに分割し、このようにして等しい帯域幅を有するサブバンド信号を提供するように構成される解析フィルタバンクであって、前記サブバンド信号の前記スペクトルが前記入力音声信号の前記スペクトルを構成する解析フィルタバンクと、
各サブバンド用のサブバンドＦＩＲフィルタであって、このようにしてそれぞれのフィルタリングされたサブバンド信号を提供するサブバンドＦＩＲフィルタと、
前記フィルタリングされたサブバンド信号を受信し、広帯域出力音声信号を提供するためにそれらを結合するように構成される合成フィルタバンクと、
を備え、
遅延が各サブバンド信号に関連付けられ、前記サブバンド信号の内の１つの前記遅延が前記解析フィルタバンクの上流の前記広帯域入力音声信号に適用され、前記残留遅延が前記解析フィルタバンクの下流の前記残りのサブバンド信号に適用される、
音声システム。
【請求項２】
前記遅延されていない入力音声信号を受信し、追加サブバンド信号として、前記解析フィルタバンクでタッピングされていない最低周波数サブバンドの前記サブバンド信号を提供する低域通過フィルタと、
前記それぞれのサブバンドＦＩＲフィルタの上流又は加硫で前記最低周波数サブバンド信号に適用される遅延と、
をさらに備え、
前記広帯域入力音声信号に適用される前記遅延が、前記解析フィルタバンクでタッピングされる前記サブバンドと関連付けられる最短の遅延に一致する、
請求項１に記載のシステム。
【請求項３】
前記サブバンドＦＩＲフィルタのフィルタ係数の前記数が、銭サブバンドフィルタの前記周波数分解能が、少なくともほぼ前記人間の耳の前記周波数分解能に相当するように設定される、
請求項１又は２に記載のシステム。
【請求項４】
少なくとも第１の音声チャネルと、第２の音声チャネルとを備え、
各広帯域入力音声信号が、前記解析フィルタバンクでタッピングされる全てのチャネルの全てのサブバンド信号の最短の遅延によって与えられる等しい遅延で遅延される、
請求項１から３の１つに記載のシステム。
【請求項５】
前記解析フィルタバンク（複数の場合がある）及び前記合成フィルタバンク（複数の場合がある）が、前記入力音声信号をブロック単位で処理する偶数で積み重ねられるＧＤＴＦフィルタであり、前記サブバンド信号が、前記ブロックサイズに対応する係数で、前記広帯域信号に比較してダウンサンプリングされ、
前記処理されたサブバンド信号の前記最低周波数及び前記最高周波数のサブバンド信号が実数値であり、前記他のサブバンド信号が複素数値である、
請求項１から４の１つに記載のシステム。
【請求項６】
前記低域通過フィルタが、前記解析フィルタバンク（複数の場合がある）によって使用されるプロトタイプフィルタに相当する、請求項２に記載のシステム。
【請求項７】
音声システムの少なくとも１つの音声出力チャネル内に提供される少なくとも１つの音声信号をフィルタリングするための方法であって、各チャネルのために、
所与の周波数バンド内にあるスペクトルを有する、フィルタリングされる広帯域入力音声信号を受信することと、
前記入力音声信号の前記周波数バンドを、複数のサブバンドに分割し、このようにしてサブバンド信号を提供することであって、前記サブバンド信号の前記スペクトルが前記入力音声信号の前記スペクトルを構成する、提供することと、
各サブバンド信号をフィルタリングし、このようにしてサブバンドごとにそれぞれのＦＩＲ又はＩＩＲフィルタを使用してそれぞれのフィルタリングされたサブバンド信号を提供することと、
広帯域出力音声信号を提供するために、前記フィルタリングされたサブバンド信号を結合することと、
を含み、
遅延が各サブバンドに関連付けられ、前記サブバンドの内の１つの前記遅延が、前記周波数バンドを分割する前に前記広帯域入力音声信号に適用され、前記残留遅延が前記周波数バンドを分割後に前記残りのサブバンド信号に適用される、
方法。
【請求項８】
前記周波数バンドを分割することが、前記入力信号をブロック単位で処理する第１の偶数で積み重ねられたＧＤＴＦ解析フィルタバンクを使用して実行され、
前記周波数バンドを分割することが、前記広帯域信号と比較して、前記ブロックサイズに相当する係数分、前記サブバンド信号を前記ダウンサンプリングすることをさらに含み、
前記フィルタリングされたサブバンド信号を結合することが、前記第１のフィルタバンクに対応し、前記広帯域出力信号をブロック単位で提供する第２の偶数で積み重ねられたＧＤＴＦ合成フィルタバンクを使用して実行される、
請求項７に記載の方法。
【請求項９】
前記遅延していない入力音声信号を低域通過フィルタリングし、追加サブバンド信号として、前記最低周波数サブバンドの前記サブバンド信号を提供することと、
それぞれのＦＩＲ又はＩＩＲフィルタを使用して前記最低周波数サブバンド信号をフィルタリングし、このようにして対応するフィルタリングされたサブバンド信号を提供することと、
ＦＩＲフィルタリングの前又は後に前記最低周波数サブバンド信号を表示することと、
をさらに含む、請求項７又は８に記載の方法。
【請求項１０】
前記サブバンドＦＩＲフィルタのフィルタ係数の前記数が、前記サブバンドフィルタの前記周波数分解能が前記人間の耳の前記周波数分解能に少なくともほぼ一致するように設定される、請求項７から９の１つに記載の方法。
【請求項１１】
前記音声システムが、第１の音声出力チャネル及び第２の音声出力チャネルを備え、
各広帯域入力音声信号が、前記周波数バンドを除算することによって得られる全てのチャネルの全てのサブバンド信号の最短遅延によって与えられる等しい遅延で遅延される、
請求項７から１０の１つに記載の方法。
【請求項１２】
前記低域通過フィルタリングが、低域通過フィルタとして、前記周波数バンドを分割するために使用される前記ＧＤＦＴフィルタバンク（複数の場合がある）によって使用されるプロトタイプフィルタに相当するフィルタを使用することを含む、
請求項９に記載の方法。
【請求項１】
少なくとも１つのデジタルフィルタが各チャネルに実装されるデジタル音声プロセッサを備える少なくとも１つの音声出力チャネルのある音声システムであって、各チャネルの前記デジタルフィルタが、
広帯域入力音声信号を受信し、前記受信した音声信号を複数のサブバンドに分割し、このようにして等しい帯域幅を有するサブバンド信号を提供するように構成される解析フィルタバンクであって、前記サブバンド信号の前記スペクトルが前記入力音声信号の前記スペクトルを構成する解析フィルタバンクと、
各サブバンド用のサブバンドＦＩＲフィルタであって、このようにしてそれぞれのフィルタリングされたサブバンド信号を提供するサブバンドＦＩＲフィルタと、
前記フィルタリングされたサブバンド信号を受信し、広帯域出力音声信号を提供するためにそれらを結合するように構成される合成フィルタバンクと、
を備え、
遅延が各サブバンド信号に関連付けられ、前記サブバンド信号の内の１つの前記遅延が前記解析フィルタバンクの上流の前記広帯域入力音声信号に適用され、前記残留遅延が前記解析フィルタバンクの下流の前記残りのサブバンド信号に適用される、
音声システム。
【請求項２】
前記遅延されていない入力音声信号を受信し、追加サブバンド信号として、前記解析フィルタバンクでタッピングされていない最低周波数サブバンドの前記サブバンド信号を提供する低域通過フィルタと、
前記それぞれのサブバンドＦＩＲフィルタの上流又は加硫で前記最低周波数サブバンド信号に適用される遅延と、
をさらに備え、
前記広帯域入力音声信号に適用される前記遅延が、前記解析フィルタバンクでタッピングされる前記サブバンドと関連付けられる最短の遅延に一致する、
請求項１に記載のシステム。
【請求項３】
前記サブバンドＦＩＲフィルタのフィルタ係数の前記数が、銭サブバンドフィルタの前記周波数分解能が、少なくともほぼ前記人間の耳の前記周波数分解能に相当するように設定される、
請求項１又は２に記載のシステム。
【請求項４】
少なくとも第１の音声チャネルと、第２の音声チャネルとを備え、
各広帯域入力音声信号が、前記解析フィルタバンクでタッピングされる全てのチャネルの全てのサブバンド信号の最短の遅延によって与えられる等しい遅延で遅延される、
請求項１から３の１つに記載のシステム。
【請求項５】
前記解析フィルタバンク（複数の場合がある）及び前記合成フィルタバンク（複数の場合がある）が、前記入力音声信号をブロック単位で処理する偶数で積み重ねられるＧＤＴＦフィルタであり、前記サブバンド信号が、前記ブロックサイズに対応する係数で、前記広帯域信号に比較してダウンサンプリングされ、
前記処理されたサブバンド信号の前記最低周波数及び前記最高周波数のサブバンド信号が実数値であり、前記他のサブバンド信号が複素数値である、
請求項１から４の１つに記載のシステム。
【請求項６】
前記低域通過フィルタが、前記解析フィルタバンク（複数の場合がある）によって使用されるプロトタイプフィルタに相当する、請求項２に記載のシステム。
【請求項７】
音声システムの少なくとも１つの音声出力チャネル内に提供される少なくとも１つの音声信号をフィルタリングするための方法であって、各チャネルのために、
所与の周波数バンド内にあるスペクトルを有する、フィルタリングされる広帯域入力音声信号を受信することと、
前記入力音声信号の前記周波数バンドを、複数のサブバンドに分割し、このようにしてサブバンド信号を提供することであって、前記サブバンド信号の前記スペクトルが前記入力音声信号の前記スペクトルを構成する、提供することと、
各サブバンド信号をフィルタリングし、このようにしてサブバンドごとにそれぞれのＦＩＲ又はＩＩＲフィルタを使用してそれぞれのフィルタリングされたサブバンド信号を提供することと、
広帯域出力音声信号を提供するために、前記フィルタリングされたサブバンド信号を結合することと、
を含み、
遅延が各サブバンドに関連付けられ、前記サブバンドの内の１つの前記遅延が、前記周波数バンドを分割する前に前記広帯域入力音声信号に適用され、前記残留遅延が前記周波数バンドを分割後に前記残りのサブバンド信号に適用される、
方法。
【請求項８】
前記周波数バンドを分割することが、前記入力信号をブロック単位で処理する第１の偶数で積み重ねられたＧＤＴＦ解析フィルタバンクを使用して実行され、
前記周波数バンドを分割することが、前記広帯域信号と比較して、前記ブロックサイズに相当する係数分、前記サブバンド信号を前記ダウンサンプリングすることをさらに含み、
前記フィルタリングされたサブバンド信号を結合することが、前記第１のフィルタバンクに対応し、前記広帯域出力信号をブロック単位で提供する第２の偶数で積み重ねられたＧＤＴＦ合成フィルタバンクを使用して実行される、
請求項７に記載の方法。
【請求項９】
前記遅延していない入力音声信号を低域通過フィルタリングし、追加サブバンド信号として、前記最低周波数サブバンドの前記サブバンド信号を提供することと、
それぞれのＦＩＲ又はＩＩＲフィルタを使用して前記最低周波数サブバンド信号をフィルタリングし、このようにして対応するフィルタリングされたサブバンド信号を提供することと、
ＦＩＲフィルタリングの前又は後に前記最低周波数サブバンド信号を表示することと、
をさらに含む、請求項７又は８に記載の方法。
【請求項１０】
前記サブバンドＦＩＲフィルタのフィルタ係数の前記数が、前記サブバンドフィルタの前記周波数分解能が前記人間の耳の前記周波数分解能に少なくともほぼ一致するように設定される、請求項７から９の１つに記載の方法。
【請求項１１】
前記音声システムが、第１の音声出力チャネル及び第２の音声出力チャネルを備え、
各広帯域入力音声信号が、前記周波数バンドを除算することによって得られる全てのチャネルの全てのサブバンド信号の最短遅延によって与えられる等しい遅延で遅延される、
請求項７から１０の１つに記載の方法。
【請求項１２】
前記低域通過フィルタリングが、低域通過フィルタとして、前記周波数バンドを分割するために使用される前記ＧＤＦＴフィルタバンク（複数の場合がある）によって使用されるプロトタイプフィルタに相当するフィルタを使用することを含む、
請求項９に記載の方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７ａ】
【図７ｂ】
【図８】
【図９】
【図１０】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７ａ】
【図７ｂ】
【図８】
【図９】
【図１０】
【公開番号】特開２０１２−２４９２９４（Ｐ２０１２−２４９２９４Ａ）
【公開日】平成２４年１２月１３日（２０１２．１２．１３）
【国際特許分類】
【外国語出願】
【出願番号】特願２０１２−１２１８９５（Ｐ２０１２−１２１８９５）
【出願日】平成２４年５月２９日（２０１２．５．２９）
【出願人】（５０４１４７９３３）ハーマン ベッカー オートモーティブ システムズ ゲーエムベーハー (165)
【Ｆターム（参考）】
【公開日】平成２４年１２月１３日（２０１２．１２．１３）
【国際特許分類】
【出願番号】特願２０１２−１２１８９５（Ｐ２０１２−１２１８９５）
【出願日】平成２４年５月２９日（２０１２．５．２９）
【出願人】（５０４１４７９３３）ハーマン ベッカー オートモーティブ システムズ ゲーエムベーハー (165)
【Ｆターム（参考）】
[ Back to top ]