State Estimation for Fuzzy Sound Environment System with Finite Amplitude Fluctuation *

Internal physical mechanism of actual sound environment system is often difficult to recognize analytically, and it contains unknown structural characteristics. Furthermore, the observation data often contain fuzziness due to several causes and exhibit level saturation owing to the existence of a finite dynamic range. Therefore, it is necessary to propose a new state estimation method by considering fuzziness and finite amplitude fluctuation of observation data. In this paper, a method for estimating the specific signal for sound environment system with unknown structure is proposed in an appropriate form for the finite level range of the measured fuzzy observation data by introducing an expansion expression of probability distribution with Bata distribution in the first term and new type of membership function. The effectiveness of the proposed theoretical method is confirmed by applying it to the actual problem in the sound environment.


Introduction
The internal physical mechanism of actual sound environment system having a complicated relation to various factors is often difficult to recognize analytically, and it contains unknown structure.Furthermore, the stochastic process observed in the actual phenomenon exhibits complex fluctuation pattern and there are potentially various nonlinear correlations in addition to the linear correlation between input and output time series [1].
Furthermore, it is necessary to pay our attention on the fact that the observation data in the sound environment system often contain fuzziness due to several causes, for example, the permissible error of the accuracy in measurements, the quantized error in the digitization of observation data, and the existence of confidence limitation in measuring instruments.In our previous studies, state estimation methods for sound environment system based on fuzzy observation have been proposed by considering the standard Gaussian type membership function [2] [3].On the other hand, owing to a finite range of the amplitude fluctuation and/or a definite dynamic range of measurement instrument, only the information defined within a finite range usually can be obtained as reliable data in the measurement of real sound environment [4] [5].In this situation, in order to evaluate more precisely the objective sound environment system, it is necessary to introduce a fuzzy theory to estimate the waveform fluctuation of the specific signal based on the observed data with fuzziness and amplitude saturation.
As a typical method in the state estimation problem, the Kalman filtering theory and its extended filter are well known [6] [7] [8] [9] [10].Furthermore, several state estimation methods for stochastic system with non-Gaussian fluctuations have previously been proposed [11] [12] [13] [14].These theories are originally based on the additive model of the specific signal and an external noise.The actual sound environment systems often contain unknown characteristics in the relationship between the state variable and the observation.Furthermore, the observation data often contain fuzziness and amplitude saturation.
In this study, based on the fuzzy observations, a Bayesian filter for estimating the specific signal of sound environment systems with unknown structural characteristic is theoretically proposed in an appropriate form for the finite amplitude range of the measured data.More specifically, complex sound environment systems which have to be treated as the systems with unknown characteristics are paid our attention.By introducing the orthogonal expansion expression of the probability distribution with a Beta distribution as the first term, and a new type of membership function, which are suitable for the finite amplitude fluctuation range of the signal and observation, a method to estimate the waveform fluctuation of the specific signal based on the fuzzy observation data is proposed.
After adopting a previously reported expansion expression of the conditional probability distribution as the correlation information between the specific signal and observation [1], by regarding the expansion coefficients as unknown parameters and introducing probability measure of fuzzy events [15], these parameters are estimated simultaneously with the specific signal based on the fuzzy observation.
Finally, the effectiveness of the proposed Bayesian filter focusing on the relationship between variables is confirmed experimentally too by applying the theory to the estimation of sound level based on the observation data containing fuzziness and amplitude saturation.

Formulation of Sound Environment System with Fuzzy Observation
Let k x and k y be the input and output signals at a discrete time k for a sound environment system.It is assumed that there are complex nonlinear relationships between k x and k y , which are difficult to find a fundamental relation- ship between them.Since the system characteristics are unknown, a system model in the form of a conditional probability is adopted.More precisely, attention is focused on the joint probability distribution function ( ) These orthonormal polynomials can be decomposed by using Schmidt's orthogonalization [16].From Equation (1), the conditional probability distribution function ( ) Though Equation ( 5) is originally infinite series expansion, finite expansion series with r R ≤ and s S ≤ is adopted because only finite expansion coefficients are available and the consideration of the expansion coefficients from the first few terms is usually sufficient in practice.The expansion coefficients defined by Equation (2) satisfy the following conditions: ( ) Furthermore, let k z be the fuzzy observation obtained from the output k y with a finite fluctuation range [ ] Since the objective system contains an unknown specific signal and unknown structure, the expansion coefficients rs A expressing hierarchically the correla- tion relationship between k x and k y must be estimated on the basis of the fuzzy observation k z .Considering the expansion coefficients rs A as unknown parameter vector a : , , , , , , , , , , , the following simple dynamical models are introduced for the simultaneous estimation of the parameters with the specific signal k x : where ( ) is the number of unknown expansion coefficients to be estimated.
On the other hand, the following time transition model for the input signal is generally established.
where k u is the random input with mean 0 and variance 2 u σ .Two parameters F and G are estimated by using an auto-correlation technique [16].
A method to estimate k x adaptively based on the fuzzy observation k z is derived in this study by introducing probability measure of fuzzy events [15] and the expansion series expressions of the conditional probability distribution function in Equation ( 5).

Derivation of Estimation Algorithm Based on Fuzzy Observation with Amplitude Saturation
In order to derive an estimation algorithm for a specific signal k x , based on the successive observations of fuzzy data k z , we focus our attention on Bayes' theorem for the conditional probability distribution [12].Since the parameter k a is also unknown, the conditional probability distribution of k x and k a is considered.
, , , is a set of fuzzy observation data up to a time k.
After applying probability measure of fuzzy events [15] to the right side of Equation (10), expanding it in a general form of the statistical orthogonal expansion series [12], the conditional probability density function ( ) with , , , .
The functions ( ) ( ) ϕ m a and ( ) ( ) are the orthogonal polynomials of degrees l , m and n with weighting functions ( ) | k k P y Z − , which can be artificially chosen as the proba- bility density functions describing the dominant parts of ( ) P y Z − .These three functions must satisfy the following orthonormal relationships: ) ( ) ( ) ( ) ( ) ( ) Based on Equation (11), the recurrence algorithm for estimating an arbitrary ( ) x and k a can be derived as follows: where is the expansion coefficient determined by the equality: In order to make the general theory for estimation algorithm more concrete, the well-known Gaussian distribution is adopted as ( ) Furthermore, the Bata distribution [17] suitable for the random variable with a finite fluctuation range [ ] , a b is adopted as ( ) where ( ) Γ ⋅ is the Gamma function.Then, the orthonormal functions with the weighting probability density functions in Equations ( 19) and ( 21) can be given by H ⋅ is the Hermite polynomial with lth order [18] and Jacobi polynomial with nth order [4].As the membership function, the following function matching to the Beta distribution is newly introduced.
Accordingly, Equation ( 12) can be given by The fuzzy data k z are reflected in , ; By considering the orthonormal condition of Jacobi polynomial [4], Equation ( 27) can be expressed as follows: where a few concrete expressions of 0 n d in Equation (30) can be expressed as follows: In two special cases when ( ) , , , estimates related to mean and variance of the specific signal are expressed as follows: { } ( ) ( ) Using the property of conditional expectation and Equation ( 5), the two variables * k y and k Ω in Equation ( 22) can be expressed in functional forms on predictions of k x and k a at a discrete time 1 k − (i.e. the expectation value of arbitrary functions of k x and k a conditioned by where T denotes the transpose of a matrix.The coefficients 1s e and 2s e in Equations ( 35) and ( 36) are determined in advance by expanding k y and ( ) in the following orthogonal series forms: Furthermore, using Equation ( 5) and the orthonormal condition of Equation ( 4), each expansion coefficient l n B m defined by Equation ( 13) can be obtained through the similar calculation process to Equations ( 35) and ( 36), as follows: In the above, the expansion coefficient l n B m can be given by the predictions of k x and k a .
Finally, by considering Equation ( 9), the prediction step which is essential to perform the recurrence estimation can be given by By replacing k with 1 k + , the recurrence estimation of k x can be achieved.

Experimental Consideration
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 shown by the frame in Figure 1.
The rock music was selected as an input signal by considering the aggravation of "Karaoke" noise pollution problem.
Figure 1.A schematic drawing of the experimental set up in sound insulation system.
After generating the music sound inside the house, the indoor and outdoor sound pressure levels were regarded as the input signal k x and the output ob- servation k y .The data were measured with a sampling interval of 1 s.The sta- tistics of the input signal and the output signal used in the experiment are shown in Table 1 and Table 2. Applying the proposed estimation method to actually observed output data quantized roughly with 1 dB, 2 dB and 3 dB widths as examples of fuzzy observation, the fluctuation wave form of the input signal was estimated.Comparison of the accuracy for the estimated probability distribution between the case considering a finite fluctuation range of the observation data and the case without consideration of the fluctuation range was discussed in our previous study [5].are 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 pressure level.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 [8] is also applied to the observation data as a trail by introducing the following observation model.
where k ε denotes the quantized noise.A uniform distribution within [ ] 2, 2 q q − (q: the quantized width) is assumed as the probability distribution of k ε .Since the parameters k a and k b are unknown, these parameters are also estimated simultaneously with the input signal k x by introducing the fol- lowing time transition models in addition to Equation ( 9): 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).For these differences on estimated results, the following reasons can be considered: 1) The standard method assumes the simple observation model in Equation (42).On the other hand, the proposed method introduces the conditional probability distribution in Equation ( 5) as the observation model, which can be considered the whole fluctuation of input and output signals; 2) The standard method is not considered the finite amplitude fluctuation of the observation data; 3) To consider the quantized observation by introducing fuzzy theory is more useful than the standard method by introducing the quantized noise.
The squared sums of the estimation error are shown in Tables 3-5.It can be found numerically that the proposed method is more useful than the extended Kalman filter.

Conclusions
In this study, based on the observed data with fuzziness and the finite level range, a new adaptive method for estimating the input signal for sound environment systems with unknown structure has been proposed.The proposed estimation method has been realized by introducing a system model of conditional probability type and the probability measure of fuzzy events.The proposed method has been applied to the estimation for the input signal of an actual sound environment system, and it has been experimentally verified that better results are obtained as compared with the standard estimation method without considering fuzzy theory.
The proposed approach is quite different from the traditional standard techniques.However, we are still in an early stage of development, and a number of practical problems are yet to be investigated in the future.These include: 1) application to a diverse range of sound signals in actual noise environment; 2) extension to cases with multi-noise sources, and 3) finding an optimal number of expansion terms for the expansion-based probability expressions adopted.

, a b
owing to a definite dynamic range of measurement instrument.The fuzziness of k z is characterized by a member- probability density function is the most standard one.
the expansion coefficients in the equality on Jacobi poly- nomial: 0 , ;

Figures 2 -
Figures 2-4 show the estimation results of the fluctuation wave form of the input signal for Data 1 in a typical case of widely fluctuating signals, by applying the proposed algorithm to the quantized observation data with 1 dB, 2 dB and 3 dB widths.Furthermore, Figures 5-7 show the estimation results for Data 2 in a typical case of narrowly fluctuating signals.In these estimations, the finite number of expansion coefficients

Figure 2 .
Figure 2. Estimation results for Data 1 of the input signal by applying the proposed method based on the quantized output observation data with 1 dB width.

Figure 3 .
Figure 3. Estimation results for Data 1 of the input signal by applying the proposed method based on the quantized output observation data with 2 dB width.

Figure 4 .
Figure 4. Estimation results for Data 1 of the input signal by applying the proposed method based on the output observation data with 3 dB width.
results by the proposed method Estimated results by the extended Kalman filter

Figure 5 .
Figure 5. Estimation results for Data 2 of the input signal by applying the proposed method based on the quantized output observation data with 1 dB width.

Figure 6 .
Figure 6.Estimation results for Data 2 of the input signal by applying the proposed method based on the quantized output observation data with 2 dB width.

Figure 7 .
Figure 7. Estimation results for Data 2 of the input signal by applying the proposed method based on the quantized output observation data with 3 dB width.
results by the proposed method Estimated results by the extended Kalman filter

Table 1 .
Mean and standard deviation of the input signal (in dB).

Table 2 .
Mean and standard deviation of the output signal (in dB).

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 1 dB width (in dB).

Table 4 .
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).

Table 5 .
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 3 dB width (in dB).