Dual Hahn polynomials

From HandWiki

In mathematics, the dual Hahn polynomials are a family of orthogonal polynomials in the Askey scheme of hypergeometric orthogonal polynomials. They are defined on a non-uniform lattice x(s)=s(s+1) and are defined as

wn(c)(s,a,b)=(a−b+1)n(a+c+1)nn!3F2(−n,a−s,a+s+1;a−b+a,a+c+1;1)

for n=0,1,...,N−1 and the parameters a,b,c are restricted to −12<a<b,|c|<1+a,b=a+N.

Note that (u)k is the rising factorial, otherwise known as the Pochhammer symbol, and 3F2(⋅) is the generalized hypergeometric functions

Roelof Koekoek, Peter A. Lesky, and René F. Swarttouw (2010, 14) give a detailed list of their properties.

Orthogonality

The dual Hahn polynomials have the orthogonality condition

∑s=ab−1wn(c)(s,a,b)wm(c)(s,a,b)ρ(s)[Δx(s−12)]=δnmdn2

for n,m=0,1,...,N−1. Where Δx(s)=x(s+1)−x(s),

ρ(s)=Γ(a+s+1)Γ(c+s+1)Γ(s−a+1)Γ(b−s)Γ(b+s+1)Γ(s−c+1)

and

dn2=Γ(a+c+n+a)n!(b−a−n−1)!Γ(b−c−n).

Numerical instability

As the value of n increases, the values that the discrete polynomials obtain also increases. As a result, to obtain numerical stability in calculating the polynomials you would use the renormalized dual Hahn polynomial as defined as

w^n(c)(s,a,b)=wn(c)(s,a,b)ρ(s)dn2[Δx(s−12)]

for n=0,1,...,N−1.

Then the orthogonality condition becomes

∑s=ab−1w^n(c)(s,a,b)w^m(c)(s,a,b)=δm,n

for n,m=0,1,...,N−1

Relation to other polynomials

The Hahn polynomials, hn(x,N;α,β), is defined on the uniform lattice x(s)=s, and the parameters a,b,c are defined as a=(α+β)/2,b=a+N,c=(β−α)/2. Then setting α=β=0 the Hahn polynomials become the Chebyshev polynomials. Note that the dual Hahn polynomials have a q-analog with an extra parameter q known as the dual q-Hahn polynomials.

Racah polynomials are a generalization of dual Hahn polynomials.

References