JP4402040B2 - リスニング領域内での、バスサウンドの音質を改善するためのイコライゼーションシステム - Google Patents

Description

本発明は一般に、リスニング領域内での、一つ以上のラウドスピーカによって生成されたバスサウンドの音質を改善することを対象とする。より詳細には、本発明は、その領域内における応答が、所望される周波数範囲内で、実質的に一定およびフラットであるように、リスニング領域内で、少なくとも一つのラウドスピーカによって生成される応答をイコライズすることを対象とする。

サウンドシステムは通常、電気信号を音響信号に変換するラウドスピーカを含む。ラウドスピーカは、高周波数信号、中周波数信号、低周波数信号などの音響信号の範囲を生成する一つ以上の変換器を含み得る。ラウドスピーカの一つのタイプは、２０Ｈｚから１００Ｈｚの範囲において低周波数信号を生成する低周波数変換器を含み得るサブウーハーである。

サウンドシステムは、様々なリスニング環境において音響信号を生成し得る。リスニング環境の例は、住居のリスニングルーム、ホームシアター、映画劇場、コンサートホール、車内、レコーディングスタジオなどを含むが、それらに限定されない。通常、リスニング環境は、ラウドスピーカによって生成される音響信号を聴く、個人または複数の人のための、単一または多数のリスニング位置を含む。リスニング位置は、ホームシアター環境におけるソファなどのような、座った位置、または、コンサートホールにおいて指揮者が立ち得る場所などのような、立ち位置であり得る。

リスニング環境は、リスニング位置での、低周波数信号、中周波数信号、および／または高周波数信号を含む音響信号に影響を与え得る。部屋の種類、ならびに、部屋におけるリスナーの位置および部屋におけるラウドスピーカの位置に依存し、音の大きさは異なる周波数によって変化し得る。これは特に、低周波数に当てはまり得る。低周波数は、音楽、映画、およびオーディオ環境の他の多くの形式を楽しむことにとって重要であり得る。ホームシアターの例において、壁、カーテン（ｄｒａｐｅｒｙ）、家具、備え付け家具などを含む部屋の境界は、ラウドスピーカからリスニング位置へ移動する音響信号に影響を与え得る。

リスニング位置で受信された音響信号は計測され得る。部屋を特徴付ける一つの方法は、リスニング領域にあるマイクロフォンへのラウドスピーカのインパルス応答である。インパルス応答は、ラウドスピーカから放たれた短いサウンドバーストのために、マイクロフォンによって計測される音響信号である。インパルス応答は、単一の周波数、多数の別個の周波数、または、周波数の一つの範囲での、振幅および／または位相を含む、音響信号の様々な性質を計測させ得る。

振幅応答は、当該周波数での音の大きさの計測である。一般に、音の大きさまたは振幅は、デシベル（ｄＢ）で計測される。振幅の偏差は、指定された目的値に関連する正のデシベル値または負のデシベル値として表現され得る。リスニング位置で計測された振幅値が目的値に近ければ近いほど、その振幅応答は良い。目的からの偏差は、部屋の境界と相互作用する場合に、音響信号において生じた変化を反映する。ピーク（ｐｅａｋ）は、目的からの正の振幅偏差を表し、一方で、窪み（ｄｉｐ）は、目的からの負の振幅偏差を表す。

振幅応答におけるこれらの偏差は、サブウーハー、サブウーハーの位置、およびリスナーの位置で再生された音響信号の周波数に依存し得る。リスナーは、サウンドトラックや映画のような記録媒体上に録音されたのと同じような低周波数を聴き得るのではなく、その代わり、部屋の境界によって歪められた低周波数を聞き得る。このように、部屋は、サブウーハーによって再生された音響信号を変化させ得、サウンドシステムの低周波数の実行に逆に影響を与え得る。一つの例として、図１が長方形の部屋におけるサウンドシステムの設定を示す。サウンドシステムは、四つのサブウーハーに接続された受信器を含み、それぞれの受信器が部屋の隅にある。その部屋は、四つのサブウーハーによって生成される低周波数音およびバスサウンドに影響を与え得る四つの壁によって境界を定められる。部屋内で、座席領域は、一人以上の個人に、四つのサブウーハーのそれぞれによって生成される組み合わされたバスサウンドを聴かせるように提供される。上で検討されたように多くの要因は、個人がほんの数フィート離れて座っている他の人よりも、より大きな音を聴き得るように、リスニングエリア内で音質に影響を与え得る。部屋のインパルス応答を計測する目的のために、受信器は、対数周波数掃引出力信号を、所定の時間の間、四つのサブウーハーに送信し得る。部屋のインパルス応答は、次いで、その部屋のリスニング領域内での異なる位置に設置されたマイクロフォンＰ１、Ｐ２、Ｐ３、およびＰ４によって受信される。

図２は、個々に四つのマイクロフォンの位置Ｐ１、Ｐ２、Ｐ３、およびＰ４において期待し得る計測されたインパルス応答に対応して、四つの周波数応答曲線Ｆ１、Ｆ２、Ｆ３、およびＦ４を示す。以前に検討したように、サブウーハーは通常、２０Ｈｚから１００Ｈｚの間の低周波数範囲にて動作する。図２は、およそ４８Ｈｚで、バスサウンドの大きさがリスニング領域内に個人が位置する場所に依存するように、バスサウンドの大きさは、広い範囲において変化する。例えば、曲線Ｆ２は、バスの大きさのレベルが、およそ４８Ｈｚでのおよそ０デシベルであることを示し、一方で、曲線Ｆ３は、バスの大きさのレベルが、同じ周波数位置での、およそ−１８デシベルであることを示す。これは、Ｐ２の位置に座っている人が、その人のちょうど後ろに座っているＰ３の位置の人よりも、４８Ｈｚでの、より大きなバスサウンドを聴くことを意味する。つまり、それぞれの人が異なるバスサウンドの音質を経験するように、サウンドレベルは部屋のリスニング領域内での異なる位置では同じではない。付け加えて、図２は、当該周波数範囲内で変動することを示す。これは、リスニングされることを意図されているにもかかわらず、個人がそのバスサウンドを聴くことができないように、あるバスサウンドが欠落することを意味する。例えば、曲線Ｆ４は、およそ４８Ｈｚから５５Ｈｚの間で、およそ５２Ｈｚのバスの大きさのレベルにおいて、相当な欠落があることを示す。これは、Ｐ３の位置に座っている人が、４８Ｈｚでのバスサウンドを聴くが、５２Ｈｚでのバスサウンドにおいて唐突な欠落、および、再び５５Ｈｚでの唐突なピークに気付くことを意味する。バスサウンドレベルにおけるそのような変動は、リスニング経験を害し得る。

多くのイコライゼーション技術が、過去に、リスニング領域内での振幅偏差を低減し、または取り除くように利用されてきた。その技術の一つは、多数のリスニング位置の平均振幅応答を計算し、次いで、システムにおける全てのサブウーハーのイコライゼーションを等しく実施する、空間的平均化である。しかしながら、空間的平均化は、現実には存在しない単一の「平均リスニング位置」を是正するのみである。このように、空間的平均化の技術を利用した場合でさえ、一部のリスニング位置は他の位置よりもより良い低周波数のパフォーマンスを有するが、他の位置では、深刻に影響を与え得る。例えば、空間的平均化は、不等なパフォーマンスと比較して、一部のリスニング位置ではパフォーマンスをより悪くし得る。さらに、単一の位置のために振幅応答をイコライズしフラットにすることを企てることは、潜在的に問題を生み出す。ピークが平均リスニング位置で低減され得る一方で、窪みが生じる周波数を増幅する企ては、サブウーハーから追加的な音声出力を相当に要求し、その結果、システムの最大音声出力を低減し、部屋の他の領域における大きなピークを潜在的に生み出してしまう。

他の既知のイコライゼーション技術は、「モードキャンセル」（ｍｏｄｅ ｃａｎｃｅｌｉｎｇ）配置において、多数のサブウーハーを位置させることである。リスニングルーム内で、対称に多数のラウドスピーカを位置させることによって、定在波は、相殺的および建設的干渉を引き出すことによって低減させ得る。しかしながら、対称の「モードキャンセル」構成は、理想的な部屋（例えば、寸法上および音響上、対称な）を想定しており、様々な形の部屋や備え付け家具を含む実際の部屋の特徴を考慮しない。さらに、ラウドスピーカを対称に位置させることは、特定の部屋のセッティングの、現実的な、または所望される構成ではあり得ない。

さらに別のイコライゼーション技術は、数学的解析を利用し、振幅偏差を低減するために、オーディオシステムを構成する。そのような一つの数学的解析は、その部屋のデータに基づいた部屋における定在波を促す。例えば、部屋の長さ、幅、高さのような部屋の寸法が入力され、様々なアルゴリズムが入力されたデータに基づき、サブウーハーを位置させる場所を予測する。しかしながら、この数学的方法は、部屋の家具、備え付け家具、構成などの音響的性質を考慮しない。例えば、石造外装を有する内壁は、木製骨組みの壁とは音響上、全く異なるように作用し得る。さらに、この数学的方法は、部分的に閉じられた部屋を効果的に埋め合わせることができず、部屋が長方形でない場合も、計算的に、難儀になり得る。

部屋におけるインパルス応答をイコライズしようとする多くの他の方法があるが、イコライゼーションの正確性は、信号に与えられた遅延およびゲインなどのような、あるパラメータを決定することを含む推測による、偶然であり得ることが多い。正確なイコライゼーションを解決するために、多大な計算出力がかかる。さらに、バスレベルの大きさがそれぞれの座席位置で一貫しているだけでなく、所望される低周波数範囲全体で、実質的に一定またはフラットでもあるように、所望される低周波数範囲内での、フラットな周波数応答になるイコライゼーションを、これらの方法は提供しない。それゆえ、所定の場所における一つ以上のリスニング位置の音声パフォーマンスが改善されるように、オーディオシステムの構成を正確に決定するシステムが、長年必要とされている。

本発明は、リスニングルームにおける低周波数イコライゼーションの、既知の問題に取り組む。本発明は、好ましいリスニング領域内での、部屋の様々な位置のインパルス応答を計測するための、一つ以上のマイクロフォンを利用する、周波数イコライゼーションシステムを対象にする。次いで、この情報は、周波数応答が、妥当な周波数範囲を介し、マイクロフォン計測ポイントにて、および、所望されるリスニング領域内にて、実質的にフラットであるように、バス応答を改善するために、その部屋におけるサブウーハーへ送信された音響信号をフィルタするために利用される。

本発明は、周波数応答が、妥当な周波数範囲を介し、リスニング領域内にて実質的にフラットであるように、それぞれの対応するサブウーハーのためのフィルタをデザインする係数を計算する、部屋のインパルス応答を使用する。一般に、部屋の応答の反転（ｉｎｖｅｒｓｅ）は、部屋によって追加される音色（ｃｏｌｏｒａｔｉｏｎ）を取り消すように決定される。その反転は、サブウーハーが扱い得る許容可能なゲインを超過し得る唐突なゲインが最小化または取り除かれるように、平滑化される。本発明はまた、サブウーハーが最適に動作する所望される周波数範囲に、イコライゼーションが与えられるように、反転上に、目的関数（ｔａｒｇｅｔ ｆｕｎｃｔｉｏｎ）を与え得る。修正された反転は、次いで、その個々のサブウーハーに送信されたそれぞれの音響信号のためのフィルタ係数を決定するために使用される。デジタル信号プロセッサ（ＤＳＰ）などのようなプロセッサは、フィルタ係数に基づいて、音響信号をフィルタするために利用され得る。

本発明は、さらに以下の手段を提供する。

（項目１）
リスニング領域における周波数範囲内で、周波数応答を実質的にイコライズするための、一つ以上のフィルタをデザインする方法であって、該方法は、
少なくとも一つの音響変換器によって生成されたリスニング領域における該周波数応答を計測すること（５０２）と、
第１のステージのイコライゼーションフィルタスペクトルを決定するために、該周波数応答を反転すること（５０４）と、
第２のステージの近似の（ａｐｐｒｏｘｉｍａｔｅ）イコライゼーションフィルタスペクトルを決定するために、該第１のステージのイコライゼーションフィルタスペクトルを平滑化すること（５０６）と、
該第２のステージの近似のイコライゼーションフィルタを、該計測された周波数応答に適合した後に、その結果である周波数応答の組み合わせから広範囲周波数応答を決定すること（５０８）と、
広範囲イコライゼーションフィルタを決定するために、該広範囲周波数応答を反転すること（５１０）と、
最終イコライゼーションフィルタを決定するために、該広範囲イコライゼーションフィルタと該第２のステージの近似のイコライゼーションフィルタとを結合すること（１６０２）と
を包含する、方法。

（項目２）
上記周波数応答を計測することは、
上記リスニング領域内に位置された少なくとも一つのマイクロフォンを介して、上記少なくとも一つの音響変換器によって生成された、上記リスニング領域のインパルス応答を受信すること（６０２）と、
該リスニング領域における該インパルス応答から任意の共通する時間遅延を取り除くこと（６０４）と、
該リスニング領域における該インパルス応答を、該リスニング領域における該周波数応答に変換すること（６０６）と
を包含する、項目１に記載の方法。

（項目３）
所望される周波数帯域内で、ゲインを制限するために、上記最終イコライゼーションフィルタの振幅応答をクリップすること（１６０４）を含む、項目１または２に記載の方法。

（項目４）
上記周波数応答を計測することは、
上記リスニング領域内に位置された少なくとも一つのマイクロフォンを介して、上記少なくとも一つの音響変換器によって生成された、上記リスニング領域のインパルス応答を受信すること（７０２）と、
上記インパルス応答を計測するために、上記リスニング領域における上記少なくとも一つの音響変換器の数と同じ数の、部屋における少なくとも一つのマイクロフォンを使用すること（７０６）と
を包含する、項目１から３のうちのいずれか一項に記載の方法。

（項目５）
複素平滑化方法を介して、上記周波数応答の反転がなされる（５０６）、項目１から４のうちのいずれか一項に記載の方法。

（項目６）
上記複素平滑化が、二つの異なる平滑化インデックス値を有する、二つの別個の周波数帯域にてなされる（５０６）、項目５に記載の方法。

（項目７）
上記少なくとも一つの上記音響変換器の数が四つであり、上記少なくとも一つのマイクロフォンの数が四つである（７０６）、項目４に記載の方法。

（項目８）
上記周波数応答を計測することは、
上記リスニング領域内に位置された少なくとも一つのマイクロフォンを介して、上記少なくとも一つの音響変換器によって生成された、上記リスニング領域のインパルス応答を受信すること（７０２）を包含し、
該インパルス応答を計測するために利用される少なくとも一つのマイクロフォンの数が、少なくとも一つの音響変換器の数と等しくない場合、該周波数応答の反転が、擬似反転法を介してなされる（７０８）、項目１から７のうちのいずれか一項に記載の方法。

（項目９）
上記リスニング領域における上記少なくとも一つの音響変換器のそれぞれが動作する周波数範囲を制限するために、目的関数を上記周波数応答に適合すること（１１０６）を包含する、項目１から８のうちのいずれか一項に記載の方法。

（項目１０）
上記少なくとも一つの音響変換器がサブウーハーである、項目９または１０に記載の方法。

（項目１１）
部屋内でのリスニング領域における周波数範囲内の周波数応答を実質的にイコライズするための、一つ以上のフィルタをデザインするシステムであって、
出力信号を生成するための、少なくとも一つの音響変換器（５０２）と、
該リスニング領域における周波数応答を計測するための、少なくとも一つのマイクロフォン（５０２）と、
該音響変換器および該マイクロフォンに接続されるプロセッサであって、該音響変換器に、該出力信号を生成させ、該マイクロフォンを介して部屋の音響応答を計測させるために、テスト信号を送信することができ、ならびに、該音響変換器への該フィルタされた出力信号が、部屋の該リスニング領域における周波数範囲内の、実質的にフラットな周波数応答を生成するように、該計測された音声応答に基づき、フィルタ係数を計算することができる、プロセッサと、
を備える、システム。

（項目１２）
上記部屋における上記少なくとも一つの音響変換器の数が、該部屋における上記少なくとも一つのマイクロフォンの数と等しい（７０６）、項目１１に記載のシステム。

（項目１３）
上記部屋における上記少なくとも一つの音響変換器の数が四つであり、該部屋における上記少なくとも一つのマイクロフォンが四つである（７０６）、項目１１または１２に記載のシステム。

（項目１４）
上記プロセッサがデジタル信号プロセッサである、項目１１から１３のうちのいずれか一項に記載のシステム。

（項目１５）
上記プロセッサが上記フィルタ係数を計算するために命令を実行することができ、該命令が、
リスニング領域における、少なくとも一つの音響変換器によって生成された周波数応答を計測するための命令（５０２）と、
第１のステージのイコライゼーションフィルタスペクトルを決定するために、上記周波数応答を反転するための命令（５０４）と、
第２のステージの近似のイコライゼーションフィルタスペクトルを決定するために、該第１のステージのイコライゼーションフィルタスペクトルを平滑化するための命令（５０６）と
該第２のステージの近似のイコライゼーションフィルタを該計測された周波数応答に適合した後に、その結果としての周波数応答の組み合わせから、広範囲周波数応答を決定するための命令（５０８）と、
広範囲イコライゼーションフィルタを決定するために、該広範囲周波数応答を反転するための命令（５１０）と、
最終イコライゼーションフィルタを決定するために、該広範囲イコライゼーションフィルタと、該第２のステージの近似のイコライゼーションフィルタとを結合するための命令（１６０２）と
を含む、項目１１から１４のうちのいずれか一項に記載のシステム。

（摘要）
周波数イコライゼーションシステムは、リスニング領域における周波数応答が、所望される周波数範囲内にて、実質的に一定でフラットであるように、リスニング領域内で、少なくとも一つのラウドスピーカによって生成された部屋の周波数応答を実質的にイコライズする。周波数イコライゼーションシステムは、部屋のインパルス応答を計測するために多数のマイクロフォンを使用し、および、妥当な周波数範囲を介して、フラットである改善されたバス応答を達成するために、一つ以上のサブウーハーの音響信号を処理するフィルタをデザインするインパルス応答を使用する。システムは、閉形式で、非反復の、数学的解であるアルゴリズムを用い、非常に短い計算時間を特徴とする。

本発明の他のシステム、方法、特徴、および利点は、以下の図面および詳細な記載の検討によって、当業者にとっては明らかであり、または、明らかになる。全ての追加的なシステム、方法、特徴、および利点は、この記載内に含まれ、本発明の範囲内であり、および、請求項によって保護されることが意図される。

本発明は、以下の図面および記載を参照し、より良く理解され得る。図における構成要素は、必然的に寸法を合わせてはおらず、その代わり、本発明の原理を示すことに強調が置かれる。さらに、図において、類似する参照番号は、異なる図を通して、同様の対応する部分を示す。

図３は、本発明に従うイコライゼーションシステム３００を図示するブロック図を示し、部屋の所望されるリスニング領域内での所定の低周波数範囲にわたってフラットである、部屋内の一つ以上のサブウーハーからの改善されたバス応答を達成するようにデザインされている。イコライゼーションシステム３００は、それぞれの部屋がそれ自体独自の特徴を有する様々な部屋の周波数応答をイコライズするために利用され得る。例えば、一つの部屋が、以下の特徴の一つ以上を有し得る。（１）部屋の一つ以上の壁が開かれ得る（２）部屋の天上または壁が弧状を有し得る（３）カーテン（ｄｒａｐｅ）が一つ以上の部屋の壁を覆い得る（４）部屋の床に段差があり得る（５）部屋に一つ以上のサブウーハーがあり得る（６）サブウーハーのそれぞれの位置が部屋の任意の場所に位置され得る、など。以下で詳述されるように、イコライゼーションシステム３００は、任意の部屋の周波数応答をイコライズするために利用され得る。

これを検討するために、イコライゼーションシステム３００（ＥＱシステム３００）は、図１において示される部屋の応答をイコライズするために利用される。部屋は一般的に、長方形構成を形成する四つの壁によって遮られる。部屋内において、位置Ｐ１、Ｐ２、Ｐ３、およびＰ４によって定義される、一人以上の個人を座らせる座席領域がある。これらの座席領域は一般に、部屋のリスニング領域を定義する。受信器３０８は部屋内に位置され、サブウーハーに音響信号を送り、イコライゼーションシステム３００を組み込む。

ＥＱシステム３００は、テスト信号を生成し、部屋におけるラウドスピーカに対応するそれぞれのフィルタの係数をデザインすることができる信号ブロック３０２を含む。この例において、信号ブロック３０２は、部屋のそれぞれの隅に位置された四つのサブウーハーＳｕｂ１、Ｓｕｂ２、Ｓｕｂ３、およびＳｕｂ４とリンクされる。信号ブロック３０２は、部屋に置かれたＰ１からＰ４のマイクロフォンそれぞれへのサブウーハーのインパルス応答を計測するために、出力信号を、一度に四つのサブウーハーのそれぞれに送信し得る。信号ブロック３０２は、それぞれのサブウーハーへ、順次的に、所定の時間量の、対数周波数掃引を出力し得る。対数周波数掃引は、信号ブロック３０２に、サブウーハーを介して、当該の広い周波数スペクトルをカバーする出力信号を送信し得る。一例として、出力信号は、およそ４秒間、送信され得る。

時間が経過するとともに、サブウーハーのそれぞれが出力信号を送信し、インパルス応答は、図１におけるＰ１からＰ４によって表される位置などのような部屋の異なる領域（「リスニング位置」）に位置されたマイクロフォンによって、独立的または同時に、計測され得る。例えば、信号ブロック３０２は、マイクロフォンがリスニング領域の左上の隅において生成された信号から、部屋のインパルス信号を計測し得るように、Ｓｕｂ１を介して出力信号を送信し得る。次いで、信号ブロック３０２は、マイクロフォンがリスニング領域の右上の隅において生成された出力信号ソースによる、部屋のインパルス信号を計測し得るように、Ｓｕｂ２を介して他の出力信号を送信し得る。同様に、出力信号は、マイクロフォンが、個々に、リスニング領域の右下および左下の隅から送信された、続く個別の信号による、インパルス応答を計測し得るように、Ｓｕｂ３、または他のＳｕｂ４を介して送信され得る。この例において、長方形の部屋の四つの隅に置かれたサブウーハー、および、所望されるリスニング領域内に置かれた四つのマイクロフォンが、部屋のインパルス応答を計測するために使用される。マイクロフォンＰ１からＰ４は、音響信号を電気信号に変換する。電気信号が信号ブロック３０２に提供される前に、電気信号が、Ａ／Ｄ変換器を使用して、所定のレートでデジタル化され得る。

マイクロフォンを介して、信号ブロック３０２は、サブウーハーおよびマイクロフォンのそれぞれの組み合わせのために、秒毎の所定のインパルス応答サンプル数を取り込み得る。取り込まれたインパルス応答は、計測されたそれぞれのインパルス応答のＮサンプルを生み出すためにダウンサンプル（ｄｏｗｎ−ｓａｍｐｌｅ）され得る。四つのサブウーハーおよび四つのマイクロフォンを有し、これは、それぞれのセットがＮ数のサンプルを有する１６通りのインパルス応答という結果になる。例えば、信号ブロック３０２は、秒毎７５０サンプルのサンプルレートにて、Ｎ＝２０４８のサンプルを取り込み得る。

信号ブロック３０２は、マイクロフォンＰ１〜Ｐ４から、計測されたインパルス応答を受信する。信号ブロック３０２は、以下で述べるように、部屋のインパルス応答に基づいて、フィルタ係数を計算する。信号ブロック３０２は、プロセッサブロック３０４にリンクされ、プロセッサブロック３０４は、四つのサブウーハーによって生成されるサウンドによる室内周波数応答を実質的にイコライズするために、対応するサブウーハーへ送信される音響信号のそれぞれを修正するように、本発明に従って計算されるように意図されたフィルタを実施する。この例において、プロセッサブロック３０４は、図３において示されるように、ＦＩＲ１、ＦＩＲ２、ＦＩＲ３、およびＦＩＲ４によって表される四つの音響信号をフィルタし得、個々にサブウーハーＳｕｂ１、Ｓｕｂ２、Ｓｕｂ３、およびＳｕｂ４のそれぞれに対応する。ＴＶ、ＤＶＤプレーヤ、音声受信器などのような様々なソース３０８によって提供される音響信号入力３０６は、その個々のサブウーハーへ送信される出力信号が、部屋の周波数応答をイコライズするためのフィルタ係数に従ってフィルタされるように、対応するフィルタＦＩＲ１からＦＩＲ４を介して、プロセッサブロック３０４によって処理される。プロセッサブロック３０４は、デジタル信号プロセッサ（ＤＳＰ）などのような様々なプロセッサであり得、フィルタはＦｉｎｉｔｅ Ｉｍｐｕｌｓｅ Ｒｅｓｐｏｎｓｅ（ＦＩＲ）フィルタであり得る。一つのプロセッサに、信号ブロック３０２およびプロセッサブロック３０４によってなされる機能を実行させることは、本発明の範囲内であることに注意したい。

図４は、サブウーハーへの出力信号が本発明に従った応答をイコライズするようにフィルタした後、図１において示された部屋の周波数応答を示す。図４は、その結果としての振幅応答がそれぞれ関連する低周波数範囲において実質的に一貫していることを示す。これは、リスニング領域内のそれぞれの個人がそれぞれの周波数のポイントで、実質的に同様の大きさのレベルを提供されていることを意味する。付け加えて、振幅レベルは、サウンドレベルの欠落が実質的に最小であるように、４０Ｈｚから１００Ｈｚの間の所望される低周波数レベルにわたって、実質的に一定またはフラットである。図４と図２とを比較し、部屋の振幅応答は実質的に改善されている。以下で、本発明に従ったサブウーハーのそれぞれのために、フィルタがどのようにデザインされるのかを詳細に検討する。

以下の検討は、図１において示されるように、部屋内で、四つのサブウーハーおよび四つのマイクロフォン、例えばｎ_ｓｕｂ＝４、および、ｎ_ｍｉｃ＝４、の特定の場合である。しかしながら、本発明は、部屋におけるサブウーハーおよびマイクロフォンの任意の組み合わせのために使用され得る。一つ以上のサブウーハーへ送信される音響信号は、以下の記述に従ってフィルタされ得る。

図５は、部屋の周波数応答をイコライズするための、フィルタデザインの手順の概観を有するフローチャートである。ブロック５０２において、入力信号は、部屋の周波数応答を実質的にイコライズするように用意され得る。一般に、データを用意することは、部屋のインパルス応答を計測し、それらを周波数領域に変換することを含む。ブロック５０４において、周波数応答のそれぞれの反転（ｉｎｖｅｒｓｅ）が決定され得る。反転のそれぞれは、部屋の壁によって追加された音色（ｃｏｌｏｒａｔｉｏｎ）を事実上取り消す。つまり、その個々の反転を有する音響信号をフィルタすること、および、それら個々のサブウーハーへフィルタされた信号を送信することは、理想的な周波数応答を生成する。しかしながら、反転は、局地的で唐突なピークおよび窪みを有し得、そのような唐突なゲインは、サブウーハーが扱い得る許容可能なゲインを上回り得る。ブロック５０６において、反転における局地的なピークおよび窪みは、以下でさらに記載される複素平滑化（ｃｏｍｐｌｅｘ ｓｍｏｏｔｈｉｎｇ）法を用いて、平滑化され得る。これは、部屋の周波数応答の近似の反転を提供する。

ブロック５０８において、低周波数帯域端または高周波数帯域端での変遷を記載する目的関数が接近され得るように、近似の反転フィルタの後の結果に、広範囲イコライゼーションが適合される。広範囲イコライゼーションはまた、以下に記載されるように、ピークおよび窪みをアドレスする平滑化方法を使用する。サブウーハーは一般に１００Ｈｚ以下で動作するが、ブロック５１０において、２０Ｈｚ以下および／または１００Ｈｚ以上などのように、サブウーハーを保護する所望される低周波数範囲の外でサブウーハーに適合され得るゲインに、制限が置かれ得る。ブロック５１２において、広範囲イコライゼーションの反転は、次いで、部屋の周波数応答を実質的にイコライズするために、それぞれのサブウーハーに送信される音響信号のそれぞれを処理するためのフィルタを決定するように使用される。

図６は、図５におけるブロック５０２において表された部屋の入力データを用意する、さらなる詳細を示したフローチャート６００である。入力データを用意することは、上で検討されたように、部屋のインパルス応答を計測するブロック６０２を含む。ブロック６０２において、一度、インパルス応答が計測されると、ブロック６０４において、インパルス応答からの任意の共通の時間遅延が取り除かれ得る。これは、以下で検討される複素平滑化の数学的問題の解答を与えるためになされる。例えば、図１において示されるように、部屋の左上の隅に配置されたＳｕｂ１によって送信された出力信号に関して、マイクロフォンＰ１は、Ｓｕｂ１に最も近い。それゆえ、マイクロフォンＰ１は、他のマイクロフォンの前に、Ｓｕｂ１からの出力信号を受信する。Ｓｕｂ１からの出力信号がマイクロフォンＰ１に到達するのにかかる時間は、他のマイクロフォンＰ２からＰ４と共通している。この時間は、Ｓｕｂ１によって送信される出力信号が四つのマイクロフォンＰ１からＰ４によって計測されるインパルス応答に関して、共通の時間遅延として定義され得る。同様に、対応する共通の時間遅延は、他のＳｕｂ２からＳｕｂ４のそれぞれによって送信される出力信号のために計測され得る。例えば、Ｓｕｂ３によって送信された出力信号のための共通の時間遅延は、Ｓｕｂ３からの出力信号がその最も近いマイクロフォンＰ４に到達するのにかかる時間である。全ての計測されたインパルス応答の最小の遅延が共通の時間遅延である。共通の時間遅延は、四つのマイクロフォンによって計測された全てのインパルス応答からオフセットされ得、または、差し引かれ得る。

ブロック６０６において、部屋の時間領域インパルス応答の入力データは、高速フーリエ変換（ＦＦＴ）を使用し、周波数領域へと変換され得る。例えば図１において、１６通りのインパルス応答（それぞれのセットがＮ数のサンプルを有する）が計測され得るように、四つのマイクロフォンおよび四つのラウドスピーカがある。この例において、それぞれの計測されたインパルス応答のＮ複素値を生み出すＮポイントＦＦＴが用いられる。それゆえ、｛ｎ_ｍｉｃｘｎ_ｓｕｂ｝×Ｎ複素ＦＦＴポイントの結果は、Ｎ数のｎ_ｍｉｃ×ｎ_ｓｕｂマトリクスＡ_ｉとして表される。ここで、ｉ＝１・・・Ｎである。それぞれのｉまたは周波数ポイントにて、ＦＦＴは、振幅および位相を提供する。

図７は、図５におけるブロック５０４によって表された周波数応答の反転を決定する方法をさらに詳述するフローチャート７００である。ブロック７０２において、インパルス応答を計測するために使用されるマイクロフォンｎ_ｍｉｃの数、および、部屋におけるサブウーハーｎ_ｓｕｂの数が決定される。決定ブロック７０４において、ｎ_ｍｉｃ＝ｎ_ｓｕｂである場合、次いで、ブロック７０６において、正確なマトリクス反転（ｍａｔｒｉｘ ｉｎｖｅｒｓｉｏｎ）方法が、インパルス応答の正確な反転を見出すために使用され得る。他方で、ｎ_ｍｉｃ＞ｎ_ｓｕｂである場合、次いで、ブロック７０８において、擬似反転（ｐｓｅｕｄｏ−ｉｎｖｅｒｓｅ）方法が、インパルス応答の反転を見出すために使用され得る。図１において、正確なマトリクス反転方法がその反転を計算するために使用されるように、四つのマイクロフォンおよび四つのサブウーハーがインパルス応答を計測するよう、使用される。ブロック６０４において周波数領域へ変換されるインパルス応答を有し、反転マトリクスは、周波数ポイントでの理想的なイコライゼーションを決定するために、その周波数ポイントのそれぞれにおいて計算され得る。これの関連で、Ｎ数の反転マトリクスＢ_ｉが決定され得る。ここで、ｉ＝１・・・Ｎである。これは、Ｎ複素値のマトリクスＢ_ｉという結果となり、Ａ_ｉＢ_ｉ＝１となる。

ｎ_ｍｉｃ＞ｎ_ｓｕｂの場合において、擬似反転方法はＢ_ｉを計算するために使用され得る。擬似反転の既知の方法は、所望される結果と実際の結果との間の平均平方誤差を最小化する。数学的に表現すると、Ｂ_ｉは、（１−Ａ_ｉＢ_ｉ）^＊ｘ（１−Ａ_ｉＢ_ｉ）が最小化されるように計算される（^＊は複素共役演算を示す）。

ブロック７１０において、一度、反転マトリクスが決定されると、目的関数は、マイクロフォンの位置Ｐ１からＰ４のそれぞれの周波数位置のそれぞれのために選択され得る。目的関数は、それぞれのリスニング位置での所望される周波数応答である。目的関数は、ｎ_ｍｉｃ要素Ｔ_ｉ（ｉ＝１・・・Ｎ）を含む複素値ベクトルであり得る。四つのマイクロフォンのこの例において、Ｔ_ｉは、周波数ポイント毎の四つの複素値要素を含む。特定の周波数ポイントｉ（ｉ＝１・・・Ｎ）でのｎ_ｓｕｂフィルタを記載するベクトルＦ_ｉは、次いで、マトリクス乗算Ｆ_ｉ＝Ｂ_ｉＴ_ｉとして計算される。ベクトルＦ_ｉは、正確な反転（理想的なイコライゼーション）を実行する、特定の周波数ポイントｉ（ｉ＝１・・・Ｎ）でのフィルタを記載する。Ａ_ｉＦ_ｉ＝Ａ_ｉＢ_ｉＴ_ｉ＝Ｔ_ｉを乗算することが理想的なイコライゼーションという結果になるように、ベクトルＦ_ｉは、部屋の壁によって追加された音色を事実上取り消す。

図８は、マトリクス反転の後に得られたフィルタＦ（ｋ）（ｋ＝１・・ｎ_ｓｕｂ＝４）の対数振幅を示すグラフである。この例において利用される目的関数は、単一のベクトルＴ_ｉ＝[１１１１]であり、ｉ＝１・・Ｎである。周波数軸ｆはｆ＝（１・・・Ｎ／２）／Ｎ^＊ｆａであり、ＮはＦＦＴの長さであり、ｆａ＝７５０Ｈｚはサンプル周波数である。図８は、例えば、マークＡ、Ｂ、Ｃ、およびＤによって示されたような、唐突なピークおよび窪みがあることを示す。部屋内で周波数応答をイコライズするために、Ｓｕｂ１からＳｕｂ４へ送信された出力信号にフィルタＦ（ｋ）を直接的に適合させることは、サブウーハーに損傷を与え得る。というのも、ある周波数でのピークは、扱うサブウーハーにとって高すぎ得るそれらの周波数での相当なゲインを適合させることを要求するからである。つまり、ベクトルＦ（ｋ）は、サブウーハーが扱い得る最大ゲインを上回り得る周波数でのゲインを押し付け得る。

周波数範囲全体にわたる平滑化は、時間領域におけるその結果のフィルタの長さを限定するためになされ得、それは、平滑化の後、より迅速にゼロに集結することが知られている。以下は、図５におけるブロック５０６によって表されたマトリクスの反転を平滑化することについてのさらなる検討である。周波数応答ベクトルＦ（ｋ）における唐突なピークおよび窪みとともに、理想的なイコライゼーションは、サブウーハーに送信される出力信号へ直接に適合され得ない。しかしながら、ベクトルＦ（ｋ）におけるピークおよび窪みは、周波数を介した複素値ベクトルＦ（ｋ）を平滑化することによって最小化され得る。これは、２０００年４月、Ｊ．Ａｕｄｉｏ Ｅｎｇ．Ｓｏｃ．， Ｖｏｌ．４８， Ｎｏ．４， ｐｐ２５９−２８０にて出版された、Ｐａｎａｇｉｏｔｉｓ Ｄ． ＨａｔｚｉａｎｔｏｎｉｏｕおよびＪｏｈｎ Ｎ． Ｍｏｕｒｊｏｐｏｕｌｏｓによる、「Ｇｅｎｅｒａｌｉｚｅｄ Ｆｒａｃｔｉｏｎａｌ−Ｏｃｔａｖｅ Ｓｍｏｏｔｈｉｎｇ ｏｆ Ａｕｄｉｏ ａｎｄ Ａｃｏｕｓｔｉｃ Ｒｅｓｐｏｎｓｅｓ」と題された論文において記載された方法を介して達成され得る。特に、複素値ベクトルＦ（ｋ）の平滑化は、周波数依存のスライドウィンドウに沿って、実際の部分と架空の部分の平均値を別個に計算することによって実行され得、その結果Ｆｓ（ｋ）を得る。例えば、１．０と２．０との間の平滑化インデックスｑが使用され得、ｉ^＊（ｑ−１／ｑ）は、周波数依存のスライドウィンドウの幅を指す。ＨａｎｎｉｎｇウィンドウまたはＷｅｌｃｈウィンドウなどのようなスライドウィンドウが使用され得る。それぞれの周波数帯域に異なる値を使用することによって、二つ以上の別個の周波数帯域における平滑化を実行することは有益であり得る。より高い周波数では、部屋における空間および周波数を介した変動は、通常、より大きいため、より高いインデックスｑが利用され得る。サブウーハーは主に８０Ｈｚ以下で動作するゆえ、その周波数より上の反転フィルタの高い正確性は必要ではあり得ず、所望され得ない。なぜならば、迅速な変動のために、それは一貫して全体のリスニング領域に適合し得ないからである。

図９は、Ｓｕｂ２へ送信された出力信号に適合され得るフィルタＦ（２）の非平滑（ｕｎｓｍｏｏｔｈｅｄ）スペクトルの振幅、および、上で検討された方法でのフィルタＦ（２）の平滑化されたバージョンを表す曲線Ｆｓ（２）の振幅を示す。曲線Ｆｓ（２）において、局地的なピークおよび窪みは曲線Ｆ（２）におけるよりも、より平滑化されており、曲線Ｆ（２）に存在する唐突なピークおよび窪みの多くは、曲線Ｆｓ（２）において、より段階的である。それゆえ、イコライゼーションが局地的な超過ゲインなしでＳｕｂ２への出力信号に適合され得るために、曲線Ｆｓ（２）は、複素値フィルタＦ（２）の近似値である。同様に、図１０は、例えば、Ｆｓ（１）、Ｆｓ（２）、Ｆｓ（３）、およびＦｓ（４）などのような平滑化後の、四つ全てのフィルタの振幅応答の曲線を表す。

図１１は、図５におけるブロック５０８によって表される広範囲イコライゼーションを決定する方法をさらに詳述するフローチャート１１００を示す。Ｆ（１）からＦ（４）の複素値フィルタのそれぞれの複素平滑化は、ピークおよび窪みのローカルな変動を取り除くが、極端なゲインはまだ存在し得る。例えば、サブウーハーは、一般に、およそ１５ｄｂからおよそ２０ｄｂの最大ゲインを扱うようにデザインされる。図９は、２０Ｈｚ以下のおよそ３０ｄｂのゲイン、および、１００Ｈｚ以上のおよそ６０ｄｂのゲインを示す。そのような極端なゲインは、サブウーハーでは扱い得ない。

ゲインを管理するために、広範囲イコライゼーション（ＥＱ）が実行され得る。広範囲ＥＱを計算する一つの方法は、図１１に記載された方法を介すことである。ブロック１１０２において、マイクロフォンの位置または座席Ｆｙ（ｊ）（ｊ＝１・・・ｎ_ｓｅａｔ）のそれぞれでの実際の応答が、Ｆｓを伴う元のマトリクスＡを乗算することによって計算され得る（上記の平滑化方法において計算される）。つまり、Ｆｙ＝Ａ^＊Ｆｓである。図１２は、図１０の（中間の）フィルタが適合された後、四つのマイクロフォンの位置（リスナー席）での応答を示す。ブロック１１０４において、上部の曲線Ｆｙｍａｘは、それぞれの周波数ポイントの最大振幅Ｍａｘ｛Ｆｙ（１・・・ｎ_ｓｅａｔ）｝を取ることによって決定され得る。それゆえ、座席での全ての応答は、曲線Ｆｙｍａｘより下である。図１２は、Ｆｙｍａｘ曲線をより良く示すために、１０ｄＢによって生じた曲線Ｆｙｍａｘを示す。これは、任意の周波数ポイントに沿って、曲線Ｆｙｍａｘより大きい応答がないことを意味する。

曲線Ｆｙｍａｘは、０Ｈｚからサンプルレートの半分の周波数範囲全体内でのｄＢの最大振幅を示す。しかしながら、サブウーハーは、上記の周波数範囲より、より限定された範囲において最適に動作するような設定である。それゆえ、ブロック１１０６において、上部の曲線Ｆｙｍａｘは、サブウーハーに最適な周波数範囲で動作させる所定の周波数範囲内に限定され得る。この関連において、広範囲ＥＱフィルタＦｒは、目的関数ＴをＦｙｍａｘで割る、または、Ｆｒ＝Ｔ／Ｆｙｍａｘによって、所定の周波数範囲において動作するように計算され得る。目的関数Ｔは、個々の変換器（サブウーハ−）が最適に動作する周波数範囲を特徴付ける高域通過（ｈｉｇｈ ｐａｓｓ）フィルタおよび低域通過（ｌｏｗ ｐａｓｓ）フィルタの振幅周波数応答を有する実数値である。通常のフィルタは、ｎ＝２・・４（コーナー周波数２０・・４０Ｈｚ）の順のＢｕｔｔｅｒｗｏｒｔｈ高域通過、および、ｎ＝２・・４（コーナー周波数８０・・１５０Ｈｚ）の順のＢｕｔｔｅｒｗｏｒｔｈ低域通過である。

図１３は、広範囲ＥＱフィルタＦｒのログ振幅応答を示す。図１３は、その応答が音質を干渉し得るピークを有することを示す。この関連において、ブロック１１０８において、曲線Ｆｒにおけるピークは、以下の方法を介して取り除かれ得る。上に記載された平滑化の方法は、Ｆｒの平滑化バージョンである中間応答Ｆｒｓを決定するために使用され得る。Ｆｒにおけるピークは実質的に、ＦｒｓとＦｒとの間の差異を計算し、その差異を修正することによって、「削り取り（ｓｈａｖｅｄ ｏｆｆ）」得る。図１４は、ＦｒｓとＦｒとの間の差異（１０ｄｂ分生じた）を表す上の曲線、および、修正された差異（１０ｄｂ分下げられた）を表す、下の曲線を示す。次いで、図１５に示されるように、広範囲ＥＱフィルタＦｒｓｆの最終周波数応答は、修正された差異から、図１３において示された非平滑フィルタである元のフィルタＦｒを差し引くことによって得られ得る。図１５において示される最終Ｆｒｓｆは、窪みと、数が減少したピークを示す。無用のピークは、窪みが元の応答において生じる周波数を増幅することを企て、サブウーハーからの相当な追加的音声出力を要求し、その結果、システムの最大音響信号を低減し、部屋の他の領域における大きなピークを潜在的に生み出してしまう。

図１６は、図５におけるブロック５１０によって表されるように、広範囲ＥＱ曲線上での最大ゲインを制限する方法をさらに詳述するフローチャート１６００を示す。ブロック１６０２において、最終ＥＱスペクトルＦｅｑは、上に記載されたように、個々に広範囲、実数値振幅スペクトルＦｒｓｆを用いて、個々のＥＱフィルタの複素スペクトルＦｓを乗算することによって計算される。図１７は、複素平滑化および広範囲ＥＱの後に得られたＥＱフィルタを示す。図１７は、実質的に、２００Ｈｚより上のゲインおよびおよそ２０Ｈｚより下のゲインが存在することを示す。これは、所望されるような最終のゲインを制限するのに十分ではない選択された目的関数のためであり得る。それゆえ、ブロック１６０４において、制限は、所定の低周波数および所定の高周波数を下回るゲインに与えられ得る。例えば、最大ゲイン上の制限は、複素値Ｆｅｑを置き換えることによって適合され得るため、最大振幅は、その位相を変化させることなしに、「Ｍａｘｇａｉｎ」にクリップされる。Ｍａｘｇａｉｎは、特定のサブウーハーの能力に依存する、ユーザによって定められた値である。好ましくは、Ｍａｘｇａｉｎの異なる値は、異なる周波数帯域において適合され得る。その結果としてのフィルタは、最大ゲインが０ｄＢを超過しないように、縮小され得る。図１８は、上記で示されたように、Ｍａｘｇａｉｎおよび０ｄＢへの正常化を適合した後の、フィルタＥＱスペクトルを示す。ＥＱスペクトルは、平均ゲインを最大化するために０ｄＢへと正常化される。

ブロック１６０６において、最終ＥＱフィルタ周波数応答は、反転ＦＦＴを利用することによって、時間領域へと変換され得、その結果、Ｆｉｎｉｔｅ Ｉｍｐｕｌｓｅ Ｒｅｓｐｏｎｓｅ（ＦＩＲ）フィルタの係数となる。時間窓は、フィルタの長さを制限するためにその係数に適合され得る。図１９は、得られたＦＩＲフィルタ（フィルタＦＩＲ１）の一つのインパルス応答を示す。図２０は、その結果としてのフィルタＦＩＲ１からＦＩＲ４の振幅応答を示す。上で検討されたように、図４は、得られたＥＱフィルタを適合した後に、四つの座席Ｐ１からＰ４での、その結果としての応答を示す。およそ４０Ｈｚと８０Ｈｚとの間のような目的周波数範内において、応答は、その周波数範囲内での実質的なイコライゼーションを提供するゆえに一貫し、フラットであることに注意したい。これは、Ｐ１からＰ４の位置の任意の一つに着席した個人が、実質的に、バスサウンドの類似した大きさのレベルを聴くことを意味する。つまり、サウンドレベルは、それぞれの個人が同じバスサウンドの音質を経験するように、部屋のリスニング領域内での異なる位置においても実質的には同じである。付け加えて、図４は、曲線が当該周波数範囲内で、実質的にフラットであることを示す。これは、所望される周波数範囲内のバスサウンドが、最小であるか、可能である場合、なくしてしまうために、バスサウンドがその所望される周波数範囲内で、実質的に一貫していることを意味する。

上で記載されたイコライゼーションシステムは、少なくとも一つのサブウーハ−を有し、異なる構成を有する様々な部屋に利用され得る。部屋は、ラウドスピーカが置かれる任意の空間のタイプを含み得る。その空間は、ドアが閉められた部屋または車内などのような十分に閉じられた境界、または、連結された廊下を有する部屋、開けっ放しの戸口または、開いた壁などのような部分的に閉じられた境界、あるいは、サンルーフの開いた車両、などを有し得る。付け加えて、部屋は、野外や、天上が閉じられている、または開いているスタジアムなどのような開かれた領域であり得る。空間における低周波数のパフォーマンスは、明細書および添付された請求項における部屋に関連して記載されるが、しかしながら、車内、レコーディングスタジオ、住居のリビング空間、コンサートホール、映画劇場、部分的に閉じられた空間などもまた含まれることは理解されよう。部屋の壁などの境界は、部分的にか十分に部屋を囲い込む間仕切りを含む。部屋の境界は、石膏、木材、コンクリート、ガラス、皮、繊維、プラスチックなどのような任意の素材から作られ得る。住居においては、部屋の境界は、しばしば、石膏、石造、または繊維から作られ得る。境界は、壁、カーテン、家具、備え付け家具などを含み得る。部屋の境界は、音を反射、拡散、および吸収する様々な性能を有する。部屋の境界の音声特徴は、音響信号に影響を与え得る。

ラウドスピーカは、様々な形およびサイズであり得る。例えば、ラウドスピーカは、変換器を覆う箱型の構成において同梱され得る。ラウドスピーカはまた、筐体の全てまたは部分として、壁または車両の一部を利用し得る。ラウドスピーカは、低周波から高周波までの音声周波数の全ての範囲を提供し得る。多くのラウドスピーカは、筐体の中の多数の変換器を有する。多数の変換器がラウドスピーカ筐体において利用される場合、個々の変換器が異なる周波数帯域において、さらに効果的に動作することは一般的である。ラウドスピーカまたはラウドスピーカの一部は、低周波数などのような、音声周波数の特定の範囲を提供するために最適化され得る。ラウドスピーカは、専用アンプ、ゲイン制御、イコライザーなどを含み得る。ラウドスピーカは、いくつかの、または追加的な構成要素を含む他の構成要素を有し得る。

低周波数を生成するように最適化される変換器を含むラウドスピーカ、またはラウドスピーカの一部は、通常、サブウーハーと呼ばれる。サブウーハーは、低周波数を生成することができる任意の変換器を含み得る。低周波数を生成することができるラウドスピーカは、明細書および添付された請求項においてサブウーハーという言葉によって参照され得るが、低周波数を生成することができ、通常の電気信号に応答することができる任意のラウドスピーカまたは任意のラウドスピーカの一部が含まれる。

マイクロフォンなどのような計測デバイスは、部屋の様々な部分における音響信号を計測するために、信号ブロック３０２などのような他の電気信号と通信し得る。異なるリスニング位置のための異なるラウドスピーカの位置からの、計測された音響信号出力は、外部ディスク上などにおいて格納され得る。外部ディスクは、コンピュータデバイスへ入力され得る。コンピュータデバイスは、他のコンピュータ環境であり得、計測デバイスに関連し、上で記載された要素の多くまたは全てを含み得る。コンピュータデバイスは、部屋だけでなく、異なる位置におけるインパルス応答を処理するために、部屋内に位置された、または、遠隔に位置された、音声／映像受信器に組み込まれ得る。

本発明の様々な実施形態が記載されている一方で、さらに多くの実施形態および実施は、本発明の範囲内において可能であることは、当業者にとって明らかである。従って、本発明は、添付された請求項およびそれらの均等物による以外に、限定されるものではない。

サブウーハーが部屋のそれぞれの隅にある長方形の部屋、および、Ｐ１からＰ４によって遮られるリスニング領域における、通常のサウンドシステムの設定を示す。 個々に、四つのマイクロフォンの位置Ｐ１、Ｐ２、Ｐ３、およびＰ４での期待し得る、計測されたインパルス応答に対応する、四つのスペクトルＦ１、Ｆ２、Ｆ３、およびＦ４を示す。 本発明に従ったイコライゼーションシステムを図示するブロック図を示す。 対応するサブウーハーへの入力信号が周波数応答をイコライズするようにフィルタされた後の部屋の、本発明にしたがった周波数応答を示す。 部屋の周波数応答をイコライズするためのフィルタデザイン手順の概観を有するフローチャートである。 図５における入力データステップを用意する、さらなる詳細を示すフローチャートである。 図５における周波数応答の反転を決定する、さらなる詳細を示すフローチャートである。 周波数応答の反転を表す曲線を示す。 Ｆ（２）の平滑化バージョンを表す、本発明に従った曲線Ｆｓ（２）を示す。 個々に、図８における曲線Ｆ（１）、Ｆ（２）、Ｆ（３）、およびＦ（４）の平滑化バージョンを表す四つの曲線Ｆｓ（１）、Ｆｓ（２）、Ｆｓ（３）、およびＦｓ（４）を示す。 図５における広範囲イコライゼーションステップを決定する、さらなる詳細を示すフローチャートである。 図１０において示された、個々に、曲線Ｆｓ（１）、Ｆｓ（２）、Ｆｓ（３）、およびＦｓ（４）に従ったフィルタが適合された後、四つのマイクロフォン位置Ｐ１、Ｐ２、Ｐ３、およびＰ４での周波数応答を示す。 反転された広範囲イコライゼーションフィルタを示す。 図１３において、１０ｄＢによって生じた、平滑化または非平滑された周波数応答との間の差異を表す、上の曲線、および、修正された差異（１０ｄＢによって下げられた）を表す、下の曲線を示す。 広範囲イコライゼーションフィルタの最終周波数応答を示す。 図５において示されたように、広範囲イコライゼーションフィルタにおける最大ゲインを制限するステップをさらに詳述するフローチャートを示す。 図１０において示される曲線Ｆｓ（１）、Ｆｓ（２）、Ｆｓ（３）、およびＦｓ（４）の複素平滑化、ならびに、図１５において示される広範囲イコライゼーションフィルタを平滑化された曲線Ｆｓ（１）、Ｆｓ（２）、Ｆｓ（３）、およびＦｓ（４）に適合した後の、サブウーハーのそれぞれのイコライゼーションフィルタを示す。 上で示されたように、Ｍａｘｇａｉｎおよび０ｄＢへの正常化を適合した後の、フィルタＥＱスペクトルを示す。 フィルタＦＩＲ１のために得られた、対応するイコライズされたインパルス応答を示す。 フィルタＦＩＲ１、ＦＩＲ２、ＦＩＲ３、およびＦＩＲ４のための振幅応答を示す。

符号の説明

３００ イコライゼーションシステム（ＥＱシステム）
３０２ 信号ブロック
３０４ プロセッサブロック
３０６ 音声入力信号
３０８ 受信器

Claims (14)

1. リスニング領域における周波数範囲内で周波数応答を実質的にイコライズするように一つ以上のフィルタをデザインする方法であって、該方法は、
   少なくとも一つのマイクロフォンによってリスニング領域における周波数応答を計測すること（５０２）であって、該周波数応答は、少なくとも一つのサブウーハーに含まれる少なくとも一つの音響変換器によって生成される、ことと、
   該周波数応答を反転することにより、第１のステージのイコライゼーションフィルタスペクトルを決定すること（５０４）と、
   該第１のステージのイコライゼーションフィルタスペクトルを平滑化することにより、第２のステージの近似のイコライゼーションフィルタスペクトルを決定すること（５０６）と、
   該少なくとも一つのマイクロフォンの最大の実際の応答を計算すること（１１０２、１１０４）であって、該計算することは、該第２のステージの近似のイコライゼーションフィルタを該計測された周波数応答に適用することを含む、ことと、
   該第２のステージの近似のイコライゼーションフィルタスペクトルにおける極端なゲインを平滑化する広範囲イコライゼーションフィルタを決定すること（１１０６）であって、該広範囲イコライゼーションフィルタを決定することは、該計算された最大の実際の応答を反転することを含む、ことと、
   該広範囲イコライゼーションフィルタスペクトルにおけるピークの数を低減するために、該広範囲イコライゼーションフィルタと該第２のステージの近似のイコライゼーションフィルタとを結合することにより、最終イコライゼーションフィルタを決定すること（１６０２）と
   を含む、方法。
2. 前記周波数応答を計測することは、
   前記リスニング領域内に位置された前記少なくとも一つのマイクロフォンを介して、前記少なくとも一つのサブウーハーに含まれる前記少なくとも一つの音響変換器によって生成された、前記リスニング領域のインパルス応答を受信すること（６０２）と、
   該リスニング領域における該インパルス応答から任意の共通する時間遅延を取り除くこと（６０４）と、
   該リスニング領域の該インパルス応答を該リスニング領域における該周波数応答に変換すること（６０６）と
   を含む、請求項１に記載の方法。
3. 前記最終イコライゼーションフィルタの振幅応答をクリップすることにより、所望される周波数帯域内にゲインを制限すること（１６０４）を含む、請求項１または２に記載の方法。
4. 前記周波数応答を計測することは、
   前記リスニング領域内に位置された前記少なくとも一つのマイクロフォンを介して、前記少なくとも一つのサブウーハーに含まれる前記少なくとも一つの音響変換器によって生成された、前記リスニング領域のインパルス応答を受信すること（７０２）と、
   前記リスニング領域における前記少なくとも一つの音響変換器の数と同じ数の、部屋における少なくとも一つのマイクロフォンを使用することにより、前記インパルス応答を計測すること（７０６）と
   を含む、請求項１から３のうちのいずれか一項に記載の方法。
5. 複素平滑化方法を介して、前記周波数応答の反転がなされる（５０６）、請求項１から４のうちのいずれか一項に記載の方法。
6. 前記複素平滑化は、二つの異なる平滑化インデックス値を有する二つの別個の周波数帯域にてなされる（５０６）、請求項５に記載の方法。
7. 前記少なくとも一つの前記音響変換器の数は四つであり、前記少なくとも一つのマイクロフォンの数は四つである（７０６）、請求項４に記載の方法。
8. 前記周波数応答を計測することは、
   前記リスニング領域内に位置された前記少なくとも一つのマイクロフォンを介して、前記少なくとも一つのサブウーハーに含まれる前記少なくとも一つの音響変換器によって生成された、前記リスニング領域のインパルス応答を受信すること（７０２）
   を含み、
   該インパルス応答を計測するために利用される少なくとも一つのマイクロフォンの数が、少なくとも一つの音響変換器の数と等しくない場合、該周波数応答の反転が、擬似反転法を介してなされる（７０８）、請求項１から７のうちのいずれか一項に記載の方法。
9. 目的関数を前記周波数応答に適用することにより、前記リスニング領域における前記少なくとも一つの音響変換器のそれぞれが動作する周波数範囲を制限すること（１１０６）を含む、請求項１から８のうちのいずれか一項に記載の方法。
10. 部屋内のリスニング領域における周波数範囲内で周波数応答を実質的にイコライズするように一つ以上のフィルタをデザインするシステムであって、
    該システムは、
    出力信号を生成する（５０２）少なくとも一つの音響変換器を含む少なくとも一つのサブウーハーと、
    該リスニング領域における周波数応答を計測する（５０２）少なくとも一つのマイクロフォンと、
    該音響変換器と該マイクロフォンとに結合されたプロセッサ（３０２）と
    を含み、
    該プロセッサは、テスト信号を送信することにより、該音響変換器が、該出力信号を生成することと、該マイクロフォンを介して該部屋の音響応答を計測することとを可能にすることが可能であり、該プロセッサは、該音響変換器への該フィルタされた出力信号が、該部屋の該リスニング領域における周波数範囲内の実質的にフラットな周波数応答を生成するように、該計測された音響応答に基づいて、フィルタ係数を計算することが可能であり、
    該プロセッサは、
    該計測された音響応答を反転することにより、第１のステージのイコライゼーションフィルタスペクトルを決定すること（５０４）と、
    該第１のステージのイコライゼーションフィルタスペクトルを平滑化することにより、第２のステージの近似のイコライゼーションフィルタスペクトルを決定すること（５０６）と、
    該少なくとも一つのマイクロフォンの最大の実際の応答を計算すること（１１０２、１１０４）であって、該計算することは、該第２のステージの近似のイコライゼーションフィルタを該計測された周波数応答に適用することを含む、ことと、
    該第２のステージの近似のイコライゼーションフィルタスペクトルにおける極端なゲインを平滑化する広範囲イコライゼーションフィルタを決定すること（１１０６）であって、該広範囲イコライゼーションフィルタを決定することは、該計算された最大の実際の応答を反転することを含む、ことと、
    広範囲イコライゼーションフィルタスペクトルにおけるピークの数を低減するために、該広範囲イコライゼーションフィルタと該第２のステージの近似のイコライゼーションフィルタとを結合することにより、最終イコライゼーションフィルタを決定すること（１６０２）と
    を実行することによって、該フィルタ係数を計算する、システム。
11. 前記部屋における前記少なくとも一つの音響変換器の数は、該部屋における前記少なくとも一つのマイクロフォンの数と等しい（７０６）、請求項１０に記載のシステム。
12. 前記部屋における前記少なくとも一つの音響変換器の数は四つであり、該部屋における前記少なくとも一つのマイクロフォンは四つである（７０６）、請求項１０または１１に記載のシステム。
13. 前記プロセッサはデジタル信号プロセッサである、請求項１０から１２のうちのいずれか一項に記載のシステム。
14. 前記プロセッサは、前記フィルタ係数を計算するために命令を実行することが可能であり、
    該命令は、
    少なくとも一つのマイクロフォンによってリスニング領域における周波数応答を計測する（５０２）ための命令であって、該周波数応答は、少なくとも一つのサブウーハーに含まれる少なくとも一つの音響変換器によって生成される、命令と、
    前記周波数応答を反転することにより、第１のステージのイコライゼーションフィルタスペクトルを決定する（５０４）ための命令と、
    該第１のステージのイコライゼーションフィルタスペクトルを平滑化することにより、第２のステージの近似のイコライゼーションフィルタスペクトルを決定する（５０６）ための命令と、
    該少なくとも一つのマイクロフォンの最大の実際の応答を計算する（１１０２、１１０４）ための命令であって、該計算することは、該第２のステージの近似のイコライゼーションフィルタを該計測された周波数応答に適用することを含む、命令と、
    該第２のステージの近似のイコライゼーションフィルタスペクトルにおける極端なゲインを平滑化する広範囲イコライゼーションフィルタを決定する（１１０６）ための命令であって、該広範囲イコライゼーションフィルタを決定することは、該計算された最大の実際の応答を反転することを含む、命令と、
    広範囲イコライゼーションフィルタスペクトルにおけるピークの数を低減するために、該広範囲イコライゼーションフィルタと該第２のステージの近似のイコライゼーションフィルタとを結合することにより、最終イコライゼーションフィルタを決定する（１６０２）ための命令と
    を含む、請求項１０から１３のうちのいずれか一項に記載のシステム。

Images