JP6955186B2 - 音響信号処理装置、音響信号処理方法および音響信号処理プログラム - Google Patents

音響信号処理装置、音響信号処理方法および音響信号処理プログラム Download PDF

Info

Publication number
JP6955186B2
JP6955186B2 JP2020503603A JP2020503603A JP6955186B2 JP 6955186 B2 JP6955186 B2 JP 6955186B2 JP 2020503603 A JP2020503603 A JP 2020503603A JP 2020503603 A JP2020503603 A JP 2020503603A JP 6955186 B2 JP6955186 B2 JP 6955186B2
Authority
JP
Japan
Prior art keywords
acoustic signal
coordinates
speaker
signal processing
focal
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.)
Active
Application number
JP2020503603A
Other languages
English (en)
Other versions
JPWO2019168083A1 (ja
Inventor
公孝 堤
公孝 堤
賢一 野口
賢一 野口
高田 英明
英明 高田
羽田 陽一
陽一 羽田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nippon Telegraph and Telephone Corp
Original Assignee
Nippon Telegraph and Telephone Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nippon Telegraph and Telephone Corp filed Critical Nippon Telegraph and Telephone Corp
Publication of JPWO2019168083A1 publication Critical patent/JPWO2019168083A1/ja
Application granted granted Critical
Publication of JP6955186B2 publication Critical patent/JP6955186B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04RLOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R3/00Circuits for transducers, loudspeakers or microphones
    • H04R3/12Circuits for transducers, loudspeakers or microphones for distributing signals to two or more loudspeakers
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04SSTEREOPHONIC SYSTEMS 
    • H04S7/00Indicating arrangements; Control arrangements, e.g. balance control
    • H04S7/30Control circuits for electronic adaptation of the sound field
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04RLOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R1/00Details of transducers, loudspeakers or microphones
    • H04R1/20Arrangements for obtaining desired frequency or directional characteristics
    • H04R1/32Arrangements for obtaining desired frequency or directional characteristics for obtaining desired directional characteristic only
    • H04R1/40Arrangements for obtaining desired frequency or directional characteristics for obtaining desired directional characteristic only by combining a number of identical transducers
    • H04R1/403Arrangements for obtaining desired frequency or directional characteristics for obtaining desired directional characteristic only by combining a number of identical transducers loud-speakers
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04RLOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R2201/00Details of transducers, loudspeakers or microphones covered by H04R1/00 but not provided for in any of its subgroups
    • H04R2201/40Details of arrangements for obtaining desired directional characteristic by combining a number of identical transducers covered by H04R1/40 but not provided for in any of its subgroups
    • H04R2201/403Linear arrays of transducers
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04RLOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R2499/00Aspects covered by H04R or H04S not otherwise provided for in their subgroups
    • H04R2499/10General applications
    • H04R2499/13Acoustic transducers and sound field adaptation in vehicles
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04RLOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R27/00Public address systems
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04SSTEREOPHONIC SYSTEMS 
    • H04S2400/00Details of stereophonic systems covered by H04S but not provided for in its groups
    • H04S2400/07Generation or adaptation of the Low Frequency Effect [LFE] channel, e.g. distribution or signal processing
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04SSTEREOPHONIC SYSTEMS 
    • H04S2400/00Details of stereophonic systems covered by H04S but not provided for in its groups
    • H04S2400/11Positioning of individual sound objects, e.g. moving airplane, within a sound field
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04SSTEREOPHONIC SYSTEMS 
    • H04S2420/00Techniques used stereophonic systems covered by H04S but not provided for in its groups
    • H04S2420/07Synergistic effects of band splitting and sub-band processing

Landscapes

  • Health & Medical Sciences (AREA)
  • Otolaryngology (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Signal Processing (AREA)
  • General Health & Medical Sciences (AREA)
  • Circuit For Audible Band Transducer (AREA)
  • Obtaining Desirable Characteristics In Audible-Bandwidth Transducers (AREA)
  • Stereophonic System (AREA)

Description

本発明は、入力音響信号を、仮想音源を実現するための複数のスピーカを並べて構成したスピーカアレイの各スピーカへの出力音響信号に変換する音響信号処理装置、音響信号処理方法および音響信号処理プログラムに関する。
パブリックビューイングやコンサートでは、上映会場に設置した複数のスピーカから音声や音楽などを再生する。近年、仮想的な音源を上映空間に作り出すことにより、これまで以上に臨場感のある音響再生を実現する取り組みが行われている。特に直線状に多数のスピーカを並べてできるスピーカアレイを用いて、スピーカより前面、客席近くまで飛び出る仮想音源を生成することで高い臨場感を実現するといったことが行なわれている。
また、一般に、楽器や人間の声は方向によって放射されるパワーが異なるため、上映空間に仮想的な音源を生成する際に方向による音響信号のパワーの違い(指向性)を再現することで、さらに臨場感の高い音響コンテンツを実現することが期待されている。
上映空間に仮想的な音源を作り出す音響再生技術に対し、波面合成と呼ばれる方法がある(特許文献1)。特許文献1に基づく方法は、音響信号を収録する地点の音響信号を複数地点に設置したマイクロフォンで収音した上で、上下左右方向の音響信号の到来方向を分析し、上映空間中に設置した複数のスピーカを用いて収録会場の音響信号を物理的に再現する。
想定する仮想音源に吸込み型音源(acoustic sink)を仮定し、第1種レイリー積分から導出される駆動信号をスピーカアレイに与えることにより、スピーカより前面に仮想音像を作り出すことができる技術がある(非特許文献1)。また、直線状スピーカアレイを用いて上映空間に生成する仮想的な音源にダイポールなどの原始的な指向性を実現できる技術がある(非特許文献2)。
スピーカから放射される音の指向性を制御する方法として、多重極音源がある(非特許文献3)。多重極音源は、音の指向性をダイポール、クアドラポールといった原始的な指向性の組み合わせで表現する手法であり、原始的な指向性それぞれは互いに近接した極性の異なる無指向性の点音源(モノポール音源)の組み合わせで実現される。非特許文献3は、指向性の向きを回転させるには、強度の異なる原始的な指向性を重ねあわせることを開示する。
特開2011−244306号公報
Sascha Spors, Hagen Wierstorf, Matthias Gainer, and Jens Ahrens, "Physical and Perceptual Properties of Focused Sources in Wave Field Synthesis," in 127th Audio Engineering Society Convention paper 7914, 2009, October. J. Ahrens, and S. Spors, "Implementation of Directional Sources in Wave Field Synthesis," Proceeding of IEEE Workshop on Applications of Signal Processing to Audio and Acoustics, pp. 66-69, 2007. 羽田陽一,古家賢一,島内末廣,"球調和関数展開に基づく多重極音源を用いた指向性合成",日本音響学会誌 69巻 11号 pp577-588 2013.
しかしながらいずれの文献も、多重極子を重ね合わせて任意の指向特性を実現する技術については触れられていないので、いずれの文献からも、多重極子を用いて、楽器などの音源が持つ指向性をモデル化することは難しい。
従って本発明の目的は、多重極子を重ね合わせて任意の指向特性を実現する音響信号処理装置、音響信号処理方法および音響信号処理プログラムを提供することである。
上記課題を解決するために、本発明の第1の特徴は、入力音響信号を、仮想音源を実現するための複数のスピーカを並べて構成したスピーカアレイの各スピーカへの出力音響信号に変換する音響信号処理装置に関する。本発明の第1の特徴に係る音響信号処理装置は、複数の初期焦点座標と、仮想音源の座標および指向性の方向を取得して、複数の初期焦点座標のうち極性の異なる一対の初期焦点座標について、仮想音源の座標に基づいて、初期焦点座標に指向性の方向から特定される回転行列をかけて、焦点座標を決定する焦点位置決定部と、円調和係数から、焦点座標を含む多重極子に与える重みを算出する円調和係数変換部と、スピーカアレイの各スピーカについて、焦点座標と、焦点座標の極性と、多重極子に与える重みから、スピーカに与える重み付き駆動関数を演算するフィルタ係数演算部と、スピーカアレイの各スピーカについて、入力音響信号に、スピーカに対応する重み付き駆動関数を畳み込んで、スピーカへの出力音響信号を出力する畳み込み演算部を備える。
ここで円調和係数変換部は、式(1)により、多重極子に与える重みを算出しても良い。
Figure 0006955186
フィルタ係数演算部は、焦点座標のそれぞれを用いて駆動関数を算出し、多重極子のそれぞれについて、多重極子を構成する各焦点座標の極性と駆動関数から算出した合成駆動関数と、多重極子に与える重みから、スピーカに与える重み付き駆動関数を演算しても良い。
フィルタ係数演算部は、多重極子に含まれる各焦点座標について、焦点座標の極性と駆動関数をかけた関数を加算して、多重極子の合成駆動関数を算出しても良い。
フィルタ係数演算部は、多重極子のそれぞれについて算出した合成駆動関数に、多重極子に与える重みをかけて加算して、重み付き駆動関数を算出しても良い。
本発明の第2の特徴は、入力音響信号を、仮想音源を実現するための複数のスピーカを並べて構成したスピーカアレイの各スピーカへの出力音響信号に変換する音響信号処理方法に関する。本発明の第2の特徴に係る音響信号処理方法は、複数の初期焦点座標と、仮想音源の座標および指向性の方向を取得するステップと、複数の初期焦点座標のうち極性の異なる一対の初期焦点座標について、仮想音源の座標に基づいて、初期焦点座標に指向性の方向から特定される回転行列をかけて、焦点座標を決定するステップと、円調和係数から、焦点座標を含む多重極子に与える重みを算出するステップと、スピーカアレイの各スピーカについて、焦点座標と、焦点座標の極性と、多重極子に与える重みから、スピーカに与える重み付き駆動関数を演算するステップと、スピーカアレイの各スピーカについて、入力音響信号に、スピーカに対応する重み付き駆動関数を畳み込んで、スピーカへの出力音響信号を出力するステップを備える。
本発明の第3の特徴は、コンピュータに、第1の特徴に記載の音響信号処理装置として機能させるための音響信号処理プログラムに関する。
本発明によれば、多重極子を重ね合わせて任意の指向特性を実現する音響信号処理装置、音響信号処理方法および音響信号処理プログラムを提供することができる。
本発明の実施の形態に係る音響信号処理装置のブロック図である。 本発明の実施の形態において多重極子を重ね合わせて実現される指向特性を説明する図である。 本発明の実施の形態に係る音響信号処理装置の焦点位置決定処理を説明するフローチャートである。 本発明の実施の形態に係る音響信号処理装置の焦点位置決定処理において、初期焦点座標を説明する図である。 本発明の実施の形態に係る音響信号処理装置の焦点位置決定処理においてもちいられる回転行列の一例を説明する図である。 本発明の実施の形態に係る音響信号処理装置の焦点位置決定処理において、指向性が考慮された焦点座標を説明する図である。 本発明の実施の形態に係る音響信号処理装置の円調和係数変換処理を説明するフローチャートである。 本発明の実施の形態に係る音響信号処理装置のフィルタ係数演算処理を説明するフローチャートである。 本発明の実施の形態に係る音響信号処理装置のフィルタ係数演算処理で算出される各関数の一例を説明する図である。 本発明の実施の形態に係る音響信号処理装置の畳み込み演算処理を説明するフローチャートである。
次に、図面を参照して、本発明の実施の形態を説明する。以下の図面の記載において、同一または類似の部分には同一または類似の符号を付している。
(音響信号処理装置)
図1を参照して、本発明の実施の形態に係る音響信号処理装置1を説明する。音響信号処理装置1は、処理装置(図示せず)、メモリ10などを備える一般的なコンピュータである。一般的なコンピュータが音響信号処理プログラムを実行することにより図1に示す機能を実現する。
本発明の実施の形態に係る音響信号処理装置1は、図2に示すような、複数のスピーカを直線状に並べた直線状スピーカアレイを用いて、多重極子に重みをかけることにより、スピーカよりも前面に飛び出し、かつ、指向性を有する仮想音源を実現する。本発明の実施の形態において、スピーカアレイを構成する各スピーカは、直線状に並ぶ場合を説明するが、これに限られない。スピーカアレイは、複数のスピーカで構成されればよく、複数のスピーカが直線状に並ばなくても良い。
本発明の実施の形態においては、仮想音源を実現するために、互いに近接する位置に極性の異なる2以上の焦点音源を生成することで、多重極音源を実現する。焦点音源は、極性の異なる無指向性の点音源(モノポール音源)の組み合わせである。本発明の実施の形態において焦点音源は、2つの多重極子を備え、一つの多重極子は、1つのモノポール音源で構成され、もう一方の多重極子は、極性の異なる2つのモノポール音源で構成される場合を説明するが、これに限られない。
本発明の実施の形態では、図2(a)に示す多重極子M1と多重極子M2を重ねて、図2(b)に示す指向特性を実現する。多重極子M1は、プラス極性を持つ焦点P1を有し、多重極子M2は、マイナス極性を持つ焦点P2と、プラス極性を持つ焦点P3を有する。本発明の実施の形態は、多重極子M1および多重極子M2にそれぞれ重みをかけて重ね合わせることにより、図2(b)に示す多重極音源の指向特性を実現する。図2(b)に示すように、様々な指向特性を持つ多重極子を重ね合わせることで、所望の範囲で所望の指向特性を実現することが可能になる。
音響信号処理装置1は、このような仮想音源を実現するために、入力音響信号Iを、直線状スピーカアレイの各スピーカへの出力音響信号Oに変換する。
図1に示すように音響信号処理装置1は、メモリ10、焦点位置決定部12、円調和係数変換部13、フィルタ係数演算部14、畳み込み演算部15、入出力インタフェース(図示せず)等を備える。入出力インタフェースは、入力音響信号を音響信号処理装置1に入力し、各スピーカへの出力音響信号を出力するためのインタフェースである。入出力インタフェースは、音響信号処理装置1が実現する仮想音源の座標および指向性の方向の各情報、さらに円調和係数を、音響信号処理装置1に入力する。
メモリ10は、焦点データ11を記憶する。焦点データ11は、仮想音源を実現するための複数の焦点の座標と、各焦点の極性を対応づける。本発明の実施の形態において、焦点データ11に記憶される焦点を、初期焦点と称し、初期焦点の座標を初期焦点座標と称する。
焦点位置決定部12は、仮想音源の位置、指向性の方向の情報、および対象周波数の各情報を受け取り、指向性を考慮した、必要な数の焦点に関する座標を出力する。焦点位置決定部12は、複数の初期焦点座標と、仮想音源の座標および指向性の方向を取得して、複数の初期焦点座標のうち極性の異なる一対の初期焦点座標について、仮想音源の座標に基づいて、初期焦点座標に指向性の方向から特定される回転行列をかけて、焦点座標を決定する。焦点位置決定部12は、仮想音源の座標に対する初期焦点座標の相対座標に回転行列をかけ、回転行列をかけて得られた座標に、仮想音源の座標を加算して、指向性を考慮した焦点座標を決定する。なお、仮想音源は、焦点座標の中心となる。
焦点位置決定部12は、複数の初期焦点座標のうち、対をなさない初期焦点座標について、何ら変換することなく、初期焦点座標を焦点座標として決定する。図2に示す例において、プラス極性の焦点を有する多重極子M1について、焦点位置決定部12は、プラス極性の初期焦点座標を焦点座標として出力する。プラス極性の焦点とマイナス極性の焦点を有する多重極子M2について、焦点位置決定部12は、これらの初期焦点座標に回転をかけた座標を、焦点座標として出力する。
焦点位置決定部12は、メモリ10から、極性の異なる一対以上の初期焦点座標を取得するとともに、音響信号処理装置1が実現する特性として、外部入力等により、仮想音源の座標および指向性の方向を取得する。焦点位置決定部12は、取得した指向性の方向から、初期焦点座標に対してかける回転方向θを特定する。
焦点位置決定部12は、一対の初期焦点座標を、
Figure 0006955186
とした場合、X軸方向に対してθ方向を指定すると、この方向から特定できる回転行列Gは、式(1)で求まるため、回転後のモノポールの座標は、式(2)で決定できる。
Figure 0006955186
焦点位置決定部12は、メモリから読み出した所望の特性に対応する1対以上の初期焦点座標に対し、指向性の方向から特定できる回転行列を座標毎にかけた上で、仮想音源の座標を座標毎に加算することで、全ての焦点座標を計算する。
焦点位置決定部12は、多重極子の識別子と、この多重極子を構成する焦点座標と、各焦点座標の極性とを対応づけて出力する。
なお、クアドラポールなど、2より多数のモノポール音源からなる多重極音源についても、回転行列で回転させて新たな座標を算出することで指向性の回転に対応したモノポール音源の座標を計算する。
図3を参照して、本発明の実施の形態に係る焦点位置決定部12による焦点位置決定処理を説明する。焦点位置決定部12は、極性の異なる一対以上の初期焦点座標について図3の処理を行い、それ以外の初期焦点座標については、初期焦点座標を焦点座標として出力する。
まずステップS11において焦点位置決定部12は、仮想音源の座標と指向性の方向の情報を取得し、ステップS12において、メモリから所望の特定に対応する1以上の初期焦点の情報を読み出す。
次に、ステップS12で読み出した各初期焦点について、焦点位置決定部12は、ステップS13およびステップS14の処理を繰り返す。ステップS13において焦点位置決定部12は、処理対象の対象焦点座標に、ステップS11で取得した指向性の方向から特定される回転行列をかける。ここで用いられる対象焦点座標は、仮想音源に対する相対座標である。ステップS14において焦点位置決定部12は、仮想音源の座標に、ステップS13により回転行列をかけた後の座標を加算して、指向性を考慮した焦点座標を決定する。
ステップS12で読み出した各初期焦点について、ステップS13およびステップS14の処理が終了すると、焦点位置決定部12は処理を終了する。
なお、ステップS13ないしステップS14の処理は、各焦点に対して行われればよく、どのような順序で行われても良い。
図4ないし図6を参照して、焦点位置決定部12の処理のシミュレーション結果を説明する。図4は、直線状スピーカアレイと、初期焦点を示す。直線状スピーカアレイは、(−2,0)から、(2,0)に配設され、一対の初期焦点座標は、(0,1−0.0345)および(0,1+0.0345)である。このとき、仮想音源の座標は、(0,1)である。この際の音場は、図4に示すように、左右対称に形成され、指向性がない。
焦点位置決定部12は、このような初期焦点座標に対して、式(1)で特定される回転行列かける。図5に示すように、初期焦点座標(1,1.0345)の仮想音源座標(0.0,1.0)に対する相対座標は、(0.0,0.0345)となる。焦点位置決定部12は、初期焦点座標の仮想音源座標に対する相対座標に対して、回転行列をかけ、仮想音源座標を加算することにより、回転後の座標(0.0172,1.0299)を得る。もう一方の初期焦点座標(0,1−0.0345)に対しても同様に処理することにより、焦点位置決定部12は、回転後の座標(−0.0172,0.9701)を得る。
図6は、図5の計算によって得られた回転後の座標における音場を示す。各モノポール座標は、図4と比べて時計回りに回転され、指向性が実現されている。
焦点位置決定部12によって、各初期焦点について、指向性を考慮した焦点座標が算出されると、フィルタ係数演算部14により、処理される。
円調和係数変換部13は、円調和係数から、焦点座標を含む多重極子に与える重みを算出する。
円調和係数変換部13は、円調和級数を解析的に変換することで各焦点音源に与える重みを決定し、現実に存在する音源が持つ指向特性を有する仮想音像の生成を実現する。円調和係数変換部13は、焦点位置決定部12が出力した各焦点座標を含む各多重極子に与える重みを、算出する。
円調和係数変換部13は、式(3)により、多重極子に与える重みを算出する。
Figure 0006955186
式(3)におけるm,nは、それぞれ、x軸方向およびy軸方向に関する音場の偏微分の次数であるが、mとnの組み合わせは重複しないので、単なるインデックスとして扱っても良い。
円調和係数変換部13は、円調和係数を、適宜取得する。例えば、円調和係数を、外部プログラムから受け取るものでも良いし、指向性を測定する対象となる音源を中心とする円上に配置された複数のマイクロフォンで観測して得られたものでも良い。また円調和係数は、別途設けられたメモリに事前に格納され、必要に応じて円調和係数変換部13によって読み出されても良い。
ここで、円調和係数から多重極子の重みを出力するための式(3)の導出を説明する。まず、xy平面上の原点に、任意の指向性を持つ音源を仮定し、この音源が生成する音場をS(x)とする。この音場を原点でテイラー展開した時、単位円上の点x=(cosα,sinα)における音場は次式で与えられる。
Figure 0006955186
一方、任意の音場は円調和展開により式(5)で表現できる。
Figure 0006955186
複素正弦波についてオイラーの公式を適用した後、νに関する二項展開を行うことで次式のように変形する。
Figure 0006955186
さらに、式(4)と式(6)の係数を比較することにより、重み係数を、式(3)の通りに算出することができる。
図7を参照して、円調和係数変換部13による円調和係数変換処理を説明する。
円調和係数変換部13は、焦点位置決定部12が出力した各多重極子について、ステップS21の処理を行う。ステップS21において円調和係数変換部13は、式(3)に従って、円調和係数から多重極子の重みを算出する。
フィルタ係数演算部14は、スピーカアレイの各スピーカについて、焦点座標と、焦点座標の極性と、多重極子に与える重みから、スピーカに与える重み付き駆動関数を演算する。フィルタ係数演算部14は、直線状スピーカアレイの各スピーカについて、焦点位置決定部12により決定された焦点座標のそれぞれから、入力音響信号Iに畳み込む重み付き駆動関数を算出する。フィルタ係数演算部14は、焦点座標のそれぞれを用いて駆動関数を算出し、多重極子のそれぞれについて、多重極子を構成する各焦点座標の極性と駆動関数から算出した合成駆動関数と、多重極子に与える重みから、スピーカに与える重み付き駆動関数を演算する。ここでフィルタ係数演算部14は、多重極子に含まれる各焦点座標について、焦点座標の極性と駆動関数をかけた関数を加算して、多重極子の合成駆動関数を算出する。またフィルタ係数演算部14は、多重極子のそれぞれについて算出した合成駆動関数に、多重極子に与える重みをかけて加算して、重み付き駆動関数を算出する。
フィルタ係数演算部14は、まず、所定のスピーカについて重み付き駆動関数を算出する際、各焦点について、式(7)により駆動関数を算出する。
Figure 0006955186
次に、フィルタ係数演算部14は、所定の多重極子について、この多重極子に属する焦点音源の極性と、式(7)で算出された各焦点の駆動関数から、式(8)による合成駆動関数を算出する。
Figure 0006955186
フィルタ係数演算部14は、多重極子のそれぞれについて、式(8)で算出した合成駆動関数に、円調和係数変換部13で算出した重みをかけて、式(9)により、重み付き駆動関数を算出する。
Figure 0006955186
次に図8を参照して、フィルタ係数演算部14によるフィルタ係数演算処理を説明する。ここで、図2に示す多重極子および焦点が与えられた場合の計算式について、図9を参照して説明する。
まずステップS31において、フィルタ係数演算部14は、焦点位置決定処理で決定された各焦点座標を取得する。このときフィルタ係数演算部14は、各焦点の極性と、多重極子を構成する焦点座標の関係を、併せて取得する。
フィルタ係数演算部14は、ステップS32からステップS37の処理を繰り返して、各スピーカに対する重み付き駆動関数を算出する。ステップS32においてフィルタ係数演算部14は、対象スピーカの重み付き駆動関数をゼロで初期化する。
フィルタ係数演算部14は、各焦点について、ステップS33の処理を繰り返す。ステップS33においてフィルタ係数演算部14は、対象焦点の座標を用いて駆動関数を算出する。図9に示す例において、フィルタ係数演算部14は、式E11ないしE13を、各焦点の駆動関数を算出する。
フィルタ係数演算部14は、各多重極子について、ステップS34ないしステップS36の処理を繰り返して、各多重極子の合成駆動関数を算出する。ステップS34においてフィルタ係数演算部14は、処理対象の多重極子の合成駆動関数を初期化する。
フィルタ係数演算部14は、処理対象の多重極子に含まれる各焦点について、ステップS35の処理を行う。ステップS35においてフィルタ係数演算部14は、対象焦点の極性を用いて、ステップS33で算出した対象焦点の駆動関数を、合成駆動関数に加算する。図9に示す例において、フィルタ係数演算部14は、多重極子M1について、式E21を算出し、多重極子M2について、式E22を算出する。
ステップS36においてフィルタ係数演算部14は、ステップS35で算出された合成駆動関数に、円調和係数変換部13で算出された重みをかけて、重み付き駆動関数を算出する。図9に示す例において、フィルタ係数演算部14は、多重極子M1について算出した式E21に、多重極子M1の重みをかけた関数と、多重極子M2について算出した式E22に、多重極子M2の重みをかけた関数と加算して、式E31の重み付き駆動関数を算出する。
ステップS37においてフィルタ係数演算部14は、各多重極子について計算した後に得られた重み付き駆動関数を、対象スピーカに与える重み付き駆動関数として出力する。
フィルタ係数演算部14により、直線状スピーカアレイの各スピーカに対する重み付き駆動関数が算出されると、畳み込み演算部15が、入力音響信号Iに、重み付き駆動関数を畳み込むことにより、各スピーカに与える出力音響信号Oを算出する。
畳み込み演算部15は、直線状スピーカアレイの各スピーカについて、入力音響信号Iに、スピーカに対応する重み付き駆動関数を畳み込んで、スピーカへの出力音響信号Oを出力する。畳み込み演算部15は、所定のスピーカについて、このスピーカに対応する重み付き駆動関数を、入力音響信号Iに畳み込むことにより、このスピーカに対する出力音響信号Oを得る。畳み込み演算部15は、各スピーカについて同様の処理を繰り返し、各スピーカに対する出力音響信号Oを得る。
図10を参照して、畳み込み演算部15による畳み込み演算処理を説明する。
畳み込み演算部15は、ステップS41およびステップS42の処理を、直線状スピーカアレイの各スピーカに対して繰り返す。ステップS41において畳み込み演算部15は、フィルタ係数演算部14から、処理対象の対象スピーカの重み付き駆動関数を取得する。ステップS42において入力音響信号Iに、ステップS41で取得した重み付き駆動関数を畳み込み、出力音響信号Oを取得する。
各スピーカについてステップS41ないしステップS42の処理が終了すると、畳み込み演算部15は、処理を終了する。なお、ステップS41ないしステップS42の処理は、各スピーカに対して行われればよく、どのような順序で行われても良い。
本発明の実施の形態に係る音響信号処理装置1は、予め、初期焦点座標に回転をかけて、所望の指向性を実現する焦点座標を算出して、各焦点座標に対して、各スピーカに対応する重み付き駆動関数を算出する。音響信号処理装置1は、入力音響信号Iに対して、各スピーカに対応する重み付き駆動関数を畳み込むことにより、各スピーカへの出力音響信号Oを得る。この重み付き駆動関数は、各多重極子毎に、円調和係数から変換された重みが与えられているので、円調和係数を適宜設定することにより、各スピーカへの出力音響信号Oを、任意に調節することが可能になる。このように本発明の実施の形態に係る音響信号処理装置1は、楽器などの音源が持つ指向性をモデル化し、多重極子を重ね合わせて任意の指向特性を実現することができる。
(その他の実施の形態)
上記のように、本発明の実施の形態によって記載したが、この開示の一部をなす論述および図面はこの発明を限定するものであると理解すべきではない。この開示から当業者には様々な代替実施の形態、実施例および運用技術が明らかとなる。
本発明はここでは記載していない様々な実施の形態等を含むことは勿論である。従って、本発明の技術的範囲は上記の説明から妥当な請求の範囲に係る発明特定事項によってのみ定められるものである。
1 音響信号処理装置
10 メモリ
11 焦点データ
12 焦点位置決定部
13 円調和係数変換部
14 フィルタ係数演算部
15 畳み込み演算部
I 入力音響信号
O 出力音響信号

Claims (7)

  1. 入力音響信号を、仮想音源を実現するための複数のスピーカを並べて構成したスピーカアレイの各スピーカへの出力音響信号に変換する音響信号処理装置であって、
    複数の初期焦点座標と、仮想音源の座標および指向性の方向を取得して、
    前記複数の初期焦点座標のうち極性の異なる一対の初期焦点座標について、前記仮想音源の座標に基づいて、前記初期焦点座標に前記指向性の方向から特定される回転行列をかけて、焦点座標を決定する焦点位置決定部と、
    円調和係数から、前記焦点座標を含む多重極子に与える重みを算出する円調和係数変換部と、
    前記スピーカアレイの各スピーカについて、前記焦点座標と、前記焦点座標の極性と、前記多重極子に与える重みから、前記スピーカに与える重み付き駆動関数を演算するフィルタ係数演算部と、
    前記スピーカアレイの各スピーカについて、前記入力音響信号に、前記スピーカに対応する重み付き駆動関数を畳み込んで、前記スピーカへの出力音響信号を出力する畳み込み演算部
    を備えることを特徴とする音響信号処理装置。
  2. 前記円調和係数変換部は、式(1)により、多重極子に与える重みを算出する
    ことを特徴とする請求項1に記載の音響信号処理装置。
    Figure 0006955186
  3. 前記フィルタ係数演算部は、前記焦点座標のそれぞれを用いて駆動関数を算出し、前記多重極子のそれぞれについて、前記多重極子を構成する各焦点座標の極性と駆動関数から算出した合成駆動関数と、前記多重極子に与える重みから、前記スピーカに与える重み付き駆動関数を演算する
    ことを特徴とする請求項1に記載の音響信号処理装置。
  4. 前記フィルタ係数演算部は、前記多重極子に含まれる各焦点座標について、前記焦点座標の極性と駆動関数をかけた関数を加算して、前記多重極子の前記合成駆動関数を算出する
    ことを特徴とする請求項3に記載の音響信号処理装置。
  5. 前記フィルタ係数演算部は、前記多重極子のそれぞれについて算出した合成駆動関数に、前記多重極子に与える重みをかけて加算して、前記重み付き駆動関数を算出する
    ことを特徴とする請求項3に記載の音響信号処理装置。
  6. 入力音響信号を、仮想音源を実現するための複数のスピーカを並べて構成したスピーカアレイの各スピーカへの出力音響信号に変換する音響信号処理方法であって、
    複数の初期焦点座標と、仮想音源の座標および指向性の方向を取得するステップと、
    前記複数の初期焦点座標のうち極性の異なる一対の初期焦点座標について、前記仮想音源の座標に基づいて、前記初期焦点座標に前記指向性の方向から特定される回転行列をかけて、焦点座標を決定するステップと、
    円調和係数から、前記焦点座標を含む多重極子に与える重みを算出するステップと、
    前記スピーカアレイの各スピーカについて、前記焦点座標と、前記焦点座標の極性と、前記多重極子に与える重みから、前記スピーカに与える重み付き駆動関数を演算するステップと、
    前記スピーカアレイの各スピーカについて、前記入力音響信号に、前記スピーカに対応する重み付き駆動関数を畳み込んで、前記スピーカへの出力音響信号を出力するステップ
    を備えることを特徴とする音響信号処理方法。
  7. コンピュータに、請求項1ないし5のいずれか1項に記載の音響信号処理装置として機能させるための音響信号処理プログラム。
JP2020503603A 2018-03-01 2019-02-28 音響信号処理装置、音響信号処理方法および音響信号処理プログラム Active JP6955186B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2018036186 2018-03-01
JP2018036186 2018-03-01
PCT/JP2019/007754 WO2019168083A1 (ja) 2018-03-01 2019-02-28 音響信号処理装置、音響信号処理方法および音響信号処理プログラム

Publications (2)

Publication Number Publication Date
JPWO2019168083A1 JPWO2019168083A1 (ja) 2021-02-04
JP6955186B2 true JP6955186B2 (ja) 2021-10-27

Family

ID=67806286

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020503603A Active JP6955186B2 (ja) 2018-03-01 2019-02-28 音響信号処理装置、音響信号処理方法および音響信号処理プログラム

Country Status (4)

Country Link
US (1) US11122363B2 (ja)
EP (1) EP3761665B1 (ja)
JP (1) JP6955186B2 (ja)
WO (1) WO2019168083A1 (ja)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6865440B2 (ja) * 2017-09-04 2021-04-28 日本電信電話株式会社 音響信号処理装置、音響信号処理方法および音響信号処理プログラム
US11570543B2 (en) 2021-01-21 2023-01-31 Biamp Systems, LLC Loudspeaker polar pattern creation procedure

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5346321B2 (ja) 2010-05-20 2013-11-20 日本電信電話株式会社 音場収音再生装置、方法及びプログラム
JP5679304B2 (ja) * 2011-02-15 2015-03-04 日本電信電話株式会社 多重極スピーカ群とその配置方法と、音響信号出力装置とその方法と、その方法を用いたアクティブノイズコントロール装置と音場再生装置と、それらの方法とプログラム
DE102012200512B4 (de) * 2012-01-13 2013-11-14 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Vorrichtung und Verfahren zum Berechnen von Lautsprechersignalen für eine Mehrzahl von Lautsprechern unter Verwendung einer Verzögerung im Frequenzbereich
WO2016074734A1 (en) * 2014-11-13 2016-05-19 Huawei Technologies Co., Ltd. Audio signal processing device and method for reproducing a binaural signal

Also Published As

Publication number Publication date
JPWO2019168083A1 (ja) 2021-02-04
EP3761665A4 (en) 2021-12-01
WO2019168083A1 (ja) 2019-09-06
EP3761665B1 (en) 2022-05-18
EP3761665A1 (en) 2021-01-06
US11122363B2 (en) 2021-09-14
US20210006892A1 (en) 2021-01-07

Similar Documents

Publication Publication Date Title
JP7551795B2 (ja) オーディオ再生のためのオーディオ音場表現のデコードのための方法および装置
JP6865440B2 (ja) 音響信号処理装置、音響信号処理方法および音響信号処理プログラム
Zotter et al. A beamformer to play with wall reflections: The icosahedral loudspeaker
EP2609759A1 (en) Method and device for enhanced sound field reproduction of spatially encoded audio input signals
JP6955186B2 (ja) 音響信号処理装置、音響信号処理方法および音響信号処理プログラム
JP6970366B2 (ja) 音像再現装置、音像再現方法及び音像再現プログラム
Murillo et al. Evaluation of ambisonics decoding methods with experimental measurements
Tsutsumi et al. Directivity synthesis with multipoles comprising a cluster of focused sources using a linear loudspeaker array
WO2018211984A1 (ja) スピーカアレイ、および信号処理装置
Hur et al. Techniques for synthetic reconfiguration of microphone arrays
WO2021024806A1 (ja) スピーカアレイ、信号処理装置、信号処理方法および信号処理プログラム
Zotter Virtual Acoustics under Variable Ambisonic Perspectives
JP2018121225A (ja) 音響再生装置
Becker Franz Zotter, Markus Zaunschirm, Matthias Frank, and Matthias Kronlachner
AU2020201419A1 (en) Method and device for decoding an audio soundfield representation
Murillo Gomez et al. Evaluation of ambisonics decoding methods with experimental measurements
JP2009139615A (ja) 音響再生装置、音響再生方法、音響再生プログラム、及び音響再生システム

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200828

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200828

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20210831

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210913

R150 Certificate of patent or registration of utility model

Ref document number: 6955186

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150