AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent...

39
1 Mintav´ etelit¨orv´ eny AD ´ atalak´ ıt´ ok + s´ avlimit´ alt jel τ id˝ok¨ oz¨ onk´ ent mintav´ etel Mintav´ etel: szorz´ as az id˝o tartom´ anyban 1orfrekvenci´ ank´ ent ism´ etl˝ odik - konvol´ uci´ o a frekvenciat´ erben.

Transcript of AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent...

Page 1: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

1

Mintaveteli torveny

AD atalakıtok + savlimitalt jel

τ idokozonkent mintavetel

Mintavetel: szorzas az ido tartomanyban

1/τ korfrekvenciankent ismetlodik - konvolucio a frekvenciaterben.

Page 2: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

2

Nem fednek at:

Mintaveteli torveny

fmintavetel ≥ 2fmax

fmintavetel: Nyquist frekvencia

Page 3: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

3

Page 4: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

4

Ez semmi mas, mint a lebeges jelensege!

Page 5: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

5

• a jel legmagasabb frekvenciaju komponensebol is legalabb ket mintat kell

vennunk.

• az ilyen modon mintavett jelekbol az eredeti jel tokeletesen re-

konstrualhato, ha a mintavett impulzusokat egy - jol megvalasztott - felso

hatarfrekvenciaju idealis alulatereszton bocsatjuk keresztul (34 abra).

Page 6: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

6

• savlimitalas→ az AD konverter bemenete elott valodi alulatereszto szuro

Page 7: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

7

Teljes DSP folyamat:

A DA atalakıto kimenete altalaban”szogletes”, a beepıtett mintatarto

aramkor miatt:

sin(x)/x korrekcio a frekvenciaterben!

Page 8: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

8

Page 9: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

9

CCD:

Page 10: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

10

Page 11: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

11

Page 12: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

12

TV videojel:

Page 13: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

13

Digitalis szurok

Teljes DSP folyamat:

Szuro:

kesleltetoelem + osszegezes:

xi bemeneti adatok, ezeket eltaroljuk/kesleltejuk, es sulyozva osszeadjuk:

Page 14: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

14

Lehet maskepp is elkepzelni:

Page 15: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

15

Page 16: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

16

Minden elem ugyanakkor kesleltetesnek felel meg:

gyurus tarolo

Page 17: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

17

Vizsgalati modszer:

sulyfuggveny: Dirac-delta be→ hi sulyok jelennek meg a kimeneten:

y[n] = h0x[n] + h1x[n− 1] + h2x[n− 2] + h3x[n− 3] + · · ·

y[n] =

NX

i=−M

hix[n− i]

y[n] =

NX

i=−M

aix[n− i]

y[n] =

∞X

i=−∞aix[n− i]

FIR (Finite Impulse Response, veges impulzusvalaszu): mindig veges

hosszusagu kimenet.

Page 18: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

18

A kimeneten a hn sulyok jelennek meg!

Pl.

Masik pelda:

Furcsa szorzas - polinomok szorzasa! A polinomok egyutthatoi az xn es hn

adatok.

Ez nem mas, mint a konvolucio!

Page 19: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

19

Legyen τ kesleltetes egy elemben. Ekkor az x(t) bemeno fuggvenyen egy

elem kesleltetese eit/τ szorzofaktort hoz be. Azaz

y(ti) =

NX

i=−M

hix(ti − τi)(ejt/τ

)i

=

NX

i=−M

x(ti − τi)hi(ejt/τ

)i

=

NX

i=−M

x(ti − τi)z−i

Itt z = hieit/τ es j =

√−1.

Ez az u.n. z-transzformacio - lathatoan kozeli a kapcsolata a

Fourier-transzformacioval!

Page 20: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

20

Elozo pelda:

z-transzformacioval:

A z-transzformaltakra igaz, hogy :

z(Y ) = z(H)z(X)

A z-transzformacio kozeli kapcsolatban all a Fourier-transzformacioval!

Page 21: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

21

Rekurzıv szurok

Most mar a korabbi y[i] ertekek is szerepet jatszanak:

y[n] = a0x[n] + a1x[n− 1] + a2x[n− 2] + a3x[n− 3] + · · ·+b0y[n] + b1y[n− 1] + b2y[n− 2] + b3y[n− 3] + · · ·

de logikusan b0 = 0 !

y[n] = a0x[n] +a1x[n− 1] + a2x[n− 2] + a3x[n− 3] + · · ·+b1y[n− 1] + b2y[n− 2] + b3y[n− 3] + · · ·

Page 22: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

22

IIR (Infinite Impulse Response, vegtelen impulzusvalaszu): altalaban nem

tudjuk, hogy meddig tart - vegtelen hosszu is lehet (oszcillacio)!

Page 23: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

23

Egyszeru pelda:

Sulyfuggveny: Dirac-delta a bemeneten

Pl. ha B1 = 0.9 , akkor a kimenet: 1.0, 0.9, 0.81 · · ·

H(z) =

∞Xm=0

z−m

Bm

=

∞Xm=0

(B

z)m

=1

1− Bz−1

Page 24: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

24

Ez csak Bz−1 < 1 -re konvergens.

Page 25: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

25

Maskepp:

yn = xn + Byn−1

yn = xn + Bynz−1

yn =1

1− Bz−1x

n

Ez pontosan olyan, mint egy kvaziintegralo aramkor!

Page 26: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

26

Hogyan lehet kvazidifferencialo aramkor jellegu szurohoz jutni?

Page 27: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

27

Page 28: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

28

z-transzformalt:

yn = A0xn + A1xnz−1 − Bynz

−1

yn =A0 + A1z

−1

1− Bz−1x

n

Tortfuggveny - hol lesz 0 a nevezo?

Ha A0 = −A1 = 1, akkor

Page 29: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

29

Ez felulatereszto jellegu:

Page 30: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

30

Page 31: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

31

Page 32: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

32

Page 33: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

33

Altalanos alak:

H(z) =A0 + A1z

−1 + A2z−2 + · · ·+ Akz

−k

1 + B1z−1 + B2z−2 + · · ·+ Bmz−m

Polinom gyokok - stabilitasvizsgalat - van-e oszcillacio?

Parcialis gyokokre bontas segıthet

Page 34: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

34

Legyen a kesleltetes valodi (muvonal, uvegszal, stb.) !

Bemenetre Dirac-delta← τ kesleltetes← e−iωτ impulzusatvitel.

Amplitudomenet:

|H(ω)| = | 1

1− Be−iωτ|

1√1 + B2 − 2B cos ωτ

Page 35: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

35

Page 36: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

36

Fazismenet: FIR: linearis

IIR: nem-linearis

Page 37: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

37

Page 38: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

38

DSP elonyok

Alkalmazasi teruletek:

• hıradastechnikai jelek (kommunikacios halozatok, TV, radio)

• kepfeldolgozas (kamera)

• hangfrekvencias technika (GSM, beszedfelismeres, szintetizalas, stb.)

• meresadat feldolgozas (fizikus, vegyesz, geologus, biologus, stb.)

Celprocesszorok (on-line) es szamıtogepes programok (off-line)

Mire hasznaljuk?

• Zavarok, zajok csokkentese: savlimitalas, Wiener szuro

• Dekonvolucios / korrekcios halozatok: meres felbontokepessegenek

novelese

Miert jobbak az analog eszkozoknel?

• nagyobb pontossag

• akarmilyen szuro (lehet, hogy ilyen nem epıtheto valos RLC elemekbol)

Page 39: AD ´atalak´ıt´ok + s´avlimit´alt jel id˝ok¨oz¨onk´ent …itl7.elte.hu/~zsolt/Oktatas/JX/j5.pdfAD ´atalak´ıt´ok + s´avlimit´alt jel τ id˝ok¨oz¨onk´ent mintav´etel

39

• a szurokarakterisztika egyszeruen modosıthato

• nem erzekenyek a kornyezetre

• kisebb suly, olcsosag