Megpróbálom megoldani a következő problémát, miközben még mindig nincs határozott elképzelésem arról, hogy mit jelent a „frekvenciafelbontás”:

Tegyük fel, hogy folytonos időjelet veszünk mintából, amelynek Ts = 1/2000 mintavételi periódusa van, majd használunk egy 1000 hosszúságú ablakot a kapott diszkrét időjelre. Ha egy 2000 pontos DFT segítségével átalakítjuk, mekkora lenne a frekvenciafelbontása?

Tudna valaki segíteni ebben?

Megjegyzések

  • Potenciális diagramfelbontást szeretne interpolációval, csúcshely becslési felbontást S / N esetén, eredménytartály-elválasztást vagy csúcselválasztási felbontást elválasztási kritériumokkal? Mindezek különböző frekvenciafelbontásokat eredményeznek ugyanolyan hosszúságú DFT esetén.
  • @ hotpaw2 Érdekelne, ha beszélhetne ezekről a felbontásokról ebben vagy egy másik informatív kérdésben.

Válasz

Szerkesztés:

Rájöttem, hogy a " frekvenciafelbontás " alábbi definícióm teljesen teljes téves (valamint az OP kérdése). A Frekvenciafelbontás mennyire hasonlít az ablak függvényének a frekvenciatérben a Dirac delta függvényéhez. Ennek oka, hogy az ablak és az időtartományban lévő jel szorzata konvolúcióvá válik a frekvenciatartományban ( és a Dirac delta függvényes konvolúció egy olyan mintavétel, amely tökéletes frekvenciafelbontást eredményezne.) Minél zsírosabb a mainlobe (annak szórásával számszerűsítve), és minél nagyobbak az oldalsó görbék, annál rosszabb a frekvenciafelbontás. Ezenkívül az Időfelbontás számszerűsíthető az ablak függvény varianciája az időtartományban.


A frekvenciafelbontás nem bináris felbontás / szélesség. Az alábbi grafikonon vegye figyelembe, hogy a lebenyek nem kerülnek közelebb (frekvenciafelbontás), annak ellenére, hogy a tároló szélessége csökken.

Hitel: Dan Boschen

A frekvenciafelbontás inkább a téglalap alakú függvény (azaz a sinc függvény) Fourier-transzformációjának tulajdonsága.

A Fourier-transzformációkkal való együttműködéshez ablakot kell kapnunk (még elméletileg is). Ennek eredményeként mindig a $ f (t) w (t) $ -val dolgozunk, nem pedig a $ f (t ) Maga a $ (itt $ w (t) $ egy téglalap alakú függvény). A konvolúciós tétel szerint az ablakos függvény Fourier-transzformációja mindig $ \ hat {f} $ konvolúciója $ \ kalap {w} = $ sinc. Különösen, ha a $ f $ szinuszos, a $ \ hat {f} $ egy Dirac delta függvény lesz, és a konvolúció csak egy sinc függvény mintavétele lesz. Így időszakosan teljesen elveszítjük a frekvenciákat ablakoláskor, ennek a veszteségnek a periodicitása a frekvenciafelbontás .

Mivel az ablakos függvényeknél a DTFT a CTFT periodikus közelítése, megszerzi ezeket a tulajdonságokat is.

A zavar azért merül fel, mert amikor nem adunk nullákat a DFT-hez (azaz csak minta $ f (t) w (t) $ ahol $ w (t) = 1 $ ), a tároló szélessége megegyezik a Frekvencia felbontással.

Ugyanakkor nullákat is kitölthetünk (azaz a mintát is $ f (t) w (t) $ ahol $ w (t) = 0 $ ), és ez azt eredményezi, hogy a DTF jobban interpolálja a span class = “math-container”> $ f (t) w (t) $ . Konferencia az első grafikonnal.


Hogy megtudja, miért a téglalap alakú függvény Fourier-transzformációja a a sinc funkció , nézze meg ezt a videót és vegye fontolóra a szinuszos függvények felcsavarodását (bár ez eléggé érintett)


Az OP példájának megválaszolásához a bin felbontása $$ \ frac {F_s} {N} = \ frac {2000} {2000} = 1 $$ ahol $ F_s = 2000 $ Hz a mintavételi arány, és $ N $ a DFT mérete.

A frekvenciafelbontás az, ami a bin felbontása lenne, ha csak mintát vennénk az ablakban (nincs nulla kitöltés)

$$ \ frac {F_s} {M} = \ frac {1} {T} = 2 $$ ahol $ M $ az ablakban található minták száma, $ T $ a minta időtartama, és $ F_s = M / T $ .

Megjegyzések

  • Kedves válasz Tom.Ha nem világos, akkor gyakran nem ' nem használunk téglalap alakú ablakot, hanem más ablakokat, amelyek elkeskenyednek, és amelyek jelentősen csökkentik az oldalsó görbéket (javítják a dinamikatartományt) a rontás rovására. a frekvenciafelbontás tovább. Az egyik kedvenc klasszikus írásom erről és a DFT alkalmazásairól általában Fred Harris. Szerintem ' nagyon élvezed, ha ' még nem láttad: web.mit.edu/xiphmont/Public/windows.pdf
  • @TomHuntington Szép, kár, hogy ' kétszer is fel tudok szavazni!
  • A @TomHuntington Wikipedia nyilván nem tud ' nem tud a képleteimről vagy a technikáimról. Továbbra is nehézségeim vannak az intrabin felbontással (a zaj és az egyenletek érzékenysége miatt), de a közeli frekvenciák feloldhatók iteratív becsléssel és eltávolítással. A nagy hang eltávolításakor a kisebb becsülhető. Ha eltávolítja a kicsi hangot, jobban olvashatja a nagyot. És így tovább, akár többféle hanggal is. Bármilyen ablak bonyolítja a matematikát.
  • Ha két közel azonos amplitúdójú, de nagyon frekvenciájú sinusoidája van, használhatja a beat jelenséget az időtartományban. A jel látszólagos frekvenciája (nulla keresztezéssel) a két frekvencia átlaga, a burok frekvenciája (ha teljes ciklust veszünk, pl. Két lebeny) a fele a frekvenciák különbségének.
  • Ezenkívül a felbontás meghatározza az Ön pontosságát a mérésben. Nem mond semmit a pontosságról.

Válasz

Kicsit attól függ, hogy mit akar elérni.

Ha $ N $ hosszúságú FFT-t hajt végre a mintából vett mintából $ F_s $ , akkor sokan azt mondanák, hogy a frekvencia felbontása $ \ frac {F_s} {N} $ . Hogy ez helyes-e vagy sem, az valóban attól függ, hogy pontosan hogyan határozza meg a frekvenciafelbontást és mit tervez vele csinálni.

Ami valójában az történik, hogy mintavételezéssel mintavételez egy frekvenciatartomány-függvényt $ \ frac {F_s} {N} $ intervalluma. Amint kiválaszt egy FFT méretet, mindkét tartományban mintavételt végez, a mintavételi intervallumok $ \ frac {1} {F_s} $ időben és $ \ frac {F_s} {N} $ gyakoriságban.

A frekvenciatartomány-mintavételnek ugyanazok a tulajdonságai, követelményei és problémái vannak, mint az időtartomány-mintavételnek, álnévhez juthat, interpolálhat, feltételezhető a periodicitás a másik tartományban, stb.

A mintavételi tétel egyszerű alkalmazásával azt állíthatjuk, hogy a jel teljes jellemzéséhez szükséges frekvenciafelbontás egyszerűen a hossza az időtartományban. Ez jól működik olyan jelek esetében, amelyek eredendően időhöz kötöttek, például egy LTI-rendszer impulzus-válasza.

Azonban ez nem praktikus hosszú folyamatos jelek esetén. Ebben az esetben olyan frekvencia felbontást kell választania, amely “elég jó” az alkalmazásához, és ez valóban a készülék követelményeitől és céljától függ. konkrét alkalmazás.

Válasz

A mintavételt $ {T} _ {s} adja = \ frac {1} {2000} $ [sec].
Az ablak hossza 1000 minta.
Mivel az ablak hosszának meg kell egyeznie az adat hosszával, arra következtetünk, hogy az adat hossza 1000 minta ami azt jelenti, hogy a mintavételi idő $ 0,5 $ [sec].

A BFT felbontása a DFT-ben a mintavételi intervallum és a DFT minták, ami ebben az esetben 2000. Ezért a bin felbontása $ \ frac {1} {4000} $ [Hz].

Válasz

Az FFT binw szélessége vagy az újrarendezés felbontása, ahogy szeretem hívni, Fs / N, ahol N az FFT mérete. A tényleges felbontás a használt ablaktól és az ablak hosszától függ.

Például: egy téglalap alakú ablak maximális felbontást, de kevésbé dinamikus tartományt biztosít. Más simább ablakok kisebb felbontást biztosítanak dinamikusabb tartományban vagy az alsó oldalsó lebenyek mellett.

Vélemény, hozzászólás?

Az email címet nem tesszük közzé. A kötelező mezőket * karakterrel jelöltük