WO2014108365A1 - A sound-field control method using a planarity measure - Google Patents
A sound-field control method using a planarity measure Download PDFInfo
- Publication number
- WO2014108365A1 WO2014108365A1 PCT/EP2014/050081 EP2014050081W WO2014108365A1 WO 2014108365 A1 WO2014108365 A1 WO 2014108365A1 EP 2014050081 W EP2014050081 W EP 2014050081W WO 2014108365 A1 WO2014108365 A1 WO 2014108365A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- sound
- zone
- target
- acc
- planarity
- Prior art date
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04S—STEREOPHONIC SYSTEMS
- H04S7/00—Indicating arrangements; Control arrangements, e.g. balance control
- H04S7/30—Control circuits for electronic adaptation of the sound field
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R2201/00—Details of transducers, loudspeakers or microphones covered by H04R1/00 but not provided for in any of its subgroups
- H04R2201/40—Details 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/401—2D or 3D arrays of transducers
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R3/00—Circuits for transducers, loudspeakers or microphones
- H04R3/005—Circuits for transducers, loudspeakers or microphones for combining the signals of two or more microphones
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R3/00—Circuits for transducers, loudspeakers or microphones
- H04R3/12—Circuits for transducers, loudspeakers or microphones for distributing signals to two or more loudspeakers
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04S—STEREOPHONIC SYSTEMS
- H04S2400/00—Details of stereophonic systems covered by H04S but not provided for in its groups
- H04S2400/15—Aspects of sound capture and related signal processing for recording or reproduction
Definitions
- the present invention relates to a method of controlling sound output to a space and in particular sound output to a space having a target zone and a dark zone, where a predetermined audio signal is desired in the target zone but no sound (or another audio signal) is desired in the dark zone.
- At least two kinds of region must be created by the loudspeaker array - the target zone, where the sound pressure reaches a certain target level, and the dark zone, a region of cancellation where the audio programme delivered to the target zone is attenuated.
- an independent audio system then can be created by superposition.
- the energy in the zones can be controlled, either via a beam forming approach, or using an energy cancellation based optimization approach.
- brightness control is an optimized beam former for focusing the energy in a particular direction
- acoustic contrast control (ACC) is an energy cancellation method creating a region of significant attenuation.
- ACC acoustic contrast control
- Various applications of this work have been investigated for the personal audio scenario, including for PC users, aircraft passengers, and users of mobile devices.
- An alternative cancellation method known as acoustic energy difference maximization (AEDM) was proposed with a modified cost function to avoid the matrix inversion and allow for adjustment of the array control effort.
- AEDM acoustic energy difference maximization
- the energy cancellation methods can produce excellent acoustic contrast (cancellation) between the zones, offering great potential for sound zone reproduction, yet the phase in the target zone is uncontrolled. Consequently, multiple plane wave components can impinge on the zone from various directions, which may create highly self-cancelling waves or other undesirable audio artefacts.
- the synthesis approaches are able to resolve this issue, but often at the cost of some contrast performance and with high array effort. Furthermore, in being strictly specified, these methods are subject to spatial aliasing problems as frequency increases, with an upper frequency bound for accurate reproduction (which is required for good cancellation).
- Plane wave reproduction has commonly been regarded as the best way to create a target zone with these properties (due to the potential of reproducing any given sound field by superposition), and synthesis approaches have been adopted for sound zone problems where this is a particular requirement.
- Other methods have considered the manipulation of intensity in a single zone (with no corresponding cancellation region). For instance, the intensity based on adjacent microphone responses in a zone may be spatially averaged, optimizoptimized and controlled, or a plane wave may be reproduced by focusing the plane wave energy towards a point in the wavenumber domain.
- the intensity is estimated using a superdirective beam forming approach that could be applied to any microphone array (rather than depending on a particular geometry), is spatially averaged, and applied to the target zone meanwhile a dark zone is also created.
- the present invention generally relates to a novel cost function 'planarity control' for sound zone optimization, where the incoming plane wave direction with respect to the target zone is constrained over a range of angles, rather than a single one (or unconstrained). I n this way, a planar sound field can be reproduced (alongside excellent cancellation) but the optimization is free to find the best plane wave direction within the specified range.
- the invention relates to a method of controlling sound output to a space by a plurality of sound generators, the method comprising:
- the properties optimized may be an optimization of a contrast between the dark zone and the target zone.
- another parameter may be that no propagation takes place via the dark zone.
- Another parameter may be the control effort which is a measure of the power required by the loudspeakers in order to create the desired sound. Excessive sound power leads to degradations in reproduced performance due to non-linear transduction, an increased effect of room reflections and increased sensitivity to errors in the system.
- a Direction of Arrival is a direction seen from the target zone.
- a direction is determined by an angle from a predetermined direction, such as a direction perpendicular to a line through a centre of the target zone and a centre of the dark zone.
- step 1) a) comprises for each of a multiple of frequencies requiring that the sound at that frequency in the target zone is planar and homogenous across the target zone, and requiring that the plane waves at the individual frequencies are at least substantially aligned to each other.
- a plurality of Directions of Arrival are defined and steps b) and c) comprise solving where: p_b is the pressure at each of M microphones in the target zone p_d is the pressure at each of N microphones in the dark zone (( is the vector of complex source weights
- the range of Directions of Arrival is 120-240 degrees to a direction perpendicular to an axis intersecting centres of the target and dark zones.
- the range of Directions of Arrival is ⁇ 20 degrees of a target direction.
- the same ⁇ is used for all frequencies.
- step c) further requires that no sound must propagate across the dark zone.
- figure 1 illustrates a geometry of a sound zone system
- - figure 4 illustrates energy distributions over azimuth for PC (top), ACC (middle) and ACC-PM
- figure 5 illustrates target vs. achieved energy distribution over azimuth
- figure 6 illustrates target vs. achieved energy distributions.
- a sound zone system comprises an array of loudspeakers and a number of microphones sampling the sound field in each zone.
- the vectors of pressures at the microphones in each zone can likewise be written.
- a and H we consider two zones, A and H: p.!
- zone A the plant matrix containing the transfer functions between the loudspeakers and the microphones in zone A is defined as G , and the equivalent notation is used for G;,.
- Acoustic contrast is a spatially averaged summary measure for sound zone performance, and is commonly used in the cancellation literature to describe system performance.
- zone A defined by M microphones
- the spatially averaged squared pressure is
- the planarity of the sound field is a physical measure for assessing the extent to which a reproduced sound field resembles a plane wave.
- the reproduction error often used in the sound field synthesis literature to quantify the performance of sound field synthesis methods, may rate a highly planar sound field very poorly if the plane wave direction does not coincide with the specified sound field.
- the absolute angle of the incoming plane wave is not important and the planarity property has been designed to test each plane wave component impinging on the microphone array.
- ii is the unit vector associated with the ith component's direction
- control effort is the energy that the loudspeaker array requires in order to achieve the reproduced sound field. It is defined as the total array energy (sum of squared source weights) relative to a single monopole qr producing the same pressure in the target zone, and expressed in decibels as
- the ACC cost function where the ratio of the spatially averaged sound pressure levels between the bright zone and the dark zone is maximised, represents the energy cancellation approach.
- the cost function can be written as a constrained optimization problem based on minimising the dark zone pressure and constrained by the bright zone pressure and control effort:
- the minimum can be found by taking the eigenvector corresponding to the minimum eigenvalue of * . * i d * ; which is equivalent to taking the eigenvector corresponding to the maximum eigenvalue of .
- the regularization term therefore regularizes both the control effort and the numerical conditioning of the inversion of GTM G d .
- any phase distribution can be specified for I'M.
- a complex pressure is specified at each microphone; in this case, a plane wave is specified propagating through the target zone, and a pressure amplitude of zero is specified for the dark zone positions.
- the proposed cost function optimizes the acoustic planarity by modification of the ACC cost function stated in Eq. (7).
- the elements of H from Eq. (4) can be written in full, with respect to the microphones in the target zone, as where him is the steering vector between the zth incident angle with respect to the mth microphone in the zone.
- h, irgmi « (16)
- each row of H* is populated by the corresponding hi.
- 11 ⁇ ' ⁇ represents a mapping between the complex pressures at the microphones and the reproduced plane wave energy distribution over azimuth, as previously introduced in Eq. (4). Therefore, it presents us with an opportunity to include it in the cost function for the sound zone optimization, and achieve some control of the plane wave energy in the target zone.
- a weighting must be applied based on the acceptable range of incoming plane wave directions. Such a weighting can be specified in terms of the desired normalised energy distribution over DOA by means of a diagonal matrix ⁇ comprising weights, typically set between zero and one:
- the vector is populated with zeros apart from a single target direction, a plane wave impinging from the specified direction should be reproduced.
- the approach is somewhat similar to that employed in the prior art, where a mapping matrix was used to deactivate certain loudspeakers for efficient 3-D reproduction based on a non-uniform array, yet here the loudspeakers are not expressly deactivated as they may contribute to the cancellation region as well as the target plane wave direction. Nevertheless, the solution is highly efficient.
- the system is free to maximise the energy under this constraint, which is best achieved by the generation of a planar sound field, and thus the planarity is optimized.
- ⁇ is kept identical over frequency, similarity between adjacent frequency bins can be achieved.
- Figure 1 illustrates a geometry of the sound zone system, where zone A is the target zone and zone B is the dark zone.
- the outer (dashed) circle represents the loudspeaker array, and the inner circle the reproduction radius with respect to the aliasing limit of the synthesis methods.
- the directions of plane wave incidence with respect to zone A are indicated.
- the test geometry comprised a circular array with 48 equally spaced loudspeakers, of radius 1.2m (see Fig. 1), and 156 omnidirectional microphones in each zone spaced at 2.1cm and arranged to sample 30cm diameter circles.
- the microphones used for calculating the sound zone filters (setup) and those for obtaining predictions (playback) were kept spatially distinct or mismatched in order to asse s a slightly wider spatial region than the specific points sampled for setup (becoming more independent with increasing frequency).
- the target sound pressure level was set to 76dB SPL (achieved by scaling the prototype source weight vector q), which has been shown to be a comfortable listening level and has been used during listening tests based on the sound zone interference situation. This imposes an upper limit on the achievable contrast scores as we do not allow sound pressure levels below OdB, although we consider OdB to be measurable (there is no noise floor imposed).
- the minimum regularization parameter component X min was set to enforce a maximum matrix condition number of 10 10 , and the effort regularization parameter component X e g- adjusted, where necessary, to enforce a maximum effort of 20dB, with reference to a single monopole on the radius of the circle (qr, Eq. (6)).
- planarity control method was applied to the array and the results obtained under the evaluation metrics introduced in section 2.1.
- Figure 2 shows the method's performance over frequency, alongside those obtained for ACC, PM and ACC-PM under the same conditions.
- the contrast performance is very good and very consistent across the extended midrange band of 50-7000Hz.
- the term responsible for cancellation in the proposed planarity control (Eq. (18)) is unchanged from that in the ACC cost function (Eq. (7)) and the dark zone creation is therefore similar in each case, resulting in perfect cancellation as for ACC, and outperforming
- planarity scores are similar to PM and ACC-PM, and greatly improved from ACC, as the DOA constraint has removed the self cancellation artefacts from the reproduced sound field.
- the planarity control method is its robustness as a function of frequency. Where PM and ACC-PM suffer from well documented limitations to the upper frequency of accurate reproduction, depending on the loudspeaker spacing and array radius, the planarity control is able to operate well above this limit.
- the aliasing problems for PM and ACC-PM can be observed in relation to each of the evaluation metrics: from the contrast the effect of aliasing lobes passing through the dark zone can be observed, and the corresponding control effort response noted.
- the planarity response is interesting. because a planar target field is still reproduced. Even under this metric, however, these methods falter around the aliasing frequency. As there is little to distinguish between the performance characteristics of PM and ACC- PM apart from the slight improvement in contrast, ACC-PM is taken forward for further simulations.
- the target (left) and dark zones are indicated by the white circles.
- the properties of the sound field reproduced by the planarity control method are of some interest to potential users.
- ACC-PM The synthesis adopted in ACC-PM can be seen to successfully place the plane wave propagation to the specified direction, with a wider lobe at low frequency due to the poor beam former resolution (c.f. planarity scores for PM at low frequency in Fig. 2), and the higher frequency aliasing effects noticeable as side lobes.
- ACC produces plane wave energy from a wide range of azimuths as well as self-cancellation patterns. It is likely that such a field would result in an unpleasant listening experience.
- the distribution of plane wave energy directions over frequency for planarity control can be noted to conform, for the most part, to the target range, with side lobes emerging at higher frequencies above the array aliasing limit.
- Figure 4 illustrates energy distributions over azimuth for PC (top), ACC (middle) and ACC-PM (bottom), plotted at 100Hz intervals from 100-7000Hz.
- the bold dot-dash line in the uppermost plot indicates the specified window along the diagonal of ⁇ , and the directions 90_ and 180_ correspond to incoming plane wave directions of west-east and north- south, respectively, in relation to Fig. 3.
- Figure 5 illustrates target vs. achieved energy distribution over azimuth at 1kHz, using planarity control to specify the DOA, for 90° (west-east) (top), 180° (north- south) (middle) and 146° (optimal) (bottom). Maximum contrast is achieved in each case. Energy reproduced by PM is included for reference (dot-dash line)
- the window was set to allow a single azimuth (with a raised-cosine weighting to smooth the transition), and the direction varied.
- Three significant results are plotted in Fig. 5 for specified directions of 90', 146" (the optimal case for this frequency) and 180", at 1 kHz.
- the planarity control method can be seen to accurately place the plane wave to arrive from the required direction (corresponding to north- south in Fig. 3), and for the optimal case this is achieved with additional side lobe suppression, although the width of the energy lobe for PM is slightly narrower.
- Figure 6 illustrates target vs. achieved energy distributions over azimuth with lines plotted over frequency, for low (top), mid (middle) and high (bottom) frequency bands, using planarity control to constrain the DOA to a window around the optimal angle of 146 .
- Maximum contrast is achieved in each case.
- the compounding of poor beam former resolution for both setup and evaluation results in very wide lobes, at mid frequencies up to the spatial aliasing limit (approximately 2kHz) the placement is satisfactory, and at high frequencies the behaviour is rather similar to that of ACC-PM, where side lobes emerge. Even so, the main energy components remain close to the specified window and good contrast and planarity are still achieved.
- a method for optimizing the planarity in the target zone, as well as producing significant cancellation between zones, has been proposed.
- the method has been shown to be comparable to the well-established acoustic control method in terms of contrast and control effort, and superior for creating a planar field in the target zone. It also outperforms the pressure matching approach and a state of the art hybrid between pressure matching and acoustic contrast control in terms of contrast and control effort, and particularly in terms of its ability to produce a good cancellation region above the spatial aliasing region, and a planar field around this limit.
- the resolution of the microphone array beam former limits planarity performance at low frequencies below 400Hz. Definition of the weighting matrix ⁇ is very important for good performance.
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Signal Processing (AREA)
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
- Circuit For Audible Band Transducer (AREA)
Abstract
A method of deriving parameters for sound output into a space having a target zone and a dark zone. Sound may be desired output into the target zone but not into the dark zone. The parameters are optimized so that the sound within the target zone is planar and homogeneous, but the direction of arrival is allowed to vary within a predetermined range, so that the parameters may be optimized in a more suitable manner.
Description
A SOUND-FIELD CONTROL METHOD USING A PLANARITY MEASURE
The present invention relates to a method of controlling sound output to a space and in particular sound output to a space having a target zone and a dark zone, where a predetermined audio signal is desired in the target zone but no sound (or another audio signal) is desired in the dark zone.
In the mid 1990s, Druyvesteyn proposed the concept of personal sound zones. Since then, array signal processing for the creation of sound zones has emerged as a key sub-topic of sound field control.
Fundamentally, at least two kinds of region must be created by the loudspeaker array - the target zone, where the sound pressure reaches a certain target level, and the dark zone, a region of cancellation where the audio programme delivered to the target zone is attenuated. I desired, an independent audio system then can be created by superposition.
The development of sound zones in the academic literature has seen the emergence of techniques which broadly fall in to two categories. One technique, with its heritage in wave-field synthesis, is to precisely specify the sound field controlled by the array. I n this manner, a target sound field can be defined and the dark zone created by attenuating the sound pressure magnitudes over a region. Such control has been investigated both analytically based on sound field coefficient translation in 2D (using line sources) and 2.5D (using point sources), and by an optimized pressure matching (PM) to directly minimise the error between a discretised desired sound field and the field reproduced by the array. These methods often require many loudspeakers, although some recent work has given attention to this. Typically, a plane-wave is specified as the desired field, although any sound field could in principle be synthesised.
Alternatively, the energy in the zones can be controlled, either via a beam forming approach, or using an energy cancellation based optimization approach. For instance, brightness control is an optimized beam former for focusing the energy in a particular direction, and acoustic contrast control (ACC) is an energy cancellation method creating a region of significant attenuation. Various applications of this work have been investigated for the personal audio scenario, including for PC users, aircraft passengers, and users of mobile devices. An alternative cancellation method known as acoustic energy difference maximization (AEDM) was proposed with a modified cost function to avoid the matrix inversion and allow for adjustment of the array control effort.
Owing to the nature of the respective cost functions, they have distinctive performance characteristics. The energy cancellation methods can produce excellent acoustic contrast (cancellation) between the zones, offering great potential for sound zone reproduction, yet the phase in the target zone is uncontrolled. Consequently, multiple plane wave components can impinge on the zone from various directions, which may create highly self-cancelling waves or other undesirable audio artefacts. The synthesis approaches are able to resolve this issue, but often at the cost of some contrast performance and with high array effort. Furthermore, in being
strictly specified, these methods are subject to spatial aliasing problems as frequency increases, with an upper frequency bound for accurate reproduction (which is required for good cancellation).
Accordingly, recent advances in the literature have included hybrid methods that attempt to recreate a plane wave in the target zone whilst using an energy cancellation approach for the dark zone. A constrained optimization approach has been devised based on PM to constrain the dark zone energy without specifying a desired field for that zone, and similar formulations have been proposed based on hybridisations of PM with ACC and A EDM. Such hybrids have been shown to be effective at relatively low frequencies for reproducing planar sound fields with good contrast between the zones. I n each case, the target sound field must be specified by the designer. Whilst in some cases reproduction of specific plane wave components may be necessary (e.g. for spatial audio), specifying a plane wave is by no means the only satisfactory propagation pattern that the array could achieve in the context of personal audio. However, plane wave sound fields exhibit other advantageous properties such as good homogeneity of sound level across the zone and the avoidance of self cancellation properties.
Plane wave reproduction has commonly been regarded as the best way to create a target zone with these properties (due to the potential of reproducing any given sound field by superposition), and synthesis approaches have been adopted for sound zone problems where this is a particular requirement. Other methods have considered the manipulation of intensity in a single zone (with no corresponding cancellation region). For instance, the intensity based on adjacent microphone responses in a zone may be spatially averaged, optimizoptimized and controlled, or a plane wave may be reproduced by focusing the plane wave energy towards a point in the wavenumber domain.
Here, the intensity is estimated using a superdirective beam forming approach that could be applied to any microphone array (rather than depending on a particular geometry), is spatially averaged, and applied to the target zone meanwhile a dark zone is also created.
The present invention generally relates to a novel cost function 'planarity control' for sound zone optimization, where the incoming plane wave direction with respect to the target zone is constrained over a range of angles, rather than a single one (or unconstrained). I n this way, a planar sound field can be reproduced (alongside excellent cancellation) but the optimization is free to find the best plane wave direction within the specified range.
The invention relates to a method of controlling sound output to a space by a plurality of sound generators, the method comprising:
1) optimizing properties of a target zone in the space, wherein also a dark (attenuated) zone is present, by:
a) for at least one frequency requiring that the sound at that frequency in the target zone is planar across the target zone, and b) selecting the properties while requiring the planar sounds each are constrained within a range of Directions of Arrival, and
2) outputting a sound into the space by the sound generators, the sound being output on the basis of the optimized properties.
In this respect, the properties optimized may be an optimization of a contrast between the dark zone and the target zone.
As is mentioned below, another parameter may be that no propagation takes place via the dark zone.
Another parameter may be the control effort which is a measure of the power required by the loudspeakers in order to create the desired sound. Excessive sound power leads to degradations in reproduced performance due to non-linear transduction, an increased effect of room reflections and increased sensitivity to errors in the system.
A Direction of Arrival is a direction seen from the target zone. Usually, a direction is determined by an angle from a predetermined direction, such as a direction perpendicular to a line through a centre of the target zone and a centre of the dark zone.
In one embodiment, step 1) a) comprises for each of a multiple of frequencies requiring that the sound at that frequency in the target zone is planar and homogenous across the target zone, and requiring that the plane waves at the individual frequencies are at least substantially aligned to each other.
In a preferred embodiment, a plurality of Directions of Arrival are defined and steps b) and c) comprise solving
where: p_b is the pressure at each of M microphones in the target zone p_d is the pressure at each of N microphones in the dark zone
(( is the vector of complex source weights
where him is a steering vector between the ith angle of Direction of Arrival with respect to the mth microphone,
In one embodiment, the range of Directions of Arrival is 120-240 degrees to a direction perpendicular to an axis intersecting centres of the target and dark zones.
I n another embodiment, the range of Directions of Arrival is ±20 degrees of a target direction. In a preferred embodiment, the same Γ is used for all frequencies.
I n a preferred embodiment, step c) further requires that no sound must propagate across the dark zone.
In the following, preferred embodiments of the invention will be described with reference to the drawing, wherein: figure 1 illustrates a geometry of a sound zone system,
figure 2 illustrates the performance of planarity control (PC) with respect to ACC, PM and ACC-PM (a = 0:9), under the metrics of contrast (top), effort (middle) and planarity (bottom), figure 3 illustrates a sound pressure level (top row) and phase (bottom row) maps for PC, ACC and ACC-PM (a = 0:9), - figure 4 illustrates energy distributions over azimuth for PC (top), ACC (middle) and ACC-PM
(bottom), plotted at 100Hz intervals from 100-7000Hz, figure 5 illustrates target vs. achieved energy distribution over azimuth, and figure 6 illustrates target vs. achieved energy distributions.
A sound zone system comprises an array of loudspeakers and a number of microphones sampling the sound field in each zone. For a single frequency, the source weight vector is written as q = [ qi cp ... qi ]T , where there are L sources and qi describes the /th loudspeaker's complex source strength. The vectors of pressures at the microphones in each zone can likewise be written. Here, we consider two zones, A and H: p.! = [pi pi ...pu]T and pn = [pip 2 ...p Ν]γ , where there are microphones in zone A and N in zone , pm is the complex pressure at the mth microphone in zone A and and pn is the complex pressure at the «th microphone in zone H. For zone A, the plant matrix containing the transfer functions between the loudspeakers and the microphones in zone A is defined as G ,, and the equivalent notation is used for G;,. The pressure vectors for each zone are populated by the summation of the contribution of each loudspeaker at each microphone, written in vector notation as p ( =GA<? and p/< =G;, for zones A and B, respectively.
2.1. Evaluation measures In this section, the three evaluation metrics used for evaluation of the novel cost function are introduced.
These measure the achieved zone separation, the extent to which the target zone sound field exhibits characteristics of a plane wave, and the physical cost of cancellation.
2.1.1. Acoustic contrast
Acoustic contrast is a spatially averaged summary measure for sound zone performance, and is commonly used in the cancellation literature to describe system performance. For zone A defined by M microphones, the spatially averaged squared pressure is
and can be more suitably expressed in decibels as the sound pressure level relative to the threshold of hearing, pref= 2 xlO~5Pa:
PSPL„
Likewise, the pressures can be obtained. The acoustic contrast between target zone A and dark zone B is defined as the ratio of spatially averaged pressures in each zone due to the reproduction of program A: contrast AB = PsPLA ~ !¾#·½ ·
2.1.2. Planarity
The planarity of the sound field is a physical measure for assessing the extent to which a reproduced sound field resembles a plane wave. The reproduction error, often used in the sound field synthesis literature to quantify the performance of sound field synthesis methods, may rate a highly planar sound field very poorly if the plane wave direction does not coincide with the specified sound field. For sound zone reproduction at a single frequency, the absolute angle of the incoming plane wave is not important and the planarity property has been designed to test each plane wave component impinging on the microphone array. The energy distribution at the microphone array over each incoming plane wave direction w = [wi ... wi] is given by wi = ^ ψί * ψί, where
.*denotes the complex conjugate, Ψ = [Ψχ ... Ψ are the plane wave components at the z'th angle, related to the observed microphone pressures by the steering matrix I I whose elements are determined by super-directive beam forming about the microphone array, if = tip, (4)
The elements of H could alternatively be calculated using a spatial Fourier decomposition approach. The planarity metric can now be defined as the ratio between the energy due to the largest plane wave component and the total energy flux of plane wave components:
planarityji = (5)
where ii; is the unit vector associated with the ith component's direction, u is the sum of all components in the "zth direction = arg max,- wi and . denotes the inner product.
2.1.3. Control effort The control effort is the energy that the loudspeaker array requires in order to achieve the reproduced sound field. It is defined as the total array energy (sum of squared source weights) relative to a single monopole qr producing the same pressure in the target zone, and expressed in decibels as
It is a necessity in any practical system to achieve a suitably low control effort. On the one hand, it is physically related to whether a set of source weights is realizable through loudspeakers. Yet in addition, limiting the control effort results in there being less sound energy overall in the environment, leading to improved robustness to reflections in reverberant rooms, and limits the white noise gain of the system, improving robustness to other kinds of errors such as measurement noise and non-linear distortion.
2.2. Existing approaches To facilitate a comparison between the proposed cost function and existing sound field control methods, ACC, PM and their hybrid are formally introduced in the following sections.
2.2.1. Acoustic contrast control
The ACC cost function, where the ratio of the spatially averaged sound pressure levels between the bright zone and the dark zone is maximised, represents the energy cancellation approach. The cost function can be written as a constrained optimization problem based on minimising the dark zone pressure and constrained by the bright zone pressure and control effort:
= P?Prf +*, (rf pt - * : · }>2. tfq - E), (7)
where the subscripts .d and .b denote assignment of the pressure vectors with respect to the dark and bright (target) zones, respectively, B is the target sound pressure in the bright zone, and E is the maximum allowed control effort.
Taking the derivative of J and setting to zero, we obtain:
which can be rearranged as an eigenvalue problem of the form λ;(| = Aq:
I . ■ . λ,1)€Ι
The minimum can be found by taking the eigenvector corresponding to the minimum eigenvalue of * . * i d * ; which is equivalent to taking the eigenvector corresponding to the maximum eigenvalue of
. The regularization term therefore regularizes both the control effort and the numerical conditioning of the inversion of G™ Gd. In order to ensure the latter over all frequencies, the regularization parameter is split such that λι = Xmin +λφ where lmin is first set to ensure that the condition number of (^1,; 1 '<; -J) is suitably controlled to avoid numerical errors and .^is subsequently adjusted, if necessary, to ensure that E does not exceed the specified value.
2.2.2. Pressure matching
As a sound field synthesis method, any phase distribution can be specified for I'M. A complex pressure is specified at each microphone; in this case, a plane wave is specified propagating through the target zone, and a pressure amplitude of zero is specified for the dark zone positions. The optimization cost function is written to minimise the error e=Gq-d between the desired sound field d and reproduced sound field, with a control effort constraint which acts as Tikhonov regularization: i = fWf I k \ n"q - E). (10)
The solution can then be found for the optimal q:
(i l l where G = [G.'G/*] ' is the complete system plant matrix and λ = Xmin +^is split as above.
2.2.3. PM and ACC hybrid
For the hybrid solution combining PM and ACC, the PM portion of the cost function is restricted to the reproduction of the bright zone, whilst the contrast control formulation is used for cancellation. Here, for consistency with the other methods, we introduce Tikhonov regularization instead of using the pseudo-inverse as in the original work:
J = «pj prf + ( I - a)ejV„ - A ι t|"q - £), · ! 2 1 where the error eb = Gbq-db is now between the desired sound field and the reproduced field in the target zone only. The weighting a provides a tuning parameter between the pure ACC solution and the pure target PM solution, with the standard pressure matching approach (Eq. (10)) being equivalent to a = 0:5. The solution can be determined by finding the gradient of Eq. (12) and rearranging for the source weights:
and as above, λ = min +λ
3. PLANARITY CONTROL OPTIMIZATION
The proposed cost function optimizes the acoustic planarity by modification of the ACC cost function stated in Eq. (7). The elements of H from Eq. (4) can be written in full, with respect to the microphones in the target zone, as
where him is the steering vector between the zth incident angle with respect to the mth microphone in the zone. Using the super-directive (ACC) beam forming approach, Hb can be determined for each steering angle by grouping the plane wave components c in each direction (based on the plane wave Green's function, gi,c = dkr "ilM) in to a passband P and stopband S:
where p denotes passband range centred upon the z'th angle and s denotes the stopband range. We can then obtain h, = irgmi«
(16)
and each row of H* is populated by the corresponding hi.
11·'· represents a mapping between the complex pressures at the microphones and the reproduced plane wave energy distribution over azimuth, as previously introduced in Eq. (4). Therefore, it presents us with an opportunity to include it in the cost function for the sound zone optimization, and achieve some control of the plane wave energy in the target zone. In order to do this, a weighting must be applied based on the acceptable range of incoming plane wave directions. Such a weighting can be specified in terms of the desired normalised energy distribution over DOA by means of a diagonal matrix Γ comprising weights, typically set between zero and one:
where ji is the weighting applied for the z'th steering angle. The planarity optimization cost function can now be introduced:
and deriving the solution in the identical manner to Eqs. (7 - 9) above, the optimal source weights can be found to be the eigenvector corresponding to the maximum eigenvalue of
The optimization is thus constrained to maximise the sound energy in the target zone from among the potential incoming azimuths allowed by Γ. The selection of Γ is clearly a significant factor. If the vector is filled with ones, then the cost function in Eq. (18) is no different from the contrast control formulation in Eq. (7) and identical performance is achieved. If, on the other hand, the vector is populated with zeros apart from a single target direction, a plane wave impinging from the specified direction should be reproduced. The approach is somewhat similar to that employed in the prior art, where a mapping matrix was used to deactivate certain loudspeakers for efficient 3-D reproduction based on a non-uniform array, yet here the loudspeakers are not expressly deactivated as they may contribute to the cancellation region as well as the target plane wave direction. Nevertheless, the solution is highly efficient. When the range of allowable angles is suitably designed, the system is free to maximise the energy under this constraint, which is best achieved by the generation of a planar sound field, and thus the planarity is optimized. Furthermore, if Γ is kept identical over frequency, similarity between adjacent frequency bins can be achieved.
4. SIMULATIONS
The operation and performance of the planarity optimization algorithm is demonstrated in the following by means of simulations. Figure 1 illustrates a geometry of the sound zone system, where zone A is the target zone and zone B is the dark zone. The outer (dashed) circle represents the loudspeaker array, and the inner circle the reproduction radius with respect to the aliasing limit of the synthesis methods. The directions of plane wave incidence with respect to zone A are indicated.
4.1. Method The simulations were conducted in Matlab, simulating a free-field lossless anechoic environment, with each source modelled as an ideal monopole. The free-field Green's function was used to populate the plant matrices, giving the transfer function at each microphone due to a loudspeaker at distance r.
where // = 1 :21kg/m3, c = 342m/s, and k is the wavenumber ω/c.
The test geometry comprised a circular array with 48 equally spaced loudspeakers, of radius 1.2m (see Fig. 1), and 156 omnidirectional microphones in each zone spaced at 2.1cm and arranged to sample 30cm diameter circles. The microphones used for calculating the sound zone filters (setup) and those for obtaining predictions (playback) were kept spatially distinct or mismatched in order to asse s a slightly wider spatial region than the specific points sampled for setup (becoming more independent with increasing frequency). The target sound pressure level was set to 76dB SPL (achieved by scaling the prototype source weight vector q), which has been shown to be a comfortable listening level and has been used during listening tests based on the sound zone interference situation. This imposes an upper limit on the achievable contrast scores as we do not allow sound pressure levels below OdB, although we consider OdB to be measurable (there is no noise floor imposed).
Figure 2 illustrates the performance of planarity control (PC) with respect to ACC, PM and ACC-PM (a = 0:9), under the metrics of contrast (top), effort (middle) and planarity (bottom)
To set the regularization conditions, the minimum regularization parameter component Xmin was set to enforce a maximum matrix condition number of 1010, and the effort regularization parameter component Xeg- adjusted, where necessary, to enforce a maximum effort of 20dB, with reference to a single monopole on the radius of the circle (qr, Eq. (6)).
The plane wave for the PM and ACC-PM hybrid (herein simply referred to as ACC-PM) approaches was specified to travel from north to south (ψ = 180 , marked on Fig. 1), and the weighting matrix Γ was set to constrain the incoming plane wave components between 120" and 240". The weighting on the diagonal of Γ is indicated in Fig. 4 (top). The weighting a for ACC-PM was set to 0.9 to encourage good contrast performance.
4.2. Planarity optimization performance
The planarity control method was applied to the array and the results obtained under the evaluation metrics introduced in section 2.1. Figure 2 shows the method's performance over frequency, alongside those obtained for ACC, PM and ACC-PM under the same conditions. The contrast performance is very good and very consistent across the extended midrange band of 50-7000Hz. The term responsible for cancellation in the proposed planarity control (Eq. (18)) is unchanged from that in the ACC cost function (Eq. (7)) and the dark zone creation is therefore similar in each case, resulting in perfect cancellation as for ACC, and outperforming
PM and ACC-PM.
Likewise, the control effort performance tends towards that of ACC, which gives preferable performance by a small margin across the whole range, outperforming the planarity control by up to 6dB at the lowest frequencies but generally being within 3dB. Nonetheless, the effort is below OdB for much of the frequency range, and it is consistently preferable to PM and ACC-PM under the same conditions. Finally, there is a good planarity performance across frequency. Under this metric, the synthesis metrics PM and ACC-PM naturally produce optimal scores for significant portions of the frequency range. However, with the exception of the low frequency performance (due to poor resolution of the planarity steering matrix in this region) and a narrow notch at 3.6kHz, the planarity scores are similar to PM and ACC-PM, and greatly improved from ACC, as the DOA constraint has removed the self cancellation artefacts from the reproduced sound field. Perhaps the most striking characteristic of the planarity control method is its robustness as a function of frequency. Where PM and ACC-PM suffer from well documented limitations to the upper frequency of accurate reproduction, depending on the loudspeaker spacing and array radius, the planarity control is able to operate well above this limit. In fact, the aliasing problems for PM and ACC-PM can be observed in relation to each of the evaluation metrics: from the contrast the effect of aliasing lobes passing through the dark zone can be observed, and the corresponding control effort response noted. The planarity response is interesting. because a planar target field is still reproduced. Even under this metric, however, these methods falter around the aliasing frequency. As there is little to distinguish between the performance characteristics of PM and ACC- PM apart from the slight improvement in contrast, ACC-PM is taken forward for further simulations.
The optimal contrast and planarity performance obtained using planarity control can be further clarified by studying the sound pressure level and phase maps shown in Fig. 3. Figure 3 illustrates a sound pressure level (top row) and phase (bottom row) maps for PC, ACC and ACC-PM (a = 0:9). The target (left) and dark zones are indicated by the white circles. For the sound pressure level maps, white indicates a high sound pressure and black a low sound pressure at high frequencies for PM and ACC-PM, We can now confirm that the planarity control produces an ACC-like dark zone, yet replaces the north- south self-cancellation (visible across the whole of the bottom-middle plot) in the target zone with a planar field (indicated by the sharp transition in the phase response), and reduces the overall sound pressure in the environment as a consequence of the low effort score with relation to PM and ACC-PM (visible by comparing the amount of bright white in the top row of Fig. 3).
4.3. Target sound field properties
The properties of the sound field reproduced by the planarity control method are of some interest to potential users.
First, we consider the energy distribution over azimuth (with respect to the target zone) obtained for the window function used for the simulations in section 4.2. We have seen from the planarity scores (Fig. 2, bottom) and the phase distributions in the enclosure (Fig. 3, bottom) that the planarity control method is capable of creating
highly planar fields in the target zone, for single frequencies. However, these plots do not give us an indication of the range of incoming plane wave directions as a function of frequency. Therefore, in Fig. 4 the normalised energy distributions for multiple frequencies have been plotted across azimuth for planarity control, ACC and ACC-PM. This gives us a useful insight in to the planarity control's performance in relation to the existing methods. The synthesis adopted in ACC-PM can be seen to successfully place the plane wave propagation to the specified direction, with a wider lobe at low frequency due to the poor beam former resolution (c.f. planarity scores for PM at low frequency in Fig. 2), and the higher frequency aliasing effects noticeable as side lobes. Conversely, ACC produces plane wave energy from a wide range of azimuths as well as self-cancellation patterns. It is likely that such a field would result in an unpleasant listening experience. The distribution of plane wave energy directions over frequency for planarity control can be noted to conform, for the most part, to the target range, with side lobes emerging at higher frequencies above the array aliasing limit.
Figure 4 illustrates energy distributions over azimuth for PC (top), ACC (middle) and ACC-PM (bottom), plotted at 100Hz intervals from 100-7000Hz. The bold dot-dash line in the uppermost plot indicates the specified window along the diagonal of Γ, and the directions 90_ and 180_ correspond to incoming plane wave directions of west-east and north- south, respectively, in relation to Fig. 3.
Figure 5 illustrates target vs. achieved energy distribution over azimuth at 1kHz, using planarity control to specify the DOA, for 90° (west-east) (top), 180° (north- south) (middle) and 146° (optimal) (bottom). Maximum contrast is achieved in each case. Energy reproduced by PM is included for reference (dot-dash line)
To test the ability of the planarity control to reproduce a specific incoming plane wave direction, the window was set to allow a single azimuth (with a raised-cosine weighting to smooth the transition), and the direction varied. Three significant results are plotted in Fig. 5 for specified directions of 90', 146" (the optimal case for this frequency) and 180", at 1 kHz. In the middle plot (180'), the planarity control method can be seen to accurately place the plane wave to arrive from the required direction (corresponding to north- south in Fig. 3), and for the optimal case this is achieved with additional side lobe suppression, although the width of the energy lobe for PM is slightly narrower. Yet for directions perpendicular to this (west-east propagation shown), which would require a beam to be placed across the dark zone, a highly self-cancelling pattern is instead reproduced and the peak in this direction is unsatisfactory. There is no variation in the contrast between these cases and the effort difference is minimal, yet if PM had been applied, the cancellation would have been poor and the effort very high, albeit with the specified plane wave component reproduced. An interesting property of the planarity control cost function is therefore exposed: that producing high contrast is the priority of the optimization, and that where specification of the incident direction does not conflict with contrast performance, the energy can be placed precisely in the desired direction. The behaviour over frequency is clarified by Fig. 6 for a constrained window (146' ±20" with a raised cosine weighting). Figure 6 illustrates target vs. achieved energy distributions over azimuth with lines plotted over frequency, for low (top), mid (middle) and high (bottom) frequency bands, using planarity control to constrain the DOA to a window around the optimal angle of 146 . Maximum contrast is achieved in each case.
At low frequencies, the compounding of poor beam former resolution for both setup and evaluation results in very wide lobes, at mid frequencies up to the spatial aliasing limit (approximately 2kHz) the placement is satisfactory, and at high frequencies the behaviour is rather similar to that of ACC-PM, where side lobes emerge. Even so, the main energy components remain close to the specified window and good contrast and planarity are still achieved.
5. CONCLUSIONS
A method for optimizing the planarity in the target zone, as well as producing significant cancellation between zones, has been proposed. The method has been shown to be comparable to the well-established acoustic control method in terms of contrast and control effort, and superior for creating a planar field in the target zone. It also outperforms the pressure matching approach and a state of the art hybrid between pressure matching and acoustic contrast control in terms of contrast and control effort, and particularly in terms of its ability to produce a good cancellation region above the spatial aliasing region, and a planar field around this limit. The resolution of the microphone array beam former limits planarity performance at low frequencies below 400Hz. Definition of the weighting matrix Γ is very important for good performance. The ability of Γ to constrain incident plane wave directions over frequency has been demonstrated, and furthermore under the condition that it does not require propagation across the dark zone, a precise plane wave direction can be specified. The method therefore presents a compelling cost function for sound zones where the self-cancellation artefacts of energy cancellation approaches can be reduced whilst allowing more flexibility over the incident plane wave specification, yet with the potential to reproduce a wave from a single direction if required.
Claims
1. A method of controlling sound output to a space by a plurality of sound generators, the method comprising:
1) optimizing properties of a target zone in the space, wherein also a dark zone is present, by: a) for at least one frequency requiring that the sound at that frequency in the target zone is planar and homogenous across the target zone, and b) selecting the properties while requiring the planar sounds each are constrained within a range of Directions of Arrival, and
2) outputting a sound into the space by the sound generators, the sound being output on the basis of the optimized properties.
2. A method according to claim 1 , wherein step 1) a) comprises for each of a multiple of frequencies requiring that the sound at that frequency in the target zone is planar and homogenous across the target zone, and requiring that the plane waves at the individual frequencies are at least substantially aligned to each other.
3. A method according to claim 1 or 2, wherein a plurality of Directions of Arrival are defined and steps b) and c) comprise solving
J |#Prf +A,(tfH?n¼|i,, ~ /J . i A;>c ii~J¾ (18) where: p is the pressure at each of M microphones in the target zone
(( is the pressure at each of N microphones in the dark zone
4. A method according to any of the preceding claims, wherein the range of angles is 120-240 degrees to a direction perpendicular to an axis intersecting centres of the target and dark zones..
5. A method according to any of claims 1-3, wherein the range of Directions of Arrival are ±20 degrees of a target direction.
6. A method according to any of the preceding claims, wherein the same Γ is used for all frequencies.
7. A method according to any of the preceding claims, wherein step c) further requires that no sound must propagate across the dark zone.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DK14700441.0T DK2944098T3 (en) | 2013-01-10 | 2014-01-06 | Method of controlling sound emitted in a room |
EP14700441.0A EP2944098B1 (en) | 2013-01-10 | 2014-01-06 | A sound-field control method using a planarity measure |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DKPA201300015 | 2013-01-10 | ||
DKPA201300015 | 2013-01-10 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2014108365A1 true WO2014108365A1 (en) | 2014-07-17 |
Family
ID=49958444
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/EP2014/050081 WO2014108365A1 (en) | 2013-01-10 | 2014-01-06 | A sound-field control method using a planarity measure |
Country Status (3)
Country | Link |
---|---|
EP (1) | EP2944098B1 (en) |
DK (1) | DK2944098T3 (en) |
WO (1) | WO2014108365A1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112224142A (en) * | 2020-09-03 | 2021-01-15 | 中国汽车技术研究中心有限公司 | Vehicle-mounted independent acoustic space system and vehicle |
-
2014
- 2014-01-06 WO PCT/EP2014/050081 patent/WO2014108365A1/en active Application Filing
- 2014-01-06 EP EP14700441.0A patent/EP2944098B1/en active Active
- 2014-01-06 DK DK14700441.0T patent/DK2944098T3/en active
Non-Patent Citations (1)
Title |
---|
MÖLLER MARTIN BO ET AL: "A Hybrid Method Combining Synthesis of a Sound Field and Control of Acoustic Contrast", AES CONVENTION 132; APRIL 2012, AES, 60 EAST 42ND STREET, ROOM 2520 NEW YORK 10165-2520, USA, 26 April 2012 (2012-04-26), XP040574591 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112224142A (en) * | 2020-09-03 | 2021-01-15 | 中国汽车技术研究中心有限公司 | Vehicle-mounted independent acoustic space system and vehicle |
Also Published As
Publication number | Publication date |
---|---|
DK2944098T3 (en) | 2017-05-22 |
EP2944098A1 (en) | 2015-11-18 |
EP2944098B1 (en) | 2017-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Coleman et al. | Personal audio with a planar bright zone | |
Coleman et al. | Acoustic contrast, planarity and robustness of sound zone methods using a circular loudspeaker array | |
KR102024284B1 (en) | A method of applying a combined or hybrid sound -field control strategy | |
US20140270245A1 (en) | Polyhedral audio system based on at least second-order eigenbeams | |
KR102353871B1 (en) | Variable Acoustic Loudspeaker | |
Coleman et al. | Optimizing the planarity of sound zones | |
JP2006508404A (en) | Loudspeaker system for virtual sound synthesis | |
US10469945B2 (en) | Sound wave field generation based on a desired loudspeaker-room-microphone system | |
US20190014430A1 (en) | Loudspeaker-room system | |
CN104980845A (en) | Sound Wave Field Generation | |
Zotter et al. | A beamformer to play with wall reflections: The icosahedral loudspeaker | |
EP3400713B1 (en) | Loudspeaker array | |
US10547943B2 (en) | Adaptive filtering audio signals based on psychoacoustic constraints | |
Poletti et al. | Design of a prototype variable directivity loudspeaker for improved surround sound reproduction in rooms | |
Betlehem et al. | Two dimensional sound field reproduction using higher order sources to exploit room reflections | |
Rosen et al. | FIR-based symmetrical acoustic beamformer with a constant beamwidth | |
Simon Galvez et al. | Loudspeaker arrays for transaural reproduction | |
US10244317B2 (en) | Beamforming array utilizing ring radiator loudspeakers and digital signal processing (DSP) optimization of a beamforming array | |
US20200169808A1 (en) | Loudspeaker assemblies and associated methods | |
Zaunschirm et al. | An interactive virtual icosahedral loudspeaker array | |
EP2755405A1 (en) | Zonal sound distribution | |
Gao et al. | Weighted loudspeaker placement method for sound field reproduction | |
Dick et al. | A comparison of measured room acoustics metrics using a spherical microphone array and conventional methods | |
EP2944098B1 (en) | A sound-field control method using a planarity measure | |
Rasumow et al. | The impact of the white noise gain (WNG) of a virtual artificial head on the appraisal of binaural sound reproduction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 14700441 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
REEP | Request for entry into the european phase |
Ref document number: 2014700441 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2014700441 Country of ref document: EP |