In the real sound environment, the observation data are usually contaminated by additional external noise (i.e., background noise) of arbitrary distribution type. In order to estimate several evaluation quantities for specific signal, like (( ) percentile level), (averaged energy on decibel scale) and peak value, based on the observed noisy data, it is fundamental to estimate the momentarily fluctuating wave form of the specific signal.
Up to now, many methodological studies have been reported on the state estimation for stochastic systems    . However, many standard estimation methods proposed previously in a study of stochastic systems are restricted only to the Gaussian distribution   . Several state estimation methods for nonlinear system have been also proposed by assuming the Gaussian distribution of system and observation noises      . The real sound environment often shows an intricate fluctuation pattern rather than the standard Gaussian distribution. In our previous studies      , several state estimation methods for a sound environment system with non-Gaussian fluctuations have been proposed on the basis of expansion expressions for the probability distribution. Furthermore, state estimation methods for stochastic systems with complex characteristics and/or unknown structure have been proposed by using Bayes theorem on probability distribution    . Especially, though the unscented Kalman filter (UKF) and particle filter are useful for non-linear systems, UKF considers only the mean and variance of variables, and the particle filter needs very complicated algorithm based on Monte Carlo simulation   .
On the other hand, in the actual case contaminated by the background noise, some signal processing methods utilizing a digital computer are indispensable for estimating precisely the latent specific signal based on the noisy observations. Therefore, the observed data in an analogue form have to be transformed into a digital one at discrete time. However, many standard estimation methods proposed previously in stochastic systems are restricted only to a continuous level form of the observation.
In this study, a state estimation algorithm of an object-oriented type matched to the above real problem and the computer technique with successive observation for the complicated sound environment system is derived. First, in the real situation where the sound environment system is affected by background noise of arbitrary probability distribution, a stochastic process model with a quantized observation is established in Section 2. Next, by considering the linear and higher order correlation information between the specific signal and observation, two types of the recursive algorithms of Bayesian filter to estimate the unknown specific signal are theoretically proposed in the quantized level form matched to the signal processing by use of a digital computer. More specifically, in Section 3.1, an estimation algorithm in functional forms is derived by introducing the forward difference operator instead of usual differential operation. Furthermore, in Section 3.2, a method combining the analytical formula with Monte Carlo simulation is derived by introducing the idea of particles. Finally, the effectiveness of the proposed theory is experimentally confirmed too by applying it to the estimation problem of real sound environment in Section 4.
2. Formulation of Sound Environment System with Quantized Observation
Let us consider an arbitrary stochastic environment system with the power state variables of arbitrary distribution type, and express the system equation as:
where denotes the specific signal at a discrete time k, is the random input with known statistics. Here, and are statistically independent of each other. Two parameters F and G are estimated by using an auto-correlation technique  . Furthermore, the observation model is established by considering the additive property of power variables and the quantized observation in decibel scale, as follows:
where is the noisy observation contaminated by the additive background noise . Though is decibel variable with continuous level, the observation data are measured in a quantized level form suitable for the signal processing by use of a digital computer through A/D converter. The function denotes a nonlinear function expressing the quantization mechanism and is the quantized observation. Therefore, denotes a nonlinear function combining the nonlinearity of decibel observation with the quantized observation mechanism. In this study, a Bayesian filter to estimate the specific signal is proposed on the basis of the quantized observation contaminated by the background noise .
3. Establishment of Bayesian Filter with Quantized Observation
3.1. General Expression of Bayesian Filter in Expansion Series Form
In order to express explicitly the effect of successive observation on the estimated probability density function by use of various types of linear and/or nonlinear correlation between and , the well-known Bayes’ theorem is introduced.
where is a set of observation data up to time k. By expanding the conditional probability density function in a statistical orthogonal expansion series on the basis of the well-known standard probability distributions describing the dominant part of the actual fluctuation, the following expression is derived   .
The above two functions and are orthonormal polynomials of degrees m and n with weighting functions and . Based on Equation (5), the estimate of the polynomial function of with Mth order can be derived as follows.
where is an appropriate constant satisfying the following equality:
3.2. Estimation Algorithm by Introducing Difference Operation
In order to make the general theory for estimation algorithm more concrete, Gaussian distribution is considered as an example of standard probability functions for the specific signal:
Furthermore, as the fundamental probability function on the level-quantized observation, the generalized binomial distribution  with level difference interval can be chosen:
where is the minimum level of observations. The orthonormal polynomials with two weighting probability distributions in Equations (9) and (10) can be determined as
where denotes the Hermite polynomial with mth order, and is the jth order factorial function defined by 
Since the function in Equation (3) is not differentiable in general, the following expansion expression of discrete type is introduced with two arbitrary constants and .
where ∆ is the forward difference operator defined as:
After substituting Equations (2) and (3) into the definition of two parameters of and in Equation (10), by applying Equation (14), the following expressions can be derived (in the case of ).
Furthermore, the expansion coefficients defined by Equation (6) can be expressed as follows:
The above Equations (16)-(18) can be obtained from the statistics of the background noise and the predictions of at a discrete time ; i.e., the expectation values of arbitrary functions of conditioned by .
Therefore, the estimation algorithm in Equation (7) can be expressed explicitly as follows:
Especially, the estimates for mean and variance can be obtained as follows:
3.3. Estimation Algorithm by Introducing Particles
Though the particle filter is useful for the state estimation problem of non-linear systems, this filter needs very complicated algorithm and a large number of computational times based on Monte Carlo simulation and the resampling procedure  . In this section, a hybrid algorithm combining the analytical formula for state estimation with Monte Carlo simulation by use of particles is proposed.
The well-known Gaussian distribution is adopted as and , because this probability density function is the most standard one.
Then, the orthonormal functions with two weighting probability density functions in (23) can be given in the Hermite polynomial:
Accordingly, the estimation algorithm of the specific signal in Equation (7) can be given by
Furthermore, the estimates for mean and variance can be obtained as follows:
Thus, two parameters , and the expansion coefficients are expressed as follows:
The integrals in Equations (29)-(31) are evaluated by use of particles with mean , variance , higher order statistics and statistics of the background noise.
3.4. Prediction Algorithm
By considering Equation (1), the prediction step essential to perform the recurrence estimation can be given by
By replacing k with , the recurrence estimation can be achieved.
4. Application to Sound Environment
In order to examine the practical usefulness of the proposed Bayesian filter based on the quantized observation, the proposed method is applied to the actual sound environmental data. The road traffic noise is adopted as an example of a specific signal with a complex fluctuation form. Applying the proposed estimation method to actually observed data contaminated by background noise and quantized with 2 dB width roughly, the fluctuation wave form of the specific signal is estimated. The statistics of the specific signal and the background noise used in the experiment are shown in Table 1 and Table 2 respectively.
Figure 1 and Figure 2 show the estimation results of the fluctuation wave form of the specific signal for Data 1 and Data 2 by applying the algorithm proposed in Secttion 3.2. In this estimation, the finite number of expansion
Figure 1. Estimation results for Data 1 of the specific signal by applying the proposed method in Section 3.2 based on the quantized observation data with 2 dB width.
Table 1. Mean and standard deviation of the specific signal (in W/m2).
Table 2. Mean and standard deviation of the background noise (in W/m2).
Figure 2. Estimation results for Data 2 of the specific signal by applying the proposed method in Section 3.2 based on the quantized observation data with 2 dB width.
coefficients is used for the simplification of the estimation algorithm. In these figures, the horizontal axis shows the discrete time k, of the estimation process, and the vertical axis expresses the sound level taking a logarithmic transformation of power-scaled variables, because the actual sound environment usually is evaluated on dB scale. For comparison, the estimation results calculated by using the usual method are also shown in these figures. Since Kalman’s filtering theory is widely used in the field of stochastic system   , the extended Kalman filter is also applied to the observation data as a trail by introducing the following observation model.
where denotes the quantized noise and a uniform distribution within (q: the quantized width) is assumed as the probability distribution of . The results estimated by the proposed method show good agreement with the true values. On the other hand, there are great discrepancies between the estimates based on the standard type dynamical estimation method (i.e., extended Kalman filter), particularly in the estimation of the lower level values of the fluctuation.
Furthermore, the estimation algorithm proposed in Section 3.3 is applied to the observation data. The estimated results of two cases by applying the proposed algorithm to the quantized data with 2 dB width are shown in Figure 3 and Figure 4 respectively.
The squared sums of the estimation error are shown in Table 3. It can be found numerically that the proposed method is more useful than the extended Kalman filter.
In this study, state estimation method for sound environment system with
Figure 3. Estimation results for Data 1 of the specific signal by applying the proposed method in Section 3.3 based on the quantized observation data with 2 dB width.
Figure 4. Estimation results for Data 2 of the specific signal by applying the proposed method in Section 3.3 based on the quantized observation data with 2 dB width.
Table 3. Comparison between the proposed method and the extended Kalman filter for root-mean squared error of the estimation based on the quantized observation data with 2 dB width (in dB).
quantized level observation has been theoretically proposed on the basis of Bayes’ theorem. More specifically, two types of the recursive algorithm of Bayesian filter to estimate the specific unknown signal have been derived based on the quantized level observation matched for the signal processing by use of a digital computer. Furthermore, the validity and effectiveness of the proposed theory have been experimentally confirmed by applying it to the real environmental noise data in sound environment.
The proposed approach is still at the early of study, and there are left a number of practical problems to be continued in the future. For example, the proposed method has to be applied to many other actual data of sound environment. Furthermore, the proposed theory has to be extended to more complicated situations involving multi-signal sources, and an optimal number of expansion terms in the proposed estimation algorithm of expansion type has to be found.
The author is grateful for valuable suggestions at Inter-Noise 2016  . This work was supported in part by fund from the Grant-in-Aid for Scientific Research No.15K06116 from the Ministry of Education, Culture, Sports, Science and Technology-Japan.
*New type method is proposed based on the Baysian filter by use of the quantized observation in sound environment system.