Stochastic Signal Processing for Sound Environment System with Decibel Evaluation and Energy Observation

In real sound environment system, a specific signal shows various types of probability distribution, and the observation data are usually contaminated by external noise (e.g., background noise) of non-Gaussian distribution type. Furthermore, there potentially exist various nonlinear correlations in addition to the linear correlation between input and output time series. Consequently, often the system input and output relationship in the real phenomenon cannot be represented by a simple model using only the linear correlation and lower order statistics. In this study, complex sound environment systems difficult to analyze by using usual structural method are considered. By introducing an estimation method of the system parameters reflecting correlation information for conditional probability distribution under existence of the external noise, a prediction method of output response probability for sound environment systems is theoretically proposed in a suitable form for the additive property of energy variable and the evaluation in decibel scale. The effectiveness of the proposed stochastic signal processing method is experimentally confirmed by applying it to the observed data in sound environment systems.


Introduction
A specific signal in real sound environment system usually exhibits multifarious and complex characteristics such as non-Gaussian distribution and nonlinear property relating to natural, social, or human factors.Furthermore, the observation data usually are contaminated by external noise (e.g., background noise) with complex statistical properties.In this situation, in order to evaluate the sound environment system, precise estimation of the system characteristics of the sound environment is required by considering the contaminated observation data.
Furthermore, the internal physical mechanism of the real sound environment system is often difficult to recognize analytically, and it contains unknown structural characteristics.In our previous study, it was found that complex sound environment systems are difficult to analyze by using usual structural methods based on the physical mechanism [1].Therefore, a nonlinear system model was derived in the expansion series form reflecting various types of correlation information from the lower order to the higher order between input and output variables [2].The conditional probability density function contains the linear and nonlinear correlations in the expansion coefficients and these correlations play an important role as the statistical information for the input and output relationship of sound environment system.
On the other hand, in considering the relationship between the evaluation from top-down viewpoint and the countermeasure from bottom-up viewpoint in the sound environment system, noise evaluation quantities in decibel scale like   ((100 − ) percentile level) and  Aeq (averaged energy on decibel scale) and some countermeasure methods in energy scale are widely used.Since there is a certain scale transform between decibel and energy variables, a unified general consideration without losing their mutual relationship has to be derived.
In this study, a general type of complex sound environment systems is considered.A stochastic signal processing method for predicting the output response probability distribution in decibel scale based on the input observations is proposed for complex sound environment systems.
More specifically, an expansion expression of the conditional probability distribution in decibel scale is adopted as the system characteristics.Next, a method to estimate the system parameters reflecting several orders of correlation information between the input and output variables is derived by considering the additive property of energy variables under existence of external noise.Furthermore, a prediction method for the output probability distribution in decibel scale is also considered.
The effectiveness of the proposed theory is confirmed experimentally by applying it to real data of a sound insulation system and the road traffic noise environment measured around a national road in Hiroshima city.

Evaluation of Sound Environment System
under Existence of External Noise In order to find explicitly the various correlation properties between  and , let us expand the joint probability density function (, ) into an orthogonal polynomial series [3], as follows: (1)   ()  (2)   () ,   ≡ ⟨ (1)   ()  (2)   ()⟩ , where ⟨⟩ denotes the averaging operation. 0 () and  0 () can be chosen arbitrarily as the probability density functions describing the dominant parts of the actual fluctuation pattern.Two functions  (1)   () and  (2)   () are orthogonal polynomials with the weighting functions  0 () and  0 ().The information on the various types of linear and/or nonlinear correlations between  and  is reflected hierarchically in each expansion coefficient   .In this section, the Gaussian distribution suitable for the random variables in decibel scale is adopted as  0 () and  0 () with Thus, orthogonal polynomials  (1)   () and  (2)   () are given by the Hermite polynomial [3]: Substituting (1) into the definition of the conditional probability, ( | ) can be expressed in an expansion series form as follows: . (5)

Estimation of Correlation Information Based on Energy
Observation.In the measurement of the sound environment, the effects by external noise (e.g., background noise) are inevitable.Then, based on the additive property of energy variable, the observed sound intensity   at a discrete time  is expressed as where   and V  are sound intensities of the output signal for the sound environment system and external noise.We assume that the statistics of the external noise are known.In this section, an estimation method for the expansion coefficients   in (1), reflecting the correlation information between  and , is derived on the basis of the observed data   .There are relationships between energy variables   ,   and decibel variables   ,   for the input and output signals, as Therefore, from (7), the following relationships can be obtained with Next, considering the expansion coefficients   as unknown parameter vector a, a ≡ ( 1 ,  2 ,  3 , . ..) ≡ (a (1) , a (2) , . ..) , the simple dynamical model, is naturally introduced for the successive estimation of the parameter.
In order to derive the estimation algorithm of the parameter, attention is focused on Bayes' theorem for the conditional probability distribution: where   ≡ { 1 ,  2 , . . .,   } is a set of observation data up to time .Based on (12), using the similar calculation process to the previously reported paper [4], the estimate for an arbitrary polynomial function  M (a  ) of a  with Mth order can be derived as follows (cf.Appendix): with Two functions  (1)  m (a  ) and  (2)   (  ) are orthonormal polynomials with the weighting functions  0 (a  |  −1 ) and  0 (  |  −1 ).Furthermore,  Mm is the coefficient when the function  M (a  ) is expanded as (1)  m (a  ) .
As the concrete expression on the fundamental probability function for the parameter a  fluctuating in both positive and negative range, a standard Gaussian distribution is adopted.Furthermore, a gamma distribution is adopted as the probability function for the sound intensity   : with Therefore, the orthogonal polynomials with the weighting functions of ( 16) are given by Hermite polynomial and Laguerre polynomial [3]: By considering (6) and independence of   and V  , two parameters  *  and Ω  in (17) can be given by Considering ( 5) and the property of conditional expectation, the first terms of the right sides in the above equations are expressed as follows: The integrals in (20) can be calculated by using the relationship between energy and decibel variables in (8) and expansion expression in (1), as follows: where and  () 0 are coefficients satisfying the following relationship: By using the orthonormal condition of Hermite polynomial, the function   in (21) can be calculated as Therefore, (20) can be expressed as follows: with Φ (  ) ≡ ( (1)  0 (  ) ,  (1)  1 (  ) , . ..)  , A (), ≡ a * (), , ( = 1, 2, . ..) , A (0), ≡ ( 00 ,  10 ,  20 , . ..) , where  denotes the transpose of a matrix.Furthermore, using the definition of Laguerre polynomial and (1), the expansion coefficient  m can be calculated as follows: From ( 19) and ( 26)-( 28), it can be found that the parameters  *  , Ω  and the expansion coefficient  m are given by the predictions of unknown parameter a  , the statistics of the external noise V  , and the input observations   .
By considering (11), the prediction to perform the recurrence estimation can be given for an arbitrary polynomial function  N (a +1 ) with Nth order of a +1 can be expressed as

Prediction of Output Probability Distribution for Sound
Environment System.Because the conditional probability density function ( | ) can be considered as an invariant system characteristic, reflecting mainly the proper correlation relationship between the two sound pressure levels  and  in the sound environment system, the output probability distribution   () in decibel scale can be predicted, as   () = ⟨( | )⟩  .Thus, based on (1) and using the estimated parameter â ≡ ( Â11 , Â12 , . ..), the output probability density function   () can be predicted from the observed input data , as follows: 2)   () .(30)

Application to Real Sound Environment System
3.1.Application to Sound Insulation System.In order to confirm the effectiveness of the proposed method, it was applied to real data observed in a sound environment system.Acoustic signals observed by two microphones in indoors and outdoors for a house were adopted as input and output data for the sound insulation system.The schematic drawing of the sound environment system is shown in Figure 1.The rock music was selected as an input signal by considering the aggravation of "Karaoke" noise pollution problem, and white noise was adopted as a background noise.The statistics of the input, output signals and the background noise are shown in Table 1.The input and output fluctuation data simultaneously measured with every sampling interval 1 s.Based on the 500 data, the expansion coefficients   in (1) were estimated on the basis of the input signal   and observation   under existence of the background noise.
Based on the estimated expansion coefficients, the output response probability distribution excited by an arbitrary input signal was predicted.The 200 sampled data following the data used for the evaluation of expansion coefficients were adopted for predicting the output response probability distribution.Figure 2 shows the comparison between theoretically predicted curves and experimentally sampled points on the output probability distribution.The cumulative distributions of sound level are shown in this figure .The "theoretical curves" in this figure were obtained by predicting the probability density function of the output level  based on the observed data of the input level  by use of the theoretical expression.In this figure, "1st-3rd approximations" considered the expansion coefficients Â11 , Â12 , Â21 , and Â22 successively in (30).The "experimental values" represent the frequency distributions obtained directly from the observed data of the output level .
For comparison, the prediction results of the output probability distribution by introducing the standard regression     The regression coefficients   ,   ,   , and   in (31) were determined by applying the extended Kalman filter [5] after introducing the following observation equation: The theoretically predicted curves based on the proposed method show better agreement with the experimentally sampled values than the results by applying the extended Kalman filter based on the standard regression models.
From the cumulative distributions in Figures 2 and 3, noise evaluation quantities   ((100 − ) percentile level) can be evaluated.Several noise evaluation quantities obtained from these figures are shown in Table 2.It is obvious that the proposed method provides a more accurate prediction than the results based on the standard regression models.

Application to Road Traffic Noise Environment.
The effectiveness of the proposed method was confirmed experimentally by applying it to real road traffic noise data observed in the complicated sound environment near a national road.In order to evaluate the sound environment around the main line, the sound level at an evaluation point has to be predicted on the basis of the observation at a reference point.After regarding the sound levels at a reference point and an evaluation point as system input  and output , respectively, the probability distribution in decibel scale at the evaluation point connected with several evaluation quantities of the sound environment was predicted on the basis of the observation at the reference point.The reference point and the evaluation point were chosen at the positions being 1 m and 25 m apart from one side of the road as shown in Figure 4.The statistics of the input, output signals and the background noise are shown in Table 3.By applying the proposed method, the probability density function of the sound level at the evaluation point was predicted on the basis of the observation at the reference point.Road traffic noise was measured by the use of the sound level meter at every 0.2 s.Through the same procedure in Section 3.1, the expansion coefficients   in (1) were first estimated.
Based on the estimates of the expansion coefficients, the probability density function of road traffic noise at the evaluation point was predicted by measuring the road traffic noise data at the reference point.The predicted results are shown in Figures 5 and 6.Several noise evaluation quantities evaluated from Figures 5 and 6 are shown in Table 4. From these results, it can be seen that the theoretically predicted

Figure 2 :
Figure 2: Comparison between experimentally sampled values and theoretically predicted curves by the proposed method on the output probability distribution for the sound insulation system.

Figure 3 :
Figure 3: Comparison between experimentally sampled values and theoretically predicted curves by the standard regression models on the output probability distribution for the sound insulation system.

Figure 5 :
Figure 5: Comparison between experimentally sampled values and theoretically predicted curves by the proposed method on the probability distribution at the evaluation point for the road traffic noise environment near a national road.

Table 1 :
Statistics of the input, output signals and the background noise in sound insulation system.

Table 2 :
Comparison between the experimental values and theoretically predicted values for several noise evaluation quantities in dB evaluated from Figures2 and 3.

Table 3 :
Statistics of the input, output signals and the background noise in road traffic noise environment.