adaptive overview talk

22
Page 1 An Overview of Optimal and Adaptive Digital Signal Processing Dr. Michael T. Johnson Marquette University Department of Electrical and Computer Engineering

Upload: hdavuluri

Post on 08-Apr-2018

221 views

Category:

Documents


0 download

TRANSCRIPT

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 1/22

Page 1

An Overview of Optimal and Adaptive

Digital Signal Processing

Dr. Michael T. JohnsonMarquette University

Department of Electrical and Computer Engineering

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 2/22

Page 2

Marquette University

Optimal filtering

Basic purpose of filtering: get rid of undesired noise, keep the desired signal.³Optimal´ filter: A filter whose parametersminimize a selected error criterion (MMSE)

OptimalFilter 7

+y

In put

d

e

u -

Error

Refere nce

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 3/22

Page 3

Marquette University

Adaptive filtering

When signal or noise components arenon-stationary, the optimal filter isconstantly changing.³Adaptive´ filter: A filter which has time-varying parameters (like an FIR filter withweights that slowly change over time

Ad aptiveFilter 7

+yIn putd

e

u -

Error

Refere nce

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 4/22

Page 4

Marquette University

Classes of adaptive filter problems

System Identification

Inverse Modeling

u=filter input, d=desired signal, e=error, y=output

Unkn ownSystem

Ad aptiveFilter

7-

+

y

In put d

e

u

Dela y

Ad aptive

Filter 7 -

+

y

In put d

e

uUnkn own

System

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 5/22

Page 5

Marquette University

Prediction

Noise suppression

Ad aptiveFilter 7

+yIn put

d

e

uDela y-

Pred

iction error Pred icted signal

Ad aptiveFilter 7+yIn put

d

e

u -

Output

Refere nce

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 6/22

Page 6

Marquette University

Room acoustic identification

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 7/22

Page 7

Marquette University

Echo cancellation

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 8/22

Page 8

Marquette University

Noise cancellation

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 9/22

Page 9

Marquette University

Channel equalization

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 10/22

Page 10

Marquette University

Active Noise Control

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 11/22

Page 11

Marquette University

Narrowband supression

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 12/22

Page 12

Marquette University

Notation

u(n)=[u[n] u[n-1] « u[n-M+1] ] T , the input window vector r(m) = R XX[m]=E{u(n) u(n-m)}, the m th autocorrelation coefficientr =[r *(1) r *(2) « r *(M)]T , the autocorrelation vector of the inputR = E{u(n) u(n) H} = the autocorrelation matrix of the input

(0) (1) ( 1)

( 1) (0)

( 1) (0)

r r r M

r r R

r M r

« »¬ ¼¬ ¼! ¬ ¼¬ ¼ ½

L

M O

OptimalFilter 7

+yIn put

d

e

u -

Error

Refere nce

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 13/22

Page 13

Marquette University

Optimality: Wiener filters

_ a*

1 12 * *

0 0

arg mi n [ ] [ ] [ ] [ ] [ ]M M

o i i

i i

J E e n E d n w u n i d n w u n i

! !

¨ ¸® ¾¨ ¸¨ ¸± ±© ¹! ! ! ¯ ¿© ¹© ¹© ¹ª º ª º± ±° Àª º§ §

ww

10 !w R p

where w 0 is the vector of M filter weightsR is the correlation matrixp is the cross-correlation vector between u and d

OptimalFilter 7

+yIn put

d

e

u -

Error

Refere nce

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 14/22

Page 14

Marquette University

Basic Adaptation Theory

Can analyze error of optimal filter

And use concept of gradient descent

To determine a theoretical update formula

_ a*

1 12 * *

0 0

2

[ ] [ ] [ ] [ ] [ ] M M

i i

i i

d

J e n d n u n i d n u n i

W

! !

® ¾¨ ¸¨ ¸! ! ¯ ¿© ¹© ¹ª º ª º À

!

§ §H H Hw p p w w Rw

_ a

12

*

( 1) ( ) ( ( ))

( ) [ ] [ ]

( ) ( )

n n J n

n E n e n

n n

Q

Q

Q

!

!

!

w w w

w u

w p Rw

1( 1) ( ) ( )

2n n n Q ! w w g

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 15/22

Page 15

Marquette University

L east Mean Square ( L MS)

A pp roxim ate p =E {u (n )d* (n )} and R=

E {u (n )u* (n )} us ing o ne wi nd ow of the s ignal:

To g iv e the LMS update r ule:

*

Ö

( ) ( ) ( )Ö ( ) ( ) ( )

H n n n

n n d n

!!

R u u

p u

*

*

Ö Ö( 1) ( ) ( ) ( ) ( ) ( )Ö ( ) ( ) ( )

H n n n d n n n

n n e n

Q

Q

!

! w w u u w

w u

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 16/22

Page 16

Marquette University

Update process

Step-by-step:1. Filter output2. Error signal3. Weight-update

Ö( ) ( ) ( ) H y n n n! w u

( ) ( ) ( )e n d n y n! *Ö Ö( 1) ( ) ( ) ( )n n n e n Q!w w u

z - 1I

I

u (n) u H(n)

7 w [n]w [n+ 1]

7Q -

d *

(n)

u(n)

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 17/22

Page 17

Marquette University

Analysis of adaptive filtersL

MS converges to optimal Wiener solution plusrandom component that acts as Brownian motionCan control learning rate and window size M

Stability

Maximum learning rate a function of dataConvergence

Convergence speed is directly proportional to

MisadjustmentRandom misadjustment is directly proportional to

and to M

Must balance convergence and misadjustment

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 18/22

Page 18

Marquette University

Adaptive equalizer example

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 19/22

Page 19

Marquette University

As function of data characteristics

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 20/22

Page 20

Marquette University

As function of learning rate

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 21/22

Page 21

Marquette University

L MS variants

Many variants to L MSNormalized L MS (N L MS)Regularized L MS ( -N L MS)L eaky L MS

Sign-error L MSFX-L MS (compensation for active noise control)Affine Projection Algorithm (APA)Block L MS in frequency domain

DCT- L MSSubband- L MSRecursive L east Squares (R L S)

8/6/2019 Adaptive Overview Talk

http://slidepdf.com/reader/full/adaptive-overview-talk 22/22

Page 22

Marquette University

Example results