# §28.33 Physical Applications

## §28.33(i) Introduction

Mathieu functions occur in practical applications in two main categories:

• Boundary-values problems arising from solution of the two-dimensional wave equation in elliptical coordinates. This yields a pair of equations of the form (28.2.1) and (28.20.1), and the appropriate solution of (28.2.1) is usually a periodic solution of integer order. See §28.33(ii).

• Initial-value problems, in which only one equation (28.2.1) or (28.20.1) is involved. See §28.33(iii).

## §28.33(ii) Boundary-Value Problems

Physical problems involving Mathieu functions include vibrational problems in elliptical coordinates; see (28.32.1). We shall derive solutions to the uniform, homogeneous, loss-free, and stretched elliptical ring membrane with mass $\rho$ per unit area, and radial tension $\tau$ per unit arc length. The wave equation

 28.33.1 $\frac{{\partial}^{2}W}{{\partial x}^{2}}+\frac{{\partial}^{2}W}{{\partial y}^{% 2}}-\frac{\rho}{\tau}\frac{{\partial}^{2}W}{{\partial t}^{2}}=0,$

with $W(x,y,t)=e^{\mathrm{i}\omega t}V(x,y)$, reduces to (28.32.2) with $k^{2}=\omega^{2}\rho/{\tau}$. In elliptical coordinates (28.32.2) becomes (28.32.3). The separated solutions $V_{n}(\xi,\eta)$ must be $2\pi$-periodic in $\eta$, and have the form

 28.33.2 $V_{n}(\xi,\eta)=\left(c_{n}\mathop{{\mathrm{M}^{(1)}_{n}}\/}\nolimits\!\left(% \xi,\sqrt{q}\right)+d_{n}\mathop{{\mathrm{M}^{(2)}_{n}}\/}\nolimits\!\left(\xi% ,\sqrt{q}\right)\right)\mathop{\mathrm{me}_{n}\/}\nolimits\!\left(\eta,q\right),$

where $q=\tfrac{1}{4}c^{2}k^{2}$ and $a_{n}(q)$ or $b_{n}(q)$ is the separation constant; compare (28.12.11), (28.20.11), and (28.20.12). Here $c_{n}$ and $d_{n}$ are constants. The boundary conditions for $\xi=\xi_{0}$ (outer clamp) and $\xi=\xi_{1}$ (inner clamp) yield the following equation for $q$:

 28.33.3 $\mathop{{\mathrm{M}^{(1)}_{n}}\/}\nolimits\!\left(\xi_{0},\sqrt{q}\right)% \mathop{{\mathrm{M}^{(2)}_{n}}\/}\nolimits\!\left(\xi_{1},\sqrt{q}\right)-% \mathop{{\mathrm{M}^{(1)}_{n}}\/}\nolimits\!\left(\xi_{1},\sqrt{q}\right)% \mathop{{\mathrm{M}^{(2)}_{n}}\/}\nolimits\!\left(\xi_{0},\sqrt{q}\right)=0.$ Symbols: $\mathop{{\mathrm{M}^{(\NVar{j})}_{\NVar{\nu}}}\/}\nolimits\!\left(\NVar{z},% \NVar{h}\right)$: modified Mathieu function, $q=h^{2}$: parameter, $n$: integer and $\xi$: variable Referenced by: §28.33(ii) Permalink: http://dlmf.nist.gov/28.33.E3 Encodings: TeX, pMML, png See also: Annotations for 28.33(ii)

If we denote the positive solutions $q$ of (28.33.3) by $q_{n,m}$, then the vibration of the membrane is given by $\omega^{2}_{n,m}=\ifrac{4q_{n,m}\tau}{(c^{2}\rho)}$. The general solution of the problem is a superposition of the separated solutions.

For a visualization see Gutiérrez-Vega et al. (2003), and for references to other boundary-value problems see:

• McLachlan (1947, Chapters XVI–XIX) for applications of the wave equation to vibrational systems, electrical and thermal diffusion, electromagnetic wave guides, elliptical cylinders in viscous fluids, and diffraction of sound and electromagnetic waves.

• Meixner and Schäfke (1954, §§4.3, 4.4) for elliptic membranes and electromagnetic waves.

• Daymond (1955) for vibrating systems.

• Troesch and Troesch (1973) for elliptic membranes.

• Alhargan and Judah (1995), Bhattacharyya and Shafai (1988), and Shen (1981) for ring antennas.

• Alhargan and Judah (1992), Germey (1964), Ragheb et al. (1991), and Sips (1967) for electromagnetic waves.

More complete bibliographies will be found in McLachlan (1947) and Meixner and Schäfke (1954).

## §28.33(iii) Stability and Initial-Value Problems

If the parameters of a physical system vary periodically with time, then the question of stability arises, for example, a mathematical pendulum whose length varies as $\mathop{\cos\/}\nolimits\!\left(2\omega t\right)$. The equation of motion is given by

 28.33.4 $w^{\prime\prime}(t)+\left(b-f\mathop{\cos\/}\nolimits\!\left(2\omega t\right)% \right)w(t)=0,$ Symbols: $\mathop{\cos\/}\nolimits\NVar{z}$: cosine function, $w(z)$: Mathieu’s equation solution and $b$: fixed Permalink: http://dlmf.nist.gov/28.33.E4 Encodings: TeX, pMML, png See also: Annotations for 28.33(iii)

with $b$, $f$, and $\omega$ positive constants. Substituting $z=\omega t$, $a=\ifrac{b}{\omega^{2}}$, and $2q=\ifrac{f}{\omega^{2}}$, we obtain Mathieu’s standard form (28.2.1).

As $\omega$ runs from $0$ to $+\infty$, with $b$ and $f$ fixed, the point $(q,a)$ moves from $\infty$ to $0$ along the ray $\mathcal{L}$ given by the part of the line $a=(2b/f)q$ that lies in the first quadrant of the $(q,a)$-plane. Hence from §28.17 the corresponding Mathieu equation is stable or unstable according as $(q,a)$ is in the intersection of $\mathcal{L}$ with the colored or the uncolored open regions depicted in Figure 28.17.1. In particular, the equation is stable for all sufficiently large values of $\omega$.

For points $(q,a)$ that are at intersections of $\mathcal{L}$ with the characteristic curves $a=\mathop{a_{n}\/}\nolimits\!\left(q\right)$ or $a=\mathop{b_{n}\/}\nolimits\!\left(q\right)$, a periodic solution is possible. However, in response to a small perturbation at least one solution may become unbounded.

References for other initial-value problems include:

• McLachlan (1947, Chapter XV) for amplitude distortion in moving-coil loud-speakers, frequency modulation, dynamical systems, and vibration of stretched strings.

• Vedeler (1950) for ships rolling among waves.

• Meixner and Schäfke (1954, §§4.1, 4.2, and 4.7) for quantum mechanical problems and rotation of molecules.

• Aly et al. (1975) for scattering theory.

• Hunter and Kuriyan (1976) and Rushchitsky and Rushchitska (2000) for wave mechanics.

• Fukui and Horiguchi (1992) for quantum theory.

• Jager (1997, 1998) for relativistic oscillators.

• Torres-Vega et al. (1998) for Mathieu functions in phase space.