Hermitovy polynomy - Hermite polynomials

V matematice jsou Hermitovy polynomy klasickou ortogonální polynomickou sekvencí .

Polynomy vznikají v:

Hermitské polynomy byly definovány Pierrem-Simonem Laplaceem v roce 1810, i když ve sotva rozpoznatelné formě, a podrobně je studoval Pafnuty Chebyshev v roce 1859. Chebyshevova práce byla přehlédnuta a byly pojmenovány později podle Charlese Hermita , který napsal o polynomech v roce 1864, popisovat je jako nové. V důsledku toho nebyly nové, ačkoli Hermite byl první, kdo definoval vícerozměrné polynomy ve svých pozdějších 1865 publikacích.

Definice

Stejně jako ostatní klasické ortogonální polynomy lze Hermitovy polynomy definovat z několika různých výchozích bodů. Hned od začátku, že existují dvě různé standardizace v běžném používání, jedna pohodlná metoda je následující:

  • V „pravděpodobnostních je hermitovy polynomy“ jsou dány
  • zatímco „fyzikální Hermitovy polynomy“ jsou dány znakem

Tyto rovnice mají formu Rodriguesova vzorce a lze je také zapsat jako

Tyto dvě definice nejsou přesně totožné; každý je změnou měřítka druhého:

Jedná se o Hermitovy polynomické sekvence různých variant; viz materiál o odchylkách níže.

Zápis He a H je ten, který se používá ve standardních referencích. Polynomy He n jsou někdy označeny H n , zvláště v teorii pravděpodobnosti, protože

je funkce hustoty pravděpodobnosti pro normální rozdělení s očekávanou hodnotou 0 a standardní odchylkou 1.

Prvních šest pravděpodobnostních Hermitových polynomů He n ( x )
  • Prvních jedenáct hermitských polynomů pravděpodobnosti je:
Prvních šest (fyzikálních) hermitských polynomů H n ( x )
  • Hermitovy polynomy prvních jedenácti fyziků jsou:

Vlastnosti

N tého řádu Hermiteův polynom je polynom stupně n . Pravděpodobnostní verze He n má vedoucí koeficient 1, zatímco fyzikova verze H n má vedoucí koeficient 2 n .

Ortogonalita

H n ( x ) a on n ( x ) jsou n th-stupeň polynomy pro n = 0, 1, 2, 3, ... . Tyto polynomy jsou ortogonální vzhledem k váhové funkci ( míra )

nebo

tj. máme

Kromě toho,

nebo

kde je Kroneckerova delta .

Pravděpodobnostní polynomy jsou tedy ortogonální vzhledem ke standardní normální hustotě pravděpodobnosti.

Úplnost

Tyto hermitovy polynomy (pravděpodobnostních je nebo fyzika) tvoří ortogonální základnu na Hilbertově prostoru funkcí, který by splňoval

ve kterém je vnitřní produkt dán integrálem

včetně Gaussovy váhové funkce w ( x ) definované v předchozí části

Ortogonální základ pro L 2 ( R , w ( x ) dx ) je kompletní ortogonální systém . U ortogonálního systému je úplnost ekvivalentní skutečnosti, že funkce 0 je jedinou funkcí fL 2 ( R , w ( x ) dx ) ortogonální ke všem funkcím v systému.

Protože lineární rozpětí Hermitových polynomů je prostorem všech polynomů, je třeba ukázat (v případě fyzika), že pokud f splňuje

pro každé n ≥ 0 , pak f = 0 .

Jedním z možných způsobů, jak toho dosáhnout, je ocenit celou funkci

zmizí stejně. Skutečnost, pak, že F ( to ) = 0 pro každé reálné t prostředky, že Fourierova transformace z f ( x ), e - x 2 je 0, tedy f je 0 téměř všude. Varianty výše uvedeného důkazu úplnosti platí pro jiné váhy s exponenciálním úpadkem.

V případě Hermite je také možné prokázat explicitní identitu, která implikuje úplnost (viz část o vztahu Úplnost níže).

Ekvivalentní formulace skutečnosti, že hermitské polynomy jsou ortogonálním základem pro L 2 ( R , w ( x ) dx ), spočívá v zavedení hermitských funkcí (viz níže) a v tom, že hermitské funkce jsou ortonormálním základem pro L 2 ( R ) .

Hermitova diferenciální rovnice

Hermitovy polynomy s pravděpodobností jsou řešením diferenciální rovnice

kde λ je konstanta. Uložením okrajové podmínky, že u by mělo být polynomiálně ohraničeno v nekonečnu, má rovnice řešení pouze v případě, že λ je nezáporné celé číslo a řešení je jednoznačně dáno vztahem , kde označuje konstantu.

Přepis diferenciální rovnice jako problém vlastní hodnoty

Hermitovy polynomy lze chápat jako vlastní funkce diferenciálního operátoru . Tento problém vlastních čísel se nazývá Hermitova rovnice , i když se tento výraz používá také pro úzce související rovnici

jehož řešení je jednoznačně dáno fyzikálními hermitskými polynomy ve formě , kde označuje konstantu, po uložení okrajové podmínky, že u by mělo být polynomiálně ohraničeno v nekonečnu.

Obecná řešení výše uvedených diferenciálních rovnic druhého řádu jsou ve skutečnosti lineární kombinace obou Hermitových polynomů a splývajících hypergeometrických funkcí prvního druhu. Například pro fyzikovu Hermitovu rovnici

obecné řešení má formu

kde a jsou konstanty, jsou fyzikální Hermitovy polynomy (prvního druhu) a fyzikální Hermitovy funkce (druhého druhu). Posledně jmenované funkce jsou kompaktně znázorněny, kde jsou konfluentní hypergeometrické funkce prvního druhu . Konvenční Hermitovy polynomy mohou být také vyjádřeny pomocí splývajících hypergeometrických funkcí, viz níže.

S obecnějšími okrajovými podmínkami lze Hermitovy polynomy zobecnit, aby získaly obecnější analytické funkce pro komplexně hodnocené λ . Je také možný explicitní vzorec Hermitových polynomů z hlediska obrysových integrálů ( Courant & Hilbert 1989 ).

Vztah opakování

Pořadí pravděpodobnostních Hermitových polynomů také splňuje relaci recidivy

Jednotlivé koeficienty souvisejí podle následujícího rekurzního vzorce:

a a 0,0 = 1 , a 1,0 = 0 , a 1,1 = 1 .

U fyzikálních polynomů za předpokladu

my máme

Jednotlivé koeficienty souvisejí podle následujícího rekurzního vzorce:

a a 0,0 = 1 , a 1,0 = 0 , a 1,1 = 2 .

Hermitovy polynomy tvoří Appellovu sekvenci , tj. Jsou to polynomické sekvence splňující identitu

Ekvivalentně tím, že Taylor expanduje ,

Tyto umbrální identity jsou samozřejmé a jsou zahrnuty v níže uvedené diferenciální reprezentaci operátora ,

V důsledku toho platí pro m th deriváty následující vztahy:

Z toho vyplývá, že Hermitovy polynomy také splňují relaci recidivy

Tyto poslední vztahy spolu s počátečními polynomy H 0 ( x ) a H 1 ( x ) lze v praxi použít k rychlému výpočtu polynomů.

Turánovy nerovnosti jsou

Navíc platí následující věta o násobení :

Explicitní výraz

Fyzikovy hermitské polynomy lze zapsat výslovně jako

Tyto dvě rovnice lze spojit do jedné pomocí funkce floor :

Pravděpodobnostní Hermitovy polynomy podobné vzorce, které z nich lze získat nahrazením síly 2 x odpovídající sílou 2 x a vynásobením celého součtu 2 - n/2:

Inverzní explicitní výraz

Inverzní výše uvedených explicitních výrazů, to znamená ty, pro monomials z hlediska Hermiteovými polynomů pravděpodobnostních je mu jsou

Odpovídající výrazy pro fyzikální Hermitovy polynomy H následují přímo správným změnou měřítka:

Generující funkce

Hermitovy polynomy jsou dány funkcí exponenciálního generování

Tato rovnost platí pro všechny komplexní hodnoty x a t a lze ji získat zápisem Taylorovy expanze na x celé funkce ze - z 2 (v případě fyzika). Lze také odvodit (fyzikovu) generující funkci pomocí Cauchyova integrálního vzorce pro zápis Hermitových polynomů jako

Použití v součtu

lze vyhodnotit zbývající integrál pomocí počtu zbytků a dospět k požadované generující funkci.

Očekávané hodnoty

Pokud X je náhodná veličina s normálním rozložením se standardní odchylkou 1 a očekávanou hodnotou μ , pak

Momenty standardní normály (s očekávanou hodnotou nula) lze odečíst přímo ze vztahu pro sudé indexy:

kde (2 n - 1) !! je dvojí faktoriál . Všimněte si, že výše uvedený výraz je zvláštním případem znázornění pravděpodobnostních Hermitových polynomů jako momentů:

Asymptotická expanze

Asymptoticky, jako n → ∞ , expanze

platí. V určitých případech týkajících se širšího rozsahu hodnocení je nutné zahrnout faktor pro změnu amplitudy:

které pomocí Stirlingovy aproximace lze dále v limitu dále zjednodušit na

Tato expanze je potřeba vyřešit wavefunction o quantum harmonický oscilátor tak, že souhlasí s klasickým aproximace v limitu principu korespondence .

Lepší přiblížení, které odpovídá odchylkám ve frekvenci, je dáno pomocí

Jemnější aproximace, která bere v úvahu nerovnoměrné rozestupy nul poblíž okrajů, využívá substituci

se kterým má jednotnou aproximaci

Podobné aproximace platí pro monotónní a přechodové oblasti. Konkrétně pokud

pak

zatímco pro

s t komplexním a ohraničeným, aproximace je

kde Ai je vzdušná funkce prvního druhu.

Zvláštní hodnoty

Hermitovy polynomy fyzika hodnocené s nulovým argumentem H n (0) se nazývají hermitská čísla .

které splňují rekurzivní vztah H n (0) = −2 ( n - 1) H n - 2 (0) .

Z hlediska polynomů pravděpodobnosti se to překládá do

Vztahy k dalším funkcím

Laguerreovy polynomy

Hermitovy polynomy lze vyjádřit jako speciální případ Laguerrových polynomů :

Vztah ke splývajícím hypergeometrickým funkcím

Hermitovy polynomy fyzika lze vyjádřit jako speciální případ funkcí parabolických válců :

v pravé polorovině , kde U ( a , b , z ) je Tricomiho splývající hypergeometrická funkce . Podobně,

kde 1 F 1 ( a , b ; z ) = M ( a , b ; z ) je Kummerova konfluentní hypergeometrická funkce .

Reprezentace diferenciálního operátoru

Pravděpodobnostní Hermitovy polynomy uspokojují identitu

kde D představuje diferenciaci vzhledem k x a exponenciál je interpretován jejím rozšířením jako mocninová řada . Když tato řada funguje na polynomech, neexistují žádné delikátní otázky konvergence, protože téměř všechny termíny zmizí.

Protože koeficienty mocninných řad exponenciálu jsou dobře známy a deriváty vyšších řádů monomiální x n lze explicitně zapsat, tato reprezentace diferenciálního operátoru dává vznik konkrétnímu vzorci pro koeficienty H n, které lze použít. k rychlému výpočtu těchto polynomů.

Protože formální výraz pro Weierstrassovu transformaci W je e D 2 , vidíme, že Weierstrassova transformace ( 2 ) n He n (X/2) je x n . Weierstrassova transformace tak v zásadě mění řadu Hermitových polynomů na odpovídající Maclaurinovu řadu .

Existence nějaké formální mocninové řady g ( D ) s nenulovým konstantním koeficientem, jako je He n ( x ) = g ( D ) x n , je dalším ekvivalentem tvrzení, že tyto polynomy tvoří Appellovu sekvenci . Vzhledem k tomu, že jsou sekvence Appell, jsou tím spíše Sheffer sekvence .

Obrysově integrální zobrazení

Ze zobrazení generující funkce výše vidíme, že Hermitovy polynomy mají reprezentaci z hlediska obrysového integrálu , jako

s obrysem obklopujícím původ.

Zobecnění

Hermitovy polynomy pravděpodobnosti definované výše jsou ortogonální s ohledem na standardní normální rozdělení pravděpodobnosti, jehož funkce hustoty je

který má očekávanou hodnotu 0 a rozptyl 1.

Při škálování lze analogicky hovořit o generalizovaných hermitských polynomech

rozptylu α , kde α je jakékoli kladné číslo. Ty jsou pak ortogonální vzhledem k normálnímu rozdělení pravděpodobnosti, jehož hustotní funkce je

Jsou dány

Teď když

potom polynomická posloupnost, jejíž n -tý člen je

se nazývá umbrální složení dvou polynomických sekvencí. Může být prokázáno, že uspokojuje identity

a

Poslední identita je vyjádřena tím, že tato parametrizovaná rodina polynomických sekvencí je známá jako křížová sekvence. (Viz výše uvedený oddíl o sekvencích Appell a o reprezentaci diferenciálního operátoru , což vede k jeho snadné derivaci. Tato binomická identita typu pro α = β =1/2, již bylo zaznamenáno ve výše uvedené části o vztazích #Recursion .)

"Negativní rozptyl"

Vzhledem k tomu, že polynomické sekvence tvoří skupinu pod působením umbrálního složení , lze je označit pomocí

posloupnost, která je inverzní k té, která je podobně označována, ale bez znaménka mínus, a tedy hovoří o hermitských polynomech negativní variance. Pro α> 0 jsou koeficienty hodnoty jen absolutní hodnoty odpovídajících koeficientů .

Ty vznikají jako momenty normálního rozdělení pravděpodobnosti: n -tý moment normálního rozdělení s očekávanou hodnotou μ a rozptylem σ 2 je

kde X je náhodná proměnná se zadaným normálním rozložením. Zvláštní případ identity křížové sekvence to pak říká

Aplikace

Hermitské funkce

Hermitovy funkce (často nazývané hermitovsko-gaussovské funkce) lze definovat z fyzikálních polynomů:

Tím pádem,

Protože tyto funkce obsahují druhou odmocninu váhové funkce a byly vhodně upraveny, jsou ortonormální :

a tvoří ortonormální základ L 2 ( R ) . Tato skutečnost je ekvivalentní odpovídajícímu tvrzení pro Hermitovy polynomy (viz výše).

Funkce Hermite úzce souvisí s funkcí Whittaker ( Whittaker & Watson 1996 ) D n ( z ) :

a tím k dalším parabolickým funkcím válce .

Funkce Hermite splňují diferenciální rovnici

Tato rovnice je ekvivalentní Schrödingerově rovnici pro harmonický oscilátor v kvantové mechanice, takže tyto funkce jsou vlastní funkce .

Hermitské funkce: 0 (modrá, plná), 1 (oranžová, čárkovaná), 2 (zelená, tečkovaná), 3 (červená, tečkovaná), 4 (fialová, plná) a 5 (hnědá, přerušovaná)
Hermitské funkce: 0 (modrá, plná), 2 (oranžová, čárkovaná), 4 (zelená, tečkovaná) a 50 (červená, plná)

Rekurzivní vztah

Po rekurzivních vztazích hermitských polynomů se funkce hermitů řídí

a

Rozšíření prvního vztahu na libovolné m th deriváty pro jakékoli kladné celé číslo m vede k

Tento vzorec lze použít ve spojení s relacemi opakování pro He n a ψ n pro efektivní výpočet jakékoli derivace Hermitových funkcí.

Cramérova nerovnost

Pro skutečné x splňují funkce Hermite následující mez kvůli Haraldovi Cramérovi a Jacku Indritzovi:

Hermit funguje jako vlastní funkce Fourierovy transformace

The Hermitovy FUNKCE ln n ( x ) je sada vlastních funkcí v kontinuální Fourierovy transformace F . Chcete -li to vidět, vezměte fyzickou verzi generující funkce a vynásobte e -1/2x 2 . To dává

Fourierova transformace levé strany je dána vztahem

Fourierova transformace pravé strany je dána vztahem

Rovnice stejných mocnin t v transformovaných verzích levé a pravé strany nakonec přináší

Hermitské funkce ψ n ( x ) jsou tedy ortonormálním základem L 2 ( R ) , který diagonalizuje operátor Fourierovy transformace .

Wignerovy distribuce funkcí Hermite

Distribuční funkce Wigner z n funkce tého řádu Hermiteho souvisí s n tého řádu Laguerre polynomu . Laguerreovy polynomy jsou

vedoucí k funkcím oscilátoru Laguerre

Pro všechna přirozená celá čísla n je to snadné vidět

kde Wignerovo rozdělení funkce xL 2 ( R , C ) je definováno jako

Toto je zásadní výsledek kvantového harmonického oscilátoru, který objevil Hip Groenewold v roce 1946 ve své disertační práci. Je to standardní paradigma kvantové mechaniky ve fázovém prostoru .

Mezi oběma rodinami polynomů existují další vztahy .

Kombinatorická interpretace koeficientů

V Hermitově polynomu He n ( x ) rozptylu 1 je absolutní hodnota koeficientu x k počet (neuspořádaných) oddílů n -prvku nastaveného na k singletů an - k/2(neuspořádané) páry. Ekvivalentně je to počet involucí n -prvkové množiny s přesně k pevnými body, nebo jinými slovy, počet shod v kompletním grafu na n vrcholech, který ponechá k vrcholů odkrytých (ve skutečnosti jsou Hermitovy polynomy shodné polynomy těchto grafů). Součet absolutních hodnot koeficientů udává celkový počet přepážek na singletony a páry, takzvaná telefonní čísla

1, 1, 2, 4, 10, 26, 76, 232, 764, 2620, 9496, ... (sekvence A000085 v OEIS ).

Tato kombinatorická interpretace může být vztažena k úplným exponenciálním Bellovým polynomům jako

kde x i = 0 pro všechny i > 2 .

Tato čísla mohou být také vyjádřena jako speciální hodnota Hermitových polynomů:

Vztah úplnosti

Christoffel-Darboux vzorec pro Hermiteovými polynomy čte

Následující identita úplnosti výše uvedených funkcí Hermite navíc platí ve smyslu distribucí :

kde δ je Diracova delta funkce , ψ n Hermitské funkce a δ ( x - y ) představuje Lebesgueovu míru na přímce y = x v R 2 , normalizovanou tak, že její projekce na horizontální osu je obvyklou Lebesgueovou mírou.

Tato distribuční identita následuje podle Wienera (1958) tím, že vezme u → 1 v Mehlerově vzorci , platném při −1 < u <1 :

který je často ekvivalentně uveden jako oddělitelné jádro,

Funkce ( x , y ) → E ( x , y ; u ) je bivariační hustota Gaussovy pravděpodobnosti na R 2 , což je, když u je blízko 1, velmi koncentrovaná kolem přímky y = x a velmi rozložená na ten řádek. Z toho vyplývá, že

když jsou f a g spojité a kompaktně podporované.

To poskytuje, že f může být vyjádřeno v Hermitových funkcích jako součet řady vektorů v L 2 ( R ) , jmenovitě,

Za účelem prokázání výše rovnosti E ( x , y ; u ) je Fourierova transformace z Gaussova funkce se používá opakovaně:

Hermitův polynom je pak reprezentován jako

S touto reprezentací pro H n ( x ) a H n ( y ) je evidentní, že

a to dává požadované rozlišení výsledku identity, opět za použití Fourierovy transformace Gaussových jader

Viz také

Poznámky

Reference

externí odkazy