[go: up one dir, main page]

10 Bessel FunctionsModified Bessel Functions

§10.41 Asymptotic Expansions for Large Order

Contents
  1. §10.41(i) Asymptotic Forms
  2. §10.41(ii) Uniform Expansions for Real Variable
  3. §10.41(iii) Uniform Expansions for Complex Variable
  4. §10.41(iv) Double Asymptotic Properties
  5. §10.41(v) Double Asymptotic Properties (Continued)

§10.41(i) Asymptotic Forms

§10.41(ii) Uniform Expansions for Real Variable

As ν→∞ through positive real values,

10.41.3 Iν⁡(ν⁢z) ∼eν⁢η(2⁢π⁢ν)12⁢(1+z2)14⁢∑k=0∞Uk⁡(p)νk,
10.41.4 Kν⁡(ν⁢z) ∼(π2⁢ν)12⁢e−ν⁢η(1+z2)14⁢∑k=0∞(−1)k⁢Uk⁡(p)νk,
10.41.5 Iν′⁡(ν⁢z) ∼(1+z2)14⁢eν⁢η(2⁢π⁢ν)12⁢z⁢∑k=0∞Vk⁡(p)νk,

uniformly for 0<z<∞. Here

10.41.7 η=(1+z2)12+ln⁡z1+(1+z2)12,
10.41.8 p=(1+z2)−12,

where the branches assume their principal values. Also, Uk⁡(p) and Vk⁡(p) are polynomials in p of degree 3⁢k, given by U0⁡(p)=V0⁡(p)=1, and

10.41.9 Uk+1⁡(p) =12⁢p2⁢(1−p2)⁢Uk′⁡(p)+18⁢∫0p(1−5⁢t2)⁢Uk⁡(t)⁢dt,
Vk+1⁡(p) =Uk+1⁡(p)−12⁢p⁢(1−p2)⁢Uk⁡(p)−p2⁢(1−p2)⁢Uk′⁡(p),
k=0,1,2,….

For k=1,2,3,

10.41.10 U1⁡(p) =124⁢(3⁢p−5⁢p3),
U2⁡(p) =11152⁢(81⁢p2−462⁢p4+385⁢p6),
U3⁡(p) =14 14720⁢(30375⁢p3−3 69603⁢p5+7 65765⁢p7−4 25425⁢p9),
10.41.11 V1⁡(p) =124⁢(−9⁢p+7⁢p3),
V2⁡(p) =11152⁢(−135⁢p2+594⁢p4−455⁢p6),
V3⁡(p) =14 14720⁢(−42525⁢p3+4 51737⁢p5−8 83575⁢p7+4 75475⁢p9).

For U4⁡(p), U5⁡(p), U6⁡(p), see Bickley et al. (1952, p. xxxv).

For numerical tables of η=η⁡(z) and the coefficients Uk⁡(p), Vk⁡(p), see Olver (1962, pp. 43–51).

§10.41(iii) Uniform Expansions for Complex Variable

The expansions (10.41.3)–(10.41.6) also hold uniformly in the sector |ph⁡z|≤12⁢π−δ (<12⁢π), with the branches of the fractional powers in (10.41.3)–(10.41.8) extended by continuity from the positive real z-axis.

Figures 10.41.1 and 10.41.2 show corresponding points of the mapping of the z-plane and the η-plane. The curve E1⁢B⁢E2 in the z-plane is the upper boundary of the domain 𝐊 depicted in Figure 10.20.3 and rotated through an angle −12⁢π. Thus B is the point z=c, where c is given by (10.20.18).

For derivations of the results in this subsection, and also error bounds, see Olver (1997b, pp. 374–378). For extensions of the regions of validity in the z-plane and extensions to complex values of ν see Olver (1997b, pp. 378–382).

See accompanying text
Figure 10.41.1: z-plane. Magnify
See accompanying text
Figure 10.41.2: η-plane. Magnify

For expansions in inverse factorial series see Dunster et al. (1993).

§10.41(iv) Double Asymptotic Properties

The series (10.41.3)–(10.41.6) can also be regarded as generalized asymptotic expansions for large |z|. Thus as z→∞ with ℓ (≥1) and ν (>0) both fixed,

Similarly for (10.41.5) and (10.41.6).

In the case of (10.41.13) with positive real values of z the result is a consequence of the error bounds given in Olver (1997b, pp. 377–378). Then by expanding the quantities η, (1+z2)−14, and Uk⁡(p), k=0,1,…,ℓ−1, and rearranging, we arrive at an expansion of the right-hand side of (10.41.13) in powers of z−1. Moreover, because of the uniqueness property of asymptotic expansions (§2.1(iii)) this expansion must agree with (10.40.2), with z replaced by ν⁢z, up to and including the term in z−(ℓ−1). It also enjoys the same sector of validity.

To establish (10.41.12) we substitute into (10.34.3), with m=0 and z replaced by ν⁢z, by means of (10.41.13) observing that when |z| is large the effect of replacing z by z⁢e±π⁢i is to replace η, (1+z2)14, and p by −η, ±i⁢(1+z2)14, and −p, respectively.

§10.41(v) Double Asymptotic Properties (Continued)

Similar analysis can be developed for the uniform asymptotic expansions in terms of Airy functions given in §10.20. We first prove that for the expansions (10.20.6) for the Hankel functions Hν(1)⁡(ν⁢z) and Hν(2)⁡(ν⁢z) the z-asymptotic property applies when z→±i⁢∞, respectively. This is a consequence of the error bounds associated with these expansions. We then extend the validity of this property from z→±i⁢∞ to z→∞ in the sector −π+δ≤ph⁡z≤2⁢π−δ in the case of Hν(1)⁡(ν⁢z), and to z→∞ in the sector −2⁢π+δ≤ph⁡z≤π−δ in the case of Hν(2)⁡(ν⁢z). This is done by re-expansion with the aid of (10.20.10), (10.20.11), and §10.41(ii), followed by comparison with (10.17.5) and (10.17.6), with z replaced by ν⁢z. Lastly, we substitute into (10.4.4), again with z replaced by ν⁢z. The final results are:

10.41.14 Jν⁡(ν⁢z) =(4⁢ζ1−z2)14⁢(Ai⁡(ν23⁢ζ)ν13⁢(∑k=0ℓAk⁡(ζ)ν2⁢k+O⁡(1ζ3⁢ℓ+3))+Ai′⁡(ν23⁢ζ)ν53⁢(∑k=0ℓ−1Bk⁡(ζ)ν2⁢k+O⁡(1ζ3⁢ℓ+1))),
10.41.15 Yν⁡(ν⁢z) =−(4⁢ζ1−z2)14⁢(Bi⁡(ν23⁢ζ)ν13⁢(∑k=0ℓAk⁡(ζ)ν2⁢k+O⁡(1ζ3⁢ℓ+3))+Bi′⁡(ν23⁢ζ)ν53⁢(∑k=0ℓ−1Bk⁡(ζ)ν2⁢k+O⁡(1ζ3⁢ℓ+1))),

as z→∞ in |ph⁡z|≤π−δ, or equivalently as ζ→∞ in |ph⁡(−ζ)|≤23⁢π−δ, for fixed ℓ (≥0) and fixed ν (>0).

It needs to be noted that the results (10.41.14) and (10.41.15) do not apply when z→0+ or equivalently ζ→+∞. This is because Ak⁡(ζ) and ζ−12⁢Bk⁡(ζ),k=0,1,…, do not form an asymptotic scale (§2.1(v)) as ζ→+∞; see Olver (1997b, pp. 422–425).