[go: up one dir, main page]

18 Orthogonal PolynomialsClassical Orthogonal Polynomials

§18.15 Asymptotic Approximations

Contents
  1. §18.15(i) Jacobi
  2. §18.15(ii) Ultraspherical
  3. §18.15(iii) Legendre
  4. §18.15(iv) Laguerre
  5. §18.15(v) Hermite
  6. §18.15(vi) Other Approximations

§18.15(i) Jacobi

With the exception of the penultimate paragraph, we assume throughout this subsection that α, β, and M (=0,1,2,…) are all fixed.

18.15.1 (sin⁡12⁢θ)α+12⁢(cos⁡12⁢θ)β+12⁢Pn(α,β)⁡(cos⁡θ)=π−1⁢22⁢n+α+β+1⁢B⁡(n+α+1,n+β+1)⁢(∑m=0M−1fm⁡(θ)2m⁢(2⁢n+α+β+2)m+O⁡(n−M)),

as n→∞, uniformly with respect to θ∈[δ,π−δ]. Here, and elsewhere in §18.15, δ is an arbitrary small positive constant. Also, B⁡(a,b) is the beta function (§5.12) and

18.15.2 fm⁡(θ)=∑ℓ=0mCm,ℓ⁡(α,β)ℓ!⁢(m−ℓ)!⁢cos⁡θn,m,ℓ(sin⁡12⁢θ)ℓ⁢(cos⁡12⁢θ)m−ℓ,

where

18.15.3 Cm,ℓ⁡(α,β)=(12+α)ℓ⁢(12−α)ℓ⁢(12+β)m−ℓ⁢(12−β)m−ℓ,

and

18.15.4 θn,m,ℓ=12⁢(2⁢n+α+β+m+1)⁢θ−12⁢(α+ℓ+12)⁢π.

When α,β∈(−12,12), the error term in (18.15.1) is less than twice the first neglected term in absolute value, in which one has to take cos⁡θn,m,ℓ=1. See Hahn (1980), where corresponding results are given when x is replaced by a complex variable z that is bounded away from the orthogonality interval [−1,1].

The case M=1 of (18.15.1) goes back to Darboux. It reads:

18.15.4_5 (sin⁡12⁢θ)α+12⁢(cos⁡12⁢θ)β+12⁢Pn(α,β)⁡(cos⁡θ)=π−12⁢n−12⁢cos⁡(12⁢(2⁢n+α+β+1)⁢θ−14⁢(2⁢α+1)⁢π)+O⁡(n−32),
α,β∈ℝ,

as n→∞, uniformly with respect to θ∈[δ,π−δ].

Next, let

18.15.5 ρ=n+12⁢(α+β+1).

Then as n→∞,

18.15.6 (sin⁡12⁢θ)α+12⁢(cos⁡12⁢θ)β+12⁢Pn(α,β)⁡(cos⁡θ)=Γ⁡(n+α+1)212⁢ρα⁢n!⁢(θ12⁢Jα⁡(ρ⁢θ)⁢∑m=0MAm⁡(θ)ρ2⁢m+θ32⁢Jα+1⁡(ρ⁢θ)⁢∑m=0M−1Bm⁡(θ)ρ2⁢m+1+εM⁡(ρ,θ)),

where Jν⁡(z) is the Bessel function (§10.2(ii)), and

18.15.7 εM⁡(ρ,θ)={θ⁢O⁡(ρ−2⁢M−(3/2)),c⁢ρ−1≤θ≤π−δ,θα+(5/2)⁢O⁡(ρ−2⁢M+α),0≤θ≤c⁢ρ−1,

with c denoting an arbitrary positive constant. Also,

18.15.8 A0⁡(θ) =1,
θ⁢B0⁡(θ) =14⁢g⁡(θ),
A1⁡(θ) =18⁢g′⁡(θ)−1+2⁢α8⁢g⁡(θ)θ−132⁢(g⁡(θ))2,

where

18.15.9 g⁡(θ)=(14−α2)⁢(cot⁡(12⁢θ)−(12⁢θ)−1)−(14−β2)⁢tan⁡(12⁢θ).

For higher coefficients see Baratella and Gatteschi (1988), and for another estimate of the error term in a related expansion see Wong and Zhao (2003). For large β, fixed α, and 0≤n/β≤c, Dunster (1999) gives asymptotic expansions of Pn(α,β)⁡(z) that are uniform in unbounded complex z-domains containing z=±1. These expansions are in terms of Whittaker functions (§13.14). This reference also supplies asymptotic expansions of Pn(α,β)⁡(z) for large n, fixed α, and 0≤β/n≤c. The latter expansions are in terms of Bessel functions, and are uniform in complex z-domains not containing neighborhoods of 1. For a complementary result, see Wong and Zhao (2004). By using the symmetry property given in the second row of Table 18.6.1, the roles of α and β can be interchanged.

For an asymptotic expansion of Pn(α,β)⁡(z) as n→∞ that holds uniformly for complex z bounded away from [−1,1], see Elliott (1971). The first term of this expansion also appears in Szegő (1975, Theorem 8.21.7).

§18.15(ii) Ultraspherical

For fixed λ∈(0,1) and fixed M=0,1,2,…,

18.15.10 Cn(λ)⁡(cos⁡θ)=22⁢λ⁢Γ⁡(λ+12)π12⁢Γ⁡(λ+1)⁢(2⁢λ)n(λ+1)n⁢(∑m=0M−1(λ)m⁢(1−λ)mm!⁢(n+λ+1)m⁢cos⁡θn,m(2⁢sin⁡θ)m+λ+O⁡(1nM)),

as n→∞ uniformly with respect to θ∈[δ,π−δ], where

18.15.11 θn,m=(n+m+λ)⁢θ−12⁢(m+λ)⁢π.

For a bound on the error term in (18.15.10) see Szegő (1975, Theorem 8.21.11).

Asymptotic expansions for Cn(λ)⁡(cos⁡θ) can be obtained from the results given in §18.15(i) by setting α=β=λ−12 and referring to (18.7.1). See also Szegő (1933) and Szegő (1975, Eq. (8.21.14)).

§18.15(iii) Legendre

For fixed M=0,1,2,…,

18.15.12 Pn⁡(cos⁡θ)=(2sin⁡θ)12⁢∑m=0M−1(−12m)⁢(m−12n)⁢cos⁡αn,m(2⁢sin⁡θ)m+O⁡(1nM+12),

as n→∞, uniformly with respect to θ∈[δ,π−δ], where

18.15.13 αn,m=(n−m+12)⁢θ+(n−12⁢m−14)⁢π.

Also, when 16⁢π<θ<56⁢π, the right-hand side of (18.15.12) with M=∞ converges; paradoxically, however, the sum is 2⁢Pn⁡(cos⁡θ) and not Pn⁡(cos⁡θ) as stated erroneously in Szegő (1975, §8.4(3)).

For these results and further information see Olver (1997b, pp. 311–313). Another expansion follows from (18.15.10) by taking λ=12; see Szegő (1975, Theorem 8.21.5).

For asymptotic expansions of Pn⁡(cos⁡θ) and Pn⁡(cosh⁡ξ) that are uniformly valid when 0≤θ≤π−δ and 0≤ξ<∞ see §14.15(iii) with μ=0 and ν=n. These expansions are in terms of Bessel functions and modified Bessel functions, respectively.

§18.15(iv) Laguerre

In Terms of Elementary Functions

For fixed M=0,1,2,…, and fixed α,

18.15.14 Ln(α)⁡(x)=n12⁢α−14⁢e12⁢xπ12⁢x12⁢α+14⁢(cos⁡θn(α)⁡(x)⁢(∑m=0M−1am⁡(x)n12⁢m+O⁡(1n12⁢M))+sin⁡θn(α)⁡(x)⁢(∑m=1M−1bm⁡(x)n12⁢m+O⁡(1n12⁢M))),

as n→∞, uniformly on compact x-intervals in (0,∞), where

18.15.15 θn(α)⁡(x)=2⁢(n⁢x)12−(12⁢α+14)⁢π.

The leading coefficients are given by

18.15.16 a0⁡(x) =1,
a1⁡(x) =0,
b1⁡(x) =148⁢x12⁢(4⁢x2−12⁢α2−24⁢α⁢x−24⁢x+3).

See also Deaño et al. (2013).

In Terms of Bessel Functions

Define

18.15.17 ν=4⁢n+2⁢α+2,
18.15.18 ξ=12⁢(x−x2+arcsin⁡(x)),
0≤x≤1.

Then for fixed M=0,1,2,…, and fixed α,

18.15.19 Ln(α)⁡(ν⁢x)=e12⁢ν⁢x2α⁢x12⁢α+14⁢(1−x)14⁢(ξ12⁢Jα⁡(ν⁢ξ)⁢∑m=0M−1Am⁡(ξ)ν2⁢m+ξ−12⁢Jα+1⁡(ν⁢ξ)⁢∑m=0M−1Bm⁡(ξ)ν2⁢m+1+ξ12⁢env⁡Jα⁡(ν⁢ξ)⁢O⁡(1ν2⁢M−1)),

as n→∞ uniformly for 0≤x≤1−δ. Here Jν⁡(z) denotes the Bessel function (§10.2(ii)), env⁡Jν⁡(z) denotes its envelope (§2.8(iv)), and δ is again an arbitrary small positive constant. The leading coefficients are given by A0⁡(ξ)=1 and

18.15.20 B0⁡(ξ)=−12⁢(1−4⁢α28+ξ⁢(1−xx)12⁢(4⁢α2−18+14⁢x1−x+524⁢(x1−x)2)).

In Terms of Airy Functions

Again define ν as in (18.15.17); also,

18.15.21 ζ =−(34⁢(arccos⁡(x)−x−x2))23,
0≤x≤1,
ζ =(34⁢(x2−x−arccosh⁡(x)))23,
x≥1.

Then for fixed M=0,1,2,…, and fixed α,

18.15.22 Ln(α)⁡(ν⁢x)=(−1)n⁢e12⁢ν⁢x2α−12⁢x12⁢α+14⁢(ζx−1)14×(Ai⁡(ν23⁢ζ)ν13⁢∑m=0M−1Em⁡(ζ)ν2⁢m+Ai′⁡(ν23⁢ζ)ν53⁢∑m=0M−1Fm⁡(ζ)ν2⁢m+envAi⁡(ν23⁢ζ)⁢O⁡(1ν2⁢M−23)),

as n→∞ uniformly for δ≤x<∞. Here Ai denotes the Airy function (§9.2), Ai′ denotes its derivative, and envAi denotes its envelope (§2.8(iii)). The leading coefficients are given by E0⁡(ζ)=1 and

18.15.23 F0⁡(ζ)=−548⁢ζ2+(x−1x⁢ζ)12⁢(12⁢α2−18−14⁢xx−1+524⁢(xx−1)2),
0≤x<∞.

§18.15(v) Hermite

Define

18.15.24 μ=2⁢n+1,
18.15.25 λn={Γ⁡(n+1)/Γ⁡(12⁢n+1),n⁢ even,Γ⁡(n+2)/(μ12⁢Γ⁡(12⁢n+32)),n⁢ odd,

and

18.15.26 ωn,m⁡(x)=μ12⁢x−12⁢(m+n)⁢π.

Then for fixed M=0,1,2,…,

18.15.27 Hn⁡(x)=λn⁢e12⁢x2⁢(∑m=0M−1um⁡(x)⁢cos⁡ωn,m⁡(x)μ12⁢m+O⁡(1μ12⁢M)),

as n→∞, uniformly on compact x-intervals on ℝ. The coefficients um⁡(x) are polynomials in x, and u0⁡(x)=1, u1⁡(x)=16⁢x3.

For more powerful asymptotic expansions as n→∞ in terms of elementary functions that apply uniformly when 1+δ≤t<∞, −1+δ≤t≤1−δ, or −∞<t≤−1−δ, where t=x/2⁢n+1 and δ is again an arbitrary small positive constant, see §§12.10(i)–12.10(iv) and 12.10(vi). And for asymptotic expansions as n→∞ in terms of Airy functions that apply uniformly when −1+δ≤t<∞ or −∞<t≤1−δ, see §§12.10(vii) and 12.10(viii). With μ=2⁢n+1 the expansions in Chapter 12 are for the parabolic cylinder function U⁡(−12⁢μ2,μ⁢t⁢2), which is related to the Hermite polynomials via

18.15.28 Hn⁡(x)=214⁢(μ2−1)⁢e12⁢μ2⁢t2⁢U⁡(−12⁢μ2,μ⁢t⁢2);

compare (18.11.3).

For an error bound for the first term in the Airy-function expansions see Olver (1997b, p. 403).

See also Geronimo et al. (2004).

§18.15(vi) Other Approximations

The asymptotic behavior of the classical OP’s as x→±∞ with the degree and parameters fixed is evident from their explicit polynomial forms; see, for example, (18.2.7) and the last two columns of Table 18.3.1.

For asymptotic approximations of Jacobi, ultraspherical, and Laguerre polynomials in terms of Hermite polynomials, see López and Temme (1999a). These approximations apply when the parameters are large, namely α and β (subject to restrictions) in the case of Jacobi polynomials, λ in the case of ultraspherical polynomials, and |α|+|x| in the case of Laguerre polynomials. See also Dunster (1999), Atia et al. (2014) and Temme (2015, Chapter 32).