Digital Library of Mathematical Functions
About the Project
NIST
28 Mathieu Functions and Hill’s EquationHill’s Equation

§28.29 Definitions and Basic Properties

Contents

§28.29(i) Hill’s Equation

A generalization of Mathieu’s equation (28.2.1) is Hill’s equation

with

28.29.2Q(z+\pi)=Q(z),

and

Q(z) is either a continuous and real-valued function for z\in\Real or an analytic function of z in a doubly-infinite open strip that contains the real axis. \pi is the minimum period of Q.

§28.29(ii) Floquet’s Theorem and the Characteristic Exponent

Let \nu be a real or complex constant satisfying (without loss of generality)

28.29.6-1<\realpart{\nu}\leq 1

throughout this section. Then (28.29.1) has a nontrivial solution w(z) with the pseudoperiodic property

iff e^{{\pi i\nu}} is an eigenvalue of the matrix

Equivalently,

This is the characteristic equation of (28.29.1), and \mathop{\cos\/}\nolimits\!\left(\pi\nu\right) is an entire function of \lambda. Given \lambda together with the condition (28.29.6), the solutions \pm\nu of (28.29.9) are the characteristic exponents of (28.29.1). A solution satisfying (28.29.7) is called a Floquet solution with respect to \nu (or Floquet solution). It has the form

28.29.10F_{{\nu}}(z)=e^{{i\nu z}}P_{{\nu}}(z),

where the function P_{{\nu}}(z) is \pi-periodic.

If \nu (\neq 0,1) is a solution of (28.29.9), then F_{\nu}(z), F_{{-\nu}}(z) comprise a fundamental pair of solutions of Hill’s equation.

If \nu=0 or 1, then (28.29.1) has a nontrivial solution P(z) which is periodic with period \pi (when \nu=0) or 2\pi (when \nu=1). Let w(z) be a solution linearly independent of P(z). Then

where c is a constant. The case c=0 is equivalent to

The solutions of period \pi or 2\pi are exceptional in the following sense. If (28.29.1) has a periodic solution with minimum period n\pi, n=3,4,\dots, then all solutions are periodic with period n\pi.

Furthermore, for each solution w(z) of (28.29.1)

A nontrivial solution w(z) is either a Floquet solution with respect to \nu, or w(z+\pi)-e^{{i\nu\pi}}w(z) is a Floquet solution with respect to -\nu.

In the symmetric case Q(z)=Q(-z), w_{{\mbox{\tiny I}}}(z,\lambda) is an even solution and w_{{\mbox{\tiny II}}}(z,\lambda) is an odd solution; compare §28.2(ii). (28.29.9) reduces to

The cases \nu=0 and \nu=1 split into four subcases as in (28.2.21) and (28.2.22). The \pi-periodic or \pi-antiperiodic solutions are multiples of w_{{\mbox{\tiny I}}}(z,\lambda),w_{{\mbox{\tiny II}}}(z,\lambda), respectively.

For details and proofs see Magnus and Winkler (1966, §1.3).

§28.29(iii) Discriminant and Eigenvalues in the Real Case

Q(x) is assumed to be real-valued throughout this subsection.

The function

is called the discriminant of (28.29.1). It is an entire function of \lambda. Its order of growth for |\lambda|\to\infty is exactly \tfrac{1}{2}; see Magnus and Winkler (1966, Chapter II, pp. 19–28).

For a given \nu, the characteristic equation \bigtriangleup(\lambda)-2\mathop{\cos\/}\nolimits\!\left(\pi\nu\right)=0 has infinitely many roots \lambda. Conversely, for a given \lambda, the value of \bigtriangleup(\lambda) is needed for the computation of \nu. For this purpose the discriminant can be expressed as an infinite determinant involving the Fourier coefficients of Q(x); see Magnus and Winkler (1966, §2.3, pp. 28–36).

To every equation (28.29.1), there belong two increasing infinite sequences of real eigenvalues:

28.29.16\lambda_{n},\;n=0,1,2,\dots,\mbox{ with $\bigtriangleup(\lambda_{n})=2$},
28.29.17\mu_{n},\;n=1,2,3,\dots,\mbox{ with $\bigtriangleup(\mu_{n})=-2$}.

In consequence, (28.29.1) has a solution of period \pi iff \lambda=\lambda_{n}, and a solution of period 2\pi iff \lambda=\mu_{n}. Both \lambda_{n} and \mu_{n}\to\infty as n\to\infty, and interlace according to the inequalities

28.29.18\lambda_{0}<\mu_{1}\leq\mu_{2}<\lambda_{1}\leq\lambda_{2}<\mu_{3}\leq\mu_{4}<%
\lambda_{3}\leq\lambda_{4}<\cdots.

Assume that the second derivative of Q(x) in (28.29.1) exists and is continuous. Then with

we have for m\to\infty

28.29.20
\mu_{{2m-1}}-(2m-1)^{2}-\dfrac{N}{(4m)^{2}}=\mathop{o\/}\nolimits\!\left(m^{{-%
2}}\right),
\mu_{{2m}}-(2m-1)^{2}-\dfrac{N}{(4m)^{2}}=\mathop{o\/}\nolimits\!\left(m^{{-2}%
}\right),
28.29.21
\lambda_{{2m-1}}-(2m)^{2}-\dfrac{N}{(4m)^{2}}=\mathop{o\/}\nolimits\!\left(m^{%
{-2}}\right),
\lambda_{{2m}}-(2m)^{2}-\dfrac{N}{(4m)^{2}}=\mathop{o\/}\nolimits\!\left(m^{{-%
2}}\right).

If Q(x) has k continuous derivatives, then as m\to\infty

see Hochstadt (1963).

For further results, especially when Q(z) is analytic in a strip, see Weinstein and Keller (1987).