Scientific Research

An Academic Publisher

**A Bayesian Filter for Sound Environment System with Quantized Observation ^{*}** ()

Share and Cite:

^{*}.

*Intelligent Information Management*,

**10**, 87-98. doi: 10.4236/iim.2018.103007.

1. Introduction

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 ${L}_{x}$ (( $100-x$ ) percentile level), ${L}_{eq}$ (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 [1] [2] [3] . However, many standard estimation methods proposed previously in a study of stochastic systems are restricted only to the Gaussian distribution [4] [5] . Several state estimation methods for nonlinear system have been also proposed by assuming the Gaussian distribution of system and observation noises [6] [7] [8] [9] [10] . The real sound environment often shows an intricate fluctuation pattern rather than the standard Gaussian distribution. In our previous studies [11] [12] [13] [14] [15] , 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 [16] [17] [18] . 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 [10] [18] .

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:

${x}_{k+1}=F{x}_{k}+G{u}_{k}$ (1)

where ${x}_{k}$ denotes the specific signal at a discrete time k, ${u}_{k}$ is the random input with known statistics. Here, ${x}_{k}$ and ${u}_{k}$ are statistically independent of each other. Two parameters F and G are estimated by using an auto-correlation technique [1] . Furthermore, the observation model is established by considering the additive property of power variables and the quantized observation in decibel scale, as follows:

${y}_{k}=10{\mathrm{log}}_{10}\left\{\left({x}_{k}+{v}_{k}\right)/{10}^{-12}\right\}$ (2)

${z}_{k}=Q\left({y}_{k}\right)\equiv g\left({x}_{k}+{v}_{k}\right)$ (3)

where ${y}_{k}$ is the noisy observation contaminated by the additive background noise ${v}_{k}$ . Though ${y}_{k}$ 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 $Q\left({y}_{k}\right)$ denotes a nonlinear function expressing the quantization mechanism and ${z}_{k}$ is the quantized observation. Therefore, $g\left(v\right)$ 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 ${x}_{k}$ is proposed on the basis of the quantized observation contaminated by the background noise ${v}_{k}$ .

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 ${z}_{k}$ on the estimated probability density function $P\left({x}_{k}|{Z}_{k}\right)$ by use of various types of linear and/or nonlinear correlation between ${x}_{k}$ and ${z}_{k}$ , the well-known Bayes’ theorem is introduced.

$P\left({x}_{k}|{Z}_{k}\right)=P\left({x}_{k},{z}_{k}|{Z}_{k-1}\right)/P\left({z}_{k}|{Z}_{k-1}\right)$ (4)

where $\left({Z}_{k}\equiv \left\{{z}_{1},{z}_{2},\cdots ,{z}_{k}\right\}\right)$ is a set of observation data up to time k. By expanding the conditional probability density function $P\left({x}_{k}|{Z}_{k}\right)$ 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 [11] [12] .

$P\left({x}_{k}|{Z}_{k}\right)=\frac{{P}_{0}\left({x}_{k}|{Z}_{k-1}\right){\displaystyle \underset{m=0}{\overset{\infty}{\sum}}{\displaystyle \underset{n=0}{\overset{\infty}{\sum}}{A}_{mn}}}{\varphi}_{m}^{\left(1\right)}\left({x}_{k}\right){\varphi}_{n}^{\left(2\right)}\left({z}_{k}\right)}{{\displaystyle \underset{n=0}{\overset{\infty}{\sum}}{A}_{0n}{\varphi}_{n}^{\left(2\right)}\left({z}_{k}\right)}}$ (5)

with

${A}_{mn}\equiv \langle {\varphi}_{m}^{\left(1\right)}\left({x}_{k}\right){\varphi}_{n}^{\left(2\right)}\left({z}_{k}\right)|{Z}_{k-1}\rangle $ (6)

The above two functions ${\varphi}_{m}^{\left(1\right)}\left({x}_{k}\right)$ and ${\varphi}_{n}^{\left(2\right)}\left({z}_{k}\right)$ are orthonormal polynomials of degrees m and n with weighting functions ${P}_{0}\left({x}_{k}|{Z}_{k-1}\right)$ and ${P}_{0}\left({z}_{k}|{Z}_{k-1}\right)$ . Based on Equation (5), the estimate of the polynomial function ${f}_{M}\left({x}_{k}\right)$ of ${x}_{k}$ with Mth order can be derived as follows.

$\begin{array}{c}{\stackrel{^}{f}}_{M}\left({x}_{k}\right)\equiv \langle {f}_{M}\left({x}_{k}\right)|{Z}_{k}\rangle \\ ={\displaystyle \underset{m=0}{\overset{M}{\sum}}{\displaystyle \underset{n=0}{\overset{\infty}{\sum}}{C}_{Mm}{A}_{mn}{\varphi}_{n}^{\left(2\right)}\left({z}_{k}\right)}}/{\displaystyle \underset{n=0}{\overset{\infty}{\sum}}{A}_{0n}{\varphi}_{n}^{\left(2\right)}\left({z}_{k}\right)}\end{array}$ (7)

where ${C}_{Mm}$ is an appropriate constant satisfying the following equality:

${f}_{M}\left({x}_{k}\right)={\displaystyle \underset{m=0}{\overset{M}{\sum}}{C}_{Mm}{\varphi}_{m}^{\left(1\right)}\left({x}_{k}\right)}$ (8)

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:

${P}_{0}\left({x}_{k}|{Z}_{k-1}\right)=N\left({x}_{k};{x}_{k}^{*},{\Gamma}_{{x}_{k}}\right)$

$N\left(x;\mu ,{\sigma}^{2}\right)\equiv \frac{1}{\sqrt{2\text{\pi}{\sigma}^{2}}}\mathrm{exp}\left\{-\frac{{\left(x-\mu \right)}^{2}}{2{\sigma}^{2}}\right\}$ ,

${x}_{k}^{*}\equiv \langle {x}_{k}|{Z}_{k-1}\rangle $ , ${\Gamma}_{{x}_{k}}\equiv \langle {\left({x}_{k}-{x}_{k}^{*}\right)}^{2}|{Z}_{k-1}\rangle $ (9)

Furthermore, as the fundamental probability function on the level-quantized observation, the generalized binomial distribution [19] with level difference interval ${h}_{z}$ can be chosen:

${P}_{0}\left({z}_{k}|{Z}_{k-1}\right)=\frac{\left(\frac{{N}_{k}-{z}_{M}}{{h}_{z}}\right)!}{\left(\frac{{z}_{k}-{z}_{M}}{{h}_{z}}\right)!\left(\frac{{N}_{k}-{z}_{k}}{{h}_{z}}\right)!}{p}_{k}^{\frac{{z}_{k}-{z}_{M}}{{h}_{z}}}{\left(1-{p}_{k}\right)}^{\frac{{N}_{k}-{z}_{k}}{{h}_{z}}}$

${p}_{k}\equiv \frac{{z}_{k}^{*}-{z}_{M}}{{N}_{k}-{z}_{M}}$ , ${N}_{k}\equiv \frac{\left({z}_{k}^{*}-{z}_{M}\right){h}_{z}{z}_{k}^{*}-{z}_{M}{\Omega}_{{z}_{k}}}{\left({z}_{k}^{*}-{z}_{M}\right){h}_{z}-{\Omega}_{{z}_{k}}}$

${z}_{k}^{*}\equiv \langle {z}_{k}|{Z}_{k-1}\rangle $ , ${\Omega}_{{z}_{k}}\equiv \langle {\left({z}_{k}-{z}_{k}^{*}\right)}^{2}|{Z}_{k-1}\rangle $ (10)

where ${z}_{M}$ is the minimum level of observations. The orthonormal polynomials with two weighting probability distributions in Equations (9) and (10) can be determined as

${\varphi}_{m}^{\left(1\right)}\left({x}_{k}\right)=\frac{1}{\sqrt{m!}}{H}_{m}\left(\frac{{x}_{k}-{x}_{k}^{*}}{\sqrt{{\Gamma}_{{x}_{k}}}}\right)$ (11)

$\begin{array}{c}{\varphi}_{n}^{\left(2\right)}\left({z}_{k}\right)={\left\{{\left(\frac{{N}_{k}-{z}_{M}}{{h}_{z}}\right)}^{\left(n\right)}n!\right\}}^{-1/2}{\left(\frac{1-{p}_{k}}{{p}_{k}}\right)}^{n/2}\frac{1}{{h}_{z}^{n}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\cdot {\displaystyle \underset{j=0}{\overset{n}{\sum}}\left(\begin{array}{c}n\\ j\end{array}\right){\left(-1\right)}^{n-j}{\left(\frac{{p}_{k}}{1-{p}_{k}}\right)}^{n-j}}{\left({N}_{k}-{z}_{k}\right)}^{\left(n-j\right)}{\left({z}_{k}-{z}_{M}\right)}^{\left(j\right)}\end{array}$ (12)

where ${H}_{m}(\cdot )$ denotes the Hermite polynomial with mth order, and ${z}^{\left(j\right)}$ is the jth order factorial function defined by [19]

${z}^{\left(n\right)}=z\left(z-{h}_{z}\right)\left(z-2{h}_{z}\right)\cdot \cdot \cdot \left(z-\left(n-1\right){h}_{z}\right),\text{\hspace{0.17em}}{z}^{\left(0\right)}=1$ (13)

Since the function $g(\cdot )$ in Equation (3) is not differentiable in general, the following expansion expression of discrete type is introduced with two arbitrary constants ${d}_{v}$ and ${h}_{v}$ .

$g\left(x+v\right)={\displaystyle \underset{n=0}{\overset{\infty}{\sum}}\frac{1}{n!}{\left(\frac{{d}_{v}}{{h}_{v}}\right)}^{n}{x}^{\left(n\right)}{\Delta}^{n}g\left(v\right)}$ (14)

where ∆ is the forward difference operator defined as:

$\Delta g\left(v\right)\equiv \frac{1}{{d}_{v}}\left\{g\left(v+{h}_{v}\right)-g\left(v\right)\right\}$ (15)

After substituting Equations (2) and (3) into the definition of two parameters of ${z}_{k}$ and ${\Omega}_{{z}_{k}}$ in Equation (10), by applying Equation (14), the following expressions can be derived (in the case of ${d}_{v}={h}_{v}$ ).

$\begin{array}{c}{z}_{k}^{*}=\langle g\left({x}_{k}+{v}_{k}\right)|{Z}_{k-1}\rangle \\ =\langle g\left({v}_{k}\right)+{x}_{k}\Delta g\left({v}_{k}\right)+\frac{1}{2}{x}_{k}\left({x}_{k}-{h}_{v}\right){\Delta}^{2}g\left({v}_{k}\right)+\cdot \cdot \cdot |{Z}_{k-1}\rangle \end{array}$ (16)

$\begin{array}{c}{\Omega}_{{z}_{k}}=\langle {\left\{g\left({x}_{k}+{v}_{k}\right)-{z}_{k}^{*}\right\}}^{2}|{Z}_{k-1}\rangle \\ =\langle {\left\{g\left({v}_{k}\right)+{x}_{k}\Delta g\left({v}_{k}\right)+\frac{1}{2}{x}_{k}\left({x}_{k}-{h}_{v}\right){\Delta}^{2}g\left({v}_{k}\right)+\cdot \cdot \cdot \right\}}^{2}|{Z}_{k-1}\rangle -{z}_{k}^{*}\end{array}$ (17)

Furthermore, the expansion coefficients defined by Equation (6) can be expressed as follows:

$\begin{array}{l}{A}_{mn}={\left\{{\left(\frac{{N}_{k}-{z}_{M}}{{h}_{z}}\right)}^{\left(n\right)}m!n!\right\}}^{-1/2}{\left(\frac{1-{p}_{k}}{{p}_{k}}\right)}^{n/2}\frac{1}{{h}_{z}^{n}}{\displaystyle \underset{j=0}{\overset{n}{\sum}}\left(\begin{array}{c}n\\ j\end{array}\right){\left(-1\right)}^{n-j}{\left(\frac{{p}_{k}}{1-{p}_{k}}\right)}^{n-j}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\cdot \langle {H}_{m}\left(\frac{{x}_{k}-{x}_{k}^{*}}{\sqrt{{\Gamma}_{{x}_{k}}}}\right){\left({N}_{k}-g\left({x}_{k}+{v}_{k}\right)\right)}^{\left(n-j\right)}{\left(g\left({x}_{k}+{v}_{k}\right)-{z}_{M}\right)}^{\left(j\right)}|{Z}_{k-1}\rangle \end{array}$ (18)

The above Equations (16)-(18) can be obtained from the statistics of the background noise ${v}_{k}$ and the predictions of ${x}_{k}$ at a discrete time $k-1$ ; i.e., the expectation values of arbitrary functions of ${x}_{k}$ conditioned by ${Z}_{k-1}$ .

Therefore, the estimation algorithm in Equation (7) can be expressed explicitly as follows:

$\begin{array}{l}{\stackrel{^}{f}}_{M}\left({x}_{k}\right)={\displaystyle \underset{m=0}{\overset{M}{\sum}}{\displaystyle \underset{n=0}{\overset{\infty}{\sum}}{C}_{Mm}{A}_{mn}}{\left\{{\left(\frac{{N}_{k}-{z}_{M}}{{h}_{z}}\right)}^{\left(n\right)}n!\right\}}^{-1/2}{\left(\frac{1-{p}_{k}}{{p}_{k}}\right)}^{n/2}\frac{1}{{h}_{z}^{n}}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\cdot {\displaystyle \underset{j=0}{\overset{n}{\sum}}\left(\begin{array}{c}n\\ j\end{array}\right){\left(-1\right)}^{n-j}{\left(\frac{{p}_{k}}{1-{p}_{k}}\right)}^{n-j}}{\left({N}_{k}-{z}_{k}\right)}^{\left(n-j\right)}{\left({z}_{k}-{z}_{M}\right)}^{\left(j\right)}/I\end{array}$ (19)

with

$\begin{array}{l}I={\displaystyle \underset{n=0}{\overset{\infty}{\sum}}{A}_{0n}{\left\{{\left(\frac{{N}_{k}-{z}_{M}}{{h}_{z}}\right)}^{\left(n\right)}n!\right\}}^{-1/2}{\left(\frac{1-{p}_{k}}{{p}_{k}}\right)}^{n/2}\frac{1}{{h}_{z}^{n}}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\cdot {{\displaystyle \underset{j=0}{\overset{n}{\sum}}\left(\begin{array}{c}n\\ j\end{array}\right){\left(-1\right)}^{n-j}\left(\frac{{p}_{k}}{1-{p}_{k}}\right)}}^{n-j}{\left({N}_{k}-{z}_{k}\right)}^{\left(n-j\right)}{\left({z}_{k}-{z}_{M}\right)}^{\left(j\right)}\end{array}$ (20)

Especially, the estimates for mean and variance can be obtained as follows:

$\begin{array}{l}{\stackrel{^}{x}}_{k}\equiv \langle {x}_{k}|{Z}_{k}\rangle \\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}={\displaystyle \underset{m=0}{\overset{1}{\sum}}{\displaystyle \underset{n=0}{\overset{\infty}{\sum}}{C}_{1m}{A}_{mn}}{\left\{{\left(\frac{{N}_{k}-{z}_{M}}{{h}_{z}}\right)}^{\left(n\right)}n!\right\}}^{-1/2}{\left(\frac{1-{p}_{k}}{{p}_{k}}\right)}^{n/2}\frac{1}{{h}_{z}^{n}}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\cdot {\displaystyle \underset{j=0}{\overset{n}{\sum}}\left(\begin{array}{c}n\\ j\end{array}\right){\left(-1\right)}^{n-j}{\left(\frac{{p}_{k}}{1-{p}_{k}}\right)}^{n-j}}{\left({N}_{k}-{z}_{k}\right)}^{\left(n-j\right)}{\left({z}_{k}-{z}_{M}\right)}^{\left(j\right)}/I\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left({C}_{10}={x}_{k}^{*},{C}_{11}=\sqrt{{\Gamma}_{{x}_{k}}}\right)\end{array}$ (21)

$\begin{array}{l}{P}_{k}\equiv \langle {\left({x}_{k}-{\stackrel{^}{x}}_{k}\right)}^{2}|{Z}_{k}\rangle \\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}={\displaystyle \underset{m=0}{\overset{2}{\sum}}{\displaystyle \underset{n=0}{\overset{\infty}{\sum}}{C}_{2m}{A}_{mn}}{\left\{{\left(\frac{{N}_{k}-{z}_{M}}{{h}_{z}}\right)}^{\left(n\right)}n!\right\}}^{-1/2}{\left(\frac{1-{p}_{k}}{{p}_{k}}\right)}^{n/2}\frac{1}{{h}_{z}^{n}}}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\cdot {\displaystyle \underset{j=0}{\overset{n}{\sum}}\left(\begin{array}{c}n\\ j\end{array}\right){\left(-1\right)}^{n-j}{\left(\frac{{p}_{k}}{1-{p}_{k}}\right)}^{n-j}}{\left({N}_{k}-{z}_{k}\right)}^{\left(n-j\right)}{\left({z}_{k}-{z}_{M}\right)}^{\left(j\right)}/I\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left({C}_{20}={\Gamma}_{{x}_{k}}+{\left({x}_{k}^{*}-{\stackrel{^}{x}}_{k}\right)}^{2},{C}_{21}=2\sqrt{{\Gamma}_{{x}_{k}}}\left({x}_{k}^{*}-{\stackrel{^}{x}}_{k}\right),{C}_{22}=\sqrt{2}{\Gamma}_{{x}_{k}}\right)\end{array}$ (22)

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 [16] . 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 ${P}_{0}\left({x}_{k}|{Z}_{k-1}\right)$ and ${P}_{0}\left({z}_{k}|{Z}_{k-1}\right)$ , because this probability density function is the most standard one.

${P}_{0}\left({x}_{k}|{Z}_{k-1}\right)=N\left({x}_{k};{x}_{k}^{*},{\Gamma}_{{x}_{k}}\right)$ , ${P}_{0}\left({z}_{k}|{Z}_{k-1}\right)=N\left({z}_{k};{z}_{k}^{*},{\Omega}_{{z}_{k}}\right)$ (23)

Then, the orthonormal functions with two weighting probability density functions in (23) can be given in the Hermite polynomial:

${\varphi}_{m}^{\left(1\right)}\left({x}_{k}\right)=\frac{1}{\sqrt{m!}}{H}_{m}\left(\frac{{x}_{k}-{x}_{k}^{*}}{\sqrt{{\Gamma}_{{x}_{k}}}}\right)$ , ${\varphi}_{n}^{\left(2\right)}\left({z}_{k}\right)=\frac{1}{\sqrt{n!}}{H}_{n}\left(\frac{{z}_{k}-{z}_{k}^{*}}{\sqrt{{\Omega}_{{z}_{k}}}}\right)$ (24)

Accordingly, the estimation algorithm of the specific signal in Equation (7) can be given by

${\stackrel{^}{f}}_{M}\left({x}_{k}\right)={\displaystyle \underset{m=0}{\overset{M}{\sum}}{\displaystyle \underset{n=0}{\overset{\infty}{\sum}}{C}_{Mm}{A}_{mn}\frac{1}{\sqrt{n!}}{H}_{n}\left(\frac{{z}_{k}-{z}_{k}^{*}}{\sqrt{{\Omega}_{{z}_{k}}}}\right)/J}}$ (25)

with

$J={\displaystyle \underset{n=0}{\overset{\infty}{\sum}}{A}_{0n}\frac{1}{\sqrt{n!}}{H}_{n}\left(\frac{{z}_{k}-{z}_{k}^{*}}{\sqrt{{\Omega}_{{z}_{k}}}}\right)}$ (26)

Furthermore, the estimates for mean and variance can be obtained as follows:

${\stackrel{^}{x}}_{k}={\displaystyle \underset{m=0}{\overset{1}{\sum}}{\displaystyle \underset{n=0}{\overset{\infty}{\sum}}{C}_{1m}{A}_{mn}\frac{1}{\sqrt{n!}}{H}_{n}\left(\frac{{z}_{k}-{z}_{k}^{*}}{\sqrt{{\Omega}_{{z}_{k}}}}\right)/J}}$ (27)

${P}_{k}={\displaystyle \underset{m=0}{\overset{2}{\sum}}{\displaystyle \underset{n=0}{\overset{\infty}{\sum}}{C}_{2m}{A}_{mn}\frac{1}{\sqrt{n!}}{H}_{n}\left(\frac{{z}_{k}-{z}_{k}^{*}}{\sqrt{{\Omega}_{{z}_{k}}}}\right)/J}}$ (28)

Thus, two parameters ${z}_{k}^{*}$ , ${\Omega}_{{z}_{k}}$ and the expansion coefficients ${A}_{mn}$ are expressed as follows:

$\begin{array}{c}{z}_{k}^{*}=\langle g\left({x}_{k}+{v}_{k}\right)|{Z}_{k-1}\rangle \\ ={\displaystyle \iint g\left({x}_{k}+{v}_{k}\right)P\left({x}_{k}|{Z}_{k-1}\right)P\left({v}_{k}\right)\text{d}{x}_{k}\text{d}{v}_{k}}\end{array}$ (29)

$\begin{array}{c}{\Omega}_{{z}_{k}}=\langle {\left(g\left({x}_{k}+{v}_{k}\right)-{z}_{k}^{*}\right)}^{2}|{Z}_{k-1}\rangle \\ ={\displaystyle \iint {\left(g\left({x}_{k}+{v}_{k}\right)-{z}_{k}^{*}\right)}^{2}P\left({x}_{k}|{Z}_{k-1}\right)P\left({v}_{k}\right)\text{d}{x}_{k}\text{d}{v}_{k}}\end{array}$ (30)

$\begin{array}{c}{A}_{mn}=\langle \frac{1}{\sqrt{m!}}{H}_{m}\left(\frac{{x}_{k}-{x}_{k}^{*}}{\sqrt{{\Gamma}_{{x}_{k}}}}\right)\frac{1}{\sqrt{n!}}{H}_{n}\left(\frac{g\left({x}_{k}+{v}_{k}\right)-{z}_{k}^{*}}{\sqrt{{\Omega}_{{z}_{k}}}}\right)|{Z}_{k-1}\rangle \\ =\frac{1}{\sqrt{m!}}\frac{1}{\sqrt{n!}}{\displaystyle \iint {H}_{m}\left(\frac{{x}_{k}-{x}_{k}^{*}}{\sqrt{{\Gamma}_{{x}_{k}}}}\right){H}_{n}\left(\frac{g\left({x}_{k}+{v}_{k}\right)-{z}_{k}^{*}}{\sqrt{{\Omega}_{{z}_{k}}}}\right)P\left({x}_{k}|{Z}_{k-1}\right)P\left({v}_{k}\right)\text{d}{x}_{k}\text{d}{v}_{k}}\end{array}$ (31)

The integrals in Equations (29)-(31) are evaluated by use of particles with mean ${x}_{k}^{*}$ , variance ${\text{\Gamma}}_{{x}_{k}}$ , higher order statistics ${A}_{m0}$ 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

${x}_{k+1}^{*}\equiv \langle {x}_{k+1}|{Z}_{k}\rangle =F{\stackrel{^}{x}}_{k}+G\langle {u}_{k}\rangle $ (32)

${\Gamma}_{{x}_{k+1}}\equiv \langle {\left({x}_{k+1}-{x}_{k+1}^{*}\right)}^{2}|{Z}_{k}\rangle ={F}^{2}{P}_{k}+{G}^{2}\langle {\left({u}_{k}-\langle {u}_{k}\rangle \right)}^{2}\rangle $ (33)

By replacing k with $k+1$ , 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/m^{2}).

Table 2. Mean and standard deviation of the background noise (in W/m^{2}).

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 ${A}_{mn}\left(m,n\le 2\right)$ 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 [11] [12] , the extended Kalman filter is also applied to the observation data as a trail by introducing the following observation model.

${z}_{k}=10{\mathrm{log}}_{10}\left({x}_{k}+{v}_{k}\right)/{10}^{-12}+{\epsilon}_{k}$ (34)

where ${\epsilon}_{k}$ denotes the quantized noise and a uniform distribution within $\left[-q/2,q/2\right]$ (q: the quantized width) is assumed as the probability distribution of ${\epsilon}_{k}$ . 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.

5. Conclusions

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.

Acknowledgements

The author is grateful for valuable suggestions at Inter-Noise 2016 [20] . 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.

NOTES

*New type method is proposed based on the Baysian filter by use of the quantized observation in sound environment system.

Conflicts of Interest

The authors declare no conflicts of interest.

[1] | Eykhoff, P. (1974) System Identification: Parameter and State Estimation. Wiley, New York. |

[2] |
Young, P. (1984) Recursive Estimation and Time-Series Analysis. Springer-Verlag, Berlin. https://doi.org/10.1007/978-3-642-82336-7 |

[3] | Gremal, M.S. and Andrews, A.P. (1993) Kalman Filtering-Theory and Practice. Prentice-Hall, New Jersey. |

[4] |
Kalman, R.E. (1960) A New Approach to Linear Filtering and Prediction Problems. ASME, Series D, Journal of Basic Engineering, 82, 35-45.
https://doi.org/10.1115/1.3662552 |

[5] |
Kalman, R.E. and Buch, R.S. (1961) New Results in Linear Filtering and Prediction Theory. Transactions on ASME, Series D, Journal of Basic Engineering, 83, 95-108.
https://doi.org/10.1115/1.3658902 |

[6] |
Kushner, H.J. (1967) Approximations to Optimal Nonlinear Filter. IEEE Transactions on Automatic Control, 12, 546-556.
https://doi.org/10.1109/TAC.1967.1098671 |

[7] |
Bell, B.M. and Cathey, F.W. (1993) The Iterated Kalman Filter Update as a Gauss-Newton Methods. IEEE Transactions on Automatic Control, 38, 294-297.
https://doi.org/10.1109/9.250476 |

[8] |
Nishiyama, K. (1997) A Nonlinear Filter for Estimating a Sinusoidal Signal and Its Parameter: On the Case of a Single Sinusoid. IEEE Transactions on Signal Processing, 45, 970-981. https://doi.org/10.1109/78.564185 |

[9] | Vincent, T.L. and Khargonekar, P.P. (1999) A Class of Nonlinear Filtering Problems Arising from Drift Sensor Gains. IEEE Transactions on Automatic Control, 44, 509-520. |

[10] |
Julier, S.J. (2002) The Scaled Unscented Transformation. Proceedings of the American Control Conference, 6, 4555-4559. https://doi.org/10.1109/ACC.2002.1025369 |

[11] | Ohta, M. and Yamada, H. (1984) New Methodological Trials of Dynamical State Estimation for the Noise and Vibration Environmental System—Establishment of General Theory and Its Application to Urban Noise Problems. Acustica, 55, 199-212. |

[12] | Ikuta, A. and Ohta, M. (1992) A State Estimation Method of Impulsive Signal Using Digital Filter under the Existence of External Noise and Its Application to Room Acoustics. IEICE Transactions on Fundamentals of Electronics, Communications and Computer Sciences, E75, 988-995. |

[13] |
Ikuta, A., Ohta, M. and Masuike, H. (2006) A Countermeasure for an External Noise in the Measurement of Sound Environment and Its Application to the Evaluation for Traffic Noise at Main Line. IEEJ Transactions on EIS, 126, 63-71.
https://doi.org/10.1541/ieejeiss.126.63 |

[14] |
Orimoto, H., Ikuta, A. and Xiao, Y. (2016) Method for Evaluating the Statistical Relationship between Sound Pressure Level and Noise Annoyance Based on a Nonlinear Time Series Regression Model and an Experiment. International Journal of Acoustics and Vibration, 21, 145-151. https://doi.org/10.20855/ijav.2016.21.2403 |

[15] |
Ikuta, A., Orimoto, H. and Gallagher, G. (2017) State Estimation for Fuzzy Sound Environment System with Finite Amplitude Fluctuation. Journal of Software Engineering and Applications, 10, 625-638. https://doi.org/10.4236/jsea.2017.107034 |

[16] |
Habbibi, B., Sayed, A.H. and Kailath, T. (1996) Linear Estimation in Krein Spaces—Parts I Theory and II Applications. IEEE Transactions on Automatic Control, 41, 18-33, 34-49. https://doi.org/10.1109/9.481605 |

[17] |
Ikuta, A., Masuike, H. and Ohta, M. (2007) State Estimation for Sound Environment System with Unknown Structure by Introducing Fuzzy Theory. The Institute of Electrical Engineers of Japan, 127, 770-777.
https://doi.org/10.1541/ieejeiss.127.770 |

[18] | Kitagawa, G. (1996) Monte Carlo Filter and Smoother for Non-Gaussian Nonlinear State Space Models. Journal of Computational and Graphical Statistics, 5, 1-25. |

[19] |
Ikuta, A., Ohta, M. and Ogawa, H. (1997) Various Regression Characteristics with Higher Order among Light, Sound and Electromagnetic Waves Leaked from VDT—Measurement and Signal Processing in the Actual Working Environment. International Measurement Confederation, 21, 25-33.
https://doi.org/10.1016/S0263-2241(97)00041-9 |

[20] | Ikuta, A. (2016) A Bayesian Filter for Sound Environment System with Quantized Observation. INTER-NOISE, Hamburg, 21-24 August2016, 3184-3192. |

Copyright © 2020 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.