Acessibilidade / Reportar erro

The fundamental elements in certain inverse acoustic problems: their roles and interactions

Abstract

Acoustic holography and holophony, wave field synthesis and active noise control are based on common elements which are causality, model, objective, and regularization. In the frequency domain (putting causality aside), a simple formulation states the influence - not the interaction - of errors of the model and objective and of regularization of the results. However, it does not give either an understanding or any relation of cause to effect. When the objective can be reached using the available model, regularization is not needed and the information liable to be extracted from this determined problem is poor, unlike in the over-determined case when the model does not allow the objective to be reached. The geometrical interpretation of the over-determined problem written in the least-mean square sense could be a tool to enlighten the influences and interactions in question. After having shown the interest of the geometrical interpretation, a pseudo-analytical inverse problem in spherical holophony and a numerical problem in plane holography provide particular illustrations. From among the properties accessible, one is highlighted: in the case of a perfect objective but inaccurate model, its adaptation brings a decrease in the amount of regularization required and an improvement in the results.

inverse problems; regularization; acoustics


TECHNICAL PAPERS

CNRS Institut Jean le Rond d'Alembert - UMR CNRS/UPMC 7190 - 75005 Paris, France vincent.martin@upmc.fr

ABSTRACT

Acoustic holography and holophony, wave field synthesis and active noise control are based on common elements which are causality, model, objective, and regularization. In the frequency domain (putting causality aside), a simple formulation states the influence - not the interaction - of errors of the model and objective and of regularization of the results. However, it does not give either an understanding or any relation of cause to effect. When the objective can be reached using the available model, regularization is not needed and the information liable to be extracted from this determined problem is poor, unlike in the over-determined case when the model does not allow the objective to be reached. The geometrical interpretation of the over-determined problem written in the least-mean square sense could be a tool to enlighten the influences and interactions in question. After having shown the interest of the geometrical interpretation, a pseudo-analytical inverse problem in spherical holophony and a numerical problem in plane holography provide particular illustrations. From among the properties accessible, one is highlighted: in the case of a perfect objective but inaccurate model, its adaptation brings a decrease in the amount of regularization required and an improvement in the results.

Keywords: inverse problems, regularization, acoustics

Introduction

Acoustic holography (Metherel et al., 1967; Maynard, Williams and Lee, 1985, for the early works), wave field synthesis (Berkhout, de Vries and Vogel, 1993, for the early work), holophony, and active noise control (Nelson and Elliott, 1992; and Elliott, 2000, textbooks presenting a synthesis) are inverse acoustic problems sharing common elements such as causality, transfer function or model, objective or goal to be reached and regularization. In the frequency domain, the question of causality is put aside.

In these problems, one or a few sources are required to radiate a given acoustic pressure field (for an extensive overview of the past three decades with 170 references, see Wu, 2008). Indeed, sound field synthesis consists in determining driving signals to feed several loudspeakers in order to radiate a particular wavefront, a particular directivity in holophony; active noise control is linked to sound synthesis indirectly as the goal is to oppose a so-called acoustic primary field at the control microphone locations; holography is concerned with reconstructing a field in a whole domain from pressure measurements at a finite number of locations, this reconstruction often resting on the finding of the source strength at the origin of the measured field.

With less sources to control than given pressures to reach, the problem is overdetermined and rests on the least-square (LS) method, while with more sources than pressures leads to underdetermined problems with the least-norm (LN) method (the smallest solution from the infinity due to underdetermination, i.e. dealing with the LS method plus a constraint). Emphasized by Hald (2009), these solutions seem to coincide when regularization is involved. Limiting the investigation to overdetermined problems, the geometrical interpretation accompanies naturally the LS method with the Euclidian distance. The use of this interpretation has proved to be a good tool for investigating some properties such as the guaranteed quality of the procedure for reaching the objective (the given pressure field).

The geometrical view is also closely linked to the description of inverse problems with the vocabulary associated with sets, as presented by Kirsh (1996). For example, the given pressure field could be inaccessible for various physical reasons, all related to the pressure field not belonging to the set of fields the sources can generate. In active control, the secondary sources are not merged with the primary ones at the origin of the given pressure field and they cannot radiate the opposite field (except in the very particular 1D case of controlling plane waves in ducts); in holography primary and secondary sources are often the same, but here an erroneous measured field has no reason to be accessible. The best we can do is to reach the projection of the given field on the set of fields the secondary sources can radiate, leading to the geometrical point of view. The same will occur when the model is erroneous, and therefore, not liable to generate the objective. To generalize, it will be said that the model is not exact when incapable of radiating the exact objective (which does not mean that the radiation itself is erroneous), while the definition of an erroneous or an exact objective is kept natural.

So much for the objective space. Now what about the solution space (that of driving signals or source velocities)? To go from one space to the other, the model and its inverse are needed, and it is in the nature of things that the inverse of the direct operator is ill-conditioned (to the point that the differentiation between the two is through conditioning) which means that errors, whatever their reason for being, are greatly amplified in the solution and in the reconstructed field. The regularization acts on the inversion to limit the amplification and also plays a preponderant role in the spatial resolution of the solution obtained, as was underlined by Nelson and Yoon (2000), by Yoon and Nelson (2000), by Nelson (2001), by Kim and Nelson (2003). Regularization acts first in the solution space; in the objective space, only the consequences can be observed and this will be proved. Besides, we have noticed that the adaptation of the model to better radiate the perfect objective (the measured pressure without any errors) is accompanied by a need for less regularization and a better reconstructed field due to a better determination of the solution (driving signal or source velocities). This could help in finding the relation between solution and objective spaces through regularization and perhaps make use of the geometrical interpretation.

Having this research in mind the formulated framework is given, and then the geometrical interpretation in the objective space is described. Its usefulness will be shown through an application. Contrary to the demonstration of a probable impossibility to describe the model regularization in the objective space, a demonstration of the relation between adaptation of a model and the amount of regularization has still not been found and the difficulties arising therein will be shown. Applications in far-field acoustic holophony and in nearfield acoustic holography illustrate the above properties. Open questions conclude this paper, devoted to an attempt to gather a variety of acoustic inverse problems with their uncertainties.

Nomenclature

Framework

The source or sources whose velocities are sought are called "secondary source(s)" to differentiate their role from the so-called primary source(s) at the origin of the radiated pressure. They differ in sound field synthesis, holophony and active noise control while they are merged in holography. Thus, in synthesis and control, secondary sources cannot radiate the primary field accurately in the whole 3D domain (contrary to cases related to 1D propagation).

As the primary field does not belong to the space of the fields liable to be generated by the secondary sources, it is said that the radiation model cannot be exact, with a clear sense in mathematical set terms. In contrast, there might be such an exact model in holography. Here a perturbed propagator, which is thus not exact, arises not from conceptual reasons but because physical data like sound speed, source and microphone locations are only approximated and/or because the measurements and calculations are not exact.

Concerning the objective, it is said to be exact if it perfectly describes the field to be generated. On a microphone array, it is necessary to have a sufficient number of sensors with a well-chosen distance between them according to the sound spectrum under study, and also strictly no measurement errors on the recorded and signal-processed pressures.

The above problems always rest on equation

where the vector p is the objective made up of the measured pressures at the M microphones of the antenna, the velocity or command vector v is sought at N points or on N sources, and E is the transfer matrix or propagator or model of dimensions M x N (M N) with rank E N. Quite often, the inverse problem is solved in the least-mean square sense to lead to the optimal velocity (or driving signals). Writing , where the asterisk indicates the transpose-conjugate, we have:

which would be the exact solution insofar as E and p are exact and also with H well-conditioned.

Very simple manipulations of Eq. (2) result in the perturbation due to perturbation on the exact model E, on the exact objective p and on the numerical inversion of H. From now on, v is the exact vector solution.

Let us consider the perturbed model and the perturbed objective . As long as is perfectly invertible i.e., as long as the strict equality holds, then according to Eq. (2). This leads to

where approximates the "inverse" of the rectangular matrix . Equation (3) is simple but with the drawback that it requires previous knowledge of E and v, making this equation uninteresting in the real world but a helpful first approach from the theoretical point of view. In fact Eq. (3) will become an equation about the error when prior additional information on and on gives an upper bound.

When matrix is either not invertible or poorly inverted, because ill-conditioned, the pseudo-inverse of is approached by a matrix where the parameter can be adjusted to solve at best the problem while still avoiding unacceptable solutions. The matrix or operator is called regulator, the most famous being that of Tikhonov. If the numerical calculation error is negligible, we have:

which is similar to Eq. (3).

How tends towards p or or, when a calculation is carried out with the adjusted value of named , how tends towards p or or how depends on or , are the questions addressed in this paper.

The geometrical interpretation of the inverse problems dealt with previously by the author could constitute a tool to contribute to answering the questions.

Geometrical Interpretation

The raw material of this section is to be found in Martin and Cariou (1997), in Le Bourdon (2009), and in Martin, Le Bourdon and Arruda (2012).

Inverse problem and projection

In presence of perturbation on the objective and on the model, solution in Eq. (2) results from the algorithm

with the L2-norm in the objective space. This Euclidian (Hermitian norm in the complex field) norm leads naturally to a geometrical point of view as this operation can be seen as the projection of the objective on the hyper-plane spanned by the columns vector of matrix . The minimization leads to the minimal distance which is the square root of the residual part of after the operation. Let us define as the value of when the velocity is zero, that is the part of that has been attained thanks to the calculation of . It appears that . These definitions allow us to read symbolically the projection operation in Fig. 1 where the objective does not belong to the hyper-plane due to a perturbation in p, in E, or in both. In these conditions there exists a non-zero angle (written for the exact p) such that


differing from the general cross-validation parameter defined as (see Hansen, 1998):

.

The efficiency is defined by

Therefore, given and , the die is cast for the optimal value of the velocity and for the efficiency. To improve the efficiency there is no choice but to reduce angle (through ) hoping to better approach the true value v under some assumptions. This angle has proved to be efficient for classifying various propagators (transfer matrices) in the automotive industry by Le Bourdon, Picard, Martin (2009).

Error of the objective

In general, the vector of the acoustic pressure is measured data and, for this reason, inevitably erroneous, while being the only data available. It is said "of reference" or "nominal" and written . The true pressure is somewhere around this nominal value and is of the form . For the time being and for the sake of general notation let us write with seen as a perturbed field relatively to , underlining that the unknown true pressure results from a variation in the available nominal pressure.

The error in the objective is defined by

As emphasized by Martin and Gronier (1998) and by Martin, Le Bourdon and Arruda (2012), this quantity may constitute a vicinity around the reference field containing all the perturbed fields such that their relative error is less than or equal to a given value of . Ideally, a vicinity around should correspond to a vicinity around the efficiency in order to establish a connection between efficiency and error in the objective and, by so doing, to be able to see how the efficiency diverges from that of reference when the error in the reference field increases. In fact, the previous definition of the vicinity around , expressed by , is such that to a given efficiency corresponds an infinite number of different values of spread over the interval . Indeed the efficiency defined by the angle associated with is the same for an infinite number of collinear fields , and, therefore, to an infinite number of values of . On the contrary, for all these collinear fields sharing the same efficiency, a single minimum value belonging to interval was defined by Martin and Gronier (1998). The following definition is thus preferred:

Moreover it has been shown that a given value of corresponds to a single value of the guaranteed efficiency when the perturbed field belongs to the vicinity defined by .

It appears that is the sinus of angle shown in Fig. 2 and therefore, is a measurement of the aperture of the cone centered on the reference objective (Fig. 3).



A point of coordinates can be associated to each perturbed field in the vicinity of the reference field . Let us notice that where is the projection of on . For a set of perturbed fields sharing the same minimal relative error, the values and (see Fig. 3) bound the range of for any field considered in that set. Figure 4 illustrates the theoretical result.


Guaranteed quality of the inverse process in the objective space

For a given reference objective , the reference or nominal quality of the process is defined by the quantity

Were we working in active noise control, would be closely linked to the predicted attenuation. Now let us envisage a perturbation of the nominal objective such that it has a minimal error, of less than or equal to a value . Geometrically speaking, with , this means that, as seen in Fig. 3, the true pressure is located somewhere close to the measured one, in or on the cone centered on the reference objective and of aperture angle . If, in the geometrical representation, is the projection of the reference pressure on the edge of the cone located farthest from the model hyper-plane, the quality of inverse process is always better than that defined by

said to be the guaranteed quality. Therefore, it is straightforward that the guaranteed quality is given by

Equation (12) was found by Martin and Gronier (1998) but via an analytical approach far more complicated than this simple geometrical way. The analytical demonstration is still interesting but no longer as for the result more directly obtained today.

Graphs showing, in logarithmic scale, the guaranteed quality against the reference quality for various values of are given in Fig. 5. The main point here is that the graphs present a plateau for high nominal quality leading us to conclude that, for a given value of , an increase in the nominal quality is not rewarded by a significant increase in the guaranteed quality. When the nominal quality is obtained by adapting a model to make the hyperplane of draw nearer to , this information gives the limit of the effort beyond which there is no need to go further.


It has thus been shown in this section that the geometrical interpretation has proved to be a useful tool for presenting the guaranteed quality of the process of inversion, and also to reveal a cost function for adapting the model. However this tool operates only in the objective space.

Model Adaptation versus Regularization: Observations, Generalization, Difficulties

When the model and/or the objective is/are not exact, there is a need for regularization such as Tikhonov's or of another type, for example by discarding the small eigenvalues of the propagation matrix. When only the model is not exact, the operation of adapting the model consists in reducing the angle to the exact objective. Thanks to case studies, it has been observed on the one hand that the regularization has no effect on and on the other hand that the smaller the angle , the less optimal regularization needed and the better the results. While we have found a demonstration for the first observation, there is still research to be done to find how regularization and errors are related.

Error of the objective and of the model

In very general terms and before speaking of exact or perturbed objective and model, the question is to see if the writing of is meaningful and if there is a solution. Let be the solution space and let be the image of through the model .

When , the writing is meaningful and the LS solution leads to the exact value . Indeed .

When , the writing has no sense for (if , then ) and either exists outside , or it does not exist. If the LS solution , it would generate , necessarily different from . When but where is not located, the same remark can be made.

In the case of errors in the objective only, the problem to solve is more precisely written which has no sense as it is unlikely for . For a well-conditioned matrix H the optimal solution differs from the exact solution v associated with the exact pressure by

as stated by Eq. (3) with E in place of . When H is ill-conditioned the regularized Tikhonov solution is

which differs from the exact solution by

in conformity with Eq. (4) with E in place of .

In the case of errors in the model only, the problem to solve is written . The regularized Tikhonov solution is

It differs from the exact solution by (see Eq. (4) with ).

It is clear that if there exists, when (in a sense to be defined), an optimal value of that tends towards zero in , then (in a sense to be defined). So far, the adaptation of the model carried out to obtain rested on the minimization of angle . And, effectively, the process was accompanied by the observation of . It was tempting to assume that the regularization modified the model in order to reduce . The observation countered this assumption and we demonstrate below in four steps the independence of angle from the coefficient with the consequence that the effect of regularization cannot be described easily in the objective space (Martin and Le Bourdon, 2010), perhaps a natural conclusion in retrospect since is not acted upon, contrarily to its "inverse" now written .

Independence of the angle from

Without jeopardizing generality, suppose only the model is not exact.

1. Let us consider with and where . The space V could contain for example only those variables the norms of which are bounded by a given maximal value. With in place of , Fig. 1 above has given the Euclidian representation with and , and in the case where .

The singular value decomposition (SVD) of is such that with and . Matrix is diagonal in its upper part and zero in its lower part.

Matrix is made up of base vectors of the objective space where the data is located, and matrix is formed with the base vectors of space of variable . When , one cannot find (one cannot invert ) and is replaced by the term called pseudo-inverse of that is . This operation leads to the solution such that it generates the nearest of in the sense of the Hermitian norm in .

As , it results in with . We have also leading to

or

.

2. Now let us consider the functional of Tikhonov which helps in approximating the inverse matrix and is a regularization process: with and .

For a given , arises from minimizing with .

3. Let us write not yet having defined (but the access to which is not difficult) and let us consider , the minimization of which for a given provides . It appears . So, for the time being, functionals and share the same objective and the same optimal value of their variables.

4. But we have

In these conditions, functionals and share the same objective and the same minimal distance and this whatever . Therefore angle such that does not depend on the Tikhonov coefficient . In the case of the Euclidean configuration shown in Fig. 1 this is true for both functionals.

5. Finally, the Tikhonov regularization does not seem to have an influence on angle in space of the objectives, but it plays a role in space V of the solutions. Therefore, we cannot expect to represent easily (if indeed it is possible) the regularization process in the objective space.

Optimal regularization and illustration on a simple example

Going back to , let us notice that, when solving with the L2-norm in the solution space, the results observed in examples (therefore without generalization yet) is whatever and therefore, cannot be used for finding (leaving aside that v is not known). Were we to accept negative values of , there could be a value to make zero. Indeed when i.e. . This can be true only if is collinear with v, a very strong constraint. Its LS value would be where it is true that the dependence on is made explicit, but this value does not minimize and moreover, what has been observed concerns , not.

In fact, the optimal value of is such that the objective is best reached while the amplitude of the velocity is kept reasonable. This qualitative definition has received a heuristic representation through the parametric so-called L-curve which gives the value of against (Kirsh, 1996; Hansen, 1998). The optimal value corresponds to the compromise located at the "corner" of the L-curve where the curvature is maximum. However, the idealized L-curve has not always so clear a form and the location for obtaining should be determined with care.

Let us say that has been defined by Steiner and Hald (2001), as depending on the signal to noise ratio (SNR) seen at the objective and due to all sorts of reasons (measurements, errors in the model, calculation approximations). The SNR is defined by rule of thumb and is not seen as model dependent. However, Hald (2003) wrote with the present notations that "re-scaling of the matrix in turn necessitates a re-scaling of the regularization parameter ". And again to emphasize the state of the art concerning the dependence of on the model, let us refer to Gomes and Hansen (2008) where, according to the type of model chosen - which could be that of Statistically Optimized Nearfied Acoustic Holography (SONAH) to replace the first version of NAH from Maynard et al. (1985), or that of the Boundary Integral Equation Method (BIEM) as developed by Bai (1992) or described by Augusztinovicz and Tournour (1999), or that of Equivalent Source Method (ESM) - , there is a change in the best regularization parameter. The work of Gomes and Hansen is an echo of the earlier work carried out by Williams (2001).

Is the observation concerning the decreasing of with that of made when solving physical inverse problems closely linked to the underlying physics, or is it a more general property? To try to answer that question, let us take a very arbitrary matrix and an arbitrary perfect velocity . In these conditions, the perfect objective is . Let us perturb the term (1,2) of E to obtain with a law such that , i.e. with unstructured uncertainties. First, calculate the angle (calculated through ) , then draw the L-curve, and in the case of a relatively visible "corner", find the , and finally calculate . It must be said that the value of retained is such that the quantity is minimal, which is never far from the maximum curvature but not always at exactly the same location. Figure 6 shows two types of L-curve retained during the process and Fig. 7 shows , and against .



This simple numerical experiment shows that the decreasing of with and that of with could be a general property rather than an observation deriving from the physical nature of the inverse problem dealt with in acoustics. However in this arbitrary example, we see a jump of and a corresponding void in the curve of error in the velocity. Moreover, if several components of matrix E are submitted to variations, the graphs are of far greater complexity as clouds take the place of the curves.

Beyond the observation made on this very simple calculation experiment, how can it be demonstrated that decreases with the angle ? What are the difficulties?

Indeed, let us write

where it can be demonstrated (via the SVD of E ) that when , where it is natural that when . But how can the dependency on be expressed? To this end, the expression of will be needed. By accepting that is the solution of , the expression in terms of such that does not seem an easy task.

While the following case studies as well as the above calculation example have shown that the adaptation of a model using the parameter could be an interesting direction to follow, the definitive interest ought to depend on the demonstration sought and this question seems to be still open.

Case Studies

Semi-analytical spherical holophony in the far field

In Pasqual and Martin, 2012, a primary field due to a source (or sources) in a sphere of centre C is given on the spherical surface (Fig. 8). This area is located in the far field relative to the source.


The given pressure is expanded in a series of P independent functions composed of spherical harmonics considering the coordinate center O. The number of « modes » to describe the primary pressure with a good approximation depends on the coordinate center O. The minimum value of P occurs when O is at the so-called acoustic centre, in general unknown, even when the geometrical centre of the primary source is given.

A secondary source is required to radiate the given primary pressure as well as possible. This source, within the sphere but also far from its surface, is constrained to radiate a controllable given number S of « modes » with ; some of them not radiating efficiently, leaving efficient controllable modes.

Where should it be located to reproduce the primary field as well as possible? The problem is equivalent to that of finding the location of the secondary source centred at O. It should be said that the problem of reproducing a field by a spherical source has also been investigated by Peleg and Rafaely (2011).

We are facing an inverse acoustic problem where the model of the secondary radiation has to be adapted through the secondary source location to reproduce at best the primary field. Due to the secondary source constraints, there are « modes » missing to reach the goal, a number that depends on O. Often, but not always, the minimal difference occurs when O is at the acoustic centre of the primary source(s). This difference will play the role of regularization.

The primary harmonic pressure is developed in the form

where are the Hankel functions and are the spherical harmonics. There are terms; however, the axisymmetry of the pressure field reduces the number of terms to as . For example, the field shown in Fig. 9 is such that there is a need for to radiate it, whatever O (i.e., ). The coefficients can be determined almost exactly providing is sufficiently large. In the example chosen, , therefore in general and here. In other words, the number of modes with non-zero (or negligible) coefficients depends on the location of point O.


After linear indexing such that and with the vector made up of components and vector with components , the primary pressure is also expressed by . As the number of components in the vectors depends on the location of point O, their dimension is kept the same but the number of non-zero components varies. The identification of the coefficients in vector is obtained through the LS-method, but to overcome a difficulty due to the ill-conditioned matrix in the inversion, it is practical to deal with coefficients such that where components of diagonal matrix H are . In these conditions:

These new coefficients make up the spherical wave spectrum of the primary pressure evaluated at .

The identification of leads to their value and, in the present situation, is sufficient to ensure a good if not exact reconstruction of , whatever O.

Similarly the secondary harmonic field is expressed by

written

Simply enough, our goal would be perfectly achieved with if the secondary source were located at the acoustic centre of the primary source and able to radiate modes. However, the secondary source constraints with have been mentioned and moreover, from among the modes, a number do not radiate efficiently.

Indeed, the radiation efficiency of a compact spherical loudspeaker array (CSLA) has the form shown in Fig. 10. Asking of such a source to radiate in the low frequency range with a large n leads to far too high a level of driving signal or velocity. It was decided to discard these disturbing high values of n and consider only those n associated with a reasonable minimum value of the efficiency. To illustrate, working with ka = 0.3 and considering only those modes with , only and are left, therefore, and or 2 for axisymmetric pressure distribution. Let us recall that in the case studied. It appears that 4 sources of the CSLA, placed at the so-called extremal points for hyperinterpolation distributed over the spherical box can be driven to radiate efficiently 4 modes (and 2 sources for 2 modes).


This discarding plays the role of the regularization as it facilitates access to the solution of the inverse problem, and the best solution in the LS-sense for leads to in the following way:

Let us define the regularization quality factor by

which is the inverse of the amount of regularization. Clearly, when thanks to seeking the best location for the secondary source, then the amount of regularization decreases as well as the residue .

In the numerical simulation, the primary field of Fig. 9 was chosen because it cannot be produced by a multipolar source of order . It was obtained by using Eq. (19) with and the wave spectra = 0.2, = 1, = 1. Due to the axisymmetry of , we only deal with points O located on the axis. Therefore only positions given by are considered here. The values of the parameters are those given in the text above added to and thus since . Moreover, the spherical area was sampled with 300 elements and each elementary area weighted in order not to privilege densely sampled regions on the sphere. The position of the secondary source is assumed to be within the range to ensure far-field conditions. Figure 11 shows the norm of the optimal radial velocity without and with regularization; only graph is considered here. The results are given as function of . It is clear that regularization is needed to prevent the loudspeakers from being overloaded. Figure 12 shows the normalized residual norm and the regularization quality factor as functions of (only graph is of interest here). It appears that because the primary field cannot be accurately decomposed into spherical harmonics of order . However is greatly reduced if the secondary source is placed at the optimal location.



Figure 13. Sound pressure field produced by a spherical secondary source placed at the optimal position for reproduction of the primary field in Fig. 9. The shape and color maps indicate respectively the amplitude and phase (in degrees).


The results presented show, in an almost natural way, that the adaptation of the model (here through the source location) is accompanied by a reduction in the regularization needed and improves the reconstruction. The synthesized field is given in Fig. 13.

Numerical nearfield holography

In Martin, Le Bourdon and Pasqual (2011), a thick vibrating disc (Fig. 14) radiates sound in the unbounded 3D space. The acoustic pressure is measured on a plane array of microphones parallel to one side of the disc, called the front face.


From the inverse procedure called here holography, the type of vibration on the front face is deduced from the measured pressure at the microphone array. It is supposed that the disc is embedded in a total front plane made up of and its complementary part . How to take into account the radiation behind this front plane due to the vibration of the rear face and the thickness? This problem originates from the acoustic holography of a rolling wheel embedded in a car body. Despite the slightly different goals, let us quote the work of Schuhmacher et al. (2003), on a quite similar subject.

Presently the radiation model is obtained through the boundary integral method that has the continuous form

leading to the discrete form

where and are the unknown admittances (in a particular sense) respectively on and .

For the numerical simulation, the given modal forms on the faces and on the thickness are

Admittances and are local reactions and depend on the point . Therefore, there are an infinite number of variables to identify and it is impossible to find the true propagator in these conditions. However, it has been observed that the admittance on tends toward zero when the radial distance tends towards infinity. Beyond a certain distance, the admittance is negligible and makes finite meshing of the source plane possible. The number of unknowns in the propagation model is directly linked to the number of points in the meshing of the source plane. It appears unrealistic for computation time reasons to identify each local admittance. To reduce the number of unknowns, the admittance is approached by a known function based on its expected form, and now only a small number of coefficients of the function must be identified. In the case of the thick disc, the following approximation has been chosen:

where the coefficients . Adapting the propagator consists now in identifying the 10 coefficients. However this simplification is at the price of information about the function.

Using the genetic algorithm with the cost function made up of the angle , the coefficients have been obtained approximately and the results are given in Fig. 15.


Table 1 shows the vibrating form obtained for the front face and Fig. 16 presents the same results in another form, namely the optimal value of the Tikhonov regularization parameter against angle , itself depending on the admittance on the front face , as well as the error in the reconstructed velocity against angle .


These results illustrate again the interaction between the adaptation of the model and the value of the regularization parameter with the consequence on the reconstructed velocity.

Conclusion and Open Questions

The present paper has given very general equations linking errors in models as well as in objectives and errors in the solution of certain inverse problems. Some interactions between models, objectives, errors, and regularizations are not visible in these equations while their consequence on the solution is well described. The current research attempts to contribute in finding the actions and interactions at play.

The geometrical interpretation has proved useful to show very easily the guaranteed quality of the inverse procedure, and also to define a cost function (angle ) to adapt the model in cases where the objective is exact. It has been demonstrated that the regularization does not act on in the objective space (it acts in the solution space) and that a smaller cost function leads to a lesser optimal regularization parameter and gives a better solution. These links have been seen in a very simple calculation example outside all physical meaning. However no demonstration has yet been found.

The question of error in the objective with an exact model has not been addressed here, except through the guaranteed quality. Nevertheless angle exists for an exact model and a perturbed objective. However, the aim here, not worked on yet, is to adapt the pressure to reach the model, not the opposite.

Now it must be underlined that the adaptation of the model has been carried out via the reduction of angle from the exact objective. In the presence of perturbation on the objective, the same procedure would be followed by reducing angle without being able to distinguish from . If the adaptation works well, the model will reach the perturbed objective and thus be far from the exact objective by angle . If , the cure is worse than the disease! How should the problem, clearly defined in geometrical terms, be tackled?

Another question concerns the geometrical representation of the best location of the sensors and their robustness. It was at the center of preoccupations in active control in the 2000's, as the paper by Baek and Elliott can testify. In this field, Martin and Gronier (2001) showed, first, that it is always possible to obtain the nominal quality of the inverse process by filtering a small number of components of the reference objective and, second, that the robustness of the filtered sensor configurations was not identical. It could be that the first property may be understood thanks to geometrical representation in the objective space (moreover with an understanding of one of the constraints associated with the assertion) and that the robustness may be understood by looking in the direction of the projection of the cone on the subspace to which the filtered components of the objective belong.

Acknowledgements

The author thanks J.R. Arruda, Professor at the University of Campinas, Brazil, and R. Sampaio, Professor at the Pontifical Catholic University, Rio de Janeiro, Brazil, for the interest they showed to the present work.

Paper received 21 April 2012.

Paper accepted 23 August 2012.

  • Augusztinovicz, F., Tournour, M., 1999, "Reconstruction of source strength distribution by inversing the boundary element method" in O. Von Estorff (Ed), Boundary Elements in Acoustics: Advances and Applications, The Bristish Library, The world's knowledge, pp. 243-284 (Chapter 8).
  • Baek, K.H., Elliott, S.J., 2000, "The effect of plant and disturbance uncertainties in active control systems on the placement of transducers", Journal of Sound and Vibration, 230(2), pp. 261-289.
  • Bai, M.R., 1992, "Application of BEM (boundary element method)-based acoustic holography to radiation analysis of sound sources with arbitrarily shaped geometries", Journal of the Acoustical Society of America, 92, pp. 199-209.
  • Berkhout, A.J., de Vries, D., Vogel, P., 1993, "Acoustic control by wave field synthesis", Journal of the Acoustical Society of America, 93(5), pp. 2764-2778.
  • Elliott, S.J., 2000, "Signal Processing for Active Control", Academic Press.
  • Gomes, J., Hansen, P.C., 2008, "A study on regularization parameter choice in nearfield acoustical holography", in the Proceedings of Euro Noise, Paris-France, pp. 895-900.
  • Hald, J., 2003, "Patch near-field acoustical holography using a new statistically optimal method", in the Proceedings of Internoise 2003, Korea, pp. 2203-22210.
  • Hald, J., 2009, "Theory and properties of statistically optimized near-field acoustical holography", Journal of the Acoustical Society of America, 125(4), pp. 2105-2120.
  • Hansen, P.C., 1998, "Rank-deficient and discrete ill-posed problems", SIAM.
  • Kim, Y., Nelson, P.A., 2003, "Spatial limits for the reconstruction of acoustic source strength by inverse methods", Journal of Sound and Vibration, 265, pp. 583-608.
  • Kirsh, A., 1996, "An introduction to the mathematical theory of inverse problems", Springer.
  • Le Bourdon, T., 2009, "Interprétation geometrique de l'holographie acoustique et usage pour l'adaptation du propagateur", Thèse UPMC.
  • Le Bourdon, T., Picard, C., Martin, V., 2009, "A parameter to classify different propagators for identifying acoustic source strengths", in the Proceedings of Novem 09 (Noise and Vibration: Emerging Methods), Oxford-UK.
  • Martin, V., Cariou, C., 1997, "Improvement of active attenuation by modifying the primary field", in the Proceedings of the Congress Active 97, Budapest, Hungary.
  • Martin, V., Gronier, C., 1998, "Minimum attenuation guaranteed by an active control system in presence of errors in the spatial distribution of the primary field", Journal of Sound and Vibration, 217(5), pp. 827-852.
  • Martin, V., Gronier, C., 2001, "Sensor configuration efficiency and robustness against spatial error in the primary field for active sound control", Journal of Sound and Vibration, 246(4), pp. 679-704.
  • Martin, V., Le Bourdon, T., 2010, "Indépendance de l'adaptation du propagateur et de la régularisation en holographie acoustique", in the Proceedings of the 10ième Congrès Français d'Acoustique, Lyon-F.
  • Martin, V., Le Bourdon, T., Pasqual, A.M., 2011, "Numerical simulation of acoustic holography with propagator adaptation; application to a 3D disc", Journal of Sound and Vibration, 330, pp. 4233-4249.
  • Martin, V., Le Bourdon, T., Arruda, R., 2012, "Geometrical interpretation of acoustic holography: adaptation of the propagator and minimum quality guaranteed in the presence of errors", Journal of Sound and Vibration, 331, pp. 3493-3508.
  • Maynard, J.D., Williams, E.G., Lee, Y., 1985, "Nearfield acoustic holography: 1. Theory of generalized holography and the development of NAH", Journal of the Acoustical Society of America, 78(4), pp. 1395-1412.
  • Metherel, A.F., El-Sum, H.M.A., Dreher, J.J., Larmore, L., 1967, "Introduction to acoustical holography", Journal of the Acoustical Society of America, 4, pp. 733-742.
  • Nelson, P.A, Elliott, S.J., 1992, "Active control of sound", Academic Press, London.
  • Nelson, P.A., 2001, "A review of some inverse problems for acoustics", International Journal of Acoustics and Vibration, 6(3), pp. 118-134.
  • Nelson, P.A., Yoon, S.H., 2000, "Estimation of acoustic source strength by inverse methods - part I: conditioning of the inverse problem", Journal of Sound and Vibration, 233(4), pp. 643-668.
  • Pasqual, A.M., Martin, V., 2012, "Optimal secondary source position in exterior spherical acoustical holophony", Journal of Sound and Vibration, 331, pp. 785-797.
  • Peleg, T., Rafaely, B., 2011, "Investigation of spherical loudspeaker array for local active control of sound", Journal of the Acoustical Society of America, 130(4), pp. 1926-1935.
  • Schuhmacher, A., Hald, J., Rasmussen, K.B., Hansen, P.C.., 2003, "Sound source reconstruction using inverse boundary element calculations", Journal of the Acoustical Society of America, 113(1), pp. 114-127.
  • Steiner, R., Hald, J., 2001, "Near-field acoustical holography without the errors and limitations caused by the use of spatial DFT", International Journal of Acoustics and Vibration, 6(2), pp. 83-89.
  • Williams, E.G., 2001, "Regularization methods for nearfield acoustical holography", Journal of the Acoustical Society of America, 110(4), pp. 1976-1988.
  • Wu, S.F., 2008, "Methods for reconstructing acoustic quantities based on acoustic pressure measurements", Journal of the Acoustical Society of America, 124(5), pp. 2680-2697.
  • Yoon, S.H., Nelson, P.A., 2000, "Estimation of acoustic source strength by inverse methods - part II: experimental investigation of methods for choosing regularization parameters", Journal of Sound and Vibration, 233(4), pp. 669-705.
  • The fundamental elements in certain inverse acoustic problems: their roles and interactions

    Vincent Martin
  • Publication Dates

    • Publication in this collection
      24 July 2013
    • Date of issue
      2012

    History

    • Received
      21 Apr 2012
    • Accepted
      23 Aug 2012
    Associação Brasileira de Engenharia e Ciências Mecânicas - ABCM Av. Rio Branco, 124 - 14. Andar, 20040-001 Rio de Janeiro RJ - Brazil, Tel.: +55 21 2221-0438, Fax: +55 21 2509-7129 - Rio de Janeiro - RJ - Brazil
    E-mail: abcm@abcm.org.br