Bessel polynomials

From The Right Wiki
(Redirected from Bessel polynomial)
Jump to navigationJump to search

In mathematics, the Bessel polynomials are an orthogonal sequence of polynomials. There are a number of different but closely related definitions. The definition favored by mathematicians is given by the series[1]: 101 

yn(x)=k=0n(n+k)!(nk)!k!(x2)k.

Another definition, favored by electrical engineers, is sometimes known as the reverse Bessel polynomials[2]: 8 [3]: 15 

θn(x)=xnyn(1/x)=k=0n(n+k)!(nk)!k!xnk2k.

The coefficients of the second definition are the same as the first but in reverse order. For example, the third-degree Bessel polynomial is

y3(x)=15x3+15x2+6x+1

while the third-degree reverse Bessel polynomial is

θ3(x)=x3+6x2+15x+15.

The reverse Bessel polynomial is used in the design of Bessel electronic filters.

Properties

Definition in terms of Bessel functions

The Bessel polynomial may also be defined using Bessel functions from which the polynomial draws its name.

yn(x)=xnθn(1/x)
yn(x)=2πxe1/xKn+12(1/x)
θn(x)=2πxn+1/2exKn+12(x)

where Kn(x) is a modified Bessel function of the second kind, yn(x) is the ordinary polynomial, and θn(x) is the reverse polynomial .[2]: 7, 34  For example:[4]

y3(x)=15x3+15x2+6x+1=2πxe1/xK3+12(1/x)

Definition as a hypergeometric function

The Bessel polynomial may also be defined as a confluent hypergeometric function[5]: 8 

yn(x)=2F0(n,n+1;;x/2)=(2x)nU(n,2n,2x)=(2x)n+1U(n+1,2n+2,2x).

A similar expression holds true for the generalized Bessel polynomials (see below):[2]: 35 

yn(x;a,b)=2F0(n,n+a1;;x/b)=(bx)n+a1U(n+a1,2n+a,bx).

The reverse Bessel polynomial may be defined as a generalized Laguerre polynomial:

θn(x)=n!(2)nLn2n1(2x)

from which it follows that it may also be defined as a hypergeometric function:

θn(x)=(2n)n(2)n1F1(n;2n;2x)

where (−2n)n is the Pochhammer symbol (rising factorial).

Generating function

The Bessel polynomials, with index shifted, have the generating function

n=02πxn+12exKn12(x)tnn!=1+xn=1θn1(x)tnn!=ex(112t).

Differentiating with respect to t, cancelling x, yields the generating function for the polynomials {θn}n0

n=0θn(x)tnn!=112tex(112t).

Similar generating function exists for the yn polynomials as well:[1]: 106 

n=0yn1(x)tnn!=exp(112xtx).

Upon setting t=zxz2/2, one has the following representation for the exponential function:[1]: 107 

ez=n=0yn1(x)(zxz2/2)nn!.

Recursion

The Bessel polynomial may also be defined by a recursion formula:

y0(x)=1
y1(x)=x+1
yn(x)=(2n1)xyn1(x)+yn2(x)

and

θ0(x)=1
θ1(x)=x+1
θn(x)=(2n1)θn1(x)+x2θn2(x)

Differential equation

The Bessel polynomial obeys the following differential equation:

x2d2yn(x)dx2+2(x+1)dyn(x)dxn(n+1)yn(x)=0

and

xd2θn(x)dx22(x+n)dθn(x)dx+2nθn(x)=0

Orthogonality

The Bessel polynomials are orthogonal with respect to the weight e2/x integrated over the unit circle of the complex plane.[1]: 104  In other words, if nm, 02πyn(eiθ)ym(eiθ)ieiθdθ=0

Generalization

Explicit form

A generalization of the Bessel polynomials have been suggested in literature, as following:

yn(x;α,β):=(1)nn!(xβ)nLn(12nα)(βx),

the corresponding reverse polynomials are

θn(x;α,β):=n!(β)nLn(12nα)(βx)=xnyn(1x;α,β).

The explicit coefficients of the yn(x;α,β) polynomials are:[1]: 108 

yn(x;α,β)=k=0n(nk)(n+k+α2)k_(xβ)k.

Consequently, the θn(x;α,β) polynomials can explicitly be written as follows:

θn(x;α,β)=k=0n(nk)(2nk+α2)nk_xkβnk.

For the weighting function

ρ(x;α,β):=1F1(1,α1,βx)

they are orthogonal, for the relation

0=cρ(x;α,β)yn(x;α,β)ym(x;α,β)dx

holds for mn and c a curve surrounding the 0 point. They specialize to the Bessel polynomials for α = β = 2, in which situation ρ(x) = exp(−2/x).

Rodrigues formula for Bessel polynomials

The Rodrigues formula for the Bessel polynomials as particular solutions of the above differential equation is :

Bn(α,β)(x)=an(α,β)xαeβx(ddx)n(xα+2neβx)

where a(α, β)
n
are normalization coefficients.

Associated Bessel polynomials

According to this generalization we have the following generalized differential equation for associated Bessel polynomials:

x2d2Bn,m(α,β)(x)dx2+[(α+2)x+β]dBn,m(α,β)(x)dx[n(α+n+1)+mβx]Bn,m(α,β)(x)=0

where 0mn. The solutions are,

Bn,m(α,β)(x)=an,m(α,β)xα+meβx(ddx)nm(xα+2neβx)

Zeros

If one denotes the zeros of yn(x;α,β) as αk(n)(α,β), and that of the θn(x;α,β) by βk(n)(α,β), then the following estimates exist:[2]: 82 

2n(n+α1)αk(n)(α,2)2n+α1,

and

n+α12βk(n)(α,2)n(n+α1)2,

for all α2. Moreover, all these zeros have negative real part. Sharper results can be said if one resorts to more powerful theorems regarding the estimates of zeros of polynomials (more concretely, the Parabola Theorem of Saff and Varga, or differential equations techniques).[2]: 88 [6] One result is the following:[7]

22n+α23αk(n)(α,2)2n+α1.

Particular values

The Bessel polynomials yn(x) up to n=5 are[8]

y0(x)=1y1(x)=x+1y2(x)=3x2+3x+1y3(x)=15x3+15x2+6x+1y4(x)=105x4+105x3+45x2+10x+1y5(x)=945x5+945x4+420x3+105x2+15x+1

No Bessel polynomial can be factored into lower degree polynomials with rational coefficients.[9] The reverse Bessel polynomials are obtained by reversing the coefficients. Equivalently, θk(x)=xkyk(1/x). This results in the following:

θ0(x)=1θ1(x)=x+1θ2(x)=x2+3x+3θ3(x)=x3+6x2+15x+15θ4(x)=x4+10x3+45x2+105x+105θ5(x)=x5+15x4+105x3+420x2+945x+945

See also

References

  1. 1.0 1.1 1.2 1.3 1.4 Krall, H. L.; Frink, O. (1948). "A New Class of Orthogonal Polynomials: The Bessel Polynomials". Trans. Amer. Math. Soc. 65 (1): 100–115. doi:10.2307/1990516.
  2. 2.0 2.1 2.2 2.3 2.4 Grosswald, E. (1978). Bessel Polynomials (Lecture Notes in Mathematics). New York: Springer. ISBN 978-0-387-09104-4.
  3. Berg, Christian; Vignat, Christophe (2008). "Linearization coefficients of Bessel polynomials and properties of Student-t distributions" (PDF). Constructive Approximation. 27: 15–32. doi:10.1007/s00365-006-0643-6. Retrieved 2006-08-16.
  4. Wolfram Alpha example
  5. Dita, Petre; Grama, Nicolae (May 14, 1997). "On Adomian's Decomposition Method for Solving Differential Equations". arXiv:solv-int/9705008.
  6. Saff, E. B.; Varga, R. S. (1976). "Zero-free parabolic regions for sequences of polynomials". SIAM J. Math. Anal. 7 (3): 344–357. doi:10.1137/0507028.
  7. de Bruin, M. G.; Saff, E. B.; Varga, R. S. (1981). "On the zeros of generalized Bessel polynomials. I". Indag. Math. 84 (1): 1–13.
  8. *Sloane, N. J. A. (ed.). "Sequence A001498 (Triangle a(n,k) (n >= 0, 0 <= k <= n) of coefficients of Bessel polynomials y_n(x) (exponents in increasing order).)". The On-Line Encyclopedia of Integer Sequences. OEIS Foundation.
  9. Filaseta, Michael; Trifinov, Ognian (August 2, 2002). "The Irreducibility of the Bessel Polynomials". Journal für die Reine und Angewandte Mathematik. 2002 (550): 125–140. CiteSeerX 10.1.1.6.9538. doi:10.1515/crll.2002.069.

External links