Abstract

Generation interior vehicle personal sound zone (PSZ) has great application prospects since it can meet the listening requirements of different people without causing interference to other people. Acoustic contrast, signal distortion, and array effort are often used to measure the performance of the PSZ system. However, these three indicators cannot reach their optimal values simultaneously. To obtain optimal performance in one aspect, it is necessary to reduce the other two indicators at a cost. Therefore, when designing the PSZ system, these three indicators must be balanced so that the comprehensive performance of the PSZ system can be optimized. Combining the power constraints of the speaker array, based on the Non-dominated Genetic Algorithm II (NSGA-II), a multiobjective optimal control method PSZ-MOGA is proposed for the generation interior vehicle PSZ, which is used to optimize the three performance indices comprehensively. In addition, the mixed speaker array is used for the reproduction of the sound field. The results show that the proposed PSZ-MOGA method improves the reproduction accuracy of the bright zones compared to the acoustic contrast control methods and also enhances the acoustic contrast (AC) while achieving lower array power consumption compared to the pressure matching method.

1 Introduction

Through the reasonable arrangement of a speaker array in the car cabin, the personal sound zone (PSZ) system meets different people's listening requirements without causing interference to others. When the speaker arrangement and external interference remain unchanged, the control algorithm determines the effect of PSZ. At present, many PSZ control algorithms have been proposed. For example, the acoustic contrast control (ACC) [1], the pressure matching (PM) [2], the ACC–PM fusion algorithm [3], the variable span trade-off filters [46], the mode matching method [7], and the amplitude matching method [8]. To evaluate the effectiveness of the PSZ system, most scholars use acoustic contrast (AC) to indicate the energy difference between bright and dark zones, signal distortion (SD) to indicate the error between the reproduced sound field and the desired sound field, and array effort (AE) to indicate the efficiency of the sound field reproduction system [9,10].

Many methods for solving the balancing of two indicators have been developed. For example, the proposed least-squares method [11] optimizes the SD while constraining the AC. In addition, Hu et al. [12] proposed a constraint optimization method to optimize the AC under the condition that the SD is constrained. However, these proposed constrained optimization methods can only be applied to achieve a balance between two indicators under some specific constraints. These methods also do not consider the problem that the actual speaker power will be constrained, thus failing to obtain the better effect of PSZ.

Cheng et al. [13] proposed the weighted sum method (WSM) and ε-constraint method (ECM) to further study the balance of AC, SD, and AE. However, this method can only find the Pareto-optimal point that is located on the boundary of the convex envelope of the target set, and there may be a situation where the optimal solution cannot be found. Cheer and Elliott [14] used a mixed speaker array to generate PSZ but did not consider the effect of speaker position or used it to optimize the overall performance of the PSZ. Peng et al. [15] proposed a method for position and number optimization of mixed speaker arrays, but the effect on PSZ is unclear. The method proposed by Zhao and Burnett [16] and Gao et al. [17] cannot further optimize the interior vehicle PSZ system with the fixed speaker position and quantities.

Stuart et al. [18] and Thipmaungprom [19] pointed out that in full-range speakers, driving multiple frequencies through the same speaker unit simultaneously tends to produce intermodulation distortion. This distortion alters the signal spectrum and waveform, affecting the quality and accuracy of the signal. The mixed speaker system can reduce the distortion by assigning different frequencies to different speaker units, resulting in clearer sound.

In this article, PSZ-MOGA is proposed to achieve the dynamic balance of AC, SD, and AE, and the accuracy of sound field reproduction is further improved by using the mixed speaker array.

Section 2 presents the theory and methodology of PSZ control. Section 3 is the method and theory proposed in this article. Section 4 presents the simulation part. An experimental part is presented in Sec. 5. Section 6 is the conclusion.

2 Principle of Personal Sound Zone Control and Evaluation Indicators

2.1 Principle of Personal Sound Zone Control.

The schematic diagram of the generation interior vehicle PSZ based on the speaker array is shown in Fig. 1. In Fig. 1, the main driver's area is the bright area, which will reproduce the target sound field in this area. The right side of the back row is the dark area, which will not be affected by the reproduction sound field of the bright area.

Fig. 1
Schematic diagram of generation interior vehicle PSZ
Fig. 1
Schematic diagram of generation interior vehicle PSZ
Close modal
The sound pressure in a single control area at a given time p(t) can be expressed as
(1)
where pi(t)(i=1M) indicates the sound pressure signal at each control point, M is the number of control points in a single control area. gB(t) and gD(t) are, respectively, the acoustic transfer function matrices from the speakers to the control points of the bright and dark zones.
g(t) is the set of transfer functions from a single control point to the speakers:
(2)
where M indicates the number of individual zone control points and L is the number of speakers:
(3)
where * denotes the time-domain convolution and s(t) is the time-domain representation of the vector of the driving signals of L speakers:
(4)
The sound pressure P in a single control region in the frequency domain is
(5)

2.2 Evaluation Indicators.

AC is the acoustic energy contrast between bright and dark zones, which is used to evaluate the sound source isolation effect of PSZ, and its expression from Ref. [1] is
where RB and RD are the correlation matrices, RB=GBHGB; RD=GDHGD, εB, and εD denote the acoustic potential energies of the bright and dark zones, respectively; ()H denote the Hermitian transpose of the matrices; and MB and MD denote the number of control points in the bright and dark zones, respectively.
SD, which represents the normalized mean square error between the reproduced sound field and the desired sound field in the bright zone, reflects the accuracy of the reproduced sound field and is expressed as Eq. (7) in Ref. [13]
(7)
where Pdes is the target reproduced sound field, Pdes=[PBdes;0]; PBdes is the ideal sound pressure in the bright zone; and α is the adjustment factor. Take the partial derivative of α in Eq. (7) and make the partial derivative zero to get Eq. (8):
(8)
Substituting Eq. (8) into Eq. (7):
(9)
The simplification operation of Eq. (9) yields Eq. (10):
(10)
AE represents the ratio of the speaker array power to the power of a reference source. The value of AE is from Ref. [20]:
(11)
where |Sr|2 indicates the source intensity of a reference source that produces the same acoustic energy as the speaker array in the bright zone.

3 Multiobjective Optimal Control Method for Generation Personal Sound Zone

3.1 Determine the Objective Function.

Considering the three optimization objectives AC, SD, and AE, higher AC often requires higher speaker power and poorer reconstruction accuracy. At some frequencies, better reconstruction accuracy also means that the speaker power consumption needs to be increased; thus, these three evaluation indicators are contradictory. The desired effect is to make the AC as large as possible and the SD and array power consumption AE as small as possible.

So the objective function F(S) is
(12)
where F1=1/fAC, F2=fSD, and F3=fAE.
Define the independent variable as the weight vector of the speaker S is
(13)
In practice, the power of the speaker is limited, which means that S is not infinite. Therefore, based on the power limitation of the speaker, S is constrained under the following condition:
(14)
where η is the constraint weight, 0<η1, and V is the rated input voltages of individual speakers.
The multiobjective optimization model can be established as follows:
(15)

3.2 Multiobjective Optimal Control Algorithm PSZ-MOGA.

As shown in Eq. (5), under the condition that the acoustic transfer function does not change (i.e., the position of the loudspeaker and the external interference do not change), the speakers' weight vectors play a decisive role in the reproduction effect of the sound field. PSZ-MOGA is proposed to achieve a better effect of PSZ. Under the constraint of the speaker's power, the algorithm uses continuous iteration to narrow the range and then finds the speakers' weight vectors that can make AC, SD, and AE reach Pareto optimality. The flow of the PSZ-MOGA is shown in Table 1.

Table 1

The flow of PSZ-MOGA

Inputs: the objective function F(S), the constraints on loudspeaker power |S|2, the target sound field Pdes, the number of loudspeakers L, the number of microphones in the bright and dark zones MB and MD, and the acoustic transfer functions in the bright and dark zones GB and GD.
Initialization: Initial population size nPop, maximum number of iterations Gen, crossover probability Pc, variation probability Pm.
1: Independent variable S=[S1,S2,S3,,SL]T
2: for i = 1 to maximum frequency
3: Finding the optimal speaker weight vector S using improved NSGA-II (see Sec. 3.2.1 for specific processes).
4: Calculate the Pareto optimal solution curve using the best S.
5: Selection of the optimal solution using the Topsis decision-making method (see Sec. 3.2.2 for the process).
6: end
Output: The values of AC, SD, and AE in the wide band
Inputs: the objective function F(S), the constraints on loudspeaker power |S|2, the target sound field Pdes, the number of loudspeakers L, the number of microphones in the bright and dark zones MB and MD, and the acoustic transfer functions in the bright and dark zones GB and GD.
Initialization: Initial population size nPop, maximum number of iterations Gen, crossover probability Pc, variation probability Pm.
1: Independent variable S=[S1,S2,S3,,SL]T
2: for i = 1 to maximum frequency
3: Finding the optimal speaker weight vector S using improved NSGA-II (see Sec. 3.2.1 for specific processes).
4: Calculate the Pareto optimal solution curve using the best S.
5: Selection of the optimal solution using the Topsis decision-making method (see Sec. 3.2.2 for the process).
6: end
Output: The values of AC, SD, and AE in the wide band

3.2.1 Improved Non-Dominated Genetic Algorithm II.

Non-dominated Genetic Algorithm II (NSGA-II) [21] has better performance in solving multiobjective optimization problems, but it also suffers from slow convergence and may fall into local optimal solutions.

In this article, we make the following three improvements to the NSGA, and the flow of the improved NSGA-II is shown in Fig. 2.

  • (1) Further optimization of the generated new populations

Fig. 2
Flowchart of the improved NSGA-II
Fig. 2
Flowchart of the improved NSGA-II
Close modal

The PSZ-MOGA aims to achieve a balance between AC, SD, and AE under the constraints of the loudspeaker array power, so the PSZ-MOGA should at least obtain an SD superior to the ACC algorithm, an AC superior to the PM algorithm, and an AE superior to both. As shown in Fig. 4, the average value of SD and AE in the calculation results of the ACC and PM algorithm in this article is about 0.4 and 0, respectively, so AC should be more than 11 dB [22], and SD and AE need to be less than 0.4 and 0, respectively. This limits the search space of the algorithm, improves its accuracy and convergence speed, and reduces the number of iterations.

  • (2) Dynamically update parameters that cannot be changed once selected

    • The optimal part of the merged population needs to be selected each time to generate a new population. However, the number of selectable populations will decrease since the populations that do not meet the AC restriction are removed. Therefore, the crossover and mutation probabilities are dynamically updated to ensure the diversity of the populations.

    • The solutions at the end of each sort are not necessarily uniformly distributed, and dynamic updating of crossover and mutation probabilities improves population coverage and avoids the local optimization problem.

The crossover probability Pgc and mutation probability Pgm are improved to address the above problem as follows:
(16)
where n is the number of current iterations, Gen is the maximum number of iterations, nPop is the population size, λ is the step of change in crossover probability, and γ is the step of change in mutation probability.
  • (3) The formula of the traditional crowding distance is
    (17)
Equation (17) only considers the distance between neighboring individuals, which reduces the diversity of the population when the solutions are unevenly distributed and may also lead to premature convergence of the algorithm. Therefore, it is proposed to use the average absolute difference distance to calculate the crowding distance to improve the distribution density of the solutions, as shown in Eq. (18):
(18)
where k is the number of objective functions, i and w are the ordinal numbers of the solutions, and R is the number of solutions in the Pareto front.

3.2.2 Topsis Decision-Making Methodology.

The Topsis algorithm is an algorithm for multiattribute decision analysis, which determines the best choice by comparing the combined scores of different alternatives under a set of evaluation indicators in Ref. [23]. The higher the score, the more superior the solution. The final score for each solution is calculated by Eq. (19):
(19)
where i is the number of solutions, j is the number of objective functions, 1jk, wj is the expected weight of the j, and eij is the normalized value of the i objective of the j solution.

The Topsis algorithm evaluates the advantages and disadvantages of the alternatives by calculating the distance of each alternative from the ideal and least ideal solutions and ultimately determines the best choice.

4 Simulation and Result Analysis

4.1 Simulation Parameter Setting.

Simulations are carried out to verify the effectiveness of the proposed PSZ-MOGA. The PSZ-MOGA is applicable to any arrangement of bright and dark zones and speakers in the car cabin. The speaker arrangement scheme used in this simulation is shown in Fig. 3. Taking the common symmetrical distribution of 10 speakers as an example, with the main driver's seat as the bright zone and the right side of the rear row as the dark zone, 18 microphones are arranged in the bright and dark zones, totaling 36 control points, to construct the PSZ system. The No. 5 speaker is selected as the reference speaker to generate the reference sound field for calculating the AE. The ideal sound field is set to have zero sound pressure in dark zones. Set the frequency response range to 60–4000 Hz and the sampling frequency to 10,240 Hz, and set the constraints, |S|2η×3×|Sr|2, 0<η1, for the loudspeaker array with the single speaker signal not exceeding 3 V. The initial number of populations is set to 80, the maximum number of iterations is 100, and the initial crossover and mutation probabilities are 0.8 and 0.5, respectively. The target weights of AC, SD, and AE in the Topsis algorithm are set to 0.5, 0.3, and 0.2, respectively.

Fig. 3
The speakers’ arrangement scheme used for the simulation: (a) speakers on the left side door, (b) speakers on the right side door, (c) speakers at the front and top of the car, and (d) top view of all the speakers
Fig. 3
The speakers’ arrangement scheme used for the simulation: (a) speakers on the left side door, (b) speakers on the right side door, (c) speakers at the front and top of the car, and (d) top view of all the speakers
Close modal

4.2 Simulation Results Analysis.

The PSZ-MOGA method proposed is compared with the ACC method and the PM method, using the same number and position of speakers. Using 10 speakers, the AC, SD, and AE values are shown in Fig. 4. From Fig. 4, it can be seen that the ACC method achieves the best AC, but this comes at the cost of significant SD. The PM method exhibits less SD and better sound quality, but its acoustic contrast effect is less pronounced. From Refs. [1,2], it can be seen that S for the ACC and PM methods is selected with the determination of the transfer function, and therefore, no constraints can be imposed on the AE, so they are not ideal for AE. In contrast, the PSZ-MOGA method can find the Pareto-optimal curves of AC, SD, and AE in fewer than 100 iterations. The optimal solution is determined by the preset weights to improve AC while keeping SD as low as possible under the speaker power constraints.

Fig. 4
Control effect indicators of PSZ-MOGA and ACC, PM methods: (a) AC, (b) AE, and (c) SD
Fig. 4
Control effect indicators of PSZ-MOGA and ACC, PM methods: (a) AC, (b) AE, and (c) SD
Close modal

Two multiobjective optimization methods, WSM and ECM, have also been proposed in Ref. [13], and since both methods can achieve the same optimization results, the WSM method is used to compare with the PSZ-MOGA method. The results are shown in Fig. 5. Although the WSM method has better SD around 1500–2500 Hz, this comes at the expense of AE, and PSZ-MOGA gives higher AC at a similar SD and AE.

Fig. 5
Control effect indicators of PSZ-MOGA and WSM methods: (a) AC, (b) AE, and (c) SD
Fig. 5
Control effect indicators of PSZ-MOGA and WSM methods: (a) AC, (b) AE, and (c) SD
Close modal

According to Ref. [13], the WSM method can only find the optimal solution on the convex Pareto front, but the Pareto front of the multiobjective optimization problem of AC, SD, and AE is not always convex in all frequency bands, so there is a situation that the optimal solution cannot be found. Moreover, the selection of weights only reflects the trend of the optimization objective, and the relationship with the AC, SD, and AE values is not clear in advance. The performance of the ECM method relies on the setting of ε constraint, and unreasonable ε constraint values can lead to a less-than-optimal solution or an infeasible optimization problem, whereas the constraint values of AC, SD, and AE for arbitrary environments are unknown beforehand, and repeated attempts to determine the constraint values are required. The selection of weights for the PSZ-MOGA is directly related to the value of the optimization objective, and a more comprehensive Pareto optimal solution can be found without constraining the optimization objective. And when the number of optimization objectives increases, the PSZ-MOGA method is still applicable. However, PSZ-MOGA increases the computation time compared to these two methods.

4.3 Sound Field Control Effects Based on Mixed Speaker Array.

In the process of sound field reproduction in a car cabin, different speaker arrays yield varying reproduction effects. To achieve more effective multiobjective optimization control of PSZ, a sound field control method based on a mixed speaker array is explored. The accuracy of the reproduced sound field in different zones by separating the target frequency bands and controlling them independently with different speaker arrays is optimized.

In this study, the frequency range of interest for sound field reproduction is divided into three frequency bands, using low, medium, and high speakers with distinct frequency response ranges for sound field reproduction. The frequency bands were divided based on the crossover method proposed in Ref. [15]:
(20)
(21)

As shown in Fig. 6, fLi and fHi(i=1,2,3), fiL and fiH(i=1,2,3), and fL and fH represent the upper and lower bounds of the target optimization frequency range of different frequency bands, working frequency range of different speakers, and needed optimization frequency range, respectively. The different frequency bands should not interfere with each other to avoid affecting the sound field reproduction of other bands. In this article, the frequency bands are defined as 60–100 Hz for low frequency, 100–500 Hz for middle frequency, and 500–4000 Hz for high frequency.

Fig. 6
Frequency band division of three kinds of loudspeakers
Fig. 6
Frequency band division of three kinds of loudspeakers
Close modal

After completing the frequency band division, a total of 2 low-frequency speakers, 4 mid-frequency speakers, and 4 high-frequency speakers (10 speakers in total) were used for sound field reproduction.

Inspired by Ref. [15], the gain Gl2 of the transfer function of the loudspeaker indicates the controllability of the lth candidate loudspeaker position to the loudspeaker array, the larger Gl2 indicates that the loudspeaker in this position has a greater influence on the sound field reproduction should be retained, and the controllability of the wide band loudspeaker candidate position is the sum of the single frequency controllability of the position.

The controllability Cl is calculated as shown in Eq. (22):
(22)

The controllability of the 10 candidate positions in the low-, mid-, and high-frequency bands was calculated separately, and the positions of the 2 low-frequency loudspeakers, 4 mid-frequency loudspeakers, and 4 high-frequency loudspeakers were determined according to the magnitude of Cl. The results are shown in Table 2, and the arrangement of the mixed speaker array is shown in Fig. 7.

Fig. 7
The arrangement of the mixed speaker array
Fig. 7
The arrangement of the mixed speaker array
Close modal
Table 2

Controllability of different loudspeaker candidate positions

FrequencyC1C2C3C4C5C6C7C8C9C10
60–10016.358.2519.865.6832.7628.675.253.589.1135.58
100–500178.38302.54358.6156.0385.61160.41296.57288.9881.7890.21
500–4000609.16526.03374.91662.07506.29735.32563.08218.8710.63533.72
FrequencyC1C2C3C4C5C6C7C8C9C10
60–10016.358.2519.865.6832.7628.675.253.589.1135.58
100–500178.38302.54358.6156.0385.61160.41296.57288.9881.7890.21
500–4000609.16526.03374.91662.07506.29735.32563.08218.8710.63533.72

According to this mixed speaker array arrangement, the PSZ-MOGA is used for multiobjective optimization, and the optimized reproduction results are shown in Fig. 8. Inspired by Refs. [20,24], the AC is reduced in the low- and mid-frequency parts due to the limitation of the number of loudspeakers, while at high frequencies the directivity of the loudspeakers is increased, so that the AC is not affected significantly. The use of mixed loudspeaker arrays offers significant advantages in the high-frequency section, with only four speakers, and gives control that approximates the PSZ of a full-range loudspeaker array and has a better SD overall.

Fig. 8
Comparison of control effects between full-range and mixed speaker arrays: (a) AC, (b) AE, and (c) SD
Fig. 8
Comparison of control effects between full-range and mixed speaker arrays: (a) AC, (b) AE, and (c) SD
Close modal

5 Experimental

The experiments were conducted in a car within a semianechoic chamber. The main driver's headrest was designated as the acoustically bright area, and the rear right seat headrest as the acoustically dark area. The mixed speaker array shown in Fig. 7 is used for sound field reproduction. Eighteen microphones (PCB-378B02) were placed in each area to measure the impulse response. Each measurement lasted 10 s. The input and output signals were processed using the BEHRINGER X32 RACK, and the LMS SCADAS system was used to collect the time-domain sound pressure signals. Frequency-domain sound pressure signals were obtained via Fourier transform to calculate the AC between the bright and dark zones. The control frequency range was 60–4000 Hz, the sampling frequency was 8192 Hz, and the frequency resolution was set to 1 Hz. The process of experimental validation is shown in Fig. 9.

Fig. 9
The process of experimental validation
Fig. 9
The process of experimental validation
Close modal

The proposed method is compared with the WSM method in Ref. [13], and the parameters are set the same as in the simulation. Figure 10 shows that the PSZ-MOGA method achieves better SD and AE, and AC is almost the same.

Fig. 10
Experimental results of PSZ-MOGA method and WSM method for AC, SD, and AE at different frequencies: (a) AC, (b) AE, and (c) SD
Fig. 10
Experimental results of PSZ-MOGA method and WSM method for AC, SD, and AE at different frequencies: (a) AC, (b) AE, and (c) SD
Close modal

6 Conclusion

To achieve better sound field reproduction while generating interior vehicle PSZ, the balance of AC, SD, and AE is studied. The PSZ-MOGA multiobjective optimization method outperforms existing algorithms and effectively addresses the balance problem among AC, SD, and AE. The performance of the reproduced sound field can be effectively controlled by simply increasing and decreasing the weights, which provides better flexibility. Using a mixed speaker array optimizes the frequency response of each band, reduces intermodulation distortion, improves overall sound quality, and provides better sound field reproduction performance compared to traditional speaker array. The proposed method is not restricted by the number of speakers, arrangement, and the configuration of bright and dark zones.

Funding Data

  • The Project of National Natural Science Foundation of China (No. 51905331).

  • Young Scientific Research Team Cultivation Program of SUES (QNTD202112).

  • Technical Service Platform for Vibration and Noise Evaluation and Control of New Energy Vehicles (No. 18DZ2295900).

Conflict of Interest

There are no conflicts of interest.

Data Availability Statement

The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.

References

1.
Choi
,
J. W.
, and
Kim
,
Y. H.
,
2002
, “
Generation of an Acoustically Bright Zone With an Illuminated Region Using Multiple Sources
,”
J. Acoust. Soc. Am.
,
111
(
4
), pp.
1695
1700
.
2.
Poletti
,
M.
,
2008
, “
An Investigation of 2-d Multizone Surround Sound Systems
,”
Audio Engineering Society Convention 125
,
San Francisco, CA
.
3.
Chang
,
J. H.
, and
Jacobsen
,
F.
,
2012
, “
Sound Field Control With a Circular Double-Layer Array of Loudspeakers
,”
J. Acoust. Soc. Am.
,
131
(
6
), pp.
4518
4525
.
4.
Shi
,
L.
,
Lee
,
T.
,
Zhang
,
L.
,
Nielsen
,
J. K.
, and
Christensen
,
M. G.
,
2020
, “
A Fast Reduced-Rank Sound Zone Control Algorithm Using the Conjugate Gradient Method
,”
Proceedings of the 2020 IEEE ICASSP
,
Barcelona, Spain
,
May 4–8
.
5.
Shi
,
L.
,
Lee
,
T.
,
Zhang
,
L.
,
Nielsen
,
J. K.
, and
Christensen
,
M. G.
,
2021
, “
Generation of Personal Sound Zones With Physical Meaningful Constraints and Conjugate Gradient Method
,”
IEEE/ACM Trans. Audio, Speech, Lang. Process.
,
29
(
1
), pp.
823
837
.
6.
Brunnström
,
J.
,
Koyama
,
S.
, and
Moonen
,
M.
,
2022
, “
Variable Span Trade-off Filter for Sound Zone Control With Kernel Interpolation Weighting
,”
ICASSP 2022-2022 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP)
,
Singapore
,
May 23–27
, pp.
1071
1075
.
7.
Zhang
,
W.
,
Zhang
,
J.
,
Abhayapala
,
T. D.
, and
Zhang
,
L.
,
2018
, “
2.5 D Multizone Reproduction Using Weighted Mode Matching
,”
2018 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP)
,
Calgary, AB, Canada
,
Apr. 15–20
, pp.
476
480
.
8.
Koyama
,
S.
,
Amakasu
,
T.
,
Ueno
,
N.
, and
Saruwatari
,
H.
,
2021
, “
Amplitude Matching: Majorization–Minimization Algorithm for Sound Field Control Only With Amplitude Constraint
,”
ICASSP 2021-2021 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP)
,
Toronto, ON, Canada
,
June 6–11
, pp.
411
415
.
9.
Coleman
,
P.
,
Jackson
,
P. J.
,
Olik
,
M.
,
Møller
,
M.
,
Olsen
,
M.
, and
Abildgaard Pedersen
,
J.
,
2014
, “
Acoustic Contrast, Planarity and Robustness of Sound Zone Methods Using a Circular Loudspeaker Array
,”
J. Acoust. Soc. Am.
,
135
(
4
), pp.
1929
1940
.
10.
Lee
,
T.
,
Shi
,
L.
,
Nielsen
,
J. K.
, and
Christensen
,
M. G.
,
2020
, “
Fast Generation of Sound Zones Using Variable Span Trade-off Filters in the DFT-Domain
,”
IEEE/ACM Trans. Audio, Speech, Lang. Process.
,
29
, pp.
363
378
.
11.
Cai
,
Y.
,
Wu
,
M.
, and
Yang
,
J.
,
2014
, “
Sound Reproduction in Personal Audio Systems Using the Least-Squares Approach With Acoustic Contrast Control Constraint
,”
J. Acoust. Soc. Am.
,
135
(
2
), pp.
734
741
.
12.
Hu
,
M.
,
Zou
,
H.
,
Lu
,
J.
, and
Christensen
,
M. G.
,
2022
, “
Maximizing the Acoustic Contrast With Constrained Reconstruction Error Under a Generalized Pressure Matching Framework in Sound Zone Control
,”
J. Acoust. Soc. Am.
,
151
(
4
), pp.
2751
2759
.
13.
Cheng
,
J.
,
Gao
,
H.
,
Pan
,
K.
,
Feng
,
X.
, and
Shen
,
Y.
,
2023
, “
Multi-Objective Optimization for Generation of Personal Sound Zone
,”
J. Acoust. Soc. Am.
,
153
(
1
), pp.
149
158
.
14.
Cheer
,
J.
, and
Elliott
,
S.
,
2013
, “
Design and Implementation of a Personal Audio System in a Car Cabin
,”
Proceedings of Meetings on Acoustics (Vol. 19, No. 1)
,
Montreal, Canada
,
June 2–7
.
15.
Peng
,
B.
,
Zheng
,
S.
,
Liao
,
X.
, and
Lian
,
X.
,
2018
, “
On the Optimization of a Mixed Speaker Array in an Enclosed Space Using the Virtual-Speaker Weighting Method
,”
Mech. Syst. Signal Process.
,
102
(
2
), pp.
214
229
.
16.
Zhao
,
S.
, and
Burnett
,
I. S.
,
2022
, “
Evolutionary Array Optimization for Multizone Sound Field Reproduction
,”
J. Acoust. Soc. Am.
,
151
(
4
), pp.
2791
2801
.
17.
Gao
,
H.
,
Feng
,
X.
, and
Shen
,
Y.
,
2022
, “
Weighted Loudspeaker Placement Method for Sound Field Reproduction
,”
IEEE/ACM Trans. Audio, Speech, Lang. Process.
,
30
(
3
), pp.
1263
1276
.
18.
Stuart
,
J. R.
,
Hollinshead
,
R.
, and
Capp
,
M.
,
2019
, “
Is High-Frequency Intermodulation Distortion a Significant Factor in High-Resolution Audio?
,”
J. Audio Eng. Soc.
,
67
(
5
), pp.
310
318
.
19.
Thipmaungprom
,
Y.
,
2023
, “
Intermodulation Distortions From a Chain of Coupled Nonlinear Oscillators Near a Hopf Bifurcation
,” Dissertation,
Chulalongkorn University
,
Bangkok, Thailand
, pp.
22
41
.
20.
Zhu
,
M.
, and
Zhao
,
S.
,
2021
, “
An Iterative Approach to Optimize Loudspeaker Placement for Multi-Zone Sound Field Reproduction
,”
J. Acoust. Soc. Am.
,
149
(
5
), pp.
3462
3468
.
21.
Murugan
,
P.
,
Kannan
,
S.
, and
Baskar
,
S.
,
2009
, “
NSGA-II Algorithm for Multi-Objective Generation Expansion Planning Problem
,”
Electr. Power Syst. Res.
,
79
(
4
), pp.
622
628
.
22.
Druyvesteyn
,
W. F.
, and
Garas
,
J.
,
1997
, “
Personal Sound
,”
Audio Eng. Soc. Convent.
,
45(9)
, pp.
685
701
.
23.
Sałabun
,
W.
,
Watróbski
,
J.
, and
Shekhovtsov
,
A.
,
2020
, “
Are MCDA Methods Benchmarkable? A Comparative Study of TOPSIS, VIKOR, COPRAS, and PROMETHEE II Methods
,”
Symmetry
,
12
(
9
), p.
1549
.
24.
Lopez
,
J. J.
,
Gutierrez-Parera
,
P.
, and
Savioja
,
L.
,
2020
, “
Effects and Applications of Spatial Acuity in Advanced Spatial Audio Reproduction Systems With Loudspeakers
,”
Appl. Acoust.
,
161
(
4
), p.
107179
.