Weighted pressure matching based on kernel interpolation
for sound field reproduction
Shoichi KOYAMA and Kazuyuki ARIKAWA(1)
(1)Graduate School of Information Science and Technology, The University of Tokyo, Japan, koyama.shoichi@ieee.org
ABSTRACT
A sound field reproduction method called weighted pressure matching is proposed. Sound field reproduction is
aimed at synthesizing the desired sound field using multiple loudspeakers inside a target region. Optimization-
based methods are derived from the minimization of errors between synthesized and desired sound fields, which
enable the use of an arbitrary array geometry in contrast with integral-equation-based methods. Pressure match-
ing is widely used in the optimization-based sound field reproduction methods because of its simplicity of
implementation. Its cost function is defined as the synthesis errors at multiple control points inside the target
region; then, the driving signals of the loudspeakers are obtained by solving a least-squares problem. However,
in pressure matching, the region between the control points is not taken into consideration. We define the cost
function as the regional integration of the synthesis error over the target region. On the basis of the kernel in-
terpolation of the sound field, this cost function is represented as the weighted square error of the synthesized
pressures at the control points. Experimental results indicate that the proposed weighted pressure matching
outperforms conventional pressure matching.
Keywords: Sound field reproduction, Kernel interpolation, Weighted pressure matching
1 INTRODUCTION
Sound field reproduction is aimed at synthesizing spatial sound using multiple loudspeakers (or secondary
sources). Such a technique can be applied to virtual/augmented reality audio, generation of multiple sound
zones for personal audio, and noise cancellation in a spatial region.
Sound field reproduction methods can be classified into two major categories: integral-equation-based and
optimization-based methods. The integral-equation-based methods are developed from the boundary integral
representations derived from the Helmholtz equation, such as wave field synthesis and higher-order ambison-
ics [3,19,17,1,24,12]. The optimization-based methods are derived from the minimization of a certain
cost function defined for synthesized and desired sound fields in a target region, such as pressure matching and
mode matching [16,9,7,17,4,21,11]. Many integral-equation-based methods require the array geometry of
loudspeakers to have a simple shape, such as a sphere, plane, circle, or line, and driving signals are obtained
from a discrete approximation of an integral equation. In optimization-based methods, the loudspeakers can be
placed arbitrarily, and driving signals are generally derived as a closed-form least-squares solution. In particular,
pressure matching is widely used among the optimization-based methods because of its simplicity of imple-
mentation. Pressure matching is based on synthesizing the desired pressures at a discrete set of control points
placed over the target region.
An issue of pressure matching is that the region between the control points is not taken into consideration
because of the discrete approximation. Therefore, its reproduction accuracy can deteriorate when the distribution
of the control points is not sufficiently dense. However, the smaller the number of control points, the better in
practice, because the transfer functions between the loudspeakers and control points are generally measured in
advance. We propose an optimization-based sound field reproduction method called weighted pressure matching.
We define the cost function as the regional integration of the synthesis error over the target region. On the
basis of the kernel interpolation of the sound field [20,22], this cost function is represented by the pressures
at the control points with the regional integration of kernel functions. When the same kernel functions are used
for interpolating primary and secondary sound fields, i.e., the desired field and the sound field of each loud-
speaker, respectively, the driving signal is obtained as the solution of weighted least squares problem with the
arXiv:2210.14711v1 [eess.AS] 26 Oct 2022