Scribd
Upload
237 views

LMS Algorithm

The least mean squares (LMS) filter is an adaptive filter used to minimize the error between the desired and actual signals by adjusting the filter coefficients. It was developed in 1960 and uses a stochastic gradient descent method to update the filter coefficients based only on the current error. The LMS filter is related to the Wiener filter but does not rely on correlations and converges to the Wiener solution. It provides a simple yet effective way to adaptively filter signals.

Uploaded by

menilanjan89nL
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
237 views

LMS Algorithm

The least mean squares (LMS) filter is an adaptive filter used to minimize the error between the desired and actual signals by adjusting the filter coefficients. It was developed in 1960 and uses a stochastic gradient descent method to update the filter coefficients based only on the current error. The LMS filter is related to the Wiener filter but does not rely on correlations and converges to the Wiener solution. It provides a simple yet effective way to adaptively filter signals.

Uploaded by

menilanjan89nL
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 8

Least mean squares filter

From Wikipedia, the free encyclopedia


Least mean squares (LMS) algorithms are a class of adaptive filter used to mimic a desired filter by
finding the filter coefficients that relate to producing the least mean squares of the error signal (difference
between the desired and the actual signal). It is a stochastic gradient descent method in that the filter is only
adapted based on the error at the current time. It was invented in 1960 by Stanford University professor
Bernard Widrow and his first Ph.D. student, Ted Hoff.

Contents

1 Problem formulation
o 1.1 Relationship to the least mean squares filter
o 1.2 Definition of symbols
2 Idea
3 Derivation
4 Simplifications
5 LMS algorithm summary
6 Convergence and stability in the mean
7 Normalised least mean squares filter (NLMS)
o 7.1 Optimal learning rate
o 7.2 Proof
8 See also
9 References
10 External links

Problem formulation

Relationship to the least mean squares filter


The realization of the causal Wiener filter looks a lot like the solution to the least squares estimate, except in
the signal processing domain. The least squares solution, for input matrix and output vector is

The FIR least mean squares filter is related to the Wiener filter, but minimizing the error criterion of the
former does not rely on cross-correlations or auto-correlations. Its solution converges to the Wiener filter

solution. Most linear adaptive filtering problems can be formulated using the block diagram above. That is,
an unknown system

is to be identified and the adaptive filter attempts to adapt the filter

it as close as possible to
and

, while using only observable signals

and

; but

to make
,

are not directly observable. Its solution is closely related to the Wiener filter.

Definition of symbols
is the number of the current input sample
is the filter order
(Hermitian transpose or conjugate transpose)

estimated filter; interpret as the estimation of the filter coefficients after n samples

Idea
The basic idea behind LMS filter is to approach the optimum filter weights
, by updating the filter
weights in a manner to converge to the optimum filter weight. The algorithm starts by assuming a small
weights (zero in most cases), and at each step, by finding the gradient of the mean square error, the weights
are updated. That is, if the MSE-gradient is positive, it implies, the error would keep increasing positively, if
the same weight is used for further iterations, which means we need to reduce the weights. In the same way,
if the gradient is negative, we need to increase the weights. So, the basic weight update equation is :
,
where represents the mean-square error. The negative sign indicates that, we need to change the weights in
a direction opposite to that of the gradient slope.
The mean-square error, as a function of filter weights is a quadratic function which means it has only one
extrema, that minimises the mean-square error, which is the optimal weight. The LMS thus, approaches
towards this optimal weights by ascending/descending down the mean-square-error vs filter weight curve.

Derivation
The idea behind LMS filters is to use steepest descent to find filter weights
function. We start by defining the cost function as

where

is the error at the current sample n and

which minimize a cost

denotes the expected value.

This cost function (


) is the mean square error, and it is minimized by the LMS. This is where the
LMS gets its name. Applying steepest descent means to take the partial derivatives with respect to the
individual entries of the filter coefficient (weight) vector

where

is the gradient operator

Now,

is a vector which points towards the steepest ascent of the cost function. To find the

minimum of the cost function we need to take a step in the opposite direction of
mathematical terms

. To express that in

where is the step size(adaptation constant). That means we have found a sequential update algorithm
which minimizes the cost function. Unfortunately, this algorithm is not realizable until we know
.
Generally, the expectation above is not computed. Instead, to run the LMS in an online (updating after each
new sample is received) environment, we use an instantaneous estimate of that expectation. See below.

Simplifications
For most systems the expectation function
the following unbiased estimator

where

must be approximated. This can be done with

indicates the number of samples we use for that estimate. The simplest case is

For that simple case the update algorithm follows as

Indeed this constitutes the update algorithm for the LMS filter.

LMS algorithm summary


The LMS algorithm for a
Parameters:
Initialisation:
Computation: For

th order algorithm can be summarized as

filter order
step size

Convergence and stability in the mean


As the LMS algorithm does not use the exact values of the expectations, the weights would never reach the
optimal weights in the absolute sense, but a convergence is possible in mean. That is, even though the
weights may change by small amounts, it changes about the optimal weights. However, if the variance with
which the weights change, is large, convergence in mean would be misleading. This problem may occur, if
the value of step-size is not chosen properly.
If is chosen to be large, the amount with which the weights change depends heavily on the gradient
estimate, and so the weights may change by a large value so that gradient which was negative at the first
instant may now become positive. And at the second instant, the weight may change in the opposite
direction by a large amount because of the negative gradient and would thus keep oscillating with a large
variance about the optimal weights. On the other hand if is chosen to be too small, time to converge to the
optimal weights will be too large.

Thus, an upper bound on


where

is needed which is given as

is the greatest eigenvalue of the autocorrelation matrix

condition is not fulfilled, the algorithm becomes unstable and

. If this
diverges.

Maximum convergence speed is achieved when

where
is the smallest eigenvalue of R. Given that is less than or equal to this optimum, the
convergence speed is determined by
, with a larger value yielding faster convergence. This means that
faster convergence can be achieved when
is close to
, that is, the maximum achievable
convergence speed depends on the eigenvalue spread of
.
A white noise signal has autocorrelation matrix
where
is the variance of the signal. In this
case all eigenvalues are equal, and the eigenvalue spread is the minimum over all possible matrices. The
common interpretation of this result is therefore that the LMS converges quickly for white input signals, and
slowly for colored input signals, such as processes with low-pass or high-pass characteristics.
It is important to note that the above upperbound on
coefficients of
practical bound is

only enforces stability in the mean, but the

can still grow infinitely large, i.e. divergence of the coefficients is still possible. A more

where
denotes the trace of
. This bound guarantees that the coefficients of
do not diverge (in
practice, the value of should not be chosen close to this upper bound, since it is somewhat optimistic due
to approximations and assumptions made in the derivation of the bound).

Normalised least mean squares filter (NLMS)


The main drawback of the "pure" LMS algorithm is that it is sensitive to the scaling of its input
. This
makes it very hard (if not impossible) to choose a learning rate that guarantees stability of the algorithm
(Haykin 2002). The Normalised least mean squares filter (NLMS) is a variant of the LMS algorithm that
solves this problem by normalising with the power of the input. The NLMS algorithm can be summarised
as:
Parameters:

filter order
step size

Initialization:
Computation: For

Optimal learning rate


It can be shown that if there is no interference (
algorithm is

and is independent of the input


with interference (

), then the optimal learning rate for the NLMS

and the real (unknown) impulse response

. In the general case

), the optimal learning rate is

The results above assume that the signals


the case in practice.

and

are uncorrelated to each other, which is generally

Proof

Let the filter misalignment be defined as


misalignment for the next sample as:

, we can derive the expected

Let

and

Assuming independence, we have:

The optimal learning rate is found at

, which leads to:

See also

Recursive least squares


For statistical techniques relevant to LMS filter see Least squares.
Similarities between Wiener and LMS
Multidelay block frequency domain adaptive filter
Zero forcing equalizer
Kernel adaptive filter

References

Monson H. Hayes: Statistical Digital Signal Processing and Modeling, Wiley, 1996, ISBN 0-47159431-8
Simon Haykin: Adaptive Filter Theory, Prentice Hall, 2002, ISBN 0-13-048434-2
Simon S. Haykin, Bernard Widrow (Editor): Least-Mean-Square Adaptive Filters, Wiley, 2003,
ISBN 0-471-21570-8
Bernard Widrow, Samuel D. Stearns: Adaptive Signal Processing, Prentice Hall, 1985, ISBN 0-13004029-0
Weifeng Liu, Jose Principe and Simon Haykin: Kernel Adaptive Filtering: A Comprehensive
Introduction, John Wiley, 2010, ISBN 0-470-44753-2
Paulo S.R. Diniz: Adaptive Filtering: Algorithms and Practical Implementation, Kluwer Academic
Publishers, 1997, ISBN 0-7923-9912-9

External links

LMS Algorithm in Adaptive Antenna Arrays www.antenna-theory.com


LMS Noise cancellation demo www.advsolned.com

Categories:

Digital signal processing


Filter theory
Stochastic algorithms

Navigation menu

Create account
Log in

Article
Talk

Read
Edit
View history

Main page
Contents
Featured content
Current events
Random article
Donate to Wikipedia
Wikimedia Shop

Interaction

Help
About Wikipedia
Community portal
Recent changes
Contact page

Tools

What links here


Related changes
Upload file
Special pages
Permanent link
Page information
Wikidata item
Cite this page

Print/export

Create a book
Download as PDF

Printable version

Languages

Catal
Deutsch
Espaol

Italiano

Edit links

You might also like