Hilbertova transformace - Hilbert transform

V matematice a při zpracování signálu je Hilbertova transformace specifický lineární operátor, který přebírá funkci, u ( t ) skutečné proměnné a vytváří další funkci skutečné proměnné H ( u ) ( t ) . Tento lineární operátor je dán konvolucí s funkcí (viz § Definice ). Hilbertova transformace má obzvláště jednoduchou reprezentaci ve frekvenční oblasti : Přináší fázový posun ± 90 ° ( π2  radiány) každé frekvenční složce funkce, znaménko posunu v závislosti na znaménku frekvence (viz § Vztah s Fourierovou transformací ). Hilbertova transformace je důležitá při zpracování signálu, kde je součástí analytické reprezentace signálu s reálnou hodnotou u ( t ) . Hilbertovu transformaci poprvé představil David Hilbert v tomto prostředí, aby vyřešil speciální případ problému Riemann – Hilbert pro analytické funkce.

Definice

Hilbertovy transformace u si lze představit jako konvoluce o u ( t ) se funkce h ( t ) = 1/ π t, známý jako jádro Cauchy . Protože 1 / t není integrovatelný napříč t = 0 , integrál definující konvoluci nemusí vždy konvergovat. Místo toho je Hilbertova transformace definována pomocí hlavní hodnoty Cauchy (zde zde označena pv ). Hilbertova transformace funkce (nebo signálu) u ( t ) je explicitně dána vztahem

za předpokladu, že tento integrál existuje jako hlavní hodnota. To je právě konvoluce u s temperovaným rozdělením p.v.1/π t. Alternativně lze změnou proměnných integrál hlavní hodnoty zapsat explicitně jako

Když se Hilbertova transformace aplikuje dvakrát za sebou na funkci u , výsledek je záporný u :

za předpokladu, že integrály definující obě iterace konvergují ve vhodném smyslu. Zejména je inverzní transformace −H . Tuto skutečnost lze nejsnadněji zjistit zvážením účinku Hilbertovy transformace na Fourierovu transformaci u ( t ) (viz § Vztah s Fourierovou transformací , níže).

Pro analytickou funkci v horní polorovině popisuje Hilbertova transformace vztah mezi skutečnou částí a imaginární částí hraničních hodnot. To znamená, že pokud f ( z ) je analytický v horní poloviční komplexní rovině { z  : Im { z }> 0} , a u ( t ) = Re { f ( t + 0 · i )} , pak Im { f ( t + 0 · i )} = H ( u ) ( t ) až do aditivní konstanty za předpokladu, že tato Hilbertova transformace existuje.

Zápis

Při zpracování signálu je Hilbertova transformace u ( t ) běžně označována . V matematice je však tento zápis již široce používán k označení Fourierovy transformace u ( t ) . Občas může být Hilbertova transformace označena . Kromě toho mnoho zdrojů definuje Hilbertovu transformaci jako negativum zde definovaného.

Dějiny

Hilbertova transformace vznikla v Hilbertově práci z roku 1905 na problému, který Riemann představoval v souvislosti s analytickými funkcemi, kterému se začalo říkat problém Riemann – Hilbert . Hilbertova práce se zabývala hlavně Hilbertovou transformací pro funkce definované na kruhu. Některé z jeho dřívějších prací souvisejících s diskrétní transformací Hilberta pocházejí z přednášek, které měl v Göttingenu . Výsledky později publikoval Hermann Weyl ve své disertační práci. Schur vylepšil Hilbertovy výsledky o diskrétní Hilbertově transformaci a rozšířil je na integrální případ. Tyto výsledky byly omezeny na mezery L 2 a 2 . V roce 1928, Marcel Riesz prokázáno, že Hilbertova transformace mohou být definovány pro u v ( L p prostoru ) pro 1 < p <∞ , že Hilbertova transformace je ohraničený operátor na dobu 1 < p <∞ , a že podobné výsledky držet pro Hilbertova transformace na kruhu, stejně jako diskrétní Hilbertova transformace. Hilbertova transformace byla motivujícím příkladem pro Antonia Zygmunda a Alberta Calderóna při studiu singulárních integrálů . Jejich vyšetřování hrálo zásadní roli v moderní harmonické analýze. Různé zobecnění Hilbertovy transformace, jako jsou bilineární a trilineární Hilbertovy transformace, jsou dnes stále aktivní oblastí výzkumu.

Vztah k Fourierově transformaci

Hilbertova transformace je operátor multiplikátoru . Násobitel H je σ H ( ω ) = - i sgn ( ω ) , kde sgn je funkce signum . Proto:

kde označuje Fourierovu transformaci . Protože sgn ( x ) = sgn (2 π x ) , vyplývá, že tento výsledek platí pro tři společné definice .

Tím, Eulerovy rovnice ,

Proto H ( u ) ( t ) má účinek posunutí fáze negativních frekvenčních složek u ( t ) o +90 ° ( π2  radiány) a fáze pozitivních frekvenčních složek o -90 °, a i · H ( u ) ( t ) má za následek obnovení pozitivních frekvenčních složek při posunutí negativních frekvenčních složek o dalších +90 °, což má za následek jejich negaci (tj. násobení −1).

Když se Hilbertova transformace aplikuje dvakrát, fáze záporné a kladné frekvenční složky u ( t ) se posunou o +180 ° a -180 °, což jsou ekvivalentní množství. Signál je negován; tj. H (H ( u )) = - u , protože

Tabulka vybraných Hilbertových transformací

V následující tabulce je parametr frekvence skutečný.

Signál
Hilbertova transformace


(viz Dawsonova funkce )
Funkce zinku
Funkce Dirac delta
Charakteristická funkce

Poznámky

  1. ^ Někteří autoři (např. Bracewell) používajíjako definici dopředné transformacenaše −H . Důsledkem je, že pravý sloupec této tabulky bude negován.
  2. ^ a b Hilbertovu transformaci funkcí sin a cos lze definovat tak, že vezmeme hlavní hodnotu integrálu v nekonečnu. Tato definice souhlasí s výsledkem distribuční definice Hilbertovy transformace.

K dispozici je rozsáhlá tabulka Hilbertových transformací. Všimněte si, že Hilbertova transformace konstanty je nulová.

Definiční doména

V žádném případě není zřejmé, že Hilbertova transformace je vůbec dobře definována, protože nesprávný integrál definující ji musí ve vhodném smyslu konvergovat. Hilbertova transformace je však dobře definována pro širokou třídu funkcí, jmenovitě pro 1 < p <∞ .

Přesněji, pokud u je pro 1 < p <∞ , pak mez definující nevlastní integrál

existuje téměř pro každé t . Mezní funkce je také in a je ve skutečnosti mezní hodnotou v průměru nesprávného integrálu. To znamená,

jako ε → 0 v normě L p , stejně jako bodově téměř všude, podle Titchmarshovy věty .

V případě p = 1 se Hilbertova transformace stále bodově sbíhá téměř všude, ale sama o sobě nemusí být integrovatelná, a to ani lokálně. Zejména konvergence průměru v tomto případě obecně nenastává. Hilbertova transformace funkce L 1 se však sbíhá v L 1 -slabá a Hilbertova transformace je omezený operátor od L 1 do L 1, w . (Zejména, protože Hilbertova transformace je také operátorem multiplikátoru na L 2 , Marcinkiewiczova interpolace a argument duality poskytují alternativní důkaz, že H je ohraničeno na L p .)

Vlastnosti

Ohraničenost

Pokud 1 < p <∞ , pak Hilbertova transformace je omezený lineární operátor , což znamená, že existuje konstantní C p taková, že

pro všechny .

Nejlepší konstanta je dána vztahem

Snadný způsob, jak najít to nejlepší pro sílu 2, je prostřednictvím takzvané Cotlarovy identity, která pro všechny skutečné hodnotné f . Stejné nejlepší konstanty platí pro periodickou Hilbertovu transformaci.

Ohraničenost Hilbertovy transformace předpokládá konvergenci operátoru symetrického částečného součtu

na f v .

Anti-self adjpointness

Hilbertova transformace je anti- self adjoint operátor ve vztahu k dualitnímu párování mezi a duálním prostorem , kde p a q jsou Hölderovy konjugáty a 1 < p , q <∞ . Symbolicky,

pro a .

Inverzní transformace

Hilbertova transformace je antiinvoluce , to znamená

za předpokladu, že každá transformace je dobře definována. Protože H zachovává prostor , znamená to zejména, že Hilbertova transformace je invertibilní na , a to

Složitá struktura

Protože H 2 = −I („ I “ je operátor identity ) ve skutečném Banachově prostoru funkcí s reálným hodnocením , definuje Hilbertova transformace lineární komplexní strukturu v tomto Banachově prostoru. Zejména, když p = 2 , Hilbertova transformace dává Hilbertovi prostor skutečných funkcí ve struktuře komplexního Hilbertova prostoru.

(Komplexní) vlastní stavy Hilbertovy transformace připouštějí reprezentace jako holomorfní funkce v horní a dolní polorovině v Hardyho prostoru H 2 podle Paley-Wienerovy věty .

Diferenciace

Formálně je derivátem Hilbertovy transformace Hilbertova transformace derivátu, tj. Tyto dva lineární operátory dojíždějí:

Iterace této identity,

To je přesně pravda, jak je uvedeno za předpokladu, že u a jeho prvních k derivátů patří . Lze to snadno zkontrolovat ve frekvenční oblasti, kde se z diferenciace stává násobení ω .

Konvoluce

Hilbertovu transformaci lze formálně realizovat jako konvoluci s temperovanou distribucí

Formálně tedy

Nicméně, a priori, to může být definována pouze u distribuci kompaktním nosičem . Je možné s tím pracovat poněkud přísně, protože kompaktně podporované funkce (což jsou distribuce a fortiori ) jsou v L p husté . Alternativně lze použít skutečnost, že h ( t ) je distribuční derivát funkce log | t |/ π ; vtipkovat

Pro většinu operačních účelů lze Hilbertovu transformaci považovat za konvoluci. Například ve formálním smyslu je Hilbertova transformace konvoluce konvolucí Hilbertovy transformace aplikované pouze na jeden z těchto faktorů:

To platí přesně, pokud u a v jsou kompaktně podporované distribuce, protože v takovém případě

Při přechodu na vhodnou mez je tedy také pravdivé, pokud uL p a vL q za předpokladu, že

z věty kvůli Titchmarshovi.

Neměnnost

Hilbertova transformace má následující vlastnosti invariance .

  • Dojíždí s překlady. To znamená, že dojíždí s operátory T a f ( x ) = f ( x + a ) pro všechny a in
  • Dojíždí pozitivními dilatacemi. To znamená, že dojíždí s operátory M λ f ( x ) = f ( λ x ) pro všechny λ > 0 .
  • To anticommutes s odrazem R f ( x ) = f (- x ) .

Až do multiplikativní konstanty je Hilbertova transformace jediným omezeným operátorem na L 2 s těmito vlastnostmi.

Ve skutečnosti existuje širší soubor operátorů, kteří dojíždějí s Hilbertovou transformací. Skupina působí pomocí unitárních operátorů U g na prostor podle vzorce

Tato unitární reprezentace je příkladem hlavní sériovou zastoupení všech v tomto případě je redukovat, rozdělení jako ortogonální součet dvou invariantních podprostorů, Hardy prostoru a jeho konjugátu. Jedná se o prostory hraničních hodnot L 2 holomorfních funkcí na horní a dolní polorovině. a jeho konjugát sestává přesně z těch funkcí L 2 s Fourierovými transformacemi mizícími na záporných a kladných částech skutečné osy. Vzhledem k tomu, Hilbertova transformace se rovná H = - i (2 P - I) , přičemž P je kolmý průmět z na a I operátor identity , to znamená, že i jeho kolmé jsou eigenspaces z H pro vlastní čísla ± i . Jinými slovy, H dojíždí s operátory U g . Omezení operátorů U g to a jeho konjugátu poskytují neredukovatelné reprezentace -takzvaného limitu reprezentací diskrétních sérií .

Rozšíření definiční domény

Hilbertova transformace distribucí

Dále je možné rozšířit Hilbertovu transformaci do určitých distribučních prostorů ( Pandey 1996 , kapitola 3). Vzhledem k tomu, Hilbertova transformace dojíždí s diferenciací, a je ohraničený operátor na L p , H omezuje, aby kontinuální transformace na inverzní limitu z Soboleových prostorů :

Hilbertovu transformaci pak lze definovat na duálním prostoru , označeném , sestávajícím z distribucí L p . Toho je dosaženo pomocí duality párování: U , definujeme:

Je možné definovat Hilbertovu transformaci na prostoru temperovaných distribucí také přístupem díky Gel'fandovi a Shilovovi, ale vzhledem k singularitě v integrálu je zapotřebí podstatně větší péče.

Hilbertova transformace ohraničených funkcí

Hilbertovu transformaci lze definovat i pro funkce , ale vyžaduje určité úpravy a upozornění. Správně pochopen, Hilbertovy transformace mapy do Banachova prostoru z ohraničené střední oscilace (BMO) tříd.

Hilbertova transformace ohraničené funkce, vykládaná naivně, je jasně špatně definována. Například s u = sgn ( x ) se integrál definující H ( u ) rozchází téměř všude na ± ∞ . Ke zmírnění těchto obtíží je proto Hilbertova transformace funkce L definována následující regularizovanou formou integrálu

kde jak výše h ( x ) =1/πx a

Upravená transformace H souhlasí s původní transformací funkcí kompaktní podpory z obecného výsledku Calderóna a Zygmunda. Výsledný integrál navíc bodově konverguje téměř všude a s ohledem na normu BMO na funkci omezené střední oscilace.

Hluboký výsledek z Fefferman práce je, že funkce je ohraničené střední kmitání pouze v případě, že má formu f + H ( g ) pro některé .

Funkce konjugátu

Hilbertovu transformaci lze chápat jako dvojici funkcí f ( x ) a g ( x ) tak, že funkce

je hraniční hodnota holomorfní funkce F ( z ) v horní polorovině. Za těchto okolností, pokud jsou f a g dostatečně integrovatelné, pak jeden je Hilbertova transformace druhého.

Předpokládejme, že pak, podle teorie Poissonova integrálu , f připouští jedinečný harmonické rozšíření do horní polorovině, a toto rozšíření je dána

což je konvoluce f s Poissonovým jádrem

Kromě toho existuje v horní polorovině jedinečná harmonická funkce v definovaná tak, že F ( z ) = u ( z ) + iv ( z ) je holomorfní a

Tato harmonická funkce je získána z f konvolucí s konjugovaným Poissonovým jádrem

Tím pádem

Skutečné a imaginární části jádra Cauchy skutečně jsou

takže F = u + iv je holomorfní podle Cauchyho integrálního vzorce .

Funkce v získaný z u tímto způsobem se nazývá harmonické konjugát z u . (Netangenciální) hraniční limit v ( x , y ) jako y → 0 je Hilbertova transformace f . Stručně tedy

Titchmarshova věta

Titchmarshova věta (pojmenovaná po EC Titchmarshovi, který ji zahrnul do své práce z roku 1937) upřesňuje vztah mezi hraničními hodnotami holomorfních funkcí v horní polorovině a Hilbertovou transformací. Poskytuje nezbytné a dostatečné podmínky, aby komplexně oceněná funkce integrovatelná do čtverce F ( x ) na reálné přímce byla hraniční hodnotou funkce v Hardyho prostoru H 2 ( U ) holomorfních funkcí v horní polorovině U .

Věta uvádí, že následující podmínky pro komplexně integrovanou funkci čtvercově integrovatelné jsou ekvivalentní:

  • F ( x ) je limita zx holomorfní funkce F ( z ) v horní polorovině tak, že
  • Skutečné a imaginární části F ( x ) jsou Hilbertovy transformace navzájem.
  • Fourierova transformace zmizí pro x <0 .

Slabší výsledek platí pro funkce třídy L p pro p > 1 . Konkrétně, pokud F ( z ) je holomorfní funkce taková, že

pro všechna y pak existuje komplexně hodnocená funkce F ( x ) v tom smyslu, že F ( x + iy ) → F ( x ) v L p normě jako y → 0 (stejně jako bodové držení téměř všude ). Kromě toho,

kde f je funkce se skutečnou hodnotou in a g je Hilbertova transformace (třídy L p ) f .

To neplatí v případě p = 1 . Ve skutečnosti, Hilbertova transformace z L 1, funkce f potřeby nekonverguje v střední k jinému L 1 funkci. Nicméně Hilbertova transformace f konverguje téměř všude na konečnou funkci g tak, že

Tento výsledek je přímo analogický výsledku Andreyho Kolmogorova pro Hardyho funkce na disku. Ačkoli se obvykle nazývá Titchmarshova věta, výsledek agreguje mnoho práce ostatních, včetně Hardyho, Paleyho a Wienera (viz Paley -Wienerova věta ), stejně jako práce Riesze, Hille a Tamarkina

Problém Riemann – Hilbert

Jedna forma problému Riemann – Hilbert se snaží identifikovat páry funkcí F + a F - takové, že F + je holomorfní v horní polorovině a F - je holomorfní v dolní polorovině, takže pro x podél skutečné osa,

kde f ( x ) je nějaká daná funkce s reálnou hodnotou . Levou stranu této rovnice lze chápat buď jako rozdíl mezí F ± od příslušných polorovin, nebo jako distribuci hyperfunkce . Dvě funkce této formy jsou řešením problému Riemann – Hilbert.

Formálně, pokud F ± vyřeší problém Riemann – Hilbert

potom je Hilbertova transformace f ( x ) dána vztahem

Hilbertova transformace na kruhu

Pro periodickou funkci f je definována kruhová Hilbertova transformace:

Kruhová Hilbertova transformace se používá při charakterizaci Hardyho prostoru a při studiu konjugované funkce ve Fourierově řadě. Jádro,

je známý jako Hilbertovo jádro, protože v této podobě byla původně studována Hilbertova transformace.

Hilbertovo jádro (pro kruhovou Hilbertovu transformaci) lze získat tak, že Cauchyho jádro bude 1 / x periodické. Přesněji pro x ≠ 0

Mnoho výsledků o kruhové Hilbertově transformaci lze odvodit z odpovídajících výsledků pro Hilbertovu transformaci z této korespondence.

Další přímější spojení poskytuje Cayleyova transformace C ( x ) = ( x - i ) / ( x + i ) , která nese skutečnou přímku na kruh a horní polovinu roviny na disk jednotky. Vyvolává unitární mapu

z L 2 ( T ) na Provozovatel U nese Hardyho prostor H 2 ( T ) na Hardy prostoru .

Hilbertova transformace při zpracování signálu

Bedrosianova věta

Bedrosianova věta uvádí, že Hilbertova transformace součinu nízkoprůchodového a vysokoprůchodového signálu s nepřekrývajícími se spektry je dána součinem nízkoprůchodového signálu a Hilbertovy transformace horního průchodu, popř.

kde f LP a f HP jsou signály dolního a vysokého průchodu. Kategorie komunikačních signálů, na které se to vztahuje, se nazývá model úzkopásmového signálu. Členem této kategorie je amplitudová modulace vysokofrekvenčního sinusového „nosiče“:

kde u m ( t ) je tvar vlny „zprávy“ s úzkou šířkou pásma, například hlas nebo hudba. Potom podle Bedrosianovy věty:

Analytická reprezentace

Specifický typ konjugované funkce je :

známý jako analytické reprezentace z The jméno ukazuje jeho matematický povolnost, převážně v důsledku Eulerovy rovnice . Při použití Bedrosianovy věty na úzkopásmový model je analytická reprezentace :

 

 

 

 

( Rovnice 1 )

Vlastnost Fourierovy transformace naznačuje, že tato složitá heterodynní operace může posunout všechny negativní frekvenční složky u m ( t ) nad 0 Hz. V takovém případě je imaginární částí výsledku Hilbertova transformace skutečné části. Toto je nepřímý způsob produkce Hilbertových transformací.

Úhlová (fázová/frekvenční) modulace

Formulář:

se nazývá úhlová modulace , která zahrnuje fázovou i frekvenční modulaci . Okamžitý kmitočet je     pro dostatečně velké W , ve srovnání s :

a:

Modulace jednoho postranního pásma (SSB)

Když u m ( t )Eq.1 je také analytická reprezentace (tvaru vlny zprávy), to je:

výsledkem je jednostranná modulace pásma :

jehož přenášenou složkou je:

Kauzalita

Tato funkce představuje dvě výzvy praktické implementace jako konvoluci:

  • Jeho trvání je nekonečné (technicky nekonečná podpora ). Místo toho je nutné použít aproximaci konečné délky. Ale okenních délku také snižuje účinný dosah frekvenční transformace. Čím kratší je okno, tím větší jsou ztráty na nízkých a vysokých frekvencích. Viz také kvadraturní filtr .
  • Je to nekauzální filtr . Je tedy vyžadována opožděná verze . Odpovídající výstup je následně zpožděn Při vytváření imaginární části analytického signálu musí být zdroj (skutečná část) zpožděn o ekvivalentní množství.

Diskrétní Hilbertova transformace

Obrázek 1 : Filtr, jehož frekvenční odezva je omezena pásmem na přibližně 95% Nyquistovy frekvence
Obrázek 2 : Hilbertův transformační filtr s vysokofrekvenční frekvenční charakteristikou
Obrázek 3 .
Obrázek 4 . Hilbertova transformace cos ( ωt ) je sin ( ωt ) . Tento obrázek ukazuje sin (ωt) a dvě přibližné Hilbertovy transformace vypočítané funkcí knihovny MATLAB, hilbert ()
Obrázek 5 . Diskrétní Hilbertovy transformace kosinové funkce pomocí konvoluce po částech

Pro diskrétní funkce , s diskrétního Fourierova transformace (DTFT) , a diskrétní Hilbertova transformace , DTFT ze v oblasti - π <ω < π je dána vztahem:

Inverzní DTFT pomocí konvoluční věty je:

kde

což je nekonečná impulsní reakce (IIR). Pokud je konvoluce prováděna numericky, aproximace FIR je nahrazena h [ n ] , jak je znázorněno na obrázku 1 . FIR filtr s lichým počtem asymetrických koeficientů se nazývá typ III, který ve své podstatě vykazuje odezvy nulové velikosti na frekvencích 0 a Nyquist, což má v tomto případě za následek tvar pásmového filtru. Konstrukce typu IV (sudý počet asymetrických koeficientů) je znázorněna na obrázku 2 . Protože odezva velikosti na Nyquistově frekvenci nevypadává, přibližuje se ideálnímu Hilbertovu transformátoru o něco lépe než lichému kohoutkovému filtru. nicméně

  • Typická (tj. Správně filtrovaná a vzorkovaná) u [ n ] sekvence nemá žádné užitečné složky na Nyquistově frekvenci.
  • Impulsní odezva typu IV vyžaduje 1 / 2 posun vzorku v sekvenci h [ n ] . To způsobí, že se koeficienty s nulovou hodnotou stanou nenulovými, jak je vidět na obrázku 2 . Konstrukce typu III je tedy potenciálně dvakrát účinnější než typ IV.
  • Skupinové zpoždění návrhu typu III je celočíselný počet vzorků, což usnadňuje zarovnání s vytvořením analytického signálu . Skupinové zpoždění typu IV je na půli cesty mezi dvěma vzorky.

Funkce MATLAB , hilbert (u, N) , spojuje au [n] posloupnost s periodickým součtem :

   

a vrací jeden cyklus ( N vzorků) periodického výsledku v imaginární části komplexní hodnotové výstupní sekvence. Konvoluce je implementována ve frekvenční oblasti jako součin pole     se vzorky distribuce - i sgn ( ω ) (jejíž skutečné a imaginární složky jsou všechny jen 0 nebo  ± 1 ). Obrázek 3 porovnává poloviční cyklus h N [ n ] s ekvivalentní délkovou částí h [ n ] . Vzhledem k tomu, FIR aproximace pro označován náhradou za - i sgn ( omega ) vzorky vede FIR verzi konvoluce.

Reálná část výstupního sekvence je původní vstupní sekvence tak, že komplex výstup je analytický znázornění z u [ n ] . Když je vstupem segment čistého kosinu, výsledná konvoluce pro dvě různé hodnoty N je znázorněna na obrázku 4 (červené a modré grafy). Efekty na hraně zabraňují tomu, aby byl výsledek čistě sinusová funkce (zelený plot). Protože h N [ n ] není sekvence FIR, teoretický rozsah účinků je celá výstupní sekvence. Ale rozdíly od sinusové funkce se zmenšují se vzdáleností od okrajů. Parametr N je délka výstupní sekvence. Pokud překročí délku vstupní sekvence, vstup se upraví připojením prvků s nulovou hodnotou. Ve většině případů to snižuje velikost rozdílů. Jejich trvání však dominují inherentní doby vzestupu a poklesu h [ n ] impulzní odezvy.

Ocenění okrajových efektů je důležité, když se k provedení konvoluce na dlouhé u [ n ] sekvenci používá metoda zvaná overlap-save . Segmenty délky N jsou spojeny s periodickou funkcí:

Když doba trvání nenulových hodnot je je výstupní sekvence zahrnuje N - M + 1 vzorky M -1 výstupů jsou vyřazeny z každého bloku N a vstupní bloky se překrývají o tuto částku, aby se zabránilo mezerám.

Obrázek 5 je příkladem použití funkce IIR hilbert (·) a aproximace FIR. V tomto příkladu je sinusová funkce vytvořena výpočtem diskrétní Hilbertovy transformace kosinové funkce, která byla zpracována ve čtyřech překrývajících se segmentech a spojena dohromady. Jak ukazuje výsledek FIR (modrý), zkreslení patrná ve výsledku IIR (červená) nejsou způsobena rozdílem mezi h [ n ] a h N [ n ] (zelená a červená na obrázku 3 ). Skutečnost, že h N [ n ] je zúžené ( okénkové ), je v této souvislosti skutečně užitečná. Skutečným problémem je, že není dostatečně oknem. Účinně M = N , zatímco metoda překrývání-uložení potřebuje M < N .

Číslově teoretická Hilbertova transformace

Číselná teoretická Hilbertova transformace je rozšířením diskrétní Hilbertovy transformace na celá čísla modulující příslušné prvočíslo. V tomto vyplývá zobecnění diskrétní Fourierovy transformace na číselné teoretické transformace. Hilbertovu transformaci teoretických čísel lze použít ke generování sad ortogonálních diskrétních sekvencí.

Viz také

Poznámky

Citace stránek

Reference

Další čtení

externí odkazy