Diskrete orthogonale Polynome

Diskrete orthogonale Polynome sind orthogonale Polynome bezüglich eines diskreten Maßes. Solche Polynome findet man unter anderem in der Stochastik und in der statistischen Physik, wo man mit diskreten Wahrscheinlichkeitsverteilungen zu tun hat.

Beispiele sind die Meixner-Polynome, die Krawtschuk-Polynome, die diskreten Tschebyscheff-Polynome, die Hahn-Polynome und die Charlier-Polynome.

Diskrete orthogonale Polynome

Konstruktion eines diskreten Maßes mit Gewichtsfunktion

Sei

• $S\in \mathbb{N}\cup \{\mathrm{\infty }\}$, eine natürliche Zahl oder Unendlich,
• $a_0,a_1,a_2,\dots$ eine positive Folge, das heißt $a_i>0\mathrm{\forall }i$,
• $s_0,s_1,s_2,\dots$ eine Folge reeller Zahlen, welche den Träger bilden werden,
• $𝒜$ eine σ-Algebra, welche mindestens alle Singletons $\{s_i\}$ enthält,
• ${\delta }_{s_i}$ das Diracmaß bezüglich $s_i$ auf der σ-Algebra $𝒜$, das bedeutet für eine Menge $A\in 𝒜$ ist dieses Maß wie folgt definiert ${\delta }_{s_i}(A):=1_A(s_i)=\{\begin{array}{cc}1,& \text{falls }{s}_i\in A,\\ 0& \text{sonst}\end{array}$.

Nun definieren wir ein diskretes Maß auf $𝒜$ mit Hilfe der Diracmaße und der beiden Folgen

$\mu =\sum _{n=0}^S{a}_n{\delta }_{s_n}.$

Weiter soll gelten, dass $\mu$ endliche Momente hat (das bedeutet wenn eine Zufallsvariable $X\sim \mu$, dann gilt $𝔼[|X{|}^n]<\mathrm{\infty }$ für alle $n=1,2,\dots$).

Für die Folge $(a_n)$ können wir nun eine Gewichtsfunktion $\omega :ℝ\to {ℝ}_{\ge 0}$ durch

$\omega (s_n)=a_n$

definieren.

Ohne Beschränkung der Allgemeinheit wählen wir nun als Träger $s_n:=n$ für alle $n=0,1,2,\dots ,S$, dann ist die Gewichtsfunktion durch

$\omega (n)=a_n$

und das diskrete Maß durch

$\mu =\sum _{n=0}^S{a}_n{\delta }_n.$

gegeben.

Diskrete orthogonale Polynome

Eine Familie von orthogonalen Polynome $(p_n(x){)}_{n\ge 0}$ heißt diskret, wenn sie orthogonal bezüglich eines diskreten Maßes $\mu$ mit Gewichtsfunktion $\omega (x)$ sind, das heißt wenn

$\sum _{x=0}^S{p}_n(x)p_m(x)\omega (x)={\kappa }_n{\delta }_{n,m}$

erfüllt ist, wobei ${\delta }_{n,m}$ das Kronecker-Delta bezeichnet. Die linke Seite dieser Formel beschreibt ein Skalarprodukt

$⟨p_n,p_m⟩:=\sum _{x=0}^S{p}_n(x)p_m(x)\omega (x)$

Beispiele

• Meixner-Polynome (Negative Binomialverteilung):

$M_n(x;\beta ,c)={}_2{F}_1\left(\begin{array}{c}-n,-x\\ \beta \end{array}\right|1-\frac{1}{c}),$ $S=\mathrm{\infty },$ $\omega (x)=\frac{(\beta {)}_x}{x!}{c}^x$ und ${\kappa }_n=\frac{n!(1-c{)}^{-\beta }}{c^n(\beta {)}_n},$

wobei die Orthogonalität nur für $\beta >0$ und $0<c<1$ gilt.

• Charlier-Polynome (Poisson-Verteilung):

$C_n(x;a):={}_2{F}_0\left(\begin{array}{c}-n,-x\\ -\end{array}\right|-\frac{1}{a}),$ $S=\mathrm{\infty },$ $\omega (x)=\frac{a^x}{x!}$ und ${\kappa }_n=\frac{n!}{a^n}{e}^a.$

Sonstiges

Sei $u(x)$ eine Funktion definiert durch die Beziehung

$\omega (x+1)-\omega (x)=-u(x+1)\omega (x+1).$

Betrachtet man allgemeine orthogonale Polynome mit einer Gewichtsfunktion $v(x)$ und sei $f(x)$ eine Funktion definiert durch die Beziehung

$v(x)=\exp (-f(x)),$

so entspricht $u(x)$ dem diskreten Pendant der Funktion $f(x)$ respektive $-\log (v(x))$.

Differenzengleichung

Es lässt sich beweisen, dass jedes diskrete orthogonale Polynom einer Differenzengleichung zweiter Ordnung genügt, wenn das Maß einen Träger über einer Halbgeraden mit äquidistanten Punkt besitzt (d. h. ein Gitter).

Annahmen

Sei $\{p_n(x)\}$ eine Familie orthogonaler Polynome bezüglich eines diskreten Maßes $\mu$ mit Träger

$𝒯=\{s,s+1,s+2,\dots ,t\}\subset ℝ,s\in ℝ,t\in ℝ\cup \{\mathrm{\infty }\}.$

Wir nehmen an, dass $p_n$ gerade vom Grad $n$ ist und die Gewichtsfunktion $\omega (l)$ normalisiert ist, d. h. es gilt

$\sum _{l=s}^t{p}_m(l)p_n(l)\omega (l)={\kappa }_m{\delta }_{m,n}$ und $\sum _{l=s}^t\omega (l)=1.$

Weiter nehmen wir an, dass auf $ℝ\setminus 𝒯$ die Gewichtsfunktion $\omega (l)$ nicht konstant $0$ ist, aber für die Randpunkt gilt $\omega (s-1)=0$ und $\omega (t+1)=0$.

Weiter notieren wir mit $\mathrm{\Delta }$ den Differenzoperator $\mathrm{\Delta }f:=f(x+1)-f(x).$ Die Funktion $u(x)$ haben wir im vorherigen Abschnitt definiert.

Aussage des Theorems

Sei

$p_n={\gamma }_n{x}^n+\dots$

ein diskretes orthogonales Polynom, welches die vorherigen Annahmen erfüllt. Dann gilt

$\mathrm{\Delta }{p}_n(x)=A_n(x)p_{n-1}(x)-B_n(x)p_n(x)$

wobei $A_n$ und $B_n$ wie folgt definiert sind

$A_n(x)=\frac{{\gamma }_{n-1}}{{\gamma }_n{\kappa }_{n-1}}\frac{p_n(t+1)p_n(t)}{(t-x)}\omega (t)+\frac{{\gamma }_{n-1}}{{\gamma }_n{\kappa }_{n-1}}\sum _{l=s}^t{p}_n(l)p_n(l-1)\frac{u(x+1)-u(l)}{(x+1-l)}\omega (l)$

und

$B_n(x)=\frac{{\gamma }_{n-1}}{{\gamma }_n{\kappa }_{n-1}}\frac{p_n(t+1)p_{n-1}(t)}{(t-x)}\omega (t)+\frac{{\gamma }_{n-1}}{{\gamma }_n{\kappa }_{n-1}}\sum _{l=s}^t{p}_n(l)p_{n-1}(l-1)\frac{u(x+1)-u(l)}{(x+1-l)}\omega (l).$

Literatur

• Vorlage:Literatur