Chapter - 2 Wiener Filters
Chapter - 2 Wiener Filters
Wiener Filters
Dr. Saad Muhi Falih
3
Mean-Squared Error Criterion
4
Mean-Squared Error Criterion
5
Mean-Squared Error Criterion
• Figure 3.1 shows the block schematic of a linear discrete-time
filter W(z) in the context of estimating a desired signal d(n)
based on an excitation x(n).
• Here, we assume that both x(n) and d(n) are samples of infinite
length random processes.
• The filter output is y(n), and e(n) is the estimation error.
• Clearly, the smaller the estimation error, the better the filter
performance. As the error approaches zero, the output of the
filter approaches the desired signal, d(n).
• the question that arises is the following: what is the most
appropriate choice for the parameters of the filter, which
would result in the smallest possible estimation error?
6
• As we want the estimation error to be as small as
possible, a straightforward approach to the design of
the filter parameters appears to be “to choose an
appropriate function of this estimation error as a cost
function and select that set of the filter parameters
which optimizes this cost function in some sense.”
• In choosing a performance function, the following
points have to be considered:
1. The performance function must be mathematically
tractable.
2. The performance function should preferably have a
single minimum (or maximum) point.
7
• The recursive (infinite-impulse response – IIR)
filters, in general, result in performance functions that may
have many minima (or maxima) points,
• whereas the non-recursive (finite-impulse response – FIR)
filters are guaranteed to have a single global minimum (or
maximum) point if a proper performance function is used.
• Because of this, application of the IIR filters in adaptive
filtering has been very limited.
• In Wiener filters, the performance function is chosen to be:
9
• The filter input, x(n), and its desired output,
d(n), are assumed to be real-valued stationary
processes.
10
11
• This is a quadratic function of the tap-weight vector w with a
single global minimum.
• To obtain the set of tap weights, which minimizes the
performance function ξ, we need to solve the system of
equations that results from setting the partial derivatives of ξ
with respect to every tap weight to zero. That is,
12
13
14
15
16
Example 3.1: Consider the modeling problem shown in
Figure 3.3. The plant is a two-tap filter with an additive
noise, ν(n), added to its output. A two-tap Wiener filter
with tap weights w0 and w1 is used to model the plant
parameters. The same input is applied to both the plant
and Wiener filter. The input, x(n), is a stationary white
process with variance of unity. The additive noise, ν(n),
is zero-mean and uncorrelated with x(n), and its
variance is v = 0.1. We want to compute the optimum
values of w0 and w1, which minimize E[e2(n)].
17
18
19
20
21
3.3 Principle of Orthogonality
23
3.4 Normalized Performance Function
• It is appropriate if we define the ratio
24
• Noting that ζmin cannot be negative,
• we find that its value remains between 0 and
1.
• The value of ζmin is an indication of the ability
of the filter in estimating the desired output.
• A value of ζmin close to zero is an indication
of good performance of the filter, and a value
of ζmin close to 1 indicates poor performance
of the filter.
25