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.
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