The goal is to find the distribution (mean and variance) of noisy speech parameter , given the distributions for noise parameters and , and clean speech parameter . First, mean, i.e. expectation value, of is given as:

(5) |

Provided , and have Gaussian distributions, does not have any closed-form expressions.

Now, the task here is to approximate the function by some polynomial that can closely approximate it within a given range, with its low order as far as possible. The polynomial approximation can be simplified by reducing function to univariate one, by taking , where .

We choose here 2nd-order Lagrange polynomial to approximate the function
, as given by:

The points , and can be specified manually (one point at and other two chosen to minimize error in the required range), or instead, Chebyshev-Lagrange polynomial can be used that specifies the points itself.

Finally, Eq.(6) reduces to
form, where ,
and are constants. Therefore:

(7) |

The accurate value of mean is more important than that of variance. Therefore, the covariance matrix can be retained as it is for the clean speech. However, expression for adapting diagonal variance can be derived from above approximation, in terms of higher-order moments (up to 4th moment) of , and diagonal variances can be adapted as well.