Download (direct link):
For any nonnegative constant p, the differential equation
t2 X'(t) + tx (t) + (t2 — p2 )x(t) = 0
is known as Bessel’s equation of order p, and its solutions are the Bessel functions of order p. In normalized form, Bessel’s equation becomes
X'(t) + -tX(t) +
t2 - p2
x(t) = 0
From this we can see that tp(t) = 1 and t2q(t) = t2 — p2, so that tp(t) and t2q(t) are analytic at t0 = 0. Therefore zero is a regular singular point and, using equation (8), we find that the indicial equation (with P0 = 1, Q0 = — p2)
r(r — 1) + r — p2 = r2 — p2 = 0 Application of Frobenius’ Theorem yields a solution Jp given by the formula
Jp (t) = tpJ2
^ 22nn! (p + 1)( p + 2) ¦¦¦( p + n)
The function Jp(t) is called the Bessel function of order p of the first kind. The series converges and is bounded for all t. If p is not an integer, it can be shown that a second solution of Bessel’s equation is J— p(t) and that the general solution of Bessel’s equation is a linear combination of Jp(t) and
J— p (t).
For the special case p = 0, we get the function J0 (t) used in the aging spring model in the second submodule of Module 11:
2n t2 t
= 1 - - +
4 64 2304
Note that even though t = 0 is a singular point of the Bessel equation of order zero, the value of J0(0) is finite [J0(0) = 1]. See Figure 11.2.
? Check that J0 (t) is a solution of Bessel’s equation of order 0.
When p is an integer we have to work much harder to get a second solution that is linearly independent of Jp(t). The result is a function Yp(t) called the Bessel function of order p of the second kind. The general formula for Yp (t) is extremely complicated. We show only the special case Y0 (t), used in the aging spring model:
Actually ó is an unending decimal and non-repeating (or so most mathematicians believe), and 0.5772 gives the first four digits.
Y0 (t) =
/ t \
(y + ln - +
where Hn = 1 + (1/2) + (1/3) + ¦¦¦ + (1/n) and y is Euler’s constant: ó =
3(Hn — lnn) ú 0.5772.
210 Chapter 11
Figure 11.2: The graph of J0(t) [dark] looks like the graph of the decaying sinusoid *J2/ntcos(t— ˙/4) [light].
The general solution of Bessel’s equation of integer order p is
x(t) = ci Jp (t) + C2Yp (t) (10)
for arbitrary constants c1 and c2. An important thing to note here is that the value of Yp (t) at t = 0 does reflect the singularity at t = 0; in fact, Yp (t) ^ —<x as t ^ 0+, so that a solution having the form given in equation (10) is bounded only if c = 0.
Bessel functions appear frequently in applications involving cylindrical geometry and have been extensively studied. In fact, except for the functions you studied in calculus, Bessel functions are the most widely used functions in science and engineering.
¦ Transforming Bessel's Equation to the Aging Spring Equation
See “Aging Springs” in Bessel’s equation of order zero can be transformed into the aging spring equa-
M°duie n. tion x" + e-atx = 0. To do this, we take
t = (2/a) ln(2/as) (11)
where the new independent variable s is assumed to be positive. Then we can use the chain rule to find the first two derivatives of the displacement x of the
Transforming Bessel's Equation to the Aging Spring Equation
We use w in place of x in the aging spring section of Module 11.
aging spring with respect to s:
dx dx dt dx
ds dt ds dt
d2x d ds2 ds
d2 x dt
dt2 ds V as/ + Ë as2 dt2 as as dt as2
d2x 4 dx 2
dt2 (as)2 dt as2
Bessel's equation of order p = 0 is given by:
2 d2 x dx 2
s2—-r + s— + Szx=0 ds2 ds
and when we substitute in the derivatives we just found, we obtain
2 ( d2 x 4 dx 2 \ dx ( 2 \ 2
I “TT 7—ă÷ + -----9 I + s~r (---I + sx = 0
dt2 (as)2 dt as2 j dt \ as)