Sparse Plane Wave Approximation of Acoustic Modes to Address Basis Mismatch

: Low-frequency sound ﬁeld reconstruction in an enclosed space has many applications where the plane wave approximation of acoustic modes plays a crucial role. However, the basis mismatch of the plane wave directions degrades the approximation accuracy. In this study, a two-stage method combining (cid:96) 1 -norm relaxation and parametric sparse Bayesian learning is proposed to address this problem. This method involves selecting sparse dominant plane wave directions from pre-discretized directions and constructing a parameterized dictionary of low dimensionality. This dictionary is used to re-estimate the plane wave complex amplitudes and directions based on the sparse Bayesian framework using the variational Bayesian expectation and maximization method. Numerical simulations show that the proposed method can efﬁciently optimize the plane wave directions to reduce the basis mismatch and improve acoustic mode approximation accuracy. The proposed method involves slightly increased computational cost but obtains a higher reconstruction accuracy at extrapolated ﬁeld points and is more robust under low signal-to-noise ratios compared with conventional methods. Author Contributions: Conceptualization, J.X.; methodology, J.X.; software, J.X. and L.W.; validation, J.X. and L.W.; investigation, J.X.; writing—original draft preparation, J.X.; writing—review and editing, J.X., L.W. and J.Z.; visualization, J.X.; supervision, K.C.; funding K.C. J.Z.


Introduction
Recently, the reconstruction of a low-frequency sound field in an enclosed space has generated interest among researchers. It is significant for many applications, including sound field auralization [1], room equalization [2], sound field control [3][4][5], etc. The reconstruction problem usually involves interpolating the impulse responses at unmeasured positions from the impulse responses measured by a discrete set of microphones. The interpolation methodology can avoid the extensive experimental effort caused by dense point-by-point measurements. However, according to the study on the plenacoustic function [6], the Nyquist sampling theorem must be applied, requiring multiple microphonebased measurements for higher cut-off frequency. This contrasts with the goal of employing the minimum number of microphones to recover the sampled sound field faithfully.
The framework of compressive sensing (CS) [7,8] compensates for the lack of observation information by exploiting prior knowledge, so the inherent sparsity of a low-frequency enclosed sound field can help reduce the number of microphones. Under this framework, most methods formulate a physical propagation model and expand the limited set of measurements into finite-order elementary waves [9]. Commonly used elementary wave (or propagation model) includes the spherical harmonic function [10,11], plane wave (PW) [11][12][13][14][15][16][17], and free-field Green function [13]. The PW approximation is widely used because Vekua's theory [18] shows that any homogeneous sound field in an enclosed convex region can be approximated by a finite number of PWs, independent of boundary conditions [19].
The current PW approximation methods can be classified into two types. The frequencydomain method [11,[13][14][15][16] approximates the steady-state sound field by a finite sum of PWs for each frequency and uses the 1 -norm relaxation to solve the PW complex amplitudes (PWAs). In previous studies working with the frequency domain, different sparsity promoting regularizations were compared [13], and a two-step 1 -norm minimization method was developed by expressing PWs in the spherical harmonic domain to address the high mutual coherence of the sensing matrix when PWs are densely discretized [14]. The modal structure was incorporated into the PW model, and the proposed method could predict frequency response functions between the unmeasured source-receiver pairs [16]. The time-domain method [12,17] decomposes the homogeneous wave equation into the superposition of damped harmonic eigen-PWs. It first uses the matching pursuit or the rational fraction polynomials global curve fitting method to identify the acoustic modes, and then conducts the sparse PW approximation for each mode. Compared with the sound field in the frequency domain, the eigen-PWs of a single mode have substantially better sparsity, especially in rectangular enclosures. Therefore, this study focuses on the time-domain method's PW approximation of acoustic modes.
However, the PW approximation has an inherent drawback in that the spatially continuous PW directions must be discretized, and regardless of how finely the directions are discretized, there is always a mismatch between the pre-discretized direction and the actual direction, which is referred to as the basis mismatch [20]. This will inevitably reduce the approximation accuracy. A dictionary or sensing matrix comprising dense PW directions is required to address this issue. This produces a high mutual coherence for the sensing matrix and has the disadvantage of increased computational time. A standard method used to tackle the basis mismatch in acoustics is to adaptively optimize pre-discretized grids based on parametric dictionary learning. For inhomogeneous sound field recording, a grid-pruning-based Newton's method was utilized to optimize monopole positions [21]. Joint sparse recovery and a parametric dictionary learning problem were formulated based on the Bayesian framework for source localization [22]. While locating sources in a highly reverberant environment, the unknown room reflection energy ratio has been estimated. The parametric sparse Bayesian learning (SBL) has also been used to simultaneously estimate the unknown monopole positions and path loss exponent [23]. The pre-discretized grids near the sources were first selected, then a combination of the twodimensional Newton optimization and a feedback mechanism was used to approximate the actual positions for beamforming [24]. In this study, PW discretization involves adaptively optimizing the propagating directions, including inclinations and azimuths, which have rarely been focused upon in previously published literature.
This study proposes a combined two-stage method to address the problem of basis mismatch in the PW approximation of acoustic modes. First, we use 1 -norm relaxation to solve the coarse PWAs, select dominant PW directions, and prune the other directions. Then, we use the extracted dominant directions to construct a low-dimensional parameterized dictionary that considers the PW directions as unknown parameters instead of densely discretizing the direction space to construct a high-dimensional known dictionary. We solve the PWAs and optimize the PW directions by applying SBL. This is facilitated by the variational Bayesian expectation and maximization (VBEM) technique [22,25,26], where the two-dimensional Newton optimization is derived to adaptively modify the inclinations and azimuths. A similar method has not been used in the field of sound reconstruction thus far. Numerical simulations show that the proposed method can reduce the basis mismatch in the PW approximation and improve the approximation accuracy compared with conventional methods, where the computational cost is slightly increased. Additional benefits of the proposed method include a high reconstruction accuracy at extrapolated field points outside the local measurement region and robustness under low signal-to-noise ratios (SNRs).
The paper is organized as follows. The problem of the PW approximation of acoustic modes is stated in Section 2, and the proposed two-stage method is described in Section 3.
The simulation results are presented to validate the effectiveness of the proposed method in Section 4. A discussion is presented in Section 5, and conclusions are summarized in Section 6.

PW Approximation of Acoustic Modes
An acoustic mode function ψ n of order n can be approximated by a finite sum of PWs sharing the same wavenumber k n . The R-order approximation of ψ n at position r = [x, y, z] T is described as follows [12]: where the unit vector s r = s(θ r , φ r ) = [sin θ r cos φ r , sin θ r sin φ r , cos θ r ] T represents the direction of the rth PW with inclination θ r ∈ [0 • , 180 • ] and azimuth φ r ∈ [0 • , 360 • ); (·) T denotes the transpose operation; k n s r is the rth wave vector of the nth mode; and w n,r is the corresponding PWA. For a complex k n (damping wall), we assume that Equation (1) remains valid at locations far from the wall [12]. The problem is to estimate the values of {w n,r } R r=1 from limited measurements. Assuming that M samples are employed, Equation (1) can be converted to a matrix form as follows: whereH n is the M × R matrix, with its (m, r)th entry being h n,m (s r ) = exp(−ik n s T r r m ), and w n is the R × 1 vector of the PWAs. The subscript (·) n is omitted below for simplicity of notation.
The discretized PW directions, S = {s r } R r=1 , are uniformly sampled on the surface of a unit sphere. Generally, R > M, andw can be solved as where (·) H denotes the conjugate transpose operation, which is the analytic solution of the underdetermined least-squares (LS) evaluated using the Tikhonov regularization. However, this solution tends to generate spatial aliasing artifacts. An alternative solution based on the CS framework exploits the prior knowledge thatw is sparse. In this case, 1 -norm relaxation can be used to formulate an 1 -norm-regularized LS, also known as Lasso [27]. Its solution,w where || · || 2 and || · || 1 represent the 2 -norm and 1 -norm, respectively, has few non-zero values and can be solved using iterative algorithms [28]. In our study, we applied the iterative reweighted least squares (IRLS) algorithm [29]. R is usually set to be considerably large to reduce the basis mismatch and obtain sufficiently dense PW directions. Thus, a high-dimensional dictionary was constructed. However, implementing this method is computationally expensive and even leads to failure of approximation.

Proposed Combined Two-Stage Method
In this section, we describe the proposed two-stage method combining 1 -norm relaxation and parametric SBL for the PW approximation of modes with reduced basis mismatch.

Stage 1: Selection of Dominant PW Directions Using 1 -Norm Relaxation
We used 1 -norm relaxation to estimate the PWAsw with respect to the R discretized PW directions in S. We then selected the dominant directions S ⊂ S whose corresponding absolute PWAs were greater than a certain percentage of the maximum, namely, where β is the percentage, called the selection threshold, and card(·) denotes the cardinality of the set, that is, the number of selected directions. Owing to the sparsity ofw L1 , R R. We also used S to form a new dictionary H of low dimensionality, M × R , which is addressed in Stage 2. This stage is intended to reduce the calculation cost of the parametric SBL.

Parameterized Dictionary
To avoid dense discretized PW directions, we explicitly modeled the dictionary with unknown inclinations and azimuths. The high dimensionality of the known dictionary can be reduced at the expense of introducing unknown parameters; therefore, we can apply the low-dimensional dictionary H here. Let h(s r ) be the rth atom of the dictionary H corresponding to the direction s r . We assumed that the actual PW direction was s l / ∈ S for some l ∈ {1, . . . , L} and that s r l ∈ S was the direction nearest to s l . In this case, we can approximate the unknown h(s l ) using the multivariable Taylor expansion in the pre-discretized direction s r l [30]: where h θ r l (s r l ) = ∂h(s r l )/∂θ r l , h φ r l (s r l ) = ∂h(s r l )/∂φ r l , ϑ l = θ l − θ r l , and ϕ l = φ l − φ r l . Extending Equation (6), we construct the parameterized dictionary as follows: where The problem in Equation (2) is reformulated as The task is to estimate the PWAs w and unknown parameters ϑ and ϕ. This problem can be solved using a sparse Bayesian framework in a probabilistic manner.

Probabilistic Model for SBL
From the SBL perspective, we attempt to find a function that maps arbitrary inputs r (field point positions) to targets ψ (acoustic mode values) based on training data {ψ(r m )} M m=1 (observations at samplings) to solve Equation (8) [26]. This function depends on a set of parameters w (PWAs).
Within the sparse Bayesian framework, we first developed a probabilistic model. We assumed a complex, Gaussian-distributed additive noise with zero mean and precision τ according to Equation (8). The likelihood function of the model ψ = H(ϑ, ϕ)w + then becomes where CN (·) indicates the complex Gaussian distribution and I M is an identity matrix of size M. A Gamma prior Ga(·) with parameters a and b is considered over τ in this study for tractable inference of τ described as follows: The components of PWAs w are probabilistically independent of each other. Assuming that the rth entry follows a zero-mean, complex Gaussian prior with precision α r ; this leads to the equation below, Furthermore, a Gamma prior with parameters c and d is imposed on each hyperparameter α r , which yields Ga(α r |c, d).
(12) Figure 1 shows the described probabilistic model of the variables and their hierarchical relationships.

Hidden Random Variable Inference
We intended to obtain the posterior P(w, τ, α|ψ) to evaluate inferences regarding hidden variables H = {w, τ, α} with the probabilistic model. However, the closed-form expression of the posterior is generally intractable owing to the difficulty in calculating the multidimensional integral when computing the model evidence P(ψ). Therefore, we resorted to the VBEM method, which involves recursively approximating the posterior to infer H at the expectation step and optimizing ϑ and ϕ at the maximization step until convergence.
We assumed a distribution Q(H) with the factorized form in the expectation calculation step, that is, Q(H) = Q 1 (w)Q 2 (τ)Q 3 (α). The posterior P(H|ψ) can then be approximated by minimizing the KL divergence between P(H|ψ) and Q(H) for given parameters ϑ and ϕ, which is equivalent to solving the following equation where E Q(H\H j ) [·] represents the expectation of ∏ i =j Q i (H i ). Substituting the priors of H, the updates of each H at the ith iteration can be summarized as follows [26]: where denotes the component-wise division, and Tr(·) represents the trace operation. At the maximization step, we estimated the unknown parameters ϑ and ϕ separately by maximizing the expected log-likelihood for one parameter while keeping the other fixed. This is equivalent to solving the equation below, where x represents ϑ or ϕ. We used Newton's method to solve this nonlinear problem. Accordingly, the single-step update of x is derived as where {x, y} represents {ϑ, θ} or {ϕ, φ}; ∇L(x) and ∇ 2 L(x) denote the gradient and Hessian matrix of function L at x, respectively; [·] is the real part of the argument; (·) * denotes the conjugate operation; • refers to the Hadamard product; and the operation diag(diag(·)) is equivalent to setting all entries of the matrix, except the diagonals, to zero. Note that the application of Newton's method must ensure that the Hessian matrix ∇ 2 L(x) is positive definite; otherwise, Equation (19) must be modified accordingly [31]. In addition, each entry of x should be constrained within an appropriate interval to prevent excessive angle adjustments at the beginning of the iterative procedure, for example [−x lim , x lim ]. If x > x lim or x < −x lim , we set x = x lim or x = −x lim . In our study, we set each entry of x lim to be the same and equal to 2 arcsin(D/2) , where D is the average Euclidean distance between adjacent points {s r } R r=1 . The second stage of the proposed method is summarized in Algorithm 1. Note that the solutionw L1 ∈ C R×1 in Stage 1 is only used to select the dominant directions, and w ∈ C R ×1 is initially assigned to 0 R . The results include the estimated PWAs and directions, which are denoted as {w, S }. Figure 2 shows a schematic illustration of the proposed method, L1-parSBL. Algorithm 1: Stage 2: Estimation of PWs using parametric SBL (parSBL).
with the PW directions ,

Stage 2 parSBL
Select dominant PW directions to construct a new PW dictionary of low dimensionality with the optimized PW directions Estimate the PWAs using parSBL a vector of

Validation Case: A Rigid-Walled Rectangular Enclosure
In our numerical simulation, we considered a rigid-walled rectangular enclosure of 3.7 m × 1.8 m × 1.2 m. Each mode can be evaluated as a sum of eight PWs whose directions can be expressed analytically [32]. The results obtained by the proposed L1-parSBL were compared with conventional methods to examine the effectiveness and performance superiority of our method.
We employed 100 sets of M = 40 randomly distributed samplings to measure the modes to eliminate the influence of random sampling positions. As it is difficult to calculate the estimation error of the PWAs directly, we used them further to reconstruct the mode and indirectly evaluate the relative mean error (RME) between the reconstructed mode values (ψ rec ) and those desired (ψ des ) at predefined field points in the volume as follows [14]: In this case, the field points were set as 30 × 15 × 10 equally-spaced grid points in the enclosure. The selection threshold β was set to 0.1 because we clearly know that there are definite PWs objectively in this case. The discretized PW directions were obtained by solving the Thomson problem [33]. We examined two sets of pre-discretized PW directions of R = 400 and 1000, and the corresponding x lim s are 11 • and 7 • , respectively. Figure 3 shows the RMEs corresponding to the first 100 modes for different methods, with R = 400 and 1000. The methods evaluated were LS, 1 -norm relaxation (L1), L1-L1, and the proposed L1-parSBL. L1-L1 denotes the method of selecting dominant PW directions based on L1 and re-estimating the PWAs based on the low-dimensional dictionary, according to L1. We observed that L1-parSBL outperformed the other methods, and even L1-parSBL with R = 400 yielded better outcomes than the other methods with R = 1000. All methods except LS provided comparable RMEs for low-order modes. LS exhibited very poor performance and could not provide an acceptable reconstruction for modes above the 22nd order (RME > 0.3). Although a maximum of eight PW directions exist, smaller eigenwavelengths relative to the enclosure size cause the mode shapes to become more complex and more sensitive to the basis mismatch for higher-order modes. Thus, the RMEs of the other methods were increased. L1-parSBL could reduce the basis mismatch and achieved a better RME because it optimized the PW directions in Stage 2. In addition, the RMEs of different higher-order modes were considerably different, corresponding to different mode shapes. Axial and tangential modes tended to be better reconstructed than oblique modes. For oblique modes, more sample data points were needed or samples needed to be arranged properly to capture sufficient spatial information.  Specifically, we considered the high-order oblique mode (1, 3, 1) as an example to compare the approximation performance of these methods. Figure 4a shows the desired mode shape. Figure 4b-d plots the relative errors of the approximation at field points (R = 400), |ψ des − ψ rec | |ψ des |, for LS, L1, and L1-parSBL, which achieved RME values of 0.73, 0.26, and 0.09, respectively. The result of L1-L1 was similar to that of L1 and is not shown here. It was observed that L1-parSBL produced a better reconstruction of the mode shape. In Figure 5a-d, note that the numbers of estimated PWs for LS and L1, {w LS , S} and {w L1 , S}, respectively, were both R = 400; in contrast, for L1-L1, there existed only R = 23 estimated PWs, {w L1-L1 , S }, whose directions were selected in Stage 1. Moreover, L1-parSBL yielded R = 23 estimated PWs, {w L1-parSBL , S }, whose directions were selected in Stage 1 and optimized in Stage 2. LS could not identify the actual PW directions because of which spatial aliasing artifacts were produced. The PWs could not be accurately estimated because of the basis mismatch, although L1 and L1-L1 could identify approximate directions. The proposed L1-parSBL efficiently optimized the PW directions to ensure that they were closer to the actual ones in Stage 2 and L1-parSBL estimated the PWAs more precisely than other methods. Similar results were obtained for the other modes as well. Thus, outcomes for only one mode are illustrated.
The elapsed times of L1, L1-L1, and L1-parSBL for the PW approximation of the (1, 3, 1) mode are listed in Table 1. The computer included two Intel Xeon central processing units with two 2.20 GHz processors. For the same R, L1-parSBL required a slightly longer time as observed from Table 1. This was because in Stage 2, parSBL handled a very low-dimensional dictionary (R R), which was constructed in Stage 1.

Case: An Aircraft Cabin with Damping Floor
In this subsection, we considered a more realistic enclosure environment for an aircraft cabin, as shown in Figure 6a. The floor was set as a damping boundary with an absorption coefficient α = 0.1. In practice, we generally focus more on local regions. Therefore, a virtual parallelepiped of 4 m × 2 m × 0.8 m with the center at [4, 0, 1.4] T was defined as the sampling region. We used the Eigenfrequency study in the Pressure Acoustics, Frequency Domain physics interface in COMSOL Multiphysics for simulation. The vertices of the meshes inside the virtual parallelepiped region were chosen as the interpolated field points for assessment, and those outside the region were chosen as the extrapolated field points. We used 100 sets of M = 40 samplings. Zero-mean, additive, Gaussian noise with an SNR of 25 dB was present during each sampling. We examined 20 modes around 107.5 Hz (the first-order blade-passing frequency) whose mode shapes are plotted in Figure 6b. Figure 7 shows the effect of the selection threshold β on the reconstruction RME. We divided the interval [10 −3 , 10 −0.3 ] of β into 20 grids which were equally-spaced in a logarithmic scale. It could be observed that the RME increased as β increased, because there were fewer PW directions selected, leading to the lack of some necessary information. When β was less than 0.1, the RME hardly changed with it, for both interpolated or extrapolated field points. However, with smaller β, more trivial PW directions were selected. Thus, we selected β to be 0.01 empirically, which meant that the PW directions whose absolute PWAs were less than 1% of the maximum absolute PWA were discarded.  Figure 8 shows the RMEs of LS, L1, and L1-parSBL for the 20 modes inside and outside the virtual region, respectively, with R = 400 and 1000. We observed that L1-parSBL achieved better RMEs than LS and L1 inside and outside the virtual region for the 20 modes. The RMEs of LS were all greater than 0.5, and L1 could obtain RMEs below 0.3 under four modes for the extrapolated field points outside the region. In contrast, L1-parSBL achieved acceptable reconstruction errors (RME < 0.3) for 13 modes. We also found that the denser discretized PW directions hardly improved the RME in this case. Additionally, the RMEs of different modes were very different, especially for the extrapolated field points outside the virtual region. This was because samplings were distributed in the local region and these samplings could not obtain sufficient spatial information for some modes. For example, for modes "No. 3" and "No. 8" of Figure 6b, the virtual region mainly included their nodal planes, so it was difficult for the samplings to capture the effective information of the mode shapes. Therefore, the estimated PWs lead to poor reconstruction accuracy outside the virtual region. Specifically, we considered the mode at approximately 106.5 Hz (closest to 107.5 Hz) as an example. Corresponding desired mode shape is plotted as the 12th subgraph in Figure 6b.
The effects of M (the number of sample data points) and SNR on the reconstruction accuracy of LS, L1, and L1-parSBL were analyzed. Figure 9a shows that the RMEs inside and outside the virtual region increased as M decreased for LS, L1, and L1-parSBL, respectively. Although the RMEs inside the region for the three methods were close, L1-parSBL had the lowest reconstruction errors. Both L1 and L1-parSBL could provide acceptable outcomes with M = 10 samplings, but according to the error bar curves, specific sampling positions had significant impact. The proposed L1-parSBL achieved much better RMEs than LS and L1 when sufficient number of samples were used for the extrapolated field points. Based on Figure 9b, the RMEs inside and outside the virtual region increased as SNR decreased for LS, L1, and L1-parSBL, respectively. LS and L1 failed to reconstruct the mode shape outside the region (RME > 0.5) when SNR is below 25 dB. However, L1-parSBL still obtained the average RME of 0.5 at 15 dB. For interpolated field points, L1-parSBL could provide RME of 0.34 at 5 dB, whereas the RMEs of LS and L1 were basically more than 0.5 at 10 dB. Although only one mode outcome is illustrated here, similar results were obtained for the other modes as well.  Figure 10 presents the plots of the estimated PWs of LS, L1, and L1-parSBL under different SNRs for the "No. 12" mode at approximately 106.5 Hz by using a set of randomly distributed samplings (R = 400). Note that there are no analytical PWs for this cabin case. The 12th subgraph in Figure 6b reveals that the mode shape had a significant characteristic of sparse PW superposition. There was no change in the modal values along the z-axis, regardless of the damping floor, indicating that the corresponding PWs of the mode were mainly concentrated on the equator of the unit sphere. It can be seen that LS could not identify sparse PWs. When SNR was 30 dB, the estimated PWs of LS barely exhibited four energy highlands around the equator. For lower SNRs, the estimated PWs distorted severely. A comparison between L1 and L1-parSBL shows that they both estimated PWs well at 30 dB. For lower SNRs, L1 obtained numerous PW artifacts that deviated from the equator. In contrast, the PWs estimated by L1-parSBL under different SNRs were similar; this reflected the robustness of L1-parSBL to some degree. At SNR = 10 dB, the PW estimation error of L1-parSBL was large. This could be partly attributed to the first stage estimate being wrong; thus, Stage 2 achieved no good results. L1-parSBL eliminated many artifacts selected in Stage 1 and significantly improved the reconstruction accuracy.

Discussion
The proposed L1-parSBL was used in this study to address the adverse effect of the basis mismatch in PW discretization and improve the approximation accuracy of acoustic modes compared with conventional methods. L1-parSBL suffered from smaller reconstruction errors at extrapolated field points outside the region in a local sampling region because it estimated the corresponding PWs more precisely.
Note that no sparsity was assumed in Equation (1), and the finite sum of Equation (1) was an approximation required for computation [12]. The sparsity depends on the consistency between the enclosure geometry and the physical propagation model employed. For example, a cylindrical cavity will barely have a sparse representation with PWs (except in the z-axis) because PWs do not sparsely represent Bessel functions. The performance of CS-based methods deteriorates when the PWs of an acoustic mode are not very sparse. In this case, the selection threshold β should be set with further insight [34,35].
The randomly distributed samplings must be arranged such that sufficient information of the sound field is captured. In addition, Stage 2 depended on the results of Stage 1. If the estimation was wrong in Stage 1 (for example, due to low SNR), it might be difficult to obtain good results in Stage 2. The IRLS was used here to calculate the coarse PWAs,w L1 , in the first stage. Other algorithms need to be compared to evaluate an optimal paradigm, especially in terms of robustness.

Conclusions
We proposed a two-stage method, L1-parSBL, that combines 1 -norm relaxation with parametric SBL for the PW approximation of acoustic modes to reduce the adverse effect of basis mismatch. The first stage was used to select the sparse dominant PW directions from pre-discretized directions. The second stage was performed under the Bayesian framework to reformulate a parameterized dictionary of low dimensionality to re-estimate PWAs while optimizing the PW directions to approximate the actual directions. Numerical simulations were conducted using two models of an enclosure. The proposed L1-parSBL achieved a considerably higher approximation accuracy of acoustic modes than LS, L1, and L1-L1 for the rectangular enclosure with rigid walls. Based on an analysis of its (1, 3, 1) mode, it was illustrated that L1-parSBL could efficiently optimize the PW directions to reduce the effect of basis mismatch, albeit at the cost of a slightly increased computation time. A local measurement region was examined considering an enclosed aircraft cabin with a damping boundary. The results showed that L1-parSBL outperformed LS and L1, especially for extrapolated field points outside the region. L1-parSBL also exhibited much better robustness in low-SNR cases. The proposed L1-parSBL can further assist in enhancing the performance of reconstruction for sound fields, and related experimental work will be performed as part of future research.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.