#2018 The Acoustical Society of Japan
Acoust. Sci. & Tech. 39, 3 (2018)
PAPER
Sidelobe suppression by desired directivity pattern optimization for a small circular loudspeaker array Koya Sato and Yoichi Haneda Graduate School of Informatics and Engineering, The University of Electro-Communications, 1–5–1 Chofugaoka, Chofu, 182–8585 Japan (Received 18 August 2017, Accepted for publication 25 December 2017) Abstract: Directivity control using a loudspeaker array is widely studied for various applications. Suppressing sidelobe levels is important for applications such as personal audio systems. In this paper, we propose a filter design method using a window function shape as the desired directivity pattern to reduce the sidelobe levels. The proposed method consists of three steps. The first step defines a cost function with a criterion for the directivity pattern. Next, filter coefficients for each loudspeaker are calculated and stored while changing the window function shape of the desired directivity pattern. Finally, we determine the optimum filter coefficients having the best performance of the cost function by using a full-search algorithm at each frequency. We conducted directivity experiments with a real six-element circular loudspeaker array having a radius of 0.055 m and evaluated its directivity to confirm the performance of the proposed method. The results, which were compared with those obtained from a conventional method, showed that the maximum sidelobe level improved by about 2 dB, although the beam was wide. We verified that using the window function shape as the desired directivity pattern is more effective than using the conventional method for sidelobe suppression. Keywords: Circular loudspeaker array, Directivity control, Sidelobe suppression, Window function, Least-squares method PACS number: 43.60.Fg
1.
[doi:10.1250/ast.39.243]
INTRODUCTION
Directivity and sound field control for loudspeaker arrays have been developed over a long time to be applied, for instance, in personal audio, public address, and surround sound systems [1–10]. In addition, directivity control, also known as beamforming in the field of microphone array processing, is widely studied for speech enhancement, sound source separation, and direction of arrival estimation [11–15]. Filter design methods that allow filter coefficients to be obtained are important to control the directivity or sound field in these applications. The filter coefficients are mainly obtained in two domains. First, the pressure-matching method (PMM) based on the least-squares method [16] in the frequency domain is widely used for sound field control [17–20]. The PMM is an extension of inverse filtering [21]. Second, the wavenumber domain has recently attracted attention because the filter coefficients, the radiation patterns of a loudspeaker array, and the corresponding gains can be
e-mail:
[email protected]
derived analytically [22–29]. However, the wavenumber domain can currently be used for specific array configurations, namely, linear, cylindrical, circular, and spherical configurations. Considering personal audio systems, directivity control aims to reduce sound leakage outside the listening area to guarantee privacy protection [30]. That is, it is important to suppress sidelobe levels, which appear in directions besides the look direction [31–33]. To attain this, we focus on a filter design method in the wavenumber domain using either a rectangular or Hanning window [34] as the distribution of the desired sound pressure [27,28]. Although the window function shapes are fixed in this method, we can modify the shape according to the desired directivity pattern in our proposed method. Moreover, we define a cost function with a criterion for the directivity pattern to determine the optimum filter coefficients that suppress the sidelobe levels. Hence, we propose a filter design method with a full-search algorithm based on the cost function. In this paper, we address directivity control with sidelobe suppression by using a small loudspeaker array with a small number of loudspeakers because we
243
Acoust. Sci. & Tech. 39, 3 (2018)
pð!Þ ¼ Gð!Þwð!Þ; with Gð!Þ defined as 0 G11 ð!Þ B G21 ð!Þ B Gð!Þ ¼ B .. B @ . GM1 ð!Þ
Fig. 1 Distribution of the circular array and control points.
assume that a personal audio system must be small for practical use. This paper is organized as follows. Section 2 describes the PMM based on the least-squares method as a conventional method. Section 3 presents the window function shape that corresponds to the desired directivity pattern in the PMM. Section 4 explains the three cost functions used to determine the filter coefficients and gives details of our proposed method using a full-search algorithm to suppress the sidelobe levels. Experiments on directivity patterns and their evaluation using a real six-element circular loudspeaker array are described in Sects. 5 and 6. Notes on the proposed method are provided in Sect. 7. Finally, conclusions are given in Sect. 8.
2.
PRESSURE-MATCHING METHOD (PMM)
The PMM is a typical filter design method in the frequency domain for sound field control [17–20]. In this section, we introduce the use of the PMM in directivity control. We consider the spherical coordinate system shown in Fig. 1 throughout this paper. A circular loudspeaker array consisting of L elements is placed on the x–y plane. M control points with coordinates ðr; m ; m Þ are set on the spherical surface, where r is the distance from the origin to the control points, and m and m are the elevation and horizontal angles of the mth control point, respectively. Note that the ranges of m and m are defined as 90 m 90 and 180 m 180 , respectively. When the circular array is driven by the filter coefficients in the vector wð!Þ ¼ ½w1 ð!Þ; w2 ð!Þ; . . . ; wL ð!ÞT , the observed sound pressures at each control point pð!Þ ¼ ½pð1 ; 1 ; !Þ; pð2 ; 2 ; !Þ; . . . ; pðM ; M ; !ÞT are determined by
244
ð1Þ 1
G12 ð!Þ
...
G1L ð!Þ
G22 ð!Þ .. .
... .. .
G2L ð!Þ C C C; .. C A .
GM2 ð!Þ
. . . GML ð!Þ
ð2Þ
where Gð!Þ is the transfer function matrix, Gml ð!Þ is the transfer function between the lth loudspeaker and the mth control point, and ! ¼ 2 f is the angular frequency corresponding to frequency f . Let dð!Þ ¼ ½dð1 ; 1 ; !Þ; dð2 ; 2 ; !Þ; . . . ; dðM ; M ; !ÞT denote a desired directivity pattern defined at each control point. From the PMM, we obtain the filter coefficients that minimize the sum of the square errors between the observed and desired directivity patterns [17,19]. This is expressed as the following minimization problem: ð3Þ min kdð!Þ Gð!Þwð!Þk22 ; wð!Þ P where kxk22 is defined as i jxi j2 . Then, the optimum filter coefficients are given by wð!Þ ¼
Gð!ÞH dð!Þ ; Gð!ÞH Gð!Þ
ð4Þ
where the superscript ðÞH denotes the complex conjugate transposition. Moreover, Eq. (4) can include a regularization parameter [16] to obtain the stable filter gains or white noise gain [19,20,35].
3. DESIRED DIRECTIVITY PATTERN WITH WINDOW FUNCTION SHAPE When using the conventional PMM for sound field control, the desired sound pressure is usually set to 1 in the look direction and 0 at the other control points [18,20]. In research to generate acoustically bright and dark zones for personal sound systems [36,37], filter design methods in which the desired sound pressure is modeled by a rectangular or Hanning window were recently proposed [27,28]. These methods achieved a high contrast between the bright and dark zones, and were analytically derived in the wavenumber domain. If we use the conventional PMM to control directivity, sidelobes appear because this method cannot completely reproduce the desired directivity pattern. In order to mitigate this problem, we apply a window function shape as the desired directivity pattern in consideration of the results in [27,28]. Although there are several window functions, we use the Hanning and rectangular windows in this paper. We describe the relation between the control points shown in Fig. 1 and the window function shape. When
K. SATO and Y. HANEDA: SIDELOBE SUPPRESSION FOR SMALL CIRCULAR ARRAY
where w is the window width at the horizontal angle, whose range is from 0 to 360 . The desired directivity pattern should be set at control points on a spherical surface because sound waves are radiated to the three-dimensional (3D) space in a real environment. The desired 3D directivity pattern is obtained by applying the Hanning window to the control points arranged in the vertical plane in an analogous way to how it was obtained for the horizontal angle, i.e., 8 ðm ¼ 0 ; w ¼ 0 Þ, > : 0 (otherwise),
-30
The ranges of w and w are the same as those when using the Hanning window. Figure 2 shows examples of the desired directivity patterns in the horizontal plane using the conventional method, Hanning, and rectangular window shapes. In Figs. 2(b) and 2(c), the solid lines indicate that the Hanning and rectangular windows are set at the control points arranged at horizontal angles from 90 to 90 , whereas the dashed lines indicate that the window functions are set at the control points from 30 to 30 , respectively. The corresponding window widths, w , are 180 and 60 . The main lobe of the conventional desired directivity pattern is narrow, as shown in Fig. 2(a), whereas that of the desired directivity patterns using the window functions is wide, as shown in Figs. 2(b) and 2(c). Hence, we can form a variety
Sound pressure
1
30
0.8 0.6
-60
0.8
60
0.4
0.6
0.2
-90
0.4
90
0.2 0 -180
-90
0
φ [ °]
90
180
120
-120 -150
150
180 -180
(a) Conventional desired directivity pattern. -30 1
±90° Hann. ±30° Hann.
0.8
0
1
30
0.8 0.6
-60
60
0.4
0.6
0.2
90
-90
0.4
±90° Hann. ±30° Hann.
0.2 0 -180
-90
0
φ [ °]
90
180
120
-120 -150
150
180 -180
(b) Desired directivity patterns using Hanning windows.
ð6Þ
-30 1
Sound pressure
where w is the window width at the elevation angle, whose range is from 0 to 180 . Therefore, we obtain the distribution of the desired 3D directivity pattern, dðm ; m ; !Þ, by multiplying the elements calculated from Eqs. (5) and (6). When using a rectangular window as the desired directivity pattern, the definitions are given as follows: 8 1 ðm ¼ 0 ; w ¼ 0 Þ, > > > < w w ð7Þ dðm ; w ; !Þ ¼ 1 m , > 2 2 > > : 0 (otherwise), 8 1 ðm ¼ 0 ; w ¼ 0 Þ, > > > < w w , m dðm ; w ; !Þ ¼ 1 ð8Þ > 2 2 > > : 0 (otherwise).
0
1
Sound pressure
using a Hanning window, the desired directivity pattern at the horizontal angle m is defined as follows: 8 ðm ¼ 0 ; w ¼ 0 Þ, > : 0 (otherwise), ð5Þ
±90° rect. ±30° rect.
0.8
0
1 0.8
30
0.6
-60
60
0.4
0.6
0.2
-90
0.4 0.2 0 -180
-90
0
φ [ °]
90
180
90
±90° rect. ±30° rect.
120
-120 -150
180 -180
150
(c) Desired directivity patterns using rectangular windows. Fig. 2 Sound pressures in the horizontal plane. The left figures show the sound pressure as a function of the angle. The right figures show the corresponding polar graphs.
of desired directivity patterns, with narrow to wide main lobes, by changing the window widths.
4.
FILTER DESIGN METHOD WITH FULL SEARCH
To obtain the filter coefficients, wð!Þ, that can suppress the sidelobe levels from Eq. (4), we must determine the window widths at the elevation and horizontal angles for the desired directivity pattern, dð!Þ. Modifying the window widths causes a variation of the directivity pattern and filter coefficients. Therefore, the window widths that result in the desired directivity pattern to suppress the sidelobe levels are unknown. For this reason, we need to define the best directivity pattern. Consequently, we prepare a cost function that defines a criterion for this directivity pattern. In addition, we propose a full-search algorithm method for
245
Acoust. Sci. & Tech. 39, 3 (2018)
the filter coefficients to satisfy the cost function. The search is performed over all the possible window widths. Although there are many ways of describing cost functions, we experimentally prepare three cost functions in our study. These cost functions are described as X X MS jpðm ; m ; !Þj2 T T X X J1 ðw ; w ; !Þ ¼ ; ð9Þ MT jpðm ; m ; !Þj2 S
J2 ðw ; w ; !Þ ¼
S
1 XX jpðm ; m ; !Þj2 ; MS S S
ð10Þ
J3 ðw ; w ; !Þ ¼ max jpðm ; m ; !Þj; ðm ; m 2 S ; S Þ; ð11Þ where T and T are the angles specifying the target area, S and S denote the sidelobe area whose levels we aim to suppress, and MT and MS are the numbers of control points in the target and sidelobe areas, respectively. In the proposed method, we consider the maximization of Eq. (9), that is, the maximization of the energy ratio between the target and sidelobe areas. On the other hand, we minimize Eqs. (10) and (11). These operations aim to minimize the energy and the maximum sidelobe in the sidelobe area, respectively. The proposed method is summarized in Algorithm 1 and consists of three main steps. In the first step, the window function is selected in line 1 of Algorithm 1. We use either a Hanning or rectangular window. Next, the cost function is determined from Eq. (9), Eq. (10), or Eq. (11) in line 2. and are the corresponding angle increments in the loop of the algorithm. In the second step, which comprises line 4 to line 14 of Algorithm 1, the filter coefficients and cost values are calculated while changing the window function shape for the desired directivity pattern. Finally, we retrieve the optimum filter coefficients from the cost value that satisfies the corresponding criterion in line 15. Line 4 to line 15 of this algorithm are executed for each angular frequency !.
5.
1: Choose window function: Hanning or rectangular window. 2: Determine cost function J from Eq. (9), Eq. (10), or Eq. (11). 3: Set increments of window widths, and . 4: Initialize window widths w ¼ w ¼ 0 . 5: while w 180 do 6: while w 360 do 7: Create dð!Þ using Eqs. (5) and (6), or Eqs. (7) and (8). 8: Calculate wðw ; w ; !Þ using Eq. (4). 9: Calculate pð!Þ using Eq. (1). 10: Calculate the cost value based on Jðw ; w ; !Þ. 11: w w þ 12: end while w þ 13: w 14: end while 15: Select optimum wðw ; w ; !Þ when Jðw ; w ; !Þ satisfies the condition.
DIRECTIVITY EXPERIMENTS
5.1. Measurements and Directivity Experiments To confirm the directivity of our proposed method using a real loudspeaker array, we measured the transfer functions between the loudspeaker array and the control points. Figure 3 shows the circular loudspeaker array and the measuring device. The circular array was composed of six loudspeaker units. The radius and height of the array were 0.055 m and 0.25 m, respectively. For the measuring device, the microphones were arranged at 5 intervals in the elevation angle, and the table at the center was rotated at 10 intervals to vary the horizontal angle. Note that the transfer functions at ðm ; m Þ ¼ ð90 ; m Þ, which was a point located below the array, could not be measured owing to the presence of the rotation shaft. The transfer functions
246
Algorithm 1 for the proposed method.
Fig. 3 Six-element circular loudspeaker array and measuring device.
between the loudspeaker unit and the control points were measured in an anechoic chamber. In total, we obtained the transfer functions at 1,296 points on the spherical surface with a radius of 0.5 m. After measuring the transfer functions, we conducted directivity experiments to evaluate the proposed method under the conditions listed in Table 1. In addition, we conducted the experiments using the conventional PMM to compare its results with those of the proposed method. In the conventional method, we set the desired directivity pattern as
K. SATO and Y. HANEDA: SIDELOBE SUPPRESSION FOR SMALL CIRCULAR ARRAY Table 1 Experimental conditions.
z- x , 800-3000 Hz
-90
Sound velocity, c [m/s]
340
Sampling frequency [Hz]
48,000
Frequency band [Hz]
800–3,000
Filter length
2,400
Number of elements, L
6
Number of control points, M
1,296
Look direction
ð0 ; 0 Þ
-10
-20
0 30
-30
0
Hann. rect. Conv. 60
0
-10
-60
-20
-90
90 120
-120
60
90
-150
150
180 -180
(a) Cost function is set to J1 .
T ¼ 90 –85 T ¼ 80 –80
Target area
x- y , 800-3000 Hz
Hann. -60 rect. Conv. -30
0
z- x , 800-3000 Hz
-90
Sidelobe area
S ¼ 90 –85 S ¼ 170 –90 S ¼ 90 –180
MT in Eqs. (9), (10), and (13)
612
MS in Eqs. (9), (10), and (13)
684
Increment angles
¼ 10 , ¼ 20
x- y , 800-3000 Hz
Hann. rect. Conv. -30
0
-10 -20
0 30
-30
Hann. rect. Conv. 60
0
-10
-60
-20
-90
90 120
-120
60
90
0
-150
150
180 -180
(b) Cost function is set to J2 .
dðm ; m ; !Þ ¼
1
ðm ; m Þ ¼ ð0 ; 0 Þ,
0
(otherwise).
z- x , 800-3000 Hz
ð12Þ
Owing to the small loudspeaker array, no regularization parameter was applied in the proposed and conventional methods because the pseudoinverse matrices of Gð!Þ were stable.
-90 -10 -20
0 30 90
5.2. Evaluation Indexes We defined three evaluation indexes to assess the directivity performance. One index is the target to sidelobe ratio (TSR) defined as follows: X X MS jpðm ; m Þj2 T T X X TSR ¼ 10 log10 : ð13Þ 2 MT jpð ; Þj m m S
60
-30
0
Hann. rect. Conv. 60
0
-10
-60
-20
-90
90
-120
120 -150
180 -180
150
(c) Cost function is set to J3 . Fig. 4 Results of directivity patterns using the proposed method and the conventional method. The left and right figures show the directivity patterns on the z–x and x–y planes, respectively (Fig. 1 shows the planes).
S
The TSR determines the difference in sound pressure level between the target and sidelobe areas. The other indexes are the maximum sidelobe level (MSL) and the beam width (BW). The MSL is the maximum sound pressure level in the sidelobe area, and the BW is defined as the angular width that is suppressed by not less than 3 dB with respect to the sound pressure level in the look direction.
6.
x- y , 800-3000 Hz
Hann. rect. Conv. -30
0
RESULTS
Figure 4 shows the results of the directivity experiments for the proposed method and for the conventional method. We set the sound pressure level in the look direction to 0 dB as a reference. Sidelobes for the conventional method appeared in several directions, as shown in the right column of Fig. 4, owing to the highly narrow desired directivity pattern. Nevertheless, the sound pressure levels for the conventional method around ð90 ; 0 Þ and
ð90 ; 0 Þ were suppressed slightly more than those for the proposed method as shown in the left column of Fig. 4. The proposed method with J1 and J3 showed similar directivity patterns and suppressed the sidelobe level by about 20 dB in the opposite direction to the look direction, as shown in Figs. 4(a) and 4(c). The proposed method with J2 showed less suppression of the sidelobe level at ð0 ; 180 Þ than the proposed method with the other cost functions from the results on the x–y plane. Overall, the proposed method showed wider main lobes than the conventional method. When we tested the directivity of the proposed method using the rectangular window, the results of the directivity patterns were similar to those using the Hanning window. From the directivity experiments, we concluded that the proposed method with J1 and J3 was superior to that with J2 for sidelobe suppression in the opposite direction. Moreover, there were no significant
247
Acoust. Sci. & Tech. 39, 3 (2018) Table 2 Evaluation results of the conventional and proposed methods.
180 160 140
Hanning window J1
J2
J3
5.7 7:9 20
9.2 9:3 40
8.6 9:2 30
8.6 10:1 30
120
θw [°]
TSR [dB] MSL [dB] BW [ ]
Conventional
100 80 60
J1 J2 J3
40
Rectangular window
TSR [dB] MSL [dB] BW [ ]
Conventional
J1
J2
J3
5.7 7:9 20
9.2 9:9 30
8.7 9:1 30
8.7 10:0 30
20 0
1000
1500
2000
Frequency [Hz]
2500
350
J1 J2 J3
300
248
φw [°]
250 200 150 100 50 0
1000
1500
2000
2500
3000
Frequency [Hz]
Fig. 5 Selected window widths of the Hanning window for each cost function. The top figure shows the window width w for the elevation angle. The bottom figure shows the window width w for the horizontal angle. 0
Relative level [dB]
differences between the two window functions in the results of the directivity patterns. Table 2 shows the evaluation results. The TSR values of the proposed method with the window functions improved by over 3 dB in comparison with the values obtained from the conventional method. These TSR results imply that more energy in the sidelobe area was suppressed. In the MSL, the greatest improvement of the proposed method was 2.2 dB when using J3 with the Hanning window. It seems that the suppression of the sidelobe level is small when referring only to the MSL results. However, the effectiveness of the proposed method is clear when considering the TSR results. When using the cost functions J1 and J3 , we obtained the largest TSR values from J1 and the lowest MSL values from J3 owing to the full-search algorithm. The wide BW was because we defined the desired directivity patterns having wide main lobes using the window function shapes. The selected window widths for the Hanning window in the proposed method are shown in Fig. 5. There were no large variations in the horizontal angle. For J1 and J2 , the corresponding window widths, w , of the elevation angle were extremely narrow around 1,620 Hz, whereas the window width for J3 was highly variable depending on the frequency. The frequency responses in the look direction followed similar trends as shown in Fig. 6. In this figure, we used the conventional method as a reference, and we biased each proposed method by 5 dB for ease of visualization and comparison. The amplitude level of the conventional method was within 5 dB and its frequency response was mostly flat. In contrast, the frequency responses for the methods using J1 and J2 showed notches around 1,620 Hz. Furthermore, the frequency response for the method using J3 was highly variable. Therefore, the frequency response corresponding to J3 suggests that the sound quality is degraded when performing a reproduction using directivity. In fact, when we listened to the directivity reproduction for female speech using the array, the
3000
-10 -20 -30 Conv. J1
-40 -50
J2 J3
1000
1500
2000
Frequency [Hz]
2500
3000
Fig. 6 Frequency response in the look direction when using the Hanning window for each cost function.
reproduced speech signal from the proposed method with J3 had poor sound quality. We were not able to distinguish differences in sound quality among the proposed method with J1 and J2 , and the conventional method. From the above results, we affirmed that the proposed method with J1 was the best for providing sidelobe suppression and suitable sound quality.
K. SATO and Y. HANEDA: SIDELOBE SUPPRESSION FOR SMALL CIRCULAR ARRAY
-30
0
0
x- y , 800-3000 Hz
-30
30
-10
-60
60
-20
-90 -120
90 120 Prototype Theoretical -150 180 150 -180
0
-10
-60
60
-20
-90
90 120
-120 Prototype Theoretical
x- y , 800-3000 Hz
x- y , 800-3000 Hz
0
0
-30
30
-10 -20
-90
90
-120
120
0
-20 -30
1000
1500
2000
Frequency [Hz]
2500
3000
30
-60 -20
-90
60 90 120
-120
Prototype Theoretical
Prototype Theoretical
(c) φw = 160◦
(d) φw = 180◦
Fig. 7 Directivity patterns when w is changed from 120 to 180 using the rectangular window. w is set to 0 in these results.
7.
-10
-50
-10
60
Prototype Theoretical
-40
(b) φw = 140◦
0
0
30
(a) φw = 120◦
-30 -60
0
Relative level [dB]
x- y , 800-3000 Hz
NOTES
From Figs. 5 and 6, we can see that a variation of the desired directivity pattern affects the frequency responses in the look direction. Therefore, a method that sets the same window widths regardless of frequency would be useful to avoid large variations of the frequency response in the look direction. For example, we can fix the window widths to w ¼ 180 and w ¼ 180 . Note that when using the fixed window widths with a rectangular window, the directivity in the look direction is split for window widths at a horizontal angle w above 140 . Figure 7 shows examples of this phenomenon. The transfer function of a rigid circular array has been analytically derived [27]; however, it can only be maintained when using a rigid cylindrical array with infinite length. Because of this, we used theoretical transfer functions assuming the 3D freefield Green’s function [22] in the theoretical simulation. The circular loudspeaker array and control points had the same configurations as those used in Sect. 5. The directivity was split when w was set to above 140 in Fig. 7. Meanwhile, the Hanning window has its maximum value at 1 and is 0 at its edges. Therefore, the desired directivity pattern in the look direction has a gain of 1 and gradually decreases in the other directions. Hence, using a Hanning window can avoid the directivity split in the look direction. We verified that the frequency response without the notches was accomplished by fixing the window widths
Fig. 8 Frequency response in the look direction. The window widths are set to w ¼ 180 and w ¼ 180 by using the Hanning window.
using the Hanning window from Fig. 8. Furthermore, we can reduce the calculation time by using fixed window widths. In contrast, the proposed method in Sect. 4 required a long calculation time owing to the full-search algorithm in wide frequency bands. In Sect. 6, the difference in the effectiveness between the rectangular and Hanning windows was not clear because the array had only six loudspeaker units. By increasing the number of loudspeakers, we may find a difference between these windows in the wavenumber domain. This problem has been described in detail using a linear array in [28]. Since this paper dealt with a small loudspeaker array, an investigation of the difference between these windows for a large circular array is a future work. In this paper, we present successful sidelobe suppression and directivity reproduction using a real six-element circular loudspeaker array. However, when we used a real eight-element circular array, the speeches reproduced with directivity were distorted using the proposed and conventional methods. Owing to the instability of the pseudoinverse matrix of Gð!Þ in Eq. (4), the filter gains [5] were large at low frequencies or the white noise gain [35] was small. To obtain stable inverse matrices, a regularization parameter [16] has to be introduced in Eq. (4) when using a large array. Moreover, it is known that the filter gains at low frequencies are larger with increasing number of loudspeaker units because of the mode strength at high orders [5,7,29]. Consequently, we must pay attention to the filter coefficient magnitudes in the filter design for a loudspeaker array consisting of a large number of loudspeakers.
8.
CONCLUSIONS
For the PMM based on the least-squares method for directivity control, we proposed the use of a window function shape to form the desired directivity pattern having
249
Acoust. Sci. & Tech. 39, 3 (2018)
a wide main lobe. Furthermore, we proposed a full-search algorithm with a cost function, which defines a criterion for the directivity, to obtain optimum filter coefficients for the suppression of sidelobe levels at each frequency. By using defined cost functions, we maximized the energy ratio between the target and sidelobe areas, minimized the energy in the sidelobe area, and minimized the maximum sidelobe level in the sidelobe area. We obtained the directivity for the proposed method using either a Hanning or rectangular window as the desired directivity pattern. Next, we evaluated the directivity performance in detail. From the results, we confirmed that the proposed method suppressed the maximum sidelobe level by about 2 dB, although the beam widths at 10 to 20 were wider than those for the conventional method. The use of a rectangular window was better than the conventional method for sidelobe suppression. However, there was no significant difference from the results obtained using the Hanning window. Overall, we verified that the desired directivity pattern with a wide main lobe using the window function shape is effective for sidelobe suppression.
ACKNOWLEDGEMENTS This work was supported by JSPS KAKENHI Grant Number 16K00232 and by SCAT. REFERENCES [1] S. J. Elliott, J. Cheer, H. Murfet and K. R. Holland, ‘‘Minimally radiating sources for personal audio,’’ J. Acoust. Soc. Am., 128, 1721–1728 (2010). [2] S. J. Elliott, J. Cheer, J.-W. Choi and Y. Kim, ‘‘Robustness and regularization of personal audio systems,’’ IEEE Trans. Audio Speech Lang. Process., 20, 2123–2133 (2012). [3] M. F. Simo´n Ga´lvez, S. J. Elliott and J. Cheer, ‘‘The effect of reverberation on personal audio devices,’’ J. Acoust. Soc. Am., 135, 2654–2663 (2014). [4] K. Imaizumi and Y. Haneda, ‘‘Wearable sound reproduction system using two end-fire arrays,’’ Proc. 142th Conv. Audio Eng. Soc., 8 pages (2017). [5] K. Sato and Y. Haneda, ‘‘Filter design of a circular loudspeaker array considering the three-dimensional directivity patterns reproduced by circular harmonic modes,’’ Proc. 142th Conv. Audio Eng. Soc., 10 pages (2017). [6] M. M. Boone, W.-H. Cho and J.-G. Ih, ‘‘Design of a highly directional endfire loudspeaker array,’’ J. Audio Eng. Soc., 57, 309–325 (2009). [7] F. M. Fazi, M. Shin, F. Olivieri and S. Fontana, ‘‘Low frequency performance of circular loudspeaker arrays,’’ Proc. 138th Conv. Audio Eng. Soc., 4 pages (2015). [8] M. Kolundzˇija, C. Faller and M. Vetterli, ‘‘Baffled circular loudspeaker array with broadband high directivity,’’ Proc. ICASSP 2010, pp. 73–76 (2010). [9] M. Kolundzˇija, C. Faller and M. Vetterli, ‘‘Design of a compact cylindrical loudspeaker array for spatial sound reproduction,’’ Proc. 130th Conv. Audio Eng. Soc., 10 pages (2011). [10] M. A. Poletti and T. Betlehem, ‘‘Design of a prototype variable directivity loudspeaker for improved surround sound reproduction in rooms,’’ Proc. 52nd Int. Conf. Audio Eng. Soc., 10
250
pages (2013). [11] H. Teutsch and W. Kellermann, ‘‘Acoustic source detection and localization based on wavefield decomposition using circular microphone arrays,’’ J. Acoust. Soc. Am., 120, 2724– 2736 (2006). [12] E. Mabande, A. Schad and W. Kellermann, ‘‘Design of robust superdirective beamformers as a convex optimization problem,’’ Proc. ICASSP 2009, pp. 77–80 (2009). [13] M. Crocco and A. Trucco, ‘‘Design of robust superdirective arrays with a tunable tradeoff between directivity and frequency-invariance,’’ IEEE Trans. Signal Process., 59, 2169–2181 (2011). [14] S. Gannot, E. Vincent, S. M. Golan and A. Ozerov, ‘‘A consolidated perspective on multimicrophone speech enhancement and source separation,’’ IEEE/ACM Trans. Audio Speech Lang. Process., 25, 692–730 (2017). [15] G. Huang, J. Benesty and J. Chen, ‘‘On the design of frequency-invariant beampatterns with uniform circular microphone arrays,’’ IEEE/ACM Trans. Audio Speech Lang. Process., 25, 1140–1153 (2017). [16] C. M. Bishop, Pattern Recognition and Machine Learning (Springer, New York, 2006). [17] O. Kirkeby and P. A. Nelson, ‘‘Reproduction of plane wave sound fields,’’ J. Acoust. Soc. Am., 94, 2992–3000 (1993). [18] O. Kirkeby, P. A. Nelson, F. Orduna-Bustamante and H. Hamada, ‘‘Local sound field reproduction using digital signal processing,’’ J. Acoust. Soc. Am., 100, 1584–1593 (1996). [19] M. F. Simo´n Ga´lvez, S. J. Elliott and J. Cheer, ‘‘A superdirective array of phase shift sources,’’ J. Acoust. Soc. Am., 132, 746–756 (2012). [20] F. Olivieri, F. M. Fazi, M. Shin and P. A. Nelson, ‘‘Pressurematching beamforming method for loudspeaker arrays with frequency dependent selection of control points,’’ Proc. 138th Conv. Audio Eng. Soc., 10 pages (2015). [21] M. Miyoshi and Y. Kaneda, ‘‘Inverse filtering of room acoustics,’’ IEEE Trans. Acoust. Speech Signal Process., 36, 145–152 (1988). [22] E. G. Williams, Fourier Acoustics: Sound Radiation and Nearfield Acoustical Holography (Academic Press, New York, 1999). [23] M. A. Poletti, ‘‘Three-dimensional surround sound systems based on spherical harmonics,’’ J. Audio Eng. Soc., 53, 1004– 1025 (2005). [24] A. Gupta and T. D. Abhayapala, ‘‘Three-dimensional sound field reproduction using multiple circular loudspeaker arrays,’’ IEEE Trans. Audio Speech Lang. Process., 19, 1149–1159 (2011). [25] F. M. Fazi, M. Shin, F. Olivieri, S. Fontana and L. Yue, ‘‘Comparison of pressure-matching and mode-matching beamforming methods for circular loudspeaker arrays,’’ Proc. 137th Conv. Audio Eng. Soc., 9 pages (2014). [26] F.-M. Hoffmann, F. M. Fazi and S. Fontana, ‘‘Sound field control with hemi-cylindrical loudspeaker arrays,’’ Proc. Conf. Audio Eng. Soc., 8 pages (2016). [27] T. Okamoto, ‘‘Analytical methods of generating multiple sound zones for open and baffled circular loudspeaker arrays,’’ Proc. IEEE Workshop on Appl. Signal Process. Audio Acoust. (WASPAA), 5 pages (2015). [28] T. Okamoto and A. Sakaguchi, ‘‘Experimental validation of spatial Fourier transform-based multiple sound zone generation with a linear loudspeaker array,’’ J. Acoust. Soc. Am., 141, 1769–1780 (2017). [29] D. P. Jarrett, E. A. P. Habets and P. A. Naylor, Theory and Applications of Spherical Microphone Array Processing (Springer, Switzerland, 2017).
K. SATO and Y. HANEDA: SIDELOBE SUPPRESSION FOR SMALL CIRCULAR ARRAY [30] F. Olivieri, F. M. Fazi, P. A. Nelson and S. Fontana, ‘‘Comparison of strategies for accurate reproduction of a target signal with compact arrays of loudspeakers for the generation of zones of private sound and silence,’’ J. Audio Eng. Soc., 64, 905–917 (2016). [31] C. C. Gaudes, I. Santamarı´a, J. Vı´a, E. M. Go´mez and T. S. Paules, ‘‘Robust array beamforming with sidelobe control using support vector machines,’’ IEEE Trans. Signal Process., 55, 574–584 (2007). [32] M. Møller, M. Olsen, F. Agerkvist, J. Dyreby and G. Munch, ‘‘Circular loudspeaker array with controllable directivity,’’ Proc. 128th Conv. Audio Eng. Soc., 12 pages (2010). [33] H. Nakajima, ‘‘Minimum sidelobe beamformer design based on Mini-Max criterion,’’ J. Acoust. Soc. Jpn. (J), 63, 341–352
(2007) (in Japanese). [34] D. H. Johnson and D. E. Dudgen, Array Signal Processing: Concepts and Techniques (Prentice Hall, New Jersey, 1993). [35] H. Cox, R. M. Zeskind and M. M. Owen, ‘‘Robust adaptive beamforming,’’ IEEE Trans. Acoust. Speech Signal. Process., 35, 1365–1376 (1987). [36] J.-W. Choi and Y.-H. Kim, ‘‘Generation of an acoustically bright zone with an illuminated region using multiple sources,’’ J. Acoust. Soc. Am., 111, 1695–1700 (2002). [37] M. Shin, S. Q. Lee, F. M. Fazi, P. A. Nelson, D. Kim, S. Wang, K. H. Park and J. Seo, ‘‘Maximization of acoustic energy difference between two spaces,’’ J. Acoust. Soc. Am., 128, 121– 131 (2010).
251