JP2011182135A - 三次元音場再現システム - Google Patents
三次元音場再現システム Download PDFInfo
- Publication number
- JP2011182135A JP2011182135A JP2010043335A JP2010043335A JP2011182135A JP 2011182135 A JP2011182135 A JP 2011182135A JP 2010043335 A JP2010043335 A JP 2010043335A JP 2010043335 A JP2010043335 A JP 2010043335A JP 2011182135 A JP2011182135 A JP 2011182135A
- Authority
- JP
- Japan
- Prior art keywords
- sound field
- speaker
- evaluation
- evaluation value
- reproduction system
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Stereophonic System (AREA)
Abstract
【課題】境界音場制御の原理を用いたBoSC再生システムにおいて、できるだけ再現音場の質を低下させることなく、スピーカ数を減らすことのできる方法を提供する。
【解決手段】境界音場制御(BoSC)理論に基づき構成された、聴者の頭部を囲う複数のスピーカから成る三次元音場再現システムにおいて、a)所定の基準スピーカと各評価対象スピーカの間のグラムシュミット直交化評価値を計算し、b)グラムシュミット直交化評価値の高い順に評価対象スピーカを選択し、c)基準スピーカと選択された評価対象スピーカから成るサブシステムと未選択評価対象スピーカの間のグラムシュミット直交化評価値を計算し、d)前記b)及びc)のステップを繰り返し、該グラムシュミット直交化評価値が所定の閾値を下回った時点で選択を停止する。この方法により、一般のA/Dコンバータを用いて簡単に構成することのできる24ch二次音源システムを決定した。
【選択図】図3
【解決手段】境界音場制御(BoSC)理論に基づき構成された、聴者の頭部を囲う複数のスピーカから成る三次元音場再現システムにおいて、a)所定の基準スピーカと各評価対象スピーカの間のグラムシュミット直交化評価値を計算し、b)グラムシュミット直交化評価値の高い順に評価対象スピーカを選択し、c)基準スピーカと選択された評価対象スピーカから成るサブシステムと未選択評価対象スピーカの間のグラムシュミット直交化評価値を計算し、d)前記b)及びc)のステップを繰り返し、該グラムシュミット直交化評価値が所定の閾値を下回った時点で選択を停止する。この方法により、一般のA/Dコンバータを用いて簡単に構成することのできる24ch二次音源システムを決定した。
【選択図】図3
Description
本発明は、原音場を物理的に忠実に記録・再現することのできる三次元音場再現システムである境界音場制御(Boundary Surface Control: BoSC)再生システムの音源配置の最適化技術に関する。
近年、5.1chサラウンドシステムなどの一般家庭への普及により、臨場感のある音による仮想空間を手軽に体験することが可能になった。しかし、より高度なヴァーチャル・リアリティの実現には、広い空間の音場をあるがままに記録し、再現することができる音場再現技術が必要不可欠である。
本願発明者は、境界音場制御(BoSC)の原理を用いることにより、原音場を物理的に忠実に記録・再現することのできる三次元音場再現システムとして没入型聴覚ディスプレイシステム(BoSC再生システム)を開発した。
BoSCの原理によれば、その領域を囲む境界面上の音圧と音圧勾配を制御することにより、3次元音場内の任意の領域内の音圧を制御することができる。ここで制御を厳密に行うためには、境界面における音圧・音圧勾配制御点を波長よりも十分短く離散化しておく必要がある。上記BoSC再生システムは、コンパクト化を目的として、境界面を聴者の頭部周辺とし、その内部での忠実な音場再現を目的としたものであり、64ch(チャンネル)の収音・再生が可能となっている(非特許文献1)。
BoSC再生システムは、システムの利用者(聴者)の頭部周辺の領域において音場を再現する4層の架台から成る楕円形のドーム部と、該ドーム部を支える4本の柱部で構成される。楕円形のドーム部の各層には合計62個のフルレンジスピーカが配置され、柱部には各2個のサブウーファースピーカが配置されている。境界面を構成するマイクロホンアレーはC80フラーレンを基に設計された構造を持ち、70個のマイクロホンより構成される。
伊勢史郎,「キルヒホッフ−ヘルムホルツ積分方程式と逆システム理論に基づく音場制御の原理」,音響誌, Vol. 53(1997), pp. 706-713
岡田耕介, 川口孝幸, 榎本成悟, 伊勢史郎,「コンパクトな没入型聴覚ディスプレイの試作と評価」,日本音響学会誌,第62巻(2006)第1号, pp. 32-41
上記BoSC再生システムは、境界面を聴者の頭部周辺とするものの、全体としては1.5帖(約2.5m2)の簡易防音室に62.8chのスピーカアレイを実装した、比較的大規模なものになっている。
本発明が解決しようとする課題は、BoSC再生システムにおいて、できるだけ再現音場の質を低下させることなく、スピーカ数を減らすことのできる方法、及び、そのような方法を用いた結果作製されたBoSC再生システムを提供することである。
本発明者は、実用化の鍵となるBoSC再生システムの簡素化・可搬化を目指し、音源数を減らすことによる音場再現への物理的影響及び心理的影響を調べた。その中で、グラムシュミット(Gram-Schmidt)の直交化法を応用した制御音源配置の最適化手法とその変法により、BoSC再生システムに用いる二次音源の設置数とその配置を検討し、再現精度及び音源定位実験により評価を行うことにより、本発明を想到するに至った。
すなわち、上記課題を解決するために成された本発明は、境界音場制御理論に基づき構成された、聴者の頭部を囲う複数のスピーカから成る三次元音場再現システムにおいて、
a) 所定の基準スピーカと各評価対象スピーカの間のグラムシュミット直交化評価値を計算し、
b) グラムシュミット直交化評価値の高い順に評価対象スピーカを選択し、
c) 基準スピーカと選択された評価対象スピーカから成るサブシステムと未選択評価対象スピーカの間のグラムシュミット直交化評価値を計算し、
d) 前記b)及びc)のステップを繰り返し、該グラムシュミット直交化評価値が所定の閾値を下回った時点で選択を停止する
ことを特徴とする。
a) 所定の基準スピーカと各評価対象スピーカの間のグラムシュミット直交化評価値を計算し、
b) グラムシュミット直交化評価値の高い順に評価対象スピーカを選択し、
c) 基準スピーカと選択された評価対象スピーカから成るサブシステムと未選択評価対象スピーカの間のグラムシュミット直交化評価値を計算し、
d) 前記b)及びc)のステップを繰り返し、該グラムシュミット直交化評価値が所定の閾値を下回った時点で選択を停止する
ことを特徴とする。
また、本発明に係る三次元音場再現システムは、境界音場制御理論に基づき構成された、聴者の頭部を囲う複数のスピーカから成る三次元音場再現システムにおいて、
a) 所定の基準スピーカと各評価対象スピーカの間のグラムシュミット直交化評価値を計算し、
b) グラムシュミット直交化評価値の高い順に評価対象スピーカを選択し、
c) 基準スピーカと選択された評価対象スピーカから成るサブシステムと未選択評価対象スピーカの間のグラムシュミット直交化評価値を計算し、
d) 前記b)及びc)のステップを繰り返し、該グラムシュミット直交化評価値が所定の閾値を下回った時点で選択を停止する
ことにより選択されたスピーカから成ることを特徴とする。
a) 所定の基準スピーカと各評価対象スピーカの間のグラムシュミット直交化評価値を計算し、
b) グラムシュミット直交化評価値の高い順に評価対象スピーカを選択し、
c) 基準スピーカと選択された評価対象スピーカから成るサブシステムと未選択評価対象スピーカの間のグラムシュミット直交化評価値を計算し、
d) 前記b)及びc)のステップを繰り返し、該グラムシュミット直交化評価値が所定の閾値を下回った時点で選択を停止する
ことにより選択されたスピーカから成ることを特徴とする。
BoSC理論によると、原音場で収音された境界面上の音圧と粒子速度を再現音場においてそれぞれ等しくなるように二次音源から信号を出力すれば、図1の左の原音場内の領域V内の音場が右の再現音場内の領域V'に忠実に再現される。
BoSC理論では、原理的に非常に多数のマイクロホンとスピーカを必要とするが、頭部伝達関数(Head-Related Transfer Function: HRTF)の個人差の問題や受聴者の身体動作との矛盾といった問題もなく、音源の近似や境界の形状の制約もないという特長を持つ。
本発明に係る三次元音場再現システムでは、このBoSC理論に基づくシステムの基本的特長をほとんど損なうことなく、スピーカの数を減らすことができる。従って、家庭内等で使用することのできる小規模システムでありながら、BoSC理論に沿った高忠実度の音場再生システムを構築することができる。
特に、一般家庭への普及を考えた場合、スピーカ(二次音源)の数の減少は、スピーカ自体のコストに加え、各chの信号処理のための機器、特にA/Dコンバータやアンプのコストを考慮することが重要である。その観点から、本発明に係る方法により、広く用いられている24ch音響機器を用いることのできる24chスピーカシステムを実現することにより、BoSC再生システムの一般への普及に大きな弾みがつくものと考えられる。
以下、本発明の実施形態を詳細に説明する。
1. はじめに
本発明の一つの実施形態であるBoSC音場再現システムは、図2(a)に示すような録音システムと、図2(b)に示すような再生システムから成る。図2(a)の録音システムは、音響空間を囲む仮想的な境界面を構成するマイクロホンアレイから成り、三次元音場をできるかぎり均一に収録すること、組立時や設置時の構造的な安定性とを考慮して、C80フラーレン(Fullerene)の構造を基に、底部の10個の頂点を切り取った70個の頂点を持つ、直径約46cmの球形に近い構造とした。各頂点の位置にはそれぞれ1個の無指向性マイクロホン(DPA4060-BM)を設置している。図2(b)の再生システムはマイクロホンと同数の70個のスピーカユニットから成り、その基台は、約1.2m×1.6m×2.0mの防音ボックス内に設けられた木製の4層の架台からなる楕円形のドーム部と、ドーム部を支える木製の4本の柱部から構成される。4層の各架台には上からそれぞれ6個、16個、24個、16個のフルレンジユニット(Fostex FE83E)が設置され、4本の柱部には低域を補うためのサブウーファーユニット(Fostex FW108N)が各2個ずつ設置されている。またドーム部の各架台と柱部の内部は空洞になっており、架台自体が密閉型エンクロージャの役割を果たす。また受聴者の耳の高さは,フルレンジスピーカの数が最大となる3層目の架台の高さを想定している。
本発明の一つの実施形態であるBoSC音場再現システムは、図2(a)に示すような録音システムと、図2(b)に示すような再生システムから成る。図2(a)の録音システムは、音響空間を囲む仮想的な境界面を構成するマイクロホンアレイから成り、三次元音場をできるかぎり均一に収録すること、組立時や設置時の構造的な安定性とを考慮して、C80フラーレン(Fullerene)の構造を基に、底部の10個の頂点を切り取った70個の頂点を持つ、直径約46cmの球形に近い構造とした。各頂点の位置にはそれぞれ1個の無指向性マイクロホン(DPA4060-BM)を設置している。図2(b)の再生システムはマイクロホンと同数の70個のスピーカユニットから成り、その基台は、約1.2m×1.6m×2.0mの防音ボックス内に設けられた木製の4層の架台からなる楕円形のドーム部と、ドーム部を支える木製の4本の柱部から構成される。4層の各架台には上からそれぞれ6個、16個、24個、16個のフルレンジユニット(Fostex FE83E)が設置され、4本の柱部には低域を補うためのサブウーファーユニット(Fostex FW108N)が各2個ずつ設置されている。またドーム部の各架台と柱部の内部は空洞になっており、架台自体が密閉型エンクロージャの役割を果たす。また受聴者の耳の高さは,フルレンジスピーカの数が最大となる3層目の架台の高さを想定している。
この再生システムにおける音場再生のためには、70個のスピーカアンプ及びオーディオインタフェース、再生用コンピュータとソフトウェアからなる大がかりな装置が必要であり、そのままでは一般家庭への普及には障害となる。一般家庭への普及も視野にいれた実用化には、システムのサイズ、形状、及びコストなどが課題となる。
このサイズ・形状に起因する導入への煩わしさやシステムの導入・維持にかかるコストは、音場再現に用いる二次音源数を減らし、その配置を単純化することで大幅に削減できる。そこで、音場再現に用いる二次音源数及びその配置について検討し、二次音源の減数及び配置の単純化によりBoSC音場再現システムの導入を容易にすることを目指し、BoSC音場再現システムに用いる二次音源配置の定量的評価法を検討した。
以下の説明では、まず、音場再現に用いる二次音源をグラムシュミットの直交化を用いて最適化する手法について述べ、この手法により選択された二次音源配置による音場再現への物理的影響及び心理的影響を計算機シミュレーションと主観評価実験により調べた結果について述べる。
2. 三次元音場再現システムの基礎
BoSC音場再現システムでは、収録に用いられるマイクロホンアレイとスピーカアレイのサイズ及び形状が異なるため、スピーカアレイから制御点までの伝達関数の逆フィルタ群を収録信号に畳み込んで出力することで、制御点において収録点と同じ音圧を再現する。
BoSC音場再現システムでは、収録に用いられるマイクロホンアレイとスピーカアレイのサイズ及び形状が異なるため、スピーカアレイから制御点までの伝達関数の逆フィルタ群を収録信号に畳み込んで出力することで、制御点において収録点と同じ音圧を再現する。
ここから、二次音源数M、制御点数Nの多チャンネル−多点制御逆システム(以下逆システム)を周波数領域で設計する手法について述べる。ここで、逆システムとはM×N個の逆フィルタ群の総称である。
二次音源iから制御点jまでの伝達関数をGji(ω)、入力信号をXj(ω)、観測信号をYj(ω)とすると、これらの関係は
Y(ω)=G(ω)H(ω)X(ω)
と表すことができる。ただしi(=1、2、…、M)は二次音源番号、j(=1、2、…、N)は制御点番号、ωは周波数であり、H(ω)は逆システムである。このとき、
Y(ω)=X(ω)
とするためには、INをN次の単位行列として、
G(ω)H(ω)=IN (1)
を満たす安定な逆システムH(ω)が必要である。
二次音源数をM、制御点数をNとすると、M<Nのとき、以下のように最小二乗解を計算することで、実際の解をG(ω)部分空間に射影した近似的な解を計算することができる。
H(ω)=(GT(ω)G(ω))-1GT(ω) (2)
GT(ω)はG(ω)の複素共役転置行列である。
二次音源iから制御点jまでの伝達関数をGji(ω)、入力信号をXj(ω)、観測信号をYj(ω)とすると、これらの関係は
Y(ω)=G(ω)H(ω)X(ω)
と表すことができる。ただしi(=1、2、…、M)は二次音源番号、j(=1、2、…、N)は制御点番号、ωは周波数であり、H(ω)は逆システムである。このとき、
Y(ω)=X(ω)
とするためには、INをN次の単位行列として、
G(ω)H(ω)=IN (1)
を満たす安定な逆システムH(ω)が必要である。
二次音源数をM、制御点数をNとすると、M<Nのとき、以下のように最小二乗解を計算することで、実際の解をG(ω)部分空間に射影した近似的な解を計算することができる。
H(ω)=(GT(ω)G(ω))-1GT(ω) (2)
GT(ω)はG(ω)の複素共役転置行列である。
また本発明では単純正則化法により、逆フィルタの不安定性の緩和を行った。単純正則化法による逆フィルタの計算式を式(3)に示す。
H^reg(ω)={GT(ω)G(ω)+β(ω)IN}-1GT(ω) (3)
ここで、H^reg(ω)は目的の逆システムを表す。また、βは正則化パラメータと呼ばれ、各計算周波数ごとに最適化する必要がある。
H^reg(ω)={GT(ω)G(ω)+β(ω)IN}-1GT(ω) (3)
ここで、H^reg(ω)は目的の逆システムを表す。また、βは正則化パラメータと呼ばれ、各計算周波数ごとに最適化する必要がある。
3. 二次音源配置の最適化
3.1 グラムシュミットの直交化法を応用した音源配置最適化アルゴリズム
グラムシュミットの直交化法を応用した音源配置選択アルゴリズムについて、総体法と逐次法について述べる。
3.1 グラムシュミットの直交化法を応用した音源配置最適化アルゴリズム
グラムシュミットの直交化法を応用した音源配置選択アルゴリズムについて、総体法と逐次法について述べる。
3.1.1 基本となる最適化アルゴリズム
まず、どちらの手法においても基本となる、Gram-Schmidtの直交化法を応用した音源配置選択アルゴリズムについて述べる。空間にNcヶ所の音源配置候補がある場合、この空間に何個の音源をどう配置して原音場を再現するかを考える。各音源配置候補の伝達インピーダンスベクトルziを
zi=[z1i…zMi]T (4)
とする。zji(i=1、2、…、Nc、j=1、2、…、M)はi番目のスピーカからj番目の制御点までの伝達インピーダンスであり、Mはセンサ数を、Tは共役転置を表す。また、Nc個の音源候補に対する伝達インピーダンスベクトルの集合を
τ={z1、…、zNc} (5)
とする。
まず、どちらの手法においても基本となる、Gram-Schmidtの直交化法を応用した音源配置選択アルゴリズムについて述べる。空間にNcヶ所の音源配置候補がある場合、この空間に何個の音源をどう配置して原音場を再現するかを考える。各音源配置候補の伝達インピーダンスベクトルziを
zi=[z1i…zMi]T (4)
とする。zji(i=1、2、…、Nc、j=1、2、…、M)はi番目のスピーカからj番目の制御点までの伝達インピーダンスであり、Mはセンサ数を、Tは共役転置を表す。また、Nc個の音源候補に対する伝達インピーダンスベクトルの集合を
τ={z1、…、zNc} (5)
とする。
まず、第1ステップとして、第1個目の音源を任意に選択する。このステップで、未使用の音源候補の部分集合をτ1とする。初期状態では、τ1=τである。この音源候補の集合τ1の中から、任意のzi∈τ1を選び、このベクトルに対応する音源を第1番目の音源とし、そのベクトルをz^1とする。また、第1番目の音源に対する正規直交ベクトルv1は
v1=z^1/|z^1| (6)
とする。
v1=z^1/|z^1| (6)
とする。
続いて、第2ステップ以下第nステップにおける音源の選択法について述べる。第n−1ステップまでの過程で、すでにn−1個の音源が選択されており、これらの音源に対応する正規直交ベクトルも求まっているものとする。この正規直交ベクトルの集合ν(n−1)を
ν(n−1)=v1、…、vn−1 (7)
とする。一つ前のステップである第n−1ステップにおいて未使用である伝達インピーダンスベクトルの集合をτn−1とすると、第nステップではτn−1の要素のうち、これまで決定された正規直交ベクトルと直交する成分が一番大きいベクトルを選択する。集合τn−1中の任意のベクトルziの正規直交ベクトルvj∈ν(n−1)上への射影pjは、
pj=(vj Hzi)vj (8)
と表される。これを用いて、複数の正規直交ベクトル{v1、…、vn−1}によって張られる平面へのziの射影pは、次式のように書ける。
p=Σ[j=1〜n-1]pj (9)
以上より、ベクトルziの成分のうち、この平面に直交する成分rは、
r=zi−p (10)
となる。第nステップでは、rのノルムが最大となるよう第n番目のベクトルz^nを決定する。すなわち、
z^n=arg max[zi∈τn−1]J(zi) (11)
ここで最大化される評価関数J(zi)は、次式で定義される。
J(zi)=|r| (12)
また、第n番目の正規直交ベクトルvnは次式により定義される。
vn=r/|r| (13)
最大化された評価関数は、次式のように書くものとする。
J^n=J(z^n) (14)
ν(n−1)=v1、…、vn−1 (7)
とする。一つ前のステップである第n−1ステップにおいて未使用である伝達インピーダンスベクトルの集合をτn−1とすると、第nステップではτn−1の要素のうち、これまで決定された正規直交ベクトルと直交する成分が一番大きいベクトルを選択する。集合τn−1中の任意のベクトルziの正規直交ベクトルvj∈ν(n−1)上への射影pjは、
pj=(vj Hzi)vj (8)
と表される。これを用いて、複数の正規直交ベクトル{v1、…、vn−1}によって張られる平面へのziの射影pは、次式のように書ける。
p=Σ[j=1〜n-1]pj (9)
以上より、ベクトルziの成分のうち、この平面に直交する成分rは、
r=zi−p (10)
となる。第nステップでは、rのノルムが最大となるよう第n番目のベクトルz^nを決定する。すなわち、
z^n=arg max[zi∈τn−1]J(zi) (11)
ここで最大化される評価関数J(zi)は、次式で定義される。
J(zi)=|r| (12)
また、第n番目の正規直交ベクトルvnは次式により定義される。
vn=r/|r| (13)
最大化された評価関数は、次式のように書くものとする。
J^n=J(z^n) (14)
この過程を逐次繰り返すことにより、各ステップで最も線形独立性の高いベクトルを持つ音源が一つずつ選択されていく。
各ステップにおいて採用された音源に対する評価関数J^nは伝達インピーダンスベクトルの線形独立性を示す尺度である。そして、採用される音源が増えるに従い、線形独立性は減少する、つまりJ^nは減少する。そこで、このJ^nにあらかじめ閾値J^thrを設けておき、ステップ数の増加に伴いJ^nが減少して、この閾値J^thrを下回った時点で音源選択を止め、この時点で選択されていた音源を制御に必要な音源とする。
各ステップにおいて採用された音源に対する評価関数J^nは伝達インピーダンスベクトルの線形独立性を示す尺度である。そして、採用される音源が増えるに従い、線形独立性は減少する、つまりJ^nは減少する。そこで、このJ^nにあらかじめ閾値J^thrを設けておき、ステップ数の増加に伴いJ^nが減少して、この閾値J^thrを下回った時点で音源選択を止め、この時点で選択されていた音源を制御に必要な音源とする。
3.1.2 BoSC音場再現システムへの適用
前節で述べたアルゴリズムは基本的なものであり、BoSC音場再現システムへ実装する場合はこれを拡張して考える必要がある。
総体法ではスピーカ選択のためのパラメータとして、式(15)によるJ^avgを用いる。
J^avg=(1/K)(a1|rω1|+…+aK|rωK|) (15)
ただし、ω1、…、ωKはそれぞれ制御対象とする離散周波数を、a1、…、aKはω1、…、ωKに対する重み係数を表す。本発明ではa1=…=aK=1とした。すなわち、式(11)は式(16)となる。
z^n=arg max[zi∈τn^1]Javg(zi) (16)
また、スピーカ選択を終了するパラメータとしては、式(17)で表されるJ^minを用いる。
J^min=min(a1|rω1|、…、aK|rωK|) (17)
すなわち、J^min<J^thrとなったときにスピーカ選択を終了する。
前節で述べたアルゴリズムは基本的なものであり、BoSC音場再現システムへ実装する場合はこれを拡張して考える必要がある。
総体法ではスピーカ選択のためのパラメータとして、式(15)によるJ^avgを用いる。
J^avg=(1/K)(a1|rω1|+…+aK|rωK|) (15)
ただし、ω1、…、ωKはそれぞれ制御対象とする離散周波数を、a1、…、aKはω1、…、ωKに対する重み係数を表す。本発明ではa1=…=aK=1とした。すなわち、式(11)は式(16)となる。
z^n=arg max[zi∈τn^1]Javg(zi) (16)
また、スピーカ選択を終了するパラメータとしては、式(17)で表されるJ^minを用いる。
J^min=min(a1|rω1|、…、aK|rωK|) (17)
すなわち、J^min<J^thrとなったときにスピーカ選択を終了する。
対して、今回比較対象とする逐次法では、制御対象とする離散周波数の低いものから前節で述べられたアルゴリズムを適用してスピーカを選択していく。まずは制御対象とする最低離散周波数であるω1において基本アルゴリズムを適用し、J^n<J^thrを満たすスピーカを二次音源として選択する。次のステップでは、前ステップで選択されたスピーカは既に選択されているものとして、ω2においてJ^n<J^thrを満たすスピーカを選択する。こうしてωKを対象とするステップまでスピーカ選択を実行する。
ここで問題となるのは、閾値J^thrの設定である。この閾値はシステムのダイナミックレンジを考慮して決定する必要がある。本発明では、両手法において同じ二次音源数で比較をするために閾値を調整した。
3.2 二次音源配置の選択
上で述べたアルゴリズムにより、音場再現に用いる二次音源数及びその配置を検討した。本発明では、真正面に位置するスピーカを初期選択音源とし、制御対象とする離散周波数はオクターブ刻みで80.5120Hzとした。本発明で用いる二次音源配置は総体法及び逐次法により選択された24ch配置、40ch配置の、計4配置とした。
上で述べたアルゴリズムにより、音場再現に用いる二次音源数及びその配置を検討した。本発明では、真正面に位置するスピーカを初期選択音源とし、制御対象とする離散周波数はオクターブ刻みで80.5120Hzとした。本発明で用いる二次音源配置は総体法及び逐次法により選択された24ch配置、40ch配置の、計4配置とした。
このようにして選択された二次音源配置を図3に示す。この図は音場再現システムの俯瞰図で、図の下方向が被験者にとって正面方向を表している。また、白丸が音場再現システムに設置されている全スピーカを、二重丸が初期選択スピーカを、黒丸が選択されたスピーカを表している。
図3右側上下の40chシステムでは、総体法と逐次法で大きく変わるところはなく、二次音源(スピーカ)は次のような配置となっている。すなわち、境界音場制御理論に基づき構成された、楕円形のドーム部に配置されたフルレンジスピーカと該ドーム部を支える4本の柱部に配置されたサブウーファースピーカから成る三次元音場再現システムにおいて、
a) 前記楕円形ドーム部の第1層の6個と、
b) 第1層の下部にある第2層の16個と、
c) 第2層の下部である第3層及び第4層の40個のうち聴者の正面の1個、背面の1個及び左右の16個と
の計40個のスピーカから成る。
a) 前記楕円形ドーム部の第1層の6個と、
b) 第1層の下部にある第2層の16個と、
c) 第2層の下部である第3層及び第4層の40個のうち聴者の正面の1個、背面の1個及び左右の16個と
の計40個のスピーカから成る。
また、図3左側上下の24chシステムでも、総体法と逐次法で大きく変わるところはなく、二次音源(スピーカ)は次のような配置となっている。すなわち、境界音場制御理論に基づき構成された、楕円形のドーム部に配置されたフルレンジスピーカと該ドーム部を支える4本の柱部に配置されたサブウーファースピーカから成る三次元音場再現システムにおいて、
a) 前記楕円形ドーム部の第1層の6個と、
b) 第1層の下部にある第2層の16個のうち左右各3個のそれぞれ、及び前後5個又は6個の計11個又は12個と、
c) 第2層の下部である第3層及び第4層のうち聴者の正面の1個、及び左右の6個又は5個の計7個又は6個と
の計24個のスピーカから成る。
a) 前記楕円形ドーム部の第1層の6個と、
b) 第1層の下部にある第2層の16個のうち左右各3個のそれぞれ、及び前後5個又は6個の計11個又は12個と、
c) 第2層の下部である第3層及び第4層のうち聴者の正面の1個、及び左右の6個又は5個の計7個又は6個と
の計24個のスピーカから成る。
4. 音場再現への影響の定量的評価
4.1 物理的影響の評価
前章で選択された各々の二次音源配置による音場再現への物理的な影響を評価するため、物理的指標として、式(18)による正確度を定義して、計算機シミュレーションでこれを求め、評価を行った。
Accuracy=log10(Σ|Pp(ω)|2/Σ|Pp(ω)−Pr(ω)|2) (18)
ただし、Pp(ω)は原音場の音圧、Pr(ω)は再現音場の音圧を表す。各二次音源配置におけるAccuracyの計算結果を図4に示す。計算対象周波数は20Hz刻みで20.20kHzまでとし、計算点はマイクロホンアレイの中心点、また実音源は正面方向にあるものとして求めた。
4.1 物理的影響の評価
前章で選択された各々の二次音源配置による音場再現への物理的な影響を評価するため、物理的指標として、式(18)による正確度を定義して、計算機シミュレーションでこれを求め、評価を行った。
Accuracy=log10(Σ|Pp(ω)|2/Σ|Pp(ω)−Pr(ω)|2) (18)
ただし、Pp(ω)は原音場の音圧、Pr(ω)は再現音場の音圧を表す。各二次音源配置におけるAccuracyの計算結果を図4に示す。計算対象周波数は20Hz刻みで20.20kHzまでとし、計算点はマイクロホンアレイの中心点、また実音源は正面方向にあるものとして求めた。
選択された二次音源数ごとに見ると低域での再現精度の差に顕著な差が見られるが、総体法と逐次法の間ではあまり差は見られない。
また、行列の縦ベクトルの線形独立性を示す指標として、「条件数」(Cond)が知られている[例えば、武者利光,岡本良夫,「逆問題とその解き方」, p. 145,オーム社,1992]。行列Aの条件数Cond(A)は式(19)で求められる。
Cond(A)=μmax(A)/μmin(A) (19)
ただし、μmax、μminはそれぞれ行列Aの最大特異値及び最小特異値を表す。条件数は、逆システム設計時に収録した伝達関数に含まれる雑音信号の拡大倍率と捉えることができ、条件数が大きいほどその逆システムでは雑音信号が拡大されて再生されるおそれがある。つまり、条件数が大きいほど不安定である、といえる。
Cond(A)=μmax(A)/μmin(A) (19)
ただし、μmax、μminはそれぞれ行列Aの最大特異値及び最小特異値を表す。条件数は、逆システム設計時に収録した伝達関数に含まれる雑音信号の拡大倍率と捉えることができ、条件数が大きいほどその逆システムでは雑音信号が拡大されて再生されるおそれがある。つまり、条件数が大きいほど不安定である、といえる。
選択された二次音源配置による伝達関数行列を比較するため、制御対象となる離散周波数ごとに各二次音源配置における伝達関数行列の条件数を求めた。結果を図5に示す。選択される二次音源数が多いほど、条件数が大きいことがわかる。どちらの手法においても、グラムシュミットの直交化によって線形独立性の高いスピーカから二次音源として選択されるため、選択された二次音源数が多くなるにつれて線形独立性は低下するためである。
また、総体法と逐次法とを比較すると、総体法のほうが条件数は小さく、またその差は選択された二次音源数が少ないほど顕著である。よって、選択された二次音源数が少ないほど、両手法による差は大きくなると考えられる。
さらに、図4と図5を見比べると、Accuracyと条件数が良い対応を示していることがわかる。
さらに、図4と図5を見比べると、Accuracyと条件数が良い対応を示していることがわかる。
4.2 心理的影響の評価
前章で選択された各々の二次音源配置による音場再現への心理的影響について評価を行うため、20歳代の聴力が正常な男女12名(男性8名、女性4名)を対象として、定位実験及び聴感実験を行った。
前章で選択された各々の二次音源配置による音場再現への心理的影響について評価を行うため、20歳代の聴力が正常な男女12名(男性8名、女性4名)を対象として、定位実験及び聴感実験を行った。
4.2.1 定位実験
音源として3秒間のピンクノイズを各二次音源配置に対応する逆システムで畳み込んだものを用い、水平面360度方向からの音を15度刻みで提示して、回答させた。提示順序は、ラテン方格を用いた実験計画法[秋田剛,「実験研究におけるサンプルサイズ−実験計画法と分散分析の適用−」,日本建築学会「建築空間における感覚・知覚心理シンポジウム(第7回)被験者は何人必要か?心理実験・調査研究におけるサンプリング」資料,別1-別4, 2009]に従い、順序効果を打ち消すよう無作為化を行い、また被験者1名につき本試験を2回ずつ行うことで、本試験を計24回行った。
音源として3秒間のピンクノイズを各二次音源配置に対応する逆システムで畳み込んだものを用い、水平面360度方向からの音を15度刻みで提示して、回答させた。提示順序は、ラテン方格を用いた実験計画法[秋田剛,「実験研究におけるサンプルサイズ−実験計画法と分散分析の適用−」,日本建築学会「建築空間における感覚・知覚心理シンポジウム(第7回)被験者は何人必要か?心理実験・調査研究におけるサンプリング」資料,別1-別4, 2009]に従い、順序効果を打ち消すよう無作為化を行い、また被験者1名につき本試験を2回ずつ行うことで、本試験を計24回行った。
被験者全員の結果をまとめたものを図6に示す。横軸は提示角度を、縦軸は被験者の回答角度を表し、図中の円の面積は回答した回数を表し、円が実線上にある場合は正解を、点線上にある場合は前後誤判定を表している。
左と右とを間違えるような大きな間違いは無いことが見てわかる。
左と右とを間違えるような大きな間違いは無いことが見てわかる。
定位実験の結果から、各二次音源配置ごとに正答率、前後誤判定率、回答誤差のRMS値を算出した。回答誤差のRMS値は、前後誤判定に大きく左右される(例えば、真正面と真後ろを誤認すると、その回答誤差は180度にもなる)。そのため、前後誤判定を取り除いた回答誤差のRMS値についても算出した。また、前後誤判定は個人差が大きいパラメータであることが知られているので、正答率についても前後誤判定を取り除いたものを算出した。
これらについて有意水準を0.05として統計解析を行った。統計解析の手順としては、はじめに等分散性を検定するためにバートレット(Bartlett)検定を行う。ここで各水準の母分散に有意差がなければ分散分析を、有意差があればフリードマン(Friedman)検定により、母平均の差の検定を行う。
バートレット検定の結果、等分散性が仮定できたので、パラメトリック検定により母平均の差の検定を行った。本発明では、定位能力には個人差があるとの見解のもと、反復測定による一要因分散分析を行った。
有意水準0.05においては、どの値においても、比較した5群の母平均に有意な差は見られなかった。このことから、比較した二次音源数(24ch、40ch)まで二次音源数を減らしても、定位精度には大きな影響を与えないということがわかった。
つぎに、正答率、前後誤判定を取り除いた正答率、前後誤判定率、回答誤差のRMS値、前後誤判定を取り除いた回答誤差のRMS値それぞれの標本平均を図7に示す。エラーバーは95%信頼区間を表す。標本平均で比較すると、正答率は、40ch再現では24ch再現に比べて改善され、62ch再現とほぼ同じ水準となっている。また選択手法で比較すると、逐次法のほうがわずかではあるが良い結果が得られた。ただし先に述べたように、有意差は見られなかった。
回答誤差のRMS値については62ch再現における値が大きく、また総体法で比較すると24ch再現よりも40ch再現のほうが大きい値をとっており、これらは推測される結果と異なる。しかし、前後誤判定を取り除いた差の二乗平均平方根(RMS of Difference(FBC))を見ると、その差はほとんど無くなっていることがわかる。このことから左のグラフ(RMS of Difference)に見られる差は前後誤判定率の差の寄与が大きいと推測できる。
そこで、二次音源配置と被験者を要因として前後誤判定率の二要因分散分析を行ったところ、二次音源配置間では有意差は見られなかったが、被験者間では有意差が見られた。このことから、前後誤判定率に見られるわずかな差は、二次音源配置に起因するものではなく、個人差によるものであると言える。
4.2.2 聴感実験
聴感実験は、ITU-RBS.1116[ITU-R Recommendation BS.1116, "Methods for the subjective assessment of small impairments in audio systems including multichannel sound systems", Geneva 1994](以下ITU勧告)を参考に、隠れ基準付き三刺激二重盲検法に基づいて行った。隠れ基準付き三刺激二重盲検法では、基準音(Ref.)に対する、評価音(A, B)の違いを評価する。ITU勧告では、評価音AとBのうち、基準音と同一であると感じた方に5.0、もう一方に劣化の度合いによって評点するとされているが、本実験では、基準音と同一と感じた方を回答した上で、他方を評価させた。評点の基準としてはITU勧告の評価基準を用いた。実験に用いる音源としては、音楽、音声、環境音の3種類、各15秒程度を抽出して用いた。これらに各二次音源配置に対応する逆システムを畳み込み、再生用の信号を作成した。
聴感実験は、ITU-RBS.1116[ITU-R Recommendation BS.1116, "Methods for the subjective assessment of small impairments in audio systems including multichannel sound systems", Geneva 1994](以下ITU勧告)を参考に、隠れ基準付き三刺激二重盲検法に基づいて行った。隠れ基準付き三刺激二重盲検法では、基準音(Ref.)に対する、評価音(A, B)の違いを評価する。ITU勧告では、評価音AとBのうち、基準音と同一であると感じた方に5.0、もう一方に劣化の度合いによって評点するとされているが、本実験では、基準音と同一と感じた方を回答した上で、他方を評価させた。評点の基準としてはITU勧告の評価基準を用いた。実験に用いる音源としては、音楽、音声、環境音の3種類、各15秒程度を抽出して用いた。これらに各二次音源配置に対応する逆システムを畳み込み、再生用の信号を作成した。
基準音と同じ評価音を正確に判断できた確率を正答率として、各音源ごとに正答率の標本平均を求めた。結果を図8に示す。二次音源数で比較すると、40ch再現の場合に比べて24ch再現のほうが音の違いを認知しやすく、正答率が高くなっている。特に、音楽を用いた場合の24ch再現における正答率は95%以上と著しく高い。
また、音源で比較しても違いが見られる。音楽、音声では24ch再現の場合と40ch再現の場合とで顕著な差が見られるが、環境音ではほとんど差が見られない。
次に、評価音の点数から隠れ基準音の点数を引いた値を評価音の差度(Xi)として、音源ごとに被験者間で式(20)を用いて点数の正規化を行った。X~si、Ssiは被験者の1試験(Session)ごとの平均値と標準偏差、X~s、Ssは音源ごとの全被験者の平均値と標準偏差である。
Zi=((Xi−X~si)/Ssi)Ss+X~s (20)
Zi=((Xi−X~si)/Ssi)Ss+X~s (20)
24回の実験結果を図9に示す。図の縦軸は各音源、二次音源配置に対して得られた差度を正規化した点数の平均値である。エラーバーは95%信頼区間を示す。いずれの項目においても、40ch再現では標本平均が−1.0以上であり、基準音である62ch再現との違いはあまり気にならない程度であると言える。
この結果について統計的な分析を行うために、以下の手順で検定を行った。まず、分散分析を行うには、要因内の水準間で等分散性を示す必要がある。本分析においてはバートレット検定によって等分散性を確認した。等分散性が仮定できる場合には、分散分析により母平均を、有意差が見られた場合にはフリードマン検定によって母集団の中央値を比較する。バートレット検定の結果、音源を要因とした水準間の分散には有意水準α=0.05において有意差が見られた。そこで、音源ごとに二次音源配置選択アルゴリズムと二次音源数を要因とした群間でそれぞれ等分散性を検定したところ、等分散性が仮定できた。そこで、アルゴリズムと二次音源数を要因とし、二要因とも対応あり・繰り返しありの二要因分散分析を行った。その結果、二次音源数の主効果には有意差があり、選択アルゴリズムの主効果と二次音源数とアルゴリズムの交互作用には有意差がないことが分かった。また、40ch再現の標本平均を比較すると有意差はないものの、音楽と音声では逐次法よりも総体法の方が良く、環境音では逐次法の方が僅かに良かった。
評価音について評点をつける際に、音のどの部分に違いを感じたか、またはどのような違いを感じたかを併せて自由に記述させた。その結果、「高音域の音が足りない、こもって聞こえる」などの高音域への影響を示唆する回答、「上方向からの音が減る」、「鳥の鳴き声が異なる、音の聞こえてくる方向が変わる」などの上下方向の定位への影響を及ぼす回答が得られた。しかし、二次音源配置の違いを明確にするような回答は得られなかった。また、音楽については、「音声や環境音よりも違いを見つけやすかった」との感想を多く得た。これは、コンテンツ間での正答率の差や差度の値がより小さいこととも一致し、音源の種類によって評価の感度が変わることを示している。
5. 結論
本発明では、BoSC音場再現システムにおける二次音源配置について、グラムシュミットの直交化を応用した制御音源配置の最適化手法とその変法を用いて4通りの音源配置を選択し、二次音源配置の違いが再現音場にどのような物理的影響及び心理的影響を与えるかを計算機シミュレーション及び主観評価実験によって調べた。
本発明では、BoSC音場再現システムにおける二次音源配置について、グラムシュミットの直交化を応用した制御音源配置の最適化手法とその変法を用いて4通りの音源配置を選択し、二次音源配置の違いが再現音場にどのような物理的影響及び心理的影響を与えるかを計算機シミュレーション及び主観評価実験によって調べた。
物理的影響については、式(18)により定義した正確度、及び式(19)により定義される条件数を計算機シミュレーションにより求め、考察した。その結果、正確度及び条件数では、特に低域で、二次音源数の違いが顕著にあらわれた。しかし、今回用いた2つの選択手法、総体法と逐次法には顕著な差は見られなかった。また、周波数ごとに見たときの正確度と条件数はよく対応していることがわかった。
心理的影響については、定位実験と聴感実験を行った。
定位実験では正答率・前後誤判定率・回答誤差のRMS値によって音場再現への心理的影響を評価した。分散分析の結果、これらの値には有意差は見られなかったが、標本平均で正答率を62ch再現の場合と比較すると、40ch再現ではその差はほとんどないが、24ch再現では低かった。また、回答誤差のRMS値では62ch再現における値が大きかったが、これは前後誤判定率の差によるものであると考えられる。総体法と逐次法による音源配置には差はほとんど見られなかった。
定位実験では正答率・前後誤判定率・回答誤差のRMS値によって音場再現への心理的影響を評価した。分散分析の結果、これらの値には有意差は見られなかったが、標本平均で正答率を62ch再現の場合と比較すると、40ch再現ではその差はほとんどないが、24ch再現では低かった。また、回答誤差のRMS値では62ch再現における値が大きかったが、これは前後誤判定率の差によるものであると考えられる。総体法と逐次法による音源配置には差はほとんど見られなかった。
聴感実験では、隠れ基準つき三刺激二重盲検法により音場再現への心理的影響を評価した。実験の結果より、聴感上には特に24ch再現の場合において劣化が感じられることが分かった。総体法と逐次法による音源配置にはほとんど差は見られなかった。さらに、自由記述回答の結果について考察したところ、高音域の再現性や再現音場の上下方向の定位への影響が示唆された。
結論としては、正確度と条件数から調べた物理的影響についてはその差は大きかった。しかし、主観評価実験によって調べた心理的影響の観点から述べると、定位精度への影響よりは聴感上の影響の方が大きく、40chまでは二次音源を減らしても音場再現への影響は少ないと言えることがわかった。ただし、高音域の再現性や再現音場の上下方向の定位・拡がりへの影響が示唆された。また、二次音源配置の最適化手法として比較した総体法と逐次法については、どちらの手法がより適しているかということについては決定できなかった。
Claims (9)
- 境界音場制御理論に基づき構成された、聴者の頭部を囲う複数のスピーカから成る三次元音場再現システムにおいて、
a) 所定の基準スピーカと各評価対象スピーカの間のグラムシュミット直交化評価値を計算し、
b) グラムシュミット直交化評価値の高い順に評価対象スピーカを選択し、
c) 基準スピーカと選択された評価対象スピーカから成るサブシステムと未選択評価対象スピーカの間のグラムシュミット直交化評価値を計算し、
d) 前記b)及びc)のステップを繰り返し、該グラムシュミット直交化評価値が所定の閾値を下回った時点で選択を停止する
ことを特徴とする三次元音場再現システムの音源配置最適化法。 - グラムシュミット直交化評価値を計算する際、複数の制御対象周波数を全て評価関数に含めておくことを特徴とする請求項1に記載の三次元音場再現システムの音源配置最適化法。
- グラムシュミット直交化評価値を計算する際、評価関数に含まれる複数の制御対象周波のそれぞれに重み付け係数を付与することを特徴とする請求項2に記載の三次元音場再現システムの音源配置最適化法。
- グラムシュミット直交化評価値を計算する際、複数の制御対象周波数を低い周波数から順に評価関数に含め、各周波数において前記選択を行うことを特徴とする請求項1に記載の三次元音場再現システムの音源配置最適化法。
- 境界音場制御理論に基づき構成された、聴者の頭部を囲う複数のスピーカから成る三次元音場再現システムにおいて、
a) 所定の基準スピーカと各評価対象スピーカの間のグラムシュミット直交化評価値を計算し、
b) グラムシュミット直交化評価値の高い順に評価対象スピーカを選択し、
c) 基準スピーカと選択された評価対象スピーカから成るサブシステムと未選択評価対象スピーカの間のグラムシュミット直交化評価値を計算し、
d) 前記b)及びc)のステップを繰り返し、該グラムシュミット直交化評価値が所定の閾値を下回った時点で選択を停止する
ことにより選択されたスピーカから成ることを特徴とする三次元音場再現システム。 - 境界音場制御理論に基づき構成された、楕円形のドーム部に配置されたフルレンジスピーカと該ドーム部を支える4本の柱部に配置されたサブウーファースピーカから成る三次元音場再現システムにおいて、
a) 所定の基準スピーカと各評価対象スピーカの間のグラムシュミット直交化評価値を計算し、
b) グラムシュミット直交化評価値の高い順に評価対象スピーカを選択し、
c) 基準スピーカと選択された評価対象スピーカから成るサブシステムと未選択評価対象スピーカの間のグラムシュミット直交化評価値を計算し、
d) 前記b)及びc)のステップを繰り返し、該グラムシュミット直交化評価値が所定の閾値を下回った時点で選択を停止する
ことにより選択された40個のスピーカから成ることを特徴とする三次元音場再現システム。 - 境界音場制御理論に基づき構成された、楕円形のドーム部に配置されたフルレンジスピーカと該ドーム部を支える4本の柱部に配置されたサブウーファースピーカから成る三次元音場再現システムにおいて、
a) 所定の基準スピーカと各評価対象スピーカの間のグラムシュミット直交化評価値を計算し、
b) グラムシュミット直交化評価値の高い順に評価対象スピーカを選択し、
c) 基準スピーカと選択された評価対象スピーカから成るサブシステムと未選択評価対象スピーカの間のグラムシュミット直交化評価値を計算し、
d) 前記b)及びc)のステップを繰り返し、該グラムシュミット直交化評価値が所定の閾値を下回った時点で選択を停止する
ことにより選択された24個のスピーカから成ることを特徴とする三次元音場再現システム。 - 境界音場制御理論に基づき構成された、楕円形のドーム部に配置されたフルレンジスピーカと該ドーム部を支える4本の柱部に配置されたサブウーファースピーカから成る三次元音場再現システムにおいて、
a) 前記楕円形ドーム部の第1層の6個と、
b) 第1層の下部にある第2層の16個と、
c) 第2層の下部である第3層及び第4層の40個のうち聴者の正面の1個、背面の1個及び左右の16個と
の計40個のスピーカから成ることを特徴とする三次元音場再現システム。 - 境界音場制御理論に基づき構成された、楕円形のドーム部に配置されたフルレンジスピーカと該ドーム部を支える4本の柱部に配置されたサブウーファースピーカから成る三次元音場再現システムにおいて、
a) 前記楕円形ドーム部の第1層の6個と、
b) 第1層の下部にある第2層の16個のうち左右各3個のそれぞれ、及び前後5個又は6個の計11個又は12個と、
c) 第2層の下部である第3層及び第4層のうち聴者の正面の1個、及び左右の6個又は5個の計7個又は6個と
の計24個のスピーカから成ることを特徴とする三次元音場再現システム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010043335A JP2011182135A (ja) | 2010-02-26 | 2010-02-26 | 三次元音場再現システム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010043335A JP2011182135A (ja) | 2010-02-26 | 2010-02-26 | 三次元音場再現システム |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2011182135A true JP2011182135A (ja) | 2011-09-15 |
Family
ID=44693188
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010043335A Pending JP2011182135A (ja) | 2010-02-26 | 2010-02-26 | 三次元音場再現システム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2011182135A (ja) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012085035A (ja) * | 2010-10-08 | 2012-04-26 | National Institute Of Information & Communication Technology | 音場共有システムおよび最適化方法 |
JP2012109643A (ja) * | 2010-11-15 | 2012-06-07 | National Institute Of Information & Communication Technology | 音再現システム、音再現装置および音再現方法 |
GB2495623A (en) * | 2011-10-10 | 2013-04-17 | Korea Electronics Telecomm | Preventing excessive output of energy when reproducing a three dimensional (3D) sound field |
KR20130038790A (ko) * | 2011-10-10 | 2013-04-18 | 한국전자통신연구원 | 3차원 음장 표현 방법 및 장치 |
JP2013172236A (ja) * | 2012-02-20 | 2013-09-02 | Nippon Telegr & Teleph Corp <Ntt> | 音場収音再生装置、方法及びプログラム |
JP2014007543A (ja) * | 2012-06-25 | 2014-01-16 | Nippon Telegr & Teleph Corp <Ntt> | 音場収音再生装置、方法及びプログラム |
JP2014150416A (ja) * | 2013-02-01 | 2014-08-21 | Nippon Telegr & Teleph Corp <Ntt> | 音場収音再生装置、方法及びプログラム |
JP2014150415A (ja) * | 2013-02-01 | 2014-08-21 | Nippon Telegr & Teleph Corp <Ntt> | 音場収音再生装置、方法及びプログラム |
JP2016082443A (ja) * | 2014-10-17 | 2016-05-16 | 学校法人 中央大学 | スピーカ配置選択装置、スピーカ配置選択方法及び音場制御システム |
JP7445975B2 (ja) | 2020-10-30 | 2024-03-08 | 株式会社カスクアコースティクス | マイクロホンアレイ装置、音場制御システム、音場収録システム、及び音場再生システム |
-
2010
- 2010-02-26 JP JP2010043335A patent/JP2011182135A/ja active Pending
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012085035A (ja) * | 2010-10-08 | 2012-04-26 | National Institute Of Information & Communication Technology | 音場共有システムおよび最適化方法 |
JP2012109643A (ja) * | 2010-11-15 | 2012-06-07 | National Institute Of Information & Communication Technology | 音再現システム、音再現装置および音再現方法 |
GB2495623A (en) * | 2011-10-10 | 2013-04-17 | Korea Electronics Telecomm | Preventing excessive output of energy when reproducing a three dimensional (3D) sound field |
KR20130038790A (ko) * | 2011-10-10 | 2013-04-18 | 한국전자통신연구원 | 3차원 음장 표현 방법 및 장치 |
GB2495623B (en) * | 2011-10-10 | 2015-10-21 | Korea Electronics Telecomm | Method and apparatus for representing three-dimensional (3D) sound field |
JP2013172236A (ja) * | 2012-02-20 | 2013-09-02 | Nippon Telegr & Teleph Corp <Ntt> | 音場収音再生装置、方法及びプログラム |
JP2014007543A (ja) * | 2012-06-25 | 2014-01-16 | Nippon Telegr & Teleph Corp <Ntt> | 音場収音再生装置、方法及びプログラム |
JP2014150416A (ja) * | 2013-02-01 | 2014-08-21 | Nippon Telegr & Teleph Corp <Ntt> | 音場収音再生装置、方法及びプログラム |
JP2014150415A (ja) * | 2013-02-01 | 2014-08-21 | Nippon Telegr & Teleph Corp <Ntt> | 音場収音再生装置、方法及びプログラム |
JP2016082443A (ja) * | 2014-10-17 | 2016-05-16 | 学校法人 中央大学 | スピーカ配置選択装置、スピーカ配置選択方法及び音場制御システム |
JP7445975B2 (ja) | 2020-10-30 | 2024-03-08 | 株式会社カスクアコースティクス | マイクロホンアレイ装置、音場制御システム、音場収録システム、及び音場再生システム |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2011182135A (ja) | 三次元音場再現システム | |
Spors et al. | Spatial sound with loudspeakers and its perception: A review of the current state | |
CN109417678A (zh) | 声场形成装置和方法以及程序 | |
Lee et al. | Effect of vertical microphone layer spacing for a 3D microphone array | |
Kobayashi et al. | The effects of spatialized sounds on the sense of presence in auditory virtual environments: a psychological and physiological study | |
CA2953808A1 (en) | Error model-based multi-zone sound reproduction method and device | |
Lindau | Binaural resynthesis of acoustical environments: technology and perceptual evaluation | |
Kamekawa et al. | Evaluation of recording techniques for three-dimensional audio recordings: Comparison of listening impressions based on difference between listening positions and three recording techniques | |
Kaplanis et al. | Perception and preference of reverberation in small listening rooms for multi-loudspeaker reproduction | |
Riedel et al. | Surrounding line sources optimally reproduce diffuse envelopment at off-center listening positions | |
US11510013B2 (en) | Partial HRTF compensation or prediction for in-ear microphone arrays | |
Tregonning et al. | The vertical precedence effect: Utilizing delay panning for height channel mixing in 3D audio | |
Rämö et al. | Validating a perceptual distraction model using a personal two-zone sound system | |
Zaunschirm et al. | Perceptual evaluation of variable-orientation binaural room impulse response rendering | |
Lindau | Binaural resynthesis of acoustical environments | |
Michaud et al. | Objective characterization of perceptual dimensions underlying the sound reproduction of 37 single loudspeakers in a room | |
Salmon et al. | The influence of the sound source on perceived differences between binaurally rendered sound spaces | |
Mittag et al. | Development and evaluation of methods for the synthesis of binaural room impulse responses based on spatially sparse measurements in real rooms | |
Kiyan et al. | Towards predicting immersion in surround sound music reproduction from sound field features | |
Li et al. | Evaluation of headphone phase equalization on sound reproduction | |
Hughes et al. | Dual frequency band amplitude panning for multichannel audio systems | |
Paine et al. | Perceptual investigation into envelopement, spatial clarity, and engulfment in reproduced multi-channel audio | |
Zhang et al. | Intelligent cloud computing platform for three‐dimensional sound reproduction | |
Werner et al. | Context-dependent quality parameters and perception of auditory illusions | |
Clapp et al. | Investigations of room acoustics with a spherical microphone array |