Hermitovy polynomy

Technology
12 hours ago
8
4
2
Avatar
Author
Albert Flores

Hermitovy polynomy jsou v matematice klasická posloupnost ortogonálních polynomů.

Hermitovy polynomy se objevují: * ve zpracování signálu jako Hermitovské vlnky pro analýzu vlnkovou transformací * v pravděpodobnosti, jako například Edgeworthova řada nebo v souvislosti s Brownovým pohybem; * v kombinatorice, jako příklad Appellovy posloupnosti, která řídí stínový počet; * v numerické matematice jako Gaussovo kvadraturní pravidlo; * ve fyzice, kde popisují kvantové stavy kvantového harmonického oscilátoru; * v teorii systémů ve spojení s nelineárními operacemi na Gaussovském šumu. * v náhodných maticích v Gaussovských náhodných maticích.

Hermitovy polynomy definoval (i když v sotva rozpoznatelném formě) Pierre-Simon Laplace v roce 1810; detailně je zkoumal Pafnutij Lvovič Čebyšev v roce 1859. Čebyševova práce však byla přehlížena a polynomy byly později pojmenovány po Charlesu Hermitovi, který je popsal jako nové v roce 1864. +more Hermite tyto polynomy tedy neobjevil jako první, ale ve svém pozdějším díle z roku 1865 definoval vícerozměrné polynomy.

Definice

Hermitovy polynomy je možné definovat stejně jako jiné klasické ortogonální polynomy několika způsoby. Je třeba si uvědomit, že existují dvě různé definice, přičemž obvyklejší je tato:

* „pravděpodobnostní Hermitovy polynomy“ jsou dány vzorcem :: \mathit{He}_n(x) = (-1)^n e^{\frac{x^2}{2}}\frac{d^n}{\mbox{d}x^n}e^{-\frac{x^2}{2}},

* zatímco „fyzikální Hermitovy polynomy“ jsou dány vzorcem :: H_n(x) = (-1)^n e^{x^2}\frac{d^n}{\mbox{d}x^n}e^{-x^2}.

Tyto rovnice mají tvar Rodriguesova vzorce a zapisují se také jako

: \mathit{He}_n(x) = \left(x - \frac{\mbox{d}}{\mbox{d}x} \right)^n \cdot 1, \quad H_n(x) = \left(2x - \frac{\mbox{d}}{\mbox{d}x} \right)^n \cdot 1.

tyto dvě definice ovšem nejsou přesně identické; liší se použitím jiného měřítka:

:H_n(x)=2^\frac{n}{2} \mathit{He}_n\left(\sqrt{2} \,x\right), \quad \mathit{He}_n(x)=2^{-\frac{n}{2}} H_n\left(\frac {x}{\sqrt 2} \right).

Tyto posloupnosti Hermitových polynomů se liší rozptylem; viz výklad o variancích níže.

Obvykle se pravděpodobnostní a fyzikální Hermitovy polynomy rozlišují značením a , v teorii pravděpodobnosti se však často místo používá , protože : \frac{1}{\sqrt{2\pi}}e^{-\frac{x^2}{2}} je hustota pravděpodobnosti pro normální rozdělení se střední hodnotou 0 a směrodatnou odchylkou 1.

Prvních šest pravděpodobnostních Hermitových polynomů

* Prvních jedenáct pravděpodobnostních Hermitových polynomů: :\begin{align} \mathit{He}_0(x) &= 1, \\ \mathit{He}_1(x) &= x, \\ \mathit{He}_2(x) &= x^2 - 1, \\ \mathit{He}_3(x) &= x^3 - 3x, \\ \mathit{He}_4(x) &= x^4 - 6x^2 + 3, \\ \mathit{He}_5(x) &= x^5 - 10x^3 + 15x, \\ \mathit{He}_6(x) &= x^6 - 15x^4 + 45x^2 - 15, \\ \mathit{He}_7(x) &= x^7 - 21x^5 + 105x^3 - 105x, \\ \mathit{He}_8(x) &= x^8 - 28x^6 + 210x^4 - 420x^2 + 105, \\ \mathit{He}_9(x) &= x^9 - 36x^7 + 378x^5 - 1260x^3 + 945x, \\ \mathit{He}_{10}(x) &= x^{10} - 45x^8 + 630x^6 - 3150x^4 + 4725x^2 - 945. \end{align}

Prvních šest (fyzikálních) Hermitových polynomů

* Prvních jedenáct fyzikálních Hermitových polynomů: :\begin{align} H_0(x) &= 1, \\ H_1(x) &= 2x, \\ H_2(x) &= 4x^2 - 2, \\ H_3(x) &= 8x^3 - 12x, \\ H_4(x) &= 16x^4 - 48x^2 + 12, \\ H_5(x) &= 32x^5 - 160x^3 + 120x, \\ H_6(x) &= 64x^6 - 480x^4 + 720x^2 - 120, \\ H_7(x) &= 128x^7 - 1344x^5 + 3360x^3 - 1680x, \\ H_8(x) &= 256x^8 - 3584x^6 + 13440x^4 - 13440x^2 + 1680, \\ H_9(x) &= 512x^9 - 9216x^7 + 48384x^5 - 80640x^3 + 30240x, \\ H_{10}(x) &= 1024x^{10} - 23040x^8 + 161280x^6 - 403200x^4 + 302400x^2 - 30240. \end{align}

Vlastnosti

Hermitův polynom -tého řádu je polynom stupně . Pravděpodobnostní verze má vedoucí koeficient 1, zatímco fyzikální verze má úvodní koeficient .

Symetrie

Z Rodriguesova vzorce uvedeného výše je vidět, že a jsou sudé nebo liché funkce podle : : H_n(-x)=(-1)^nH_n(x),\quad \mathit{He}_n(-x)=(-1)^n\mathit{He}_n(x).

Ortogonalita

a jsou polynomy -tého stupně pro . Tyto polynomy jsou ortogonální vzhledem k váhové funkci (míře)

: w(x) = e^{-\frac{x^2}{2}} \quad (\text{pro }\mathit{He})

nebo

: w(x) = e^{-x^2} \quad (\text{pro } H),

tj. máme : \int_{-\infty}^\infty H_m(x) H_n(x)\, w(x) \,\mbox{d}x = 0 \quad \text{pro každé }m \neq n.

a navíc

: \int_{-\infty}^\infty \mathit{He}_m(x) \mathit{He}_n(x)\, e^{-\frac{x^2}{2}} \,\mbox{d}x = \sqrt{2 \pi}\, n!\, \delta_{nm},

nebo

:\int_{-\infty}^\infty H_m(x) H_n(x)\, e^{-x^2} \,\mbox{d}x = \sqrt{\pi}\, 2^n n!\, \delta_{nm},

kde \delta_{nm} je Kroneckerovo delta.

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

Úplnost

Hermitovy polynomy (jak pravděpodobnostní tak fyzikální) tvoří ortogonální bázi Hilbertova prostoru funkcí, které splňují : \int_{-\infty}^\infty \bigl|f(x)\bigr|^2\, w(x) \,\mbox{d}x ve kterém je vnitřní součin definován integrálem : \langle f,g\rangle = \int_{-\infty}^\infty f(x) \overline{g(x)}\, w(x) \,\mbox{d}x, kde je gaussovská váhová funkce definovaná v předchozí části.

Ortogonální báze pro tvoří úplný ortogonální systém. Pro ortogonální systém je úplnost ekvivalentní se skutečností, že nulová funkce je jedinou funkcí , která je ortogonální se všemi funkcemi v systému.

Protože lineárním obalem Hermitových polynomů je prostor všech polynomů, pro důkaz úplnosti (pro fyzikální polynomy) stačí dokázat, že pokud splňuje : \int_{-\infty}^\infty f(x) x^n e^{- x^2} \,\mbox{d}x = 0 pro každé , pak .

Jedním ze způsobů, jak to udělat, je uvědomit si, že celá funkce : F(z) = \int_{-\infty}^\infty f(x) e^{z x - x^2} \,\mbox{d}x = \sum_{n=0}^\infty \frac{z^n}{n. } \int f(x) x^n e^{- x^2} \,\mbox{d}x = 0 bude mít nulovou hodnotu identicky. +more Skutečnost, že pak bude pro každé reálné znamená, že Fourierova transformace je 0, a tedy že je 0 skoro všude. Varianty výše uvedeného důkazu úplnosti platí i pro jiné váhy s exponenciálním poklesem.

V Hermitově případě je možné dokázat i explicitní identitu, která implikuje úplnost (viz část na #Relace úplnosti|Relace úplnosti níže).

Ekvivalentně lze fakt, že Hermitovy polynomy jsou ortogonální bází pro , formulovat zavedením Hermitových funkcí (viz níže) a ukázáním, že jsou ortonormální bází pro .

Hermitova diferenciální rovnice

Pravděpodobnostní Hermitovy polynomy jsou řešením diferenciální rovnice : \left(e^{-\frac12 x^2}u'\right)' + \lambda e^{-\frac12 x^2}u = 0, kde je konstanta. Zavedením okrajové podmínky, že funkce musí být v nekonečnu polynomiálně omezená, má rovnice řešení, pouze pokud je nezáporné celé číslo, a pak je řešení jednoznačně dáno u(x) = C_1 He_\lambda(x) , kde C_{1} je konstanta.

Přepsáním diferenciální rovnice jako problém vlastních hodnot : L[u] = u - x u' = -\lambda u, Hermitovy polynomy He_\lambda(x) je možné chápat jako vlastní funkce diferenciálního operátoru L[u] . Tento problém vlastní hodnoty se nazývá 'Hermitova rovnice', i když tento termín se také používá pro blízce příbuzné rovnice : u - 2xu' = -2\lambda u. +more jejichž řešení lze, po stanovení okrajové podmínky, že musí být polynomiálně omezená v nekonečnu, jednoznačně vyjádřit pomocí fyzikálních Hermitových polynomů ve tvaru u(x) = C_1 H_\lambda(x) , kde C_{1} je konstanta.

Obecné řešení výše uvedené diferenciální rovnice druhého řádu je vlastně lineární kombinací obou Hermitových polynomů a konfluentní hypergeometrické funkce prvního druhu. Například pro fyzikální Hermitovu rovnici : u - 2xu' + 2\lambda u = 0, má obecné řešení tvar : u(x) = C_1 H_\lambda(x) + C_2 h_\lambda(x), kde C_{1} a C_{2} jsou konstanty, H_\lambda(x) jsou fyzikální Hermitovy polynomy (prvního druhu) a h_\lambda(x) jsou fyzikální Hermitovy funkce (druhého druhu). +more Tyto funkce lze kompaktně reprezentovat jako h_\lambda(x) = {}_1F_1(-\tfrac{\lambda}{2};\tfrac{1}{2};x^2) kde {}_1F_1(a;b;z) jsou konfluentní hypergeometrické funkce prvního druhu. Obvyklé Hermitovy polynomy lze také vyjádřit pomocí konfluentních hypergeometrických funkcí, viz níže.

S obecnějšími okrajovými podmínkami je možné zobecnit Hermitovy polynomy pro získání obecnějších analytických funkcí pro komplexní . Explicitní vzorec Hermitových polynomů pomocí křivkových integrálů je také možné.

Rekurentní vzorec

Posloupnost pravděpodobnostních Hermitových polynomů také vyhovuje diferenční rovnici : \mathit{He}_{n+1}(x) = x \mathit{He}_n(x) - \mathit{He}_n'(x). Vztah jednotlivých koeficientů popisuje následující rekurentní vzorec: :a_{n+1,k} = \begin{cases} - n a_{n-1,k} & k = 0, \\ a_{n,k-1} - n a_{n-1,k} & k > 0, \end{cases} a, .

Pro fyzikální polynomy, předpokládá : H_n(x) = \sum^n_{k=0} a_{n,k} x^k, máme : H_{n+1}(x) = 2xH_n(x) - H_n'(x). Vztah jednotlivých koeficientů popisuje následující rekurentní vzorec: : a_{n+1,k} = \begin{cases} - a_{n,k+1} & k = 0, \\ 2 a_{n,k-1} - (k+1)a_{n,k+1} & k > 0, \end{cases} a, .

Hermitovy polynomy tvoří Appellovu posloupnost, tj. jsou posloupností polynomů, která vyhovuje rovnici : \begin{align} \mathit{He}_n'(x) &= n\mathit{He}_{n-1}(x), \\ H_n'(x) &= 2nH_{n-1}(x). +more \end{align} Ekvivalentně, podle Taylorova rozvoje, :\begin{align} \mathit{He}_n(x+y) &= \sum_{k=0}^n \binom{n}{k}x^{n-k} \mathit{He}_{k}(y) &&= 2^{-\frac n 2} \sum_{k=0}^n \binom{n}{k} \mathit{He}_{n-k}\left(x\sqrt 2\right) \mathit{He}_k\left(y\sqrt 2\right), \\ H_n(x+y) &= \sum_{k=0}^n \binom{n}{k}H_{k}(x) (2y)^{(n-k)} &&= 2^{-\frac n 2}\cdot\sum_{k=0}^n \binom{n}{k} H_{n-k}\left(x\sqrt 2\right) H_k\left(y\sqrt 2\right). \end{align} Tyto identity stínového počtu jsou evidentní a #Zobecnění|obsažené v reprezentaci diferenciálním operátorem rozebrané níže :\begin{align} \mathit{He}_n(x) &= e^{-\frac{D^2}{2}} x^n, \\ H_n(x) &= 2^n e^{-\frac{D^2}{4}} x^n. \end{align}.

V důsledku pro -tou derivaci platí: :\begin{align} \mathit{He}_n^{(m)}(x) &= \frac{n. }{(n-m). +more} \mathit{He}_{n-m}(x) &&= m. \binom{n}{m} \mathit{He}_{n-m}(x), \\ H_n^{(m)}(x) &= 2^m \frac{n. }{(n-m). } H_{n-m}(x) &&= 2^m m. \binom{n}{m} H_{n-m}(x). \end{align}.

Odtud plyne, že Hermitovy polynomy také vyhovují diferenční rovnici :\begin{align} \mathit{He}_{n+1}(x) &= x\mathit{He}_n(x) - n\mathit{He}_{n-1}(x), \\ H_{n+1}(x) &= 2xH_n(x) - 2nH_{n-1}(x). \end{align}

Tyto poslední vzorce se spolu s počátečními polynomy a používají v praxi pro rychlé vyčíslení hodnoty polynomů.

Platí Turánovy nerovnosti: : \mathit{H}_n(x)^2 - \mathit{H}_{n-1}(x) \mathit{H}_{n+1}(x) = (n-1)! \sum_{i=0}^{n-1} \frac{2^{n-i}}{i!}\mathit{H}_i(x)^2 > 0.

a následující multiplikační věta: : \begin{align} H_n(\gamma x) &= \sum_{i=0}^{\left\lfloor \tfrac{n}{2} \right\rfloor} \gamma^{n-2i}(\gamma^2 - 1)^i \binom{n}{2i} \frac{(2i). }{i. +more} H_{n-2i}(x), \\ \mathit{He}_n(\gamma x) &= \sum_{i=0}^{\left\lfloor \tfrac{n}{2} \right\rfloor} \gamma^{n-2i}(\gamma^2 - 1)^i \binom{n}{2i} \frac{(2i). }{i. }2^{-i} \mathit{He}_{n-2i}(x). \end{align}.

Explicitní vyjádření

Fyzikální Hermitovy polynomy je možné psát explicitně jako : H_n(x) = \begin{cases} \displaystyle n. \sum_{l = 0}^{\frac{n}{2}} \frac{(-1)^{\tfrac{n}{2} - l}}{(2l). +more \left(\tfrac{n}{2} - l \right). } (2x)^{2l} & \text{pro sudé } n, \\ \displaystyle n. \sum_{l = 0}^{\frac{n-1}{2}} \frac{(-1)^{\frac{n-1}{2} - l}}{(2l + 1). \left (\frac{n-1}{2} - l \right ). } (2x)^{2l + 1} & \text{pro liché } n. \end{cases}.

Tyto dvě rovnice je možné zkombinovat do jedné pomocí funkce celá část: : H_n(x) = n. \sum_{m=0}^{\left\lfloor \tfrac{n}{2} \right\rfloor} \frac{(-1)^m}{m. +more(n - 2m). } (2x)^{n - 2m}.

Pravděpodobnostní Hermitovy polynomy mají podobné vzorce, které je možné získat z těchto nahrazením mocniny odpovídající mocninou a znásobením celého součtu výrazem : : He_n(x) = n. \sum_{m=0}^{\left\lfloor \tfrac{n}{2} \right\rfloor} \frac{(-1)^m}{m. +more(n - 2m). } \frac{x^{n - 2m}}{2^m}.

Inverzní explicitní výraz

Inverzí výše uvedených explicitních výrazů, tj. výrazů pro jednočleny v členech pravděpodobnostních Hermitových polynomů jsou : x^n = n. +more \sum_{m=0}^{\left\lfloor \tfrac{n}{2} \right\rfloor} \frac{1}{2^m m. (n-2m). } He_{n-2m}(x).

Odpovídající výrazy pro fyzikální Hermitovy polynomy zjistíme přímo správnou změnou měřítka takto: : x^n = \frac{n. }{2^n} \sum_{m=0}^{\left\lfloor \tfrac{n}{2} \right\rfloor} \frac{1}{m. +more(n-2m). } H_{n-2m}(x).

Vytvořující funkce

Hermitovy polynomy lze zadat vytvořující funkcí

:\begin{align} e^{xt - \frac12 t^2} &= \sum_{n=0}^\infty \mathit{He}_n(x) \frac{t^n}{n!}, \\ e^{2xt - t^2} &= \sum_{n=0}^\infty H_n(x) \frac{t^n}{n!}. \end{align}

Tato rovnost je platná pro všechny komplexní hodnoty a a je možné ji získat zapsáním Taylorova rozvoje v celé funkce (ve fyzikálním případě). Můžeme také odvodit (fyzikální) vytvořující funkce pomocí Cauchyův vzorec zapsat Hermitovy polynomy jako : H_n(x) = (-1)^n e^{x^2} \frac{\mbox{d}^n}{\mbox{d}x^n} e^{-x^2} = (-1)^n e^{x^2} \frac{n. +more}{2\pi i} \oint_\gamma \frac{e^{-z^2}}{(z-x)^{n+1}} \,dz.

Dosazením do součtu :\sum_{n=0}^\infty H_n(x) \frac {t^n}{n!}, je možné vyhodnotit zbývající integrál pomocí reziduového počtu a tak získat požadovanou vytvořující funkci.

Střední hodnoty

Pokud je náhodná veličina s normálním rozdělením se standardní odchylkou 1 a střední hodnotou , pak : \operatorname{\mathbb E}\left[\mathit{He}_n(X)\right] = \mu^n.

Momenty standardního normálního rozdělení (se střední hodnotou nula) je možné číst přímo z relace pro sudé indexy: : \operatorname{\mathbb E}\left[X^{2n}\right] = (-1)^n \mathit{He}_{2n}(0) = (2n-1)!!,

kde je dvojitý faktoriál. Pamatujte, že výše uvedený výraz je speciálním případem reprezentace pravděpodobnostních Hermitových polynomů jako momentů:

: \mathit{He}_n(x) = \frac{1}{\sqrt{2\pi}} \int_{-\infty}^\infty (x + iy)^n e^{-\frac{y^2}{2}} \,dy.

Asymptotický rozvoj

Asymptoticky pro , lze použít rozvoj : e^{-\frac{x^2}{2}}\cdot H_n(x) \sim \frac{2^n}{\sqrt \pi}\Gamma\left(\frac{n+1}2\right) \cos \left(x \sqrt{2 n}- \frac{n\pi}{2} \right) Pro určité případy zabývající se širším rozsahem vyhodnocování je nutné zahrnout faktor pro změnu amplitudy: : e^{-\frac{x^2}{2}}\cdot H_n(x) \sim \frac{2^n}{\sqrt \pi}\Gamma\left(\frac{n+1}2\right) \cos \left(x \sqrt{2 n}- \frac{n\pi}{2} \right)\left(1-\frac{x^2}{2n+1}\right)^{-\frac14}=\frac{2 \Gamma(n)}{\Gamma\left(\frac{n}2\right)} \cos \left(x \sqrt{2 n}- \frac{n\pi}{2} \right)\left(1-\frac{x^2}{2n+1}\right)^{-\frac14}, což lze, pomocí Stirlingova vzorce dále zjednodušit; v limitě na : e^{-\frac{x^2}{2}}\cdot H_n(x) \sim \left(\frac{2n}{e}\right)^{\frac{n}{2}} \sqrt{2} \cos \left(x \sqrt{2n}- \frac{n\pi}{2} \right)\left(1-\frac{x^2}{2n+1}\right)^{-\frac14}.

Toto rozvoj je potřebný pro řešení vlnové funkce kvantového harmonického oscilátoru tak, že souhlasí s klasickou aproximací v limitě principu korespondence.

Lepší aproximaci, která odpovídá za variaci ve frekvenci, popisuje vztah : e^{-\frac{x^2}{2}}\cdot H_n(x) \sim \left(\frac{2n}{e}\right)^{\frac{n}{2}} \sqrt{2} \cos \left(x \sqrt{2n+1-\frac{x^2}{3}}- \frac {n\pi}{2} \right)\left(1-\frac{x^2}{2n+1}\right)^{-\frac14}.

Jemnější aproximace, která bere v úvahu nestejný odstup kořenů blízko hrany, používá substituci : x = \sqrt{2n + 1}\cos(\varphi), \quad 0 se kterou máme rovnoměrnou aproximaci : e^{-\frac{x^2}{2}}\cdot H_n(x) = 2^{\frac{n}{2}+\frac14}\sqrt{n!}(\pi n)^{-\frac14}(\sin \varphi)^{-\frac12} \cdot \left(\sin\left(\frac{3\pi}{4}+\left(\frac{n}{2}+\frac{1}{4}\right)\left(\sin 2\varphi-2\varphi\right) \right)+O\left(n^{-1}\right) \right).

Podobná aproximace platí pro monotonní a přechod oblasti. Konkrétně pokud : x = \sqrt{2n+1} \cosh(\varphi), \quad 0 pak : e^{-\frac{x^2}{2}}\cdot H_n(x) = 2^{\frac{n}{2}-\frac34}\sqrt{n. +more}(\pi n)^{-\frac14}(\sinh \varphi)^{-\frac12} \cdot e^{\left(\frac{n}{2}+\frac{1}{4}\right)\left(2\varphi-\sinh 2\varphi\right)}\left(1+O\left(n^{-1}\right) \right), zatímco pro : x = \sqrt{2n + 1} + t s komplerxním a omezeným je aproximace : e^{-\frac{x^2}{2}}\cdot H_n(x) =\pi^{\frac14}2^{\frac{n}{2}+\frac14}\sqrt{n. }\, n^{-\frac{1}{12}}\left( \operatorname{Ai}\left(2^{\frac12}n^{\frac16}t\right)+ O\left(n^{-\frac23}\right) \right), kde je Airyho funkce prvního druhu.

Speciální hodnoty

Fyzikální Hermitovy polynomy vyčíslené v bodě nula se nazývají Hermitova čísla.

:H_n(0) = \begin{cases} 0 & \text{pro lichá }n, \\ (-2)^\frac{n}{2} (n-1)!! & \text{pro sudé }n, \end{cases} což vyhovuje rekurentnímu vzorci .

V členech pravděpodobnostních polynomů se převádí na :He_n(0) = \begin{cases} 0 & \text{pro liché }n, \\ (-1)^\frac{n}{2} (n-1)!! & \text{pro sudé }n. \end{cases}

Vztahy k jiným funkcím

Laguerrovy polynomy

Hermitovy polynomy lze vyjádřit jako speciální případ Laguerrových polynomů: : \begin{align} H_{2n}(x) &= (-4)^n n. L_n^{\left(-\frac12\right)}(x^2) &&= 4^n n. +more \sum_{k=0}^n (-1)^{n-k} \binom{n-\frac12}{n-k} \frac{x^{2k}}{k. }, \\ H_{2n+1}(x) &= 2(-4)^n n. x L_n^{\left(\frac12\right)}(x^2) &&= 2\cdot 4^n n. \sum_{k=0}^n (-1)^{n-k} \binom{n+\frac12}{n-k} \frac{x^{2k+1}}{k. }. \end{align}.

Vztah k konfluentním hypergeometrickým funkcím

Fyzikální Hermitovy polynomy je možné vyjádřit jako speciální případ parabolických válcových funkcí: : H_n(x) = 2^n U\left(-\tfrac12 n, \tfrac12, x^2\right) v pravé polorovině, kde je Tricomiho konfluentní hypergeometrické funkce. Podobně : \begin{align} H_{2n}(x) &= (-1)^n \frac{(2n). +more}{n. } \,_1F_1\big(-n, \tfrac12; x^2\big), \\ H_{2n+1}(x) &= (-1)^n \frac{(2n+1). }{n. }\,2x \,_1F_1\big(-n, \tfrac32; x^2\big), \end{align} kde je Kummerova konfluentní hypergeometrické funkce.

Reprezentace diferenciálním operátorem

Pravděpodobnostní Hermitovy polynomy vyhovují vztahu : \mathit{He}_n(x) = e^{-\frac{D^2}{2}}x^n, kde reprezentuje derivaci podle a exponenciální funkce je interpretována svým rozvojem na mocninnou řadu. O konvergenci této řady aplikované na polynomy není pochyb, protože všechny členy až na konečný počet zanikají.

Protože koeficienty mocninné řady exponenciální funkce jsou známé a derivace vyššího řádu jednočlenu je možné zapsat explicitně, tato reprezentace diferenciálním operátorem dává konkrétní vzorec pro koeficienty , který lze použít pro rychlý výpočet těchto polynomů.

Protože formální výraz pro Weierstrassovu transformaci je , vidíme, že Weierstrassova transformace je . Weierstrassova transformace tedy v zásadě převádí řadu Hermitových polynomů na odpovídající Taylorovu řadu.

Existence nějaké formální mocninné řady s nenulovým konstantním koeficientem, takové, že , je dalším ekvivalentem tvrzení, že tyto polynomy tvoří Appellovu posloupnost. Protože jsou Appellovou posloupností, jsou také Shefferovou posloupností.

Reprezentace křivkovým integrálem

Z reprezentace generující funkce uvedené výše, vidíme, že Hermitovy polynomy lze reprezentovat pomocí křivkový integrál, protože :\begin{align} \mathit{He}_n(x) &= \frac{n. }{2\pi i} \oint_C \frac{e^{tx-\frac{t^2}{2}}}{t^{n+1}}\,dt, \\ H_n(x) &= \frac{n. +more}{2\pi i} \oint_C \frac{e^{2tx-t^2}}{t^{n+1}}\,dt, \end{align} s křivkou obkružující počátek souřadnicového systému.

Zobecnění

Pravděpodobnostní Hermitovy polynomy definované výše jsou ortogonální vůči standardnímu normálnímu rozdělení pravděpodobnosti, jehož hustota je : \frac{1}{\sqrt{2\pi}} e^{-\frac{x^2}{2}}, a které má střední hodnotu 0 a rozptyl 1.

Při změně měřítka je možné obdobně mluvit o zobecněných Hermitových polynomech : \mathit{He}_n^{[\alpha]}(x) s rozptylem , kde je jakékoli kladné číslo. Tyto jsou pak ortogonální vzhledem k normální rozdělení pravděpodobnosti, jejíž hustota je : (2\pi\alpha)^{-\frac12} e^{-\frac{x^2}{2\alpha}}. +more Jsou daný : \mathit{He}_n^{[\alpha]}(x) = \alpha^{\frac{n}{2}}\mathit{He}_n\left(\frac{x}{\sqrt{\alpha}}\right) = \left(\frac{\alpha}{2}\right)^{\frac{n}{2}} H_n\left( \frac{x}{\sqrt{2 \alpha}}\right) = e^{-\frac{\alpha D^2}{2}} \left(x^n\right).

Nyní, pokud : \mathit{He}_n^{[\alpha]}(x) = \sum_{k=0}^n h^{[\alpha]}_{n,k} x^k, pak posloupnost polynomů, jejíž -tý člen je : \left(\mathit{He}_n^{[\alpha]} \circ \mathit{He}^{[\beta]}\right)(x) \equiv \sum_{k=0}^n h^{[\alpha]}_{n,k}\,\mathit{He}_k^{[\beta]}(x) se nazývá stínová kompozice dvou posloupností polynomů. Lze dokázat, že vyhovuje identitám : \left(\mathit{He}_n^{[\alpha]} \circ \mathit{He}^{[\beta]}\right)(x) = \mathit{He}_n^{[\alpha+\beta]}(x) a : \mathit{He}_n^{[\alpha+\beta]}(x + y) = \sum_{k=0}^n \binom{n}{k} \mathit{He}_k^{[\alpha]}(x) \mathit{He}_{n-k}^{[\beta]}(y). +more Poslední vztah lze vyjádřit tím, že řekneme, že tato parametrizovaná rodina posloupností polynomů je známá jako křížová posloupnost. (Viz výše uvedená část o Appellových posloupnostech a o #Reprezentace diferenciálním operátorem|reprezentaci diferenciálním operátorem, která vede k její připravené derivaci. Tento vztah identity binomického typu pro jsme již zaznamenali ve výše uvedené části o #Rekurentní vzorec|rekurentních vzorcích. ).

„Záporný rozptyl“

Protože posloupnosti polynomů tvoří grupu s operací stínové kompozice, je možné pomocí : \mathit{He}_n^{[-\alpha]}(x) zapsat posloupnost, která je inverzní k podobně označené posloupnosti, ale bez znaménka minus, proto mluvíme o Hermitových polynomech se záporným rozptylem. Pro jsou koeficienty \mathit{He}_n^{[-\alpha]}(x) pouze absolutními hodnotami odpovídajících koeficientů \mathit{He}_n^{[\alpha]}(x).

Tyto koeficienty se objevují jako momenty normálního rozdělení pravděpodobnosti: -tý moment normálního rozdělení se střední hodnotou a rozptylem je : E[X^n] = \mathit{He}_n^{[-\sigma^2]}(\mu), kde je náhodná proměnná s uvedeným normálním rozdělením. Speciální případ křížové posloupnosti identit pak říká, že : \sum_{k=0}^n \binom{n}{k} \mathit{He}_k^{[\alpha]}(x) \mathit{He}_{n-k}^{[-\alpha]}(y) = \mathit{He}_n^{[0]}(x + y) = (x + y)^n.

Aplikace

Hermitovy funkce

Z fyzikálních polynomů je možnédefinovat Hermitovy funkce (často nazývané Hermitovy-gaussovy funkce): : \psi_n(x) = \left (2^n n. \sqrt{\pi} \right )^{-\frac12} e^{-\frac{x^2}{2}} H_n(x) = (-1)^n \left (2^n n. +more \sqrt{\pi} \right)^{-\frac12} e^{\frac{x^2}{2}} \frac{d^n}{\mbox{d}x^n} e^{-x^2}. tedy : \sqrt{2(n+1)}~~\psi_{n+1}(x)= \left ( x- {d\over \mbox{d}x}\right ) \psi_n(x).

Protože tyto funkce obsahují druhou odmocninu váhové funkce a jejich měřítko bylo vhodným způsobem upravené, jsou ortonormální: : \int_{-\infty}^\infty \psi_n(x) \psi_m(x) \,\mbox{d}x = \delta_{nm}, a tvoří ortonormální bázi prostoru . Tato skutečnost je ekvivalentní se stejným tvrzením pro Hermitovy polynomy (viz výše).

Hermitovy funkce úzce souvisí s Whittakerovou funkcí : : D_n(z) = \left(n! \sqrt{\pi}\right)^{\frac12} \psi_n\left(\frac{z}{\sqrt 2}\right) = (-1)^n e^\frac{z^2}{4} \frac{d^n}{dz^n} e^\frac{-z^2}{2} a díky tomu i s dalšími parabolickými válcovými funkcemi.

Hermitovy funkce vyhovují diferenciální rovnici : \psi_n(x) + \left(2n + 1 - x^2\right) \psi_n(x) = 0. Tato rovnice je ekvivalentní se Schrödingerovou rovnicí pro harmonický oscilátor v kvantové mechanice, a tyto funkce jsou tedy vlastní funkce.

Hermitovy funkce: 0 (plná modrá), 1 (čárkovaná oranžová), 2 (čerchovaná zelená), 3 (tečkovaná červená), 4 (plná fialová), a 5 (čárkovaná hnědá) :\begin{align} \psi_0(x) &= \pi^{-\frac14} \, e^{-\frac12 x^2}, \\ \psi_1(x) &= \sqrt{2} \, \pi^{-\frac14} \, x \, e^{-\frac12 x^2}, \\ \psi_2(x) &= \left(\sqrt{2} \, \pi^{\frac14}\right)^{-1} \, \left(2x^2-1\right) \, e^{-\frac12 x^2}, \\ \psi_3(x) &= \left(\sqrt{3} \, \pi^{\frac14}\right)^{-1} \, \left(2x^3-3x\right) \, e^{-\frac12 x^2}, \\ \psi_4(x) &= \left(2 \sqrt{6} \, \pi^{\frac14}\right)^{-1} \, \left(4x^4-12x^2+3\right) \, e^{-\frac12 x^2}, \\ \psi_5(x) &= \left(2 \sqrt{15} \, \pi^{\frac14}\right)^{-1} \, \left(4x^5-20x^3+15x\right) \, e^{-\frac12 x^2}. +more \end{align} Hermitovy funkce: 0 (plná modrá), 2 (čárkovaná oranžová), 4 (čerchovaná zelená) 50 (plná červená).

Rekurentní vzorec

Podle rekurentního vzorce pro Hermitovy polynomy platí pro Hermitovy funkce : \psi_n'(x) = \sqrt{\frac{n}{2}}\,\psi_{n-1}(x) - \sqrt{\frac{n+1}{2}}\psi_{n+1}(x) a : x\psi_n(x) = \sqrt{\frac{n}{2}}\,\psi_{n-1}(x) + \sqrt{\frac{n+1}{2}}\psi_{n+1}(x).

Rozvoj prvního vzorce pro libovolnou -tou derivaci pro jakékoli kladné celé číslo vede k : \psi_n^{(m)}(x) = \sum_{k=0}^m \binom{m}{k} (-1)^k 2^\frac{m-k}{2} \sqrt{\frac{n. }{(n-m+k). +more}} \psi_{n-m+k}(x) \mathit{He}_k(x).

Tento vzorec může být používán ve spojení s rekurentními vzorci pro a pro efektivní výpočet jakékoli derivace Hermitovy funkce.

Cramérova nerovnost

Pro reálné vyhovují Hermitovy funkce následujícímu omezení, které dokázal Harald Cramér a Jack Indritz: : \bigl|\psi_n(x)\bigr| \le \pi^{-\frac14}.

Hermitovy funkce jako vlastní funkce Fourierovy transformace

Hermitovy funkce jsou sadou vlastních funkcí Fourierovy transformace \mathcal{F}. Pro ověření použijeme fyzikální verzi vytvořující funkce a znásobíme ji . +more Tím dostaneme : e^{-\frac12 x^2 + 2xt - t^2} = \sum_{n=0}^\infty e^{-\frac12 x^2} H_n(x) \frac{t^n}{n. }.

Fourierovu transformaci levé strany popisuje vzorec : \begin{align} \mathcal{F} \left\{ e^{ -\frac12 x^2 + 2xt - t^2 } \right\}(k) &= \frac{1}{\sqrt{2 \pi}}\int_{-\infty}^\infty e^{-ixk}e^{-\frac12 x^2 + 2xt - t^2}\, \mbox{d}x \\ &= e^{-\frac12 k^2 - 2kit + t^2 } \\ &= \sum_{n=0}^\infty e^{ -\frac12 k^2 } H_n(k) \frac{(-it)^n}{n. }. +more \end{align}.

Fourierovu transformaci pravé strany pak vzorec : \mathcal{F} \left\{ \sum_{n=0}^\infty e^{-\frac12 x^2} H_n(x) \frac {t^n}{n. } \right\} = \sum_{n=0}^\infty \mathcal{F} \left \{ e^{-\frac12 x^2} H_n(x) \right\} \frac{t^n}{n. +more}.

Srovnáním stejných mocnin v transformované verzi levé a pravé strany dostáváme : \mathcal{F} \left\{ e^{-\frac12 x^2} H_n(x) \right\} = (-i)^n e^{-\frac12 k^2} H_n(k).

Hermitovy funkce jsou tedy ortonormální bází prostoru , která diagonalizuje Fourierův transformační operátor. V tomto případě byla použita unitární verze Fourierovy transformace, takže vlastní čísla jsou . +more Další pak slouží k definici mocnin (včetně racionálních) Fourierovy transformace pro nápadité získání zobecnění zlomkové Fourierovy transformace, resp. Mehlerova jádra.

Wignerova distribuce Hermitovy funkce

Wignerova distribuční funkce -tého řádu Hermitovy funkce souvisí s Laguerrovými polynomy -tého řádu. Laguerrovy polynomy jsou : L_n(x) := \sum_{k=0}^n \binom{n}{k} \frac{(-1)^k}{k. +more}x^k, které vedou k oscilátorovým Laguerrovým funkcím :l_n (x) := e^{-\frac{x}{2}} L_n(x). Pro všechna přirozená čísla je zřejmé, že :W_{\psi_n}(t,f) = (-1)^n l_n \big(4\pi (t^2 + f^2) \big), kde Wignerova funkce rozdělení je definována jako : W_x(t,f) = \int_{-\infty}^\infty x\left(t + \frac{\tau}{2}\right) \, x\left(t - \frac{\tau}{2}\right)^* \, e^{-2\pi i\tau f} \,d\tau. To je základní výsledek pro kvantový harmonický oscilátor, který v roce 1946 objevil Hilbrand J. Groenewold a publikoval ve své disertační práci. Jedná se o standardní paradigma kvantové mechaniky ve fázovém prostoru.

Mezi těmito dvěma rodinami polynomů existují další vztahy.

Kombinatorická interpretace koeficientů

V Hermitově polynomu s rozptylem 1 je absolutní hodnota koeficientu rovna počtu (neuspořádaných) dělení -prvkové množiny na singletonů a (neuspořádaných) dvojic. Ekvivalentně je to počet involucí -prvkové množiny s právě pevnými body, což je počet párování v úplném grafu s vrcholy, které ponechávají vrcholů nepokrytých (skutečně, Hermitovy polynomy jsou polynomy párování těchto grafů). +more Součet absolutních hodnot koeficientů dává celkový počet dělení na singletony a dvojice, tak zvaná telefonní čísla : 1, 1, 2, 4, 10, 26, 76, 232, 764, 2620, 9496,. .

Tato kombinatorická interpretace je příbuzná s kompletními exponenciálními Bellovými polynomy jako : \mathit{He}_n(x) = B_n(x, -1, 0, \ldots, 0), kde pro všechna .

Tato čísla je možné také vyjádřit jako speciální hodnotu Hermitových polynomů: : T(n) = \frac{\mathit{He}_n(i)}{i^n}.

Relace úplnosti

Christoffelův-Darbouxův vzorec pro Hermitovy polynomy má tvar :\sum_{k=0}^n \frac{H_k(x) H_k(y)}{k. 2^k} = \frac{1}{n. +more2^{n+1}}\,\frac{H_n(y) H_{n+1}(x) - H_n(x) H_{n+1}(y)}{x - y}.

Navíc pro výše uvedené Hermitovy funkce platí následující identita úplnosti ve smyslu distribucí: : \sum_{n=0}^\infty \psi_n(x) \psi_n(y) = \delta(x - y), kde je Diracovo delta, jsou Hermitovy funkce a reprezentuje Lebesgueovu míru na přímce v normalizovanou tak, že její projekce na horizontální osu je obvyklá Lebesgueova míra.

Použitím limity v Mehlerově vzorci, který platí pro , z této distribuční identity podle plyne : E(x, y; u) := \sum_{n=0}^\infty u^n \, \psi_n (x) \, \psi_n (y) = \frac{1}{\sqrt{\pi (1 - u^2)}} \, \exp\left(-\frac{1 - u}{1 + u} \, \frac{(x + y)^2}{4} - \frac{1 + u}{1 - u} \, \frac{(x - y)^2}{4}\right), což bývá často uváděno ekvivalentně jako separabilní jádro, : \sum_{n=0}^\infty \frac{H_n(x) H_n(y)}{n!} \left(\frac u 2\right)^n = \frac{1}{\sqrt{1 - u^2}} e^{\frac{2u}{1 + u}xy - \frac{u^2}{1 - u^2}(x - y)^2}.

Funkce je gaussovská hustota pravděpodobnosti funkce dvou proměnných na , která je při přiblížení k 1, velmi zahuštěná kolem přímky a velmi roztažená dále od této přímky. Odtud plyne, že : \sum_{n=0}^\infty u^n \langle f, \psi_n \rangle \langle \psi_n, g \rangle = \iint E(x, y; u) f(x) \overline{g(y)} \,\mbox{d}x \,\mbox{d}y \to \int f(x) \overline{g(x)} \,\mbox{d}x = \langle f, g \rangle pokud jsou funkce a spojité a mají kompaktní support.

Odtud je možné vyjádřit v Hermitově funkci jako sumu řady vektorů v , jmenovitě, : f = \sum_{n=0}^\infty \langle f, \psi_n \rangle \psi_n.

Pro důkaz této rovnosti pro použijeme opakovaně Fourierovu transformaci Gaussových funkcí: : \rho \sqrt{\pi} e^{-\frac{\rho^2 x^2}{4}} = \int e^{isx - \frac{s^2}{\rho^2}} \,ds \quad \text{for }\rho > 0.

Hermitův polynom je pak reprezentován jako : H_n(x) = (-1)^n e^{x^2} \frac {\mbox{d}^n}{\mbox{d}x^n} \left( \frac {1}{2\sqrt{\pi}} \int e^{isx - \frac{s^2}{4}} \,ds \right) = (-1)^n e^{x^2}\frac{1}{2\sqrt{\pi}} \int (is)^n e^{isx - \frac{s^2}{4}} \,ds.

Z této reprezentace a je zřejmé, že : \begin{align} E(x, y; u) &= \sum_{n=0}^\infty \frac{u^n}{2^n n. \sqrt{\pi}} \, H_n(x) H_n(y) e^{-\frac{x^2+y^2}{2}} \\ &= \frac{e^{\frac{x^2+y^2}{2}}}{4\pi\sqrt{\pi}}\iint\left( \sum_{n=0}^\infty \frac{1}{2^n n. +more} (-ust)^n \right ) e^{isx+ity - \frac{s^2}{4} - \frac{t^2}{4}}\, ds\,dt \\ & =\frac{e^{\frac{x^2+y^2}{2}}}{4\pi\sqrt{\pi}}\iint e^{-\frac{ust}{2}} \, e^{isx+ity - \frac{s^2}{4} - \frac{t^2}{4}}\, ds\,dt, \end{align} což, opět pomocí Fourierovy transformace gaussovských jader při substituci, dává požadovanou identitu : s = \frac{\sigma + \tau}{\sqrt 2}, \quad t = \frac{\sigma - \tau}{\sqrt 2}.

Odkazy

Poznámky

Reference

Literatura

[url=https://gallica. bnf. +morefr/ark:/12148/bpt6k77600r/f362. image. r=oeuvres%20completes%20de%20laplace]Oeuvres complètes 12, pp. 357-412[/url], [url=http://cerebro. xu. edu/math/Sources/Laplace/defint. pdf]anglický překlad/translace[/url] . * - 2000 referencí of Bibliography na Hermitovy polynomy.

Související články

Hermitova transformace * Legendrovy polynomy * Mehlerovo jádro * Parabolická válcová funkce * Romanovského polynomy * Turánovy nerovnosti

Externí odkazy

[url=https://www. gnu. +moreorg/software/gsl/]GNU Scientific Library[/url] - knihovna pro jazyk C implementující Hermitovy polynomy, funkce, jejich derivace a kořeny (viz také GNU Scientific Library).

Kategorie:Ortogonální polynomy Kategorie:Polynomy Kategorie:Speciální hypergeometrické funkce

5 min read
Share this post:
Like it 8

Leave a Comment

Please, enter your name.
Please, provide a valid email address.
Please, enter your comment.
Enjoy this post? Join Cesko.wiki
Don’t forget to share it
Top