(57)【要約】【課題】自動車内の対称軸外に配列された聴取位置における両耳間の時間差を最小化することを可能にする方法を提供すること。【解決手段】聴取室内の少なくとも１つの聴取位置（１０）において音響ローカライゼーションを最適化する方法であって、該方法は、該第２のラウドスピーカ（４）に供給する少なくとも該オーディオチャネルのための位相イコライゼーションフィルタのフィルタ係数を計算することであって、該位相イコライゼーションフィルタの位相応答は、該少なくとも１つの聴取位置（１０）における両耳の位相差（Δφｍｎ）、または２つ以上の聴取位置（１０、１１）の上に平均される平均的な両耳の位相差（ｍΔφｍｎ）が所定の周波数範囲内で最小化されるように設計される、ことと、該位相イコライゼーションフィルタを該それぞれのオーディオチャネルに適用することとを包含する、方法。【選択図】図３
2011097561.tif 000002
【特許請求の範囲】
【請求項１】
聴取室内の少なくとも１つの聴取位置（１０）において音響ローカライゼーションを最適化する方法であって、音場は、該少なくとも１つの聴取位置（１０、１１）に割り当てられた一群のラウドスピーカ（２、４）によって生成され、該一群のラウドスピーカは、第１のラウドスピーカ（２）および少なくとも第２のラウドスピーカ（４）を備え、該スピーカの各々は、オーディオチャネルを介してオーディオ信号が供給され、該方法は、
該第２のラウドスピーカ（４）に供給する少なくとも該オーディオチャネルのための位相イコライゼーションフィルタのフィルタ係数を計算することであって、該位相イコライゼーションフィルタの位相応答は、該少なくとも１つの聴取位置（１０）における両耳の位相差（Δφｍｎ）、または２つ以上の聴取位置（１０、１１）の上に平均される平均的な両耳の位相差（ｍΔφｍｎ）が所定の周波数範囲内で最小化されるように設計される、ことと、
該位相イコライゼーションフィルタを該それぞれのオーディオチャネルに適用することと
を包含する、方法。
【請求項２】
前記位相イコライゼーションフィルタの係数を計算するステップは、
少なくとも１つのオーディオチャネルに適用可能な周波数および位相シフトに依存する位相差のアレイ内の最小限の探索を実行することであって、該最小限の探索は周波数（ｆｍ）の関数として最適な位相シフト（φＸ）を表す最適位相関数φＸ，ＦＩＬＴ（ｆｍ）を生み出す、ことを包含する、請求項１に記載の方法。
【請求項３】
前記位相イコライゼーションフィルタの係数を計算するステップは、
各聴取位置（１０、１１）に対して、それぞれの聴取位置（１０、１１）に割り当てられた群の各ラウドスピーカ（２、４）に対する両耳の伝達特性を決定することと、
所定の周波数範囲から一組の周波数（ｆｍ）を選択し、所定の位相範囲から一組の位相シフト（φｎ）を選択することと、
各聴取位置（１０、１１）、該一組の周波数の各周波数（ｆｍ）、および該一組の位相シフトの各位相シフト（φｎ）に対して両耳の位相差（Δφｍｎ）を計算し、該計算に関して、オーディオ信号が各ラウドスピーカ（２、４）に供給されると仮定することであって、該少なくとも１個の第２のラウドスピーカ（４）に供給される該オーディオ信号は、該第１のラウドスピーカ（２）に供給される該オーディオ信号に対してそれぞれの位相シフト（φｎ）によって位相シフトされ、従って、該それぞれの聴取位置（１０、１１）に対する両耳の位相差（Δφｍｎ）のアレイを提供する、ことと、
該少なくとも１つの聴取位置（１０、１１）において、該両耳の位相差（Δφｍｎ）の重み付け平均を計算することによって、平均的な両耳の位相差（ｍΔφｍｎ）のアレイを提供することと、
該平均的な両耳の位相差（ｍΔφｍｎ）のアレイにおいて、各周波数（ｆｍ）に対する最適な位相シフト（φｎ）を探索することであって、該最適な位相シフト（φＸ）は該平均的な両耳の位相差（ｍΔφｍｎ）の最小値を生み出し、従って、結果として、周波数（ｆｍ）の関数として該最適な位相シフト（φＸ）を表す最適位相関数φＸ，ＦＩＬＴ（ｆｍ）をもたらす、ことと
を包含する、請求項１に記載の方法。
【請求項４】
各考えられる聴取位置（１０、１１）において前記両耳の位相差（Δφｍｎ）を計算するステップは、
前記一組の周波数の各周波数（ｆｍ）および前記一組の位相シフト（φｎ）に対して、各聴取位置（１０、１１）において重なりスペクトル値を計算することと、
各計算された重なりスペクトル値に対する該重なりスペクトルの位相を計算することであって、該重なりスペクトルの位相は、該それぞれの聴取位置（１０、１１）において該両耳の位相差（Δφｍｎ）を表す、ことと
を包含する、請求項３に記載の方法。
【請求項５】
前記最適位相関数φＸ，ＦＩＬＴ（ｆｍ）を近似する位相応答を提供するように設計されるデジタル位相イコライゼーションフィルタを提供するステップをさらに包含する、請求項２〜４のうちの１項に記載の方法。
【請求項６】
前記両耳の伝達特性を決定するステップは、
広帯域試験信号を各ラウドスピーカ(２、４、３)に順次に供給することと、
各聴取位置（１０、１１）に到達する、結果として生じる音響信号を両耳で測定することと、
ラウドスピーカ（２、４、３）と聴取位置（１０、１１）との各対に対して対応する両耳の伝達特性を計算することと
を包含する、請求項３〜５のうちの１項に記載の方法。
【請求項７】
前記位相イコライジングフィルタの位相応答を計算する前に、前記最適位相関数φＸ，ＦＩＬＴ（ｆｍ）のスムージングを行うステップをさらに包含する、請求項２〜６のうちの１項に記載の方法。
【請求項８】
前記スムージングを行うステップは、非線形の複素スムージングフィルタを用いて行われる、請求項７に記載の方法。
【請求項９】
前記スムージングステップは、周波数の増加にともない動的応答が減少するスムージングフィルタを用いて行われる、請求項７または８に記載の方法。
【請求項１０】
前記位相イコライゼーションフィルタのフィルタ係数を計算するステップは、
所定の周波数範囲から一組の周波数（ｆｍ）を選択し、所定の位相範囲から一組の位相シフト（φｎ）を選択することと、
音場を生成するために、各選択された周波数（ｆｍ）に対して、該それぞれの周波数（ｆｍ）を有するオーディオ信号を各ラウドスピーカ（２、４）に供給することであって、該少なくとも１個の第２のラウドスピーカ（４）に供給される該オーディオ信号は、該第１のラウドスピーカ（２）に供給される該オーディオ信号に対して該それぞれの位相シフト（φｎ）によって位相シフトされる、ことと、
位相シフト（φｎ）と周波数（ｆｍ）との各組み合わせに対して、各聴取位置（１０、１１）に到達する、結果として生じる音響信号を両耳で測定することと、
該それぞれの両耳で測定された音響信号から各聴取位置（１０、１１）に対する両耳の位相差（Δφｍｎ）を計算し、従って、位相シフト（φｎ）と周波数（ｆｍ）との各組み合わせに対する両耳の位相差値を含む該各聴取位置（１０、１１）に対する両耳の位相差（Δφｍｎ）のアレイを提供することと、
該少なくとも１つの聴取位置（１０、１１）において該両耳の位相差（Δφｍｎ）の重み付けされた平均を計算することによって、平均的な両耳の位相差（ｍΔφｍｎ）のアレイを提供することと、
該平均的な両耳の位相差（ｍΔφｍｎ）のアレイにおいて、各周波数（ｆｍ）に対する最適な位相シフト（φｎ）を探索することであって、該最適な位相シフト（φＸ）は該平均的な両耳の位相差（ｍΔφｍｎ）の最小値を生み出し、従って、結果として周波数（ｆｍ）の関数として該最適な位相シフト（φＸ）を表す最適位相関数φＸ，ＦＩＬＴ（ｆｍ）をもたらす、ことと、
該最適な位相φＸ，ＦＩＬＴ（ｆｍ）関数に近似する位相イコライゼーションフィルタに対する位相応答を計算することと
を包含する、請求項１に記載の方法。
【請求項１１】
聴取室内の少なくとも１つの聴取位置（１０）において音響ローカライゼーションを最適化するシステムであって、該システムは、
音場を生成するために該少なくとも１つの聴取位置（１０、１１）に割り当てられたラウドスピーカ（２、４）の群であって、該ラウドスピーカの群は、第１のラウドスピーカ（２）および少なくとも第２のラウドスピーカ（４）を含む、ラウドスピーカ（２、４）の群と、
オーディオ信号を各ラウドスピーカにそれぞれのオーディオチャネルを介して提供する信号源と、
該第２のラウドスピーカ（４）に供給する少なくとも該オーディオチャネル適用される位相イコライゼーションフィルタのフィルタ係数を計算するように構成される信号処理ユニットであって、該位相イコライゼーションフィルタの位相応答は、該少なくとも１つの聴取位置（１０）における両耳の位相差（Δφｍｎ）または複数の聴取位置（１０、１１）の上に平均される平均的な両耳の位相差（ｍΔφｍｎ）が所定の周波数範囲内で最小化されるように設計される、信号処理ユニットと
を備えている、システム。
【請求項１２】
前記位相イコライゼーションフィルタの係数を計算するために、前記信号処理ユニットは、少なくとも１つのオーディオチャネルに適用可能な周波数および位相シフトに依存する位相差のアレイ内の最小限の探索を実行するように構成され、該最小限の探索は、周波数（ｆｍ）の関数として最適な位相シフト（φＸ）を表す最適位相関数φＸ，ＦＩＬＴ（ｆｍ）を生み出す、請求項１１に記載のシステム。
【請求項１３】
位相イコライゼーションフィルタの係数を計算するために、信号処理ユニットは、
各聴取位置（１０、１１）に対して、該それぞれの聴取位置（１０、１１）に割り当てられた群の各ラウドスピーカ（２、４）に対する両耳の伝達特性を決定することと、
所定の周波数範囲から一組の周波数（ｆｍ）を選択し、所定の位相範囲から一組の位相シフト（φｎ）を選択することと、
各聴取位置（１０、１１）、該一組の周波数の各周波数（ｆｍ）、および該一組の位相シフトの各位相シフト（φｎ）に対して両耳の位相差（Δφｍｎ）を計算し、該計算に関して、オーディオ信号が各ラウドスピーカ（２、４）に供給されると仮定することであって、該少なくとも１個の第２のラウドスピーカ（４）に供給される該オーディオ信号は、該第１のラウドスピーカ（２）に供給される該オーディオ信号に対してそれぞれの位相シフト（φｎ）によって位相シフトされ、従って、該それぞれの聴取位置（１０、１１）に対して両耳の位相差（Δφｍｎ）のアレイを提供する、ことと、
該少なくとも１つの聴取位置（１０、１１）において、該両耳の位相差（Δφｍｎ）の重み付け平均を計算することによって、平均的な両耳の位相差（ｍΔφｍｎ）のアレイを提供することと、
該平均的な両耳の位相差（ｍΔφｍｎ）のアレイにおいて、各周波数（ｆｍ）に対する最適な位相シフト（φｎ）を探索することであって、該最適な位相シフト（φＸ）は該平均的な両耳の位相差（ｍΔφｍｎ）の最小値を生み出し、従って、結果として、周波数（ｆｍ）の関数として該最適な位相シフト（φＸ）を表す最適位相関数φＸ，ＦＩＬＴ（ｆｍ）をもたらす、ことと、
前記最適な位相φＸ，ＦＩＬＴ（ｆｍ）関数に近似する位相イコライゼーションフィルタに対する位相応答を計算することと
を行うように構成される、請求項１１または１２に記載のシステム。
【請求項１４】
前記位相イコライジングフィルタの位相応答を計算する前に、前記最適位相関数φＸ，ＦＩＬＴ（ｆｍ）のスムージングを行うように構成されるスムージングフィルタをさらに備えている、請求項１２または１３に記載のシステム。
【請求項１５】
前記スムージングフィルタは、非線形の複素スムージングフィルタであり、動的応答が、周波数が増加するにつれ動的応答が減少する、請求項１４に記載のシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
（技術分野）
本発明は、オーディオシステムにおける位相イコライゼーションのための方法に関し、特に、自動車の客室内の任意の聴取位置におけるステレオ信号の両耳間の時間差を最小化する方法に関する。
【背景技術】
【０００２】
（背景）
特に高級リムジンにおける先進の車両サウンドシステムは、通常、車両の客室内の異なる位置に置かれた複数の単一のスピーカとそれらのアレイとを備えている高度に複雑な構成を有し、これらの単一のスピーカおよびそれらのアレイは、典型的には、多様な周波数帯域専用である（例えば、サブウーハ、ウーハ、中間域スピーカおよびトウィータスピーカなど）。
【０００３】
そのような先行技術のサウンドシステムは、各々の場合に関係する車両のタイプのために音響技術者によって手動で調整、すなわち最適化され、主として彼らの経験に基づき、そして彼らの訓練された聴覚に基づいて、主観的に所望の音響品質を達成する。この目的のために、彼らは、典型的には、双二次フィルタ（例えば、ハイパスフィルタ、バンドパスフィルタ、ローパスフィルタ、オールパスフィルタ）、双一次フィルタ、デジタルディレイライン、クロスオーバフィルタなどの公知の信号処理アセンブリと、信号動的応答を変化させるアセンブリ（コンプレッサ、リミッタ、エクスパンダ、ノイズゲートなど）とを利用し、クロスオーバフィルタ、ディレイラインおよび周波数応答の大きさの関連するカットオフ周波数パラメータを規定し、その結果、最終的に、動力車両におけるサウンドシステムの音の印象は、サウンドシステムのスペクトルバランス（すなわち、調性、調性の優秀性）および周囲（すなわち、空間バランス、音の空間性）に関して最適化され、達成される。
【０００４】
そのような調整の目的は、すべての聴取位置、換言すると車両の客室におけるすべての着席位置（聴取位置）において、最適化されるような音を達成することであり、従って、動力車両のサウンドシステムの調整において複雑さを実質的にさらに追加する。オーディオ信号が周囲に知覚される仕方およびオーディオ信号がステレオ音響でローカライズされる仕方に大いに影響を及ぼすのは、特に、動力車両における様々な聴取位置または着席位置における両耳間の時間差である。
【発明の概要】
【発明が解決しようとする課題】
【０００５】
車両の客室内の任意の聴取位置、特に自動車内の対称軸外に配列された聴取位置における両耳間の時間差を最小化することを可能にする方法に対する一般的なニーズがある。
【課題を解決するための手段】
【０００６】
（概要）
聴取室内の少なくとも１つの聴取位置において音響ローカライゼーションを最適化する方法が開示される。音場は、少なくとも１つの聴取位置に割り当てられたラウドスピーカの群によって生成され、ラウドスピーカの群は、第１のラウドスピーカおよび少なくとも第２のラウドスピーカを備え、スピーカの各々は、オーディオチャネルを介してオーディオ信号が供給される。方法は、第２のラウドスピーカに供給する少なくともオーディオチャネルのための位相イコライゼーションフィルタのフィルタ係数を計算することを包含し、位相イコライゼーションフィルタの位相応答は、少なくとも１つの聴取位置における両耳の位相差または複数の聴取位置の上に平均される平均的な両耳の位相差が所定の周波数範囲内で最小化され、位相イコライゼーションフィルタをそれぞれのオーディオチャネルに適用するように設計される。
【０００７】
例えば、本発明は以下の項目を提供する。
（項目１）
聴取室内の少なくとも１つの聴取位置（１０）において音響ローカライゼーションを最適化する方法であって、音場は、該少なくとも１つの聴取位置（１０、１１）に割り当てられた一群のラウドスピーカ（２、４）によって生成され、該一群のラウドスピーカは、第１のラウドスピーカ（２）および少なくとも第２のラウドスピーカ（４）を備え、該スピーカの各々は、オーディオチャネルを介してオーディオ信号が供給され、該方法は、
該第２のラウドスピーカ（４）に供給する少なくとも該オーディオチャネルのための位相イコライゼーションフィルタのフィルタ係数を計算することであって、該位相イコライゼーションフィルタの位相応答は、該少なくとも１つの聴取位置（１０）における両耳の位相差（Δφｍｎ）、または２つ以上の聴取位置（１０、１１）の上に平均される平均的な両耳の位相差（ｍΔφｍｎ）が所定の周波数範囲内で最小化されるように設計される、ことと、
該位相イコライゼーションフィルタを該それぞれのオーディオチャネルに適用することと
を包含する、方法。
（項目２）
上記位相イコライゼーションフィルタの係数を計算するステップは、
少なくとも１つのオーディオチャネルに適用可能な周波数および位相シフトに依存する位相差のアレイ内の最小限の探索を実行することであって、該最小限の探索は周波数（ｆｍ）の関数として最適な位相シフト（φＸ）を表す最適位相関数φＸ，ＦＩＬＴ（ｆｍ）を生み出す、ことを包含する、上記項目のいずれかに記載の方法。
（項目３）
上記位相イコライゼーションフィルタの係数を計算するステップは、
各聴取位置（１０、１１）に対して、それぞれの聴取位置（１０、１１）に割り当てられた群の各ラウドスピーカ（２、４）に対する両耳の伝達特性を決定することと、
所定の周波数範囲から一組の周波数（ｆｍ）を選択し、所定の位相範囲から一組の位相シフト（φｎ）を選択することと、
各聴取位置（１０、１１）、該一組の周波数の各周波数（ｆｍ）、および該一組の位相シフトの各位相シフト（φｎ）に対して両耳の位相差（Δφｍｎ）を計算し、該計算に関して、オーディオ信号が各ラウドスピーカ（２、４）に供給されると仮定することであって、該少なくとも１個の第２のラウドスピーカ（４）に供給される該オーディオ信号は、該第１のラウドスピーカ（２）に供給される該オーディオ信号に対してそれぞれの位相シフト（φｎ）によって位相シフトされ、従って、該それぞれの聴取位置（１０、１１）に対する両耳の位相差（Δφｍｎ）のアレイを提供する、ことと、
該少なくとも１つの聴取位置（１０、１１）において、該両耳の位相差（Δφｍｎ）の重み付け平均を計算することによって、平均的な両耳の位相差（ｍΔφｍｎ）のアレイを提供することと、
該平均的な両耳の位相差（ｍΔφｍｎ）のアレイにおいて、各周波数（ｆｍ）に対する最適な位相シフト（φｎ）を探索することであって、該最適な位相シフト（φＸ）は該平均的な両耳の位相差（ｍΔφｍｎ）の最小値を生み出し、従って、結果として、周波数（ｆｍ）の関数として該最適な位相シフト（φＸ）を表す最適位相関数φＸ，ＦＩＬＴ（ｆｍ）をもたらす、ことと
を包含する、上記項目のいずれかに記載の方法。
（項目４）
各考えられる聴取位置（１０、１１）において上記両耳の位相差（Δφｍｎ）を計算するステップは、
上記一組の周波数の各周波数（ｆｍ）および上記一組の位相シフト（φｎ）に対して、各聴取位置（１０、１１）において重なりスペクトル値を計算することと、
各計算された重なりスペクトル値に対する該重なりスペクトルの位相を計算することであって、該重なりスペクトルの位相は、該それぞれの聴取位置（１０、１１）において該両耳の位相差（Δφｍｎ）を表す、ことと
を包含する、上記項目のいずれかに記載の方法。
（項目５）
上記最適位相関数φＸ，ＦＩＬＴ（ｆｍ）を近似する位相応答を提供するように設計されるデジタル位相イコライゼーションフィルタを提供するステップをさらに包含する、上記項目のいずれかに記載の方法。
（項目６）
上記両耳の伝達特性を決定するステップは、
広帯域試験信号を各ラウドスピーカ(２、４、３)に順次に供給することと、
各聴取位置（１０、１１）に到達する、結果として生じる音響信号を両耳で測定することと、
ラウドスピーカ（２、４、３）と聴取位置（１０、１１）との各対に対して対応する両耳の伝達特性を計算することと
を包含する、上記項目のいずれかに記載の方法。
（項目７）
上記位相イコライジングフィルタの位相応答を計算する前に、上記最適位相関数φＸ，ＦＩＬＴ（ｆｍ）のスムージングを行うステップをさらに包含する、上記項目のいずれかに記載の方法。
（項目８）
上記スムージングを行うステップは、非線形の複素スムージングフィルタを用いて行われる、上記項目のいずれかに記載の方法。
（項目９）
上記スムージングステップは、周波数の増加にともない動的応答が減少するスムージングフィルタを用いて行われる、上記項目のいずれかに記載の方法。
（項目１０）
上記位相イコライゼーションフィルタのフィルタ係数を計算するステップは、
所定の周波数範囲から一組の周波数（ｆｍ）を選択し、所定の位相範囲から一組の位相シフト（φｎ）を選択することと、
音場を生成するために、各選択された周波数（ｆｍ）に対して、該それぞれの周波数（ｆｍ）を有するオーディオ信号を各ラウドスピーカ（２、４）に供給することであって、該少なくとも１個の第２のラウドスピーカ（４）に供給される該オーディオ信号は、該第１のラウドスピーカ（２）に供給される該オーディオ信号に対して該それぞれの位相シフト（φｎ）によって位相シフトされる、ことと、
位相シフト（φｎ）と周波数（ｆｍ）との各組み合わせに対して、各聴取位置（１０、１１）に到達する、結果として生じる音響信号を両耳で測定することと、
該それぞれの両耳で測定された音響信号から各聴取位置（１０、１１）に対する両耳の位相差（Δφｍｎ）を計算し、従って、位相シフト（φｎ）と周波数（ｆｍ）との各組み合わせに対する両耳の位相差値を含む該各聴取位置（１０、１１）に対する両耳の位相差（Δφｍｎ）のアレイを提供することと、
該少なくとも１つの聴取位置（１０、１１）において該両耳の位相差（Δφｍｎ）の重み付けされた平均を計算することによって、平均的な両耳の位相差（ｍΔφｍｎ）のアレイを提供することと、
該平均的な両耳の位相差（ｍΔφｍｎ）のアレイにおいて、各周波数（ｆｍ）に対する最適な位相シフト（φｎ）を探索することであって、該最適な位相シフト（φＸ）は該平均的な両耳の位相差（ｍΔφｍｎ）の最小値を生み出し、従って、結果として周波数（ｆｍ）の関数として該最適な位相シフト（φＸ）を表す最適位相関数φＸ，ＦＩＬＴ（ｆｍ）をもたらす、ことと、
該最適な位相φＸ，ＦＩＬＴ（ｆｍ）関数に近似する位相イコライゼーションフィルタに対する位相応答を計算することと
を包含する、上記項目のいずれかに記載の方法。
（項目１１）
聴取室内の少なくとも１つの聴取位置（１０）において音響ローカライゼーションを最適化するシステムであって、該システムは、
音場を生成するために該少なくとも１つの聴取位置（１０、１１）に割り当てられたラウドスピーカ（２、４）の群であって、該ラウドスピーカの群は、第１のラウドスピーカ（２）および少なくとも第２のラウドスピーカ（４）を含む、ラウドスピーカ（２、４）の群と、
オーディオ信号を各ラウドスピーカにそれぞれのオーディオチャネルを介して提供する信号源と、
該第２のラウドスピーカ（４）に供給する少なくとも該オーディオチャネル適用される位相イコライゼーションフィルタのフィルタ係数を計算するように構成される信号処理ユニットであって、該位相イコライゼーションフィルタの位相応答は、該少なくとも１つの聴取位置（１０）における両耳の位相差（Δφｍｎ）または複数の聴取位置（１０、１１）の上に平均される平均的な両耳の位相差（ｍΔφｍｎ）が所定の周波数範囲内で最小化されるように設計される、信号処理ユニットと
を備えている、システム。
（項目１２）
上記位相イコライゼーションフィルタの係数を計算するために、上記信号処理ユニットは、少なくとも１つのオーディオチャネルに適用可能な周波数および位相シフトに依存する位相差のアレイ内の最小限の探索を実行するように構成され、該最小限の探索は、周波数（ｆｍ）の関数として最適な位相シフト（φＸ）を表す最適位相関数φＸ，ＦＩＬＴ（ｆｍ）を生み出す、上記項目のいずれかに記載のシステム。
（項目１３）
位相イコライゼーションフィルタの係数を計算するために、信号処理ユニットは、
各聴取位置（１０、１１）に対して、該それぞれの聴取位置（１０、１１）に割り当てられた群の各ラウドスピーカ（２、４）に対する両耳の伝達特性を決定することと、
所定の周波数範囲から一組の周波数（ｆｍ）を選択し、所定の位相範囲から一組の位相シフト（φｎ）を選択することと、
各聴取位置（１０、１１）、該一組の周波数の各周波数（ｆｍ）、および該一組の位相シフトの各位相シフト（φｎ）に対して両耳の位相差（Δφｍｎ）を計算し、該計算に関して、オーディオ信号が各ラウドスピーカ（２、４）に供給されると仮定することであって、該少なくとも１個の第２のラウドスピーカ（４）に供給される該オーディオ信号は、該第１のラウドスピーカ（２）に供給される該オーディオ信号に対してそれぞれの位相シフト（φｎ）によって位相シフトされ、従って、該それぞれの聴取位置（１０、１１）に対して両耳の位相差（Δφｍｎ）のアレイを提供する、ことと、
該少なくとも１つの聴取位置（１０、１１）において、該両耳の位相差（Δφｍｎ）の重み付け平均を計算することによって、平均的な両耳の位相差（ｍΔφｍｎ）のアレイを提供することと、
該平均的な両耳の位相差（ｍΔφｍｎ）のアレイにおいて、各周波数（ｆｍ）に対する最適な位相シフト（φｎ）を探索することであって、該最適な位相シフト（φＸ）は該平均的な両耳の位相差（ｍΔφｍｎ）の最小値を生み出し、従って、結果として、周波数（ｆｍ）の関数として該最適な位相シフト（φＸ）を表す最適位相関数φＸ，ＦＩＬＴ（ｆｍ）をもたらす、ことと、
上記最適な位相φＸ，ＦＩＬＴ（ｆｍ）関数に近似する位相イコライゼーションフィルタに対する位相応答を計算することと
を行うように構成される、上記項目のいずれかに記載のシステム。
（項目１４）
上記位相イコライジングフィルタの位相応答を計算する前に、上記最適位相関数φＸ，ＦＩＬＴ（ｆｍ）のスムージングを行うように構成されるスムージングフィルタをさらに備えている、上記項目のいずれかに記載のシステム。
（項目１５）
上記スムージングフィルタは、非線形の複素スムージングフィルタであり、動的応答が、周波数が増加するにつれ動的応答が減少する、上記項目のいずれかに記載のシステム。
【０００８】
(摘要）
聴取室内の少なくとも１つの聴取位置において音響ローカライゼーションを最適化する方法が開示される。音場は、少なくとも１つの聴取位置に割り当てられたラウドスピーカの群によって生成され、ラウドスピーカの群は、第１のラウドスピーカおよび少なくとも第２のラウドスピーカを備え、スピーカの各々は、オーディオチャネルを介してオーディオ信号が供給される。方法は、第２のラウドスピーカに供給する少なくともオーディオチャネルのための位相イコライゼーションフィルタのフィルタ係数を計算することを包含し、位相イコライゼーションフィルタの位相応答は、少なくとも１つの聴取位置における両耳の位相差または複数の聴取位置の上に平均される平均的な両耳の位相差が所定の周波数範囲内で最小化され、位相イコライゼーションフィルタをそれぞれのオーディオチャネルに適用するように設計される。
【０００９】
本発明は、以下の図面および説明を参照してより良く理解され得る。図における構成要素は、必ずしも一定の縮尺ではなく、代わりに、本発明の原理を例示することに強調がなされている。さらに図において同様の参照数字は、対応する部品を示す。
【図面の簡単な説明】
【００１０】
【図１】図１は、ダミーヘッドを用いて車両の中心線（対称軸）において測定される場合の両耳の位相差を例示する図である。
【図２】図２は、車両の中心線外に位置を定められる運転席の聴取位置においてダミーヘッドを用いて測定される場合の両耳の位相差を例示する図である。
【図３】図３は、例としてオーディオシステムのラウドスピーカの配列を示す、車両の客室において試験されたすべての測定位置の上面図である。
【図４】図４は、車両の客室において試験されたすべての測定位置の側部図である。
【図５】図５は、前部左チャネルに対して０°から１８０°まで、１°のステップで連続位相シフトが適用される、車両における２個の異なる着席位置における周波数の関数としての両耳の伝達関数の重なりスペクトルの位相の三次元図である。
【図６】図６は、図５に示されるような重なりスペクトルの位相の三次元図の上面図であり、両耳の重なりスペクトルの位相を最小化する、前部左チャネルに対する１周波数当たりの位相シフトを示す。
【図７】図７は、前部左チャネルに対する最適な位相シフトを例示する図であり、この最適な位相シフトは、重なりスペクトルの位相の最大の最小化を、そして従って、平均して車両の両方の前部着席位置における最適な水平ローカライゼーションという結果をもたらす。
【図８】図８は、図７に示されるように、最適な位相シフトに近似するための、前部左チャネルに対する位相イコライザの群遅延を例示する図である。
【図９】図９は、図８に示されるような前部左チャネルの位相イコライザのインパルス応答を例示する図（下図：大きさの線形表現、上図：大きさの対数表現）である。
【図１０】図１０は、図８に示されるような前部左チャネルの位相イコライザのボードダイアグラム（下図：周波数応答の大きさ、上図：位相周波数応答）である。
【図１１ａ】図１１は、位相イコライザの適用前および適用後の車両内のすべての４つの着席位置における両耳の重なりスペクトルの位相差を例示する図である。
【図１１ｂ】図１１は、位相イコライザの適用前および適用後の車両内のすべての４つの着席位置における両耳の重なりスペクトルの位相差を例示する図である。
【図１１ｃ】図１１は、位相イコライザの適用前および適用後の車両内のすべての４つの着席位置における両耳の重なりスペクトルの位相差を例示する図である。
【図１１ｄ】図１１は、位相イコライザの適用前および適用後の車両内のすべての４つの着席位置における両耳の重なりスペクトルの位相差を例示する図である。
【発明を実施するための形態】
【００１１】
（詳細な説明）
オーディオシステムを手動で調整するための音響的手段を用いることは、時代遅れであるが、その目的は、特に、位相を微調整するために、例えば、個々の増幅器チャネルの遅延を主としてイコライズするディレイラインを利用することである。位相応答を直接的に修正するために、オールパスフィルタが通常利用される。しかしながら、個々のラウドスピーカの伝達帯域を制限するために主として利用されるクロスオーバフィルタがまた、再生されるオーディオ信号の位相応答を微調整する。部分的に、勾配が異なる多様なタイプのフィルタ（Ｂｕｔｔｅｒｗｏｒｔｈ、Ｂｅｓｓｅｌ、Ｌｉｎｋｗｉｔｚ−Ｒｉｌｅｙなど）が、異なる位相遷移によって明確にサウンドを微調整するために意図的に利用される。
【００１２】
強力なデジタル信号プロセッサの利用可能性は、非常に大きなフィルタの柔軟性を（低コストで）もたらし、例えば、大きさおよび位相周波数応答の各々が、互いに別個に設定されることを可能にする。しかしながら、優先度は、ＦＩＲフィルタの利用に向けられる。なぜならば、適切な等価ＩＩＲフィルタを実現することは、ＩＩＲフィルタのより低いフィルタ次数に起因して、より安価に実装されるにもかかわらず、今もなお極めて困難であるからである。
【００１３】
ＦＩＲフィルタは、有限のインパルス応答を有することと、通常アナログ信号のサンプリング周波数によって決定される別個の時間的なステップで作動することとによって特徴付けられる。Ｎ次のＦＩＲフィルタは、次の差分方程式
【００１４】
【数１】
2011097561.tif 000003
によって記述され、ここで、ｙ（ｎ）は、時間ｎ（ｎは、サンプル番号であり、従って時間のインデックスである）の時点における開始値であり、フィルタ係数ｂｉによって重み付けされる、現時点およびＮ個の最近サンプリングされた入力値ｘ（ｎ−Ｎ−１）〜ｘ（ｎ）の合計から得られ、それによって、所望の伝達関数が、フィルタ係数ｂｉを指定することによって実現される。
【００１５】
例えば、分割（ｐａｒｔｉｔｉｏｎｅｄ）高速畳み込みなどの多様な信号処理アルゴリズムを利用することまたはフィルタバンクを用いることが、実際に任意の市販されているデジタル信号プロセッサによって達成可能であるような、十分に長いＦＩＲフィルタを実現することを可能にする。このことは、実装に含まれる問題を目立たなくさせ、音響の永続的な向上のためのオーディオ信号の位相周波数応答の正しく導かれた微調整と、特に、車両の客室内の多様な聴取位置におけるオーディオ信号のローカライゼーションとを可能にする。
【００１６】
ローカライゼーションは、両方の耳（両耳の聴覚）による聴覚の結果として、音源の水平方向および音源からの距離を認識することであると理解される。サウンドがどの側から到来するかを決定するために、人の耳の知覚の感覚は、左、真正面、右などの方向を区別することに関して、両方の耳の間の遅延およびレベルにおける差を評価する。
【００１７】
人の耳は、知覚されるサウンドがどの方向から到来するかを決定することに関して、両方の耳の間の遅延における差（「両耳間の時間差」と呼ばれ、ＩＴＤと略記する）を主として評価する。右から到来するサウンドは、左耳よりも早く右耳に到達し、区別が、低い周波数における位相遅延の評価と、高い周波数における群遅延の評価と、両方の耳の間の周波数の関数としてのレベル差（「両耳間のレベル差」と呼ばれ、ＩＬＤと略記する）の評価との間でなされる。
【００１８】
右から到来するサウンドは、右耳において左耳よりも高いレベルを有する。なぜならば、頭部が左耳における信号をさえぎるからである。これらのレベル差は、周波数の関数であり、周波数を増大させることによって増大する。約８００Ｈｚ未満の低い周波数において評価されるのは、特に、遅延における差（位相遅延または遅延における位相差）であるのに対して、約１５００Ｈｚ超の高い周波数においては、特に、レベルの差が評価される。両方の機構が役割を果たす（「トレード（ｔｒａｄｉｎｇ）」）重複領域が間に配置される。
【００１９】
低い周波数においては、耳から耳までの距離ｄ＝２１．５ｃｍを有する頭部の寸法が、０．６３ｍｓの遅延における差に対応し、サウンドの波長の半分よりも小さい。両方の耳の間の遅延における差を耳が非常に正確に評価し得るのはこの点であるが、レベル差は、かなり小さいので、どのような精度によっても評価され得ない。８０Ｈｚ未満の周波数は、もはやそれらの方向においてローカライズされ得ない。低い周波数において、頭部の寸法は、サウンドの波長よりも小さい。ここにおいて、耳は、もはや遅延における差から方向を明示的に決定し得ないが、両耳間のレベル差は、より大きくなるので、次いで、耳によって評価される。
【００２０】
そのような変数を測定することにおいて目的の結果を達成するために、人の頭部の形状および反射／回折の特性を複製した、いわゆるダミーヘッドが使用される。耳の代わりに、そのようなダミーヘッドは、様々な条件下で到着する信号を測定するために、対応して位置決めされた２個のマイクロホンを有する。例えば、そのようなダミーヘッドの位置は、聴取室において変わり得る。
【００２１】
両耳間のレベル差に加えて（より高い周波数においても同様に）、両方の耳の間の群遅延が評価され、群遅延は、新しいサウンドが発生する場合に、その方向が、両方の耳の間のサウンド発生における遅延から決定され得ることを意味している。この機構は、特に残響環境において重要である。サウンドの発生において、直接的サウンドがすでに聴取者に到達しているが、反射サウンドがまだそのうえに到達する短い時間が存在する。耳は、この時間の間隙を開始時間において用いることにより、方向を決定し、反射が原因で、もはや方向を明示的に決定することが可能でなくなるまで、測定された方向を保持する。この現象は、「ハース効果」、「先行音効果」、または「第一波面の法則」と呼ばれる。
【００２２】
音源のローカライゼーションは、いわゆる周波数群で行われる。人の聴覚範囲は、約２４の周波数群（各々１Ｂａｒｋまたは１００Ｍｅｌの幅）に分割される。方向を決定するために、人の耳は、１つの周波数群に入る共通の信号成分を評価する。
【００２３】
評価することにおいて、人の耳は、臨界周波数群または臨界帯域幅（ＣＢ）とも呼ばれる限定された周波数帯域に存在するサウンドキューを組み合わせるが、臨界周波数群または臨界帯域幅の幅は、特定の周波数帯域に存在するサウンドを、これらのサウンドから発する音響心理学的な聴覚感覚に関して、共通の聴覚感覚の中に組み合わせることができる人の耳に基づいている。単一周波数群の中に存在するサウンド事象は、様々な周波数群に存在するサウンドのそれと異なる効果を有する。例えば、ある周波数群の中で同じレベルを有する２個の音調が、様々な周波数群に存在するときよりも柔らかく聞かれる。
【００２４】
エネルギーが同じであり、かつマスカが中心周波数としてテストトーンの周波数を有する周波数帯域に入る場合には、マスカ内のテストトーンが可聴であるので、周波数群の所望の帯域幅が決定され得る。低い周波数において、周波数群は、１００Ｈｚの帯域幅を有する。５００Ｈｚ超の周波数において、周波数群は、関係する周波数群の中心周波数の約２０％に達する帯域幅を有する（Ｚｗｉｃｋｅｒ，Ｅ．；Ｆａｓｔｌ，Ｈ．Ｐｓｙｃｈｏａｃｏｕｓｔｉｃｓ−Ｆａｃｔｓ  ａｎｄ  Ｍｏｄｅｌｓ，２ｎｄ  ｅｄｉｔｉｏｎ，Ｓｐｒｉｎｇｅｒ−Ｖｅｒｌａｇ，Ｂｅｒｌｉｎ／Ｈｅｉｄｅｌｂｅｒｇ／Ｎｅｗ  Ｙｏｒｋ，１９９９）。
【００２５】
全部の聴覚範囲にわたってすべての臨界周波数群を一列に並べることは、ピッチ（単位が「Ｂａｒｋ」である）と呼ばれる、聴覚指向の非線形な周波数尺度を結果としてもたらす。それは周波数軸の歪んだ尺度化を表し、それによって、周波数群が、各々の位置（ｐｏｉｎｔ）において正確に１Ｂａｒｋの同一の幅を有することを特徴とする。周波数とピッチとの非線形関係は、その起源を基底板上の周波数／場所の変換に有している。ピッチ関数は、Ｚｗｉｃｋｅｒによって、聴取閾値および音の大きさ（ｌｏｕｄｎｅｓｓ）に関する試験に基づき、表および式の形態で定式化されている（Ｚｗｉｃｋｅｒ，Ｅ．；Ｆａｓｔｌ，Ｈ．Ｐｓｙｃｈｏａｃｏｕｓｔｉｃｓ−Ｆａｃｔｓ  ａｎｄ  Ｍｏｄｅｌｓ，２ｎｄ  ｅｄｉｔｉｏｎ，Ｓｐｒｉｎｇｅｒ−Ｖｅｒｌａｇ，Ｂｅｒｌｉｎ／Ｈｅｉｄｅｌｂｅｒｇ／Ｎｅｗ  Ｙｏｒｋ，１９９９）。ちょうど２４個の周波数群が、０ｋＨｚ〜１６ｋＨｚの可聴周波数範囲に一列に並べられ得、それによって、対応するピッチ範囲が、０Ｂａｒｋから２４Ｂａｒｋに達することが示された。１Ｂａｒｋ当たりのピッチｚは、式
【００２６】
【数２】
2011097561.tif 000004
に示され、対応する周波数群の幅ΔｆＧは、式
【００２７】
【数３】
2011097561.tif 000005
に示される。
【００２８】
閉じた室において、耳に作用するのは、サウンドシステムの方向からのサウンドだけでなく、壁から反射されたサウンドもまた作用する。しかしながら、方向を決定することにおいて、耳は、後から到着する任意の反射音ではなく、まさに最初に到着する直接的サウンドを評価し（第一波面の法則）、それによって、正しく音源の方向を決定することが可能であり続ける。この目的のために、耳は、時間と共に様々な周波数群における音の大きさの強い変化を評価する。１つ以上の周波数群において音の大きさの強い増大がある場合には、このことは、音源から新しく発生する直接的サウンドかまたは特性が変わる信号を含むあらゆる可能性を有する。この短い時間が、耳によって方向を決定するために用いられる。
【００２９】
後から到着する反射サウンドは、もはや関係する周波数群における音の大きさに追加されず、それによって、後から到着する反射サウンドが、方向の新しい決定を促し得るところまで、すなわち、一旦認識された方向が、音の大きさのより強い増大によって方向の再決定が可能にされるまで、音源の方向として維持される。２個のラウドスピーカまたはラウドスピーカアレイの間の正確に中間の聴取位置においては、中心のローカライゼーション、強いローカライゼーションが焦点に集まり、従って、対称的な周囲知覚が自動的に実現される。この考察は、信号が左側ステレオチャネルと右側ステレオチャネルとの間に同じレベルおよび同じ遅延を有する状態でいつも放出されることを想定する。
【００３０】
しかしながら、車両の客室における通常の聴取位置が一般にそうであるように、聴取位置がこの対称軸外にある場合には、ローカライゼーションの所望の品質はもはや、単にレベルをイコライズすることだけでは達成され得ない。ラウドスピーカの左側ステレオチャネルおよび右側ステレオチャネルの信号の振幅を適応させて、ラウドスピーカの放出角の差を補償することでさえも、ステレオラウドスピーカの間の対称軸上の聴取位置に対応する結果を達成することができない。
【００３１】
位相をそれぞれ調整するいかなる方法においても、遅延における信号の差が対称性を欠く着席位置によって変わることは、簡単な測定によって示され得る。耳を表す２個のマイクロホンを含むダミーヘッドを位置決めすることによって、客室内の聴取者の生理機能を、車両内に配列されたラウドスピーカの間の長手方向の中心線において正確にシミュレートし、両耳の位相差を測定することは、両方のステレオ信号が非常に高い程度に一致することを示す。約１５００Ｈｚまでの音響心理学的に関連する領域における対応する測定結果は、図１から明らかである。
【００３２】
ここで、図１を参照すると、ダミーヘッドのマイクロホンによって測定される場合の信号の位相差の曲線が例示され、左側の測定信号と右側の測定信号との間の位相差を、対数目盛の周波数の関数として度の単位で示している。実際の車両の客室内部で測定された場合の例から、２つの測定された信号の位相差が、１００Ｈｚ未満の周波数に対しては比較的わずかであり、正の方向および負の方向において４５度を超えないことは明らかである。
【００３３】
ここで、図２を参照すると、運転者の場所で位置決めされたダミーヘッドのマイクロホンによって測定される場合の信号の位相差の曲線が例示され、再び左側の測定信号と右側の測定信号との間の位相差を、対数目盛の周波数の関数として度の単位で示している。図２から、この場合において、測定された２つの信号の位相差が、１００Ｈｚ超の周波数に対して正の方向および負の方向において４５度をすでに超えていることは明らかである。３００Ｈｚ超の周波数において、位相差は、１８０度と同じ高さに現れている。図１と図２との測定結果を比較すると、ラウドスピーカの間の対称軸外の聴取位置は、この場合の運転者の位置がそうであるように、左耳および右耳に到着する信号の実質的により大きな位相差、関係するオーディオ信号のローカライゼーションの大きな損失という結果をもたらすことは明らかである。
【００３４】
動力車両のオーディオシステムを手動で調整する過程において、位相を微調整（すなわち、調整）するための上記のような手段のすべてが、良い音響を達成するために、いわゆる「ステージ」を配置し、構成するように用いられる。対照的に、周波数応答の大きさをイコライズすることは、いわゆる調性を最適化することにもっぱら役立つ。これらの目的はまた、ここに記載されるような方法に関係し、すなわち、周波数応答の大きさをイコライズすることに関して任意に所定の標的の機能を達成することに関係する。位相イコライゼーションの方法を重点的に取り扱うことは、さらに、車両内のすべての可能性のある聴取位置において、対称性および広がりをステージに与えることを向上させること、および現実的なステージ幅を維持しながらローカライゼーションの正確さを向上させることに役立つ。
【００３５】
他の研究グループは、様々なラウドスピーカの異種の位相調整によって発生させられた櫛形フィルタ効果を、調査された測定点において低減するために位相を用いることにより、室内のこの位置において、スペクトル的により接近した、従って、向上した周波数応答の大きさを発生させる。このことに加えて、最適な周波数応答の大きさが原理的にローカライゼーションの品質について結論を可能にしないので、これはたぶん主要な所望の目標ではなかったが、ローカライゼーションがまた改善され得るのはこの様態においてである。
【００３６】
位相イコライゼーションに対して従来から公知のアプローチから、所望の位相周波数応答を単に複製するこの目的のために設計されたＦＩＲオールパスフィルタは、その位相に全く影響しないが、変動する大きさの狭帯域グリッチを主として含む、周波数応答の大きさに対して特定の強い影響を有することもまた明らかである。このことは別にして、上記の観点によって生成される位相イコライザは、サウンドの知覚をだめにし得る、長いインパルス応答によって特徴付けられる。位相イコライゼーションにおいてインパルス応答を試験することは、調性の妨害と、位相イコライザの群遅延が設計される仕方との間に直接的つながりがあることを示した。
【００３７】
例えば、位相イコライザの群遅延の狭いスペクトル帯域における大きく、かつ突然の変化は、ちょうどこれらの周波数において高いＱファクタ／利得を有するフィルタに類似して、インパルス応答内の振動を結果としてもたらす。この効果は、「時間的拡散」とも呼ばれ、すなわち、調性の妨害は、ますますさらに長く持続し、従って、ますます有害であり、狭いスペクトル帯域における偏差はますます動的である。代わりに、位相イコライザの群遅延の突然の変化が非常に低い周波数帯域内にある場合には、これは、ほとんど有害ではないとして経験され、ほとんどの場合に無視して良いことさえある。しかしながら、ともかく、この文脈は、オーディオシステムの衝撃性をだめにしないように、聴覚指向のスムージングによって、例えば、位相イコライザを設計する場合に考慮される必要がある。換言すると、良い衝撃性のために、位相イコライザの群遅延は、より高い周波数に対して低減された動的応答を有する必要がある。
【００３８】
位相イコライゼーションフィルタ用のフィルタに加えて、大きさイコライゼーション用のフィルタもまた、オーディオシステムの衝撃性に影響を及ぼす。ここで、位相イコライゼーション用のフィルタ、すなわち位相イコライザの設計におけるように、聴覚指向の非線形で、複雑なスムージングが使用される。この上に、衝撃性が影響される仕方はまた、大きさイコライゼーション用のフィルタの設計に依存する。すなわち、周波数応答の大きさの所定の所望の曲線が、線形に変換されるかまたは最小の位相調整が行われるか否かに依存して、妨害は多くなるかまたは少なくなる。
【００３９】
これは、良い衝撃性に対して、大きさイコライゼーションを実現するために最小位相フィルタの独占的使用が推奨される理由であるが、最小位相フィルタは、位相イコライゼーションを実装することにおいて考慮される必要があるように、特定の最小位相応答を特徴とする。これは、まさに同じことが位相に影響を及ぼす他の構成要素（例えばディレイライン、クロスオーバフィルタなど）に対してあてはまる。このことに加えて、周波数応答の大きさをイコライズするために最小位相フィルタを用いることは、線形の位相設計に比べて、同じ所望の周波数応答の大きさを達成するために、それらがほんの半分の個数のフィルタ係数を必要とし、従って非常に効率的に実現される利点を有する。
【００４０】
以下は、位相応答を周波数の関数としてイコライズする仕方が、ローカライゼーションの著しい向上を結果としてもたらすように設計され得るかを記述する。この目的のために、対応する前の考察および実行される場合の試験がここで詳細になる。
【００４１】
基本的に、３つの因子が水平ローカライゼーションに責任があり、３つの因子とは、すなわち、上記のハース効果または先行音効果（第一波面の法則とも呼ばれる）と、両耳間の時間差（ＩＴＤ）と、両耳間のレベル差（ＩＬＤ）とである。先行音効果は、残響環境において主として有効であり、両耳間の時間差は、Ｂｌａｕｅｒｔに従うと最高約１５００Ｈｚまでのより低いスペクトル帯域において有効であり、両耳間のレベル差は、主として約４０００Ｈｚ超において有効である。
【００４２】
しかしながら、本システムによって考慮されるローカライゼーションに対して、関心のあるスペクトルの範囲は、聴取者によって知覚されるようなローカライゼーションを分析したり、修正したりする場合に、両耳間の時間差（ＩＴＤ）だけが考慮される必要のある最高約１５００Ｈｚまでの可聴周波数範囲であることが判明した。
【００４３】
従って、車両内部のすべての着席位置における各々のラウドスピーカの両耳室内インパルス応答（ＢＲＩＲ）が記録された。この目的のために、人工的な頭部（「人の頭部の耳がある位置においてマイクロホンを備えるダミーヘッド」）が、マネキンに取り付けられ、これに加えて、車両の客室内のすべての残りの座席が、所望の調整のタイプ（すなわち、運転者最適化調整、前部最適化調整、後部最適化調整、またはすべての位置に対する最適化調整）に依存して、乗客および／またはマネキンによって占められ得るかまたは占められないことがあり得る。
【００４４】
ここで、図３を参照すると、例としてオーディオシステムのラウドスピーカの配列と共に、概略的に示される車両内部１においてダミーヘッドの助けによって試験される場合のすべての位置が上面図で例示される。車両内部は、以下のラウドスピーカを備えているオーディオシステムを特徴とする。すなわち、前部左ラウドスピーカ２、前部中央ラウドスピーカ３、前部右ラウドスピーカ４、側部左ラウドスピーカ５、側部右ラウドスピーカ６、後部左ラウドスピーカ７、後部中央サブウーハ８ならびに後部右ラウドスピーカ９である。また図３から明らかなことは、ダミーヘッドの助けによってＢＲＩＲを測定する場所であり、これらの場所は、１０ａおよび１１ａ（運転者および前部乗客の前部着席位置）、１０ｂおよび１１ｂ（運転者および前部乗客の中央着席位置）、ならびに１０ｃおよび１１ｃ（運転者および前部乗客の後部着席位置）である。さらに明らかなことは、測定位置１２（後部左着席位置）および１３（後部右着席位置）である。
【００４５】
ここで、図４を参照すると、ダミーヘッドの助けによって車両客室１４において試験されるすべての測定位置の側面図が例示される。車両の２つの前部着席場所において、ダミーヘッド信号の各々は、座席の位置を決めることによって３つの位置（前部、中央、後部）において測定され、座席の位置を決めることは、結果として、客室において前部左に配列された運転者の座席および前部右に配列された前部乗客の着席位置、測定位置１０ａおよび１１ａ（運転者および前部乗客の前部着席位置）、１０ｂおよび１１ｂ（運転者および前部乗客の中央着席位置）ならびに１０ｃおよび１１ｃ（運転者および前部乗客の後部着席位置）となった。この状況において、着席位置における前部から後部への移動の他に、同時に下部から上部への移動が、小柄な人、普通の大きさの人および背の高い人を考慮して行われる。
【００４６】
図４に示される側面図からさらに明らかなことは、ダミーヘッドの高さが２つの後部着席場所（左手および右手、図３を参照されたい）において調整される仕方であり、２つの後部着席場所の各々は、信号を測定する、左および右の後部着席場所に対する３つの位置、すなわち、左および右の後部着席場所に対する位置１２ａおよび１３ａ、左および右の後部着席場所に対する中央位置１２ｂおよび１３ｂ、ならびに左および右の後部着席場所に対する下方位置１２ｃおよび１３ｃに設定される。各場合における最高位置と最低位置との間のダミーヘッド配列の高さの変更は、これも異なる大きさの人を考慮して行われる。この配列の意図は、体の大きさの差、従って車両の客室における実際の乗客の耳に関する聴覚位置における差を再現することである。
【００４７】
前部着席位置における水平ローカライゼーションに対して、前部ラウドスピーカ２、４、およびオプションで３のみが関連する。同様に、前部ラウドスピーカ２、３、および使用可能な場合３の他に後部着席位置における水平ローカライゼーションに対して、後部ラウドスピーカ７、９、ならびに側部ラウドスピーカ５および６が関連する。しかしながら、着席位置が環境（すなわち、客室）ならびにそこにあるラウドスピーカの配列に依存するローカライゼーションにどれかのラウドスピーカが関連する。換言すると、各着席位置（および従って聴取位置）に対して、ラウドスピーカの規定された群が考えられ、ラウドスピーカの各群は少なくとも２個の単一のラウドスピーカを備えている。
【００４８】
（関連する群から選ばれる）聴取位置およびラウドスピーカの各対に対して両耳室内インパルス応答（ＢＲＩＲ）を測定した後に、さらなる分析およびフィルタ統合がオフラインで実行され得る。位相を調整するためのすべての手段を考慮に入れる際に、考えられる聴取位置に関連する群の対応するラウドスピーカを重ねることは、重なりスペクトルの所望の位相周波数応答を生成する。
【００４９】
２つの前部座席の聴取位置に対する両耳間の時間差（ＩＴＤ）を最適化することは、１°のステップで０から１８０°までの位相シフトをラウドスピーカの関連する群のラウドスピーカのうちの１つに特定の周波数で供給されるオーディオ信号に賦課することによって、実行される。すなわち、特定の周波数ｆｍのオーディオ信号は、前部聴取位置、例えばラウドスピーカ２および４（中央スピーカ３が存在しない場合）に割り当てられる群のラウドスピーカに供給される。次いで後に、０°から１８０°までの位相シフトφｎがラウドスピーカ２（または、代わりにラウドスピーカ４）に供給され、このとき、他のラウドスピーカに供給されるオーディオ信号の位相は、不変のままである。このことは、例えば１００Ｈｚ〜１５００Ｈｚなどの所与の周波数範囲における異なる周波数に対してなされる。上記のように、１５００Ｈｚ未満の周波数範囲は、車両の客室などの反響する環境において水平ローカライゼーションに対して主として決定的である。
【００５０】
各考えられる聴取位置に対する測定された両耳室内インパルス応答（ＢＲＩＲ）を用いて、結果として生じる位相差Δφｍｎは、周波数ｆｍと位相シフトφｎとの各対に対して計算され得る。位相差Δφｍｎは、ダミーヘッドの２個のマイクロホン（すなわち、「耳」）に存在する音響信号の位相差であるか、または、換言すると、考えられる聴取位置に配置されるダミーヘッドの「耳」に存在する、結果として生じる音響信号から計算される重なりスペクトルの位相である。
【００５１】
本例において、左前部ラウドスピーカ２の信号の位相は変化させられが、但し、代案として、右ラウドスピーカ４の信号は同様に変化させられ得る。それに続いて、対象とする全スペクトルバンドにおいて重なりスペクトルの、結果として生じる位相Δφｍｎを得ること、次いで、結果を行列の中に入れることがなされる。複数のラウドスピーカが試験される特定の動力車両のサウンドシステムに存在する場合、３個以上のラウドスピーカの信号もまた、考えられる聴取位置に対する最適の結果を達成するために変化させられるように選ばれ得る。この場合、結果は、位相差の三次元「行列」である。しかしながら、ものごとを複雑にすることを避けるために、さらなる考察は、２個のラウドスピーカ（例えば、前部ラウドスピーカ３および４)のみを備えているラウドスピーカの群に限定し、その結果、１個のラウドスピーカの音声信号のみが位相シフトされる必要がある。
【００５２】
位相シフトを挿入し、結果として生じる位相差Δφｍｎを計算する手順は、割り当てられた関連するラウドスピーカの同じ群を有する各聴取位置に対してなされ得る。本例において、前部ラウドスピーカ２および４を備えている群が考えられた。ラウドスピーカのこの群は、車両の前部の６つの聴取位置（運転者位置：前方、中央、後方、前部乗客位置：前方、中央、後方）に割り当てられる。その結果、上記の手順を用いて、６つの行列Δφｍｎが計算され得、各行列は特定の聴取位置に属する。
【００５３】
以下の最適化に関して、各聴取位置に対して計算された位相差Δφｍｎは平均され得、平均位相差ｍΔφｍｎの行列を得る。従って、平均位相差ｍΔφｍｎの最適化は、達成され得、すべての考えられる聴取位置において良好なローカライゼーションの理由となり得る。
【００５４】
ここで、図５を参照すると、２個の前部測定位置１０および１１（例えば、中心位置１０ｂ、１１ｂ）の上の重なりスペクトルの位相の形式において上記に得られるような結果の三次元表現ｍΔφｍｎが例示され、２個の前部測定位置１０および１１において、設定された位相シフトφｎは０〜１８０°のｙ軸に入り、一方、各々の場合において、ｚ軸は重なりスペクトルの平均位相差ｍΔφｍｎをプロットし、ｘ軸は対応する周波数ｆｍをプロットする。この三次元表現における最小限の高さの線は、対応する着席位置または複数のそれぞれの着席位置に対する最小限の両耳間の時間差の意味において最適な位相シフトに対応する。位相差ｍΔφｍｎのＮ×Ｎ行列を仮定すると（ここで、周波数指数ｍは０からＭ−１に動き、位相指数ｎは０からＮ−１に動く）、周波数ｆｍにおいて最適なシフトφｘ（ｆｍ）を生む指数Ｘは、以下の関係から導かれ得る。
【００５５】
ｍΔφｍＸ＝ｍｉｎ｛ｍΔφｍｎ｝、ここでｎ＝０，１，．．．，Ｎ−１、
ここで、上記に論議される例において、Ｎ＝１８０、すなわちφｎ＝ｎ°、ここでｎ＝０、１、．．．、１７９である。例をあげると、周波数値数Ｍは、Ｍ＝１５００、すなわちｆｍ＝ｍ  Ｈｚから選択され得、ここでｍ＝１，２，．．．，１５００である。あるいは、周波数値ｆｍの代わりに対数間隔が選ばれ得る。最適な位相シフトは、最小限の位相差という結果をもたらす。
【００５６】
ここで、図６を参照すると、図５に示されるような三次元表現の上面図が例示され、図６において、横座標は、Ｈｚ単位で測定周波数ｆｍをプロットし、一方、縦座標は、ここでは左ラウドスピーカ２（図３を参照されたい）であるラウドスピーカのオーディオ信号に賦課される位相シフトΔφｎをプロットする。この上面図に重ねられるのは、位相差、従って図５に示されるように三次元表現ｍΔφｍｎからの最小値として得られる両耳間の時間差（ＩＴＤ）に対する最小限の高さの「線」（例えば、ｆｍの関数としての最適な位相シフトφｘ）である。
【００５７】
ここで、図７を参照すると、測定された結果の三次元表現から分離された、上面図における最小限の「高さ」の線（すなわち、最小限の位相差であり、図６も参照されたい）がより明快に例示される。ここで再び、横座標は、Ｈｚ単位で周波数ｆｍをプロットし、一方、縦座標は、左ラウドスピーカ２（図３を参照されたい）に対する対応する位相シフトφｎをプロットする。図７から明らかな曲線は、従って前部左チャネルに対する最適量としての（周波数依存の）最適な位相シフトφＸの曲線であり、重なりスペクトル位相の最大の最小化、および従って２個の前部着席場所の上の平均としての最適の水平ローカライゼーションの結果をもたらし、重なりスペクトル位相の最大の最小化および最適な水平ローカライゼーションの各々はまた、結果として生じる重なりスペクトルを計算するためにオプションで重み付けされ得る。結果は、図６および図７に示されるように、試験された２個の左および右の前部着席位置の同等の重み付けから得られる。しかし、計算において、より大きな重みを有する運転者の場所に入り、おそらく最も多く占有される着席位置であるものに対して両耳間の時間差を最適化することに、より多くの重みを与えることは正に可能である。
【００５８】
行列最小値を利用して、位相イコライザを直接形成することは、上記に説明されるように、非最適化の瞬間力を有するフィルタという結果をもたらすが、可能な限り最良のローカライゼーションを提供する。従って、このことは、最適なローカライゼーションと瞬間力ノイズ含有量との間の妥協を伴う。
【００５９】
この目的のために、行列最小値φＸ（ｆｍ）は、位相イコライゼーションフィルタが計算される前に、スライディング、非線形の複素スムージングフィルタの助けによってスムージングされる（複合フィルタリングの詳細については、Ｍｏｕｒｊｏｐｏｕｌｏｓ，Ｊｏｈｎ  Ｎ．；Ｈａｔｚｉａｎｔｏｎｉｏｕ，Ｐａｎａｇｉｏｔｉｓ  Ｄ．：“Ｒｅａｌ−Ｔｉｍｅ  Ｒｏｏｍ  Ｅｑｕａｌｉｚａｔｉｏｎ  Ｂａｓｅｄ  ｏｎ  Ｃｏｍｐｌｅｘ  Ｓｍｏｏｔｈｉｎｇ：Ｒｏｂｕｓｔｎｅｓｓ  Ｒｅｓｕｌｔｓ”，ＡＥＳ  Ｐａｐｅｒ  ６０７０，ＡＥＳ  Ｃｏｎｖｅｎｔｉｏｎ  １１６，２００４年５月を参照されたい）。このことは、一方において、車両内のその後の聴取試験によって確認されたように、ローカライゼーションの精度が今までと同じく良好のままであることを確実にし、他方において、再び車両内のその後の聴取試験によって確認されたように、位相イコライザの瞬間力をそれがもはや迷惑物として感じられない点まで高める。
【００６０】
スムージングされた最適位相関数φＸ，ＦＩＬＴ（ｆｍ）は、考察中のラウドスピーカ（上記に論議される例において前部左ラウドスピーカ２）に供給されるオーディオ信号の位相をイコライズする位相イコライザの設計のための基準（設計目標）として用いられる。イコライジングフィルタは、例えば、ＦＩＲフィルタまたはＩＩＲフィルタなどの任意のデジタルフィルタ技術によって実装され得る。
【００６１】
ここで、図８を参照すると、非線形の複素スムージングフィルタの適用後の位相イコライザの、結果として生じる群遅延と、Ｈｚ単位にて対数で周波数ｆｍをプロットする横座標と、周波数の関数として位相イコライザφＸ，ＦＩＬＴ（ｆｍ）の対応する群遅延をプロットする縦座標とが例示される。図８から見られ得るように、この場合の群遅延の動的応答が少なければ少ないほど、周波数は高い。既に上記に説明されたように、これは利点である。なぜなら、このような方法で時間的拡散が実質的に防げるからである。
【００６２】
ここで、図９を参照すると、前部左チャネル（図３に示されるようなラウドスピーカ２）の得られたＦＩＲ位相イコライザのインパルス応答が例として示される。図９の下図は時間の関数としてのインパルス応答の大きさの線形の表現を例示し、図９の上図は時間の関数としてのインパルス応答の大きさの対数の表現を例示する。
【００６３】
ここで、図１０を参照すると、ＦＩＲフィルタとして実装された、図９に示されるような位相イコライザφＸ，ＦＩＬＴ（ｆｍ）のボードダイアグラムが例示され、上と下の両図における横座標は対数目盛りで周波数をプロットし、図１０の下図の縦座標はｄＢ単位でレベルをプロットし、図１０の上図の縦座標は位相をプロットする。
【００６４】
このような方法で実現されるような位相イコライザは、次いで左前部ラウドスピーカ２（図３を参照されたい）の信号に適用された。本例において、全部の手順は、関連する群の他のラウドスピーカ、すなわちラウドスピーカ３および４（図３を参照されたい）に対して行われる。これらのラウドスピーカまたはそれらのそれぞれの活性化信号（ラウドスピーカに供給されるオーディオ信号）のために、対応する位相イコライザが、測定されたＢＲＩＲおよびその後の信号処理の結果から得られた。前部ラウドスピーカの位相イコライゼーションに対する最適曲線を得、それを適用した後に、後部着席位置のために着席位置最適化がさらに行われた。この目的のために、オーディオ信号のローカライゼーションは、側部左および側部右にそれぞれ配列された（図３を参照されたい）ラウドスピーカ５および６を利用する前部着席位置に関して説明された方法と同じ方法で最適化された。ダミーヘッドが対応して位置決めされた方法は、図３および図４に示される（位置１２ａ、１３ａ、１２ｂ、１３ｂ、１２ｃおよび１３ｃ）。
【００６５】
オーディオ信号のローカライゼーションは、ここで、時間的拡散の迷惑物を作ることなく、位相イコライザによる周波数応答の大きさにおける望まない変化に我慢する必要もなく、説明された方法で動力車両の客室内の４つのすべての着席位置において相当に向上し得る。
【００６６】
ここで、図１１ａ〜図１１ｄを参照すると、最適化（位相イコライザを挿入すること、とりわけすべての位相イコライズされたチャネルに対する位相関数φＸ，ＦＩＬＴ（ｆｍ））の後に、車両内の４つすべての着席位置１０、１１、１２および１３において測定された両耳の重なりスペクトルの位相周波数応答が位相イコライザの適用前に測定された両耳の重なりスペクトルの位相周波数応答と比較する仕方が例示される。図１１ａ〜図１１ｄの横座標はＨｚ単位にて対数表現で周波数をプロットし、縦座標は両耳の位相差曲線を度でプロットする。図１１ａは、車両内の左前部着席位置に対する最適化の前および後の両耳の位相差周波数応答である。図１１ｂにおいて、車両内の右前部着席位置に対する最適化の前および後の両耳の位相差周波数応答が対応して比較される。図１１ｃにおいて車両内の左後部着席位置に対する最適化の前および後の両耳の位相差周波数応答が比較され、そして図１１ｄにおいて車両内の右後部着席位置に対する最適化の前および後の両耳の位相差周波数応答が比較される。最適化の前に得られた周波数依存の両耳の位相差は、図において「Ａ」と識別され、最適化の後に得られた位相差は「Ｂ」と各々識別される。図１１ａ〜図１１ｄから明らかなことは、理想の０ラインから位相周波数応答のより少ない偏移が、特に、車両内の４つのすべての着席位置に対してより低い周波数で達成可能であり、その結果、すべての着席位置に対して車両オーディオシステム内のローカライゼーションの実質的な向上がもたらされることである。
【００６７】
本方法は、聴取室１内の少なくとも１つの聴取位置（例えば、運転者中央位置１０ｂ）における音響ローカライゼーションを最適化するために有効に用いられ得る。音場は少なくとも１つの聴取位置に割り当てられるラウドスピーカの群（例えば、前部ラウドスピーカ２および４）によって生成され、この場合、ラウドスピーカの群は、第１のラウドスピーカ（例えば、前部左ラウドスピーカ２）と、少なくとも第２のラウドスピーカ（例えば、前部右ラウドスピーカ４、およびオプションで中央ラウドスピーカ３）とを備えている。各ラウドスピーカは、オーディオチャネルを介してオーディオ信号が供給される。本方法は、第２のラウドスピーカ４に供給する少なくともオーディオチャネルのための位相イコライゼーションフィルタのフィルタ係数を計算することを包含する。位相イコライゼーションフィルタの位相応答は、少なくとも１つ聴取位置１０における両耳の位相差Δφｍｎ、または代わりに複数の聴取位置が考えられる場合、複数の聴取位置（例えば、前部位置１０ｂ、１１ｂ）の平均された平均的な両耳の位相差ｍΔφｍｎが所定の周波数範囲内で最小化されるように、設計される。方法は、位相イコライゼーションフィルタをそれぞれのオーディオチャネルに適用するステップをさらに包含する。
【００６８】
上記のように、それぞれの聴取位置（例えば、前部左位置１０および前部右位置１１、図３を参照されたい）において、１人または１人より多い聴取者によって知覚される両耳間の時間差は、本方法によって最小化され得る。各考えられる聴取位置１０、１１において位相イコライゼーションフィルタを第１に計算するステップを実行するために、両耳の伝達特性は、考えられる聴取位置１０、１１に割り当てられる群の各ラウドスピーカ２、４に対して決定され得る。これは、例えば上記のようにダミーヘッドを用いる測定によって達成され得る。
【００６９】
最適化は所定の周波数範囲内で実行され得、従って、所定の周波数範囲内から選ばれた一組の周波数ｆｍは、所定の位相範囲（例えば、φｎ＝｛１°，２°，．．．，１８０°｝）から選ばれた一組の位相シフトφｎと同様に定義される。
【００７０】
両耳の位相差Δφｍｎは、各考えられる聴取位置１０、１１において計算され得、それによって、計算は、一組の周波数の各周波数ｆｍおよび一組の位相シフトの各位相シフトφｎに対して行われる。それによって、計算（これはシミュレーションとも呼ばれ得る）の目的のために、オーディオ信号が各ラウドスピーカ２、４に供給されることが仮定され、それによって、少なくとも１個の第２のラウドスピーカ４に供給されるオーディオ信号は、第１のラウドスピーカ２に供給されるオーディオ信号に対して位相シフトφｎだけ位相シフトされる。各考えられる聴取位置１０、１１に対する両耳の位相差Δφｍｎのアレイは、このように生成される。Ｍの異なる周波数値ｆｍおよびＮの異なる位相シフトφｎの場合、結果として生じる行列は、ラウドスピーカの関連する群が２個のラウドスピーカを備えている場合、Ｍ×Ｎ行列である。３つのラウドスピーカ（例えば、追加の中央ラウドスピーカ３、図３を参照されたい）に関して、結果として生じる行列は、同じ組の位相シフトφｎが第２のラウドスピーカ３および第３のラウドスピーカ４に供給されるオーディオ信号に適用される場合、Ｍ×Ｎ×Ｎ構成要素を備えている三次元アレイである。
【００７１】
すべての考えられる聴取位置において、向上したローカライゼーションを達成するために、考えられる聴取位置１０、１１における両耳の位相差Δφｍｎの重み付けされた平均である平均的な両耳の位相差ｍΔφｍｎのアレイが計算され得る。重み付け係数は、０もしくは１または区間［０，１］であり得る。しかしながら、１つのみの聴取位置（例えば、運転者位置１０）が考えられる場合、運転者位置１０における両耳の位相差ΔφｍｎのそれぞれのアレイがアレイｍΔφｍｎとして用いられ得る。
【００７２】
実際の最適化は、平均的な両耳の位相差ｍΔφｍｎのアレイにおいて、少なくとも１個の第２のラウドスピーカ４に出力されるオーディオ信号に適用される各周波数ｆｍに対して最適位相シフトφＸを探索することによって、実行され得る。最適の位相シフトφＸは、平均的な両耳の位相差ｍΔφｍｎの最小値を生み出すように定義される。従って、位相関数φＸ，ＦＩＬＴ（ｆｍ）は、周波数ｆｍの関数として最適な位相シフトφＸを表す少なくとも１個の第２のラウドスピーカに対して得られ得る。さらなるラウドスピーカが考えられる場合（例えば、第３の中央ラウドスピーカ３、図３を参照されたい）、最適な位相シフトφＸは、第２のラウドスピーカおよび各さらなるラウドスピーカ３、４に供給されるオーディオ信号に対して最適な位相シフトを含むベクトルである。
【００７３】
両耳の位相差Δφｍｎは、各聴取位置に存在する音響信号の重なりスプペクトルの位相である。これらの重なりスペクトルは、ラウドスピーカの関連する群のラウドスピーカおよび以前に測定された対応するＢＲＩＲに供給されるオーディオ信号を考えることによって容易に計算（すなわち、シミュレート）され得る。
【００７４】
本方法は、測定された両耳室内インパルス応答（ＢＲＩＲ）用いて、音響信号をシミュレートし、その音響信号は、計算において仮定されるように、オーディオ信号がすべての関連するラウドスピーカに供給され、位相シフトが少なくとも１個の第２のラウドスピーカの供給チャネルに挿入される場合に存在する。各聴取位置におけるシミュレートされた（両耳の）信号から、対応する両耳間の位相差が引き出され得る。しかしながら、このシミュレーションは、現実の測定値によって置き換えられ得る。すなわち、シミュレーションにおいて言及されたオーディオ信号は、ラウドスピーカに実際に供給され得、聴取位置において、結果として生じる音響信号は、両耳で測定され得る。所望の両耳間の位相差は、シミュレートされた信号から引き出される方法と同じ方法で測定された信号から引き出され得、従って、結果として、シミュレーションに基づく「オフライン」方法に関して上記に論議されるように両耳間の位相差の行列と同じ行列がもたらされる。この両耳間の位相差の行列は、両方の場合において同じ方法で処理される。しかしながら、後者の場合、ラウドスピーカから放出されるオーディオ信号の周波数および位相は実際には変化させられ、第１の場合において、このことは、シミュレーション中にコンピュータにおいてのみ行われる。
【００７５】
本発明を実現する様々な例が開示されたが、本発明の精神および範囲から逸脱することなく、本発明の利点のいくつかを達成する様々な変更および修正が行われ得ることは当業者に明らかである。同じ機能を果たす他の構成要素が適切に代わりに用いられ得ることは当業者に明らかである。本発明の概念に対するそのような修正は、添付の特許請求の範囲に含まれることが意図される。さらに本発明の範囲は、自動車の用途に限定されるのではなく、例えば、ホームシネマまたは類似のものならびにシネマおよびコンサートホールまたは類似のものの場合などの任意の他の環境においても適用され得る。
【符号の説明】
【００７６】
１  車両内部
２  前部左ラウドスピーカ
３  前部中央ラウドスピーカ
４  前部右ラウドスピーカ
５  側部左ラウドスピーカ
６  側部右ラウドスピーカ
７  後部左ラウドスピーカ
８  後部中央サブウーハ
９  後部右ラウドスピーカ
【図１】
2011097561.tif 000006
【図２】
2011097561.tif 000007
【図３】
2011097561.tif 000008
【図４】
2011097561.tif 000009
【図５】
2011097561.tif 000010
【図６】
2011097561.tif 000011
【図７】
2011097561.tif 000012
【図８】
2011097561.tif 000013
【図９】
2011097561.tif 000014
【図１０】
2011097561.tif 000015
【図１１ａ】
2011097561.tif 000016
【図１１ｂ】
2011097561.tif 000017
【図１１ｃ】
2011097561.tif 000018
【図１１ｄ】
2011097561.tif 000019