„Méréstechnika szerkesztőlap” változatai közötti eltérés

A Fizipedia wikiből
(A Lock-in erősítő)
(Spektrumanalízis)
 
29. sor: 29. sor:
  
 
Könnyen belátható, hogy a diszkrét Fourier-transzformált fenti képlet szerinti kiértékelése $N$ mérési pont esetén $\sim N^2$ művelet ($\Delta \omega \approx 2\pi/N\Delta t$ frekvenciafelbontás és $\omega_\mathrm{max}=2\pi/2\Delta t$ maximális felbontható frekvencia esetén csak $\omega_\mathrm{max}/\Delta \omega\approx N$ diszkrét pontban érdemes kiértékelni a diszkrét Fourier-transzformáltat, és a definíció szerint egy adott frekvencián $~N$ művelet a szumma kiszámítása). Egy ügyes trükkel azonban jelentősen csökkenthető a számítási műveletek mennyisége. Megmutatható, hogy ha a mérési pontok száma kettő hatványa ($N=2^p$), és a frekvenciatérben $\omega_k=2\pi k/N\Delta t,\ \ \  k=0,1,...,N/2$ diszkrét körfrekvenciáknál értékeljük ki a Fourier-transzformáltat, akkor az ún. Fast Fourier Transform (FFT) algoritmus segítségével a számítási műveletek száma $N^2$-ről $N\log_2 N$-re csökken, ami nagy N esetén lényeges különbség.
 
Könnyen belátható, hogy a diszkrét Fourier-transzformált fenti képlet szerinti kiértékelése $N$ mérési pont esetén $\sim N^2$ művelet ($\Delta \omega \approx 2\pi/N\Delta t$ frekvenciafelbontás és $\omega_\mathrm{max}=2\pi/2\Delta t$ maximális felbontható frekvencia esetén csak $\omega_\mathrm{max}/\Delta \omega\approx N$ diszkrét pontban érdemes kiértékelni a diszkrét Fourier-transzformáltat, és a definíció szerint egy adott frekvencián $~N$ művelet a szumma kiszámítása). Egy ügyes trükkel azonban jelentősen csökkenthető a számítási műveletek mennyisége. Megmutatható, hogy ha a mérési pontok száma kettő hatványa ($N=2^p$), és a frekvenciatérben $\omega_k=2\pi k/N\Delta t,\ \ \  k=0,1,...,N/2$ diszkrét körfrekvenciáknál értékeljük ki a Fourier-transzformáltat, akkor az ún. Fast Fourier Transform (FFT) algoritmus segítségével a számítási műveletek száma $N^2$-ről $N\log_2 N$-re csökken, ami nagy N esetén lényeges különbség.
 
A mérőműszerek jelentős része, így a laborgyakorlaton használt digitális oszcilloszkóp is az FFT algoritmus numerikus kiértékelése alapján határozza meg a mért jel spektrumát. A legtöbb esetben a műszer nem adja meg külön a spektrum valós és képzetes részét, hanem csak a Fourier-transzformált abszolút érték négyzetét látjuk. Ezen kívül a mérőműszerek általában a frekvencia, és nem a körfrekvencia függvényében adják meg a spektrumot, erre érdemes odafigyelni a mérés kiértékelésénél.
 
  
 
A Fourier-transzformáció nem csak a jel frekvenciájának a vizsgálatára alkalmas, hanem a jel amplitúdójának a mérésére is. A fenti képletek szerint egy $A\exp(i\omega_0 t)$ harmonikus jel téglalap ablakkal vett Fourier transzformáltja az $\omega_0$ körfrekvenciánál pontosan $A$ amplitudóval rendelkezik. Az $A$ amplitúdót azonban csak akkor kapjuk vissza, ha a jel körfrekvenciája megegyezik az FFT valamelyik frekvenciapontjával ($\omega_0=2\pi n/T$), ami akkor teljesül, ha a jelnek pontosan egész számú periódusa fér el a $T$ mintavételezési ablakban. (Érdemes megjegyezni, hogy ebben az esetben az FFT összes többi frekvenciapontja a téglalapablakkal számolt Fourier-transzformált zérushelyeivel esik egybe, azaz a spektrális szivárgás is megszűnik.) Azonban ha a mérési ablak szélessége nem egyezik meg a jel periódusidejének egész számú többszörösével, akkor az FFT-algoritmus által számolt frekvenciapontok biztosan eltérnek a mért jelünk frekvenciájától, így az amplitúdóra is pontatlan értéket kapunk. Ha a jel amplitúdója érdekel minket, akkor olyan ablakfüggvényt érdemes választni, melynek a Fourier-transzormáltja a centrális csúcs körüli $2\pi/T$ szélességű tartományban minél kevesebbet változik. '''Téglalapablak esetén az amplitúdópontosság mintegy 36%, Hanning-ablak esetén pedig 15%. Speciális, ún. Flattop-ablakot választva ennél sokkal jobb, mintegy 0.01%-os amplitúdópontosság is elérhető, viszont cserébe lényegesen leromlik a frekvenciafelbontásunk.'''   
 
A Fourier-transzformáció nem csak a jel frekvenciájának a vizsgálatára alkalmas, hanem a jel amplitúdójának a mérésére is. A fenti képletek szerint egy $A\exp(i\omega_0 t)$ harmonikus jel téglalap ablakkal vett Fourier transzformáltja az $\omega_0$ körfrekvenciánál pontosan $A$ amplitudóval rendelkezik. Az $A$ amplitúdót azonban csak akkor kapjuk vissza, ha a jel körfrekvenciája megegyezik az FFT valamelyik frekvenciapontjával ($\omega_0=2\pi n/T$), ami akkor teljesül, ha a jelnek pontosan egész számú periódusa fér el a $T$ mintavételezési ablakban. (Érdemes megjegyezni, hogy ebben az esetben az FFT összes többi frekvenciapontja a téglalapablakkal számolt Fourier-transzformált zérushelyeivel esik egybe, azaz a spektrális szivárgás is megszűnik.) Azonban ha a mérési ablak szélessége nem egyezik meg a jel periódusidejének egész számú többszörösével, akkor az FFT-algoritmus által számolt frekvenciapontok biztosan eltérnek a mért jelünk frekvenciájától, így az amplitúdóra is pontatlan értéket kapunk. Ha a jel amplitúdója érdekel minket, akkor olyan ablakfüggvényt érdemes választni, melynek a Fourier-transzormáltja a centrális csúcs körüli $2\pi/T$ szélességű tartományban minél kevesebbet változik. '''Téglalapablak esetén az amplitúdópontosság mintegy 36%, Hanning-ablak esetén pedig 15%. Speciális, ún. Flattop-ablakot választva ennél sokkal jobb, mintegy 0.01%-os amplitúdópontosság is elérhető, viszont cserébe lényegesen leromlik a frekvenciafelbontásunk.'''   
 
</wlatex>
 
</wlatex>
 +
 
==A zaj==
 
==A zaj==
 
===A zaj definíciója===
 
===A zaj definíciója===

A lap jelenlegi, 2019. október 23., 15:11-kori változata

Tartalomjegyzék

Spektrumanalízis

Egy \setbox0\hbox{$F(t)$}% \message{//depth:\the\dp0//}% \box0% időfüggvény különböző frekvenciájú komponenseinek felbontását matematikailag a Fourier-transzformált segítségével adhatjuk meg:

\[f(\omega)=\int_{-\infty}^{\infty} F(t)\mathrm{e}^{-i\omega t}\mathrm{d}t.\]

Egy valós mérésnél a Fourier-transzformált függvényt csak közelítőleg tudjuk megadni, hiszen egyrészt véges ideig tart a mérésünk, másrészt a mérési adatok csak diszkrét időfelbontással álnak rendelkezésre. Először nézzük meg a véges idejű mérés hatását a Fourier-transzformáltra.

A véges idejű mérés megfelel annak, mintha az eredeti függvényt megszoroznánk a mérési intervallumnak megfelelő \setbox0\hbox{$W(t)$}% \message{//depth:\the\dp0//}% \box0% ablakfüggvénnyel, és ezen szorzatfüggvény Fourier-transzformáltját számolnánk ki:

\[f_W(\omega)=\int_{-\infty}^{\infty} W(t)\cdot F(t)\mathrm{e}^{-i\omega t}\mathrm{d}t,\]

ahol a \setbox0\hbox{$W(t)$}% \message{//depth:\the\dp0//}% \box0% függvény \setbox0\hbox{$|t|<T/2$}% \message{//depth:\the\dp0//}% \box0% esetén \setbox0\hbox{$1/T$}% \message{//depth:\the\dp0//}% \box0%, ezen \setbox0\hbox{$T$}% \message{//depth:\the\dp0//}% \box0% hosszúságú időintervallumon kívül pedig zérus. Megmutatható, hogy egy szorzatfüggvény Fourier-transzformáltja a két komponens Fourier-transzformáltjának a konvolúciója, azaz:

\[ f_W(\omega)=\int_{-\infty}^{\infty} f(\omega ')w(\omega - \omega ')\frac{\mathrm{d}\omega}{2\pi},\ \ \ \ \mathrm{ahol}\ \ \  w(\omega)=\int_{-\infty}^{\infty} W(t)\mathrm{e}^{-i\omega t}\mathrm{d}t. \]

Nézzünk egy egyszerű példát, legyen \setbox0\hbox{$F(t)=A\cdot \exp(i\omega_0 t)$}% \message{//depth:\the\dp0//}% \box0% egy harmonikus függvény, melynek a Fourier-transzformáltja egy Dirac-delta függvény: \setbox0\hbox{$f(\omega)=A\cdot 2\pi\delta(\omega-\omega_0)$}% \message{//depth:\the\dp0//}% \box0%. Véges idejű mérés esetén azonban a Fourier integrál értéke a fentiek alapján \setbox0\hbox{$f_W(\omega)=A\cdot w(\omega-\omega_0)$}% \message{//depth:\the\dp0//}% \box0%, azaz a harmonikus függvény Fourier-transzormáltjában egy valós mérés esetén a végtelenül keskeny Dirac-delta csúcs helyett az ablakfüggvény Fourier-transzormáltját látjuk az \setbox0\hbox{$\omega_0$}% \message{//depth:\the\dp0//}% \box0% középfrekvenciához eltolva. A fent definiált téglalap ablak esetén (azaz amikor \setbox0\hbox{$W(t)$}% \message{//depth:\the\dp0//}% \box0% egy \setbox0\hbox{$T$}% \message{//depth:\the\dp0//}% \box0% szélességű intervallumban konstans, azon kívül zérus, lásd 5a. ábra, kék folytonos vonal) az ablakfüggvény Fourier-transzformáltja \setbox0\hbox{$w(\omega)=(2/\omega T)\cdot \sin(\omega T/2)$}% \message{//depth:\the\dp0//}% \box0%, azaz \setbox0\hbox{$f_W(\omega)=A\cdot\left(2/(\omega-\omega_0 ) T \right)\cdot \sin\left((\omega-\omega_0 ) T/2\right)$}% \message{//depth:\the\dp0//}% \box0% (5b. ábra, kék folytonos vonal). A véges időintervallumra számolt Fourier-integrál is mutat egy határozott csúcsot az \setbox0\hbox{$\omega_0$}% \message{//depth:\the\dp0//}% \box0% középfrekvencia körül, azonban ez a csúcs véges szélességű, ráadásul a csúcstól távolabb is oszcillációkat látunk a Fourier-transzformáltban, amit spektrális szivárgásnak nevezünk. Az \setbox0\hbox{$f_W(\omega)$}% \message{//depth:\the\dp0//}% \box0% függvény \setbox0\hbox{$\omega_0$}% \message{//depth:\the\dp0//}% \box0% melletti első zérushelyeinek a távolsága \setbox0\hbox{$4\pi/T$}% \message{//depth:\the\dp0//}% \box0%, így az \setbox0\hbox{$\omega_0$}% \message{//depth:\the\dp0//}% \box0% körüli csúcs szélessége \setbox0\hbox{$\sim 2\pi/T$}% \message{//depth:\the\dp0//}% \box0%. Tehát az első fontos konklúzió, hogy véges időtartamú mérés esetén a jelünket a Fourier-térben csak véges, nagyságrendileg \setbox0\hbox{$\Delta \omega \approx 2\pi/T$}% \message{//depth:\the\dp0//}% \box0% frekvenciafelbontással látjuk!

Ablak.jpg
5. ábra. a) \setbox0\hbox{$W(t)$}% \message{//depth:\the\dp0//}% \box0% ablakfüggvény téglalap ablak (kék folytonos vonal) és Hanning ablak (piros szaggatott vonal) esetén. b) \setbox0\hbox{$F(t)=A\cdot \exp(i\omega_0 t)$}% \message{//depth:\the\dp0//}% \box0% harmonikus jel Fourier-transzformáltjának abszolút érték négyzete téglalap ablak (kék folytonos vonal) és Hanning ablak (piros szaggatott vonal) esetén. A téglalap ablakot Hanning ablakra cserélve az \setbox0\hbox{$\omega_0$}% \message{//depth:\the\dp0//}% \box0% körüli csúcs kiszélesedik, azaz romlik a frekvenciafelbontás, azonban az \setbox0\hbox{$\omega_0$}% \message{//depth:\the\dp0//}% \box0%-tól távolabbi mellékcsúcsok amplitúdója lecsökken, azaz csökken a spektrális szivárgás.

Érdemes megjegyezni, hogy a fent említett téglalap ablak helyett választhatunk más ablakfüggvényt is, például \setbox0\hbox{$W(t)=\cos^2(t\pi/T)$}% \message{//depth:\the\dp0//}% \box0% ún. Hanning-ablak esetén a mért jelben elnyomjuk a \setbox0\hbox{$|t|<T/2$}% \message{//depth:\the\dp0//}% \box0% mintavételezési időablak széleihez közeli részeket (5a. ábra, piros szaggatott vonal). Ebben az esetben az \setbox0\hbox{$\omega_0$}% \message{//depth:\the\dp0//}% \box0% körfrekvenciás jel Fourier-transzformáltjában \setbox0\hbox{$\omega_0$}% \message{//depth:\the\dp0//}% \box0% körül egy még szélesebb csúcsot látunk (azaz a frekvenciafelbontás romlik), viszont az \setbox0\hbox{$\omega_0$}% \message{//depth:\the\dp0//}% \box0%-tól távolabbi oszcillációk amplitúdója (az ún. spektrális szivárgás) lecsökken (5b. ábra, piros szaggatott vonal).

Következő lépésként nézzük meg, hogy mi a hatása annak, hogy a jelünket nem folytonosan látjuk, hanem csak diszkrét mintavételezési időpontokban. Emiatt a jel Fourier-transzformáltját a folytonos integrál helyett kénytelenek vagyunk egy diszkrét összeggel, az ún. diszkrét Fourier-transzformálttal (DFT) közelíteni:

\[ f_W(\omega)=\sum_{n=0}^{N-1} W(n\cdot \Delta t)F(n\cdot \Delta t) \mathrm{e}^{-i\omega n \Delta t}\Delta t, \]

ahol \setbox0\hbox{$\Delta t$}% \message{//depth:\the\dp0//}% \box0% a szomszédos mérési pontok közötti idő, \setbox0\hbox{$N=T/\Delta t$}% \message{//depth:\the\dp0//}% \box0% pedig a mintavételezett pontok száma. Az ún. Nyquist-Shannon mintavételezési törvény szerint \setbox0\hbox{$\Delta t$}% \message{//depth:\the\dp0//}% \box0% sűrűségű mintavételezés esetén a jelet \setbox0\hbox{$\omega_\mathrm{max}=2\pi/2\Delta t$}% \message{//depth:\the\dp0//}% \box0% maximális körfrekvenciáig tudjuk rekonstruálni.

Könnyen belátható, hogy a diszkrét Fourier-transzformált fenti képlet szerinti kiértékelése \setbox0\hbox{$N$}% \message{//depth:\the\dp0//}% \box0% mérési pont esetén \setbox0\hbox{$\sim N^2$}% \message{//depth:\the\dp0//}% \box0% művelet (\setbox0\hbox{$\Delta \omega \approx 2\pi/N\Delta t$}% \message{//depth:\the\dp0//}% \box0% frekvenciafelbontás és \setbox0\hbox{$\omega_\mathrm{max}=2\pi/2\Delta t$}% \message{//depth:\the\dp0//}% \box0% maximális felbontható frekvencia esetén csak \setbox0\hbox{$\omega_\mathrm{max}/\Delta \omega\approx N$}% \message{//depth:\the\dp0//}% \box0% diszkrét pontban érdemes kiértékelni a diszkrét Fourier-transzformáltat, és a definíció szerint egy adott frekvencián \setbox0\hbox{$~N$}% \message{//depth:\the\dp0//}% \box0% művelet a szumma kiszámítása). Egy ügyes trükkel azonban jelentősen csökkenthető a számítási műveletek mennyisége. Megmutatható, hogy ha a mérési pontok száma kettő hatványa (\setbox0\hbox{$N=2^p$}% \message{//depth:\the\dp0//}% \box0%), és a frekvenciatérben \setbox0\hbox{$\omega_k=2\pi k/N\Delta t,\ \ \  k=0,1,...,N/2$}% \message{//depth:\the\dp0//}% \box0% diszkrét körfrekvenciáknál értékeljük ki a Fourier-transzformáltat, akkor az ún. Fast Fourier Transform (FFT) algoritmus segítségével a számítási műveletek száma \setbox0\hbox{$N^2$}% \message{//depth:\the\dp0//}% \box0%-ről \setbox0\hbox{$N\log_2 N$}% \message{//depth:\the\dp0//}% \box0%-re csökken, ami nagy N esetén lényeges különbség.

A Fourier-transzformáció nem csak a jel frekvenciájának a vizsgálatára alkalmas, hanem a jel amplitúdójának a mérésére is. A fenti képletek szerint egy \setbox0\hbox{$A\exp(i\omega_0 t)$}% \message{//depth:\the\dp0//}% \box0% harmonikus jel téglalap ablakkal vett Fourier transzformáltja az \setbox0\hbox{$\omega_0$}% \message{//depth:\the\dp0//}% \box0% körfrekvenciánál pontosan \setbox0\hbox{$A$}% \message{//depth:\the\dp0//}% \box0% amplitudóval rendelkezik. Az \setbox0\hbox{$A$}% \message{//depth:\the\dp0//}% \box0% amplitúdót azonban csak akkor kapjuk vissza, ha a jel körfrekvenciája megegyezik az FFT valamelyik frekvenciapontjával (\setbox0\hbox{$\omega_0=2\pi n/T$}% \message{//depth:\the\dp0//}% \box0%), ami akkor teljesül, ha a jelnek pontosan egész számú periódusa fér el a \setbox0\hbox{$T$}% \message{//depth:\the\dp0//}% \box0% mintavételezési ablakban. (Érdemes megjegyezni, hogy ebben az esetben az FFT összes többi frekvenciapontja a téglalapablakkal számolt Fourier-transzformált zérushelyeivel esik egybe, azaz a spektrális szivárgás is megszűnik.) Azonban ha a mérési ablak szélessége nem egyezik meg a jel periódusidejének egész számú többszörösével, akkor az FFT-algoritmus által számolt frekvenciapontok biztosan eltérnek a mért jelünk frekvenciájától, így az amplitúdóra is pontatlan értéket kapunk. Ha a jel amplitúdója érdekel minket, akkor olyan ablakfüggvényt érdemes választni, melynek a Fourier-transzormáltja a centrális csúcs körüli \setbox0\hbox{$2\pi/T$}% \message{//depth:\the\dp0//}% \box0% szélességű tartományban minél kevesebbet változik. Téglalapablak esetén az amplitúdópontosság mintegy 36%, Hanning-ablak esetén pedig 15%. Speciális, ún. Flattop-ablakot választva ennél sokkal jobb, mintegy 0.01%-os amplitúdópontosság is elérhető, viszont cserébe lényegesen leromlik a frekvenciafelbontásunk.

A zaj

A zaj definíciója


Egy \setbox0\hbox{$V(t)$}% \message{//depth:\the\dp0//}% \box0% időben fluktuáló zaj jellegű (azaz széles frekvenciaspektrumú) jelet egy \setbox0\hbox{$f_0$}% \message{//depth:\the\dp0//}% \box0% középfrekvencia körüli \setbox0\hbox{$\Delta f$}% \message{//depth:\the\dp0//}% \box0% szélességű sávszűrűn keresztül vizsgálva azt találjuk, hogy a szűrt jel szórásnégyzete arányos a frekvenciatartomány \setbox0\hbox{$\Delta f$}% \message{//depth:\the\dp0//}% \box0% szélességével, az arányossági tényező a feszültségzaj spektrális sűrűsége (spectral density of noise), azaz a teljesítmény spektrális sűrűsége (Power Spectral Density - PSD):

\[\left\langle (\Delta V(t|f_0,\Delta f))^2 \right\rangle=s_V(f_0)\Delta f.\]

Ez a zaj kísérleti definíciója. A feszültség szórásnégyzete a zajsűrűség teljes frekvenciatartományra vett integráljával egyenlő.

\[\langle(\Delta V(t))^2 \rangle=\int_{0}^{\infty} df s_V(f).\]

Ez a képlet közvetlenül is használható a zajspektrum kísérleti meghatározásához, azonban emellett jó kiindulópontként szolgál más fizikai mennyiségekkel való kapcsolatának vizsgálatára is. Érdemes megvizsgálni a zajsűrűség és a feszültség korrelációs függvény kapcsolatát. Ez utóbbi a következőképpen definiálható:

\[C(\Delta t)=\langle\Delta V(t)\cdotp\Delta V(t+\Delta t) \rangle=\lim\limits_{T\rightarrow\infty}\frac{1}{T}\int_{0}^{T}dt\Delta V(t)\cdotp\Delta V(t+\Delta t).\]

Könnyen belátható, hogy a korrelációs függvény \setbox0\hbox{$\Delta t=0$}% \message{//depth:\the\dp0//}% \box0%-ban megegyezik a feszültség szórásnégyzetével.

\[C(0)=\langle(\Delta V(t))^2 \rangle.\]

Most fejezzük ki a korrelációs függvényt a Fourier-transzformáltja segítségével:

\[C(t)=\frac{1}{2\pi}\int_{-\infty}^{\infty}d\omega e^{i\omega t}c(\omega).\]

Ez alapján a feszültség szórásnégyzete:

\[\langle(\Delta V(t))^2 \rangle=C(0)=\frac{1}{2\pi}\int_{-\infty}^{\infty}d\omega c(\omega)=\frac{1}{\pi}\int_{0}^{\infty}d\omega c(\omega),\]

mivel \setbox0\hbox{$c(\omega)=c(-\omega)$}% \message{//depth:\the\dp0//}% \box0%. Másrészt a szórásnégyzet és a zaj közötti összefüggésből következik:

\[\langle(\Delta V(t))^2 \rangle=\frac{1}{2\pi} \int_{0}^{\infty} d\omega s_V(\omega).\]

Az utóbbi két egyenletből látható, hogy a zaj teljesítménysűrűsége a feszültség korrelációs függvény Fourier-transzformáltjának kétszerese.

Sikerült kapcsolatot teremteni a zajsűrűség és a feszültség korrelációs függvény között. Most vizsgáljuk meg a feszültség korrelációs függvény és a feszültség mint mérhető fizikai mennyiségek kapcsolatát. A feszültség átlagtól való eltérésének (\setbox0\hbox{$\Delta V(t)=V(t)-\langle V(t)\rangle $}% \message{//depth:\the\dp0//}% \box0%) a Fourier-transzformáltja szerint:

\[\Delta v(\omega)= \int_{-\infty}^{\infty} dt \Delta V(t)e^{-i\omega t}.\]

Vizsgáljuk ezen Fourier-transzformált abszolút érték négyzetének várható értékét:

\[\langle|(\Delta v(\omega))|^2 \rangle= \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} dt' dt \langle\Delta V(t)\Delta V(t') \rangle e^{-i\omega t} e^{i\omega t'}.\]

Az egyenlet jobb oldalán látható korrelációs függvény kifejezhető a következőképpen:

\[ \langle\Delta V(t)\Delta V(t')\rangle=\frac{1}{2\pi}\int_{-\infty}^{\infty}d\omega'c(\omega')e^{i\omega'(t-t')}.\]

Továbbá az exponenciális tagokat átcsoportosítva a következő tag integrálja egy Dirac-deltát ad:

\[\int_{-\infty}^{\infty}dt'\frac{e^{i(\omega+\omega')t'}}{2\pi}=\delta(\omega+\omega').\]

Így a feszültség Fourier-transzformált abszolút érték négyzetének várható értéke könnyen belátható, hogy arányos a korrelációs függvény Fourier-transzformáltjával, illetve az időablak szélességével:

\[\langle|(\Delta v(\omega))|^2 \rangle= c(\omega) \int_{-\infty}^{\infty} dt.\]


Innen könnyen kifejezhető a zaj teljesítménysűrűsége, felhasználva, hogy az a feszültség korrelációs függvény Fourier-transzformáltjának kétszerese.

\[s_V(\omega)=\lim\limits_{T->\infty} \frac{2}{T}\left\langle |\int_{-T/2}^{T/2} dt \Delta V(t)e^{-i\omega t}|^2\right\rangle.\]

A mérésünk során diszkrét pontokon történik a mintavételezés, így ugyanez a számolás DFT segítségével a következő módon fejezhető ki:

\[s_V(\omega)\approx\frac{2}{N\Delta t}|\sum_{n=0}^{N-1} \Delta V(n\Delta t)e^{-i\omega n\Delta t}\Delta t|^2=\frac{2\Delta t}{N}|\sum_{n=0}^{N-1} \Delta V(n\Delta t)e^{-i\omega n\Delta t}|^2.\]

A fentiekben sikerült megállapítanunk, hogy a vizsgált \setbox0\hbox{$\Delta V(t)$}% \message{//depth:\the\dp0//}% \box0% jel Fourier-transzformáltjából hogyan számolható a zaj spektrális sűrűsége. Azonban ennél a számolásnál nem vettünk figyelembe ablakfüggvényt, vagy fogalmazhatunk úgy is, hogy téglalap ablakkal számoltunk.

Nézzük meg, hogy egy tetszőleges ablakfüggvény esetén hogyan származtatható a zajsűrűség. Sajnos tetszőleges spektrumú zajra és tetszőleges ablakfüggvényre általános összefüggés nem adható, viszont ha \setbox0\hbox{$\Delta V(t)$}% \message{//depth:\the\dp0//}% \box0% fehér zaj, akkor tetszőleges ablakfüggvényre egyszerűen számolható a konverziós faktor.

Az ablakfüggvénnyel szorzott jel Fourier-transzformáltja:

\[\Delta v_W(\omega)=\int_{-\infty}^{\infty}dt\Delta V(t)W(t)e^{-i\omega t}.\]

Vizsgáljuk meg ezen Fourier-transzformált abszolút érték négyzetének várható értékét, ekkor a következőt kapjuk:

\[\langle|\Delta v(\omega)|^2\rangle=\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}dtdt'W(t)W(t')\langle\Delta V(t)\Delta V(t')\rangle e^{-i\omega t}e^{-i\omega t'}.\]

A fehér zaj jellegéből következik, hogy:

\[\langle\Delta V(t)\Delta V(t')\rangle=C(t-t')=\frac{s_0}{2}\delta(t-t').\]

Így a Fourier-transzformált abszolút érték négyzetének várható értéke a következőképpen egyszerűsödik:

\[\langle|\Delta v_W(\omega)|^2\rangle=\frac{s_0}{2}\int_{-\infty}^{\infty}W^2(t)dt.\]

Azaz a zajsűrűség számolása:

\[s_0=\frac{2\langle|v_W(\omega)|^2\rangle}{\int_{-\infty}^{\infty}W^2(t)dt}.\]

A feszültségzaj spektrális sűrűsége egy normálási faktor erejéig egyenlő a mért feszültségpontok diszkrét Fourier-transzformáltjának abszolútérték-négyzetének kétszeresével. A normálási faktor az ablakfüggvény négyzetének a mérési időablakra vett integrálja. Láttuk, hogy ez téglalap ablakkal számolva \setbox0\hbox{$T$}% \message{//depth:\the\dp0//}% \box0%. Most vizsgáljuk meg a zajméréshez leggyakrabban használt Hanning-ablakot, melynek ablakfüggvénye:

\[W(t)=2\sin^2\left(\frac{t\pi}{T}\right).\]

A zajsűrűség Hanning-ablakot használva:

\[s_0=\frac{4\langle|v_W(\omega)|^2\rangle}{3T}.\]

Diszkrét mérési pontok esetén a következőképpen módosul a kifejezés:

\[s_0\approx\dfrac{2}{\sum_{n=0}^{N-1}W^2(n\Delta t)\Delta t}|\sum_{n=0}^{N-1}W(n\Delta t) V(N\Delta t)e^{-i\omega n\Delta t}\Delta t|^2.\]

Hanning-ablakkal számolva:

\[s_0\approx\frac{4\Delta t}{3N}|\sum_{n=0}^{N-1}W(n\Delta t)\Delta V(n\Delta t)e^{-i\omega n\Delta t}|^2.\]

Így frekvenciafüggetlen zajsűrűségek esetén sikerült analitikusan meghatározni, hogy a mérési pontokból miképpen számolható ki a zaj spektrális sűrűsége.

Termikus zaj


Egy kristályban lévő elektronok termikus fluktuációi miatt külső feszültség nélkül az áram időátlaga nulla, azonban az adott időpillanatokban véletlen irányokba mutató áramokból áramfluktuációk alakulnak ki. Ezt, az elektronok termikus fluktuációjából adódó zajt nevezzük termikus zajnak. Jellegét tekintve ez a zaj egy fehér zaj, azaz a zajsűrűsége független a frekvenciától.

A Drude-modell értelmében az elektronok elektromos tér által nyert impulzusát a kristályráccsal történő ütközése (szóródás, rácsrezgés keltése) során veszíti el. A modell értelmében a fajlagos vezetőképesség (\setbox0\hbox{$\sigma$}% \message{//depth:\the\dp0//}% \box0%) arányos az elektronok átlagos sűrűségével (\setbox0\hbox{$n$}% \message{//depth:\the\dp0//}% \box0%), töltésük négyzetével (\setbox0\hbox{$e$}% \message{//depth:\the\dp0//}% \box0%) és az átlagosan két ütközés között eltelt idővel (\setbox0\hbox{$\tau$}% \message{//depth:\the\dp0//}% \box0%), továbbá fordítottan arányos az elektron tömegével (\setbox0\hbox{$m$}% \message{//depth:\the\dp0//}% \box0%):

\[\sigma=\frac{ne^2\tau}{m}.\]

Most vizsgáljuk meg ugyanezen rendszer áramzaját külső feszültség nélkül. (Tételezzük fel, hogy a vizsgált vezető mintánk két végét egy árammérővel összekötve mérjük az áramfluktuációkat.) A \setbox0\hbox{$\tau$}% \message{//depth:\the\dp0//}% \box0% karakterisztikus időnként elszenvedett ütközések miatt az áram korrelációja is elvész. Ezt egy \setbox0\hbox{$\tau$}% \message{//depth:\the\dp0//}% \box0% időállandóval jellemzett exponenciálisan lecsengő korrelációs függvénnyel modellezzük:

\[C(\Delta t)=C_0e^{-\frac{|\Delta t|}{\tau}}.\]

Ilyen korrelációs függvény esetén a zajsűrűségre a következő érték adódik:

\[s_I(\omega)=2c(\omega)=2\int_{-\infty}^{\infty}dte^{-i\omega t}C_0e^{-\frac{|\Delta t|}{\tau}}=\frac{4C_0\tau}{1+\omega^2\tau^2}.\]

A feladat tehát nem más, mint a \setbox0\hbox{$C_0$}% \message{//depth:\the\dp0//}% \box0% érték, azaz az áram teljes szórásnégyzetének a meghatározása a Drude-modell segítségével. Vizsgáljunk egy \setbox0\hbox{$A$}% \message{//depth:\the\dp0//}% \box0% keresztmetszetű, \setbox0\hbox{$L$}% \message{//depth:\the\dp0//}% \box0% hosszúságú téglalap alakú dobozt. Ebben egy elektron árama a következőképpen fejezhető ki:

\[I_1=-Anev_x(t)=-\frac{e}{L}v_x(t),\]

ahol \setbox0\hbox{$v_x(t)$}% \message{//depth:\the\dp0//}% \box0% a részecske \setbox0\hbox{$x$}% \message{//depth:\the\dp0//}% \box0% irányú, azaz a doboz hosszanti irányába mutató sebessége adott \setbox0\hbox{$t$}% \message{//depth:\the\dp0//}% \box0% időpillanatban. Innen következik egy elektron áram-áram korrelációs függvénye:

\[\langle (\Delta I_1(t))^2\rangle=\frac{e^2}{L^2}\langle v_x^2(t)^2\rangle.\]

N független elektronra a szórásnégyzetek összeadódnak:

\[\langle (\Delta I(t))^2\rangle=C_0=N\langle (\Delta I_1(t))^2\rangle=N\frac{e^2}{L^2}\langle v_x^2(t)\rangle.\]

Bővítsük a törtet \setbox0\hbox{$\tau Am$}% \message{//depth:\the\dp0//}% \box0%-mel, hogy észrevegyük benne a Drude-modell vezetőképességét:

\[C_0=\frac{N}{AL}\frac{e^2\tau}{m}\frac{A}{L}\frac{1}{\tau}m\langle v_x^2(t)\rangle =\frac{G}{\tau}m\langle v_x^2(t)\rangle.\]

Tovább egyszerűsödik a képlet a klasszikus ekvipartíciós tételt alkalmazva:

\[C_0=\frac{G}{\tau}k_BT.\]

Ezt felhasználva, továbbá élve azzal a közelítéssel, hogy a frekvencia sokkal kisebb az átlagos ütözkési idő reciprokánál, adódik a termikus áramzaj spektrális sűrűsége:

\[s_I(\omega)=\frac{4C_0\tau}{1+\omega^2\tau^2}=\frac{4k_BTG}{1+\omega^2\tau^2}\approx 4k_BTG.\]

Felhasználva, hogy \setbox0\hbox{$s_V=R^2s_I$}% \message{//depth:\the\dp0//}% \box0% a feszültségzaj spektrális sűrűsége:

\[s_V=4k_BTR.\]

Ezzel számítással is beláttuk a termikus zaj fehér zaj jellegét, továbbá azt is látjuk, hogy a zaj mértéke arányos a hőmérséklettel és az ellenállás nagyságával.

Meg kell jegyezni továbbá, hogy a számolás két helyen sem volt teljesen korrekt. Egyrészt csak a Fermi-energia körüli \setbox0\hbox{$k_BT$}% \message{//depth:\the\dp0//}% \box0% energiasávban lévő elektronok tekinthetők valamilyen szinten függetlennek, a mélyebb energiákon lévő elektronállapotok mindegyike teljesen betöltött, így az áramfluktuációjuk zérus. Másrészt a Fermi-Dirac-eloszlást követő elektronok sebességét nem írhatjuk le az ekvipartíció tétellel, hanem a zajhoz járulékot adó részlegesen betöltött elektronállapotok alapvetően a Fermi-sebességgel mozognak. Ez a két hiba kompenzálja egymást, és korrekt szilárdtestfizikai számításokkal is a végeredményben kapott képlettel azonos eredményre jutunk.

Poisson zaj


A Poisson-zaj az elektronok diszkrét töltéséből fakadó áramzajnak, azaz a sörétzaj egy speciális fajtája. A Poisson-zaj specialitását az egymástól független elektronok detektálása jelenti. Az elnevezése is innen származik, hiszen ahogyan a klasszikus példában a véletlen időközönként elsütött puskagolyók, úgy az egymástól függetlenül kilépő elektronok száma is Poisson-eloszlást mutat. A Poisson-eloszlás sörétzaj szempontjából fontos tulajdonsága, hogy a várható értéke megegyezik a szórásnégyzetével:

\[\langle N\rangle=\langle (\Delta N)^2\rangle.\]

Tegyük fel, hogy \setbox0\hbox{$\Delta t$}% \message{//depth:\the\dp0//}% \box0% idő alatt átlagosan N elektron lép ki az egyik elektródából. Ekkor az áram várható értéke a következőképpen definiálható:

\[\langle I\rangle=\langle N\rangle\frac{e}{\Delta t}.\]

Az áram szórásnégyzete pedig felhasználva a várható érték és szórásnégyzet közötti egyenlőséget:

\[\langle (\Delta I)^2\rangle=\langle (\Delta N)^2\rangle\left(\dfrac{e}{\Delta t}\right)^2=\langle I\rangle\frac{e}{\Delta t}.\]

Ebben az esetben azonban a \setbox0\hbox{$\Delta t$}% \message{//depth:\the\dp0//}% \box0% nem a fentebb definiált két diszkrét mintavételezés között eltelt időt jelenti, hanem azt az időintervallumot, amelyen belül az áram várható értéke \setbox0\hbox{$\langle I\rangle$}% \message{//depth:\the\dp0//}% \box0%, illetve szórásnégyzete \setbox0\hbox{$\langle (\Delta I)^2\rangle$}% \message{//depth:\the\dp0//}% \box0%.

Határozzuk meg ebben az időintervallumban az átlagos áramot. Ehhez a jelünket konvolváljuk össze egy \setbox0\hbox{$G(t)$}% \message{//depth:\the\dp0//}% \box0% függvénnyel, amely egyenlő \setbox0\hbox{$\dfrac{1}{\Delta t}$}% \message{//depth:\the\dp0//}% \box0%-vel, ha \setbox0\hbox{$|t|<\dfrac{\Delta t}{2}$}% \message{//depth:\the\dp0//}% \box0%, egyébként nulla.

\[I_{out}(t)=\int_{-\infty}^{\infty}dt'I_{in}(t')G(t-t')=\frac{1}{\Delta t}\int_{t-\frac{\Delta t}{2}}^{t+\frac{\Delta t}{2}}dt'I_{in}(t').\]

A konvolúció Fourier-transzformáltja a Fourier-transzformáltak szorzatával egyenlő:

\[i_{out}(\omega)=i_{in}(\omega)g(\omega).\]

Ez az összefüggés azt sejteti, hogy a konvolváló függvény szűrőként viselkedik. Mivel a zaj az áram Fourier-transzformáltjának négyzetével arányos az áramzajok aránya a szűrő karakterisztika Fourier-transzformáltjának négyzete lesz:

\[s_{out}(\omega)=s_{in}(\omega)g^2(\omega).\]

Vizsgáljuk meg a szűrő karakterisztikáját:

\[g(\omega)=\int_{-\infty}^{\infty}dtg(t)e^{-i\omega t}=\frac{1}{\Delta t}\int_{-\frac{\Delta t}{2}}^{+\frac{\Delta t}{2}}dte^{-i\omega t}=\frac{2}{\omega \Delta t}sin\left(\frac{\omega \Delta t}{2}\right).\]

A továbbiakhoz bevezetjük az ekvivalens zaj sávszélesség (Equivalent Noise Bandwidth, ENBW) fogalmát. Az ENBW az a sávszélesség, amelyen fehér zajt feltételezve ugyanakkorának adódik az áram szórásnégyzete, mint a szűrőnkön mérve.

\[\langle (\Delta I)^2\rangle=\frac{1}{2\pi}\int_{0}^{\infty}d\omega s_{out}(\omega)=\frac{s_0}{\pi \Delta t}\int_{0}^{\infty}\frac{sin^2(x)}{x^2}=\frac{s_0}{2\Delta t}=s_0f_{max}.\]

Azaz az ekvivalens zaj sávszélesség:

\[ENBW=\frac{1}{2\Delta t}.\]

Ez az összefüggés érezhetően összecseng a Nquist-Shannon mintavételezési törvénnyel. Ha egy \setbox0\hbox{$\Delta t$}% \message{//depth:\the\dp0//}% \box0% időablakra átlagolunk egy jelet, az a konvolúció egy olyan szűrőként viselkedik a frekvenciatérben, melyen keresztül fehér zajt mérve az áram szórásnégyzet egyenlő lesz egy \setbox0\hbox{$\dfrac{1}{2\Delta t}$}% \message{//depth:\the\dp0//}% \box0% tökéletes sávszűrőn mért szórásnégyzettel.

Most vizsgáljuk az áram szórásnégyzetét:

\[\langle (\Delta I)^2\rangle=\langle I\rangle\dfrac{e}{\Delta t}=\int_{0}^{\infty}s_I(f)g(f)df=\int_{0}^{\frac{1}{2\Delta t}}s_I(f)df.\]

Ez bármely \setbox0\hbox{$\Delta t$}% \message{//depth:\the\dp0//}% \box0% szélességű átlagolási időre csak akkor teljesülhet, ha a zajsűrűség független a frekvenciától, értéke pedig:

\[s_I(f)=2e\langle I\rangle.\]

Aliasing jelenség


A DFT a mért jel spektrumát \setbox0\hbox{$0$}% \message{//depth:\the\dp0//}% \box0% és \setbox0\hbox{$\omega_{max}=\frac{2\pi}{2\Delta t}$}% \message{//depth:\the\dp0//}% \box0% frekvenciák közötti diszkrét pontokon értékeli ki. Felmerül a kérdés, hogy egy \setbox0\hbox{$\omega_{max}$}% \message{//depth:\the\dp0//}% \box0%-nál nagyobb frekvenciakomponenseket tartalmazó jelnél mi történik a magas frekvenciakomponensekkel.

Vizsgáljunk egy \setbox0\hbox{$\omega_0$}% \message{//depth:\the\dp0//}% \box0% frekvenciájú tiszta szinuszos jelet:

\[F(t)=A_0e^{i\omega_0t}.\]

Vegyük ennek a jelnek a Fourier-transzformáltját:

\[f(\omega)=\int_{-\infty}^{\infty}A_0e^{i\omega_0t}e^{-i\omega t}dt=A_02\pi\delta(\omega-\omega_0).\]

A jel Fourier-transzformáltja tehát a várakozásnak megfelelően egy Dirac-delta. Most nézzük meg, hogy a DFT számolása során hogyan változik a spektrum meghatározása a gyakorlatban.

\[DFT(\omega_k)=\sum_{n=0}^{N-1}A_0e^{i\omega_0n\Delta t}W(n\Delta t)e^{-i\omega_k n\Delta t}\Delta t.\]

Vegyük észre, hogy

\[DFT(\omega_k)=DFT\left(\omega_k+\frac{2\pi  m}{\Delta t}\right),\]

ahol \setbox0\hbox{$m$}% \message{//depth:\the\dp0//}% \box0% egész szám. Másrészt

\[DFT(\omega_k)=DFT^*(-\omega_k),\]

azaz

\[\left|DFT(\omega_k)\right|=|DFT(-\omega_k)|.\]

Így belátható, hogy tetszőleges magas körfrekvenciájú jelet úgy látunk, mint ha az a \setbox0\hbox{$\left[0,\omega_{max}=\dfrac{2\pi}{2\Delta t}\right]$}% \message{//depth:\the\dp0//}% \box0% tartományban lévő jel lenne a saját \setbox0\hbox{$A_0$}% \message{//depth:\the\dp0//}% \box0% amplitúdójával.

A \setbox0\hbox{$(2m)\dfrac{2\pi}{2\Delta t}<\omega<(2m+1)\dfrac{2\pi}{2\Delta t},$}% \message{//depth:\the\dp0//}% \box0% \setbox0\hbox{$m=1,2,\ldots$}% \message{//depth:\the\dp0//}% \box0% frekvenciatartományokban lévő komponensek az \setbox0\hbox{$\omega_m=\omega-\dfrac{2\pi m}{\Delta t}$}% \message{//depth:\the\dp0//}% \box0% frekvencián jelennek meg a spektrumban, míg a \setbox0\hbox{$(2m-1)\dfrac{2\pi}{2\Delta t}<\omega<(2m)\dfrac{2\pi}{2\Delta t},$}% \message{//depth:\the\dp0//}% \box0% \setbox0\hbox{$m=1,2,\ldots$}% \message{//depth:\the\dp0//}% \box0% tartományokban lévők a \setbox0\hbox{$\omega_m=\dfrac{2\pi m}{\Delta t}-\omega$}% \message{//depth:\the\dp0//}% \box0% frekvencián.

A fentiekkel beláttuk, hogy minden \setbox0\hbox{$\omega_{max}$}% \message{//depth:\the\dp0//}% \box0%-nál nagyobb frekvenciatag megjelenik a DFT-t követően a diszkrét frekvenciákat tartalmazó \setbox0\hbox{$[0,\omega_{max}]$}% \message{//depth:\the\dp0//}% \box0% tartományban. Gondoljuk át, hogy a végeredményben kapott DFT spektrum \setbox0\hbox{$\omega_k$}% \message{//depth:\the\dp0//}% \box0% pontjában összeadódnak-e a \setbox0\hbox{$DFT_m(\omega_k)$}% \message{//depth:\the\dp0//}% \box0% tagok amplitúdói \setbox0\hbox{$m=0,1,2,\ldots$}% \message{//depth:\the\dp0//}% \box0% értékekre. Mivel ezen frekvenciatagok fázisai nem azonosak, így a különböző \setbox0\hbox{$m$}% \message{//depth:\the\dp0//}% \box0% értékekhez tartozó amplitúdók nem adhatóak össze.

Később látni fogjuk, hogy a zajsűrűség számításához a DFT abszolútérték négyzetére van szükség, amely azonban a véletlenszerű, egymástól független fázisok miatt egyenlő a különböző \setbox0\hbox{$m$}% \message{//depth:\the\dp0//}% \box0% értékekre vett amplitúdók négyzeteivel:

\[|DFT(\omega_k)|^2=\left(\sum_{m}^{} DFT_m(\omega_k)\right)\left(\sum_{m}^{} DFT^*_m(\omega_k)\right)=\sum_{m}^{}|DFT_m(\omega_k)|^2.\]


A leképezés leginkább úgy képzelhető el, mint a magas frekvenciás tartományok \setbox0\hbox{$|DFT|^2$}% \message{//depth:\the\dp0//}% \box0%-einek a \setbox0\hbox{$[0,\omega_{max}]$}% \message{//depth:\the\dp0//}% \box0% intervallumra történő visszahajtogatása. Az \setbox0\hbox{$\omega_{max}$}% \message{//depth:\the\dp0//}% \box0% fölötti részt visszahajtjuk \setbox0\hbox{$[-\infty; \omega_{max}]$}% \message{//depth:\the\dp0//}% \box0% intervallumra. Majd ezt a visszahajtott spektrumot 0 \setbox0\hbox{$Hz$}% \message{//depth:\the\dp0//}% \box0%-nél előrehajtjuk, stb. Ezt mindaddig folytatjuk, amíg minden jelentős frekvenciakomponenst be nem hajtottunk a \setbox0\hbox{$[0,\omega_{max}]$}% \message{//depth:\the\dp0//}% \box0% intervallumba.


Beláttuk, hogy a magas komponensű jelek fals spektrumot okozhatnak, így gondoskodni kell kiszűrésükről egy aluláteresztő szűrővel. Ezt a jelenséget aliasingnak, azaz magas frekvenciájú komponensek beszűrődésének nevezzük. Fontos megjegyezni, hogy az aliasing kialakulásának az oka a DFT diszkrét mintavételezése, és független a mérési pontok \setbox0\hbox{$N$}% \message{//depth:\the\dp0//}% \box0% számától. A gyakorlatban a legtöbbször az erősítők, vagy a mérőkártyák rendelkeznek anti-aliasing szűrővel, ami a maximális frekvencia fölött rendszerint élesen levágja a spektrumot.