(57)【特許請求の範囲】
【請求項１】
多チャンネルオーディオシステムにおいて、多チャンネルラウドスピーカ構成を特徴付ける方法であって、
第１のプローブ信号を生成する段階と、
前記第１のプローブ信号を第１の音響応答に変換し、かつ該音響応答を無音期間によって分離された非重ね合わせ時間スロットで音波としてリスニング環境内に順次送信するために、該リスニング環境に多チャンネル構成で位置決めされたそれぞれの電気−音響変換器に結合された複数のオーディオ出力に該第１のプローブ信号を供給する段階と、
を含み、
更に、各前記オーディオ出力に対して、
各々が前記音響応答を第１の電気応答信号に変換する少なくとも２つの非共在音響−電気変換器を含む多マイクロフォンアレイにおいて音波を受信する段階と、
各前記音響−電気変換器における前記電気−音響変換器に対する第１の室内応答を決定するよう、前記第１の電気応答信号に前記第１のプローブ信号を用いて逆畳み込みを行う段階と、
各前記音響−電気変換器における前記電気−音響変換器に対する遅延を計算してメモリに記録する段階と、
各前記音響−電気変換器における前記電気−音響変換器に対する前記遅延によってオフセットされた指定期間にわたって前記第１の室内応答をメモリに記録する段階と、
を含み、
方法が、更に、
各前記音響−電気変換器までの前記遅延に基づいて、各前記電気−音響変換器に対する距離及び少なくとも第１の角度を決定する段階と、
前記電気−音響変換器に対する前記距離及び少なくとも前記第１の角度を使用して、特定の多チャンネル構成を自動的に選択し、かつ前記リスニング環境内のその多チャンネル構成における各電気−音響変換器に対する位置を計算する段階と、
前記多チャンネルオーディオシステムのための室内補正フィルタを発生させる段階とを含み、
前記室内補正フィルタを発生させる段階は、
Ｐが整数である時に、オーディオ信号をＰ個のサブバンドに対するベースバンドまでダウンサンプリングするＰバンドオーバーサンプリング解析フィルタバンクと、該オーディオ信号を再構成するよう該Ｐ個のサブバンドをアップサンプリングするＰバンドオーバーサンプリング合成フィルタバンクとを与える段階と、
各チャンネルに対してスペクトル測度を与える段階と、
各前記スペクトル測度をチャンネルターゲット曲線と組み合わせてチャンネル毎の集合スペクトル測度を与える段階と、
少なくとも１つのチャンネルについて、前記集合スペクトル測度のうちで異なるサブバンドに対応する部分を抽出する段階と、
前記少なくとも１つのチャンネルについて、前記解析フィルタバンクの前記ダウンサンプリングを模擬するよう前記集合スペクトル測度の前記抽出部分をベースバンドに再マップする段階であって、偶数及び奇数のサブバンドに対応する前記抽出部分を再マップする前記段階は、前記抽出部分をベースバンドに周波数シフトし、該シフトされた周波数部分をマイナス及びプラス９０度でそれぞれ転換し、これにより、転換された該スペクトルの一部を急転させ、各サブバンドに対する再マップされた前記スペクトル測度における不連続部を生成する段階を含む、段階と、
前記少なくとも１つのチャンネルについて、各サブバンドに対して前記再マップされたスペクトル測度に対する自己回帰（ＡＲ）モデルを推定する段階と、
前記少なくとも１つのチャンネルについて、各前記ＡＲモデルの係数を最小位相全ゼロサブバンド補正フィルタの係数にマップする段階と、
前記解析及び合成フィルタバンクの間でＰ個のベースバンドオーディオ信号を周波数補正するＰ個のデジタル全ゼロサブバンド補正フィルタを対応する係数から構成する段階と
によって行われることを特徴とする方法。
【請求項２】
前記遅延を計算する前記段階は、
各前記第１の電気応答信号及び前記第１のプローブ信号を処理して時間シーケンスを発生させる段階と、
前記オーディオ出力が前記電気−音響変換器に結合されているか否かを示すものとして前記時間シーケンスにおける顕著なピークの存在又は不在を検出する段階と、
前記ピークの位置を前記遅延として計算する段階と、
を含む、
ことを特徴とする請求項１に記載の方法。
【請求項３】
前記第１の電気応答信号は、該第１の電気応答が前記音響−電気変換器で受信された時にブロックに区分化され、かつ前記第１のプローブ信号の区分を用いて逆畳み込みが施され、
前記遅延及び第１の室内応答は、次のプローブ信号の送信の前の前記無音期間に計算され、かつメモリに記録される、
ことを特徴とする請求項１に記載の方法。
【請求項４】
前記第１のプローブ信号の後に前記複数のオーディオ出力の各々に第２のプリ−エンファシスプローブ信号を供給し、第２の電気応答信号を記録する段階と、
前記第１のプローブ信号の前記区分を用いて前記第２の応答信号の重ね合わせブロックに逆畳み込みを行い、第２の候補室内応答のシーケンスを発生させる段階と、
前記第１のプローブ信号に対する前記遅延を使用して連続する第２の候補室内応答を付加し、第２の室内応答を形成する段階と、
を更に含むことを特徴とする請求項１に記載の方法。
【請求項５】
多チャンネルオーディオを処理するためのデバイスであって、
リスニング環境に多チャンネル構成で位置決めされてオーディオ出力に結合されたそれぞれの電気−音響変換器を駆動するための複数のオーディオ出力と、
オーディオ入力に結合された複数の音響−電気変換器から第１の電気応答信号を受信するための１つ又はそれよりも多くのオーディオ入力と、
前記複数の第１の電気応答信号を受信するよう前記１つ又はそれよりも多くのオーディオ入力に結合された入力レシーバと、
デバイスメモリと、
第１のプローブ信号を生成し、かつ
前記第１のプローブ信号を無音期間によって分離された非重ね合わせ時間スロットで前記複数のオーディオ出力の各々に供給する、
ようになったプローブ発生及び送信スケジューリングモジュール、及び
各前記オーディオ出力に対して、前記第１の電気応答信号に前記第１のプローブ信号を用いて逆畳み込みを行って各前記音響−電気変換器での第１の室内応答を決定し、各該音響−電気変換器での遅延を計算して前記デバイスメモリに記録し、かつ各該音響−電気変換器での該遅延によってオフセットされた指定期間にわたって該第１の室内応答を該デバイスメモリに記録し、
各前記電気−音響変換器に対する各前記音響−電気変換器での前記遅延に基づいて、該電気−音響変換器に対する距離及び少なくとも第１の角度を決定し、かつ
前記電気−音響変換器に対する距離及び少なくとも前記第１の角度を使用して特定の多チャンネル構成を自動的に選択し、かつ前記リスニング環境内のその多チャンネル構成における各電気−音響変換器に対する位置を計算し、
前記多チャンネルオーディオシステムのための室内補正フィルタを発生させ、
前記室内補正フィルタの発生は、
Ｐが整数である時に、オーディオ信号をＰ個のサブバンドに対するベースバンドまでダウンサンプリングするＰバンドオーバーサンプリング解析フィルタバンクと、該オーディオ信号を再構成するよう該Ｐ個のサブバンドをアップサンプリングするＰバンドオーバーサンプリング合成フィルタバンクとを与えること、
各チャンネルに対してスペクトル測度を与えること、
各前記スペクトル測度をチャンネルターゲット曲線と組み合わせてチャンネル毎の集合スペクトル測度を与えること、
少なくとも１つのチャンネルについて、前記集合スペクトル測度のうちで異なるサブバンドに対応する部分を抽出すること、
前記少なくとも１つのチャンネルについて、前記解析フィルタバンクの前記ダウンサンプリングを模擬するよう前記集合スペクトル測度の前記抽出部分をベースバンドに再マップすることであって、偶数及び奇数のサブバンドに対応する前記抽出部分を再マップすることは、前記抽出部分をベースバンドに周波数シフトし、該シフトされた周波数部分をマイナス及びプラス９０度でそれぞれ転換し、これにより、転換された該スペクトルの一部を急転させ、各サブバンドに対する再マップされた前記スペクトル測度における不連続部を生成することを含むこと、
前記少なくとも１つのチャンネルについて、各サブバンドに対して前記再マップされたスペクトル測度に対する自己回帰（ＡＲ）モデルを推定すること、
前記少なくとも１つのチャンネルについて各前記ＡＲモデルの係数を最小位相全ゼロサブバンド補正フィルタの係数にマップすること、
前記解析及び合成フィルタバンクの間でＰ個のベースバンドオーディオ信号を周波数補正するＰ個のデジタル全ゼロサブバンド補正フィルタを対応する係数から構成すること
によって行われることを特徴とするデバイス。
【請求項６】
前記室内解析モジュールは、前記第１の電気応答信号を該第１の電気応答が受信された時に重なるブロックに区分化して前記第１のプローブ信号の区分を用いて各ブロックに逆畳み込みを行い、かつ次のプローブ信号の送信の前の前記無音期間に前記遅延及び第１の室内応答を計算して記録するようになっていることを特徴とする請求項５に記載のデバイス。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、多チャンネルオーディオ再生デバイス及び方法に関し、より具体的には、多チャンネルラウドスピーカ構成を特徴付け、かつラウドスピーカ／室内の遅延、利得、及び周波数応答を補正するようになったデバイス及び方法に関する。
【背景技術】
【０００２】
ホームエンターテインメントシステムは、単純なステレオシステムからサラウンドサウンドシステム及び最近では３Ｄサウンドシステムのような多チャンネルオーディオシステム、及びビデオディスプレイを有するシステムに移行してきた。これらのホームエンターテインメントシステムは改善されたが、室内音響は、室内の面からの反響によってもたらされるサウンド歪み、及び／又はリスナーに対するラウドスピーカの不均一配置のような欠陥に依然として悩まされている。ホームエンターテインメントシステムは、広く家庭内に使用されているので、室内音響の改善は、自分達の好みのリスニング環境をより良く楽しむホームエンターテインメントシステムユーザに対する関心事である。
【０００３】
「サラウンドサウンド」は、スピーカの間に位置決めされたリスナーに模擬音源配置を与えるために複数のチャンネル及びスピーカを使用するサウンド再生システムを指すためにオーディオ工学において使用される用語である。リスナーを音源で「取り囲み」、それによってより面白い又は現実的なリスニング体験を発生させるために、サウンドは、スピーカのうちの１つ又はそれよりも多くを通じて異なる遅延を伴って異なる強度で再生することができる。従来のサラウンドサウンドシステムは、２次元スピーカ構成、例えば、前部、中央、背後、及び可能な場合は側部を含む。より最近の３Ｄサウンドシステムは、３次元スピーカ構成を含む。例えば、この構成は、高及び低前部、中央、背後、又は側部のスピーカを含むことができる。本明細書に使用する時に、多チャンネルスピーカ構成は、ステレオ、サラウンドサウンド、及び３Ｄサウンドシステムを含む。
【０００４】
多チャンネルサラウンドサウンドは、映画館及びホームシアターの用途に使用される。１つの一般的な構成では、ホームシアター内のリスナーは、従来のホームステレオシステムに使用される２つのスピーカの代わりに、５つのスピーカによって囲まれる。５つのスピーカのうちの３つは、室内の前部に配置され、残りの２つのサラウンドスピーカは、リスニング／ビューイング位置の後部又は側部（ＴＨＸ（登録商標）二重極）に位置付けられる。新しい構成は、サラウンドサウンド体験を模擬することができる複数のスピーカを含む「サウンドバー」を使用することである。現在使用されている様々なサラウンドサウンドフォーマットの中で、「Ｄｏｌｂｙ  Ｓｕｒｒｏｕｎｄ（登録商標）」は、映画館に向けて１９７０年代最初に開発されたオリジナルフォーマットである。１９９６年には「Ｄｏｌｂｙ  Ｄｉｇｉｔａｌ（登録商標）」が登場した。「Ｄｏｌｂｙ  Ｄｉｇｉｔａｌ（登録商標）」は、６つの離散オーディオチャンネルを有するデジタルフォーマットであり、４つのオーディオチャンネルを記録媒体上に格納される２つのチャンネルに組み合わせるマトリックスシステムを拠り所とする「Ｄｏｌｂｙ  Ｓｕｒｒｏｕｎｄ（登録商標）」のある一定の制約を解消する。「Ｄｏｌｂｙ  Ｄｉｇｉｔａｌ（登録商標）」は、５．１チャンネルフォーマットとも呼ばれ、映画サウンド記録に向けて数年前に世界的に採用された。現在使用されている別のフォーマットは、「Ｄｏｌｂｙ  Ｄｉｇｉｔａｌ（登録商標）」よりも高いオーディオ品質を提供する「ＤＴＳ  Ｄｉｇｉｔａｌ  Ｓｕｒｒｏｕｎｄ（登録商標）」（１，４１１，２００ビット対３８４，０００ビット毎秒）、並びに多くの異なるスピーカ構成、例えば、５．１、６．１、７．１、１１．２など、及びその変形、例えば、「フロント・ワイド」、「フロント・ハイト」、「センター・オーバーヘッド」、「サイド・ハイト」、又は「センター・ハイト」である。例えば、ＤＴＳ−ＨＤ（登録商標）は、Ｂｌｕ−Ｒａｙ（登録商標）ディスク上で７つの異なる７．１チャンネル構成をサポートする。
【０００５】
オーディオ／ビデオプリアンプ（又はＡ／Ｖコントローラ又はＡ／Ｖレシーバ）は、２チャンネルの「Ｄｏｌｂｙ  Ｓｕｒｒｏｕｎｄ（登録商標）」、「Ｄｏｌｂｙ  Ｄｉｇｉｔａｌ（登録商標）」、又は「ＤＴＳ  Ｄｉｇｉｔａｌ  Ｓｕｒｒｏｕｎｄ（登録商標）」、又はＤＴＳ−ＨＤ（登録商標）の信号をそれぞれの別々のチャンネルに復号するジョブを処理する。Ａ／Ｖプリアンプ出力は、それぞれ左、中央、右、左サラウンド、右サラウンド、及びサブウーファのチャンネルに対して６つのラインレベル信号を供給する。これらの別々の出力は、ホームシアターラウドスピーカシステムを駆動するために多チャンネル電力増幅器に供給されるか、又は一体型レシーバの場合と同じく内部で増幅される。
【０００６】
最良性能に関してＡ／Ｖプリアンプを手動で設定及び微調整するのは困難である可能性がある。ホームシアターシステムを製品マニュアルに従って接続した後に、ラウドスピーカ設定のためのプリアンプ又はレシーバを構成しなければならない。例えば、Ａ／Ｖプリアンプは、使用中の特定のサラウンドサウンドスピーカ構成を認識しなければならない。多くの場合に、Ａ／Ｖプリアンプは、デフォルト出力構成のみをサポートし、ユーザが、５．１又は７．１スピーカをそれらの場所に配置することができない場合には、ユーザは、単に運が無かったことになる。少数のハイエンドＡ／Ｖプリアンプは、複数の７．１構成をサポートし、ユーザに、室内のための適切な構成をメニューから選択させる。更に、オーディオチャンネルの各々（チャンネルの実際の数は、使用中の特定のサラウンドサウンドフォーマットによって決定される）のラウドネスは、ラウドスピーカからの音量に全体的な均衡を与えるために個々に設定しなければならない。この処理は、各スピーカからノイズの形態の「試験信号」を順次生成し、リスニング／ビューイング位置での各スピーカの音量を独立して調節することによって始まる。このタスクに推奨されるツールは、「音圧レベル（ＳＰＬ）」メーターである。これは、異なるラウドスピーカ感度、リスニング室内音響、及びラウドスピーカ配置に対する補償を与える。非対称リスニング空間及び／又は傾斜ビューイング区域、窓、アーチ通路、並びに傾斜天井のような他のファクタは、較正をより複雑にする可能性がある。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】米国特許第７，１５８，６４３号明細書
【特許文献２】米国特許公開第２００７，０１２１９５５号明細書
【非特許文献】
【０００８】
【非特許文献１】Ｔｈｏｍａｓ，Ｄ．Ｃ．（２００８年）著「音響強度及びエネルギ密度の理論及び推定（Ｔｈｅｏｒｙ  ａｎｄ  Ｅｓｔｉｍａｔｉｏｎ  ｏｆ  Ａｃｏｕｓｔｉｃ  Ｉｎｔｅｎｓｉｔｙ  ａｎｄ  Ｅｎｅｒｇｙ  Ｄｅｎｓｉｔｙ）」、科学修士論文、ブリガムヤング大学
【非特許文献２】Ｆａｈｙ，Ｆ．Ｊ．（１９９５年）著「サウンド強度（Ｓｏｕｎｄ  Ｉｎｔｅｎｓｉｔｙ）」、第２版、ロンドン、Ｅ＆ＦＮ  Ｓｐｏｎ
【非特許文献３】Ｃｏｂｏｓ，Ｍ．、Ｌｏｐｅｚ，Ｊ．Ｊ．、及びＭａｒｔｉ，Ａ．（２０１０年）著「四面体マイクロフォンアレイを用いた３Ｄ  ＤＯＡ推定における室内残響の効果に関して（Ｏｎ  ｔｈｅ  Ｅｆｆｅｃｔｓ  ｏｆ  Ｒｏｏｍ  Ｒｅｖｅｒｂｅｒａｔｉｏｎ  ｉｎ  ３Ｄ  ＤＯＡ  Ｅｓｔｉｍａｔｉｏｎ  Ｕｓｉｎｇ  Ｔｅｔｒａｈｅｄｒａｌ  Ｍｉｃｒｏｐｈｏｎｅ  Ａｒｒａｙ）」、第１２８回ＡＥＳ会議、ロンドン、ＵＫ、２０１０年５月２２〜２５日
【発明の概要】
【発明が解決しようとする課題】
【０００９】
従って、各オーディオチャンネルの周波数応答、振幅応答、及び時間応答を調節することによって多チャンネルサウンドシステムを自動的に較正するシステム及び処理を提供することが望ましいと考えられる。この処理は、リスナーを煩わすことなくサラウンドサウンドシステムの通常作動中に実施することができることがより望ましい。
【００１０】
「自動較正サラウンドシステム（Ａｕｔｏ−Ｃａｌｉｂｒａｔｉｎｇ  Ｓｕｒｒｏｕｎｄ  Ｓｙｓｔｅｍ）」という名称の米国特許第７，１５８，６４３号明細書は、サラウンドサウンドシステムの各チャンネルの周波数応答、振幅応答、及び時間応答の独立した自動較正及び自動調節を可能にする１つの手法を記載している。システムは、スピーカを通じて再生されてマイクロフォンによって記録される試験信号を生成する。システムプロセッサは、受信したサウンド信号を試験信号と相関付け、相関付けられた信号から白色応答を決定する。「室内音響補正デバイス（Ｒｏｏｍ  Ａｃｏｕｓｔｉｃｓ  Ｃｏｒｒｅｃｔｉｏｎ  Ｄｅｖｉｃｅ）」という名称の米国特許公開第２００７，０１２１９５５号明細書は、類似の手法を記載している。
【課題を解決するための手段】
【００１１】
以下は、本発明の一部の態様の基本的な理解をもたらすための本発明の要約である。この要約は、本発明の重要又は決定的な要素を明らかにすること、又は本発明の範囲を線引きするように意図したものではない。その唯一の目的は、後に提供するより詳細な説明及び決定的な特許請求の範囲の導入部として本発明のいくつかの概念を簡易形態で提供することである。
【００１２】
本発明は、多チャンネルラウドスピーカ構成を特徴付け、ラウドスピーカ／室内の遅延、利得、及び周波数応答を補正し、又はサブバンド領域補正フィルタを構成するようになったデバイス及び方法を提供する。
【００１３】
多チャンネルラウドスピーカ構成を特徴付けるための実施形態において、ブロードバンドプローブ信号は、複数のものがリスニング環境内で多チャンネル構成のラウドスピーカに結合されたＡ／Ｖプリアンプの各オーディオ出力に供給される。ラウドスピーカは、プローブ信号を無音期間によって分離された非重ね合わせ時間スロットで音波としてリスニング環境内に送信される音響応答に変換する。探査される各オーディオ出力に対して、音波が、音響応答をブロードバンド電気応答信号に変換する多マイクロフォンアレイによって受信される。次のプローブ信号の送信の前の無音期間において、プロセッサは、ブロードバンドプローブ信号を用いてブロードバンド電気応答信号に逆畳み込みを行ってラウドスピーカに対する各マイクロフォンでのブロードバンド室内応答を決定し、ラウドスピーカに対する各マイクロフォンでの遅延を計算してメモリに記録し、ラウドスピーカに対する遅延によってオフセットされた指定期間にわたって各マイクロフォンでのブロードバンド応答をメモリに記録し、オーディオ出力がラウドスピーカに結合されているか否かを決定する。オーディオ出力が結合されているか否かの決定は、各チャンネルに対する室内応答が処理されるまで延期することができる。プロセッサは、ブロードバンド電気応答信号をそれを受信した時に区分化し、区分化された信号を例えば分割ＦＦＴを用いて処理してブロードバンド室内応答を形成することができる。プロセッサは、区分化された信号から「ヒルベルト包絡線（ＨＥ）」を計算し、かつ継続的に更新することができる。ＨＥ内の顕著なピークを使用して遅延を計算し、かつオーディオ出力がラウドスピーカに結合されているか否かを決定することができる。
【００１４】
計算された遅延に基づいて、プロセッサは、各接続チャンネルに対してラウドスピーカに対する距離及び少なくとも第１の角度（例えば、方位角）を決定する。多マイクロフォンアレイが２つのマイクロフォンを含む場合には、プロセッサは、半平面において前部、いずれかの側部、又は後部に位置決めされたラウドスピーカに対する角度を分解することができる。多マイクロフォンアレイが３つのマイクロフォンを含む場合には、プロセッサは、３つのマイクロフォンによって定められる平面において前部、側部、及び後部に位置決めされたラウドスピーカに対する角度を分解することができる。多マイクロフォンアレイが４つ又はそれよりも多くのマイクロフォンを３Ｄ配置で含む場合には、プロセッサは、３次元空間に位置決めされたラウドスピーカに対する方位角と仰角の両方を分解することができる。結合されたラウドスピーカに対するこれらの距離及び角度を用いて、プロセッサは、特定の多チャンネル構成を自動的に選択し、かつリスニング環境内の各ラウドスピーカの位置を計算する。
【００１５】
ラウドスピーカ／室内周波数応答を補正するための実施形態において、ブロードバンドプローブ信号及び可能な場合はプリ−エンファシスプローブ信号は、少なくとも複数のものがリスニング環境内で多チャンネル構成のラウドスピーカに結合されたＡ／Ｖプリアンプの各オーディオ出力に供給される。ラウドスピーカは、プローブ信号を無音期間によって分離された非重ね合わせ時間スロットで音波としてリスニング環境内に送信される音響応答に変換する。探査される各オーディオ出力に対して、音波が、音響応答を電気応答信号に変換する多マイクロフォンアレイによって受信される。プロセッサは、ブロードバンドプローブ信号を用いて電気応答信号に逆畳み込みを行い、ラウドスピーカに対する各マイクロフォンでの室内応答を決定する。
【００１６】
プロセッサは、室内応答から室内エネルギ測度を計算する。プロセッサは、遮断周波数よりも高い周波数に対して、室内エネルギ測度の第１の部分を音圧の関数として計算し、遮断周波数よりも低い周波数に対して、室内エネルギ測度の第２の部分を音圧及び音速の関数として計算する。音速は、マイクロフォンアレイにわたる音圧の勾配から得られる。ブロードバンド及びプリ−エンファシスプローブ信号の両方を含む二重プローブ信号が利用される場合には、音圧だけに基づくエネルギ測度の高周波部分は、ブロードバンド室内応答から抽出され、音圧と音速の両方に基づくエネルギ測度の低周波部分は、プリ−エンファシス室内応答から抽出される。二重プローブ信号を使用して、音速成分なしで室内エネルギ測度を計算することができ、この場合に、プリ−エンファシスプローブ信号は、ノイズ整形に使用される。プロセッサは、エネルギ測度の第１及び第２の部分を混合して指定音響バンドにわたる室内エネルギ測度を与える。
【００１７】
知覚的により適切な測定値を得るために、室内応答又は室内エネルギ測度は、最も低い周波数で実質的に時間応答全体を取り込み、かつ最も高い周波数で基本的に時間応答の直接経路プラス数ミリ秒のみを取り込むように漸進的に平滑化することができる。プロセッサは、プロセッサ内のデジタル補正フィルタを構成するのに使用されるフィルタ係数を室内エネルギ測度から計算する。プロセッサは、ユーザ定義のものか又はチャンネルエネルギ測度の平滑化バージョンであるチャンネルターゲット曲線のためのフィルタ係数を計算することができ、次に、ユーザ定義のものか又はチャンネルターゲット曲線の平均とすることができる共通ターゲット曲線に対してフィルタ係数を調節することができる。プロセッサは、オーディオ信号をリスニング環境内への再生に向けてラウドスピーカまで対応するデジタル補正フィルタを通過させる。
【００１８】
多チャンネルオーディオシステムのためのサブバンド補正フィルタを発生させるための実施形態において、Ｐが整数である時に、オーディオ信号をＰ個のサブバンドに対してベースバンドまでダウンサンプリングするＰバンドオーバーサンプリング解析フィルタバンク、及びオーディオ信号を再構成するためにＰ個のサブバンドをアップサンプリングするＰバンドオーバーサンプリング合成フィルタバンクが、Ａ／Ｖプリアンプ内のプロセッサに与えられる。スペクトル測度が、各チャンネルに対して与えられる。プロセッサは、各スペクトル測度をチャンネルターゲット曲線と組み合わせてチャンネル毎の集合スペクトル測度を与える。各チャンネルに対して、プロセッサは、集合スペクトル測度のうちで異なるサブバンドに対応する部分を抽出し、スペクトル測度の抽出部分をベースバンドに再マップして解析フィルタバンクのダウンサンプリングを模擬する。プロセッサは、各サブバンドに対する再マップしたスペクトル測度に対する自己回帰（ＡＲ）モデルを計算し、かつ各ＡＲモデルの係数を最小位相全ゼロサブバンド補正フィルタの係数にマップする。プロセッサは、再マップしたスペクトル測度の逆ＦＦＴとして自己相関シーケンスを計算し、かつＡＲモデルを計算するために自己相関シーケンスにレビンソン−ダービンアルゴリズムを適用することにより、ＡＲモデルを計算することができる。レビンソン−ダービンアルゴリズムは、補正フィルタの次数を選択するのに使用することができるサブバンドに対する残存電力推定値を生成する。プロセッサは、解析及び合成フィルタバンク間でＰ個のベースバンドオーディオ信号を周波数補正するＰ個のデジタル全ゼロサブバンド補正フィルタを対応する係数から構成する。プロセッサは、ユーザ定義のものであるか又はチャンネルエネルギ測度の平滑化バージョンであるチャンネルターゲット曲線のためのフィルタ係数を計算することができ、次に、チャンネルターゲット曲線の平均とすることができる共通ターゲット曲線に対してフィルタ係数を調節することができる。
【００１９】
以下に続く好ましい実施形態の詳細説明を添付図面と合わせて理解することにより、本発明のこれら及び他の特徴及び利点は当業者に明らかであろう。
【図面の簡単な説明】
【００２０】
【図１ａ】多チャンネルオーディオ再生システム及びリスニング環境の解析モードにおける実施形態のブロック図である。
【図１ｂ】四面体マイクロフォンの実施形態の図である。
【図２】多チャンネルオーディオ再生システム及びリスニング環境の再生モードにおける実施形態のブロック図である。
【図３】解析モードにおいて決定されたラウドスピーカ／室内周波数応答のずれを補正するようになったサブバンド通過フィルタバンクの再生モードにおける実施形態のブロック図である。
【図４】解析モードの実施形態の流れ図である。
【図５ａ】全域通過プローブ信号に対する時間シーケンスを示す図である。
【図５ｂ】全域通過プローブ信号に対する周波数シーケンスを示す図である。
【図５ｃ】全域通過プローブ信号に対する自己相関シーケンスを示す図である。
【図５ｄ】全域通過プローブ信号に対する自己相関シーケンスを示す図である。
【図６ａ】プリ−エンファシスプローブ信号の時間シーケンスを示す図である。
【図６ｂ】プリ−エンファシスプローブ信号のマグニチュードスペクトルを示す図である。
【図７】同じ周波数領域信号から全域通過プローブ信号とプリ−エンファシスプローブ信号とを生成するための実施形態の流れ図である。
【図８】取得のためのプローブ信号の送信をスケジューリングするための実施形態の図である。
【図９】室内応答及び遅延を与えるためのプローブ信号の実時間取得処理における実施形態のブロック図である。
【図１０】補正フィルタを与えるための室内応答の後処理に対する実施形態の流れ図である。
【図１１】ブロードバンドプローブ信号のスペクトル測度とプリ−エンファシスプローブ信号のスペクトル測度とで混合された室内スペクトル測度の実施形態の図である。
【図１２】異なるプローブ信号及びマイクロフォン組合せに対してエネルギ測度を計算するための実施形態の流れ図である。
【図１３】周波数補正フィルタを計算するためにエネルギ測度を処理するための実施形態の流れ図である。
【図１４ａ】解析フィルタバンクのダウンサンプリングを模擬するためのエネルギ測度の抽出及びベースバンドへの再マッピングにおける実施形態を示す図である。
【図１４ｂ】解析フィルタバンクのダウンサンプリングを模擬するためのエネルギ測度の抽出及びベースバンドへの再マッピングにおける実施形態を示す図である。
【図１４ｃ】解析フィルタバンクのダウンサンプリングを模擬するためのエネルギ測度の抽出及びベースバンドへの再マッピングにおける実施形態を示す図である。
【発明を実施するための形態】
【００２１】
本発明は、多チャンネルラウドスピーカ構成を特徴付け、ラウドスピーカ／室内の遅延、利得、及び周波数応答を補正し、又はサブバンド領域補正フィルタを構成するようになったデバイス及び方法を提供する。様々なデバイス及び方法は、空間内でラウドスピーカを自動的に位置付けて、オーディオチャンネルが接続されているか否かを決定し、特定の多チャンネルラウドスピーカ構成を選択し、かつリスニング環境内で各ラウドスピーカを位置決めするようになっている。様々なデバイス及び方法は、低周波数で音圧と音速の両方を取り込み、かつ広いリスニング区域にわたって正確である知覚的に適切なエネルギ測度を抽出するようになっている。エネルギ測度は、リスニング環境内の単一の場所に配置されてデジタル補正フィルタを構成するのに使用される密に離間された非共在多マイクロフォンアレイを使用することによって収集された室内応答から導出される。様々なデバイス及び方法は、入力多チャンネルオーディオ信号の周波数応答を例えば室内応答及びラウドスピーカ応答によってもたらされるターゲット応答からのずれに対して補正するためのサブバンド補正フィルタを構成するようになっている。解析フィルタバンクのダウンサンプリングを模擬するために、スペクトル測度（室内スペクトル／エネルギ測度のような）が区分化され、ベースバンドに再マップされる。各サブバンドに対してＡＲモデルが独立して計算され、これらのモデルの係数は、全ゼロ最小位相フィルタにマップされる。注意点として、解析フィルタの形状は、再マッピングには含まれない。サブバンド通過フィルタの実施は、ＭＩＰＳ要件と、メモリ要件と、処理遅延とを均衡調整するように構成することができ、他のオーディオ処理に向けて解析／合成フィルタバンクアーキテクチャが既に存在する場合にはその上に搭載することができる。
【００２２】
多チャンネルオーディオ解析及び再生システム
ここで図面を参照すると、図１ａ〜図１ｂ、図２、及び図３は、リスニング環境１４内の多チャンネルスピーカ構成１２を探査及び解析して多チャンネルスピーカ構成を自動的に選択して室内でスピーカを位置決めし、広いリスニング区域にわたって知覚的に適切なスペクトル（例えば、エネルギ）測度を抽出し、かつ周波数補正フィルタを構成するためのかつ室内補正（遅延、利得、及び周波数）を用いた多チャンネルオーディオ信号１６の再生のための多チャンネルオーディオシステム１０の実施形態を示している。多チャンネルオーディオ信号１６は、ケーブル給送又は衛星給送によって供給することができ、又はＤＶＤ又はＢｌｕ−Ｒａｙ（登録商標）ディスクのようなストレージ媒体から読み取ることができる。オーディオ信号１６は、テレビジョン１８に供給されるビデオ信号と対にすることができる。代替的に、オーディオ信号１６は、ビデオ信号を伴わない音楽信号とすることができる。
【００２３】
多チャンネルオーディオシステム１０は、多チャンネルオーディオ信号１６を供給するためのケーブルレシーバ又は衛星レシーバ又はＤＶＤプレーヤ又はＢｌｕ−Ｒａｙ（登録商標）プレーヤのようなオーディオソース２０と、多チャンネルオーディオ信号をオーディオ出力２４における個別オーディオチャンネルに復号するＡ／Ｖプリアンプ２２と、それぞれのオーディオ出力２４に結合され、Ａ／Ｖプリアンプによって供給される電気信号をリスニング環境１４内に音波２８として送信される音響応答に変換する複数のラウドスピーカ２６（電気−音響変換器）とを含む。オーディオ出力２４は、ラウドスピーカに有線接続した端子、又はラウドスピーカに無線結合された無線出力とすることができる。オーディオ出力がラウドスピーカに結合される場合に、対応するオーディオチャンネルが接続されると言われる。ラウドスピーカは、サラウンドサウンド体験を模倣するようになった複数のスピーカを各々が含む離散２Ｄレイアウト又は離散３Ｄレイアウト又はサウンドバーで配置された個別スピーカとすることができる。システムは、１つ又はそれよりも多くのマイクロフォン３０とマイクロフォン送信ボックス３２とを含むマイクロフォンアセンブリも含む。マイクロフォン（音響−電気変換器）は、ラウドスピーカに供給されるプローブ信号に対する音波を受信し、かつ音響応答を電気信号に変換する。送信ボックス３２は、電気信号をＡ／Ｖプリアンプのオーディオ入力３４のうちの１つ又はそれよりも多くに有線又は無線の接続を通じて供給する。
【００２４】
Ａ／Ｖプリアンプ２２は、典型的に独自のプロセッサメモリが備えられた汎用「コンピュータ処理ユニット（ＣＰＵ）」又は専用「デジタル信号プロセッサ（ＤＳＰ）」チップのような１つ又はそれよりも多くのプロセッサ３６と、システムメモリ３８と、オーディオ出力２４に接続したデジタル／アナログコンバータ及び増幅器４０とを含む。一部のシステム構成では、Ｄ／Ａコンバータ及び／又は増幅器は、別々のデバイスとすることができる。例えば、Ａ／Ｖプリアンプは、補正されたデジタル信号をＤ／Ａコンバータに出力することができ、Ｄ／Ａコンバータは、アナログ信号を電力増幅器に出力する。解析及び再生の作動モードを実施するために、メモリ、プロセッサ、又はシステム内にコンピュータプログラム命令の様々な「モジュール」が格納され、１つ又はそれよりも多くのプロセッサ３６によって実行される。
【００２５】
Ａ／Ｖプリアンプ２２は、入力マイクロフォン信号を受信するために１つ又はそれよりも多くのオーディオ入力３４に接続されてプロセッサ３６に別々のマイクロフォンチャンネルを供給する入力レシーバ４２も含む。マイクロフォン送信ボックス３２と入力レシーバ４２とは対応する対である。例えば、送信ボックス３２は、マイクロフォンアナログプリアンプと、Ａ／Ｄコンバータ及びＴＤＭ（時間領域マルチプレクサ）又はＡ／Ｄコンバータと、圧縮器と、ＵＳＢトランスミッタとを含むことができ、対応する入力レシーバ４２は、アナログプリアンプ及びＡ／Ｄコンバータと、ＳＰＤＩＦレシーバ及びＴＤＭデマルチプレクサ又はＵＳＢレシーバと、解凍器とを含むことができる。Ａ／Ｖプリアンプは、各マイクロフォン信号に対するオーディオ入力３４を含むことができる。代替的に、複数のマイクロフォン信号を単一の信号に多重化して、単一のオーディオ入力３４に供給することができる。
【００２６】
作動の解析モード（図４に提供している）をサポートするために、Ａ／Ｖプリアンプには、プローブ発生及び送信スケジューリングモジュール４４と室内解析モジュール４６とが備えられる。図５ａ〜図５ｄ、図６ａ〜図６ｂ、図７、及び図８に詳描するように、モジュール４４は、ブロードバンドプローブ信号、及び可能な場合はそれと対になったプリ−エンファシスプローブ信号を発生し、これらのプローブ信号をＡ／Ｄコンバータ及び増幅器４０を通じて各オーディオ出力２４にスケジュールに従って無音期間によって分離された非重ね合わせ時間スロットで送信する。各オーディオ出力２４は、出力がラウドスピーカに結合されているか否かに関して探査される。モジュール４４は、１つ又は複数のプローブ信号及び送信スケジュールを室内解析モジュール４６に供給する。図９から図１４に詳描するように、モジュール４６は、送信スケジュールに従って多チャンネルスピーカ構成を自動的に選択し、室内にスピーカを位置決めし、広いリスニング区域にわたって知覚的に適切なスペクトル（エネルギ）測度を抽出し、周波数補正フィルタ（サブバンド周波数補正フィルタなど）を構成するようにマイクロフォン信号及びプローブ信号を処理する。モジュール４６は、ラウドスピーカ構成及びスピーカ位置とフィルタ係数とをシステムメモリ３８に格納する。
【００２７】
マイクロフォン３０の個数及びレイアウトは、多チャンネルラウドスピーカ構成を選択してラウドスピーカを位置決めし、広いリスニング区域にわたって有効である知覚的に適切なエネルギ測度を抽出する解析モジュールの機能に影響を及ぼす。これらの機能をサポートするために、マイクロフォンレイアウトは、ラウドスピーカを２次元又は３次元で「局在化」させ、かつ音速を計算するためのある一定量の多様性を与える。一般的に、これらのマイクロフォンは非共在のものであり、固定の間隔を有する。例えば、単一のマイクロフォンは、ラウドスピーカまでの距離のみを推定する段階をサポートする。１対のマイクロフォンは、ラウドスピーカまでの距離及び平面の半域（前部、後部、又はいずれかの側部）内での方位角のような角度を推定する段階と、単一の方向の音速を推定する段階とをサポートする。３つのマイクロフォンは、ラウドスピーカまでの距離及び平面全域（前部、後部、及び両側部）内での方位角を推定する段階と、３次元空間における音速を推定する段階とをサポートする。３次元ボール上に位置決めされた４つ又はそれよりも多くのマイクロフォンは、ラウドスピーカまでの距離、並びに完全な３次元空間内での方位角及び仰角を推定する段階と、３次元空間における音速を推定する段階とをサポートする。
【００２８】
四面体マイクロフォンアレイの場合の特定的に選択された座標系における多マイクロフォンアレイ４８の実施形態を図１ｂに示している。４つのマイクロフォン３０は、四面体物体（「ボール」）４９の頂点に配置される。全てのマイクロフォンは全方位性のものであり、すなわち、マイクロフォン信号は、様々な場所における圧力測定値を表すと仮定する。マイクロフォン１、２、及び３は、ｘ，ｙ平面に位置し、マイクロフォン１は、座標系の原点、及びマイクロフォン２及び３は、ｘ軸から等距離に位置する。マイクロフォン４は、ｘ，ｙ平面から外れて位置する。マイクロフォンの各々の間の距離は等しく、ｄで表している。到来方向（ＤＯＡ）は、音波到来方向（補記Ａの位置処理に使用されることになる）を示している。マイクロフォンの間隔「ｄ」は、最大で５００Ｈｚから１ｋＨｚまでの音速を正確に計算するために小さい間隔を必要とすること、及びラウドスピーカを正確に位置決めするために大きい間隔を必要とすることの兼ね合いを表している。約８．５ｃｍから９ｃｍの間隔が、両方の要件を満たす。
【００２９】
作動の再生モードをサポートするために、Ａ／Ｖプリアンプには、入力レシーバ／復号器モジュール５２とオーディオ再生モジュール５４とが備えられる。入力レシーバ／復号器モジュール５２は、多チャンネルオーディオ信号１６を別々のオーディオチャンネルに復号する。例えば、多チャンネルオーディオ信号１６は、標準の２チャンネルフォーマットで配信することができる。モジュール５２は、２チャンネルの「Ｄｏｌｂｙ  Ｓｕｒｒｏｕｎｄ」信号、「Ｄｏｌｂｙ  Ｄｉｇｉｔａｌ」信号、又は「ＤＴＳ  Ｄｉｇｉｔａｌ  Ｓｕｒｒｏｕｎｄ（登録商標）」信号又はＤＴＳ−ＨＤＲ（登録商標）信号をそれぞれの別々のオーディオチャンネルに復号するジョブを処理する。モジュール５４は、汎用フォーマット変換、並びにラウドスピーカ／室内の較正及び補正を実施するように各オーディオチャンネルを処理する。例えば、モジュール５４は、アップミキシング又はダウンミキシング、スピーカの再マッピング又は仮想化を実施し、遅延、利得、又は極性の補償を適用し、バス管理を実施し、かつ室内周波数補正を実施することができる。モジュール５４は、１つ又はそれよりも多くのデジタル周波数補正フィルタを各オーディオチャンネルに向けて構成するのに、解析モードによって生成され、システムメモリ３８に格納された周波数補正パラメータ（例えば、遅延及び利得の調節値及びフィルタ係数）を使用することができる。周波数補正フィルタは、時間領域、周波数領域、又はサブバンド領域に実施することができる。各オーディオチャンネルは、各々の周波数補正フィルタを通過し、リスニング環境内に音波として送信される音響応答を生成するようにラウドスピーカを駆動するアナログオーディオ信号に変換される。
【００３０】
サブバンド領域に実施されたデジタル周波数補正フィルタ５６の実施形態を図３に示している。フィルタ５６は、Ｐバンド複素非臨界サンプリング解析フィルタバンク５８と、Ｐ個のサブバンドに対するＰ個の最小位相ＦＩＲ（有限インパルス応答）フィルタ６２と、Ｐバンド複素非臨界サンプリング合成フィルタバンク６４とを含み、ここでＰは整数である。図示のように、室内周波数補正フィルタ６０は、汎用アップミキシング／ダウンミキシング／スピーカ再マッピング／仮想化の機能６６をサブバンド領域に実施する「ＤＴＳ  ＮＥＯ−ＸＴＭ」のような既存のフィルタアーキテクチャに追加される。サブバンド単位の室内周波数補正における計算の大部分は、解析フィルタバンク及び合成フィルタバンクの実施にある。室内補正の追加によってＮＥＯ−ＸＴＭのような既存のサブバンドアーキテクチャに課せられる処理要件の増分は最小である。
【００３１】
周波数補正は、サブバンド領域で、最初にオーディオ信号（例えば、入力ＰＣＭサンプル）をオーバーサンプリング解析フィルタバンク５８に通し、次に、適宜異なる長さの最小位相ＦＩＲ補正フィルタ６２を各バンド内で独立して適用し、最後に合成フィルタバンク６４を適用して、周波数補正された出力ＰＣＭオーディオ信号を作成することによって実施される。周波数補正フィルタは、最小位相のものであるように設計されるので、サブバンド信号は、異なる長さのフィルタを通過した後でさえも、依然としてバンド間で時間整合されている。その結果、この周波数補正手法によって導入される遅延は、解析フィルタバンク及び合成フィルタバンクのチェーン内の遅延だけによって決定される。６４バンドオーバーサンプリング複素フィルタバンクを用いた特定の実施では、この遅延は２０ミリ秒よりも小さい。
【００３２】
取得、室内応答処理、及びフィルタ構成
作動の解析モードの実施形態に対する高レベルの流れ図を図４に示している。一般的に、解析モジュールは、ブロードバンドプローブ信号、及び可能な場合はプリ−エンファシスプローブ信号を発生し、これらのプローブ信号をスケジュールに従ってラウドスピーカを通じてリスニング環境内に音波として送信し、マイクロフォンアレイにおいて検出される音響応答を記録する。これらのモジュールは、各マイクロフォンでの各ラウドスピーカに対する遅延及び室内応答、並びに各プローブ信号を計算する。この処理は、次のプローブ信号の送信の前に「実時間」で行うか、又は全てのプローブ信号を送信し終わり、マイクロフォン信号を記録し終わった後にオフラインで行うことができる。モジュールは、各ラウドスピーカに対してスペクトル（例えば、エネルギ）測度を計算し、このスペクトル測度を用いて周波数補正フィルタ及び利得調節値を計算するように室内応答を処理する。ここでもまた、この処理は、次のプローブ信号の送信の前の無音期間中に行うか、又はオフラインで行うことができる。取得及び室内応答処理が実時間又はオフラインのいずれで行われるかは、１００万命令毎秒（ＭＩＰＳ）で測定される計算と、メモリと、全体の取得時間との兼ね合いであり、特定のＡ／Ｖプリアンプのリソース及び要件に依存する。モジュールは、計算された各ラウドスピーカまでの遅延を用いて、接続した各チャンネルにおけるラウドスピーカに対する距離及び少なくとも方位角を決定し、これらの情報を用いて特定の多チャンネル構成を自動的に選択し、リスニング環境内の各ラウドスピーカに対する位置を計算する。
【００３３】
解析モードは、システムパラメータ及び解析モジュールパラメータを初期化する段階で始まる（段階７０）。システムパラメータは、利用可能なチャンネルの個数（ＮｕｍＣｈ）と、マイクロフォンの個数（ＮｕｍＭｉｃｓ）と、マイクロフォン感度、出力レベルなどに基づく出力音量設定とを含むことができる。解析モジュールパラメータは、１つ又は複数のプローブ信号Ｓ（ブロードバンド）及びＰｅＳ（プリ−エンファシス処理されたもの）と、利用可能チャンネルの各々に信号を送信するためのスケジュールとを含む。プローブ信号は、システムメモリに格納するか又は解析が開始される時に発生させることができる。スケジュールは、システムメモリに格納するか又は解析が開始される時に発生させることができる。スケジュールは、１つ又はそれよりも多くのプローブ信号を各プローブ信号が無音期間によって分離された非重ね合わせ時間スロットでスピーカによってリスニング環境内に音波として送信されるようにオーディオ出力に供給する。無音期間の幅は、処理のうちのいずれかが、次のプローブ信号の送信の前に実施されているか否かに少なくとも部分的に依存することになる。
【００３４】
第１のプローブ信号Ｓは、指定された音響バンドにわたって実質的に一定のマグニチュードスペクトルによって特徴付けられるブロードバンドシーケンスである。音響バンド内での一定のマグニチュードスペクトルからのずれは、ＳＮ比（ＳＮＲ）を犠牲にし、それによって室内及び補正フィルタの特徴付けが影響を受ける。システム仕様は、音響バンドにわたる一定量からの最大ｄＢずれを指定することができる。第２のプローブ信号ＰｅＳは、指定された音響バンドの一部分にわたって増幅されたマグニチュードスペクトルを与えるベースバンドシーケンスに適用されるプリ−エンファシス関数によって特徴付けられるプリ−エンファシスシーケンスである。プリ−エンファシスシーケンスは、ブロードバンドシーケンスから導出することができる。一般的に、第２のプローブ信号は、指定された音響バンドと部分的又は完全に重なる可能性がある一定のターゲットバンド内でのノイズの整形又は減衰に有利であるとすることができる。特定の用途では、プリ−エンファシス関数のマグニチュードは、指定された音響バンドの低周波数領域と重なるターゲットバンド内の周波数に反比例する。多マイクロフォンアレイとの組合せに使用される場合には、この二重プローブ信号は、ノイズの存在下でよりロバストな音速計算を与える。
【００３５】
プリアンプのプローブ発生及び送信スケジューリングモジュールは、スケジュールに従ってプローブ信号の送信及びマイクロフォン信号Ｐ及びＰｅＰの取り込みを開始する（段階７２）。プローブ信号（Ｓ及びＰｅＳ）及び取り込みマイクロフォン信号（Ｐ及びＰｅＰ）は、室内解析モジュールに供給され、室内応答の取得が実施される（段階７４）。この取得は、室内応答を時間領域室内インパルス応答（ＲＩＲ）又は周波数領域室内周波数応答（ＲＦＲ）のいずれかで出力し、更に各ラウドスピーカに対する各取り込みマイクロフォン信号における遅延を出力する。
【００３６】
一般的に、取得処理は、室内応答を抽出するためのプローブ信号を用いたマイクロフォン信号の逆畳み込みを含む。逆畳み込みは、ブロードバンドプローブ信号を用いてブロードバンドマイクロフォン信号に施される。プリ−エンファシスマイクロフォン信号には、プリ−エンファシスマイクロフォン信号を使用するか、又はブロードバンドプローブ信号とすることができるそのベースバンドシーケンスを用いた逆畳み込みを行うことができる。プリ−エンファシスマイクロフォン信号に、そのベースバンドシーケンスを用いて逆畳み込みを行うことにより、プリ−エンファシス関数が室内応答の上に重ね合わせられる。
【００３７】
逆畳み込みは、マイクロフォン信号のＦＦＴ（高速フーリエ変換）を計算し、プローブ信号のＦＦＴを計算し、マイクロフォン周波数応答をプローブ周波数応答で割算して室内周波数応答（ＲＦＲ）を形成することによって実施することができる。ＲＩＲは、ＲＦＲの逆ＦＦＴを計算することによって与えられる。逆畳み込みは、全体のマイクロフォン信号を記録し、全体のマイクロフォン信号及びプローブ信号に対して単一のＦＦＴを計算することによって「オフライン」に実施することができる。この逆畳み込みは、プローブ信号の間の無音期間内に行うことができるが、この計算に対応するために、無音期間の持続時間を増大する必要がある場合がある。代替的に、いずれかの処理が始まる前に、全てのチャンネルにおけるマイクロフォン信号を記録してメモリに格納することができる。逆畳み込みは、マイクロフォン信号が取り込まれる時に、それをブロックに区分化し、この区分化に基づいてマイクロフォン信号及びプローブ信号に対してＦＦＴを計算することによって「実時間」に実施することができる（図９を参照されたい）。「実時間」手法は、メモリ要件を軽減するが、取得時間を延長させる傾向を有する。
【００３８】
取得は、各ラウドスピーカに対する取り込みマイクロフォン信号の各々に対する遅延を計算する段階を有する。遅延は、プローブ信号とマイクロフォン信号とで信号の相互相関、相互スペクトル位相、又はヒルベルト包絡線（ＨＥ）のような解析包絡線を用いて計算することができる。遅延は、例えば、ＨＥ内の顕著なピーク（例えば、定められた閾値よりも大きい最大ピーク）の位置に対応する場合がある。時間領域シーケンスを発生させるＨＥのような技術は、サンプリング期間精度の数分の１のより細かい時間尺度でピークの新しい場所を計算するために、ピークの前後で内挿することができる。サンプリング期間は、受信マイクロフォン信号がサンプリングされる間隔であり、当業技術で一般に知られているように、サンプリングされる最大周波数の逆数の２分の１よりも小さいか又はそれに等しいように選ばなければならない。
【００３９】
取得は、オーディオ出力が実際にラウドスピーカに結合されるか否かを決定する段階も伴う。端子が結合されていない場合には、マイクロフォンは、いずれかの環境信号を依然として拾って記録することになるが、相互相関／相互スペクトル位相／解析包絡線は、ラウドスピーカ接続を示す顕著なピークを示すことにはならない。取得モジュールは、最大ピークを記録し、それを閾値と比較する。ピークがそのピークよりも大きい場合には、ＳｐｅａｋｅｒＡｃｔｉｖｉｔｙＭａｓｋ［ｎｃｈ］が真値に設定され、オーディオチャンネルが接続されていると見なされる。この決定は、無音期間中又はオフラインで行うことができる。
【００４０】
接続した各オーディオチャンネルに対して、解析モジュールは、各マイクロフォンでの各ラウドスピーカからの室内応答（ＲＩＲ又はＲＦＲのいずれか）及び遅延を処理し、各ラウドスピーカに対する室内スペクトル測度を出力する（段階７６）。この室内応答処理は、次のプローブ信号の送信の前の無音期間中に、又は全ての探査及び取得が終了した後にオフラインに実施することができる。最も単純な場合には、室内スペクトル測度は、可能な場合は複数のマイクロフォンにわたって平均化され、更に可能な場合は高周波数でブロードバンドＲＦＲを用い、低周波数でプリ−エンファシスＲＦＲを使用するように混合された単一のマイクロフォンに対するＲＦＲを含むことができる。室内応答の更に別の処理は、知覚的により適切なスペクトル応答及びより広いリスニング区域にわたって有効なものをもたらすことができる。
【００４１】
標準の室内（リスニング環境）には、通常の利得／距離の問題以外に、室内補正を如何に測定し、計算し、適用することができるかに影響を及ぼすいくつかの音響的な問題が存在する。これらの問題を理解するためには、知覚問題を考慮すべきである。特に、人間の聴覚において「先行音効果」としても公知の「最初の到来音」は、画像及び音質の実際の知覚にある一定の役割を果たす。無響室を除くあらゆるリスニング環境では、「直接」音質、すなわち、音源の実際の知覚音質は、最初の到来（スピーカ／機器から直接の）サウンドと、最初のいくつかの反響音とによって影響を受ける。この直接音質を聴き取った後に、リスナーは、この音質を室内で反響された後のサウンドのものと比較する。耳の全空間電力応答に対する直接音への頭部伝達関数（ＨＲＴＦ）の影響の比較は、人間が把握しており、使用することを学習するものであるので、上述の比較は、取りわけ、前部／後部曖昧性除去のような問題において役立つ。注意する点は、直接信号が重み付き間接信号よりも多い高周波数を有する場合には、この信号は、一般的に「前方」のものとして聞こえ、それに対して高周波数を欠く直接信号は、リスナーの背後に位置することになる点である。この効果は、約２ｋＨｚから上で最も強い。聴覚システムの性質に起因して、下側遮断周波数から約５００Ｈｚまでの信号が１つの方法によって位置し、それよりも大きい信号が別の方法によって位置する。
【００４２】
最初の到来音に起因する高周波数知覚効果に加えて、物理的な音響が室内補償において大きい効果を有する。殆どのラウドスピーカは、最初の到来音に関して理想的なものに確実に近づく場合であっても、全体として平坦な電力放射曲線を持たない。これは、高周波数では、リスニング環境が、低周波数で駆動されることになる場合よりも少ないエネルギで駆動されることになることを意味する。これだけを取ると、それは、補償の計算において長期のエネルギ平均を用いた場合に、直接信号に望ましくないプリ−エンファシスを適用してしまうことを意味することになる。一般的に、より高い周波数では、壁、家具、人間などがより多くのエネルギを吸収することになり、それによって室内のエネルギ累積量（すなわち、Ｔ６０）が減少し、長期測定値が、直接音質に対するより一層高い誤認関係を有することがもたらされるので、残念ながら、一般的な室内音響によって状況は悪化する。
【００４３】
その結果、本出願人の手法は、実際の蝸牛機構によって決定される直接音に関する測定を低周波数で長い測定期間を有し（蝸牛フィルタの長いインパルス応答に起因して）、高周波数で短い測定期間を有するものにする。低周波数から高周波数への遷移は滑らかに変更される。この期間は、ｔ＝２／ＥＲＢバンド幅という規則によって近似することができ、この場合ＥＲＢは、聴覚システム内の他のファクタが、時間を更に短縮すべきではないことを示す時間である数ミリ秒の下側の限界に「ｔ」が達するまでの同一矩形バンド幅である。この「段階的平滑化」は、室内インパルス応答又は室内スペクトル測度に対して実施することができる。
【００４４】
低い周波数、すなわち、長い波長では、サウンドエネルギは、様々な場所にわたって、音圧だけ又は速度のいずれかの軸だけと比較して僅かしか変化しない。非共在多マイクロフォンアレイからの測定値を用いて、モジュールは、低周波数で、音圧だけでなく、好ましくは、全方向の音速も考慮する全体エネルギ測度を計算する。それを行うことにより、モジュールは、室内の１つの点から低周波数における実際の累積エネルギを取り込む。それによってゼロ圧力が体積速度の最大値と一致することになることで測定点における圧力が過剰累積量を指定しない場合であっても、Ａ／Ｖプリアンプが、過剰累積量が存在する周波数におけるエネルギを室内に放射することを回避することが好都合に可能になる。多マイクロフォンアレイとの組合せに使用される場合には、二重プローブ信号は、ノイズの存在下でよりロバストな室内応答をもたらす。
【００４５】
解析モジュールは、室内スペクトル（例えば、エネルギ）測度を用いて、接続した各オーディオチャンネルに対して周波数補正フィルタ及び利得調節値を計算し、これらのパラメータをシステムメモリに格納する（段階７８）。ラウドスピーカ／室内周波数補正を適用するのに、時間領域フィルタ（例えば、ＦＩＲ又はＩＩＲ）、周波数領域フィルタ（例えば、重ね合わせ加算、重ね合わせ保留によって実施されたＦＩＲ）、及びサブバンド領域フィルタを使用することができる。非常に低い周波数における室内補正は、数百ミリ秒の持続時間に容易に達することができるインパルス応答を有する補正フィルタを必要とする。サイクル毎に必要とされる作動に関して、これらのフィルタを実施する最も効率的な手法は、重ね合わせ保留法又は重ね合わせ加算法を用いた周波数領域におけるものになる。必要とされるＦＦＴの大きいサイズに起因して、いくつかの庭用電化製品用途では、固有遅延及びメモリ要件が法外のものである可能性がある。分割ＦＦＴ手法が使用される場合には、遅延は、サイクル毎の作動回数の増加の代償として短縮することができる。しかし、本方法は、依然として高いメモリ要件を有する。処理がサブバンド領域に実施される場合には、サイクル毎に必要とされる作動回数と、メモリ要件と、処理遅延との間の妥協点を微調整することができる。サブバンド領域での周波数補正は、特に、非常に少ないサブバンド（非常に少ない低周波数バンドのみを有する室内補正の場合のような）内のフィルタが、全ての他のサブバンド内のフィルタよりもかなり高い次数を有する場合に、異なる周波数領域内で異なる次数のフィルタを効率良く利用することができる。取り込み室内応答が、低周波数で長い測定期間を用い、高周波数に向けて漸進的に短い測定期間を用いて処理される場合には、室内補正フィルタリングは、低周波数から高周波数にフィルタリングを行う時に更に低い次数のフィルタを必要とする。この場合に、サブバンド単位の室内周波数補正フィルタリング手法は、重ね合わせ保留法又は重ね合わせ加算法を用いた高速畳み込みと類似の計算の複雑さを与えるが、サブバンド領域手法は、この計算をかなり低いメモリ要件、並びにかなり低い処理遅延しか伴わずに与える。
【００４６】
オーディオチャンネルの全てを処理した状態で、解析モジュールは、ラウドスピーカに対して特定の多チャンネル構成を自動的に選択し、リスニング環境内で各ラウドスピーカに対する位置を計算する（段階８０）。モジュールは、各ラウドスピーカからマイクロフォンの各々までの遅延を用いて、距離及び少なくとも方位角、好ましくは、定められた３Ｄ座標系内でのラウドスピーカに対する仰角を決定する。方位角及び仰角を分解するモジュールの機能は、マイクロフォンの個数及び受信信号の多様性に依存する。モジュールは、これらの遅延をラウドスピーカから座標系の原点までの遅延に一致するように再調節する。所定のシステム電子機器の伝播遅延に基づいて、モジュールは、ラウドスピーカから原点までの空気伝播に一致する絶対遅延を計算する。この遅延とサウンドの一定の速度とに基づいて、モジュールは、各ラウドスピーカまでの絶対距離を計算する。
【００４７】
各ラウドスピーカの距離及び角度を用いて、モジュールは、最も近い多チャンネルラウドスピーカ構成を選択する。室内の物理特性、ユーザの誤操作、又はユーザプリファレンスのいずれかに起因して、ラウドスピーカ位置は、サポートされた構成と正確に対応しない可能性がある。工業規格に従って適宜指定された所定のラウドスピーカの場所のテーブルが、メモリに格納される。標準のサラウンドサウンドスピーカは、例えば、大まかにゼロの仰角の水平平面にほぼ位置し、方位角を指定する。いずれかの高さのラウドスピーカは、例えば、３０度と６０度の間の仰角を有することができる。以下は、そのようなテーブルの例である。
【００４８】
（表）
0006023796.tif 000002
【００４９】
現在の工業規格は、モノラルから５．１までの約９つの異なるレイアウトを指定している。ＤＴＳ−ＨＤ（登録商標）は、現在、以下の４つの６．１構成：
−Ｃ＋ＬＲ＋ＬsＲs＋Ｃs
−Ｃ＋ＬＲ＋ＬsＲs＋Ｏh
−ＬＲ＋ＬsＲs＋ＬhＲh
−ＬＲ＋ＬsＲs＋ＬcＲc
及び７つの７．１構成：
−Ｃ＋ＬＲ＋ＬＦＥ1＋ＬsrＲsr＋ＬssＲss
−Ｃ＋ＬＲ＋ＬsＲs＋ＬＦＥ1＋ＬhsＲhs
−Ｃ＋ＬＲ＋ＬsＲs＋ＬＦＥ1＋ＬhＲh
−Ｃ＋ＬＲ＋ＬsＲs＋ＬＦＥ1＋ＬsrＲsr
−Ｃ＋ＬＲ＋ＬsＲs＋ＬＦＥ1＋Ｃs＋Ｃh
−Ｃ＋ＬＲ＋ＬsＲs＋ＬＦＥ1＋Ｃs＋Ｏh
−Ｃ＋ＬＲ＋ＬsＲs＋ＬＦＥ1＋ＬwＲw
を指定している。
【００５０】
産業が３Ｄに向けて動く時に、更に多くの工業規格及びＤＴＳ−ＨＤ（登録商標）レイアウトが定められることになる。接続したチャンネルの個数と、これらのチャンネルにおける距離及び角度とが与えられると、モジュールは、テーブルから個々のスピーカの場所を決定し、指定された多チャンネル構成に最も近い適合を選択する。「最も近い適合」は、誤差計量処理又は論理処理によって決定することができる。誤差計量処理は、例えば、特定の構成に対する正しい適合の個数を計数するか又は特定の構成におけるスピーカの全てのものまでの距離（例えば、誤差の二乗和）を計算することができる。論理処理は、最多のスピーカ適合個数を有する１つ又はそれよりも多くの構成候補を決定し、次に、いずれかの不適合に基づいて、どの構成候補が最も有望であるかを決定することができる。
【００５１】
解析モジュールは、各オーディオチャンネルに対する遅延及び利得の調節値、並びにフィルタ係数をシステムメモリに格納する（段階８２）。
【００５２】
プローブ信号は、室内応答の効率的で正確な測定、及び広いリスニング区域にわたって有効なエネルギ測度の計算を可能にするように設計することができる。第１のプローブ信号は、指定された音響バンドにわたって実質的に一定のマグニチュードスペクトルによって特徴付けられるブロードバンドシーケンスである。指定された音響バンドにわたり「一定」からのずれは、これらの周波数でＳＮＲの失損をもたらす。一般的に、設計仕様が、指定された音響バンドにわたるマグニチュードスペクトルにおける最大ずれを指定することになる。
【００５３】
プローブ信号及び取得
第１のプローブ信号Ｓの１つのバージョンは、図５ａに示す全域通過シーケンス１００である。図５ｂに示すように、全域通過シーケンスＡＰＰのマグニチュードスペクトル１０２は、全ての周波数にわたってほぼ一定（すなわち、０ｄＢ）である。このプローブ信号は、図５ｃ及び図５ｄに示すように非常に狭いピークの自己相関シーケンス１０４を有する。ピークの狭さは、マグニチュードスペクトルが一定であるバンド幅に反比例する。自己相関シーケンスのゼロ遅延値は、あらゆる非ゼロ遅延値を遥かに上回り、繰り返さない。どれ程超えるかは、シーケンスの長さに依存する。１，０２４（２10）個のサンプルのシーケンスは、あらゆる非ゼロ遅延値を少なくとも３０ｄＢ超えるゼロ遅延値を有することになり、それに対して６５，５３６（２16）個のサンプルのシーケンスは、あらゆる非ゼロ遅延値を少なくとも６０ｄＢ超えるゼロ遅延値を有することになる。非ゼロ遅延値が低い程、ノイズ除去は大幅であり、遅延はより正確である。全域通過シーケンスは、室内応答取得処理中に、室内のエネルギが、全ての周波数で同時に累積されることになるようなものである。それによって正弦波プローブを掃引するのと比較した場合に短いプローブ長が可能になる。更に、全域通過励振は、ラウドスピーカをその公称作動モードの近くで作動させる。同時に、このプローブは、ラウドスピーカ／室内応答の正確な全バンド幅測定を可能にし、非常に迅速な全体の測定処理を可能にする。２16個のサンプルのプローブ長は、０．７３Ｈｚの周波数分解能を可能にする。
【００５４】
第２のプローブ信号は、第１のプローブ信号の指定された音響バンドに部分的又は完全に重なるとすることができる特定のターゲットバンド内でのノイズの整形又は減衰に向けて設計することができる。第２のプローブ信号は、指定された音響バンドの一部分にわたって増幅されたマグニチュードスペクトルを与えるベースバンドシーケンスに適用されるプリ−エンファシス関数によって特徴付けられるプリ−エンファシスシーケンスである。このシーケンスは、音響バンドの一部分にわたって増幅されたマグニチュードスペクトル（＞０ｄＢ）を有するので、エネルギ保存の理由で音響バンドの他の部分にわたって減衰されたマグニチュードスペクトル（＜０ｄＢ）を示すことになり、従って、第１又は１次のプローブ信号としての使用には適さない。
【００５５】
図６ａに示す第２のプローブ信号ＰｅＳの１つのバージョンは、ベースバンドシーケンスに適用されるプリ−エンファシス関数が周波数（ｃ／ωｄ）に反比例するプリ−エンファシス処理が施されたシーケンス１１０であり、ここで、ｃは、音速であり、ｄは、指定された音響バンドの低周波数領域にわたるマイクロフォンの間隔である。ｆがＨｚである時に、半径方向周波数ω＝２πｆであることに注意されたい。これらの２つは、一定のスケール係数によって表されるので、交換可能に使用される。更に、簡略化のために、周波数への関数依存性を省くことができる。図６ｂに示すように、マグニチュードスペクトル１１２は、周波数に反比例する。５００Ｈｚよりも低い周波数では、マグニチュードスペクトルは＞０ｄＢである。増幅は、２０ｄＢの位置で最低周波数として制限される。低周波数における室内スペクトル測度を計算する際の第２のプローブ信号の使用は、単一のマイクロフォンの場合に低周波数ノイズを減衰させ、多マイクロフォンアレイの場合に圧力成分中の低周波数ノイズを減衰させ、速度成分の計算を改善するという利点を有する。
【００５６】
第１のブロードバンドプローブ信号及び第２のプリ−エンファシスプローブ信号を構成するのに、多くの異なる手法が存在する。第２のプリ−エンファシスプローブ信号は、第１のプローブ信号のブロードバンドシーケンスであってもなくてもよいベースバンドシーケンスから生成される。全域通過プローブ信号及びプリ−エンファシスプローブ信号を構成する方法の実施形態を図７に示している。
【００５７】
本発明の一実施形態により、プローブ信号は、好ましくは、２nというべき乗の長さを有する−π，＋πの間の乱数シーケンスを発生させることによって周波数領域で構成される（段階１２０）。乱数シーケンスを発生させる多くの公知の技術が存在し、本発明では、メルセンヌ・ツイスタアルゴリズムに基づくＭＡＴＬＡＢ（Ｍａｔｒｉｘ  Ｌａｂｏｒａｔｏｒｙ）の「ｒａｎｄ」関数を一様に分散された擬似ランダムシーケンスを発生させるのに適宜使用することができる。乱数シーケンスには、平滑化フィルタ（例えば、重なるハイパスフィルタとローパスフィルタの組合せ）が適用される（段階１２１）。ランダムシーケンスが、全域通過振幅を仮定した周波数応答の位相（φ）として使用され、周波数領域で全域通過プローブシーケンスＳ（ｆ）が生成される（段階１２２）。全域通過振幅は、Ｓ（ｆ）＝１*ｅ(j2πφ(f))であり、ここでＳ（ｆ）は共役対称である（すなわち、負の周波部分は、正の部分の複素共役であるものとして設定される）。Ｓ（ｆ）の逆ＦＦＴが計算され（段階１２４）、時間領域で第１の全域通過プローブ信号Ｓ（ｎ）を生成するように正規化され（段階１２６）、ここでｎは、サンプルの時間桁である。周波数依存（ｃ／ωｄ）のプリ−エンファシス関数Ｐｅ（ｆ）が定められ（段階１２８）、全域通過周波数領域信号Ｓ（ｆ）に付加されてＰｅＳ（ｆ）がもたらされる（段階１３０）。最低周波数でＰｅＰ（ｆ）の限界を定めるか、又はＰｅＰ（ｆ）を制限することができる（段階１３２）。ＰｅＳ（ｆ）の逆数ＦＦＴが計算され（段階１３４）、深刻なエッジ効果が存在しないことを確実にするために精査され、制限を回避しながら高レベルを有するように正規化され（段階１３６）、時間領域で第２のプリ−エンファシスプローブ信号ＰｅＳ（ｎ）が生成される。プローブ信号は、オフラインで計算し、メモリに格納することができる。
【００５８】
図８に示すように、実施形態において、Ａ／Ｖプリアンプは、各プローブ信号が、ラウドスピーカによって音波としてリスニング環境内に無音期間によって分離された非重ね合わせ時間スロットで送信されるように、持続時間（長さ）「Ｐ」を有する１つ又はそれよりも多くのプローブ信号、全域通過プローブ（ＡＰＰ）、及びプリ−エンファシスプローブ（ＰＥＳ）を送信スケジュール１４０に従ってオーディオ出力に供給する。プリアンプは、１つのプローブ信号を１つのラウドスピーカに一度に１つずつ送る。二重探査の場合には、最初に全域通過プローブＡＰＰが単一のラウドスピーカに送られ、所定の無音期間の後に、プリ−エンファシスプローブ信号ＰＥＳが同じラウドスピーカに送られる。
【００５９】
同じスピーカへの１回目のプローブ信号の送信と２回目のプローブ信号の送信の間に、無音期間「Ｓ」が挿入される。ロバストであるが、依然として高速な取得を可能にするために、１番目のラウドスピーカと２番目のラウドスピーカの間、及びｋ番目のラウドスピーカとｋ番目＋１のラウドスピーカの間の１回目の送信と２回目の送信の間には、無音期間Ｓ1,2及びＳk,k+1がそれぞれ挿入される。無音期間Ｓの最小持続時間は、取得される最大ＲＩＲ長である。無音期間Ｓ1,2の最小持続時間は、最大ＲＩＲ長とシステムを通じて仮定される最大遅延との和である。無音期間Ｓk,k+1の最小持続時間は、（ａ）取得される最大ＲＩＲ長と、（ｂ）ラウドスピーカ間で仮定される最大相対遅延の２倍と、（ｃ）室内応答処理ブロック長の２倍との和によって加えられる。プロセッサが、無音期間内に取得処理又は室内応答処理を実施し、計算を完了するのにより長い時間を必要とする場合には、異なるラウドスピーカへのプローブの間の無音を延長することができる。第１のチャンネルは、開始時点に一度、更に全ての他のラウドスピーカの後に遅延における整合性を検査するために一度の２回適宜探査される。合計システム取得長は、Ｓｙｓ＿Ａｃｑ＿Ｌｅｎ＝２*Ｐ＋Ｓ＋Ｓ1,2＋Ｎ＿ＬｏｕｄＳｐｋｒｓ*（２＊Ｐ＋Ｓ＋Ｓk,k+1）である。６５，５３６個というプローブ長及び６つのラウドスピーカの二重プローブ試験では、合計取得時間は、３１秒未満とすることができる。
【００６０】
上述したように、非常に長いＦＦＴに基づく取り込みマイクロフォン信号の逆畳み込みの手法は、オフライン処理の状況に適している。この場合に、プリアンプは、全ての取り込みマイクロフォン信号を格納して、取り込み処理が終了した後に初めて伝播遅延及び室内応答の推定を始めるのに十分なメモリを有すると仮定する。
【００６１】
室内応答取得のＤＳＰ実施では、取得処理の必要メモリ及び必要持続時間を最小にするために、Ａ／Ｖプリアンプは、マイクロフォン信号を取り込みながら、逆畳み込み及び遅延推定を適宜実施する。遅延及び室内応答の実時間推定の手法は、メモリ要件と、ＭＩＰＳ要件と、取得時間要件との間の兼ね合いに関して、異なるシステム要件に合わせることができる。
・取り込みマイクロフォン信号の逆畳み込みは、インパルス応答が時間反転プローブシーケンスである（すなわち、６５５３６個のサンプルプローブに対して６５５３６個のタップＦＩＲフィルタを有する）対応フィルタによって実施される。複雑さの軽減のために、適合フィルタリングは周波数領域で行われ、メモリ要件及び処理遅延の軽減のために、分割ＦＦＴ重ね合わせ及び保存法が、５０％重なるように使用される。
・各ブロックにおいて、この手法は、室内インパルス応答候補の特定の時間部分に対応する周波数応答候補をもたらす。各ブロックに対して逆ＦＦＴが実施され、室内インパルス応答（ＲＩＲ）候補の新しいサンプルブロックが得られる。
・更に同じ周波数応答候補から、負の周波数に対する値をゼロにし、結果にＩＦＦＴを適用し、ＩＦＦＴの絶対値を取ることにより、室内インパルス応答候補の解析包絡線（ＡＥ）の新しいサンプルブロックが得られる。一部の実施形態において、ＡＥはヒルベルト包絡線（ＨＥ）である。
・ＡＥの広域ピーク（全てのブロックにわたる）が追跡され、その場所が記録される。
・ＡＥ広域ピークの場所の前の所定数のサンプルで始めてＲＩＲ及びＡＥが記録され、それによって室内応答処理中の伝播遅延の微調整が可能になる。
・全ての新しいブロックにおいて、ＡＥの新しい広域ピークが見つかった場合には、それまでに記録されたＲＩＲ及びＡＥ候補がリセットされ、新しいＲＩＲ及びＡＥ候補の記録が始められる。
・不正検出を低減するために、ＡＥ広域ピーク検索空間が予想領域に限定され、各ラウドスピーカに対するこれらの予想領域は、システムを通じて仮定される最大遅延と、ラウドスピーカ間で仮定される最大相対遅延とに依存する。
【００６２】
次に、図９を参照すると、特定の実施形態において、Ｎ／２個のサンプルの各連続ブロック（５０％重ね合わせた）が処理されてＲＩＲが更新される。各マイクロフォンに対する各ブロックに対してＮ点ＦＦＴが実施され、長さＮ×１の周波数応答が出力される（段階１５０）。各マイクロフォン信号（非負の周波数のみ）に対する現在のＦＦＴ区分が、長さ（Ｎ／２＋１）×１のベクトルで格納される（段階１５２）。これらのベクトルは、先入れ先出し（ＦＩＦＯ）に基づいて累積され、（Ｎ／２＋１）×Ｋ次元のＫ個のＦＦＴ区分の行列Ｉｎｐｕｔ＿ＦＦＴ＿Ｍａｔｒｉｘが作成される（段階１５４）。長さＫ*Ｎ／２のサンプルの時間反転ブロードバンドプローブ信号の１組の分割ＦＦＴ（非負の周波数のみ）が事前計算され、（Ｎ／２＋１）×Ｋ次元の行列Ｆｉｌｔ＿ＦＦＴとして格納される。Ｉｎｐｕｔ＿ＦＦＴ＿Ｍａｔｒｉｘに対して、重ね合わせ及び保存法を用いたＦｉｌｔ＿ＦＦＴ行列との高速畳み込みが実施され、現在のブロックに対してＮ／２＋１個の点の周波数応答候補が与えられる。重ね合わせ及び保存法は、Ｆｉｌｔ＿ＦＦＴ＿ｍａｔｒｉｘの各周波数ビン内の値に、Ｉｎｐｕｔ＿ＦＦＴ＿Ｍａｔｒｉｘ内の対応する値を乗算し、行列のＫ個の列にわたる値を平均化する。各ブロックに対して、Ｎ点逆ＦＦＴが、負の周波数に対して共役対称拡張を用いて実施され、室内インパルス応答（ＲＩＲ）候補のＮ／２×１個のサンプルの新しいブロックが得られる（段階１６０）、ＲＩＲ候補の連続ブロックが付加され、指定されたＲＩＲ長（ＲＩＲ＿Ｌｅｎｇｔｈ）に達するまで格納される（段階１６２）。
【００６３】
更に同じ周波数応答候補から、負の周波数に対する値をゼロにし、結果にＩＦＦＴを適用し、ＩＦＦＴの絶対値を取ることにより、室内インパルス応答候補のＨＥのＮ／２×１個のサンプルの新しいブロックが得られる（段階１６４）。到着するＮ／２個のサンプルのブロックにわたってＨＥの最大値（ピーク）が追跡され、全てのブロックにわたる広域ピークを追跡するために更新される（段階１６６）。ＨＥの広域ピーク付近のＭ個のサンプルが格納される（段階１６８）。新しい広域ピークが検出された場合には、格納されているＲＩＲ候補をフラッシュして再開する制御信号が発せられる。ＤＳＰは、ＲＩＲと、ＨＥピークの場所と、ＨＥのピーク付近のＭ個のサンプルとを出力する。
【００６４】
二重プローブ手法が使用される実施形態において、プリ−エンファシスプローブ信号が同じく処理され、ＲＩＲ＿Ｌｅｎｇｔｈに達するまで格納されるＲＩＲ候補が生成される。全域通過プローブ信号におけるＨＥの広域ピークの場所は、ＲＩＲ候補の累積を始めるために使用される。ＤＳＰは、プリ−エンファシスプローブ信号に対するＲＩＲを出力する。
【００６５】
室内応答処理
取得処理が完了した状態で、室内応答のうちの長い部分が低い周波数の位置にあると考えられ、室内応答のうちの漸進的に短くなる部分が、周波数で次第に高くなる位置にあると考えられる蝸牛機構にヒントを得た時間−周波数処理によって室内応答が処理される。この可変分解能時間−周波数処理は、時間領域ＲＩＲ又は周波数領域スペクトル測度のいずれに対しても実施することができる。
【００６６】
室内応答処理の方法の実施形態を図１０に例示している。オーディオチャンネルインジケータｎｃｈがゼロに設定される（段階２００）。ＳｐｅａｋｅｒＡｖｔｉｖｉｔｙＭａｓｋ［ｎｃｈ］が真値ではない（すなわち、それ以上ラウドスピーカが結合されていない）場合は（段階２０２）、ループ処理は終了し、全ての補正フィルタを共通のターゲット曲線に調節する最終段階までにスキップする。真値ではない場合に、処理は、任意的にＲＩＲに可変分解能時間−周波数処理を適用する（段階２０４）。ＲＩＲには、時変フィルタが適用される。時変フィルタは、ＲＩＲの開始点は全くフィルタリングを受けないが、フィルタが時間と共にＲＩＲを通じて進む時に、時間と共に漸進的に小さくなるバンド幅を有するローパスフィルタが適用されるように構成される。
【００６７】
時変フィルタを構成してＲＩＲに適用する例示的な処理は、以下の通りである。
・ＲＩＲの最初の数ミリ秒を変更しないままに残す（全ての周波数が存在する）。
・ＲＩＲの開始時点から数ミリ秒後に、ＲＩＲに時変ローパスフィルタを適用し始める。
・ローパスフィルタの時間変更は、各段によって行うことができる。
○各段は、ＲＩＲ内の特定の期間に対応する。
○この期間は、前の段における期間と比較した場合に２倍に増大することができる。
○２つの連続する段中の期間は、５０％（先行する段に対応する期間の）だけ重なることができる。
○各新しい段では、ローパスフィルタは、そのバンド幅を５０％だけ減幅することができる。
・初期段における期間は、数ミリ秒前後になる。
・時変フィルタの実施は、重ね合わせ加算の手法を用いてＦＦＴ領域で行うことができ、取りわけ、以下のように行うことができる。
○ＲＩＲのうちで現在のブロックに対応する部分を抽出する。
○抽出したＲＩＲブロックに窓関数を適用する。
○現在のブロックにＦＦＴを適用する。
○現在の段のローパスフィルタのＦＦＴと同じサイズの対応する周波数ビンで乗算する。
○結果の逆ＦＦＴを計算して出力を生成する。
○現在のブロックの出力を抽出して、保存されている前のブロックからの出力を加算する。
○次のブロックとの組合せに向けて出力の残りの部分を保存する。
○これらの段階は、ＲＩＲの「現在のブロック」が時間と共にＲＩＲを通じて前のブロックに対して５０％重ねて摺動しながら繰り返される。
○ブロックの長さは、各段において増大することができ（各段に関連付けられた期間の持続時間をマッチングながら）、ある一定の段で増大を止めることができ、又は全体を通じて一様なものとすることができる。
【００６８】
異なるマイクロフォンでの室内応答が再整合される（段階２０６）。単一のマイクロフォンの場合は再整合を必要としない。室内応答は、時間領域でＲＩＲとして与えられる場合には、各マイクロフォンでのＲＩＲの間の相対遅延が回復され、ＦＦＴが計算され、整合されたＲＦＲが得られるように再整合される。室内応答が周波数領域でＲＦＲとして与えられる場合には、再整合は、マイクロフォン信号の間の相対遅延に対応する位相シフトによって達成される。全域通過プローブ信号に対する各周波数ビンｋにおける周波数応答をＨkとし、プリ−エンファシスプローブ信号に対するものをＨk,peとし、周波数への関数依存性を割愛した。
【００６９】
現在のオーディオチャンネルに対する再整合ＲＦＲから、スペクトル測度が構成される（段階２０８）。一般的に、スペクトル測度は、マグニチュードスペクトル及びエネルギ測度を含むが、これらに限定されない幾つもの手法でＲＦＲから計算することができる。図１１に示すように、スペクトル測度２１０は、遮断周波数ビンｋtよりも低い周波数におけるプリ−エンファシスプローブ信号に対する周波数応答Ｈk,peから計算されたスペクトル測度２１２と、遮断周波数ビンｋtよりも高い周波数におけるブロードバンドプローブ信号に対する周波数応答Ｈkからのスペクトル測度２１４とを混合することができる。最も単純な場合には、スペクトル測度は、遮断周波数よりも大きいＨkを遮断周波数よりも小さいＨk,peに付加することによって混合される。代替的に、必要に応じて、異なるスペクトル測度を遮断周波数ビン付近の遷移領域２１６内の重み付き平均として組み合わせることができる。
【００７０】
段階２０４において室内応答に可変分解能時間−周波数処理が適用されなかった場合には、可変分解能時間−周波数処理をスペクトル測度に適用することができる（段階２２０）。スペクトル測度には、平滑化フィルタが適用される。平滑化フィルタは、平滑化の量が周波数と共に増加するように構成される。
【００７１】
平滑化フィルタを構成してスペクトル測度に適用するための例示的な処理は、単極ローパスフィルタ差分方程式を用いて、それを周波数ビンに適用する段階を含む。平滑化は、バンド１：０〜９３．８、バンド２：９３．８〜１８７．５、バンド３：１８７．５〜３７５、バンド４：３７５〜７５０、バンド５：７５０〜５００、バンド６：１５００〜３０００、バンド７：３０００〜６０００、バンド８：６０００〜１２０００、及びバンド９：１２０００〜２４０００という９つの周波数バンド（Ｈｚで表した）内に実施される。平滑化は、可変指数関数忘却因子を用いた前方及び後方の周波数領域平均化を使用する。指数関数忘却因子の変動性は、周波数バンドのバンド幅（Ｂａｎｄ＿ＢＷ）によって決定され、すなわち、Ｃがスケール定数である時にλ＝１−Ｃ／Ｂａｎｄ＿ＢＷである。１つのバンドから次のものに遷移する場合には、λ値は、これらの２つのバンド内のλ値の間の線形内挿によって得られる。
【００７２】
最終のスペクトル測度が生成され終わると、周波数補正フィルタを計算することができる。この計算を行うには、システムに、望ましい補正された周波数応答又は「ターゲット曲線」を設けなければならない。このターゲット曲線は、あらゆる室内補正システムの特徴的なサウンドへの主な寄与ファクタのうちの１つである。１つの手法は、いずれかのユーザプリファレンスを反映する単一の共通ターゲット曲線を全てのオーディオチャンネルに対して使用することである。図１０に反映されている別の手法は、各オーディオチャンネルに対して独特なチャンネルターゲット曲線を生成して保存し（段階２２２）、全てのチャンネルに対して共通ターゲット曲線を生成する（段階２２４）ことである。
【００７３】
ステレオ又は多チャンネルの画像を補正するために、室内補正処理は、室内のラウドスピーカの各々からのサウンドのうちで最初の到来音のマッチング（時間、振幅、及び音質における）を最初に提供しなければならない。室内スペクトル測度は、測度の傾向のみが保持されるように、非常に粗いローパスフィルタを用いて平滑化される。言い換えれば、全ての室内寄与が除外又は平滑化されるので、ラウドスピーカ応答の直接経路の傾向が保持される。これらの平滑化された直接経路ラウドスピーカ応答は、各ラウドスピーカに対する別々の周波数補正フィルタの計算中にチャンネルターゲット曲線として使用される（段階２２６）。その結果、ターゲット付近のピーク及びディップしか補正しなくてもよいので、比較的小さい次数の補正濾過しか必要とされない。オーディオチャンネルインジケータｎｃｈが１だけ増分され（段階２２８）、チャンネルの合計個数ＮｕｍＣｈに対して探査され、見込まれる全てのオーディオチャンネルを処理し終えたか否かが決定される（段階２３０）。決定が偽であった場合には、次のオーディオチャンネルに対して全処理が繰り返される。決定が真であった場合には、処理は、補正フィルタに共通ターゲット曲線に対する最終調節を加える段階に進む。
【００７４】
段階２２４では、共通ターゲット曲線は、全てのラウドスピーカにわたるチャンネルターゲット曲線の平均として生成される。チャンネルターゲット曲線上には、いずれかのユーザプリファレンス又はユーザ選択可能なターゲット曲線を重なることができる。補正フィルタへのいずれかの調節は、チャンネルターゲット曲線と共通ターゲット曲線との差を補償するように加えられる（段階２２９）。チャンネル毎のターゲット曲線と共通のターゲット曲線の間の比較的小さい変化、及び大幅に平滑化された曲線に起因して、共通ターゲット曲線によって課せられる要件は、非常に単純なフィルタを用いて実施することができる。
【００７５】
上述したように、段階２０８で計算されるスペクトル測度は、エネルギ測度とすることができる。単一のマイクロフォンの様々な組合せ又は四面体マイクロフォン及び単一のプローブ又は二重プローブにおけるエネルギ測度を計算するための実施形態を図１２に例示している。
【００７６】
解析モジュールは、１つのマイクロフォン又は４つのマイクロフォンのいずれが存在するかを決定し（段階２３０）、次に、単一プローブ室内応答又は二重プローブ室内応答のいずれが存在するかを決定する（単一のマイクロフォンの場合は段階２３２、四面体マイクロフォンの場合は段階２３４）。この実施形態を４つのマイクロフォンに対して記載し、より一般的には、方法は、あらゆる多マイクロフォンアレイに適用することができる。
【００７７】
単一のマイクロフォン及び単一プローブ室内応答Ｈkの場合には、解析モジュールは、各周波数ビンｋ内のエネルギ測度Ｅk（割愛された周波数に関数依存する）をＥk＝Ｈk*ｃｏｎｊ（Ｈk）として構成し、ここでｃｏｎｊ（*）は、共役演算子である（段階２３６）。エネルギ測度Ｅkは、音圧に対応する。
【００７８】
単一のマイクロフォンと二重プローブ室内応答Ｈk及びＨk,peとの場合には、解析モジュールは、低周波数ビンｋ＜ｋtにおいてエネルギ測度ＥkをＥk＝Ｄｅ*Ｈk,peｃｏｎｊ（Ｄｅ*Ｈk,pe）として構成し、ここでＤｅは、プリ−エンファシス関数Ｐｅに対する相補的なデ−エンファシス関数である（すなわち、全ての周波数ビンｋにおいてＤｅ*Ｐｅ＝１である）（段階２３８）。例えば、プリ−エンファシス関数は、Ｐｅ＝ｃ／ωｄであり、デ−エンファシス関数は、Ｄｅ＝ωｄ／ｃである。高周波数ビンｋ＞ｋtでは、Ｅk＝Ｈk*ｃｏｎｊ（Ｈk）である（段階２４０）。二重プローブを使用する効果は、エネルギ測度での低周波数ノイズを減衰することである。
【００７９】
四面体マイクロフォンの場合には、解析モジュールは、マイクロフォンアレイにわたる圧力勾配を計算し、そこから音速成分を抽出することができる。詳述するように、低周波数における音圧と音速の両方に基づくエネルギ測度は、より広いリスニング区域にわたってよりロバストである。
【００８０】
四面体マイクロフォン及び単一プローブ応答Ｈkの場合には、各低周波数ビンｋ＜ｋtにおいて、エネルギ測度の第１の部分は、音圧成分と音速成分とを含む（段階２４２）。音圧成分Ｐ＿Ｅkは、全てのマイクロフォンにわたって周波数応答を平均化し、ＡｖＨk＝０．２５*（Ｈk（ｍ１）＋Ｈk（ｍ２）＋Ｈk（ｍ３）＋Ｈk（ｍ４））、Ｐ＿Ｅk＝ＡｖＨkｃｏｎｊ（ＡｖＨk）を計算することによって計算することができる（段階２４４）。この「平均」は、重み付き平均のいずれかの変形として計算することができる。音速成分Ｖ＿Ｈkは、４つ全てのマイクロフォンでのＨkから圧力勾配
0006023796.tif 000003
を推定し、
0006023796.tif 000004
に周波数依存の重み付け（ｃ／ωｄ）を適用してｘ、ｙ、及びｚの座標軸に沿った速度成分Ｖk_x、Ｖk_y、及びＶk_zを取得し、更にＶ＿Ｅk＝Ｖk_xｃｏｎｊ（Ｖk_x）＋Ｖk_yｃｏｎｊ（Ｖk_y）＋Ｖk_zｃｏｎｊ（Ｖk_z）を計算することによって計算される（段階２４６）。周波数依存の重み付けの適用は、低周波数におけるノイズを増幅する効果を有することになる。エネルギ測度の低周波部分はＥK＝０．５（Ｐ＿Ｅk＋Ｖ＿Ｅk）であるが（段階２４８）、重み付き平均のいずれかの変形を使用することができる。各高周波数ビンｋ＞ｋtにおけるエネルギ測度の第２の部分は、例えば、和の二乗ＥK＝｜０．２５（Ｈk（ｍ１）＋Ｈk（ｍ２）＋Ｈk（ｍ３）＋Ｈk（ｍ４））｜2、又は二乗の和ＥK＝０．２５（｜Ｈk（ｍ１）｜2＋｜Ｈk（ｍ２）｜2＋｜Ｈk（ｍ３）｜2＋｜Ｈk（ｍ４）｜2）として計算される（段階２５０）
【００８１】
四面体マイクロフォンと二重プローブ応答Ｈk及びＨk,peとの場合には、各低周波数ビンｋ＜ｋtにおいて、エネルギ測度の第１の部分は、音圧成分と音速成分とを含む（段階２６２）。音圧成分Ｐ＿Ｅkは、全てのマイクロフォンにわたる周波数応答を平均化し、ＡｖＨk,pe＝０．２５*（Ｈk,pe（ｍ１）＋Ｈk,pe（ｍ２）＋Ｈk,pe（ｍ３）＋Ｈk,pe（ｍ４））、デ−エンファシススケーリングを適用し、Ｐ＿Ｅk＝Ｄｅ*ＡｖＨk,peｃｏｎｊ（Ｄｅ*ＡｖＨk,pe）を計算することによって計算することができる（段階２６４）。この「平均」は、重み付き平均のいずれかの変形として計算することができる。音速成分Ｖ＿Ｈk,peは、４つ全てのマイクロフォンでのＨk,peから圧力勾配
0006023796.tif 000005
を推定し、
0006023796.tif 000006
からｘ、ｙ、及びｚの座標軸に沿った速度成分Ｖk_x、Ｖk_y、及びＶk_zを推定し、Ｖ＿Ｅk＝Ｖk_xｃｏｎｊ（Ｖk_x）＋Ｖk_yｃｏｎｊ（Ｖk_y）＋Ｖk_zｃｏｎｊ（Ｖk_z）を計算することによって計算される（段階２６６）。プリ−エンファシスプローブ信号の使用は、周波数依存の重み付けを適用する段階を排除する。エネルギ測度の低周波部分は、ＥK＝０．５（Ｐ＿Ｅk＋Ｖ＿Ｅk）である（段階２６８）（あるいは、他の重み付き組合せ）。各高周波数ビンｋ＞ｋtにおけるエネルギ測度の第２の部分は、例えば、和の二乗ＥK＝｜０．２５（Ｈk（ｍ１）＋Ｈk（ｍ２）＋Ｈk（ｍ３）＋Ｈk（ｍ４））｜2、又は二乗の和ＥK＝０．２５（｜Ｈk（ｍ１）｜2＋｜Ｈk（ｍ２）｜2＋｜Ｈk（ｍ３）｜2＋｜Ｈk（ｍ４）｜2）として計算される（段階２７０）。二重プローブ、多マイクロフォンの場合は、音圧及び音速からエネルギ測度を形成する段階と、音速成分を抽出するのに周波数依存のスケーリングを回避するためにプリ−エンファシスプローブ信号を使用する段階との両方を組合せて、従って、ノイズの存在下でよりロバストな音速を与える。
【００８２】
四面体マイクロフォンアレイにおいて単一プローブ技術又は二重プローブ技術を用いてエネルギ測度、取りわけ、エネルギ測度の低周波成分を構成する手法のより精密な発展形態が続行される。この発展は、多マイクロフォンアレイの利点と、二重プローブ信号の使用との両方を示している。
【００８３】
実施形態において、低周波数で室内のサウンドエネルギ密度のスペクトル密度が推定される。この点において、瞬間的なサウンドエネルギ密度は次式で与えられる（式１）。
0006023796.tif 000007
ここで、太線で記した全ての変数はベクトル変数を表し、ｐ（ｒ，ｔ）及びｕ（ｒ，ｔ）は、それぞれ位置ベクトルｒで決定される瞬間的な音圧及び音速ベクトルであり、ｃは音速であり、ρは平均空気密度である。
0006023796.tif 000008
は、ベクトルＵのｌ２ノルムを示している。解析がフーリエ変換を通じて周波数領域で行われる場合には、次式が適用される（式２）。
0006023796.tif 000009
ここで、以下の通りである。
0006023796.tif 000010
【００８４】
場所ｒ（ｒx，ｒy，ｒz）における音速は、線形オイラー方程式を用いて圧力に関連付けられる（式３）。
0006023796.tif 000011
かつ周波数領域では次式になる（式４）。
0006023796.tif 000012
∇ｐ（ｒ，ｗ）項は、ｘ、ｙ、及びｚに沿った圧力勾配の周波数ｗにおけるフーリエ変換である。以下では、全ての解析を周波数領域で行うことにし、フーリエ変換を示すｗへの関数依存性は、前と同じく割愛する。同様に、場所ベクトルｒへの関数依存性も、表記から割愛する。
【００８５】
それによって望ましい低周波数領域内の各周波数における望ましいエネルギ測度に対する表現は、次式のように書くことができる（式５）。
0006023796.tif 000013
圧力勾配を計算するのに、複数のマイクロフォンの場所における圧力の間の差を使用する技術は、Ｔｈｏｍａｓ，Ｄ．Ｃ．（２００８年）著「音響強度及びエネルギ密度の理論及び推定（Ｔｈｅｏｒｙ  ａｎｄ  Ｅｓｔｉｍａｔｉｏｎ  ｏｆ  Ａｃｏｕｓｔｉｃ  Ｉｎｔｅｎｓｉｔｙ  ａｎｄ  Ｅｎｅｒｇｙ  Ｄｅｎｓｉｔｙ）」、科学修士論文、ブリガムヤング大学に記載されている。四面体マイクロフォンアレイの場合及び図１ｂに示す特定的に選択した座標系におけるこの圧力勾配推定技術を提供する。全てのマイクロフォンは全方位性のものであり、すなわち、マイクロフォン信号は、様々な場所における圧力測定値を表すと仮定する。
【００８６】
圧力勾配は、マイクロフォンアレイによって占有される空間領域にわたって圧力場における空間変動が小さいようにマイクロフォンが位置決めされるという仮定から得ることができる。この仮定は、周波数範囲でこの仮定を使用することができる上限を設ける。この場合に、圧力勾配は、
0006023796.tif 000014
によっていずれかのマイクロフォン対の間の圧力差に近似的に関連付けることができ、ここでＰkは、マイクロフォンｋにおいて測定される圧力成分であり、ｒklは、マイクロフォンｋからマイクロフォンｌに向くベクトルであり、すなわち、
0006023796.tif 000015
であり、Ｔは、行列転置演算子であり、・は、ベクトルドット積である。特定のマイクロフォンアレイ及び特定の座標系の選択では、マイクロフォン位置ベクトルは、
0006023796.tif 000016
0006023796.tif 000017
0006023796.tif 000018
及び
0006023796.tif 000019
である。四面体アレイ内の６つ全ての可能なマイクロフォン対に着目して、過剰決定の連立方程式を最小二乗解を用いて圧力勾配の未知の成分（ｘ、ｙ、及びｚの座標に沿った）に関して解くことができる。特に全ての式が行列形式でまとめられる場合には、次の行列方程式が得られる（式６）。
0006023796.tif 000020
ここで、
0006023796.tif 000021
0006023796.tif 000022
であり、Δは推定誤差である。最小二乗的に推定誤差を最小にする圧力勾配
0006023796.tif 000023
は、次式の通りに得られる（式７）。
0006023796.tif 000024
ここで、（ＲTＲ）-1ＲTＰは、行列Ｒの左側擬似逆行列である。行列Ｒは、選択されたマイクロフォンアレイ幾何学構成及び選択された座標系原点にしか依存しない。行列Ｒの擬似逆行列の存在は、マイクロフォンの個数が次元数よりも大きい限り確実にされる。３Ｄ空間（３つの次元）内の圧力勾配の推定では、少なくとも４つのマイクロフォンが必要である。
【００８７】
圧力勾配の実際の測定、更に最終的には音速の実際の測定への上述の方法の適用性に関しては、考察しなければならないいくつかの問題が存在する。
・本方法は、位相整合マイクロフォンを使用するが、マイクロフォンの間の距離が拡大するときに、一定の周波数における僅かな位相不整合の効果は減少する。
・マイクロフォンの間の最大距離は、圧力場における空間変動がマイクロフォンアレイによって占有される空間領域にわたって小さいという仮定によって制限を受け、すなわち、マイクロフォンの間の距離は、関連の最も高い周波数の波長λよりもかなり小さくなる。Ｆａｈｙ，Ｆ．Ｊ．（１９９５年）著「サウンド強度（Ｓｏｕｎｄ  Ｉｎｔｅｎｓｉｔｙ）」、第２版、ロンドン、Ｅ＆ＦＮ  Ｓｐｏｎにより、圧力勾配の推定において有限差近似を使用する方法において、圧力勾配における５％よりも大きい誤差を回避するために、マイクロフォン間隔が０．１３λよりも小さくなければならないことが示されている。
・実際の測定では、特に低周波数でマイクロフォン信号内にノイズが常に存在することを考えると、勾配は非常にノイズの多いものになる。ラウドスピーカから異なるマイクロフォンの場所に到着する音波に起因する圧力差は、同じマイクロフォン間隔において低周波数で非常に小さくなる。速度推定において、関連の信号が低周波数における２つのマイクロフォンの間の差であることを考えると、有効ＳＮ比は、マイクロフォン信号内の元のＳＮＲと比較すると低い。事態を更に悪化させることに、速度信号の計算中に、これらのマイクロフォン差信号は、周波数に反比例する関数によって重み付けされ、実質的にノイズ増幅がもたらされる。それによって速度推定において、離間されたマイクロフォンの間の圧力差に基づく手法を適用することができる周波数領域に下限が課せられる。
・室内補正は、マイクロフォンアレイ内の異なるマイクロフォンの間の高い位相整合を仮定することができない様々な消費者ＡＶ機器において実施しなければならない。その結果、マイクロフォン間隔は、可能な限り大きくなければならない。
【００８８】
室内補正では、室内モードが優勢な効果を有する２０Ｈｚと５００Ｈｚの間の周波数領域内で、圧力及び速度に基づくエネルギ測度を得ることに関心がある。その結果、マイクロフォンカプセルの間で、約９ｃｍ（０．１３*３４０／５００ｍ）を超えない間隔が適切である。
【００８９】
音圧マイクロフォンｋにおいて受信される信号と、そのフーリエ変換Ｐk（ｗ）とを考える。ラウドスピーカ供給信号Ｓ（ｗ）（すなわち、プローブ信号）に着目し、ラウドスピーカからマイクロフォンｋまでのプローブ信号の送信を室内周波数応答Ｈk（ｗ）によって特徴付ける。この場合に、Ｐk（ｗ）＝Ｓ（ｗ）Ｈk（ｗ）＋Ｎk（ｗ）であり、ここでＮk（ｗ）は、マイクロフォンｋにおけるノイズ成分である。以下に続く式における表記の簡略化のために、ｗへの依存性、すなわち、Ｐk（ｗ）をＰkなどで表している。
【００９０】
室内補正の目的では、目指す位置は、周波数補正フィルタの計算に対して使用することができる代表的な室内エネルギスペクトルを求めることである。理想的には、システム内にノイズが存在しない場合には、代表的な室内エネルギスペクトル（ＲｍＥＳ）は、次式のように表すことができる（式８）。
0006023796.tif 000025
実際には、システム内にはノイズが常に存在することになり、ＲｍＥＳの推定値を次式のように表すことができる（式９）。
0006023796.tif 000026
非常に低い周波数では、ラウドスピーカから密に離間されたマイクロフォンカプセルへの周波数応答の間の差の二乗、すなわち、｜Ｈk−Ｈl｜2は非常に小さい。それに対して、異なるマイクロフォン内のノイズは無相関のものであり、従って、｜Ｎk−Ｎl｜2〜｜Ｎk｜2＋｜Ｎl｜2である。それによって見出されるＳＮ比が実質的に低下し、低周波数で圧力勾配がノイズの多いものになる。マイクロフォンの間の距離を拡大することにより、見出される信号の振幅（Ｈk−Ｈl）は大きくなり、従って、有効ＳＮＲが改善される。
【００９１】
関連の全ての周波数に対する周波数重み係数
0006023796.tif 000027
は＞１であり、周波数に反比例するスケールでノイズを実質的に増幅する。それによって、
0006023796.tif 000028
内により低い周波数に向う時に上方傾斜が導入される。推定エネルギ測度
0006023796.tif 000029
におけるこの低周波数傾斜を阻止するために、低周波数における室内探査には、プリ−エンファシスプローブ信号が使用される。特に、プリ−エンファシスプローブ信号は
0006023796.tif 000030
である。更に、マイクロフォン信号から室内応答を抽出する際に、送信されたプローブ信号Ｓpeではなく、元のプローブ信号Ｓを用いて逆畳み込みが実施される。この手法で抽出された室内応答は、以下のフォーマット：
0006023796.tif 000031
を有することになる。結果的に、エネルギ測度の推定演算子の修正形態は次式になる（式１０）。
0006023796.tif 000032
ノイズ増幅に関する上（式１１）。
0006023796.tif 000033
この推定演算子を用いて、速度推定に入り込むノイズ成分は
0006023796.tif 000034
に基づいて増幅されず、更に、音圧推定に入り込むノイズ成分は
0006023796.tif 000035
によって減衰され、従って、音圧マイクロフォンのＳＮＲが改善される。上述したように、この低周波数処理は、２０Ｈｚから５００Ｈｚ前後までの周波数領域内で適用される。その目指す位置は、室内の広いリスニング区域を表すエネルギ測度を得ることである。高い周波数で目指すのは、ラウドスピーカからリスニング区域までにおける直接経路及び少数の早期の反響を特徴付けることである。これらの特性は、主にラウドスピーカの構造及び室内でのその位置に依存し、従って、リスニング区域内の異なる場所の間で大幅に変化することはない。従って、高周波数では、四面体マイクロフォン信号の単純な平均（又はより複雑な重み付き平均）が使用される。得られる全体の室内エネルギ測度は、式（１２）のように書かれる（式１２）。
0006023796.tif 000036
【００９２】
これらの式は、単一プローブマイクロフォン構成及び二重プローブ四面体マイクロフォン構成におけるエネルギ測度Ｅkを構成する場合に直接関連する。特に式８は、Ｅkの低周波成分を計算するための段階２４２に対応する。式８の１番目の項は、平均周波数応答の振幅の二乗であり（段階２４４）、２番目の項は、圧力勾配に周波数依存の重み付けを適用して速度成分を推定し、振幅の二乗を計算する（段階２４６）。式１２は、段階２６０（低周波数）及び２７０（高周波数）に対応する。式１２における１番目の項は、デ−エンファシス平均周波数応答である（段階２６４）。２番目の項は、圧力勾配から推定された速度成分の振幅の二乗である。単一プローブと二重プローブの両方の場合に、低周波数測度の音速成分は、測定された室内応答Ｈk又はＨk,peから直接計算され、圧力勾配を推定する段階と、速度成分を取得する段階とは一体的に実施される。
【００９３】
サブバンド周波数補正フィルタ
最小位相ＦＩＲ補正フィルタの構造は、上述の室内スペクトル（エネルギ）測度を独立して用いた各バンドにおけるＡＲモデル推定に基づいている。解析／合成フィルタバンクは非臨界的にサンプリングされるので、各バンドは、独立して構成することができる。
【００９４】
ここで図１３及び図１４ａ〜図１４ｃを参照すると、各オーディオチャンネル及びラウドスピーカに対してターゲット曲線が与えられる（段階３００）。上述したように、チャンネルターゲット曲線は、室内スペクトル測度に周波数平滑化を適用し、ユーザ定義のターゲット曲線を選択することにより、又はユーザ定義のターゲット曲線を周波数平滑化された室内スペクトル測度上に重なることによって計算することができる。更に、補正フィルタに対する極度に厳しい要件を防ぐために、室内スペクトル測度の限界を定めることができる。チャンネル毎の中間バンドの利得を中間バンド周波数領域にわたる室内スペクトル測度の平均として推定することができる。中間バンド利得の最大値に上限（例えば、２０ｄＢ）を加えたものと、中間バンド利得の最小値から下限（例えば、１０ｄＢ）を差し引いたものとの間で室内スペクトル測度の範囲が限定される。室内スペクトル測度が深いヌル点を有する周波数バンド内に過度のエネルギを注入するのを回避するために、上限は、典型的には下限よりも大きい値である。チャンネル毎のターゲット曲線は、限界が定められたチャンネル毎の室内スペクトル測度と組み合わされ、集合室内スペクトル測度３０３が得られる（段階３０４）。各周波数ビンにおいて、室内スペクトル測度は、対応するターゲット曲線ビンによって区分化され、集合室内スペクトル測度が与えられる。サブバンドカウンタｓｂがゼロに初期化される（段階３０６）。
【００９５】
解析フィルタバンクのダウンサンプリングを模擬するために、集合スペクトル測度のうちで異なるサブバンドに対応する部分が抽出され、ベースバンドに再マップされる（段階３０８）。集合室内スペクトル測度３０３は、オーバーサンプリングフィルタバンク内の各バンドに応答して、重なる周波数領域３１０ａ、３１０ｂ、以降同じく続く領域に区分化される。各区分は、図１４ｃ及び図１４ｂにそれぞれ示す偶数及び奇数のフィルタバンクバンドに適用される間引き規則に従ってベースバンドまでマップされる。マッピングに解析フィルタの形状が含まれないことに注意されたい。可能な限り低い次数を有する補正フィルタを得ることが望ましいので、これは重要である。解析フィルタバンクのフィルタが含まれない場合には、マップされたスペクトルは、急勾配の立ち下がりエッジを有することになる。従って、補正フィルタは、解析フィルタの形状を不要に補正するために高い次数を必要とすることになる。
【００９６】
ベースバンドへのマッピングの後に、奇数又は偶数に対応する区分は、シフトされたスペクトル部分だけでなく、急転したいくつかの他の部分も有することになる。それによって高次の周波数補正フィルタを必要とすることになるスペクトル不連続部がもたらされる可能性がある。補正フィルタ次数の不要な増加を防ぐために、急転スペクトル領域が平滑化される。それと引き換えに、平滑化された領域内のスペクトルの詳細部が変化する。しかし、急転区画は、合成フィルタが既に減衰高い減衰を有しており、その結果、区分のうちのこの部分の最終スペクトルへの寄与を無視することができる領域内に常に存在することに注意されたい。
【００９７】
再マップされた集合室内スペクトル測度に対して自己回帰（ＡＲ）モデルが推定される。室内スペクトル測度の各区分は、ベースバンドまでマップされて間引き効果を模擬した後に、何らかの同等スペクトルとして解釈される。従って、その逆フーリエ変換は、対応する自己相関シーケンスになる。この自己相関シーケンスは、所定のエネルギスペクトルに最小二乗的に最適に適合する望ましい次数のＡＲモデルを計算するレビンソン−ダービンアルゴリズムへの入力として使用される。このＡＲモデル（全極）フィルタの分母は、最小位相多項式である。各サブバンドに対する周波数補正フィルタの長さは、全体の室内エネルギ測度の作成中に着目した対応する周波数領域内で室内応答の長さによって大まかに決定される（長さは、低周波数から高周波数に移行するときに比例して短くなる）。しかし、最終的な長さは、経験的に微調整するか、又は残存電力を観察し、望ましい分解能に達した時に停止するＡＲ次数選択アルゴリズムの使用によって自動的に微調整するかのいずれかとすることができる。
【００９８】
ＡＲの係数は、最小位相全ゼロサブバンド補正フィルタの係数にマップされる。このＦＩＲフィルタは、ＡＲモデルによって得られるスペクトルの反転スペクトルに従って周波数補正を実施することになる。異なるバンドの間でフィルタを適合させるために、補正フィルタの全ては適宜正規化される。
【００９９】
サブバンドカウンタｓｂが増分され（段階３１６）、サブバンドの個数ＮＳＢと比較されて（段階３１８）、次のオーディオチャンネルに対する処理が繰り返されるか、又は補正フィルタのチャンネル毎の構成が終了する。この時点で、チャンネルＦＩＲフィルタ係数を共通ターゲット曲線に対して調節することができる（段階３２０）。調節されたフィルタ係数は、システムメモリに格納され、図３に示す各オーディオチャンネルに対するＰ個のデジタルＦＩＲサブバンド補正フィルタを実施するように、１つ又はそれよりも多くのプロセッサを構成するのに使用される。
【０１００】
補記Ａ：ラウドスピーカの局在化
完全に自動化されたシステム較正及び設定に関して、室に存在するラウドスピーカの場所及び個数の把握情報を有することが望ましい。距離は、ラウドスピーカからマイクロフォンアレイまでの推定伝播遅延に基づいて計算することができる。ラウドスピーカとマイクロフォンアレイの間の直接経路に沿って伝播する音波を平面波によって近似することができると仮定すると、マイクロフォンアレイによって定められる座標系の原点に対する到来角（ＡＯＡ）、仰角をアレイ内の異なるマイクロフォン信号の間の関係を観察することによって推定することができる。推定されたＡＯＡから、ラウドスピーカの方位角及び仰角が計算される。
【０１０１】
ＡＯＡを決定するのに、原理的にはラウドスピーカからマイクロフォンカプセルの各々までの周波数応答の各ビンにおける位相の間の比を拠り所とする周波数領域に基づくＡＯＡアルゴリズムを使用することができる。しかし、Ｃｏｂｏｓ，Ｍ．、Ｌｏｐｅｚ，Ｊ．Ｊ．、及びＭａｒｔｉ，Ａ．（２０１０年）著「四面体マイクロフォンアレイを用いた３Ｄ  ＤＯＡ推定における室内残響の効果に関して（Ｏｎ  ｔｈｅ  Ｅｆｆｅｃｔｓ  ｏｆ  Ｒｏｏｍ  Ｒｅｖｅｒｂｅｒａｔｉｏｎ  ｉｎ  ３Ｄ  ＤＯＡ  Ｅｓｔｉｍａｔｉｏｎ  Ｕｓｉｎｇ  Ｔｅｔｒａｈｅｄｒａｌ  Ｍｉｃｒｏｐｈｏｎｅ  Ａｒｒａｙ）」、第１２８回ＡＥＳ会議、ロンドン、ＵＫ、２０１０年５月２２〜２５日に示すように、室内反響の存在は、推定ＡＯＡの精度に対してかなりの影響を有する。周波数領域に基づくＡＯＡアルゴリズムの代わりに、プローブ信号と対になった解析包絡線手法を使用することによってもたらすことができる本出願人の直接経路遅延推定の精度を拠り所とする時間領域手法がＡＯＡ推定に向けて使用される。四面体マイクロフォンアレイを用いてラウドスピーカ／室内応答を測定することにより、各ラウドスピーカから各マイクロフォンカプセルまでの直接経路を推定することが可能になる。これらの遅延を比較することにより、ラウドスピーカを３Ｄ空間内に局在化させることができる。
【０１０２】
図１ｂを参照すると、方位角θ及び仰角φは、ラウドスピーカから四面体マイクロフォンアレイまで伝播する音波の推定到来角（ＡＯＡ）から決定される。ＡＯＡの推定のためのアルゴリズムは、２つのベクトルの間の角度を特徴付けるベクトルドット積に基づいている。特に、特定的に選択された座標系原点を用いて、ドット積式を以下のように書くことができる（式１３）。
0006023796.tif 000037
ここで、ｒlkは、マイクロフォンｋをマイクロフォンｌに接続するベクトルを示し、Ｔは、行列／アレイ転置演算を示し、
0006023796.tif 000038
は、平面音波の到来方向に位置合わせした単項ベクトルを表し、ｃは、音速を示し、Ｆｓは、サンプリング周波数を示し、ｔkは、マイクロフォンｋへの音波の到来時間を示し、ｔlは、マイクロフォンｌへの音波の到来時間を示している。
【０１０３】
図１ｂに示す特定のマイクロフォンアレイに対して、
0006023796.tif 000039
が成り立ち、この場合、
0006023796.tif 000040
0006023796.tif 000041
0006023796.tif 000042
及び
0006023796.tif 000043
である。全てのマイクロフォン対に対する式を集計すると、以下の行列方程式が得られる（式１４）。
0006023796.tif 000044
この行列方程式は、最小二乗法によって解くことができる過剰決定の線形連立方程式を表し、到来方向ベクトルに対して以下の表現がもたらされる（式１５）。
0006023796.tif 000045
方位角及び仰角は、正規化ベクトル
0006023796.tif 000046
の推定座標から
0006023796.tif 000047
及び
0006023796.tif 000048
として得られ、ａｒｃｔａｎ（）は、第４象限逆正接関数であり、ａｒｃｓｉｎ（）は、逆正弦関数である。
【０１０４】
時間遅延推定値を用いて最終的に達成可能なＡＯＡアルゴリズムの角精度は、遅延推定値の精度と、マイクロフォンカプセルの間の間隔とによって制限される。カプセル間の小さい間隔は、低い達成可能精度を意味する。マイクロフォンカプセルの間の間隔は、速度推定の要件、並びに最終製品の審美性によって上側から制限される。従って、望ましい角精度は、遅延推定精度を調節することによって達成される。必要とされる遅延推定精度が、サンプリング間隔の数分の１になる場合には、室内応答の解析包絡線は、室内応答の対応するピークの前後で内挿される。新しいピークの場所は、ＡＯＡアルゴリズムによって使用される新しい遅延推定値をサンプル精度の数分の１の精度で表している。
【０１０５】
本発明のいくつかの例示的な実施形態を図示して説明したが、当業者には、多くの変形及び代替実施形態が想起されるであろう。そのような変形及び代替実施形態は考えられており、かつ添付の特許請求の範囲に定める本発明の精神及び範囲から逸脱することなく作ることができる。
【符号の説明】
【０１０６】
１０  多チャンネルオーディオシステム
１２  多チャンネルスピーカ構成
１４  リスニング環境
２２  Ａ／Ｖプリアンプ
３２  マイクロフォン送信ボックス
【図１ａ】
0006023796.tif 000049
【図１ｂ】
0006023796.tif 000050
【図２】
0006023796.tif 000051
【図３】
0006023796.tif 000052
【図４】
0006023796.tif 000053
【図５ａ】
0006023796.tif 000054
【図５ｂ】
0006023796.tif 000055
【図５ｃ】
0006023796.tif 000056
【図５ｄ】
0006023796.tif 000057
【図６ａ】
0006023796.tif 000058
【図６ｂ】
0006023796.tif 000059
【図７】
0006023796.tif 000060
【図８】
0006023796.tif 000061
【図９】
0006023796.tif 000062
【図１０】
0006023796.tif 000063
【図１１】
0006023796.tif 000064
【図１２】
0006023796.tif 000065
【図１３】
0006023796.tif 000066
【図１４ａ】
0006023796.tif 000067
【図１４ｂ】
0006023796.tif 000068
【図１４ｃ】
0006023796.tif 000069