Issue 
Acta Acust.
Volume 5, 2021



Article Number  32  
Number of page(s)  13  
Section  Structural Acoustics  
DOI  https://doi.org/10.1051/aacus/2021023  
Published online  29 July 2021 
Scientific Article
Application of the LMS algorithm to identify the surface velocity responsible for the radiated sound pressure
^{1}
Applied Acoustics, Chalmers University of Technology, 41296 Göteborg, Sweden
^{2}
HEAD acoustics GmbH, 52134 Herzogenrath, Germany
^{*} Corresponding author: Wolfgang.Kropp@chalmers.se
Received:
11
March
2021
Accepted:
10
June
2021
A time domain approach based on the least mean square (LMS) algorithm is applied to reconstruct the source amplitude and source distribution on a plate. For this a numerical experiment is established. A boundary element model is used to calculate the required impulse response functions describing the pressure in near and farfield for a given volume flow at individual patches on the plate. Three different cases are considered. Firstly, a volume flow is given to a single patch. The LMS algorithm is used to reconstruct the source signal by means of receiving positions in the farfield. Secondly, the approach is used to identify the vibration pattern and source signal on a line of patches. Thirdly, a vibration pattern was given to the plate as the whole. For the reconstruction an assumption was made about the underlying vibration patterns (e.g. expansion in vibrational modes). Such an approach proved to be very time efficient and powerful. It also showed the need to place the receiving positions in the nearfield to be able to obtain correct results over the whole frequency range. However, this is not a problem of the approach based on the LMS algorithm, but just due to the underlying physics. It is not possible to deduce the nearfield from farfield observations, and therefore the vibrations mainly leading to a nearfield are simply not visible in the farfield.
Key words: Source identification / Inverse problem / Time domain
© W. Kropp et al., Published by EDP Sciences, 2021
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
For the successful design of sound and vibration properties of products and environments, experimental work will be indispensable even in the foreseeable future. One of the main tools in this context is the socalled transfer path analysis (TPA). Classical TPA dates back to the early eighties when Verheij measured indirectly the forces at the interface between vibration isolators, engine and ship hull (see e.g. [1]). Since the early eighties different approaches have been made to develop the TPA towards a tool applicable in industry. A multitude of different methods have emerged whose similarities, but also differences can easily be perceived as confusing. In [2] an attempt is made to classify the existing methods as “classical TPA”, “componentbased TPA” and “transmissibilitybased TPA” where the latter also includes the operational transfer path analysis. The authors in [2] conclude that the three families are highly interrelated and that little has been suggested about the methods’ relative success for typical industrial cases and sensitivity to various sources of error. Further research, in both numerical simulation and experimental validation, would substantiate this TPA framework in that respect.
In the special case of airborne sound source characterisation (i.e. reconstructing the velocity distribution on radiating surfaces) different methods have been developed over the recent years. Methods such as the nearfield acoustic holograpy (NAH) [3, 4], the inverse boundary element method (IBEM) (e.g. [5] and the least squared methods (LSM) (see e.g. [6, 7]) are used in this context. The NAH demands special geometries to allow the application of the Fourier transform. The IBEM is more general in geometry but is sensitive to similar errors as the TPA since an inversion has to be carried out.
Typical sources of errors are:
Insufficient signal to noise ratio which might lead to erroneous observation of field quantities (e.g. measured vibrations or sound pressures) or frequency response functions (FRFs). Measurement noise in some elements of e.g. a transfer matrix is distributed over all elements after inverting the matrix.
Nonaccurate FRFs due to e.g. small misplacement of transducers or the inherent difficulties to measure certain degrees of freedom (DoF) such as rotational degrees.
Not correct or not complete description of the interfaces as part of the transfer paths between components. Neglecting transfer paths can result in an erroneous result for the transfer paths considered.
Neglecting the presence of correlated multiple sources. Multiple correlated sources might lead to “crosstalking” at interfaces. Neglecting some of the sources will give erroneous results for those sources taken into account.
Performing the operational test in an environment different from the targeted assembly. This only leads to correct results when the mounting condition does not influence the source itself. E.g. the blockedforce [8] and freevelocity [9] methods provide two extreme boundary conditions for the description of an active component. The boundary condition in the case of the hybrid interface method [10] depends on the dynamic stiffness of the test rig. The different boundary conditions will influence the active component differently.
And finally just the limitation by lack of information due to the physics of the problem might not allow for identifying sources or transferpaths from measurements in the farfield accurately.
The scope of this paper is to demonstrate the capability of the LMS algorithm to identify the surface velocity responsible for the radiated sound pressure and to investigate the influence of different parameters such as background noise or accuracy of transfer function as well as limitations by physics on the quality of the results. In practical applications typically measurements of sound pressure signals in a multitude of receiving positions are carried out. At the same time a set of transfer functions from the anticipated source positions to the receiver positions is measured. In this way an inverse problem is established which has to be solved. Among the solving methods discussed in literature the overwhelming majority is formulated in the frequency domain and involve matrix inversion of different types. The inverse problem is well known to be illposed and has been shown to be sensitive to measurement noise and to suffer from numerical illconditioning at frequencies associated with the natural frequencies of the components [11]. Usually, the robustness of such solutions has to be improved by applying some form of regularization [12]. One usual attempt to cure problems is to create an overdetermined equation system and use truncated singular value decomposition [13].
In the following a different approach is utilised based on work originally presented in [14]. There two of the authors presented an approach based on the LMS algorithm in the time domain. They demonstrated the approach which is explained in more detail in Section 1.2 for a single input multiple output system (SIMO) where a force is exciting a structure in a known excitation position. The approach showed to be very robust when reconstructing forces acting on a beam. Although not included in the paper [14] but presented at the conference they also extended the approach to a multiple input multiple output system (MIMO system), i.e. excitation with several forces in known excitation positions. This extension was later published in [15] and applied in [16] for the indirect measurement of walking forces on wooden floors.
1.1 The setup
To investigate the possibilities to reconstruct the velocity distribution on a surface a numerical experiment is set up consisting of a 4 cm thick plate with the length of 0.78 m and the width of 0.3 m as shown in Figure 1.
Figure 1 Geometry of the plate used in the numerical experiments. 
The plate is situated at a height of 30 cm over an infinite and rigid ground. The plate surfaces are divided into quadratic shaped elements of the size of 3 cm times 3 cm. At the nodes of each element a velocity in normal direction of the surface is prescribed which determines the volume flow created by each element.
To calculate the radiation from the source distribution on the plate a half space BEM formulation is used based on the work by Brick presented in [17]. The sound pressure is evaluated at 321 receiving positions on a half sphere with a radius of 1 m around the origin and 260 positions situated directly 6 cm in front of each of the elements on the top surface of the plate (see Fig. 2).
Figure 2 The plate above the ground in combination and two sets of receiving positions. One set directly above the plate and one set distributed on a sphere. 
From the radiated pressure transfer functions between each element on the top surface and all receiving positions are calculated by normalising the calculated pressure values by the volume flow at the surface. Typical results for the transfer functions are shown in Figure 3.
Figure 3 Example for transfer functions for the radiation of one patch on the plate to a position directly above the plate (solid line) and on the half sphere (dashed line). 
The transfer functions to the positions close to the plate are relatively smooth over frequency and differ mainly in level. For positions on the sphere interference effects with strong dips in the FRF can be observed. From these FRFs the impulse response functions (IRFs) are calculated by inverse Fourier transform. In the following a maximum frequency range of 1500 Hz (i.e. a sampling rate of 3000 Hz) is considered. With 200 frequency lines the frequency resolution is 7.5 Hz.
1.2 The LMS algorithm applied to source identification
Today the LMS algorithm is a standard tool in filter design. It has been invented by Widrow and Hoff [18] and was the basis for their work with neural networks and machine learning. By means of the LMS algorithm a desired filter is approximated by determining a limited number of coefficients of a new filter that results in the least mean square of the error signal, i.e. the difference between a desired and the actual signal observed. At each time step the LMS algorithm follows the steepest gradient in order to minimise this error. The result converges to the optimal filter, the socalled Wiener–Hopf filter. Figure 4 shows a typical situation where the output a_{d}(n) of a system described by the IRF h_{o} should be reconstructed by a filter h with I coefficients.
Figure 4 Typical task for filter design with x as input signal, h_{o} the desired impulse response function of the filter, h the designed filter and e the error between desired signal a_{d} and the signal created by the filter h. 
By comparing the output a_{d}(n) of the system with the output from the filter h the error e(n) at each time sample n can be calculated as,(1)where x is the observed input into the system (x = [x(n), x(n − 1), … , x(n – I + 1)]^{T}) with T indicating the transposed form of x and h the unknown filter (h = [h(0), h(1), … , h(I − 1)]).
The expectation value of the squared error is then,(2)
The main idea of the LMS algorithm is to minimise the mean value of the quadratic error. As this is a quadratic equation the minimum is found by following the gradient of Equation (2) with respect to the filter coefficients h_{i},(3)
Starting from this Widrow and Hoff showed that instead of calculating the expectation value at each time step the gradient can be calculated as e(i) x(n −i) which leads to an adaptive process where,(4)and that this iterative process converges towards the results obtained by (3). α is the step size following the gradient in the adaptive process for calculating the filter coefficients. To ensure stability Widrow suggested in [19] as a “rule of thumb” that the step size has to obey the condition,(5)where E[x(n)^{2}] is the mean square value of the reference signal.
The frequent use of the LMS algorithm for filter design in e.g. active noise control is based on the fact,
That it does not demand a matrix inversion and, therefore the solution becomes more robust in the case of low SNR and illconditioned systems.
That the gradient is determined by the absolute error e(n). This means that the error due to strong components in the signals (e.g. due to resonances) governs the gradient rather than weak components (e.g. anti resonances). This circumstance makes the algorithm very robust with respect to background noise.
And that the algorithm is extremely simple to implement.
These properties make the LMS algorithm also an interesting approach for the reconstruction of input signals into systems described by impulse response functions. As the convolution is a commutative operation the only minor modification is the exchange of x(n) and h, i.e. make h to the input and x(n) to the function to be updated as shown in Figure 5.
Figure 5 Adaptive process to determine the unknown excitation x(n). 
The equation corresponding to Equation (4) is therefore,(6)The main difference is that h(n) is of finite length and always used completely. This has two consequences.
Firstly, the reconstructed signal x(n) is erroneous as long as n is smaller than the length of h. Consequently, the length of the desired/observed signal a_{d}(n) has to be at least twice the number of coefficients of h to capture the response to an excitation completely at least once.
Secondly, the time iterative process described by Equation (6) might be too short to converge to a sufficient accurate result. Therefore the process is repeated several times using the result of x of the previous run as start values.
The formulation in (4) is for a single input single output system (SISO). It easily can be extended for a single input multiple output (SIMO) system which would be the case in an overdetermined system. For a system with R output signals, Equation (1) yields the error at each receiving position as,(7)
a_{d,r} is the measured signal at each receiving position r and h_{r} the IRF from the excitation position to the receiving positions r. The gradient can be determined in analogy to (4) and the iterative updating procedure yields,(8)
By this the mean gradient is calculated. From the equation it is obvious that always the strongest absolute error will determine the gradient most. Extending the problem from a SIMO system to a MIMO system (i.e. multiple input and multiple output system) is a further but straight forward step. The error in this case is calculated as before as the deviation between the measured signal a_{d,r} and the predicted signal. However, as several sources are acting the predicted system at each receiving position is a superposition of the responses due to all sources, i.e.,(9)
The iterative updating is done for each source separately,(10)
The implementation of Equations (6)–(10) is rather simple resulting in just a few lines of code which also is a strength of the LMS algorithm.
2 Numerical experiments
Numerical experiments allow under controlled conditions to investigate the potential to reconstruct a given source distribution on an object by means of the approach based on the LMS algorithm as described in the previous section. Based on the numerical experiments for the plate presented in Section 1.1 a systematic study is presented in the following starting from a single point source and extended to more complex source distributions. To limit the complexity for the beginning only patches along a line as indicated in Figure 6 are considered. The patches are numbered from 1 to 24. When prescribing a velocity to a patch it is defined in the four nodes around the patch. The velocity in the normal direction to the surface in the centre of the patch is calculated as average of the velocity at the four nodes. As a consequence even neighbour patches will have a velocity as they share nodes with the patch where the velocity is prescribed.
Figure 6 Position of possible sources in front of the plate. 
2.1 A single source
In order to demonstrate the functioning of the LMS algorithm the most simple case of only one source (at patch 14) and one receiver is studied first (a SISO system). The source signal consists of white noise, but any other signal could have been chosen. As receiving position one of the positions on the half sphere is selected. It is assumed that the source position is known. For this case a reconstruction of the source signal is carried out. The results of the original time signal of the source and the reconstructed time signal after 800 repetitions of the updating process are shown in Figure 7. The iteration process could be interrupted when the error e becomes smaller than a certain limit or when the change in the error between two subsequent iteration steps is smaller than a given value. In this paper a limit of 800 iteration steps was used for all cases.
Figure 7 Original source signal (dashed line) and reconstructed source signal (solid line). 
As the length of the IRFs is 400 samples the reconstructed values for the sources signal for the first 400 samples are erroneous. After 400 samples the agreement between original and reconstructed signal seems to be fine, although such judgements are very difficult to make by means of inspecting the time records. Instead the autospectra of both signals (only for valid samples, i.e. higher than 400) are shown in Figure 8.
Figure 8 The upper curve shows the comparison of the original source spectrum (dashed line) and reconstructed source spectrum (solid line). The lower curve shows the difference between original source spectrum and reconstructed source spectrum. 
From the figure the good agreement is visible over a wide frequency range which is not too surprising as the data used in the simulation are free from noise and the source position is known. However, at very low frequencies there is a substantial deviation. The reason is the low response in the FRF (see Fig. 3). As the gradient is determined by the absolute error small values in the FRF will lead to small values in the pressure signals. As weak pressure signals hardly influence the adaptive process such frequency ranges are not updated properly. This can be cured as usual when applying the LMS algorithm by only focusing on the frequency range where a good adaptation is needed. Figure 9 shows the resulting error when applying a low pass filter with a cutoff frequency of 200 Hz to the input signal, i.e. the IRFs.
Figure 9 Error in the reconstructed source spectrum for the original IRFs (dotted line), after applying a lowpass filter to the IRFs with a cut off frequency of 200 Hz (dashed line) and after integrating the IRFs (solid line). 
Now the signal below 200 Hz controls the adaptive process. Consequently the results are improved at low frequencies but are erroneous above 200 Hz. The low response of the FRFs at low frequencies is based on the underlying physics as the sound pressure radiated from volume sources is proportional to the volume acceleration, i.e. the time derivative of the volume flow . Carrying out a time integration of the IRFs compresses the dynamics and the FRFs become more flat. Applying the integrated IRFs leads to that the reconstruction will be successful over more or less the whole frequency range as shown in Figure 9. Derivation with respect to time applied to the reconstructed source signals will correct such a choice of IRFs.
In the following the integrated IRFs are used throughout the paper.
2.1.1 The influence of measurement noise
A usual problem when using the TPA approach is that the observations (i.e. the measured pressure signals) are disturbed by measurement noise. Therefore, the effect of measurement noise on the performance of the LMS algorithm is investigated in the following.
In typical TPA approaches the influence of measurement noise is reduced by formulating an overdetermined system (SIMO), i.e. more receiving positions than sources are used for the reconstruction of the source properties. In the case here the pressure signals in several receiving positions on the half sphere are used. To demonstrate the functioning of the LMS algorithm in connection with an overdetermined system, noise is added to the simulated pressure signals at the receiving positions. The amplitude of the noise sequences is chosen so that for each receiving position a prescribed signal to noise ratio (SNR) is obtained. The number of receiving positions (Nr) in the half sphere is chosen as Nr = 2^{n} with n = 0, 1, … , 8. Figure 10 shows the results from this study in the form of the error calculated from the ratio between the energy in the reconstructed source signal Q_{r} and the energy in the original source signal Q_{o},(11)
Figure 10 Error according to Equation (11) in the prediction as function of receiving positions and signal to noise ratio shown as boxplots. The results are compared with the error predicted by Equation (12) (solid lines). The upper curve corresponds to 1 receiving position. The consequent curves correspond to 2, 4, 8, 16, 32, 64, 128 and 256 (lowest curve) receiving points. 
The error will depend on the SNR, the number Nr of receiving positions, but also on the selection of the receiving positions. Therefore it is necessary to calculate the error for a number of realisations with different receiving positions. 64 realisations are used for each combination of SNR and Nr and the results for each combination are represented as boxplot in Figure 10. The median is indicated by a notch and the edges of the box represent the 25th and 75th percentiles respectively. The whiskers are extended to the most extreme data points the algorithm considers to be not outliers, and the outliers (light grey) are plotted individually.
From the figure it is visible that the increase of receiving positions in the reconstruction process is able to suppress the negative impact of measurement noise. For a SIMO system the gradient in the adaptive process is calculated as the average over the error in all receiving positions. As the added noise is uncorrelated, averaging will suppress the influence of measurement noise on the gradient. In addition, in the time iterative procedure errors in the gradient will lead that the algorithm will march into a somewhat erroneous direction. However, in the next time step the error in the gradient will be different and therefore also the error in the direction. Therefore as long as the error is random, uncorrelated to the input signal and not too big the adaptive process will in any case converge on a random path towards the correct value. However, it will never reach the correct value exactly. How close the value will come to the true value is depending on the SNR and the degree of overdetermination.
The results show that already for Nr = 16 receiving positions the error ΔL can be kept below 0.5 dB in the case of an SNR = 0 dB. In Figure 10 analytical results are plotted describing the influence of averaging on the error ΔL as function of the SNR,(12)
There is a good agreement between the simulation results and Equation (12), indicating that the use of the LMS algorithm in connection with an overdetermined system is actually equivalent to averaging out the noise contribution and that the number of receiving positions is equivalent to the number of ensembles averaged over in Equation (12).
These results concern the overall energy in the spectrum. In the case of individual frequency components, however, there can be considerable deviations. To evaluate the similarity between the original and reconstructed source spectrum, an approach similar to the Modal Assurance Criterion (MAC) [20] is applied. The MAC value is widely used in modal analysis and is a metric that calculates the degree of linear similarity between two mode shapes. Using the same mathematical definition of the MAC we define here the Spectral Cross Correlation (SCC), which describes the similarity between two complex spectra. This enables a very robust similarity estimate including the spectrum shape and its energy content. The SCC value varies between 0 and 1 and it is defined as the normalized crosscorrelation of the two spectra. It is defined in Equation (13) (where (·)^{H} is the Hermitiantranspose operator),(13)
In Figure 11 a similar tendency for the MAC can be seen, in which the more similar the spectra are in shape, the higher the number of receiver positions and the SNR.
Figure 11 SCC between the original spectrum and the reconstructed spectrum of the source patch according to Equation (13) in the prediction as function of receiving positions and signal to noise ratio shown as boxplots. The curve indicates the median values. The lowest curve corresponds to 1 receiving position. The consequent curves correspond to 2, 4, 8, 16, 32, 64, 128 and 256 (highest curve) receiving points. 
2.1.2 Wrongly anticipated source location
The application of the LMS algorithm demands a correct assumption of the position of the source in order to establish the appropriate IRFs. In typical real case applications the source positions might often not be known exactly. The consequence might be an erroneous reconstruction of the source strengths. This is investigated in the following. Therefore, it is assumed that patch 14 (see Fig. 6) is the source while in reality the source is placed on different positions. The position was varied between patch 13 to patch 9, which corresponds to a misplacement by 3 cm to 15 cm. As expected with increasing distance the reconstruction of the source strength became increasingly erroneous over a wide range of the frequency spectrum. However, below about 500 Hz the error stayed rather small. Only one receiving position was chosen. To investigate the influence of more receiving positions the biggest distance (i.e. 15 cm misalignment) was chosen and the number of receiving positions for the reconstruction was gradually increased (Fig. 12). As the number of considered receiving positions increased the error decreases clearly.
Figure 12 Level difference between prescribed pressure signals and pressure signal which is obtained for the reconstructed source signal assumed at erroneous position. The calculations are carried out for 1 (—), 2 (…), 4 (–), and 8 (.)receiving positions. 
However, even for the case of evaluating two receiving positions only, similar good results have been obtained as for eight positions in Figure 12 for frequencies below 800 Hz, depending on the selected receiving positions utilised for the reconstruction. Closer investigations showed that it is more important to select proper positions rather than many. However, having many positions increases the chance that among these there are sufficient many relevant positions for reconstructing the source strength correctly, even if the source position is assumed 15 cm apart from the correct position which corresponds to a positioning error of more than 60% of the wavelength at the highest frequency considered here. Figure 13 shows that the spread is rather small for different randomly selected sets of receiving positions as long as there are sufficient (in this case eight) positions.
Figure 13 Level difference between prescribed pressure signals and the pressure signal which is obtained for the reconstructed source signal assumed at an erroneous position. The results are presented for seven simulation runs with different eight randomly selected receiving positions. 
From the results we can conclude that even with assuming wrong source positions it is possible to reconstruct the source signal rather accurately. However, these are also bad news as it might implicate that it could be difficult to reconstruct the spatial source distributions correctly as the space information does not seem to play a sufficiently important role at lower frequencies. This will be investigated in the next section.
2.2 Identification of vibration patterns
Up to now only a SISO or a SIMO system has been investigated. In the following this is extended to a MIMO system by including all patches indicated in Figure 6 (i.e. marked with a red line) as sources. Two cases are investigated. In the first case the patches are treated as individual sources as part of a line array. For each individual source the strength is identified by the LMS algorithm.
In the second case the individual sources are collected to global vibration patterns along the line array. The source strength of each of these patterns is then determined by the LMS algorithm. This method is shown to be faster. However, it demands some preknowledge about the source distribution to chose these patterns in an appropriate way.
2.2.1 A line array of correlated pistons
The prescribed source signals at each patch are correlated but different in amplitude. As time signal a white noise sequence is used. The velocity v_{m} at each patch m is described as,(14)where A is the amplitude of the pattern and N_{p} the number of half wavelength fitting to the whole length of the strip. N_{p} is varied between 1 and 5. In this way (14) might represent the modal pattern for e.g. a simply supported beam.
For the reconstruction of the velocity fields 40 receiving positions on the half sphere are utilised. Figures 14 and 15 show typical results for N_{p} = 3 and N_{p} = 5. The simulations show good agreement above a certain frequency. Below this frequency the results are erroneous. With increasing N_{p} this frequency is increasing as shown in Figure 16. There the spatial averaged squared velocity for the prescribed velocity pattern and the reconstructed pattern are compared. The results are expressed as relative error in decibel for N_{p} = 1 to N_{p} =5. In the figure also the socalled critical frequencies f_{c} in accordance to Table 1 are indicated as small circles. The critical frequency f_{c} as used here is defined as the frequency where the wavelength in the surrounding fluid (i.e. air) is identical with the prescribed wavelength on the plate.
Figure 14 Comparison of the original source distribution along the 24 patches (upper) and the reconstructed source distribution (lower) as function of frequency for N_{p} = 3. Reconstruction based on farfield receiving positions. 
Figure 15 Comparison of the original source distribution along the 24 patches (upper) and the reconstructed source distribution (lower) as function of frequency for N_{p} = 5. Reconstruction based on farfield receiving positions. 
Figure 16 Error in the averaged velocity between original and reconstructed source distribution. Circles indicate the critical frequencies as described in Table 1. The numbering of the curves corresponds to N_{p}. 
Inspecting Figures 14 and 15 one can argue that below f_{c} only the part of the prescribed velocity field can be observed in the farfield which is responsible for the radiated sound field. The part of the vibration resulting in a nearfield mainly cannot be observed in the farfield positions. The reconstructed source distribution show strong contributions at both ends of the strips where the actual velocity is close to zero. However, we know from theory that for instance for a radiator in a rigid baffle the radiation takes place at the edges/ends of the radiator (see e.g. [21]) at frequencies far below the critical frequency. In this way the results might correctly describe a possible source distribution explaining the radiated field although is not identical with the prescribed velocity. This is further discussed in Section 3.
If these considerations are correct it should be possible to identify the prescribed velocity exactly by using receiving positions in the nearfield directly above the plate as shown in Figure 2 and indeed the so reconstructed field is more or less identical with the prescribed velocity distribution (see Figs. 17 and 18).
Figure 17 Comparison of the original source distribution along the 24 patches (upper) and the reconstructed source distribution (lower) as function of frequency for N_{p} = 3. Reconstruction based on nearfield receiving positions. 
Figure 18 Comparison of the original source distribution along the 24 patches (upper) and the reconstructed source distribution (lower) as function of frequency for N_{p} = 5. Reconstruction based on nearfield receiving positions. 
The only bigger deviations are in the nodal lines where the velocity values are very small and therefore the LMS approach might have difficulties to converge to the correct values. However, this is hardly visible in the figures and does not play an important role for the sound radiation as the values are very small. In any case one can conclude that the vibration distribution on the plate can be correctly reconstructed when using nearfield receiver positions, but this might not be possible at frequencies below the critical frequencies of the underlying vibration patterns when using farfield positions.
2.2.2 Identification of predefined source configurations
The alternative to the identification of the source distribution for individual patches is to assume certain patterns from the very beginning. Often the vibration distribution is given by modes on the structure. Having knowledge about the modes, one could directly identify the contribution of each mode to the vibration of the structure and to the radiated sound. Instead of structural modes of course any pattern can be used. In the case here the pattern given by (14) is utilised. For this N_{p} = 1 to N_{p} = 10 are considered. The impulse response functions for each of these patterns are calculated as,(15)where v_{m} is the prescribed velocity for a modal pattern according to (14) and h_{m,r} the impulse reaponse function between patch m and receiving position r.
The reconstruction of a prescribed vibration pattern is based on two different sets of 36 receiving positions (one set in the nearfield 10 cm above the plate and one set in the farfield on the half sphere). In the following the prescribed vibration consists of a combination of pattern according to (14) and the velocity at each patch can be written as,(16)where A(N_{p}) is the amplitude for each pattern.
The reconstruction by means of farfield positions shows the same phenomena as in the previous section. Below the critical frequency given by the wavelength of the pattern N_{p} the velocity distribution cannot be reconstructed correctly as only the part of the source distribution can be observed which is responsible for the radiated sound. Consequently, a possible source distribution is reconstructed which explains the radiated sound at the farfield receiving positions. Below f_{c} the area in the middle of the radiating strip is contributing much less than the borders. However, due to the choice of the assumed base functions for the vibration pattern (i.e. sinusoidal functions), the patches at both ends always show zero velocity after reconstruction as shown in Figure 19. This also might indicate that the choice of the base of vibration patterns is not well adapted to the radiation problem although it certainly is the most appropriate to explain the vibrations on the plate.
Figure 19 Comparison of the original source distribution along the 24 patches (upper) given by (16) and the reconstructed source distribution (lower) as function of frequency for N_{p} = 5. Reconstruction based on farfield receiving positions. 
In Figure 19 the amplitudes A(N_{p}) are set to zero for all N_{p} besides for N_{p} = 5 which has unit amplitude. This is identical with the prescribed vibration patterns used in Figure 15. However, the results differ from the results shown there. This demonstrates the importance to use for the radiation relevant base functions and does not depend on the procedure of reconstructing the patterns. In the case considered here the sinusoidal function are not adequate. Similar problems might be observed when using modal patterns as base functions. The approach in 2.2.1 might be in this context clearly superior as it does not assume any pattern in advance.
The situation is again different when using nearfield receiving positions for the reconstruction process. Results show that in the absence of measurement noise at the receiving positions the amplitudes A(N_{p}) of the underlying patterns can be determined perfectly. This is true if the time signals for all patterns are correlated or uncorrelated and if their amplitudes are in the same order. Even if one single pattern is very small in amplitude in comparison to all other prescribed amplitudes the reconstruction works well over a wide frequency range as shown in Figure 20.
Figure 20 Level difference between original source strength and reconstructed source strength for pattern N_{p} = 3 for three different amplitudes: .. 10^{−4},   10^{−5}, —10^{−6} and the amplitudes for the pattern N_{p} = 1, 2, 4, 5 set to one. 
Figure 20 shows the level difference between original source strength and reconstructed source strength for the pattern N_{p} = 3. The amplitude A(3) is varied between 10^{−4} and 10^{−6} while the amplitudes for the pattern N_{p} = 1, 2, 4, 5 are set to one.
In these calculations all patterns have the same underlying time signal, i.e. all patterns are correlated but different in amplitude. In the presence of measurement noise the potential to reconstruct such weak source signals will strongly depend on the SNR. For decreasing SNR the possibility to reconstruct low amplitudes gets successively lost as expected. Similar results are obtained when assuming uncorrelated signals for all individual patterns.
3 Application of the LMS algorithm to reconstruct a complex vibration pattern
As final demonstration of the functioning of the LMS algorithm in the context of source identification the whole plate has a prescribed velocity pattern. This pattern is assumed to consist of sinusoidal functions similar to Equation (16) as,(17)where x and y are the coordinates of the patch centre and n and l the number of half wavelengths fitting to the length and width respectively.
As amplitudes A(n, l) the values were chosen randomly with the same order of magnitude (see Table 2 below).
Using nearfield receiving positions allows to reconstruct the pattern on the plate rather well. However, to ensure a good agreement at low frequencies a high number of receiving positions is needed (see Fig. 21).
Figure 21 Comparison of the original pattern at 100 Hz (left upper corner) with the reconstructed pattern for 12 receiving position (right upper corner), 48 receiving position (left lower corner) and 64 receiving positions (right lower corner). 
The differences mainly occur in nodal lines where the amplitudes are small which is typical for the LMS algorithm as it is controlled by the absolute error. However, looking at the averaged squared sound pressure level on the sphere which is proportional to the radiated sound power, the agreement between original values and the values obtained from the reconstructed source distribution are very close even if the reconstruction is based only on 12 receiving positions (see Fig. 22).
Figure 22 Averaged squared pressure on the half sphere for the original vibration pattern (solid gray line) and the reconstructed vibration pattern (dashed black line). 
Using the farfield receiving positions for the reconstruction the results are very similar to the results in Section 2.2.
Figure 23 shows the results for the reconstruction based on 128 receiving positions for different frequencies. At low frequencies the pattern differs strongly from the original pattern shown in Figure 21. With increasing frequencies the differences decrease and at high frequencies the agreement is very good. However, looking at the averaged squared pressure on the sphere, the agreement is perfect over the whole frequency range. One can again conclude that the problem suffers from the fact that the same farfield pressure fields can be created by different source compositions. The problem is not the applied reconstruction method, but the underlying physics which makes it difficult to observe the required nearfield information when the farfield is dominating.
Figure 23 Results of the reconstructed vibration pattern at 100 Hz (left upper corner), 250 Hz (right upper corner), 500 Hz (left lower corner) and 1000 Hz (right lower corner). 
4 Conclusions
A time domain approach for identifying the source distribution on a plate based on the LMS algorithm has been used. The approach has been demonstrated for SISO, SIMO and MIMO systems. A numerical experiment has been created as test case. FRFs have been calculated by means of a BE model describing the pressure at receiving positions normalised by the volume flow at discrete patches on the plate. Two sets of receiving positions were used, one in the farfield in the form of a half sphere around the centre of the plate and one in the nearfield close to the surface of the plate. From the FRFs the IRFs were calculated as input to the LMS algorithm.
For the case of a single monopole represented by one of the patches on the plate the LMS algorithm demonstrated to be very robust. As for frequency domain approaches an overdetermined system (SIMO) can be formulated to suppress efficiently measurement noise as long as the noise is not correlated at the different receiving positions. A misspositioning of the source (i.e. the assumption of the wrong source position) showed to result in correct source signals as long as relevant receiving positions are utilised for the identification. This can be achieved by just choosing enough receiving positions. In the case considered here for a misspositioning of 15 cm the resulting error was still smaller than 1 dB. Two different sets of receiving positions were utilised for the reconstruction of the vibration pattern on the plate.
The first set consisted of receiving positions in the nearfield of the source on the plate. The examples in this paper showed that in this case accurate reconstruction is achievable as long as sufficient many observation positions are utilised. It was shown that the so reconstructed source distribution creates the same farfield sound field as the original distribution.
The second set consisted of receiving positions in the farfield of the sources on the plate. In this case the reconstruction leads to source distributions which at lower frequencies clearly differ from the original source distributions. This is due to the fact that in the farfield only the part of the source pattern can be observed that is responsible for the radiated sound. At the end it is not a failure of the approach, but a consequence of the underlying physics which determines the limits of the reconstruction as the required nearfield information is hardly to be extracted when the farfield is dominating.
Referring back to the list of typical sources of errors as mentioned in the introduction of this paper one can conclude the following for the topics considered in this paper:
Insufficient signal to noise ratio which might lead to erroneous observation of field quantities (e.g. measured vibrations or sound pressures) or frequency response functions (FRFs) can be improved by an overdetermined system. The LMS algorithm is extremely robust with respect to measurement noise due to the fact that it is controlled by the strongest signals which has the highest SNR. In addition as long as the measurement noise is not correlated to the signals to be reconstructed, even erroneous gradients in the updating process are random and averaged out over the iteration process.
Nonaccurate FRFs due to e.g. small misplacement of transducers or the inherent difficulties to measure certain degrees of freedom (DoF) such as rotational degrees of freedom. The LMS algorithm showed to be capable even for this case to find correct source signals as long as sufficiently relevant receiving position are used.
Neglecting the presence of correlated multiple sources. Multiple correlated sources might lead to crosstalking at interfaces. Neglecting some of the source will give erroneous results for those sources taken into account. This has turned out not to be a problem in the case considered here. Assuming several correlated sources did not create any problems for the identification process. Even correlated sources with substantially low amplitudes could be detected as long as the background noise is not too high and therefore masking the contribution of low amplitude sources.
In practical applications it will be required either to calculate or to measure a complete set of impulse response functions between the expected sources patches and or vibration pattern (e.g. modes) and observation positions. Calculation using e.g. the boundary element method will be very efficient but often demand simple situations (e.g. anechoic conditions) to keep the computational effort acceptable. Measurement in contrary can include the complexity of the sound field. It can be expected that the complexity of a sound field in a room will improve the efficiency of the method as it will increase the difference between individual impulse response functions as long as the reverberation will not become too high.
In some way the presented approach is similarto NAH or IBEM as it uses preknowledge about the sound field to reconstruct the source distribution. The big advantage of the approach, however, is that does not demand certain geometries and no inversion is required and therefore typical problems of inverse methods are avoided. In addition the time domain formulation is perfectly for transient signals and time varying source distributions.
Finally, there is the limitation by lack of information due to the physics of the problem might not allow for identifying sources or transfer paths accurately. In the case considered here, it has been demonstrated that the observation of the sound field in the farfield is not sufficient for reconstructing the source distribution if the sources mainly create a nearfield. In this case it is essential to place observation positions in the nearfield to be able to “see” the correct source distribution.
References
 W.J. Verheij: Multipath sound transfer from resiliently mounted shipboard machinery: Experimental methods for analyzing and improving noise control. Ph.D. Thesis at Delft University of Technology, The Netherlands, 1982. [Google Scholar]
 M. Van der Seijs, D. de Klerk, D.J. Rixen: General framework for transfer path analysis: History, theory and classification of techniques. Mechanical Systems and Signal Processing 68–69 (2016) 217–244. [Google Scholar]
 G. Williams, J.D. Maynard: Numerical evaluation of the Rayleigh integral for planar radiators using the FFT. Journal of the Acoustical Society of America 72 (1982) 2020–2030. [Google Scholar]
 J.D. Maynard, E.G. Williams, Y. Lee: Nearfield acoustic holography: I. Theory of generalized holography and the development of NAH. Journal of the Acoustical Society of America 78 (1985) 1395–1413. [Google Scholar]
 G.T. Kim, B.T. Lee: 3D sound source reconstruction and field reproduction using the Helmholtz integral equation. Journal of Sound and Vibration 136 (1990) 245–261. [Google Scholar]
 Y.C. Chao: An implicit leastsquares method for the inverse problem of acoustic radiation. Journal of the Acoustical Society of America 81 (1987) 1288–1292. [Google Scholar]
 Z. Wang, S.F. Wu: Helmholtz equationleastsquares method for reconstructing the acoustic pressure field. Journal of the Acoustical Society of America 102 (1997) 2020–2032. [Google Scholar]
 A.S. Elliott, A.T. Moorhouse: Characterisation of structure borne sound sources from measurement insitu. Journal of the Acoustical Society of America 123, 5 (2008) 3176. [Google Scholar]
 J.M. Mondot, B.A.T. Petersson: Characterization of structureborne sound sources: the source descriptor and the coupling function. Journal of Sound and Vibration 114, 3 (1987) 507–518. [Google Scholar]
 D. De Klerk, D.J. Rixen: Component transfer path analysis method with compensation for test bench dynamics. Mechanical Systems and Signal Processing 24, 6 (2010) 1693–1710. [Google Scholar]
 B.J. Dobson, E. Rider: A review of the indirect calculation of excitation forces from measured structural response data. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science 204–2 (1990) 69–75. [Google Scholar]
 A.N. Thite, D.J. Thompson: The quantification of structureborne transmission paths by inverse methods. Part 2: Use of regularization techniques. Journal of Sound and Vibration 264, 2 (2003) 433–451. [Google Scholar]
 A.N. Thite, D.J. Thompson: The quantification of structureborne transmission paths by inverse methods. Part 1: Improved singular value rejection methods. Journal of Sound and Vibration 264, 2 (2003) 411–431. [Google Scholar]
 W. Kropp, K. Larsson: Force estimation in the time domain by applying an LMSalgorithm, in Proceedings NOVEM (2005). [Google Scholar]
 M. Sturm, A.T. Moorhouse, W. Kropp, T. Alber: Robust calculation of simultaneous multichannel blocked force signatures from measurements made insitu using an adaptive algorithm in time domain, in Bangkok, Thailand, 2013. [Google Scholar]
 N. Amiryarahmadi, W. Kropp, K. Larsson: Application of the LMS algorithm to measure lowfrequency transient forces from human walking. Acta Acustica United with Acustica 202, 1 (2016) 23–34. [Google Scholar]
 H. Brick: Application of the Boundary Element Method to combustion noise and halfspace problems. Ph.D. Thesis, Chalmers University of Technology, Göteborg, Sweden, 2009. [Google Scholar]
 B. Widrow, M. Hoff: Adaptive switching circuits. Proceedings IRE Wescon Convention Record, Part 4 Session 16 (1960) 96–104. [Google Scholar]
 B. Widrow, J.M. Cool, M.G. Larimore, C.R. Johnson: Stationary and nonstationary learning characteristics of the LMS adaptive filter. Proceedings of the IEEE 64, 8 (1976) 1151–1162. [Google Scholar]
 R.J. Allemang: The modal assurance criterion – twenty years of use and abuse. Sound and Vibration 37, 1, (2003) 14–23. [Google Scholar]
 F. Fahy: Sound and structural vibration. Academic Press, 1985, p. 86. [Google Scholar]
Cite this article as: Kropp W. Larsson K. Lobato T. & Sottek R. 2021. Application of the LMS algorithm to identify the surface velocity responsible for the radiated sound pressure. Acta Acustica, 5, 32.
All Tables
All Figures
Figure 1 Geometry of the plate used in the numerical experiments. 

In the text 
Figure 2 The plate above the ground in combination and two sets of receiving positions. One set directly above the plate and one set distributed on a sphere. 

In the text 
Figure 3 Example for transfer functions for the radiation of one patch on the plate to a position directly above the plate (solid line) and on the half sphere (dashed line). 

In the text 
Figure 4 Typical task for filter design with x as input signal, h_{o} the desired impulse response function of the filter, h the designed filter and e the error between desired signal a_{d} and the signal created by the filter h. 

In the text 
Figure 5 Adaptive process to determine the unknown excitation x(n). 

In the text 
Figure 6 Position of possible sources in front of the plate. 

In the text 
Figure 7 Original source signal (dashed line) and reconstructed source signal (solid line). 

In the text 
Figure 8 The upper curve shows the comparison of the original source spectrum (dashed line) and reconstructed source spectrum (solid line). The lower curve shows the difference between original source spectrum and reconstructed source spectrum. 

In the text 
Figure 9 Error in the reconstructed source spectrum for the original IRFs (dotted line), after applying a lowpass filter to the IRFs with a cut off frequency of 200 Hz (dashed line) and after integrating the IRFs (solid line). 

In the text 
Figure 10 Error according to Equation (11) in the prediction as function of receiving positions and signal to noise ratio shown as boxplots. The results are compared with the error predicted by Equation (12) (solid lines). The upper curve corresponds to 1 receiving position. The consequent curves correspond to 2, 4, 8, 16, 32, 64, 128 and 256 (lowest curve) receiving points. 

In the text 
Figure 11 SCC between the original spectrum and the reconstructed spectrum of the source patch according to Equation (13) in the prediction as function of receiving positions and signal to noise ratio shown as boxplots. The curve indicates the median values. The lowest curve corresponds to 1 receiving position. The consequent curves correspond to 2, 4, 8, 16, 32, 64, 128 and 256 (highest curve) receiving points. 

In the text 
Figure 12 Level difference between prescribed pressure signals and pressure signal which is obtained for the reconstructed source signal assumed at erroneous position. The calculations are carried out for 1 (—), 2 (…), 4 (–), and 8 (.)receiving positions. 

In the text 
Figure 13 Level difference between prescribed pressure signals and the pressure signal which is obtained for the reconstructed source signal assumed at an erroneous position. The results are presented for seven simulation runs with different eight randomly selected receiving positions. 

In the text 
Figure 14 Comparison of the original source distribution along the 24 patches (upper) and the reconstructed source distribution (lower) as function of frequency for N_{p} = 3. Reconstruction based on farfield receiving positions. 

In the text 
Figure 15 Comparison of the original source distribution along the 24 patches (upper) and the reconstructed source distribution (lower) as function of frequency for N_{p} = 5. Reconstruction based on farfield receiving positions. 

In the text 
Figure 16 Error in the averaged velocity between original and reconstructed source distribution. Circles indicate the critical frequencies as described in Table 1. The numbering of the curves corresponds to N_{p}. 

In the text 
Figure 17 Comparison of the original source distribution along the 24 patches (upper) and the reconstructed source distribution (lower) as function of frequency for N_{p} = 3. Reconstruction based on nearfield receiving positions. 

In the text 
Figure 18 Comparison of the original source distribution along the 24 patches (upper) and the reconstructed source distribution (lower) as function of frequency for N_{p} = 5. Reconstruction based on nearfield receiving positions. 

In the text 
Figure 19 Comparison of the original source distribution along the 24 patches (upper) given by (16) and the reconstructed source distribution (lower) as function of frequency for N_{p} = 5. Reconstruction based on farfield receiving positions. 

In the text 
Figure 20 Level difference between original source strength and reconstructed source strength for pattern N_{p} = 3 for three different amplitudes: .. 10^{−4},   10^{−5}, —10^{−6} and the amplitudes for the pattern N_{p} = 1, 2, 4, 5 set to one. 

In the text 
Figure 21 Comparison of the original pattern at 100 Hz (left upper corner) with the reconstructed pattern for 12 receiving position (right upper corner), 48 receiving position (left lower corner) and 64 receiving positions (right lower corner). 

In the text 
Figure 22 Averaged squared pressure on the half sphere for the original vibration pattern (solid gray line) and the reconstructed vibration pattern (dashed black line). 

In the text 
Figure 23 Results of the reconstructed vibration pattern at 100 Hz (left upper corner), 250 Hz (right upper corner), 500 Hz (left lower corner) and 1000 Hz (right lower corner). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.