JPH10509565A

JPH10509565A - 録音及び再生システム

Info

Publication number: JPH10509565A
Application number: JP8507891A
Authority: JP
Inventors: フェリープオーデュナ−バスタマンテ; オールキルケビ; 晴夫濱田; フィリップアーサーネルソン
Original assignee: Adaptive Audio Ltd
Current assignee: Adaptive Audio Ltd
Priority date: 1994-08-25
Filing date: 1995-08-24
Publication date: 1998-09-14
Anticipated expiration: 2015-08-24
Also published as: WO1996006515A1; US5862227A; JP3913775B2; EP0776592B1; DE69525163T2; AU3350495A; GB9417185D0; EP0776592A1; DE69525163D1

Abstract

(57)【要約】複数のラウドスピーカにより再生するために録音するか、又は複数のラウドスピーカにより再生するために音を処理する方法が記載され、再生された音は、聴取者にはラウドスピーカから間隔をあけて設けられている仮想音源から生じているように聞こえるものもある。フィルタ手段（Ｈ）は、録音又はラウドスピーカに伝達するために録音された信号の処理のいずれかに使用され、このフィルタ手段（Ｈ）はフィルタ設計ステップで作成される。ステップにおいて、ａ）ラウドスピーカにより録音を再生すると聴取者の意図された位置で再生される信号（ｗ）と、意図された位置での所望の信号（ｄ）との間の誤差を最少にするために技術が使用される。またｂ）聴取者の位置で作られるべき前記所望の信号（ｄ）は、仮想音源の所望位置で、音源によって、前記意図した位置における聴取者の耳（または領域）で作られるであろう信号（または予測信号）により画定される。最小自乗法は、聴取者の意図された位置で再生された信号と所望の信号との誤差を平均した時間を最少にするように使用可能であり、または周波数領域に加えることができる。

Description

【発明の詳細な説明】録音及び再生システムこの発明は録音と再生システムに関するものである。序論発明は音の収音と再生のための新しい手段を提案する。ここに述べられる手段は一般に多チャンネルディジタル信号処理技術に基づいており、これは、２つかあるいはそれ以上のラウドスピーカを用い、通常の多チャンネル再生方式による音場再生のための収音手段を改善することに直接に採用することができる。使用された技術は、多数のラウドスピーカによって再生するために一般に録音された音声信号の処理にも拡張することができ、また特別なときは録音された音声信号は単一チャンネルでもあり得る。それらの信号の再生を改善するために、録音された音声信号を処理するという再生過程で、ディジタルフィルタが使用されるということに関する一般的なアプローチは参考文献「1,2」で記述されている。不十分な位置に設置されたラウドスピーカによる既存の２チャネル録音の音声再生を補正するためにこれらの技術を利用することについては、参考文献「3,4」に記述された。その後者において、 ‘仮想’ラウドスピーカ配置の考えは紹介された。再生信号と‘所望’の信号における時間平均された自乗誤差の和が最小であることを保証するために、再生の際、録音された信号に対して作用するために、信号処理系ではフィルタが使用された。所望の信号は今度は、音場の音源によって特定された明確な位置に作り出されるものとして特定された。フィルタの動作によって、再生信号は所望の信号と良好に一致し、これにより、受聴者のなかに‘仮想’音源の位置から発生する音のイリュージョンが作られる。本発明は再び仮想音源の考えを利用する。本発明の目的は、２つ（またはそれ以上）のラウドスピーカを介して再生するための録音手段を提供するものであり、それによって、実際のラウドスピーカの位置からは隔たった特定の空間位置から到来する音のイリュージョンを受聴者に対して作り出す。再生におけるこの目的を達成するための技術を最初に記述したのはＡtal＆Ｓchr ocder「5」であり、彼らは‘２つだけのラウドスピーカによって任意に位置に配置された音像’を作り出すための手段を提案した。‘擬似音源トランスレータ’ と名付けられた彼らの発明のなかで、Ａtal＆Ｓchrocderはまた、２つのラウドスピーカに入力するのに先だって、単一（チャンネル）の信号を操作するためにフィルタネットワークを用いた。本発明における１つの局面によって、我々は複数のラウドスピーカによる再生のための収音手段、もしくは複数のラウドスピーカによる再生のための音声処理に関する手段を提案しており、そのなかで再生されるいくつかの音は受聴者に対して仮想音源から発生されているがごとく現れ、ラウドスピーカに供給するための信号を録音したり処理するためにフィルタ手段（Ｈ）を用いて成り立っており、フィルタ手段（Ｈ）はフィルタ設計の過程で作られるものであり、そのフィルタ設計過程とは以下のように特徴づけられる： (a) 収音したものをラウドスピーカを介して再生する際、意図される位置の受聴者に対して再生される信号（w）と、意図される位置における所望の信号（d）との間の誤差を最小にするように用いられる技術で、そのなかで、 (b)受聴者に対して再生される前記所望の信号(d)は、仮想音源の所望の位置での音源によって前記の意図される位置の受聴者の耳元(あるいは受聴者の領域)において作り出されるであろう信号(あるいは推定信号)で定義される。好ましくは、所望信号は、仮想音源の前記所望の位置と再生音場における特定位置、すなわち受聴者の耳元、あるいは頭部の領域、との間の伝達関数であるフィルタ（Ａ）の形式で明確に記述されることが望ましい。伝達関数を導出する方法は様々であるが、しかし、好ましくは伝達関数は実音源への入力と受聴者の‘頭部伝達関数’(Ｈead Ｒelated Ｔransfer Ｆunctions,ＨＲＴＦ)の影響をモデル化するために用いるダミーヘッドの耳元（あるいは頭部の領域）におけるマイクロフォンからの出力との間のまず第一に実施される測定によって導出される。最小自乗法は、意図される受聴者の位置における再生信号と所望の信号との時間平均誤差を最小にするように採用される。また、最小自乗法は時間領域よりはむしろ周波数領域において採用される。トランスデューサ関数は、実際の受聴者による測定か、または解析上、あるいは経験的に得られる頭部伝達関数（ＨＲＴＦ）のモデルによって導出される。好ましくは、再生のためにラウドスピーカに入力するに先だって仮想音源を処理するために使用されるフィルタは‘クロストークキャンセル’のマトリックスと共に所望の信号を特定する伝達関数を表しているディジタルフィルタの畳み込みによって導出される。単一の逆フィルタ設計手順（数値計算的に集中している）のみが必要とされる。発明の最初の局面のとおりに手法を利用した結果は、２つのラウドスピーカのみを使用した場合に、受聴者は音像が、受聴者の耳を含む平面上のほとんどすべて位置にある任意の場所に配置された仮想音源から到来するように知覚する。しかしながら、（受聴者の前方である）この平面の前方の円弧において、仮想音源をおくことが特に効果的であるようにシステムは開発された。受聴者の後方に２つの追加されたラウドスピーカを使用することにより、受聴者の横方向、あるいは後方に仮想音源を創成することが可能となる。発明の１つの利用法は、より向上した２チャンネルの音声収音を行うための手段を提供することである。前述したすべてのフィルタ設計の過程は、それ以上の必要な処理をすることはなく後に続く２つのラウドスピーカを介するトランスミッションに対してそのまま使用できる２つの収音信号を発生するために実行することができる。従って、この発明の２番目の局面は、従来の多チャンネル音場再生方式をとおして収音された音声を再生することを可能にする、マルチチャンネルの収音の手段を提案することであり、その手段は前述のフィルタ設計手順を利用する。収音された音声は、コンパクトディスクやアナログ、あるいはディジタルオーディオテープ、その他いかなる適当な通常の媒体に録音されることが可能であることは明確である。重畳により、異なる楽器、ボーカルに属性を与える目的で、異なる複数の位置に仮想音源を配置させて、このような録音は実施される。従って、２つのラウドスピーカによる再生のための録音のプロダクションが向上する。以下の図を参照しながら、例によって、この発明の様々な具体的側面について記述していく。図１は仮想音源配置のための信号処理を、(a)に概要図解で、(b)にそのブロック図で示しており、図２にクロストークキャンセルフィルターの設計行列を示す。Ｈ_x21，Ｈ_x12とＨ_x22が最小二乗法に基づいて設計される。これは、再とを確証する。従って、w₁(n)と w₂(n)は、それぞれ信号u₁(n)とu₂(n)が単純に遅延をもったものであり、図３はラウドスピーカの位置補正問題のアウトラインを図(a)に示しており(b)ではブロックダイアグラムを表している。u₁(n)とu₂(n)は、従来のステレオで録音された信号を示すことに留意する。つまりディジタルフィルタＡ₁₁,Ａ₂₁，Ａ₁₂ とＡ₂₂は‘理想的’に配置された仮想音源に対する入力から受聴者までの伝達関数を示しており、図４は仮想音源の等化領域に対する主観評価を行った際のレイアウトである。仮想音源は受聴者に対面する対のスピーカを介して模擬されている。黒いカーテンは音源が実際どこにあるか分からないようにするために用いられた。スクリーンの外に描かれた円には異なる角度の定位のために仮想と付加的な実音源が置かれる距離がマークされており、図５は無響室における電気音響システムのインパルス応答を示しており、 a)はＬchのスピーカから左耳、b)はＬchのスピーカから右耳、c)はＲchのスピーカから左耳、d)はＲchのスピーカから右耳までの応答であり、図６は無響室で測定したクロストークキャンセル用のフィルタ行列におけるインパルス応答を示しており、a)ｈ₁₁(n)、b)h₂₁(n),c)ｈ₁₂(n),d）h₂₂(n)によって示されており、図７は無響室内での電気音響システムのインパルス応答とクロストーク相殺用の逆フィルタを畳み込んだ結果を表しており、図８と図９は共に無響室内で音声信号に対する、a)仮想音源と、b)実音源の定位実験の結果をを表しており、図10 はリスニングルームでのインパルス応答を表し：a)Ｌchスピーカから左耳までの系、b)Ｌchスピーカから右耳までの系、c)Ｒchスピーカから左耳までの系、 d)Ｒchスピーカから右耳までの系であり、図11 はリスニングルーム内において用いられたクロストークキャンセルフィルタのインパルス応答を示しており、a)h₁₁(n)，b)h₁₂(n)，c)h₂₁(n)，d)h₂₂(n)，によって表され、図12 はリスニングルームでの電気音響系のインパルス応答とクロストーク相殺用逆フィルタを畳み込んだ結果を表しており、図13と図14 は共にリスニングルームでの音声信号に対する、a)仮想音源と、b)実音源の定位実験の結果を表しており、図15 は車室内での主観評価に用いられたスピーカとダミーヘッドマイクロフォンとのレイアウトを示しており、a)頭上視界、b) 側面視界であり、図16 は車室内の運転者席（左側運転車）に配置されたダミーヘッドマイクロフォンの両耳と前面に置かれた対のスピーカ間のインパルス応答を実測したものであり、図17 は車室内におけるクロストークキャンセルフィルタのインパルス応答を示しており、図18 はクロストークキャンセルフィルタの入力からダミーヘッドの耳元のマイクロフォンまでのインパルス応答を表している。これらの結果は図17で示した逆フィルタと図16での車室内でのインパルス応答を畳み込んだ結果であり、図19 は車室内実験における仮想音源配置に関する主観評価を示しており、図20 は逆システムと所望の特性のフィルタデータベースを用いた無響室内での主観評価を表している。±45度と±135度の音源が仮想音源を生成するに当たって用いられている。実音源は165度、-150度、135度を除いた上述の全ての角度に配置された。仮想音源は135、150、-165度を除いた上述のすべての音源位置に配置された。音源はＫemarモデルのダミーヘッドマイクロフォンの中心から2.2m離れた円上に配置され、図21 は仮想音源を生成するために音声信号（スピーチ）や４つの音源を用いた無響室内での定位実験の結果を表しており、a)仮想音源の結果、b)実音源の結果である。２つのラウドスピーカを用いた単一仮想音源再生のための信号処理技術離散時間信号u(n)は、受聴者に対して我々が配置させたい任意の位置での仮想音源信号を定義する。信号Ｄ₁(z)とＤ₂(z)は仮想音源によって生じる受聴者の耳元での所望の信号である。ディジタルフィルタＡ₁(z)と′Ａ₂(z)は、配置された仮想音源から受聴者の耳元までの伝達特性を定義する。Ｚ変換の領域において、これらの伝達関数は高品質なラウドスピーカへの入力（またはラウドスピーカの位置に設置された高品質のマイクロフォンにより測定された音圧）とダミーヘッドの耳元に設置された高品質のマイクロフォンの出力との間の伝達関数を測定することで導出することができる。このような実験手順は無響条件下において仮想音源位置の範囲で実行される。これにより、仮想音源位置の範囲を網羅した頭部伝達関数(ＨＲＴＦ)のデータベースが導出される。別の方法としては、このデータベースは分析的、経験的なＨＲＴＦのモデルを使うことによって得ることもできる。図１に戻って、信号Ｖ₁(n)とＶ₂(n)は再生の際用いられているラウドスピーカへの入力を定義する。これらの信号は収音信号を形成している。スピーカから受聴者の耳までの伝達過程において、収音信号はＣ₁₁(z),Ｃ₁₂(z),Ｃ₂₁(z)とＣ₂₂(z) を要素とする電気音響伝達関数行列を経て伝搬する。これらの伝達関数は信号v₁ (n)及びv₂(n)と受聴者の耳元で再生された信号w₁(n)及びw₂(n)とを関係づける。そこでｚ領域では(2)の様に書き記すことができる。伝達特性Ａ₁(z)、Ａ₂(z)などと同様に、無響室において２つのラウドスピーカへの入力とダミーヘッドの両耳のマイクロフォンの出力との伝達関数を測定することにより、伝達関数、Ｃ₁₁(z),Ｃ₁₂(z),Ｃ₂₁(z)とＣ₂₂(z)を導出することができる。さらに付記するとこれらの伝達関数を示すために他の技術を用いる場合もある。収音に適切な信号処理過程を導出することで、これらの伝達関数を示すのに用いられるフィルタは、収音されたものが再生される場合に起こるような伝達関数の良好な標本であることが保証される必要があることは明らかである。伝達特性Ｃ₁₁(z),Ｃ₂₁(z),Ｃ₁₂(z)とＣ₂₂(z)が今、理想的に表記できたと仮定すると、その際仮想音源信号u(n)を操作する逆フィルタＨ₁(z)とＨ₂(z)を導出することが可能となる。これにより２つのラウドスピーカ再生のために準備の整った録音されるべき信号v₁(n)及びv₂(n)の生成が可能となるフィルタ係数（ただし有限長インパルス応答ィジタルフィルタと仮定される）を設計するにあったって我々は最小自乗法を用いている。これについては参考文献「１、２、３、４」に概要が述べられている。これらのフィルター設計手順は十分に参考文献「1,2,3,4 」で記述されているので、ここでは繰り返さない。しかしながら、最小自乗法によるアプローチは(3)式で与えられるコスト関数Ｊを最小とするような最適なフィルタ係数を得るために用いられていることを強調しておく。ここで、Ｅ[]は期待値演算子である。この最小自乗法に基づいたと再生された信号w1(n)とw2(n)の差違で定量化される値を時間平均された自乗平均の意味で小さくすることができることを付記しておく。いくつかの条件下で逆問題を改善する意味で、式(3)で示したコスト関数に新たに別の項を加えることは有用である。この項はフィルタＨ₁(z)とＨ₂(z)で用いられている係数の絶対値の自乗によって与えられる。この設計手法は参考文献「3,4」でより詳しく述べられている。しかしながら収音技術として有用であるためには、要求されるそれぞれの仮想音源に対して逆フィルタを設計することは不可欠である。フィルタ設計手法はとても冗長なため（特に高音質における高い周波数においては）、それぞれの位置に対応するフィルタを設計するには大変時間を消費する。ここに述べる他の技術は逆フィルタのマトリックスを使って、スピーカ１から受聴者の耳２とスピーカ２から受聴者の耳１への‘クロストーク’を最小とすることを保証するフィルタ設計法である。ここでもまた、参考文献「1,2」で記述されているこの‘クロストーク相殺用’の行列を設計するために最小自乗技術が使われている。この手法は（図２に示されている）式(4)が良好な近似であることを保証するために用いられている。ここで、z^- ^ΔはΔサンプルのモデリングディレイを意味する。ひとたび、フィルタ行列Ｈ_x11(z),Ｈ_x12(z),Ｈ_x21(z)とＨ_x22(z)が最小自乗法に基づいて設計されれば、フィルタＨ₁(z)とＨ₂(z)は、要求される仮想音源のそれぞれの位置に対応する所望信号を特定するのに用いられるフィルタＡ₁(z)とＡ₂(z)のそれぞれのペアーへと容易に導出される。これは式(4)を用いた事実に基づいており、我々はつぎの近似を行うことができる。従って我々は、フィルタＨ₁(z)とＨ₂(z)をから導出し、式（７）を導くことができる．再生信号は式(8)によって与えられるため、それに従い式(6)からＨ₁(z)とＨ₂(z)は与えられ再生信号は式(9)となる．言い換えれば、再生信号はΔサンプル分の遅延を付加した所望の信号と同等の良好な近似を有していることになる．この付加的な遅延を別にすれば、信号を再生する目的と合致している。仮想音源に着目した信号の再生に帰着する。従って、はじめにクロストーク相殺用のフィルタを設計することで、フィルタＨ₁ (z)とＨ₂(z)は、与えれた仮想音源の特性を有しているフィルタＡ₁(z)とＡ₂(z) のインパルス応答とクロストーク相殺用フイルタマトリックスの適切な要素のインパルス応答との単純に畳み込みの結果設計される。インパルス応答は小文字を用いて以下のように表される．ここで記号＊は畳み込みを意味する．これらのインパルス応答を得るために必要とされる数値計算は、係数h₁(n)とh₂( n)がＨ₁(z)とＨ₂(z)の最適設計の最小自乗問題を解くことによって得るのに必要な数値計算と比較して大幅にその量が軽減される。３．ラウドスピーカ配置の補正システムへのフィルタ設計手順の拡張先に概説したフィルタ設計手順は、この発明に従って、ラウドスピーカの配置補正のための逆フィルタ設計を支援することに利用することができる。これらについては参考文献「3」と「4」で詳細に述べられている。この場合、目的は簡易的に録音された２チャンネルの信号を操作するようなフィルタマトリックスの設計となる。フィルタは従来のステレオ録音された信号が、受聴者の耳元で”仮想音源”として最良な再生となるよう設計される。このシステムに関するブロックダイアグラムを図３に示す。再び我々は、式(4)を用いて式(12)の様に表す。従って、逆フィルタ行列は次式に従い、さらにフィルタの設計手法は初めにクロストーク相殺用のフィルタ行列を設計することで単純化される。これは先に述べた同じ理由から説明することができる。しかしながら、この場合については再生信号が式(14)によって与えられ、さらに、式(13)にしたがって設計された逆フィルタにより、次式が得られる。再生された信号は所望の信号が単純に遅延されたものとなりスピーカ配置補正システムの目的と合致する。４．多チャンネルスピーカ再生方式への技術の拡張上述された仮想音源配置に関してのアプローチは２つ以上のスピーカを用いた音場再生システムとして容易に拡張することができる。Ｌ個のラウドスピーカが再生に用いられると仮定する。受聴者の頭部の領域において、Ｍ個の位置に生成される所望の信号を明示すると仮定し、これらは仮想音源の位置とダミーヘッド（あるいは分析的、経験的に明示されるもの）の頭部の位置（あるいは領域）の間の伝達関数のＭ次のベクトルを測定することにより導出される。我々は以下に与えられるようにベクトルを定義した。また所望の信号のベクトルは（頭部の領域において）Ｍ個の位置のでの再生信号のベクトルと、式(18)に示すスピーカに対する入力信号とを関連づけるところの伝達特性を測定ないし明細する。ここでw(z)とy(z)は以下で与えられるように定義される。さらに、Ｃ(z)が以下の行列で与えられるその際、逆フィルタベクトルは(22)式によって明示される。Ｍ>Ｌの時、逆フィルタは参考文献「1,2」で示された技術により導出され、このことからコスト関数を最小とするような最適なＦＩＲフィルタ係数が得られる。ここで、e_m(n)は所望の信号d_m(n)とＤＨＭがあるところのでの再生信号z_m(n)との誤差を表している。これは、しかしながら、強力な数値計算を必要とする作業である。しかし、もし測定ポイントＭがラウドスピーカの数Ｍと同じ分だけ用意されたならば、我々は上述した内容での効率的なフィルタ設計法を用いることができる。まずはじめに次式で示す良好な近似を与えるＬ×Ｌの逆フィルタについて定義したことを留意する。ここで、Ｉは単位行列である。この関係を用いれば近似(25)が得られるその際逆フィルタベクトルはとなりさらにそして再生信号はつぎのように与えられるため、次式が誘導できる。Ｍ＝Ｌの際における再生信号は所望の信号に単純に遅延をもうけたものと等化である。この手順はクロストーク相殺フィルタＨ_x(z)の設計に依存する。これは聴取者頭部の領域でＭ＝Ｌポイントにおいて作り出された信号がこれらの信号の単に遅延しているものになることを保証するために信号u(z)のＬベクトルに作用するマトリックスＨ_x(z)として定義される。換言すれば、クロストーク・キャンセレーション・フィルタを設計するために使われた所望信号はつぎで与えられ、さらにこの場合の再生信号は次式で与えられる。行列Ｈ_x(z)は、広範囲にわたって参考文献「1,2,3,4」で記述されている技術を利用して設計される。５．周波数領域フィルタ設計技術前述の参考文献「1-4」で参照されている時間領域における最小自乗法による逆フィルタ設計法に加えて、逆フィルタを周波数領域で設計することは同じく可能である。これは特にクロストークキャンセル行列を設計する場合、とりわけ多数のラウドスピーカを利用して行う再生においてはより効果的である。しかし、もし周波数領域による設計を効率よく行うためには多くの手順を必要とする。まず初めに、逆フィルタの実現不能性の問題が、時間領域におけるそれと同様にモデリングディレイの最適な選択という点において考慮されなければならない。次にインバーションの不良条件に関する問題であり、これは周波数領域において行う場合に明確に考慮に入れられなければならない。これは時間領域において最小自乗解を求めるために適応アルゴリズムが使われるとき、本質的に避けられる。周波数領域における設計の技法は、可能な不良条件問題を例証する単一チャンネルでの例を使って最も容易に説明される。もし、われわれが仮に電気音響的な伝送系Ｃ(z)を有する時、逆フィルタＣ(z)を設計することは、単純に1/Ｃ(z)を算出することである。明らかに、もしＣ(z)が最小位相（１つかそれ以上の零点をz-平面の単位円の外側に有する時）ではない場合は、1/Ｃ(z)はその後の過程に不安定である。（Ｃ(z)の零点は1/Ｃ(z)極であり、単位円の外側にあるため）しかし、前向き時間での1/Ｃ(z)の不安定な応答は後ろ向き時間における安定な応答であると解釈される。すなわち、1/Ｃ(z)は安定であるが非因果なインパルス応答を有するといえる。非因果なインパルス応答に関する問題はモデリングディレイを加味することで部分的には補正することができる。それにより、正の時間における方向のΔサンプルにより逆フィルタのインパルス応答を効果的にシフトすることで、原則としてz^- ^Δ/Ｃ(z)からＨ(z)を算出する。しかしながら、もしＣ(z)の零点の１つか単位円の外側にあり、単位円に近接している時、逆の時間におけるインパルス応答の減衰は遅いであろう（極は僅かに押さえられている）。これは、ゼロより小さい時間の値として生じる“理想的な”逆フィルタのインパルス応答において十分なエネルギーを有するという結果となる。同様に、もし単位円の中のＣ(z)の零点の１つが単位円に近接している場合、前向き時間でのインパルス応答の減衰は遅く、また要求される逆フィルタは、非常に長いインパルス応答を前向き時間において有する。この問題を軽減するのを補助するための技術として、逆フィルタの設計を“組織化(regularise)”するために新たにパラメータが提案される。これにより、逆フィルタの極は操作され、そして単位円から離れていく。従って、正と負の時間において逆フィルタのインパルス応答は短縮される。この議論は、ある特定の例題を考えることで、単一チャンネルの場合について明確に説明される。我々は誤差信号（所望の信号と再生信号の差違）のフーリエ変換の自乗絶対値を最小し、さらにこれに付加して逆フィルタからの出力信号のフーリエ変換の自乗絶対値に比例する項をもつコスト関数を定義する。我々は次のコスト関数をを最小にするように探索する。ここで、βがレギュラリゼーションパラメータであり、インバージョンを行う逆フィルタで用いる”努力(effort)”の項を重みづける。この式で使われたフーリエ変換はz＝e^j ^ωの代入を行った前に使われたz- 変換に関連している。すなわち、例えば、Ｄ(e^j ^ω)はＺ変換Ｄ(z)をもつ所望の信号のフーリエ変換を表している。所望の信号と逆フィルタの出力信号はそれぞれＤ(e^j ^ω)＝e^-j ^ωΔＵ(e^j ^ω)とＶ(e^j ^ω)＝Ｈ(e^j ^ω)Ｕ(e^j ^ω)によって逆フィルタへの入力信号（仮想音源信号）と関連しているので、コスト関数はつぎのようにまとめられる。この２次の評価関数を最小にする逆フィルタのフーリエ変換は容易に示すことができる（参考文献[6]の付録を参照）。ここに上付添字＊は共役複数を示す。我々は、この式のＺ変換についてz＝e^j ^ω を代入することで示すことができる。故に|Ｃ(e^j ^ω)|²は、Ｃ^*(e^j ^ω)Ｃ(e^j ^ω)＝Ｃ(e^-j ^ω)Ｃ(e^j ^ω)と表されるので、Ｚ変換はここで、例えば、aが実数の場合において、Ｃ(z)＝1＋az^-1で示される伝達関数のインバージョンを考えてみる。これはz＝−aにおいて単一の零をもち、|a|＞1 （すなわち零点が単位円の外側にある）の時に、最小位相となる。上で定義した評価関数を最小にする最適な逆フィルタはつぎで与えられる。この式の分母を展開して、つぎのように書くことができる。ここで、p1とp2 は逆フィルタの極である。これらは次式によって与えられる。システムの零点がインバートされて単位円の近傍に横たわっている場合は特に興味のあるケースである。このような場合、逆フィルタは非常に長いレスポンスを有する。つまり概括的にいえば、逆フィルタの極は周波数領域における大きな応答と時間領域における長いインパルス応答とに対応しながら単位円に近接するであろう。図解によって、Ｃ(z)の零点がａ＝1＋εにある時を仮定してみると、ここでεは、レギュラライゼーションパラメータβと同じ次数の小さなパラメータであるが、最初の近似の次数定める式(38)の級数展開に従うと、逆フィルタの極はつぎのように与えられる。ここで、次数βとεの項は無視されている。本式は逆フィルタが２つの極を有し、その１つは単位円の内側、もう１つは外側にあることを示す。明らかにβが増加する時、極は単位円から遠くへ移動し、逆フィルタのインパルス応答の継続長はより小さくなる。実際、逆フィルタのＺ変換の表現の部分分数展開は、我々はつぎのように書けることを示している。極であるp1 と p2は式(30)によって与えられる。二項展開により、我々はつぎのように書くことができ、さらに、対応する逆Ｚ変換は数列1,p,p²,p³,…によって与えられる。る。つまり、極が単位円の外側にあることが、インパルス応答を不安定にするのに寄与している。それにもかかわらず、上に強調されるように、前向きの時間におけるこの不安定な応答は後ろ向きの時間においては安定であるという２つの解釈を有する。以下のことによって最も容易にその有効性を確認することができる。つまり、z/(z-p)は(-z/p)/(1-(z/p))としてまた表すことができ、二項式展開を引き続いて用いることで次式を示すことができる。り大きな値はより迅速な減衰をもたらす。レギュラライゼーションパラメータβ の利用により、以下のことが保証される。すなわち、たとえシステムの零点が単位円に非常に近接してインバートされた場合でも、逆フィルタのインパルス応答は十分迅速に減衰する。最後に、式(40)のz^- ^Δがすべてのインパルス応答にΔサンプルの遅延を与えることを記述しておく。それで、もし十分に大きいβの値が選択されるなら、逆向きの時間における逆フィルタの応答はΔサンプルの中で、無視することのできる値に減衰する。これは逆フィルタの因果性を保証する。この逆フィルタのインパルス応答の所要時間の抑制は、フィルタの設計を周波数領域で行う場合に、非常に重要である。技術の基本はＤiscrete Ｆourier Ｔran sform（ＤＦＴ）と、高速フーリエ変換(Ｆast Ｆourier Ｔransform)（ＦＦＴ）アルゴリズムを使った高速な処理をよりどころとしている。適切な変換及び逆変換はつぎのように表される。ここで、数列f(n)がＦ(k)によって与えられる対応するＤＦＴを有する。この際、kは変換が行われる周波数の不連続性を意味するインデックスとして用いられる。最初にインバートされるシステムのインパルス応答c(n)を測定し、次にＦＦＴアルゴリズムを使うことによってＤＦＴＣ(k)を計算する。逆フィルタの周波数応答は次式から算出される。上に明記される逆変換の関係を利用して、対応するインパルス応答は算出される。逆フィルタのインパルス応答が、その周期において、ＤＦＴや逆ＤＦＴの計算の際に用いられるＮサンプルの”基本周期”よりも短いということは、計算のこの段階において、非常に重要となる。もしこのインパルス応答の周期がこの値より大きい場合、その計算は誤った結果ををもたらすであろう。これはもちろん扱われている信号は周期的であるというＤＦＴを使う際に暗黙のうちにたてられる仮定によるものである。実際にこの計算を行うための過程は以下のように要約されることができる。我々は逆フィルタh(n)のフィルタ係数の数をＮ_hで示し、インパルス応答c(n)の所要時間をＮ_cで表す。Ｎ_hは２の乗数でなければならない(2,4,8,16,32,...)。またＮ_hは２Ｎ_cより大きくなければならない。これらは因果性ＦＩＲ逆フィルタh(n)を算出するために必要な手順である。１．c(n)の零付加法(zero-padding)を利用することによって、インバートされる伝達経路におけるインパルス応答の周期はＮ_hサンプルであることが保証される。例えば、もしＮ_c＝256でＮ_h＝1024である場合、768個のゼロが元のインパルス応答c(n)に加えられる。２．ゼロを加えられた数列c(n)のＤＦＴの計算結果はＮ_hの等間隔なポイントの周波数応答Ｃ(k)となる。３．式Ｃ^*(k)/(Ｃ^*(k)Ｃ(k)＋β)からＮ_h周波数の逆フィルタの周波数応答を算出実際は最初を計算することはただ必要なだけである実際はこの式の最初の(Ｎ_h /2)＋1個の値を計算することだけが必要となる。何故なら、実際の数列のＤＦＴは対称であるためである。４．式Ｃ^*(k)/(Ｃ^*(k)Ｃ(k)＋β)の逆ＤＦＴの計算５．逆ＤＦＴの１番目と２番目の半分を交換することによるh(n)の探索例えば、逆ＤＦＴが数列{1,2,3,4,5,6,7,8}である場合、h(n)＝{5,6,7,8,1,2,3,4}となる。この操作（処理）はＮ_h/2＋1のモデリングディレイを実行している。モデリングディレイはそれ故に、逆フィルタのインパルス応答の長さの半分に選択される。多チャンネルケースへのこの技術の拡張は簡単である。まず次式で表されるコスト関数を最小にする逆フィルタの行列を求める。ここでe(e^j ^ω)はベクトル誤差信号（すなわち、所望信号と再生信号の間の差違で定義される信号のベクトル）のフーリエ変換のベクトルであり、v(e^j ^ω)は、逆フィルタ行列の出力信号のフーリエ変換のベクトルである。コスト関数を最小にするような逆フィルタのマトリックスは次式で与えられる。（解析の詳細は参考文献「７」を参照）Ｉが単位マトリックスである場合、所望の信号の行列はΔサンプルで遅延された収音信号の行列であると仮定されている。ここで、レギュラーライゼーションパラメータβは多チャネルケースでも、また単一チャンネルケースでも同一の役割をはたし、これを利用することはインバートされるマトリックスは良好な条件であることを保証し、また、これはこの解の大変に重要な特徴である。逆フィルタ行列を計算するためのステップは次のように要約される。１．インパルス応答の周期はＮ_hサンプルであることを保証するために行われるＣ(n)の要素の零付加法(zero-padding)を利用して、すべての電気音響的な伝送経路のインパルス応答を測定する。２．零付加されたインパルス応答のＤＦＴsを計算する。Ｎ_h等間隔を配置されたポイントにおける周波数応答Ｃ(k)を、結果として得る。３．式[Ｃ^H(k)Ｃ(k)＋βI]^-1Ｃ(k)からＮ_h周波数において逆フィルタの周波数応答を計算する。実際は、この行列でそれぞれの要素の最初の値Ｎ_h/2＋1を訃算することのみが必要である。これは実際のsequcnceのＦＦＴは対称であるという特性に原因を帰する。この式はラウドスピーカのチャンネル数や再生音場での測定の数に関わらず、採用することができる。何故ならβ>0の場合、行列Ｃ^H(k)Ｃ(k )＋βIは正則ではないからである。４．この式の逆ＤＦＴsの行列を計算する。５．逆ＤＦＴの行列のそれぞれの要素の逆ＦＦＴにおいて第１番目と２番目の半分を交換することにより、逆フィルタのインパルス応答を探す。この操作は(Ｎ_h /2)＋1サンプルのモデリング遅延を実行する。６．２つのスピーカシステムに関する主観評価実験６．１これまでの実験に関するレビュー所定の空間に位置している音源のイリュージョンを受聴者に対して作り出すことは、長い間（今まで）音響技術者の目標であった。ラウドスピーカ間のどこかに位置する"幻影（ファントム）"あるいは”仮想”音源により作り出される音像であるイリュージョンを受聴者に対して作り出すために、１対のラウドスピーカに入力される信号を操作す過程に比較的にシンプルな信号処理の技術が利用できることは長い間認識されてきた「8」。これらの技術は通常のステレオフォニックの基礎を形成し、心理音響的な基礎は”音像定位の足し合わせ”というカテゴリ下においてＢlauert「９」によって報告された。すなわち、受聴者に対して適切に配置された１対のラウドスピーカへの２つの入力信号間に単純に振幅差（あるいは時間差）を与えることで、仮想音源のイメージの位置を２つのラウドスピーカ間で移動させることができるというものである。更に洗練された信号処理についてはＡtal とＳchrocder「５」によって報告された。（彼らは一般にはその発明によって知られている）しかしながら、それ以前にＢauer[10]によってダミーヘッド録音の再生という状況のなかで類似した手順が研究されていた。ＡtalとＳchrocderは、”定位ネットワーク(localisation network)"という装置を考案した。それは１対のラウドスピーカに入力される前の信号（これは仮想音源に関連している）を操作するものである。上述されたように、技術の原理は仮想音源信号を１対のフィルタを介して処理することであり、上述されたように、技術の原理は仮想音源信号を１対のフィルタを介して処理することであり、そのフィルタは受聴者の耳元に作り出された信号が、仮想音源の所望の位置に選択されて作り出された信号とが十分に等化であることを保証するために設計された。ＡtalとＳchroederは、提案し採用したフィルタ設計手順において、仮想音源によって受聴者の両耳において作り出された信号は単純に周波数に依存しない利得と時間遅延にのみ関係していると仮定した。そして、この周波数に関して独立な受聴者の両耳での信号間の差違は、音響空間上での仮想音源の位置に依存していると仮定した。これらの仮定は分析的に扱いやすい定位ネットワークのデザインをもたらし、そのパラメータは外見上は異なった仮想音源配置を提供するために変えることができた。この技術の包括的な主観的な評価は、この発明者によっては実施されなかったが、「５」によると、この方法は水平面上に±60度以下の方位角で位置する（つまり、典型的はインテンシティーステレオにて包括される位置角度±30度の範囲の外側「9」）仮想音源のイリュージョンを受聴者に対して作り出すのに十分に効果的であることが報告された。しかしながら、±60度を越えると音像定位はその周波数依存性を増すので、はっきりとは定義されにくくなる。Ｓchrocder他「11」は、後にこの方法をダミーヘッド録音のラウドスピーカ再生に採用した。この場合、ダミーヘッドの耳元で収音された信号は、フィルタネットワークを通して処理される。このフィルタネットワークによって、１対のラウドスピーカにより受聴者の耳元で実質上は同様の信号が再生されることが保証される。このネットワークによって、右のラウドスピーカと左の耳の間、あるいはその逆において、クロストークがキャンセルされることが保証される。そして、主観実験は実施されなかったが、仮想音源は横方向のすべて、または受聴者の後方にさえ、その音像を作り出すことができると報告された。一方、同じタイプのシステム、すなわちクロストークキャンセルのためのプロセスを経て、１対のラウドスピーカから再生されるダミーヘッド録音であるが、これにおける主観評価実験の結果がＤamaske and Ｍellert[12]によって報告された。その中で、彼らは"ＴＲＡＤＩＳ"（ステレオによる全方向情報の再生）技術と名付けた。水平面並びに垂直面の両方向における定位実験の結果は、明らかに技術の効果を証明する。近年、浜田ら「13」は、このアプローチの本質はデジタルでクロストークキャンセルネットワークを利用し、これをＯrthoＳtereophoni c Ｓystem(ＯＳＳ)と呼んだ。まずダミーヘッドの耳元で作り出された信号を収音し、続いて１対のラウドスピーカによる伝送系の前段の2×2のディジタルフィルタマトリックスを介してその信号を処理することによって生成された仮想音源について、主観評価実験はその定位の十分な正確さを示した。Ｎeu他「１４」とＵrbach他「１５」によって、ダミーヘッドの両耳の位置において収音された信号を処理するためのクロストークキャンセルシステムのディジタルでの実行についての一歩進んだ主観評価実験が実施され、その結果が報告された。ここでも良好な結果が見られた。特に水平面における仮想音源の位置についてそれは顕著である。仮想音響信号を作り出すための一般的なアプローチについては、Ｃooper and Ｂauck[16]によって議論された。その中で彼らは ”トランスオーラル・ステレオ”の技術について述べ、また多数の受聴者のための再生システムの一般化についても議論した「17」。トランスオーラル・ステレオについてはＭoller[18]and by Ｋotorynski[19].によっても報告がなされた。これらすべての著者によって利用されたフィルタ設計手順は、ラウドスピーカへの入力から無響室において受聴者の両耳の位置で再生された信号とを関連づける４つの頭部伝達関数(ＨＲＴＦs)の測定あるいは分析的な叙述の形式によるクロストークキャンセルネットワークを構成するフィルタの導出を含んでいる。クロストークキャンセルマトリックスは４つのＨＲＴＦsのマトリックスの逆である。ＡtalとＳchroeder「５」によって認識されたように、もしＨＲＴＦの成分が最小位相でない時、このinversionは非現実的なクロストークキャンセルマトリックスを構築するというリスクを負っている。ＨＲＴＦの成分が最小位相でない場合に（例えば、部屋の表面による屈折などによって「20」）、上述されたフィルタ設計手順を利用することで、このフィルタ問題は処理することができる。これにより、音場再生における問題が非常に一般化した形で体系づけられる（収音と再生信号の多チャンネル化を説明するような）。またクロストークキャンセルフィルタの設計において、最小自乗法を時間領域及び周波数領域のどちらにも適用する。ここでは、仮想音源イメージシステムにおける主観評価実験の結果を示す。このシステムは、水平面上に置かれた（再生された）仮想音源のイリュージョンを受聴者に対して作り出すことができ、さらに様々な音響空間（環境）下においても効果的な再生能力を有することが観測された。しかしながら、上述されたようにＡtalとＳchroederの最初の意図に戻り、つまり信号処理技術を用いて単一の信号が１つの仮想音源に対応するように処理を施し、ダミーヘッド録音をあまりはっきりとは使用しないということである。しかし我々は暗黙にダミーヘッドを用い、ラウドスピーカ入力とダミーヘッドの両耳出力の間の伝達関数の測定によるセットを用いることとする。このダミーヘッドのＨＲＴＦsのデータベースは、空間上の規定された位置にある仮想音源によって、ダミーヘッドの両耳の位置に再生された信号を生成するために仮想音源信号にフィルタをかける際に用いられる。これらの２つの信号はクロストーク相殺用のフィルタにかけらる。このフィルタはイメージが求められる音響空間において同じダミーヘッドの両耳の位置で、これら２つの信号の再生を保証する。無響室、（ＩＥＣ仕様に建てられた）リスニングルーム、並びに車室内にて受聴者に対して行われた主観評価実験の結果を示す。主観評価実験に関する詳細についてはＤ.Ｅnglerの修士論文[21]、及びＦ.Ｏrduna-Ｂustamanteの博士号論文[22]を参照する。上述された信号処理技術の概念は様々な音響空間において仮想の音響イメージを首尾良く生成するための優れた基礎を提供する。６．２．無響条件の下の実験図４は無響室内にて実施される実験のためのクロストークキャンセルフィルタＨ_x (z)を設計するために用いられた音源とダミーヘッドの幾何学的な配置を示す。使用されたラウドスピーカは(ＫＥＦＴype Ｃ35 ＳＰ3093であり、ＫＥＭＡＲＤＢ 4004 artificial head and torsoがダミーヘッドとして使用された。これは明らかにＨＲＴＦデータベースをコンパイルする際に使用されたものと同等のものである。このデータベースは、無響室内にてダミーヘッドから半径2[m]にラウドスピーカを設置し、ラウドスピーカへの入力とダミーヘッドマイクロフォンの出力の間のインパルス応答を測定した。ダミーヘッドの水平面内の円弧上に10度毎にラウドスピーカの位置を移動して、データベースの測定は実施された。インパルス応答は、ＭＬＳＳＡシステムを使用して決定された。このＭＬＳＳＡは線形システムのインパルス応答を決定するためにＭＬＳを適用する[23]。ＨＲＴＦは７２kＨz のサンプルレートにおいて測定され、結果の応答は48[kＨz]にダウンサンプルされた。同様の技術はマトリックスＣ(z)の成分を測定する際にも採用された。このＣ(z)は再生用の２つのラウドスピーカへの入力信号とダミーヘッドマイクロフォンの出力を関連づけるものである。図（Ｆigure）５は行列Ｃ(z)の成分（要素）に対応しているインパルス応答の結果をを示す。図６に時間領域における最小自乗法により前述の手順[1-4] を用いて設計したクロストークキャンセルフィルタＨ_x(z)の成分に対応したインパルス応答を示す。これらのインパルス応答は48kＨzの標本化レートにおいて測定された。最終的に、図７がマトリックスＣ(z)とマトリックスＨ_x(z)を畳み込んだ結果である。これは、十分に効果的ななクロストークの抑圧を示し、積Ｈ_x( z)Ｃ(z)の対角線上の要素のみが重要であり、その良好な近似である式(４)は満足された。選択されたモデリングディレイは150サンプルであったことをここに示す。上に記述されるように、クロストークキャンセルフィルタを設計した後、所望の信号d₁(n)とd₂(n)を対応する選ばれた仮想音源の位置に再生するために、種々の仮想音源信号u(n)を処理するためにＨＲＴＦデータベースは用いられた。これらはクロストークキャンセルフィルタを通って、ラウドスピーカの入力信号となった。クロストークキャンセルフィルタを設計した時にダミーヘッドが設置されていた場所と受聴者の頭部の位置がラウドスピーカの位置に対して、可能な限り同じ位置であるように受聴者は席に着いた。受聴者の周りには音響学上透過なスクリーン（図４参照）が設けられ、その内側に水平面（この平面はラウドスピーカの中心と受聴者の耳の位置を含む）に沿って10度毎の間隔でマーカーがセットされた。受聴者は正面０度に対応するマークを真っ直ぐに向いているように要求され、ラウドスピーカはスクリーンの後方に受聴者に対して対称に左右方位角±30度の位置に配置された。仮想音源の呈示の終了後（すなわち入力信号u(n)と与えられる仮想音源の位置に対応するフィルタＡ₁(z)とＡ₂(z)の選択の組み合わせ）、受聴者は仮想音源の位置についてその角度を決定した。受聴者は、仮想音源の角度の決定を下すまでは正面を向いているように指示され、その後、（もし必要ならば）彼らの選択した仮想音源の位置に最も近いマークをスクリーン上から選び出すために頭部を回転させてもよい。受聴者の頭部の動きは物理的には、制限されなかった。与えられた位置での仮想音源のイリュージョンを作り出すことにおける、このシステムの効果を直接評価するために、実際のラウドスピーカを用いた聴感実験も試みられた。ラウドスピーカは受聴者の周りの半径2[ｍ]の円上の様々な位置に配置された。仮想音源を評価した実験のそれぞれのセットについて、同様の実験が実音源を用いても行われた。各被験者は実音源と仮想音源の両方の試験音のセットを呈示された。試験時間は典型的なセッションの場合で約50分間であり、まず実音源が被験者に呈示された。２日間のインターバルを空けて、仮想音源についての聴感実験が各被験者毎に実施された。実音源と仮想音源として入力された信号u(n)はスピーチ、1/3オクターブ帯域制限されたランダムノイズ（中心周波数250Ｈz,1kＨz,4kＨz）と250Ｈz,1kＨz,4 kＨzの純音で構成された。実施されたすべての実験について表１にまとめる。実音源と仮想音源の両方において、異なった方位角に位置する呈示は３つの角度の”セット”に分けられた。これらは表１に示されている。”Ｓet 0”は前後両方向における角度で構成されているのに対し、”Ｓet 1”と”Ｓet 2”は水平面の前半分の角度のみを含んでいる。表１で明記される各実験で、所定のセットからの角度は特定なシーケンスで呈示された。例えば、シーケンス「ＯＡ」が、"Ｓet 0"からの角度をある特定の順番で呈示することを示している。一方、シーケンス"1Ａ"は"Ｓet 1"からの別の角度の呈示のシーケンスである。実験に用いられた特定のシーケンスを表２に示す。所定のシーケンスでの角度の呈示の順番は、被験者が試験音の順番から学習が出できないようにランダムに選択された。加えて、主観評価において呈示の順番に起因するところのいかなるバイアスをも最小にするために、それぞれのシーケンスは逆からの順番でも呈示された。つまり、シーケンス「1Ａr」はシーケンス「 1Ａ」を逆の順番で呈示した場合を示す。表１で明記される実験のそれぞれについて３人の被験者が実験に参加し、合計12人の被験者がこの実験を受けた。被験者はすべて20代で正常な聴覚を有する。被験者はほぼ男女が半数づつで、３人で構成される被験者のグループに少なくとも１人は女性が含まれるように配慮した。この主観的な実験についての詳細はＤ.Ｅngler[21]に報告されている。システムのパフォーマンスにおいて最も配慮された点は、受聴者の後方に配置された仮想音源のイリュージョンを確かに作り出すことは一般には不可能だということである。これは図８に示される結果によって明らかである。図８は実音源と仮想音源の音像定位の比較について示したものである。図において、正方形の１辺の長さは所定の応答角（回答角）が特定の”呈示角” に対して記録された回数に直接比例している。つまり、被験者がその音源は所定の角度の位置に定位したと回答することで、与えられた試験音にに対して応答した回数である。図８に示されている結果により（これはスピーチ信号による実験結果）、実音源の音像定位はその位置が受聴者の後方にあっても十分に正確であるが、受聴者の後方に作り出された仮想音源に関しては大変しばしばそれらと等しい前方の角度の位置に”反転”される。例えば、150度の呈示角が、受聴者に対しては30度の応答角をもたらすことである。スピーチを実音源として受聴した場合、このような”前後混同”は非常に希であるが、その他の信号が受聴された場合、特に純音の場合（これらの試験のデータについては参考文献[21]を参照することを勧める）には、この前後誤判定は顕著に見られる。このシステムが受聴者の前方に仮想音源の確かなイリュージョンを作り出す能力を有することが図９よりいっそう明らかである。これは、とりわけ±60°の範囲で顕著である。しかし、時折被験者はこの角度の範囲でも前後方向に関して混乱している。±60°の外側の範囲で、被験者は音像を呈示された角度よりやや前方に音像定位する傾向にある。(すなわち、90度の呈示角の音像は60°,70°また 8 0°に知覚される）。これは中心周波数をそれぞれ250Ｈz, 1kＨz 4kＨz とした1/3オクターブで帯域制限された白色雑音を音源とした際の実験結果より、はっきりと観測できる。時折、前後方向の誤判定が見られるが、このデータは基本的にシステムの効果にはある種の周波数依存性があることを示している。4kＨｚのデータを見る[21]と、仮想音源のイメージは受聴者の横方向の本来意図される位置より前方に音源が定位されている時に、大きな程度でイメージが前方に知覚されているのが見て取れる。純音による実験結果[21]は、1/3 のオクターブバンドの白色雑音の結果と比較してかなりデータ間にばらつきが見られるが、それでも前後判定に関して類似した傾向を示した。６．３リスニングルームにおける実験ＩＥＣ規格に基づいて建てられたリスニングルームにおいて、無響室にて実施されたのと同様の評価実験を行った。ラウドスピーカ、受聴者とスクリーンの幾何学的な配置は図４に示されたそれと全く同一であった。しかしながら、インバートされる電気音響系の応答は明らかに異なり、その応答を図10に示す。図５と比較すると、ラウドスピーカに入力される信号はリスニングルームの表面（壁、床等）に起因して、ダミーヘッドの両耳に十分に強い反射を連続してもたらす。図11がクロストークキャンセルフィルタ（ここで再び、時間領域において最小自乗法により設計された[1-4]）のインパルス応答であり、図12は図10に示す測定されたインパルス応答との畳込みの結果である。再び、フィルタ設計手順はシステムを逆畳み込みして、十分なネットワークの応答を積Ｃ(z )Ｈx(z)の対角線上の項のみから得ることができるという点で、非常に効果的である。上述したように無響条件の下で行われた実験と同一の手順で実施された。表１に挙げられたすべての試験（表２で分類されたシーケンスを用いて）はリスニングルーム内において繰り返された。しかしながら、12人の被験者は異なったグループに分けられた。しかしながら、12人の被験者は異なったグループに分けられたが、実音源と仮想音源の定位に関してあくまでも同一の実験手順による主観評価実験が実行された。再び、被験者は20代の正常な聴覚を有する男女であり、男女同人数である。図13は仮想音源イメージシステムの有効性と、受聴者が実際のスピーチ音源を音像定位する能力を比較した結果を示す。再び、システムは受聴者の後方に確かなイメージを作り出すことができない、ほとんどすべての仮想音源再生においてそのイメージはそれらの”鏡像の位置”である前方に知覚される。図13はスピーチ音源による実験結果であり、結果がここには示されてないが、その他の種類の信号（純音や1/3オクターブバンド雑音）を実音源とした時の音像定位はスピーチ音源に比べて正確さに欠け、十分な数の前後誤判定が見られる[21]。ここで再び、システムは受聴者の前方、特に±60°の範囲において、正確に定位されるイメージを再生するのに非常に効果的であった。図14はこの実験結果であり、無響室内にて実施された同様の実験で観測されたもの（図9）に比べて前後の誤判定が少ないことを示している。図14はまた、このシステムが受聴者の両側において仮想音源の前方のイメージを生成する傾向にあることも示している。この傾向は1/3のオクターブバンド雑音を音源とした実験によっても観測された。特に4kＨzの中心周波数を有した音源の場合にこの傾向は顕著であった。また、2 50Ｈzの音源による実験結果から、この周波数においてはリスニングルームで受聴した場合の方が無響室内での実験に比べてその被験者の応答にばらつきがあることも興味深い挙動である。Ｅngler[21]による追加のデータは、反響（reverbe rant）な音響空間下での純音の受聴はその定位においては一般に大変に不十分であることを示す。1kＨzと4kＨzによる実験結果から、そのばらつきは無響室内で実行された実験の結果と類似しているが、250Ｈzの場合にはそのばらつきはの程度は無響室内にくらべて非常に大きい。６．４車室内の実験最後により冒険的な試みとして、確かな音源を再生するためのシステムの能力について検討するために、いくつかの簡易な実験を自動車の車室内にて行った。車は左ハンドルの”イスズＩ-ＭarkＸＳ”であり、既存ののオーディオシステムのラウドスピーカを受聴者に信号を呈示する際に使用した。これらは自動車のダッシュボードの外側に下向きに水平面に対して45度の角度を有して取り付けられ。既存ののオーディオシステムのラウドスピーカを受聴者に信号を呈示する際に使用した。これらは自動車のダッシュボードの外側に下向きに水平面に対して45度の角度を有して取り付けられ。図15におおよその実験のレイアウトを示す。ラウドスピーカは受聴者の両耳を含む水平面の下に上手く配置された。ダミーヘッドがクロストークキャンセルフィルタの設計に用いたダミーヘッド、および受聴者は自動車の左側にあるドライバーシートの同じ位置に配置された。ラウドスピーカと自動車の車室内の間のインパルス応答は、クロストークキャンセルフィルタの設計は使用できるフィルタ係数の数に限界があるために困難であり、満足に逆システムを設計することも困難であることを証明した。自動車後方のトランクに無響材質のくさびを挿入し、車内のインテリアを覆うことでこの状況を緩和するためのいくつかの試みを実行した。電気音響系の伝達関数のマトリックスを含んだインパルス応答は図16に示される。これらのインパルス応答の形状や周期は、明らかに無響室やリスニングルームにおけるものとは異なっており、インパルス応答の中に直接波の後に到来する本質的なエネルギーが存在する。これはもちろん、受聴者に近接して非常に強い反射をもつ自動車のインテリアの表面の特性によるものである。クロストーク相殺フィルタは同様に非常に長い所要時間を有し、これらのインパルス応答は図17で示される。これらは時間領域技術[1-4]を使うことによって再び設計された。インパルス応答のトランケーションにより、上述されたケースに比べてインバージョンの効果は低減する。これは逆畳み込みされたシステムの伝達関数の周波数分析によってより詳細に確認することができる。逆畳み込みされたシステムの対応するインパルス応答を図18に示す。これより、クロストークの相殺はその困難性にも関わらず、基本的には効果的に行われていることが確認できる。ここで考慮されている環境では実音源と仮想音源の両方について実験を行い、比較することは不可能なので、実験は仮想音源のみにつちえ実施された。上に記述された実験により、スピーチ信号を仮想音源とした時に、このシステムは最も効果的であったので、スピーチ信号のみを実験で使用した。本質的には、前述された実験と同様な手順に従って、これらの実験は実施され、被験者は直接正面を向くように指示され、仮想音源の位置の角度を決定し、そして車の外側の水平面上に設置されたマーカから対応する角度ｗを選択した。位置角度の判定に加えて、被験者は仮想音源が水平面に対して”上方”、”下方 ”、あるいは”同じ高さ”のいずれであるか、その”エレベーション”についても回答した。この単純な実験はは、前の実験と異なり、信号を再生するために使用されたラウドスピーカが水平面からかなり下がった位置に設置されていたことを考慮して加えられた。受聴者の両耳における所望の信号は水平面内の仮想音源であった。合計12人の正常な聴覚をもった被験者がこの実験に参加した。これらの被験者は無響室やリスニングルームで行われた実験とは異なった。合計で38個のランダムに抽出された位置角度を有する仮想音源がそれぞれの被験者に呈示された。その定位実験の結果を図19に示す。スピーチ信号を用い、先だって行われた２つの条件下での実験に比べてデータのばらつきは大きいが、大変類似した傾向もデータから見られた。例えば、中央に位置したイメージは確実に中央に定位するが、受聴者の横方向に置かれた仮想音源は”前方イメージ”となる傾向がある。また、前方イメージとは矛盾するいくつかの結果も見られる。つまり、比較的多数の実験において、受聴者の横方向に（60〜90度と-60〜-90の間に）呈示された仮想音源はすべて正確に90度と-90度の位置に定位した。これらの結果は実際は” 前後の混同”から推定され、自動車の外側に設置されたマーカーの列から選択されうる位置角度の一番端（つまり±90°度）に位置することが可能である。上昇テストの結果より、回答不能だったデータも無視はできないが、平均して被験者は仮想音源が水平面上にあると判定したことがわかる。十分な数の被験者が彼らの左側に定位した仮想音源について、水平面の下に位置すると回答したことは、比較的広い角度で被験者の下に置かれた左側のラウドスピーカについて考慮するとあまり驚くべき結果ではない。振り返ってみて、この上昇テストは、被験者が正中面方向のレンジを利用しながら仮装音源のエレベーションについて回答すれば、よりよい検討が可能であったと思われる。このデータから明確なことは、しかしながら、被験者が首尾一貫して仮想音源が水平面の下の位置をにあると判断（回答）しなかったことである。６．５議論上記の実験の結果は、記述された信号処理技術は受聴者の前方の水平面上に±60 °以上の範囲で仮想音源のイメージを確実に作り出すのに大変効果的な手段であることを明示する。さらに、これは様々な環境下（Ａcoustic property）において、ほとんどそれら音響環境の複合性に係わりなく達成され得る。この技術は正常な聴覚を有する被験者に対して一貫して効果的であることが証明され、クロストークキャンセルフィルタはその音場には依存するが、各受聴者それぞれについて設計されないといったことは強調されるべきである。ある特定の傾向はデータを通して、繰り返し見られた。例えば、システムは被験者の後方に仮想音源を作り出すことはできず、これらの呈示は一般にはそれら音像の鏡像が受聴者の前方に知覚された。受聴者の横方向に呈示された仮想音源は、そのイメージか前方によるという傾向があり、それらは水平面上の本来意図される位置よりも前方に知覚される。従って、受聴者に対して横方向に仮想音源を確実に定位させることは困難であるが、しかしそれでも±60°の外側の範囲に位置する仮想音源を生成することは可能であることが見いだされた。７．多チャンネルラウドスピーカを使った水平面の後方半分への音像の生成上に記述された２チャネルの仮想音源イメージシステムは多数の受聴者に対してその前方にイメージを生成するのに非常に効果的なシステムであり、受聴者の横、ならびに後方にイメージを作り出すための可能性について開発を進めることは明らかに重要である。参考文献[11-15]で報告されているように、過去のいくつかの実験のように、２つのラウドスピーカで受聴者の前方にそれら仮想音源を再生することは可能なことである。しかしながら、それら過去の研究は無響室内で実施され、音源としてダミーヘッド録音によるソースを使用していた。クロストークキャンセルフィルタの設計において十分な配慮と注意を払うことにより、任意の環境下において、２つのラウドスピーカのみを用いて同様の効果を得られる可能性が高い。これは、各受聴者のＨＲＴＦについてその詳細を考慮に入れるというように、個別の基準の上に達成されなければらならい。例えば、受聴者の両耳（鼓膜）での２つの信号を正確にヘッドフォンを用いて呈示することによって、任意の位置にイメージを作り出すことは可能であり、そのイメージは受聴者の後方や、あるいはその上方にさえ再生され得るであろう[24,25]。しかしながらその研究の中で、正確な信号を呈示することを保証するためには、（外耳道の応答も含んだ）各個人それぞれののＨＲＴＦを用いて逆システムを設計しなければならなかった。これはまた、ヘッドフォン再生において、受聴者は”頭の内側”に作り出されたイメージを知覚することがよく知られているために、このシステムによるイメージが受聴者の頭の外側に作り出されることを保証する必要があった。最終的に、側面と後方へイメージを作り出すための、これまでの方法は一般に頭部の回転に非常に敏感である。ここで記述されるシステムについて実験に基づいた詳細な研究は行われなかったが、このシステムによって生成されたイメージは頭部の回転に対して比較的敏感でないことが見い出されており、仮に頭部が大きく回転されて（例えば60度）そのイメージか破壊された場合も、受聴者がその頭部を本来の位置に戻すと、すぐにそのイメージもそれらがあるべき位置に定位するであろう。さらに、少しの回転に対しては（例えば30度以下）、イメージの位置は非常に安定していることがわかった。従って、多数の受聴者に対して、そのの側面と後方に確実なイメージを作り出し、しかも受聴者の頭部の回転に対して敏感でないようなシステムを構築することは、非常に興味深い。また、様々な環境下(acoustic environment)でこれを達成することが可能であることは同様に重要である。それを実現するために、ここで、上述された２チャンネルでの技術の多チャンネル一般化に基づいた新たな方法についてその概要を述べる。技術の本質は、受聴者の後方に追加のラウドスピーカを設置し、セクション４で記述された方法にて設計された逆フィルタ列を介して、仮想音源信号を処理することである。これは、２つの音源を用いたシステムに対しては十分に有効であったフィルタ設計の技術が、多チャンネル、または任意の数のラウドスピーカによるシステムに対してどのように一般化されることができるかを示す。したがって、再生に用いられるいかなる個数の音源においても、クロストークキャンセルマトリックスを構築し、これと仮想音源により作り出されるであろう信号を特定するインパルス応答のマトリックスを畳み込むことができる。明らかに、より多数のチャネルが選択されれば、再生音場においてより多くのポイントで所望の信号は複製される。概括的に話をして、それ故に、もし多くのラウドスピーカが再生のために用いられれば、受聴者へのイリュージョンはより確実となることが期待される。しかしながら、この挑戦は最小の個数のラウドスピーカで受聴者の側面と後方に仮想音源の確実なイリュージョンを作り出そうとするものである。これを実用化するにあたって、都合が良く効率的な手段は、つまり２つのラウドスピーカを前方に、後２つを後方に設置する方法であることが見出された。この構成はＱuadraphonic音場再生を用いることと同様の効果をもたらす。しかしながら、このシステムにおける信号処理に関しては、Ｑuadraphonic音場再生におけるそれらの本来の手法とは全く異なっていることをここに強調する。初期におけるこれらのシステムは、信号は従来のステレオ再生に使われる通常の音源とする通常の方式と同様の一般的な方式を用いることにより、確実なイメージを生成することができなかったため、一般的な受容を獲得できなかった。つまり、音源となる信号はしばしば、イメージの所望の位置に依存するように異なった振幅を有した信号をラウドスピーカに与えてやることで容易に決定される。Ｂla uertはいくつかの簡単な信号処理手法について議論している[9]。従って、ここで記述される信号処理手法は、従来のステレオ方式を改善するのと同様の方法で、クワドラフォニック(quadraphonic)システムの上に改善を施されるものである。すなわち、信号処理は受聴者の両耳の領域で、所望の仮想音源の振幅と位相の両方を調整する。さらに、ここで記述された特別なアプローチでは、受聴者の両耳の領域での再生音場における音の到来方向と、所望の音場での仮想音源のそれとをぴったりと合わせることを保証することが可能である。これは我々が、その位置において仮想音源信号の非常に正確な再生を保証しようと試みているポイントを、受聴者の頭部の領域で注意深く選択することで達成される。それらのポイントとは、つまり、クロストークキャンセルフィルタを設計する際に用いる応答が測定されるポイントである。明らかに、受聴者頭部の領域においてそのポイントに関して多くの組み合わせが存在するが、受聴者の頭部に対して片側に非常に近接した２点を、反対側にさらに２点をこれも近接して選ぶ方法が、最も効果的であることが理解された。これは図20に示される。参考文献[7]で記述された研究では、音場を収音するために過去に用いられていた音の波長に比較してコンパクトなマイクロフォンの集合を囲むようにいくつかのラウドスピーカが用いられ、収音された信号が最適な逆フィルタによって処理される場合、原音場における指向特性はそのマイクロフォン群の領域において良好に再生されることが報告されている。この原理はここでは、受聴者の後方へ仮想音源をシミュレーションすることは、受聴者の後方に設置されたラウドスピーカに寄与するところが非常に大きいということを保証するために用いられる。同様に、もし受聴者の前方へ仮想音源を再生することが要求されるなら、信号処理の手法は、この音場再生は受聴者の前方に設置されたラウドスピーカによって支配されているということを保証する。このアプローチを最も都合のよい状態で実行するには、上述された２チャンネルの場合と類似した方法で、ダミーヘッドを用いるのがよい。しかしながら、この場合、受聴者の頭部が意図される位置に場合のみならず、わずかに回転している状況も考慮したうえで、ダミーヘッドマイクロフォンの位置において、非常に正確な再生を保証するためにクロストークキャンセルフィルタは設計される。４つのラウドスピーカの入力信号と受聴者の頭部領域での４つの位置の関係であることろの4×4のマトリックスＣ(z)を定義するために、合計で４つの測定の位置が与えられる。4×4のクロストークキャンセルフィルタ行列Ｈx(z)は上述の式(24) を満足するように設計される。これは再び参考文献[1-4]で記述された時間領域での技術、あるいはセクション５で記述された周波数領域での技術によって成し遂げられる。明らかに、この原理は多数のラウドスピーカを用いた再生方式や、またダミーヘッドの表面上、あるいは近傍に追加されたマイクロフォンの利用に拡張されることが可能である。さらに、多チャンネルのクロストークキャンセルフィルタ行列を設計するためには、ダミーヘッドマイクロフォンの利用は必ずしも必要ではなく、回転楕円な散乱物体の表面上に取り付けられたマイクロフォンアレイの利用も可能であろう。もちろん、各受聴者の頭部の近傍に配置されたマイクロフォンを利用することも可能であろう。最後に、受聴者の頭部の領域において所望の仮想音源の信号を特定し、クロストークキャンセルフィルタ行列を設計するために、解析的、数学的、あるいは経験的なＨＲＴＦのモデルを用いて、逆フィルタを設計することは有用であろう。しかし、４つのラウドスピーカによる仮想音源再生システムを設計する上で、ダミーヘッドを回転させるという単純な技術は、実用上大変効果的であることが証明されている。図20に示すようにダミーヘッドを±5度ずつ回転させて行われた測定（つまり、まずダミーヘッドを両耳から+5度回転した状態で測定は実施され、次に-5度に回転して測定された。）の結果を用いて設計されたシステムに関する聴感実験が、無響室内にて実施された。この場合、逆システムを構成するための行列Ｃ(z)を定義するために過去の測定で得られたＨＲＴＦのデータベースを利用して（つまり、再生における実際の電気音響系はインバートされていない）、クロストークキャンセルフィルタは設計された。これらの実験により、上述された２チャンネル再生システムにおいて観測された前後誤判定の問題はこのシステムでは軽減されることが結果より明らかになった。図21はスピーチ音源による実験結果である。仮想音源のイメージが確実に受聴者の側面、および後方に受聴者の側面、および後方に再生されていることが見て取れる。頭部がより大きな角度で回転された場合（例えば±15度）についても実験が行われたが、あまり良好な結果は得られなかった。受聴者の頭部近傍の片側に配置されたマイクロフォン間の距離は、対象とする音のうちで最も高い周波数をもつ音の波長より短くなければならないことに起因すると思われる。このような間隔を置くことは、仮想音源の到来方向を正確に再生することを保証する[7]であろう。ダミーヘッドが10度回転する場合、マイクロフォン間の距離は直線で 1[cm]離れていなければならないことをここに記しておく。これは、16[kＨz]までの周波数帯域において確実な再生が得られることを明示している。回転角度が小さくなれば、この帯域は広がると思われる。より大きな頭部の回転に関しては、十分な周波数帯域での確実な再生は期待できない。最後に、ここで述べたきた実験では、セクション５で記述した周波数領域での逆システム設計技術使われたことを記しておく。式（47）を用いてクロストークキャンセルフィルタを算出する際に、パラメータβを用いてインバージョンを一般化することが特に重要であることが解った。βの値はカット＆トライによって抽出されるが、伴うコンピューテーションの速度と効率の点において考慮する際には、それほど困難ではない。反復フィルタ設計手順を設けることによってパラメータβの選択を自動化することも可能である。

【手続補正書】特許法第１８４条の８第１項【提出日】１９９６年７月８日【補正内容】請求の範囲１．複数のラウドスピーカによる再生のための収音方法、もしくは複数のラウドスピーカによる再生のための音声処理に関する方法を提案し、そのなかで再生されるいくつかの音は受聴者に対して仮想音源から発生されているがごとく現れ、ラウドスピーカに供給するための信号を録音したり処理するためにフィルタ手段 (Ｈ)を用いて成り立っており、そのフィルタ手段(Ｈ)は以下のような２つのフィルタ設計の手順で特徴づけられる、すなわち： a）仮想音源の前記所望の位置と再生音場における特定位置、すなわち受聴者の耳元、あるいは頭部の領域、の間の伝達関数であるフィルタ(Ａ)を明確に記述し、 b）前記伝達関数フィルタ(Ａ)と、スピーカの入力からある特定の地点までの伝達系(Ｃ)、あるいは電気音響的な伝送系の逆システムを得るためのクロストークキャンセルフィルタ(Ｈ_x)の畳み込みを行う。２．クロストーク相殺用のフィルタマトリクス(Ｈ_x)は、多チャンネルシステムのためのフィルタ設計手順によって、以下のように構成されるという特許請求の範囲第１項記載の方法。 a）すべての電気音響系の送信経路のインパルス応答c(n)を測定し、インパルス応答のタップ長がＮ_h点であることを保証するために要素c(n)に零値を付加し、 b）Ｎ_h点で均等に間隔をおかれた周波数応答Ｃ(k)を得るために、ゼロパディングされたインパルス応答のＤＦＴsを算出し、 c）式[Ｃ^H(k)Ｃ(k)＋βI]^-1Ｃ(k)から、周波数Ｎ_hでのフィルタの周波数応答を算出し、ここで"Ｈ"はその行列の転置、"Ｉ"は単位行列、βはレギュラライゼーションパラメータであり、 d）前記の式の逆ＤＦＴs行列を算出し、 e）逆ＤＦＴ行列の各要素の逆ＦＦＴの最初と２番目の半分を交換することによって、フィルタのインパルス応答を算出し、(Ｎ_h/2)+1 サンプルでモデリングディレイを実行する。３．インパルス応答h(n)を有する単一チャンネルの逆フィルタを設計するために、最小自乗法が周波数領域において採用され、最小自乗技術は以下のように構成されるフィルタ設計手順を用いるとする特許請求の範囲第２項記載の方法： a）フィルタh(n)におけるフィルタ係数の個数を表すためにＮ_hを用い、インパルス応答c(n)のタップ長をＮ_cとし、ここで、Ｎ_hは２の階乗(2,4,8,16,32 .... ) であり、Ｎ_hは２Ｎ_cより大きい値でなければならず、 b）逆システムとなる伝送経路のインパルス応答のタップ長がＮ_hサンプルであることを保証するためにc(n)に零値を加える手段が用いられ、 c）Ｎ_h点の均等な間隔で周波数応答Ｃ(k)を与えるため、零値を付加された数列 c(n)のＤＦＴ(Ｄiscrete Ｆourier Ｔransform)の算出し、 d）式Ｃ^*(k)/(Ｃ^*(k)Ｃ(k)＋β)から、周波数Ｎ_hにおけるフィルタの周波数応答を算出し、 e）式Ｃ^*(k)/(Ｃ^*(k)Ｃ(k)＋β)の逆ＤＦＴを算出し、ここでβはレギュラライジングパラメータであり、 f）逆ＤＦＴの最初と２番目の半分を交換することによって h(n)を計算する。４．伝達関数ＡやＣが、実音源への入力とダミーヘッドの耳元（あるいは頭部の領域）におけるマイクロフォンからの出力とのまず第一に実施される測定によって導出され、このダミーヘッドは受聴者の頭部伝達関数(ＨＲＴＦ)の影響をモデル化するために利用されれという特許請求の範囲第１、第２項、または第３項記載の方法。５．最小自乗法が、意図される受聴者の位置における再生信号(w)と所望の信号( d)との時間平均誤差を最小にするように採用されるとする特許請求の範囲第１項から第４項の何れか１項に記載の方法。６．伝達関数が、実際の受聴者によってまず第一に実施される測定によって導出されるとする特許請求の範囲第１項から第５項の何れか１項に記載の方法。７．伝達関数が、または解析上、あるいは経験的に得られる頭部伝達関数(ＨＲＴＦ)のモデルによって導出されるとする特許請求の範囲第１項から第５項の何れか１項に記載の方法。８．フィルタ設計の手段が、仮想音源を受聴者耳軸上前方に位置するように配置させることを特徴とする２つのスピーカのみを用いた特許請求の範囲第１項から第７項の何れか１項に記載の方法。９．２つのスピーカが受聴者正面に設置され、また少なくとも１つのスピーカが受聴者の後方に用いられるとする特許請求の範囲第１項から第７項の何れか１項に記載の方法。 10．２つのスピーカが受聴者後方にあるという特許請求の範囲第９項記載の方法。 11．伝達関数フィルタ設計の手順が、仮想音源の所望の位置と受聴者の両耳に近接した４つの特定の位置、その２つは片耳に、もう２つは別の耳に近接した、との間の伝達特性を決定することから成るとする特許請求の範囲第９項および第10 項記載の方法。 12．耳の位置にマイクロホンが装着されているダミーヘッドが、前記伝達関数の測定に用いられ、各耳に近接した前記２つの位置を供給するためにダミーヘッドが微少な角度で回転され、４つのラウドスピーカの入力信号と受聴者の頭部の領域の４つの位置を関連させる4×4のマトリックスＣ(z)が決定されるとする特許請求の範囲11項記載の方法。 13．多チャンネル音場再生システムを介して演奏することによって、再生されることが可能となる、多チャンネル音場収音を生み出し、フィルタ設計手順を利用する特許請求の範囲第１項から第12項の何れか１項に記載の方法。 14．複数のラウドスピーカとラウドスピーカに入力される事前に収音された信号を処理するためのフィルタ手段から成り、前記フィルタ手段が、フィルタ設計手順をふまえた方法により構築されるとする特許請求の範囲第１項から第12項の何れか１項に記載の手順であるとする音場再生システム。 15．再生のためのオーディオシステムを備えた自動車、特許請求の範囲第１項から第13項、または第15項から第17項の何れか１項に記載の方法によって構築されるフィルタ手段を用い、ラウドスピーカを使用する前記システム。

───────────────────────────────────────────────────── フロントページの続き (81)指定国ＥＰ(ＡＴ，ＢＥ，ＣＨ，ＤＥ，ＤＫ，ＥＳ，ＦＲ，ＧＢ，ＧＲ，ＩＥ，ＩＴ，ＬＵ，ＭＣ，ＮＬ，ＰＴ，ＳＥ)，ＯＡ(ＢＦ，ＢＪ，ＣＦ，ＣＧ，ＣＩ，ＣＭ，ＧＡ，ＧＮ，ＭＬ，ＭＲ，ＮＥ，ＳＮ，ＴＤ，ＴＧ)，ＡＰ(ＫＥ，ＭＷ，ＳＤ，ＳＺ，ＵＧ)，ＡＭ，ＡＴ，ＡＵ，ＢＢ，ＢＧ，ＢＲ，ＢＹ，ＣＡ，ＣＨ，ＣＮ，ＣＺ，ＤＥ，ＤＫ，ＥＥ，ＥＳ，ＦＩ，ＧＢ，ＧＥ，ＨＵ，ＩＳ，ＪＰ，ＫＥ，ＫＧ，ＫＰ，ＫＲ，ＫＺ，ＬＫ，ＬＲ，ＬＴ，ＬＵ，ＬＶ，ＭＤ，ＭＧ，ＭＫ，ＭＮ，ＭＷ，ＭＸ，ＮＯ，ＮＺ，ＰＬ，ＰＴ，ＲＯ，ＲＵ，ＳＤ，ＳＥ，ＳＧ，ＳＩ，ＳＫ，ＴＪ，ＴＭ，ＴＴ，ＵＡ，ＵＧ，ＵＳ，ＵＺ，ＶＮ (72)発明者濱田晴夫東京都千代田区神田錦町２−２東京電機大学情報通信工学科音響研究室内 (72)発明者ネルソンフィリップアーサーイギリス国エスオー２３エヌワイハンプシャーサウザンプトンバセットチェトウィンドドライブ 19

Claims

【特許請求の範囲】１．複数のラウドスピーカによる再生のための収音方法、もしくは複数のラウドスピーカによる再生のための音声処理に関する方法を提案し、そのなかで再生されるいくつかの音は受聴者に対して仮想音源から発生されているがごとく現れ、ラウドスピーカに供給するための信号を録音したり処理するためにフィルタ手段 (Ｈ)を用いて成り立っており、そのフィルタ手段(Ｈ)はフィルタ設計の過程で作られるものであり、フィルタ設計過程とは以下のように特徴づけられる： (a)収音したものをラウドスピーカを介して再生する際、意図される位置の受聴者に対して再生される信号(w)と、意図される位置における所望の信号(d)との間の誤差を最小にするように用いられる技術 (b)受聴者に対して再生される前記所望の信号(d)は、仮想音源の所望の位置での音源によって前記の意図される位置の受聴者の耳元(あるいは受聴者の領域)において作り出されるであろう信号（あるいは推定信号）で定義される。２．所望信号が、仮想音源の前記所望の位置と再生音場における特定位置、すなわち受聴者の耳元、あるいは頭部の領域、の間の伝達関数であるフィルタ(Ａ)の形式で明確に記述されることが望ましいことを特徴とする特許請求の範囲第１項記載の方法。３．伝達関数が、実音源への入力とダミーヘッドの耳元(あるいは頭部の領域)におけるマイクロフォンからの出力とのまず第一に実施される測定によって導出され、このダミーヘッドは受聴者の頭部伝達関数の影響をモデル化するために利用されることを特徴とする特許請求の範囲第２項記載の方法。４．最小自乗法が、意図される受聴者の位置における再生信号(w)と所望の信号( d)との時間平均誤差を最小にするように採用されることを特徴とする特許請求の範囲前記何れか１項に記載の方法。５．伝達関数が、実際の受聴者による測定によって導出されることを特徴とする特許請求の範囲第２項記載の方法。６．伝達関数が、または解析上、あるいは経験的に得られる頭部伝達関数（ＨＲＴＦ）のモデルによって導出されることを特徴とする特許請求の範囲第２項記載の方法。７．フィルタが、‘クロストークキャンセル’のマトリックス(Ｈ_x)と共に所望の信号を特定する伝達関数を表しているディジタルフィルタ(Ａ)の畳み込みによって導出されることを特徴とする特許請求の範囲第２項から第６項何れか１項に記載の方法。８．フィルタ設計の手段が、仮想音源を受聴者耳軸上前方に位置するように配置させることを特徴とする２つのスピーカのみを用いた特許請求の範囲第１項から第７項の何れか１項に記載の方法。９．２つのスピーカが受聴者正面に設置され、また少なくとも１つのスピーカが受聴者の後方に用いられることとする特許請求の範囲第１項から第７項の何れか１項に記載の方法。 10．２つのスピーカが受聴者後方にあることとする特許請求の範囲第９項記載の方法。 11．伝達関数フィルタ設計の手順が、仮想音源の所望の位置と受聴者の両耳に近接した４つの特定の位置、その２つは片耳に、もう２つは別の耳に近接した、との間の伝達特性を決定することから成るとする特許請求の範囲第９項および第10 項記載の方法。 12．耳の位置にマイクロホンが装着されているダミーヘッドが、前記伝達関数の測定に用いられ、各耳に近接した前記２つの位置を供給するためにダミーヘッドが微少な角度で回転され、４つのラウドスピーカの入力信号と受聴者の頭部の領域の４つの位置を関連させる₄ _×4のマトリックスＣ(z)が決定されるとする特許請求の範囲11項記載の方法。 13．多チャンネル音場再生システムを介して演奏することによって、再生されることが可能となる、多チャンネル音場収音を生み出し、フィルタ設計手順を利用する特許請求の範囲第１項から第12項の何れか１項に記載の方法。 14．複数のラウドスピーカとラウドスピーカに入力される事前に収音された信号を処理するためのフィルタ手段から成り、前記フィルタ手段が、フィルタ設計手順をふまえた方法により構築されるとする特許請求の範囲第１項から第12項の何れか１項に記載の手順であるとする音場再生システム。 15．最小自乗法が、周波数領域でフィルタを設計するために、周波数領域において採用されるとする特許請求の範囲第１項記載の方法。 16．フィルタ手段が、以下のような手順で構成されているとする特許請求の範囲第15項記載の方法： a）フィルタh(n)におけるフィルタ係数の個数を表すためにＮ_hを用い、インパルス応答c(n)のタップ長をＮ_cとし、ここで、Ｎ_hは２の階乗(2,4,8,16,32 ....) であり、Ｎ_hは２Ｎ_cより大きい値でなければならず、 b）逆システムとなる伝送経路のインパルス応答のタップ長がＮ_hサンプルであることを保証するためにc(n)に零値を加える手段が用いられ、 c）Ｎ_h点の均等な間隔で周波数応答Ｃ(k)を与えるため、零値を付加された数列c (n)のＤＦＴ(Ｄiscrete Ｆourier Ｔransform)の算出し、 d）式Ｃ^*(k)/(Ｃ^*(k)Ｃ(k)＋β)から、周波数Ｎ_hにおけるフィルタの周波数応答を算出し、 e）式Ｃ^*(k)/(Ｃ^*(k)Ｃ(k)＋β)の逆ＤＦＴを算出し、ここでβはレギュラライジングパラメータであり、 f）逆ＤＦＴの最初と２番目の半分を交換することによって h(n)を計算する。 17．フィルタマトリックスを利用し、前記マトリックスは以下のようなフィルタ設計手順で構成されているとする特許請求の範囲第16項記載の方法： a）すべての電気音響系の伝送経路のインパルス応答を測定し、インパルス応答のタップ長がＮ_h点であることを保証するために要素c(n)にゼロパディングを導入し、 b）Ｎ_h点の均等な間隔で周波数応答Ｃ(k)を与えるため、零値を付加されたインパルス応答のＤＦＴを算出し、 c）式[Ｃ^H(k)Ｃ(k)＋βI]^-1Ｃ(k)から、周波数Ｎ_hでのフィルタの周波数応答を算出し、 d）前記の式の逆ＤＦＴsマトリックスを算出し、 e）逆ＤＦＴ行列の各要素の逆ＦＦＴの最初と２番目の半分を交換することによって、フィルタのインパルス応答を算出し、(Ｎ_h/2)+1 サンプルでモデリングディレイを実行する。 18．再生のためのオーディオシステムを備えた自動車、特許請求の範囲第１項から第13項、または第15項から第17項の何れか１項に記載の方法によって構築されるフィルタ手段を用い、ラウドスピーカを使用する前記システム。