Funkce Dirac delta - Dirac delta function

Schematické znázornění Diracovy delta funkce čarou převyšovanou šipkou. Výška šipky je obvykle určena k určení hodnoty jakékoli multiplikativní konstanty, která poskytne oblast pod funkcí. Druhou konvencí je napsat oblast vedle šipky.
Funkce Diracovy delty jako limitu jako (ve smyslu distribucí ) sekvence nulových středů normálních distribucí

V matematiky je funkce Dirac delta ( δ funkce ), také známý jako impulsní symbolu, je generalizovaná funkcí nebo distribuci nad reálnými čísly , jejichž hodnota je nula všude kromě na nulu, a jehož integrální po celé reálné ose je roven do jednoho. To může také být interpretován jako lineární funkční , který mapuje každou funkci na hodnotu nula, nebo jako slabé hranice části sekvence z bump funkcí , které jsou k nule přes většinu z reálné ose, s vysokou hrotem na původu. Bump funkcím se proto někdy říká „přibližné“ nebo „rodící se“ delta funkce.

Funkci delta představil fyzik Paul Dirac jako nástroj pro normalizaci stavových vektorů. Má také využití v teorii pravděpodobnosti a zpracování signálu . Protože to není skutečná matematická funkce , někteří matematici proti tomu protestovali jako nesmysl, dokud Laurent Schwartz nerozvinul teorii rozdělení.

Funkce Kroneckerova delta , která je obvykle definována v diskrétní doméně a nabývá hodnot 0 a 1, je diskrétním analogem Diracovy delta funkce.

Motivace a přehled

Graf funkce delta je obvykle myšlenka jako po celou x aretačním kroužkem a pozitivní y aretačním kroužkem. Diracova delta se používá k modelování funkce vysokého úzkého hrotu ( impuls ) a dalších podobných abstrakcí , jako je bodový náboj , hmota bodu nebo elektronový bod. Například k výpočtu dynamiky příslušníky kulečníková koule úderu, lze aproximovat sílu dopadu pomocí delta funkce. Přitom člověk nejen zjednodušuje rovnice, ale je také schopen vypočítat pohyb koule pouze uvažováním celkového impulsu srážky bez podrobného modelu veškerého pružného přenosu energie na subatomárních úrovních (například) .

Abychom byli konkrétní, předpokládejme, že kulečníková koule je v klidu. Časem je zasažen jiným míčem, který mu dodá hybnost P , in . Výměna hybnosti není ve skutečnosti okamžitá, je zprostředkována elastickými procesy na molekulární a subatomární úrovni, ale pro praktické účely je vhodné považovat tento přenos energie za efektivně okamžitý. Síla je tedy . (Jednotky jsou .)

Chcete -li tuto situaci modelovat přísněji, předpokládejme, že síla je místo toho rovnoměrně rozložena v malém časovém intervalu . To znamená,

Integrace pak zjistí hybnost kdykoli t :

Modelová situace okamžitého přenosu hybnosti nyní vyžaduje vzít limit jako dávání

Zde jsou funkce považovány za užitečné aproximace myšlenky okamžitého přenosu hybnosti.

Funkce delta nám umožňuje sestrojit idealizovaný limit těchto aproximací. Bohužel skutečný limit funkcí (ve smyslu bodové konvergence ) je všude nula, kromě jediného bodu, kde je nekonečný. Abychom správně rozuměli funkci delta, měli bychom místo toho trvat na této vlastnosti

který platí pro všechny , by se měl i nadále držet v limitu. V rovnici se tedy rozumí, že mez je vždy vzata mimo integrál .

V aplikované matematice, jak jsme to udělali zde, je delta funkce často manipulována jako druh limitu ( slabý limit ) posloupnosti funkcí, z nichž každý člen má na počátku vysoký bod: například posloupnost Gaussovské distribuce se soustředily na počátek s rozptylem směřujícím k nule.

Navzdory svému názvu není funkce delta skutečně funkcí, alespoň ne obvyklou s doménou a rozsahem v reálných číslech . Například objekty f ( x ) = δ ( x ) a g ( x ) = 0 jsou všude stejné, s výjimkou bodu x = 0, ale mají integrály, které jsou různé. Podle Lebesgueovy integrační teorie platí , že pokud f a g jsou funkce takové, že f = g téměř všude , pak je f integrovatelné právě tehdy, když je g integrovatelné a integrály f a g jsou totožné. Přísný přístup k tomu, aby Diracova delta funkce byla sama o sobě matematickým objektem, vyžaduje teorii míry nebo teorii distribucí .

Dějiny

Joseph Fourier představil to, čemu se nyní říká Fourierova integrální věta, ve svém pojednání Théorie analytique de la chaleur ve formě:

což se rovná zavedení funkce δ ve tvaru:

Později Augustin Cauchy vyjádřil větu pomocí exponenciálů:

Cauchy poukázal na to, že za určitých okolností je pořadí integrace v tomto výsledku významné (kontrast Fubiniho věty ).

Jak je odůvodněno pomocí teorie distribucí , Cauchyovu rovnici lze přeskupit tak, aby se podobala Fourierově původní formulaci, a vystavit δ -funkci jako

kde funkce δ je vyjádřena jako

Přísná interpretace exponenciální formy a různá omezení funkce f nezbytné pro její aplikaci trvala několik století. Problémy s klasickou interpretací jsou vysvětleny následovně:

Největší nevýhodou klasické Fourierovy transformace je poměrně úzká třída funkcí (originálů), pro které ji lze efektivně vypočítat. Totiž je nutné, aby tyto funkce dostatečně rychle klesaly na nulu (v sousedství nekonečna), aby byla zajištěna existence Fourierova integrálu. Například Fourierova transformace takových jednoduchých funkcí, jako jsou polynomy, v klasickém smyslu neexistuje. Rozšíření klasické Fourierovy transformace na distribuce značně rozšířilo třídu funkcí, které bylo možné transformovat, a to odstranilo mnoho překážek.

Další vývoj zahrnoval zobecnění Fourierova integrálu, „počínaje Plancherelovou průlomovou teorií L 2 (1910), pokračující Wienerovou a Bochnerovou tvorbou (kolem roku 1930) a kulminující sloučením do teorie distribucí L. Schwartze (1945) ... “, a vedoucí k formálnímu rozvoji funkce Diracovy delty.

V textu Augustina Louise Cauchyho z roku 1827 se výslovně objevuje nekonečně malý vzorec pro nekonečně vysokou funkci jednotkové impulzní delty (nekonečně malá verze Cauchyho distribuce ) . Siméon Denis Poisson zvažoval problém v souvislosti se studiem šíření vln stejně jako Gustav Kirchhoff o něco později. Kirchhoff a Hermann von Helmholtz také představili jednotkový impuls jako limit Gaussianů , což také odpovídalo pojmu lorda Kelvina o bodovém zdroji tepla. Na konci 19. století používal Oliver Heaviside formální Fourierovu sérii k manipulaci s jednotkovým impulzem. Funkce Diracovy delty jako takové byla zavedena jako „pohodlný zápis“ Paulem Diracem v jeho vlivné knize z roku 1930 Zásady kvantové mechaniky . Říkal tomu „delta funkce“, protože ho používal jako spojitý analog diskrétní Kroneckerovy delty .

Definice

O Diracově deltě lze volně uvažovat jako o funkci na reálné přímce, která je všude nulová, kromě jejího počátku, kde je nekonečná,

a který je také omezen na uspokojení identity

Toto je pouze heuristická charakteristika. Diracova delta není funkcí v tradičním smyslu, protože žádná funkce definovaná na skutečných číslech tyto vlastnosti nemá. Funkce Dirac delta může být přísně definována buď jako distribuce, nebo jako míra .

Jako měřítko

Jedním ze způsobů, jak důsledně zachytit pojem Diracovy delta funkce, je definovat míru , nazývanou Diracova míra , která přijímá jako argument podmnožinu A skutečné linie R a vrací δ ( A ) = 1, pokud 0 ∈ A , a δ ( A ) = 0 jinak. Pokud je funkce delta koncipována jako modelování idealizované hmotný bod při 0 ° C, pak δ ( ) představuje hmotnost, které je součástí A . Lze pak definovat integrál proti δ jako integrál funkce proti tomuto rozdělení hmotnosti. Formálně Lebesgueův integrál poskytuje potřebné analytické zařízení. Lebesgueův integrál s ohledem na míru δ splňuje

pro všechny souvisle kompaktně podporované funkce f . Míra δ není absolutně spojitá s ohledem na Lebesgueovu míru - ve skutečnosti je to singulární míra . V důsledku toho delta opatření nemá derivát Radon -Nikodym (s ohledem na Lebesgueovu míru) - žádnou skutečnou funkci, pro kterou by vlastnost

drží. Výsledkem je, že tato notace je praktickým zneužitím notace a není standardním ( Riemann nebo Lebesgue ) integrálem.

Jako měřítko pravděpodobnosti na R je delta opatření charakterizováno kumulativní distribuční funkcí , což je funkce jednotkového kroku .

To znamená, že H ( x ) je integrálem funkce kumulativního indikátoru 1 (−∞, x ] vzhledem k míře δ ;

posledně jmenovaný je měřítkem tohoto intervalu; formálněji δ ((−∞, x ]) . Integraci delta funkce proti spojité funkci lze tedy správně chápat jako Riemann – Stieltjesův integrál :

Všechny vyšší momenty z delta jsou nulové. Zejména charakteristická funkce a funkce generující moment jsou obě stejné jako jedna.

Jako distribuce

V teorii distribucí není generalizovaná funkce považována za funkci sama o sobě, ale pouze za to, jak ovlivňuje jiné funkce, když je „integrována“ proti nim. V souladu s touto filozofií stačí ke správné definici funkce delta říci, co je „integrál“ funkce delta proti dostatečně „dobré“ testovací funkci  φ . Testovací funkce jsou také známé jako bump funkce . Pokud je funkce delta již chápána jako míra, pak Lebesgueův integrál testovací funkce proti této míře dodává potřebný integrál.

Typický prostor testovacích funkcí se skládá ze všech hladkých funkcí na R s kompaktní podporou, které mají tolik derivací, kolik je potřeba. Distribuce Dirac delta je lineární funkce v prostoru testovacích funkcí a je definována

 

 

 

 

( 1 )

pro každou testovací funkci .

Aby δ byla správně distribucí, musí být spojitá ve vhodné topologii na prostoru testovacích funkcí. Obecně platí, že pro lineární funkční S na prostoru testovacích funkcí k definování rozdělení je nutné a dostatečné, že pro každé kladné celé číslo N existuje celé číslo M N a konstanta C N takové, že pro každou testovací funkci φ , jeden má nerovnost

S delta distribuci, jeden má takovou nerovnost (s C N = 1) s M N = 0 pro všechny N . Tak δ je distribuce řádu nula. Je to navíc distribuce s kompaktní podporou ( podpora je {0}).

Delta distribuci lze také definovat několika ekvivalentními způsoby. Například se jedná o distribuční derivát z kroku funkce Heaviside . To znamená, že pro každou testovací funkci φ má jeden

Intuitivně, pokud by byla povolena integrace po částech , měl by se tento integrál zjednodušit na

a skutečně je pro Stieltjesův integrál povolena určitá forma integrace po částech, a v takovém případě člověk má

V kontextu teorie míry vede Diracovo opatření k distribuci integrací. Naopak rovnice ( 1 ) definuje Daniellův integrál v prostoru všech kompaktně podporovaných spojitých funkcí φ, které lze pomocí Rieszovy věty o reprezentaci reprezentovat jako Lebesgueův integrál φ týkající se nějaké radonové míry .

Obecně platí, že když je použit termín „ Diracova delta funkce “, znamená to spíše distribuce než míry, Diracova míra je jedním z několika termínů pro odpovídající pojem v teorii opatření. Některé zdroje mohou také používat termín Dirac delta distribution .

Zobecnění

Delta funkce může být definována v n -rozměrném euklidovském prostoru R n jako míra taková, že

pro každou kompaktně podporovanou spojitou funkci f . Jako měřítko je n -dimenzionální delta funkce součinem 1 -rozměrných delta funkcí v každé proměnné zvlášť. Formálně tedy s x = ( x 1 , x 2 , ..., x n ) má jeden

 

 

 

 

( 2 )

Delta funkce může být také definována ve smyslu rozdělení přesně jako výše v jednorozměrném případě. I přes široké použití v technických kontextech by ( 2 ) mělo být zacházeno opatrně, protože produkt distribucí lze definovat pouze za docela úzkých okolností.

Pojem Diracovy míry má smysl na každé sadě. Takže pokud X je množina, x 0X je označený bod a Σ je jakýkoliv sigma algebry podmnožin X , pak se míra definována na sady ∈ å podle

je míra delta nebo jednotková hmotnost koncentrovaná na x 0 .

Další běžnou generalizací funkce delta je diferencovatelný variátor, kde lze také díky distribuovatelné struktuře využít většinu jejích distribučních vlastností . Funkce delta na potrubí M se středem v bodě x 0M je definována jako následující rozdělení:

 

 

 

 

( 3 )

Pro všechny podporované kompaktně hladké reálných funkcí cp na M . Běžným zvláštním případem této konstrukce je ten, ve kterém M je otevřená množina v euklidovském prostoru R n .

Na místně kompaktním Hausdorffově prostoru X je Diracova delta míra soustředěná v bodě x je radonová míra spojená s Daniell integrálem ( 3 ) na kompaktně podporovaných spojitých funkcích φ . Na této úrovni obecnosti již není kalkul jako takový možný, nicméně jsou k dispozici různé techniky z abstraktní analýzy. Například mapování je kontinuální vkládání X do prostoru konečných radonových měr na X , vybavených jeho vágní topologií . Kromě toho je konvexní obal z obrazu X podle tohoto zakotvení je hustá v prostoru pravděpodobnostních opatření na X .

Vlastnosti

Škálování a symetrie

Funkce delta splňuje následující vlastnost škálování pro nenulovou skalární α:

a tak

 

 

 

 

( 4 )

Důkaz:

Zejména funkce delta je rovnoměrné rozdělení v tom smyslu

který je homogenní stupně −1.

Algebraické vlastnosti

Distribuční produkt z delta s x se rovná nule:

A naopak, pokud xf ( x ) = xg ( x ) , kde f a g jsou rozdělení, pak

pro nějakou konstantu c .

Překlad

Integrál časově zpožděné Diracovy delty je

Toto je někdy označováno jako vlastnost prosévání nebo vlastnost vzorkování . Funkce delta se říká, že „tříbil ven“ hodnotě na t = T .

Z toho vyplývá, že účinek konvoluce funkce f ( t ) s časově zpožděnou Diracovou deltou je na časové zpoždění f ( t ) o stejnou hodnotu:

To platí za přesné podmínky, že f je temperované rozdělení (viz diskuse o Fourierově transformaci níže ). Jako zvláštní případ například máme identitu (chápanou ve smyslu distribuce)

Kompozice s funkcí

Obecněji může být rozdělení delta složeno s hladkou funkcí g ( x ) takovým způsobem, že platí známý vzorec změny proměnných, že

za předpokladu, že g je spojitě diferencovatelná funkce s g 'nikde nula. To znamená, že existuje jedinečný způsob, jak distribuci přiřadit význam , aby tato identita platila pro všechny kompaktně podporované testovací funkce f . Proto musí být doména rozdělena, aby se vyloučil bod g ′ = 0. Toto rozdělení splňuje δ ( g ( x )) = 0, pokud g není nikde nula, a jinak pokud g má skutečný kořen v x 0 , pak

Je tedy přirozené definovat kompozici δ ( g ( x )) pro spojitě diferencovatelné funkce g pomocí

kde součet sahá přes všechny kořeny (tj. všechny různé) g ( x ), o kterých se předpokládá, že jsou jednoduché . Tak například

V integrální formě lze generalizovanou vlastnost škálování zapsat jako

Vlastnosti v n rozměrech

Distribuce delta v n -dimenzionálním prostoru splňuje místo toho následující vlastnost škálování,

takže δ je homogenní rozdělení stupně - n .

Při jakémkoli odrazu nebo rotaci ρ je funkce delta invariantní,

Stejně jako v případě jedné proměnné je možné definovat složení δ pomocí bi-Lipschitzovy funkce g : R nR n jedinečně tak, aby identita

pro všechny kompaktně podporované funkce f .

Pomocí vzorce coarea z teorie geometrických opatření lze také definovat složení delta funkce s ponořením z jednoho euklidovského prostoru do jiného s jinou dimenzí; výsledkem je typ proudu . Ve zvláštním případě kontinuálně diferencovatelné funkce g : R nR tak, aby se sklon o g nikde nula, následující identita drží

kde integrál vpravo je nad g −1 (0), ( n -1) -dimenzionální povrch definovaný g ( x ) = 0 s ohledem na míru obsahu Minkowski . Toto je známé jako integrál jednoduché vrstvy .

Obecněji řečeno, je-li S je hladký nadplochy z R n , pak můžeme spojit k S distribuce, který integruje jakýkoliv kompaktně podporováno hladké funkce g nad S :

kde σ je nadplochy opatření spojená s S . Toto zobecnění je spojena s teorie potenciálu z jednoduchých vrstev potenciálů na S . Pokud D je doména v R n s hladkým ohraničující S , pak δ S se rovná normální derivátu z funkce indikátoru části D v distribučním smyslu,

kde n je vnější normál. Důkaz najdete např. V článku o funkci povrchové delty .

Fourierova transformace

Funkce delta je temperovaná distribuce , a proto má dobře definovanou Fourierovu transformaci . Formálně člověk najde

Správně řečeno, Fourierova transformace distribuce je definována uložením sebe-sousedství Fourierovy transformace v rámci dualistického párování temperovaných distribucí se Schwartzovými funkcemi . Je tedy definováno jako uspokojivé jedinečné temperované rozdělení

pro všechny funkce Schwartz . A skutečně z toho vyplývá, že

V důsledku této identity je konvoluce funkce delta s jakoukoli jinou temperovanou distribucí S jednoduše S :

To znamená, že δ je prvek identity pro konvoluci na temperovaných distribucích a ve skutečnosti je prostor kompaktně podporovaných distribucí pod konvolucí asociativní algebrou s funkcí delta funkce. Tato vlastnost je zásadní při zpracování signálu , protože konvoluce s temperovaným rozdělením je lineární časově invariantní systém a aplikace lineárního časově invariantního systému měří jeho impulsní odezvu . Impulsní odezvu lze vypočítat na libovolný požadovaný stupeň přesnosti výběrem vhodné aproximace pro δ , a jakmile je známa, charakterizuje systém úplně. Viz teorie systému LTI § Impulsní odezva a konvoluce .

Inverzní Fourierova transformace temperovaného rozdělení f ( ξ ) = 1 je delta funkce. Formálně je to vyjádřeno

a přísněji z toho plyne od

pro všechny funkce Schwartz f .

V těchto termínech, funkce delta poskytuje sugestivní prohlášení o majetku orthogonality jádra Fourier na výzkum . Formálně jeden má

To je samozřejmě zkratka pro tvrzení, že Fourierova transformace temperované distribuce

je

což opět následuje uložením sebe-sousedství Fourierovy transformace.

Podle analytické pokračování Fourierovy transformace je Laplaceova transformace funkce delta je zjištěno, že

Distribuční deriváty

Distribuční derivát Diracovy delta distribuce je distribuce δ ′ definovaná na kompaktně podporovaných hladkých testovacích funkcích φ pomocí

První rovnost je druh integrace po částech, protože kdyby δ tehdy byla skutečná funkce

K -The derivát delta je definován stejně jako distribuce uvedené na testovacích funkcemi

Zejména δ je nekonečně odlišitelná distribuce.

První derivací funkce delta je distribuční limit rozdílových kvocientů:

Přesněji řečeno, jeden má

kde τ h je operátor překlad, definovaná na funkcemi ▼ se h cp ( x ) = cp ( x + h ) , a na distribuční S pomocí

V teorii elektromagnetismu první derivace delta funkce představuje bodový magnetický dipól umístěný na počátku. V souladu s tím se označuje jako dipólová nebo dubletová funkce .

Derivace delta funkce splňuje řadu základních vlastností, včetně:

Poslední z těchto vlastností lze snadno prokázat aplikací distribuční derivační definice, Liebnitzovy věty a linearity vnitřního produktu:

Navíc konvoluce δ ′ s kompaktně podporovanou hladkou funkcí f je

což vyplývá z vlastností distribučního derivátu konvoluce.

Vyšší rozměry

Obecněji platí, že na otevřené množině U v n -dimenzionálním euklidovském prostoru R n je Diracova delta distribuce soustředěná v bodě aU definována

Pro všechny cpS ( U ) , prostor všech hladkých kompaktně podporované funkce na U . Pokud α = ( α 1 , ..., α n ) je jakýkoliv multi-index a ∂ α označuje spojené smíšené parciální derivace operátor, pak α th derivace ∂ α ó z delta a je dána vztahem

To znamená, že α th derivace δ a je rozdělení, jehož hodnota na jakékoli testovací funkci φ je α th derivace φ na a (s příslušným kladným nebo záporným znaménkem).

První dílčí derivace delta funkce jsou považovány za dvojité vrstvy podél souřadnicových rovin. Obecněji je normální derivát jednoduché vrstvy nesené na povrchu dvojitou vrstvou nesenou na tomto povrchu a představuje laminární magnetický monopól. Vyšší deriváty delta funkce jsou ve fyzice známé jako multipoly .

Vyšší deriváty vstupují do matematiky přirozeně jako stavební kameny kompletní struktury rozdělení s bodovou podporou. Jestliže S je jakékoli rozdělení na U podporované na množině { a } sestávající z jednoho bodu, pak existuje celé číslo m a koeficienty c α takové, že

Reprezentace funkce delta

Na funkci delta lze pohlížet jako na limit posloupnosti funkcí

kde η ε ( x ) se někdy říká rodící se delta funkce. Tato hranice je míněna ve slabém smyslu: buď to

 

 

 

 

( 5 )

pro všechny spojité funkce f s kompaktní podporou , nebo že tento limit platí pro všechny hladké funkce f s kompaktní podporou. Rozdíl mezi těmito dvěma mírně odlišnými způsoby slabé konvergence je často jemný: první je konvergence ve vágní topologii opatření a druhá konvergence ve smyslu distribucí .

Přibližování identity

Rodící se delta funkce η ε může být typicky konstruována následujícím způsobem. Nechť η je absolutně integrovatelná funkce na R celkového integrálu 1, a definuj

V n dimenzích se místo toho používá škálování

Potom jednoduchá změna proměnných ukazuje, že η ε má také integrál 1. Lze ukázat, že ( 5 ) platí pro všechny spojité kompaktně podporované funkce f , a tak η ε slabě konverguje k δ ve smyslu opatření.

Takto konstruované η ε jsou známé jako přiblížení identity . Tato terminologie je dána tím, že prostor L 1 ( R ) absolutně integrovatelných funkcí je uzavřen pod operací konvoluce funkcí: fgL 1 ( R ), kdykoli jsou f a g v L 1 ( R ). V L 1 ( R ) pro konvoluční produkt však neexistuje žádná identita : žádný prvek h takový, že fh = f pro všechna f . Nicméně sekvence η ε přibližuje takové identitě v tom smyslu, že

Tento limit platí ve smyslu střední konvergence (konvergence v L 1 ). K zajištění bodové konvergence téměř všude jsou zapotřebí další podmínky týkající se η ε , například to, že se jedná o mollifier spojený s kompaktně podporovanou funkcí .

Pokud je počáteční η = η 1 sama hladká a kompaktně podepřená, pak se posloupnost nazývá mollifier . Standardní mollifier se získá volbou r být vhodně normalizovanou funkci bump , například

V některých situacích, jako je numerická analýza , je žádoucí lineární aproximace identity po částech . Toho lze dosáhnout tím, že η 1 považujeme za funkci klobouku . S touto volbou η 1 má člověk

které jsou všechny spojité a kompaktně podepřené, i když ne hladké, a tudíž ani ne mollifier.

Pravděpodobnostní úvahy

V kontextu teorie pravděpodobnosti je přirozené uložit další podmínku, že počáteční η 1 v aproximaci identity by měla být kladná, protože taková funkce pak představuje rozdělení pravděpodobnosti . Konvoluce s distribucí pravděpodobnosti je někdy příznivá, protože nevede k překročení nebo podtržení, protože výstup je konvexní kombinací vstupních hodnot, a proto spadá mezi maximum a minimum vstupní funkce. Vezmeme -li η 1 jako jakékoli rozdělení pravděpodobnosti a ponecháme -li η ε ( x ) = η 1 ( x / ε ) / ε, jak je uvedeno výše, povede to k aproximaci identity. Obecně to konverguje rychleji na delta funkci, pokud má navíc η průměr 0 a má malé vyšší momenty. Například, je-li η 1 je rovnoměrné rozdělení na [-1/2, 1/2] , také známý jako obdélníkové funkce , pak:

Další příklad je s půlkruhovou distribucí Wigner

Toto je spojité a kompaktně podporované, ale není to mollifier, protože není hladké.

Poloskupiny

Vznikající delta funkce často vznikají jako konvoluční pologrupy . To se rovná dalšímu omezení, které musí splňovat konvoluce η ε s η δ

pro všechna ε , δ > 0 . Konvoluční pologrupy v L 1, které tvoří rodící se delta funkci, jsou vždy aproximací identity ve výše uvedeném smyslu, nicméně podmínka semigroup je docela silné omezení.

V praxi vznikají poloskupiny aproximující delta funkci jako základní řešení nebo Greenovy funkce pro fyzicky motivované eliptické nebo parabolické parciální diferenciální rovnice . V kontextu aplikované matematiky vznikají pologrupy jako výstup lineárního časově invariantního systému . Abstraktně, pokud A je lineární operátor působící na funkce x , pak konvoluční poloskupina vzniká řešením počátečního problému s hodnotou

ve kterém je limit jako obvykle chápán ve slabém smyslu. Nastavením η ε ( x ) = η ( ε , x ) získáte přidruženou rodící se funkci delta.

Některé příklady fyzicky důležitých semigroup konvoluce vyplývající z takového zásadního řešení zahrnují následující.

Tepelné jádro

Teplo jádro , definovaný

představuje teplotu v nekonečném drátu v čase t > 0, pokud je jednotka tepelné energie uložena na počátku drátu v čase t = 0. Tato poloskupina se vyvíjí podle jednorozměrné tepelné rovnice :

V teorii pravděpodobnosti , r | e ( x ) je normální rozdělení z rozptylu e a střední 0. představuje hustotu pravděpodobnosti v čase t = epsilon polohy částice začíná na počátku následující standardní Brownova pohybu . V této souvislosti je podmínka poloskupiny pak výrazem Markovovy vlastnosti Brownova pohybu.

Ve vyšších dimenzích euklidovského prostoru R n je tepelné jádro

a má stejnou fyzickou interpretaci, mutatis mutandis . Představuje také rodící se delta funkci ve smyslu, že η εδ ve smyslu distribuce jako ε → 0 .

Poissonovo jádro

Kernel Poisson

je základní řešení Laplaceovy rovnice v horní polorovině. Představuje elektrostatický potenciál v semi-nekonečné desce, jejíž potenciál podél okraje je pevně držen ve funkci delta. Poissonovo jádro je také úzce spjato s distribucí Cauchy a funkcí Epanechnikovova a Gaussova jádra . Tato poloskupina se vyvíjí podle rovnice

kde operátor je přísně definován jako Fourierův multiplikátor

Oscilační integrály

V oblastech fyziky, jako je šíření vln a mechanika vln , jsou zapojené rovnice hyperbolické, a proto mohou mít více singulárních řešení. V důsledku toho jsou rodící se delta funkce, které vznikají jako základní řešení souvisejících Cauchyho problémů, obecně oscilační integrály . Jako příklad, který pochází z roztoku Euler-Tricomi rovnice z transsonických dynamiky plynů , je přepočtou funkce vzdušný

Ačkoli používáme Fourierovu transformaci, je snadné vidět, že to v jistém smyslu generuje poloskupinu - není absolutně integrovatelná, a proto nemůže definovat poloskupinu ve výše uvedeném silném smyslu. Mnoho rodících se delta funkcí konstruovaných jako oscilační integrály konverguje pouze ve smyslu distribucí (příkladem je jádro Dirichlet níže), nikoli ve smyslu opatření.

Dalším příkladem je Cauchyův problém pro vlnovou rovnici v R 1+1 :

Řešení u představuje posunutí z rovnováhy nekonečného elastického řetězce s počáteční poruchou na počátku.

Mezi další aproximace identity tohoto druhu patří funkce sinc (široce používaná v elektronice a telekomunikacích)

a Besselova funkce

Rozklad v rovinné vlně

Jeden přístup ke studiu lineární parciální diferenciální rovnice

kde L je diferenciální operátor na R n , je nejprve hledat základní řešení, což je řešení rovnice

Když je L obzvláště jednoduché, lze tento problém často vyřešit přímo pomocí Fourierovy transformace (jako v případě již zmíněného Poissonova jádra a tepelného jádra). U složitějších operátorů je někdy jednodušší nejprve zvážit rovnici formuláře

kde h je funkce rovinné vlny , což znamená, že má formu

pro nějaký vektor ξ. Taková rovnice může být vyřešen (v případě, že koeficienty L jsou analytické funkce ) pomocí Cauchy-Kovalevskaya teorém , nebo (v případě, že koeficienty L jsou konstantní) od kvadratuře. Pokud tedy lze delta funkci rozložit na rovinné vlny, pak lze v zásadě řešit lineární parciální diferenciální rovnice.

Takový rozklad funkce delta na rovinné vlny byl součástí obecné techniky, kterou nejprve zavedl v zásadě Johann Radon a poté ji v této formě vyvinul Fritz John ( 1955 ). Vyberte k tak, aby n + k bylo sudé celé číslo, a pro reálné číslo s vložte

Potom δ získáme aplikací síly Laplaciána na integrál vzhledem k jednotkové sféře měřímeg ( x · ξ ) pro ξ v jednotkové sféře S n −1 :

Laplacian je zde interpretován jako slabá derivace, takže tato rovnice je brána tak, že pro jakoukoli testovací funkci  φ ,

Výsledek vyplývá ze vzorce pro newtonovský potenciál (základní řešení Poissonovy rovnice). Toto je v podstatě forma inverzního vzorce pro radonovou transformaci, protože získává hodnotu φ ( x ) ze svých integrálů nad hyperplanes. Například pokud n je liché a k = 1 , pak integrál na pravé straně je

kde ( ξ , p ) je radonová transformace φ :

Alternativní ekvivalentní vyjádření rozkladu rovinných vln od Gelfand & Shilov (1966–1968 , I, §3.10) je

pro n sudý, a

pro n liché.

Fourierova jádra

Při studiu Fourierových řad spočívá hlavní otázka v určení, zda a v jakém smyslu Fourierova řada spojená s periodickou funkcí konverguje k funkci. N tý částečný součet Fourierovy řady funkčního f o období 2 n je definován konvolucí (v intervalu [-π, π] ) s jádrem Dirichlet :

Tím pádem,

kde

Základní výsledek elementární Fourierovy řady uvádí, že jádro Dirichlet má tendenci k násobku delta funkce jako N → ∞ . To je interpretováno ve smyslu distribuce, že

pro každou kompaktně podporovanou hladkou funkci f . Formálně tedy člověk má

na intervalu  [−π, π] .

Navzdory tomu výsledek neplatí pro všechny kompaktně podporované spojité funkce: to znamená, že D N nekonverguje slabě ve smyslu opatření. Nedostatek konvergence Fourierovy řady vedl k zavedení různých metod sčítatelnosti, které produkují konvergenci. Metoda souhrnu Cesàra vede k jádru Fejér

Na jádra Fejér tendenci delta funkce v silnějším smyslu, že

pro každou kompaktně podporovanou spojitou funkci f . Z toho vyplývá, že Fourierova řada jakékoli spojité funkce je Cesàro sčítatelná s hodnotou funkce v každém bodě.

Hilbertova vesmírná teorie

Distribuce Dirac delta je hustě definovaný neomezený lineární funkční na Hilbertově prostoru L 2 ze čtvercového integrovatelných funkcí . Hladké kompaktně podporované funkce jsou ve skutečnosti husté v L 2 a působení distribuce delta na takové funkce je dobře definováno. V mnoha aplikacích je možné identifikovat podprostory L 2 a poskytnout silnější topologii, na které funkce delta definuje ohraničenou lineární funkci .

Sobolevovy prostory

Sobolev vkládání teorém pro Soboleových prostorů na reálné ose R znamená, že jakékoliv čtverečních integrovatelná funkce f tak, že

je automaticky kontinuální a zejména uspokojuje

Tak δ je omezená lineární funkční na Sobolev prostoru H 1 . Ekvivalentně δ je prvek spojitého duálního prostoru H −1 z H 1 . Obecněji řečeno, v n rozměrech má jeden δH - s ( R n ) za předpokladu  s > n  / 2 .

Prostory holomorfních funkcí

V komplexní analýze vstupuje funkce delta prostřednictvím Cauchyho integrálního vzorce , který tvrdí, že pokud D je doména v komplexní rovině s hladkou hranicí, pak

pro všechny holomorfních funkcí F v D , které jsou kontinuální o uzavření D . Výsledkem je, že delta funkce δ z je v této třídě holomorfních funkcí reprezentována Cauchyovým integrálem:

Kromě toho, ať H 2 (∂ D ) je Hardy prostor se skládá z uzávěru v L 2 (∂ D ) všech holomorfních funkcí v D kontinuální až do hranice D . Pak se funkce v H 2 (∂ D ) jednoznačně rozšíří na holomorfní funkce v D a Cauchyův integrální vzorec nadále platí. Zejména pro zD je delta funkce δ z spojitá lineární funkční na H 2 (∂ D ). To je zvláštní případ situace z více komplexních proměnných , ve kterých, na hladké domény D je Szegö jádro hraje roli Cauchy integrálu.

Řešení identity

Vzhledem k tomu, kompletní ortonormální báze sadu funkcí { φ n } v Hilbertově prostoru, například, normalizované vektory příslušníky kompaktní operátora vlastním adjoint , jakýkoliv vektor, f může být vyjádřena jako

Koeficienty {α n } se nacházejí jako

který může být reprezentován zápisem:

forma Bracetova zápisu Diraca. Přijetím tohoto zápisu, expanze f vezme Dyadická tvar:

Nechal jsem označují operátor identity na Hilbert prostoru, výraz

se nazývá řešení identity . Když je Hilbertův prostor prostor L 2 ( D ) čtvercově integrovatelných funkcí v doméně D , množství:

je integrální operátor a výraz pro f lze přepsat

Pravá strana konverguje k f ve smyslu L 2 . Nemusí to platit v bodovém smyslu, i když f je spojitá funkce. Přesto je běžné zneužívat notaci a psát

což má za následek reprezentaci funkce delta:

Při vhodném upraveném Hilbertově prostoru (Φ, L 2 ( D ), Φ*), kde Φ ⊂ L 2 ( D ) obsahuje všechny kompaktně podporované hladké funkce, se toto součty mohou sbíhat v Φ*, v závislosti na vlastnostech báze φ n . Ve většině případů praktického zájmu pochází ortonormální základ od integrálního nebo diferenciálního operátoru, v takovém případě se řada sbíhá ve smyslu distribuce .

Nekonečně malé funkce delta

Cauchy použil nekonečně malý α k zapsání jednotkového impulsu, nekonečně vysoké a úzké delta funkce Diracového typu δ α vyhovující v řadě článků v roce 1827. Cauchy definoval v Cours d'Analyse (1827) nekonečně malé z hlediska tendence sekvence na nulu. Totiž, taková nulová sekvence se stane nekonečně malou v terminologii Cauchyho a Lazare Carnota .

Nestandardní analýza umožňuje přísně zacházet s infinitesimals. Článek Yamashity (2007) obsahuje bibliografii o moderních Diracových delta funkcích v kontextu nekonečně malého obohaceného kontinua poskytovaného hyperrealy . Zde může být Diracova delta dána skutečnou funkcí, mající tu vlastnost, že pro každou skutečnou funkci F má, jak předpokládali Fourier a Cauchy.

Hřeben Dirac

Hřeben Dirac je nekonečná řada funkcí Dirac delta rozmístěných v intervalech T

Takzvaný jednotný „sled pulsů“ měření Diracovy delty, který je známý jako Diracovský hřeben nebo Shahova distribuce, vytváří funkci vzorkování , často používanou při zpracování digitálního signálu (DSP) a při analýze diskrétního časového signálu. Hřeben Dirac je uveden jako nekonečný součet , jehož limit je chápán ve smyslu distribuce,

což je posloupnost bodových hmot na každém z celých čísel.

Až do celkové normalizační konstanty je Diracov hřeben roven vlastní Fourierově transformaci. To je důležité, protože pokud je nějaká funkce Schwartz , pak periodizace ze je dán konvolucí

Zejména,

je přesně Poissonův součtový vzorec . Obecněji platí, že tento vzorec zůstává pravdivý, pokud jde o temperovanou distribuci rychlého klesání nebo, ekvivalentně, pokud jde o pomalu rostoucí, běžnou funkci v prostoru temperovaných distribucí.

Sokhotski – Plemeljova věta

Sokhotski-Plemelj věta , významný v kvantové mechaniky, se týká funkce delta do rozvodné pv 1 / x je hlavní hodnota Cauchy funkce 1 / x , definovanou

Sokhotského vzorec to říká

Zde je limit chápán ve smyslu distribuce, že pro všechny kompaktně podporované hladké funkce f ,

Vztah k Kroneckerově deltě

Kronecker delta δ ij je množství definované

pro všechna celá čísla i , j . Tato funkce pak splňuje následující analogii prosévající vlastnosti: pokud existuje nějaká dvojnásobně nekonečná posloupnost , pak

Podobně pro jakoukoli skutečnou nebo komplexní oceňovanou spojitou funkci f na R splňuje Diracova delta prosévací vlastnost

Toto ukazuje funkci Kronecker delta jako diskrétní analog funkce Dirac delta.

Aplikace

Teorie pravděpodobnosti

V teorii pravděpodobnosti a statistikách se Diracova delta funkce často používá k reprezentaci diskrétního rozdělení nebo částečně diskrétního, částečně spojitého rozdělení pomocí funkce hustoty pravděpodobnosti (která se normálně používá k reprezentaci absolutně spojitých rozdělení). Například funkci hustoty pravděpodobnosti f ( x ) diskrétního rozdělení sestávajícího z bodů x = { x 1 , ..., x n }, s odpovídajícími pravděpodobnostmi p 1 , ..., p n , lze zapsat jako

Jako další příklad zvažte rozdělení, ve kterém 6/10 času vrací standardní normální rozdělení a 4/10 času vrací přesně hodnotu 3,5 (tj. Částečně spojité, částečně diskrétní rozdělení směsi ). Hustotní funkci této distribuce lze zapsat jako

Funkce delta se také používá k reprezentaci výsledné funkce hustoty pravděpodobnosti náhodné proměnné, která je transformována spojitě diferencovatelnou funkcí. Pokud Y = g ( X ) je spojitá diferencovatelná funkce, pak lze hustotu Y zapsat jako

Funkce delta se používá také v úplně jiným způsobem reprezentovat místní čas o difuzního procesu (jako Brownova pohybu ). Místní čas stochastického procesu B ( t ) je dán vztahem

a představuje množství času, které proces stráví v bodě x v rozsahu procesu. Přesněji, v jedné dimenzi lze tento integrál zapsat

kde 1 [ x - ε , x + ε ] je indikátorová funkce intervalu [ x - ε , x + ε ] .

Kvantová mechanika

Funkce delta je v kvantové mechanice účelná . Vlnová funkce částice dává amplitudu pravděpodobnost nalezení částice v dané oblasti prostoru. Vlnové funkce se předpokládá, že prvky Hilbertova prostoru L 2 ze čtvercových-integrovatelných funkcí , a celková pravděpodobnost nalezení částice v daném intervalu je integrál velikosti vlnové funkce na druhou přes interval. Sada { } vlnových funkcí je ortonormální, pokud jsou normalizovány

kde je Kroneckerova delta. Sada funkcí ortonormálních vln je v prostoru čtvercově integrovatelných funkcí úplná, pokud lze jakoukoli vlnovou funkci vyjádřit jako lineární kombinaci { } s komplexními koeficienty:

s . Kompletní ortonormální systémy vlnové funkce se objeví přirozeně jako eigenfunctions z hamiltoniánu (jednoho vázaného systému ) v kvantové mechanice, že opatření na snížení energetické hladiny, které se nazývají vlastní čísla. Soubor vlastních čísel je v tomto případě znám jako spektrum hamiltoniánu. V notaci bra -ket , jak je uvedeno výše , tato rovnost znamená vyřešení identity:

Zde se předpokládá, že vlastní čísla jsou diskrétní, ale množina vlastních čísel pozorovatelného může být spojitá spíše než diskrétní. Příkladem je pozorovatelná poloha , ( x ) = x ψ ( x ) . Spektrum polohy (v jedné dimenzi) je celá skutečná čára a nazývá se spojité spektrum . Na rozdíl od hamiltoniánu však operátorovi polohy chybí správné vlastní funkce. Konvenčním způsobem, jak tento nedostatek překonat, je rozšířit třídu dostupných funkcí tím, že umožníme také distribuce: to znamená nahradit Hilbertův prostor kvantové mechaniky vhodným zmanipulovaným Hilbertovým prostorem . V této souvislosti má operátor polohy kompletní sadu vlastních distribucí, označených body y skutečné čáry, dané

Vlastní funkce polohy jsou označeny v Diracově notaci a jsou známé jako vlastní stavy polohy.

Podobné úvahy platí pro vlastní stavy operátora hybnosti nebo pro jakýkoli jiný nezávislý nevázaný operátor P v Hilbertově prostoru za předpokladu, že spektrum P je spojité a neexistují zdegenerovaná vlastní čísla. V takovém případě existuje množina Ω reálných čísel (spektrum) a sbírka φ y distribucí indexovaných prvky Ω, takže

To znamená, že φ y jsou vlastní vektory P . Pokud jsou vlastní vektory normalizovány tak, že

ve smyslu rozdělení pak pro jakoukoli testovací funkci ψ,

kde

To znamená, že stejně jako v diskrétním případě existuje řešení identity

kde integrál oceňovaný operátorem je opět chápán ve slabém smyslu. Pokud má spektrum P spojité i diskrétní části, pak rozlišení identity zahrnuje součet přes diskrétní spektrum a integrál přes spojité spektrum.

Funkce delta má také mnohem více specializovaných aplikací v kvantové mechanice, například modely delta potenciálu pro jamku s jednoduchým a dvojitým potenciálem.

Strukturální mechanika

Funkci delta lze ve strukturní mechanice použít k popisu přechodových nebo bodových zatížení působících na konstrukce. Lze zapsat řídící rovnici jednoduché soustavy hmota - pružina buzené náhlým silovým impulzem I v čase t = 0

kde m je hmotnost, £ vychýlení a k na pružinovou konstantu .

Jako další příklad je rovnice řídící statickou výchylku štíhlého paprsku podle Eulerovy -Bernoulliho teorie ,

kde EI je ohybová tuhost paprsku, w k vychýlení , x prostorové souřadnici a q ( x ) na rozložení zatížení. Pokud je paprsek zatížen bodovou silou F při x = x 0 , je zapsáno rozdělení zatížení

Protože integrace delta funkce má za následek krokovou funkci Heaviside , vyplývá z toho, že statická výchylka štíhlého paprsku vystaveného vícenásobnému bodovému zatížení je popsána množinou po částech polynomů .

Bodový moment působící na paprsek lze také popsat pomocí funkcí delta. Uvažujme dvě protichůdné bodové síly F ve vzdálenosti d od sebe. Poté vytvoří moment M = Fd působící na paprsek. Nyní je vzdálenost d přiblíží limitu nula, přičemž M je udržována konstantní. Rozložení zatížení za předpokladu, že moment ve směru hodinových ručiček působí na x = 0, je zapsáno

Bodové momenty tak mohou být reprezentovány derivací funkce delta. Integrace rovnice paprsku opět vede k polynomiální výchylce po částech .

Viz také

Poznámky

Reference

externí odkazy