Perfectly diffuse sound fields play an important role in architectural acoustics and there are established theoretical characterizations of perfect diffuseness. Although sound fields in real rooms are diffuse to some extent, they are not perfectly diffuse, and therefore theories are required to describe pseudo-perfectly diffuse sound fields. Here, we aim to spatially characterize pseudo-perfect diffuseness via directional characterization of that, finite-degree spherical harmonic diffuseness. Our results show that finite-degree diffuse sound fields yield local spatial diffuseness, suggesting that spatial pseudo-perfect diffuseness is characterized using the effective radius of diffuseness.

Perfectly diffuse sound fields are theoretically characterized by using their directional and spatial properties. For example, a random-wave model1,2 composed of an infinite number of isotropic and incoherent plane waves simulates perfectly diffuse sound fields based on the directionally random properties. The random-wave model can also be defined by using random coefficients of spherical wave functions.3,4 The random-wave model also enables spatial characterization of the perfect diffuseness. The spatial perfect diffuseness is defined by exploiting the correspondence of the two-point covariance of sound pressures at a pair of positions with the imaginary part of the Green's function of the field between the positions5 or with the sinc function2 (the imaginary part of the Green's function in a free sound field). The above models and characterizations establish perfect diffuseness theoretically and provide the foundations for architectural acoustics such as the measurements of absorption coefficients using reverberation rooms.6 

In contrast, sound fields in real environments never achieve perfect diffuseness; therefore, previous studies in room acoustics have attempted to evaluate the extent of diffuseness in a room.7–10 Although they are not perfectly diffuse, real sound fields are assumed to be diffuse to some extent. This therefore requires another theory describing pseudo-perfect diffuseness.

We have recently characterized pseudo-perfect diffuseness by discretizing or truncating the directional properties of perfect diffuseness.11,12 In these studies, discrete random-wave models composed of a finite number of random plane waves simulate pseudo-perfectly diffuse sound fields, and finite-degree spherical harmonic diffuseness, which is defined using random coefficients of truncated spherical harmonics, characterizes pseudo-perfect diffuseness. Nonetheless, spatial characterization of pseudo-perfect diffuseness remain unclear.

Here, we show that pseudo-perfect diffuseness can be defined spatially as local diffuseness, spatial diffuseness valid for a restricted area, via finite-degree spherical harmonic diffuseness.

Section 2 preliminarily introduces models of perfectly diffuse sound fields, defines spatially characterized perfect diffuseness, and introduces finite-degree diffuseness, directionally characterized pseudo-perfect diffuseness. Section 3 then shows that spatial diffuseness is locally valid in sound fields with finite-degree diffuseness and suggests defining spatial pseudo-perfect diffuseness as local diffuseness with an effective radius.

There are two established methods of modeling perfectly diffuse sound fields, as shown below.

A (continuous) random-wave model simulating perfectly diffuse sound fields consists of an infinite number of isotropic and incoherent plane waves,1 which is rigorously formulated as follows:5,
(1)
where P ( k , r ) denotes sound pressure for a wavenumber k [rad/m] and a position r. A ( k , θ , ϕ ) is a complex-valued function representing plane wave amplitudes arriving from a direction indicated by a unit vector,
(2)
To simulate an isotropic and incoherent mix of plane waves, the distributional function A ( k , n ) is assumed to provide a spherical Gaussian random field with the following statistical properties5 for k 0:
(3)
where E [ · ] represents an ensemble average (i.e., the expected value for repeatedly synthesized sound fields), ( · ) ¯ denotes taking a complex conjugate, and δ S 2 ( · ) denotes the Dirac delta on the unit sphere in a 3-dimensional space. The last two equations respectively represent the pseudo covariance and the covariance of the distribution of plane wave amplitudes, although the Dirac-delta covariance may be difficult to physically understand as also indicated in an optical literature.13 Notably, the above (integral) formulation is more rigorous than a formulation including the infinite sum of plane waves,14 where the arrival directions of the plane waves are provided by words rather than equations. The latter formulation is not employed here also because the question remains whether it is equivalent to the following spherical harmonic formulation of perfectly diffuse sound fields. We hereafter assume k 0 as well and omit the argument k without notice.
The random-wave model can also be constructed using random coefficients a l , m of spherical wave functions3,4
(4)
where j l ( · ) denotes the lth-order spherical Bessel function of the first kind, Y l , m ( · ) denotes the real-valued spherical harmonics of the lth degree and the mth order, with l , m : = l = 0 m = l l. The coefficients a l , m are assumed to have Gaussian distributions with the following statistical properties:3,4,11
(5)
(5a)
(5b)
(5c)
In short, the real and the imaginary parts of each coefficient are zero means and independent, and have the same variance of 1 / 2. The above statistical properties of a l , m define infinite-degree spherical harmonic diffuseness, which represents directionally characterized perfect diffuseness.
Spatial diffuseness is defined by the statistical properties of sound pressure yielded using the random-wave model. First, sound pressure in Eq. (4), given by the linear sum of Gaussian random coefficients, has a Gaussian distribution. Second, the expected value and the two-point pseudo covariance of sound pressures are both zeros,
(6a)
(6b)
Finally, the two-point covariance results in a sinc function2,5
(6c)
where the addition theorems15 of spherical harmonics and spherical Bessel functions allow the last transformation. Notably, the sinc function corresponds to the imaginary part of a Green's function in free sound fields. In fact, Weaver and Lobkis have shown that the two-point covariance in diffuse sound fields includes the imaginary part of the Green's function between the two positions and suggested defining the spatial diffuseness by the correspondence between the covariance and the Green's function.5 The sinc covariance derived above thus represents a special case of the diffuseness in a free sound field.

Finite-degree spherical harmonic diffuseness, directionally characterized pseudo-perfect diffuseness, is defined by simply truncating the infinite-degree spherical harmonic diffuseness.11,12 Namely, Gaussian random coefficients α l , m of 0 l L and | m | l with the statistical properties in Eqs. (5a)–(5c) define the Lth-degree spherical harmonic diffuseness.

The Lth-degree spherical harmonic diffuseness is achieved by using a discrete random-wave model of the Lth-degree. The model is composed of a finite number of Gaussian random plane waves with the following statistical properties:11,12
(7a)
(7b)
(7c)
for q , q = 1 , , Q, where A q ( k ) denotes the amplitude of the qth plane wave, Q denotes the number of plane waves, and w q denotes the weight of a 2Lth-degree spherical quadrature rule.
A spherical quadrature rule of the 2Lth exactitude degree is a set of points and weights { n q ( θ q , ϕ q ) ; w q } q = 1 Q on the sphere accurately calculating the integrals of spherical functions, yielding the following equality:16,
(8)
where F 2 L ( θ , ϕ ) is a function written by the linear sum of spherical harmonics up to the 2Lth degree,
(9)
which can represent any polynomial function on the sphere of degree at most 2L.17 In addition, the weights of the quadrature rule are all positive: w q > 0, applied 4 π-normalization: q w q = 4 π , and the sampling points of the quadrature rule { n q ( θ q , ϕ q ) } q = 1 Q determine the arrival directions of the plane waves in the discrete random-wave model. Importantly, the following Lth-degree discrete orthonormality of spherical harmonics is enabled by the 2Lth-degree quadrature rules16 
(10)
which is necessary for the coefficients (13) to satisfy the spherical harmonic diffuseness up to the Lth degree.11,12
Sound pressure in this model can then be written as follows:
(11)
The sound pressure can also be represented by the linear sum of spherical wave functions
(12)
with
(13)

Note that the above models of diffuse sound fields simulate sound fields generated by multiple sound sources emitting Gaussian noise5 rather than room impulse responses.

This section shows that the statistical properties of sound pressure yielded using discrete random-wave models can characterize pseudo-perfect diffuseness spatially. First, we theoretically show that the expected value and the two-point pseudo covariance of sound pressures always result in zeros regardless of the diffuseness degree, L. Following the definition of spatial diffuseness suggested by Weaver and Lobkis,5 we then show numerically that the two-point covariance locally corresponds with the imaginary part of the Green's function between the two positions and suggest defining local diffuseness, which represents spatially characterized pseudo-perfect diffuseness, by exploiting such correspondence.

The amplitudes of the plane waves in the discrete model have Gaussian distributions, and thus sound pressure (11), written as the linear sum of the amplitudes, also has a Gaussian distribution. The expected value and the two-point pseudo covariance of P ( k , r ) are zeros,
(14)
(15)

These results show that only the zero statistical properties derived above cannot characterize spatial pseudo-perfect diffuseness because they hold for any pair of positions, r and r , regardless of the diffuseness degree, L. In addition, these zero properties are identical to those in the (continuous) random-wave model [see Eqs. (6a) and (6b)].

First, we calculated the two-point covariance in pseudo-perfectly diffuse sound fields synthesized using the Lth-degree discrete random-wave models under free conditions without scatterers or reflectors present in the field, as follows:
(16)
where O denotes the origin and τ is the index indicating the number of syntheses of sound fields with Ts = 50 000 denoting the total number of syntheses. The random amplitudes A q with the statistical properties as in Eqs. (7c) to (7b) and the arrival directions n q of the plane waves were generated by using a numerical quadrature rule of the 2Lth degree.18 The frequency and the sound speed were respectively 880 Hz and 345 m/s, yielding k 16.0 rad/m. The positions r were on 64 straight lines through the origin, whose directions were selected using random points on a unit hemisphere.

The upper panel of Fig. 1 shows the real part of the two-point covariance computed using discrete random-wave models of select degrees, (a) L = 8, (b) L = 16, (c) L = 24, and (d) L = 32. Colored lines denote the two-point covariance on randomly selected straight lines and circles denote the sinc function, sin ( k r ) / k r. The two-point covariance was scaled by 1 / 4 π. The lower panel of Fig. 1 shows the absolute error between the sinc function and the two-point covariance above.

Fig. 1.

(Upper panel) Colored lines show the real part of two-point covariance computed between the origin and the positions on randomly selected 64 straight lines thorough the origin in sound fields repeatedly synthesized using discrete random-wave models of select diffuseness degrees, (a) L = 8, (b) L = 16, (c) L = 24, and (d) L = 32. The number of sound field syntheses was 50 000. The two-point covariance was scaled by 1 / 4 π. Circles denote the sinc function, sin ( k r ) / k r , with k 16.0 rad/m. (Lower panel) Absolute errors between the sinc function and the two-point covariance above. Colors of the lines therein correspond to those in the above figures.

Fig. 1.

(Upper panel) Colored lines show the real part of two-point covariance computed between the origin and the positions on randomly selected 64 straight lines thorough the origin in sound fields repeatedly synthesized using discrete random-wave models of select diffuseness degrees, (a) L = 8, (b) L = 16, (c) L = 24, and (d) L = 32. The number of sound field syntheses was 50 000. The two-point covariance was scaled by 1 / 4 π. Circles denote the sinc function, sin ( k r ) / k r , with k 16.0 rad/m. (Lower panel) Absolute errors between the sinc function and the two-point covariance above. Colors of the lines therein correspond to those in the above figures.

Close modal
As shown in Figs. 1(e)–1(h), the two-point covariance remarkably corresponds to the sinc function for r < 1, r < 2, r < 3, and r < 4, respectively. This implies that the radius of the area of such correspondence can be connected with the diffuseness degrees L and the wavenumber k as follows:
(17)
Next, we calculated the two-point covariance in pseudo-perfectly diffuse sound fields with a sound-hard spherical scatterer of radius, r sc = 0.6 m, at the origin. (Unless otherwise specified, the calculation conditions are the same as in the above experiment.) The synthesized sound fields can be written as follows:
(18)
where j l ( · ) is the spherical Bessel functions of the first kind; h l ( · ) is the spherical Hankel functions of the first kind; and their first derivatives are represented by j l ( · ) and h l ( · ), respectively, and α l , m ( τ ) is written as
(19)
where again, τ is the index for the number of syntheses of sound fields. As the truncation degree in Eq. (18), we let L tr = 49 to accurately calculate sound fields for the positions such that r < 3. We then computed the two-point covariance between r c = ( 1 , 1 , 0 ) and the positions r on z = 0 plane,
(20)
The real parts of the two-point covariance computed using discrete random-wave models of diffuseness degrees, L = 16, L = 24 , and L = 32, are shown in the upper panel of Fig. 2. The two-point covariance was scaled by k / ( 4 π ) 2 for easier comparison with the imaginary part of the Green's function for a point source at r c shown in Fig. 2(d). The Green's function was computed as follows:
(21)
with
(22)
The lower panel of Fig. 2 shows the absolute error between the imaginary part of the Green's function and the two-point covariance above.
As shown in the lower panel of Fig. 2, the area of correspondence between the Green's function and the two-point covariance extends around the reference position of covariance, r c, as the diffuseness degree increases. Here, again, the radius of such area can be written as in Eq. (17). Thus, in the Lth-degree directionally diffuse sound fields, spatial diffuseness can hold for any position, r, such that
(23)
That is, we can define spatial pseudo-perfect diffuseness as local diffuseness with the effective radius, 2 L / k , via the Lth-degree directional diffuseness. Notably, in plane wave sound fields, the effective radius of local diffuseness is constant regardless of the reference position, r c; however, in sound fields not composed only of plane waves, the effective radius will depend on the diffuseness degree at the reference position.
Fig. 2.

(Upper panel) Real part of the two-point covariance between r c = ( 1 , 1 , 0 ) and the positions on z = 0 plane in sound fields repeatedly synthesized using discrete random-wave models of select diffuseness degrees, (a) L = 16, (b) L = 24, and (c) L = 32. The two-point covariance was scaled by k / ( 4 π ) 2. (d) Imaginary part of the Green's function for a point source at r c. (Lower panel) Absolute errors between the imaginary part of the Green's function and the two-point covariance above.

Fig. 2.

(Upper panel) Real part of the two-point covariance between r c = ( 1 , 1 , 0 ) and the positions on z = 0 plane in sound fields repeatedly synthesized using discrete random-wave models of select diffuseness degrees, (a) L = 16, (b) L = 24, and (c) L = 32. The two-point covariance was scaled by k / ( 4 π ) 2. (d) Imaginary part of the Green's function for a point source at r c. (Lower panel) Absolute errors between the imaginary part of the Green's function and the two-point covariance above.

Close modal
In addition to characterizing pseudo-perfect diffuseness spatially, the present results provide insight into the sufficient number of isotropic and incoherent sound sources for generating locally diffuse sound fields, which is not very trivial.5 Namely, the number of sampling points in the 2Lth-degree quadrature rule, which composes the Lth-degree discrete random-wave model, can be regarded as the sufficient number of random sound sources. For example, we can numerically construct the 2Lth-degree spherical quadrature rules with
(24)
sampling points19 and positive weights,18 where · denotes a ceiling function.

Note that emerging the imaginary part of the Green's function from the two-point covariance of sound pressures in diffuse sound fields is different from reproducing a sound field based on theories such as higher-order ambisonics20 (HOA). For example, HOA usually deals only with interior problems (i.e., sound fields without sound sources or scatterers). Our results should therefore not be interpreted based on HOA.

Aiming to spatially characterize pseudo-perfect diffuseness, we presented locally diffuse sound fields with the effective radius, R = 2 L / k, via the Lth-degree spherical harmonic diffuseness. Our findings connect the directional and the spatial characterizations of pseudo-perfect diffuseness.

This work was supported by JSPS KAKENHI (Grant Nos. JP19H04153, JP22H00523) and Grant-in-Aid for JSPS Fellows (Grant No. 22KJ1941).

The authors have no conflicts to disclose.

The data are available from the corresponding author upon reasonable request.

1.
C. G.
Balachandran
and
D. W.
Robinson
, “
Diffusion of the decaying sound field
,”
Acustica
19
(
5
),
245
257
(
1967
).
2.
R. K.
Cook
,
R. V.
Waterhouse
,
R. D.
Berendt
,
S.
Edelman
, and
M. C.
Thompson
, Jr.
, “
Measurement of correlation coefficients in reverberant sound fields
,”
J. Acoust. Soc. Am.
27
(
6
),
1072
1077
(
1955
).
3.
S.
Yan
,
H.
Sun
,
U. P.
Svensson
,
X.
Ma
, and
J. M.
Hovem
, “
Optimal modal beamforming for spherical microphone arrays
,”
IEEE Trans. Audio. Speech. Lang. Process.
19
(
2
),
361
371
(
2011
).
4.
B.
Rafaely
,
Fundamentals of Spherical Array Processing
, 2nd ed. (
Springer
,
Cham, Switzerland
,
2019
), Sec. 7.1.
5.
R. L.
Weaver
and
O. I.
Lobkis
, “
Diffuse fields in open systems and the emergence of the Green's function (L)
,”
J. Acoust. Soc. Am.
116
(
5
),
2731
2734
(
2004
).
6.
H.
Kuttruff
,
Room Acoustics
, 6th ed. (
CRC Press
,
Boca Raton, FL
,
2016
), Sec. 8.7.
7.
E.
Meyer
and
R.
Thiele
, “
Raumakustische Untersuchungen in zahlreichen Konzertsälen und Rundfunkstudios unter Anwendung neuerer Messverfahren
” (“Room acoustic studies in numerous concert halls and radio studios using newer measurement methods”),
Acustica
6
,
425
444
(
1956
) (in German with English summary).
8.
M.
Koyasu
and
M.
Yamashita
, “
Evaluation of the degree of diffuseness in reverberation chambers by spatial correlation techniques
,”
J. Acoust. Soc. Jpn.
27
(
3
),
132
143
(
1971
) (in Japanese with English summary).
9.
K. J.
Ebeling
, “
Statistical properties of random wave fields
,” in
Physical Acoustics
, edited by
W. P.
Mason
and
R. N.
Thurston
(
Academic Press
,
Orlando, FL
,
1984
), Vol.
17
, Secs. 5.3 and 5.4, pp.
233
310
.
10.
M.
Nolan
,
E.
Fernadez-Grande
,
J.
Brunskog
,
A.
Richard
, and
C. H.
Jeong
, “
A wavenumber approach to characterizing the diffuse field conditions in reverberation rooms
,” in
Proceedings of the 22nd International Congress on Acoustics
, paper ICA2016-578, Buenos Aires, Argentina (
2016
).
11.
T.
Tanaka
and
M.
Otani
, “
A diffuseness of a sound field in the domains of spherical harmonics and plane waves
,”
Inter-Noise
268
(
6
),
2164
2171
(
2023
).
12.
T.
Tanaka
and
M.
Otani
, “
Directionally characterized pseudo-perfect diffuseness: A detailed comparison of its theoretical formulations
,”
Acoust. Sci. Technol.
(published online 2024).
13.
M. J.
Beran
and
G. B.
Parrent
,
Theory of Partial Coherence
(
Prentice-Hall
,
New York
,
1964
), Sec. 4.4.
14.
F.
Jacobsen
and
A. R.
Molares
, “
Sound power emitted by a pure-tone source in a reverberation room
,”
J. Acoust. Soc. Am.
126
(
2
),
676
684
(
2009
).
15.
F. W. J.
Olver
,
A. B.
Olde Daalhuis
,
D. W.
Lozier
,
B. I.
Schneider
,
R. F.
Boisvert
,
C. W.
Clark
,
B. R.
Miller
,
B. V.
Saunders
,
H. S.
Cohl
, and
M. A.
McClain
, “
NIST Digital Library of Mathematical Functions
” (2023), https://dlmf.nist.gov/ (Last viewed 2 December
2023
).
16.
I. H.
Sloan
, “
Polynomial interpolation and hyperinterpolation over general regions
,”
J. Approximation Theory
83
(
2
),
238
254
(
1995
).
17.
R. S.
Womersley
and
I. H.
Sloan
, “
How good can polynomial interpolation on the sphere be?
,”
Adv. Comput. Math.
14
,
195
226
(
2001
).
18.
T.
Tanaka
, “
Efficient spherical quadrature rules with positive weights [Dataset]
” (2023), doi: 10.5281/zenodo.10042399 (Last viewed 26 October
2023
).
19.
M.
Gräf
, “
Efficient algorithms for the computation of optimal quadrature points on Riemannian manifolds
,” Ph.D. thesis,
Technische Universität Chemnitz
, Chemintz, Germany (
2013
), Chap. 6.
20.
F.
Zotter
and
M.
Frank
,
Ambisonics: A Practical 3D Audio Theory for Recording, Studio Production, Sound Reinforcement, and Virtual Reality
(
Springer Nature
,
New York
,
2019
).