Adaptive Signal
Adaptive Signal
2 MATHEMATICAL BACKGROUND
To fully understand the derivation of the Wiener filter, it is essential to review some
foundational concepts related to stochastic processes, random variables, and the mean
square error (MSE) criterion. These mathematical concepts form the basis for defining and
solving the Wiener filter problem.
2.1 Stochastic Processes and Random Variables
A stochastic process is a collection of random variables indexed by time. In Wiener filtering,
both the signal of interest s(n) and the noise v(n) are modeled as stochastic processes. This
means that for each time n, s(n) and v(n) are random variables that follow certain probability
distributions.
The signal s(n) is the underlying process that we aim to estimate. It is a random process, and
its exact values at each time instant are not known deterministically.
The noise v(n), which corrupts the signal, is also modeled as a random process. It is assumed
to be uncorrelated with the signal s(n).
Key characteristics of stochastic processes include:
Mean E[s(n)]: The average or expected value of the process at time n. If the mean of the
process does not change with time, the process is said to be stationary.
E[s(n)]= μs, where μs is constant.
Variance Var[s(n)]: The measure of how much the values of the process fluctuate around the
mean.
Autocorrelation function Rss(m): This function describes how the values of the process at
different times are related to each other.
Let:
s(n) be the true, desired signal.
s’(n) be the estimate of the signal, which is produced by applying the Wiener filter to
the noisy observation x(n).
The MSE is defined as:
MSE=E [ ( s ( n )−^s ( n ) ) ]
2
This expression represents the expected value of the squared difference between the actual
signal s(n) and its estimate ^s ( n )
Minimizing the MSE
The Wiener filter aims to minimize the MSE by finding the optimal filter coefficients that
best approximate the true signal from the observed data. Minimizing the MSE involves
taking the derivative of the MSE with respect to the filter coefficients and setting it to zero,
which leads to the Wiener-Hopf equations.
Let us assume signal ^s ( n ) is generated by filtering the observed signal x(n), which includes
both the desired signal s(n) and the noise v(n). If the filter is represented by a set of
coefficients h(k), the output of the filter is
M −1
^s ( n )= ∑ h ( k ) x ( n−k )
k=0
Where h(k) are the filter coefficients, and M is the length of the filter. The MSE can then be
expressed as a function of these coefficients, and the goal is to find the coefficients that
minimize the MSE.
2.3 Cross-Correlation and Autocorrelation Functions
To derive the Wiener filter, it is essential to use autocorrelation and cross-correlation
functions, which measure the statistical relationships between the signals and the noise.
R xx ( m) =E [ x ( n ) x ( n−m ) ]
Where E[⋅] denotes the expectation operator. This function describes how the values of x(n)
at different time instants are correlated.
Cross-correlation function Rsx(m): The cross-correlation between the desired signal s(n) and
the observed signal x(n) measures how s(n) and x(n)are related. It is defined as:
R xx ( m) =E [ s ( n ) x ( n−m) ]
This function plays a crucial role in determining how the observed noisy signal x(n) can be
used to estimate the desired signal s(n).
The Wiener filter uses these correlation functions to compute the optimal filter coefficients
that minimize the MSE. By exploiting the relationship between the autocorrelation of the
observed signal and the cross-correlation between the observed and desired signals, the
Wiener filter can find the best possible estimate of the desired signal.
3 DERIVATION OF WINER FILTER
Let x(n) represent the observed noisy signal, which is a combination of the desired signal s(n)
and noise v(n):
x(n) = s(n) + v(n)
We want to find a filter h(n) that produces an estimate ^s ( n ) of the signal s(n) from the
observation x(n) :
M −1
^s ( n )= ∑ h ( k ) x ( n−k )
k=0
Where h(k) represents the filter coefficients and M is the length of the filter.
Now we need to define the Mean Square Error (MSE), which is the expected value of the
squared difference between the true signal s(n) and the estimated signal ^s ( n ) :
MSE=E [ ( s ( n )−^s ( n ) ) ]
2
MSE=E [ s ( n )2 ]−2 E [ s ( n ) s^ ( n ) ] + E [ ^s ( n )2 ]
M −1
Using the filter output ^s ( n )= ∑ h ( k ) x ( n−k ), we express the cross-correlation between s(n)
k=0
and ^s ( n ) , as well as the autocorrelation of x(n).
Cross- correlation between s(n) and x(n) :
The cross-correlation function Rsx(m)between the desired signal s(n)s(n)s(n) and the
observed signal x(n) is:
R sx ( m )=E [ s ( n ) x ( n−m ) ]
Autocorrelation of x(n):
The autocorrelation function Rxx (m) of the observed signal x(n) is:
R xx ( m) =E [ x ( n ) x ( n−m ) ]
The Wiener-Hopf equations represent a system of linear equations that can be solved for the
filter coefficients h(k). These equations involve the cross-correlation between the desired
signal and the observed signal, Rsx(k), and the autocorrelation of the observed signal, Rxx(k)
In matrix form, these equations can be written as:
R xx ( h )=R xx
Where:
R xx is the autocorrelation matrix of the observed signal,
h is the vector of filter coefficients,
R sx is the cross-correlation vector between the desired and observed signals.
One of the primary applications of the Wiener filter is noise reduction in signals. When a
signal is corrupted by additive noise, the Wiener filter can estimate the original signal by
minimizing the error between the actual signal and the filtered output.
2. Image Restoration
In image processing, the Wiener filter is commonly used for image deblurring and noise
reduction. When an image is blurred due to motion or other factors, and noise is added, the
Wiener filter can restore the image by estimating the original pixel values.
4. Echo Cancellation
Wiener filters are also employed in echo cancellation systems, particularly in
telecommunications. Echo occurs when a portion of the transmitted signal is reflected back
to the sender, leading to a delayed and distorted version of the original signal.