21.01.2015 Views

MIK 200 Anvendt signalbehandling, 2012. Lab. 8, Notch-filter.

MIK 200 Anvendt signalbehandling, 2012. Lab. 8, Notch-filter.

MIK 200 Anvendt signalbehandling, 2012. Lab. 8, Notch-filter.

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

2 Litt teori om IIR-<strong>filter</strong> med mer.<br />

Dette er et noe stikkordspreget og kladdepreget notat over litt teorien fra<br />

<strong>signalbehandling</strong>en. Det er det en bør kunne for å få til denne øvinga.<br />

2.1 IIR-<strong>filter</strong><br />

Et IIR-<strong>filter</strong> er et <strong>filter</strong> med (teoretisk) uendelig impulsrespons. I praksis er det<br />

implementert nesten like enklet som et FIR-<strong>filter</strong>. Utgangen av et FIR-<strong>filter</strong> er<br />

vekta sum av inngangssignal ved dette og noen foregående tidspunkt.<br />

y(k) =<br />

N−1<br />

∑<br />

n=0<br />

For eksempel for et boxcar-<strong>filter</strong> med lengde 4<br />

b(n)x(k − n). (3)<br />

y(k) = x(k) + x(k − 1) + x(k − 2) + x(k − 3) (4)<br />

For et IIR-<strong>filter</strong> har en med ledd av foregående utganger også.<br />

y(k) =<br />

N−1<br />

∑<br />

n=0<br />

b(n)x(k − n) −<br />

M−1<br />

∑<br />

m=1<br />

a(m)y(k − m) (5)<br />

Merk at andre sum begynner med m = 1. Et eksempel på et IIR-<strong>filter</strong> er et<br />

AR(1)-<strong>filter</strong> (her med ρ = 0.95)<br />

y(k) = x(k) + 0.95y(k − 1) (6)<br />

Vi ser ganske enkelt at dette <strong>filter</strong>et vil få uendlig lang impulsrespons, selv<br />

om utgangsverdiene går mot null så vil de aldri bli helt like null. Dere stusser<br />

gjerne på minustegnet framfor det siste summetegnet i ligning 5. Ved å flytte<br />

summen over på venstre siden, og la a(0) = 1 kan en skrive<br />

M−1<br />

∑<br />

m=0<br />

a(m)y(k − m) = y(k) +<br />

M−1<br />

∑<br />

m=1<br />

a(m)y(k − m) =<br />

Nå kan en ta z-transform på begge sider og får<br />

som gir transferfunksjonen H(z).<br />

N−1<br />

∑<br />

n=0<br />

b(n)x(k − n) (7)<br />

A(z)Y (z) = B(z)X(z) (8)<br />

Y (z) = B(z) X(z) = H(z)X(z) (9)<br />

A(z)<br />

H(z) = Y (z)<br />

X(z) = B(z)<br />

A(z)<br />

(10)<br />

10

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!