Пређи на садржај

Краткорочна Фуријеова трансформација

С Википедије, слободне енциклопедије

Краткорочна Фуријеова трансформација је Фуријеова трансформација која се користи за одређивање синусоидне фреквенције и фазног садржаја локалних одсека сигнала како се током времена мења. [1] У пракси, поступак израчунавања ове трансформације састоји се од дељења дужих временских сигнала на краће сегменте једнаке дужине, а затим засебно израчунати Фуријеову трансформацију на сваком краћем сегменту. Ово открива Фуријеов спектар на сваком краћем сегменту. Тада обично поставља променљиви спектар као функцију времена.

Пример кратких времена Фуријеове трансформације које се користе за одређивање времена удара од аудио сигнала

Напредна краткорочна Фуријеова трансформација

[уреди | уреди извор]

Временски континуална краткорочна Фуријеова трансформација

[уреди | уреди извор]

Једноставно, у случају континуираног времена, функција која се трансформише множи се прозорском функцијом која је само кратко време не постоји. Фоуриерова трансформација (једнодимензионална функција) резултирајућег сигнала узима се док се прозор помиче дуж временске оси, резултирајући дводимензионалним приказом сигнала. Математички, ово се пише као:

где јесте функција прозора, обично Ханов прозор или Гаусов прозор који је центриран око нуле а јесте сигнал који треба трансформисати (имајте на уму разлику између функције прозора и фреквенције ). у суштини је Фуријеова трансформација , сложена функција која представља фазу и јачину сигнала током времена и фреквенције. Често се одмотавање фаза употребљава дуж једне или обе временске осе, и осне фреквенције, , да се сузбије било какав прекид скока фазног резултата краткорочне Фуријеове трансформације. Временски индекс обично се сматра „ спорим “ временом и обично се не изражава у тако високој резолуцији као време .

Краткорочна Фуријеова трансформација у дискретном времену

[уреди | уреди извор]

У дискретном временском случају, подаци који се трансформишу могу се рашчланити на делове или оквире (који се обично преклапају, да би се смањили артефакти на граници). Сваки комад се трансформише Фуријеом трансформацијом, а сложен резултат се додаје у матрицу која бележи величину и фазу за сваку тачку у времену и фреквенцији. То се може изразити као:

слично, са сигналом x[n] и прозором w[n]. У овом случају, m је дискретан и ω је континуaлан, али у већини типичних примена СТФТ се изводи на рачунару помоћу брзе Фуријеове трансформације, тако да су обе променљиве дискретне и квантизоване .

Величина у квадратури краткорочне Фуријеове трансформациеј даје спектрограмски приказ спектралне густоће снаге функције:

Погледајте и модификовану дискретну трансформацију косинуса, која је такође Фуријеова трансформација која користи прозоре који се преклапају.

Клизна дискретна Фуријеова трансформација

[уреди | уреди извор]

Ако је пожељан само мали број ω или ако се жели да се краткорочна Фуријеова трансформација процењује за сваки помак m прозора, тада се она може ефикасније проценити користећи клизни ДФТ алгоритам. [2]

Инверзна краткорочна Фуријеова трансформација

[уреди | уреди извор]

Краткорочна Фуријеова трансформација је инвертибилна, то јест, изворни сигнал се може повратити из трансформације помоћу инверзне краткорочне Фуријеове трансформације. Најприхваћенији начин претварања краткорочне Фуријеове трансформације је коришћењем методе преклапања и додавања, која такође омогућава измене у сложеном спектру краткорочне Фуријеове трансформације. Ово чини свестраним методе обраде сигнала [3] и назива се преклапањем и додавањем модификација.

Временски континуирана краткорочна Фуријеова трансформација

[уреди | уреди извор]

С обзиром на ширину и дефиницију прозорске функције w(t), у почетку захтевамо да се површина функције прозора умањи тако да

Континуална Фуријеова трансформација је

Ако заменимо х(t) одозго:

Мењање редоследа интеграције:

Дакле, Фуријеова трансформација може се посматрати као врста фазног кохерентног сума свих СТФТ-ова к ( т ). Пошто је инверзна Фоуриерова трансформација

тада се x(t) може повратити из X(τ, ω) као

или

Може се видети, да је у поређењу са горенаведеним „таласом“ прозора од х(t)

инверзна Фуријеова трансформација X(τ, ω) за τ фиксна.

Питања резолуције

[уреди | уреди извор]

Једна од замки краткорочне Фуријеове трансформације је та што има фиксну резолуцију. Ширина функције прозора односи се на то како је сигнал представљен - он одређује да ли постоји добра фреквенциона резолуција (фреквентне компоненте близу међусобно се могу раздвојити) или добра временска резолуција (време у коме се фреквенције мењају). Широки прозор даје бољу резолуцију фреквенције, али лошу временску резолуцију. Ужи прозор даје добру временску резолуцију, али слабу фреквенцију. То се назива ускопојасна и широкопојасна трансформација.

Поређење резолуције краткорочне Фуријеове трансформације . Лева има бољу временску резолуцију, а десна бољу резолуцију фреквенције.

Ово је један од разлога за стварање валне трансформације и мултиресолуцијске анализе, који могу дати добру временску резолуцију за високофреквентне догађаје и добру фреквенцијску резолуцију за догађаје са ниским фреквенцијама, комбинација која је најприкладнија за многе стварне сигнале.

Ово својство је повезано са Хеисенберговим принципом неизвесности, али није директно - видети ограничење Габора за расправу. Производ стандардног одступања у времену и фреквенцији је ограничен. Граница принципа несигурности (најбоље истовремено решење оба) се постиже Гауссовом функцијом прозора, јер Гауссиан минимизира принцип Фоуриерове несигурности . То се назива Габорова трансформација (и модификацијама за мултирезолуцију постаје Морлет- ова валутна трансформација).

Може се размотрити СТФТ за промену величине прозора као дводимензионална домена (време и фреквенција), као што је илустровано у доњем примеру, који се може израчунати променом величине прозора. Међутим, ово више није строго представљање временске фреквенције - кернел није константан током целог сигнала.

Користите следеће узорке сигнала који су састављени од скупа четири таласа синусоидног облика спојених у низу. Сваки таласни облик састоји се само од једне од четири фреквенције (10, 25, 50, 100 Hz ). Дефиниција је:

Затим се узоркује на 400 Hz. Добијени су следећи спектрограми:

25 мс прозор
125 мс прозор
Прозор од 375 мс
1000 мс прозора

Прозор од 25 мс омогућава нам да идентификујемо тачно време у коме се сигнали мењају, али прецизне фреквенције је тешко препознати. На другом крају скале, прозор од 1000 мс омогућава да се фреквенције тачно виде, али време између промена фреквенције је замагљено.

Објашњење

[уреди | уреди извор]

То се такође може објаснити референцом на фреквенцију узорковања и Никуист фреквенцију .

Погледајте прозор Н узорака из произвољног сигнала реалне вредности брзином узорковања ф с . Преузимање Фоуриерове трансформације даје Н сложене коефицијенте. Од ових коефицијената корисна је само половина (последњи Н / 2 је сложен коњугат прве Н / 2 обрнутим редоследом, јер је то стварни сигнал).

Ови Н / 2 коефицијенти представљају фреквенције од 0 до ф с / 2 (Никуист), а два узастопна коефицијента међусобно су растављена ф с / Н Хз.

Да бисте повећали фреквенцијску разлучивост прозора, размак коефицијената треба смањити. Постоје само две променљиве, али смањење ф с (и одржавање Н константног) узроковаће да се повећа величина прозора - јер сада има мање узорака по јединици времена. Друга алтернатива је повећање Н, али то опет узрокује повећање величине прозора. Стога сваки покушај повећања резолуције фреквенције узрокује већу величину прозора, а самим тим и смањење временске разлучивости, и обрнуто.

Рејлова фреквенција

[уреди | уреди извор]

Како је Никуист-ова фреквенција ограничење у максималној фреквенцији која се може смислено анализирати, тако је и Раилеигх-ова фреквенција ограничење на минималну фреквенцију.

Раилеигх-ова фреквенција је минимална фреквенција која се може ријешити временским временским оквиром. [4] [5]

С обзиром на временски прозор у трајању од Τ секунде, најмања фреквенција која се може разрешити је 1 / Τ   Хз.

Раилеигх-ова фреквенција важан је аспект у примјени краткотрајне Фоуриер-ове трансформације (СТФТ), као и било које друге методе хармоничке анализе на сигналу коначне дужине записа. [6] [7]

Апликација

[уреди | уреди извор]
Краткорочна Фуријеова трансформација се користи за анализу аудио сигнала током времена

Краткорочне Фуријеове трансформације, као и стандардне Фуријеове трансформације и други алати се често користе за анализу музике. Спектрограм, на пример, може да приказује фреквенцију на водоравној оси, а најниже фреквенције лево, а највише са десне стране. Висина сваке траке (повећана бојом) представља амплитуду фреквенција унутар тог опсега. Димензија дубине представља време, где је свака нова трака била засебна посебна трансформација. Аудио инжењери користе ову врсту визуелног материјала за прибављање информација о аудио узорку, на пример, за проналажење фреквенција специфичних звукова (нарочито када се користе са већом резолуцијом фреквенције) или за проналажење фреквенција које могу бити мање или више резонантне у простору у коме сигнал јесте снимљен. Ове информације се могу користити за изједначавање или подешавање других аудио ефеката.

Имплементација

[уреди | уреди извор]

Оригинална функција

Претварање у дискретни облик:

Претпостављајући претходно

Тада можемо уписати оригиналну функцију у

Директна примена

[уреди | уреди извор]

Ограничења

[уреди | уреди извор]

а. Никвистов критеријум (избегавање ефекта алијасинга):

, где јесте ширина опсега од

Метода заснована на брзој Фуријеовој трансформацији

[уреди | уреди извор]

Ограничење

[уреди | уреди извор]

а. , где је цели број

б.

в. Никвистов критеријум (избегавање ефекта алијасинга):

, је ширина опсега од

Рекурзивна метода

[уреди | уреди извор]

Ограничење

[уреди | уреди извор]

а. , где је цели број

б.

в. Никвистов критеријум (избегавање ефекта алијасинга):

, је ширина опсега од

г. Само за имплементацију правоуглог облика краткорочне Фуријеове трансформације

Правоугли прозор намеће ограничење

Замена даје:

Смена променљиве n-1 за n :

Израчунавање за N-тачака брзе Фуријеове трансформације:

где је

Примена рекурзивне формуле за израчунавање

Генерализација дискретне Фуријеове трансформације

[уреди | уреди извор]

Ограничење

[уреди | уреди извор]

па је

Поређење имплементација

[уреди | уреди извор]
Метод Сложеност
Директна примена
Заснована на брзој Фуријеовој трансформацији
Рекурзија
Генерализација дискретне Фуријеове трансформације
  • Процена спектралне густине
  • Временско-фреквентна репрезентација
  • Метода преноса

Остале временске фреквенције:

  • Функција дистрибуције у облику конуса
  • Константно-К трансформација
  • Фракционарна Фуријеова трансформација
  • Габорова трансформација
  • Њувландова трансформација
  • С трансформација
  • Таласна трансформација
  • Чирплет трансформација

Референце

[уреди | уреди извор]
  1. ^ Sejdić E.; Djurović I.; Jiang J. (2009). „Time-frequency feature representation using energy concentration: An overview of recent advances”. Digital Signal Processing. 19 (1): 153—183. doi:10.1016/j.dsp.2007.12.004. 
  2. ^ E. Jacobsen and R. Lyons, The sliding DFT, Signal Processing Magazine vol. 20, issue 2, pp. 74–80 (March 2003).
  3. ^ Jont B. Allen (јун 1977). „Short Time Spectral Analysis, Synthesis, and Modification by Discrete Fourier Transform”. IEEE Transactions on Acoustics, Speech, and Signal Processing. ASSP-25: 235—238. 
  4. ^ „Архивирана копија” (PDF). Архивирано из оригинала (PDF) 04. 03. 2016. г. Приступљено 14. 12. 2019. 
  5. ^ „What does "padding not sufficient for requested frequency resolution" mean? - FieldTrip toolbox”. 
  6. ^ „Biased competition through variations in amplitude of gamma-oscillations”. J Comput Neurosci. 25 (1): 89—107. 2008. PMC 2441488Слободан приступ. PMID 18293071. doi:10.1007/s10827-007-0066-2. 
  7. ^ Wingerden, Marijn van; Vinck, Martin; Lankelma, Jan; Pennartz, Cyriel M. A. (19. 5. 2010). „Theta-Band Phase Locking of Orbitofrontal Neurons during Reward Expectancy”. Journal of Neuroscience (на језику: енглески). 30 (20): 7078—7087. ISSN 0270-6474. PMID 20484650. doi:10.1523/JNEUROSCI.3860-09.2010. 

Спољашње везе

[уреди | уреди извор]