Jacobi elliptic functions

In mathematics, the Jacobi elliptic functions are a set of basic elliptic functions, and auxiliary theta functions, that are of historical importance. Many of their features show up in important structures and have direct relevance to some applications (e.g. the equation of a pendulumalso see pendulum (mathematics)). They also have useful analogies to the functions of trigonometry, as indicated by the matching notation sn for sin. The Jacobi elliptic functions are used more often in practical problems than the Weierstrass elliptic functions as they do not require notions of complex analysis to be defined and/or understood. They were introduced by Carl Gustav Jakob Jacobi (1829).

Introduction

Auxiliary rectangle construction

There are twelve Jacobian elliptic functions. Each of the twelve corresponds to an arrow drawn from one corner of a rectangle to another. The corners of the rectangle are labeled, by convention, s, c, d and n. The rectangle is understood to be lying on the complex plane, so that s is at the origin, c is at the point K on the real axis, d is at the point K + iK' and n is at point iK' on the imaginary axis. The numbers K and K' are called the quarter periods. The twelve Jacobian elliptic functions are then pq, where each of p and q is a different one of the letters s, c, d, n.

The Jacobian elliptic functions are then the unique doubly periodic, meromorphic functions satisfying the following three properties:

More generally, there is no need to impose a rectangle; a parallelogram will do. However, if K and iK' are kept on the real and imaginary axis respectively, then the Jacobi elliptic functions pq u will be real functions when u is real.

Notation

The elliptic functions can be given in a variety of notations, which can make the subject unnecessarily confusing. Elliptic functions are functions of two variables. The first variable might be given in terms of the amplitude φ, or more commonly, in terms of u given below. The second variable might be given in terms of the parameter m, or as the elliptic modulus k, where k2 = m, or in terms of the modular angle α, where m =  sin2 α. A more extensive review and definition of these alternatives, their complements, and the associated notation schemes are given in the articles on elliptic integrals and quarter period.

Definition as inverses of elliptic integrals

Elliptic Jacobi function, sn, corresponding to k=0.8, generated using a version of the Domain coloring method.[1]

The above definition, in terms of the unique meromorphic functions satisfying certain properties, is quite abstract. There is a simpler, but completely equivalent definition, giving the elliptic functions as inverses of the incomplete elliptic integral of the first kind. Let

u=\int_0^\phi \frac{\mathrm d\theta} {\sqrt {1-m \sin^2 \theta}}\, .

Then the elliptic function sn u is given by

\operatorname {sn}\; u = \sin \phi\,

and cn u is given by

\operatorname {cn}\; u = \cos \phi

and

\operatorname {dn}\; u = \sqrt {1-m\sin^2 \phi}\, .
Model of amplitude (measured along vertical axis) as a function of independent variables u and k

Here, the angle \phi is called the amplitude. On occasion, dn u = Δ(u) is called the delta amplitude. In the above, the value m is a free parameter, usually taken to be real, 0  m  1, and so the elliptic functions can be thought of as being given by two variables, the amplitude \phi and the parameter m.

The remaining nine elliptic functions are easily built from the above three, and are given in a section below.

Note that when \phi=\pi/2, that u then equals the quarter period K.

Definition as trigonometry

 \cos \theta, \sin \theta are defined on the unit circle, with r = 1. Similarly, Jacobi elliptic functions are defined on the unit ellipse, with a = 1. Let


\begin{align}
& x^2 + \frac{y^2}{b^2} = 1, \\
& m = 1 - \frac 1 {b^2} < 1, \\
& x = r \cos \theta, \quad y = r \sin \theta
\end{align}

then:

 r( \theta,m) =  \frac{1} {\sqrt {1-m \sin^2 \theta}}\, .

Replace radians on the unit circle by u on the unit ellipse, where:

 \cos \varphi= x ,
u = u(\varphi,m)=\int_0^\varphi r(\theta,m) \, d\theta

then:

\operatorname{cn}(u,m)=x,\quad \operatorname{sn}(u,m) = \frac y b,\quad \operatorname{dn}(u,m) = \frac 1 {r(\varphi,m)} .

Definition in terms of theta functions

Equivalently, Jacobi elliptic functions can be defined in terms of his theta functions. If we abbreviate \vartheta(0;\tau) as \vartheta, and \vartheta_{01}(0;\tau), \vartheta_{10}(0;\tau), \vartheta_{11}(0;\tau) respectively as \vartheta_{01}, \vartheta_{10}, \vartheta_{11} (the theta constants) then the elliptic modulus k is k=\left({\vartheta_{10} \over \vartheta}\right)^2. If we set u = \pi \vartheta^2 z, we have

Elliptic Jacobi function, cn, k=0.8
Elliptic Jacobi function, dn, k=0.8
\mbox{sn}(u; k) = -{\vartheta \vartheta_{11}(z;\tau) \over \vartheta_{10} \vartheta_{01}(z;\tau)}


\mbox{cn}(u; k) = {\vartheta_{01} \vartheta_{10}(z;\tau) \over \vartheta_{10} \vartheta_{01}(z;\tau)}


\mbox{dn}(u; k) = {\vartheta_{01} \vartheta(z;\tau) \over \vartheta \vartheta_{01}(z;\tau)}

Since the Jacobi functions are defined in terms of the elliptic modulus k(τ), we need to invert this and find τ in terms of k. We start from k' = \sqrt{1-k^2}, the complementary modulus. As a function of τ it is

k'(\tau) = \left({\vartheta_{01} \over \vartheta}\right)^2.

Let us first define

\ell = {1 \over 2} {1-\sqrt{k'} \over 1+\sqrt{k'}} =
{1 \over 2} {\vartheta - \vartheta_{01} \over \vartheta + \vartheta_{01}}.

Then define the nome q as q = \exp (\pi i \tau) and expand \ell as a power series in the nome q, we obtain

\ell = {q + q^9 + q^{25} + \cdots \over 1 + 2q^4 + 2q^{16} + \cdots}.

Reversion of series now gives

q = \ell + 2\ell^5 + 15\ell^9 + 150\ell^{13} + 1707\ell^{17} + 20910\ell^{21} + 268616\ell^{25} + \cdots.

Since we may reduce to the case where the imaginary part of τ is greater than or equal to 1/2 sqrt(3), we can assume the absolute value of q is less than or equal to exp(-1/2 sqrt(3) π) ~ 0.0658; for values this small the above series converges very rapidly and easily allows us to find the appropriate value for q.

Minor functions

Reversing the order of the two letters of the function name results in the reciprocals of the three functions above:


\begin{align}
\operatorname{ns}(u) & = \frac{1}{\operatorname{sn}(u)} \\[8pt]
\operatorname{nc}(u) & = \frac{1}{\operatorname{cn}(u)} \\[8pt]
\operatorname{nd}(u) & = \frac{1}{\operatorname{dn}(u)}
\end{align}

Similarly, the ratios of the three primary functions correspond to the first letter of the numerator followed by the first letter of the denominator:

Elliptic Jacobi function, sc, k = 0.8

\begin{align}
\operatorname{sc}(u) & = \frac{\operatorname{sn}(u)}{\operatorname{cn}(u)} \\[8pt]
\operatorname{sd}(u) & = \frac{\operatorname{sn}(u)}{\operatorname{dn}(u)} \\[8pt]
\operatorname{dc}(u) & = \frac{\operatorname{dn}(u)}{\operatorname{cn}(u)} \\[8pt]
\operatorname{ds}(u) & = \frac{\operatorname{dn}(u)}{\operatorname{sn}(u)} \\[8pt]
\operatorname{cs}(u) & = \frac{\operatorname{cn}(u)}{\operatorname{sn}(u)} \\[8pt]
\operatorname{cd}(u) & = \frac{\operatorname{cn}(u)}{\operatorname{dn}(u)}
\end{align}

More compactly, we have

\operatorname{pq}(u)=\frac{\operatorname{pr}(u)}{\operatorname{qr}(u)}

where each of p, q, and r is any of the letters s, c, d, n, with the understanding that ss = cc = dd = nn = 1.

(This notation is due to Gudermann and Glaisher and is not Jacobi's original notation.)

Addition theorems

The functions satisfy the two algebraic relations

\operatorname{cn}^2(u,k) + \operatorname{sn}^2(u,k) = 1,\,
\operatorname{dn}^2(u,k) + k^2 \ \operatorname{sn}^2(u,k) = 1.\,

From this we see that (cn, sn, dn) parametrizes an elliptic curve which is the intersection of the two quadrics defined by the above two equations. We now may define a group law for points on this curve by the addition formulas for the Jacobi functions


\begin{align}
\operatorname{cn}(x+y) & =
{\operatorname{cn}(x)\;\operatorname{cn}(y)
- \operatorname{sn}(x)\;\operatorname{sn}(y)\;\operatorname{dn}(x)\;\operatorname{dn}(y)
\over {1 - k^2 \;\operatorname{sn}^2 (x) \;\operatorname{sn}^2 (y)}}, \\[8pt]
\operatorname{sn}(x+y) & =
{\operatorname{sn}(x)\;\operatorname{cn}(y)\;\operatorname{dn}(y) +
\operatorname{sn}(y)\;\operatorname{cn}(x)\;\operatorname{dn}(x)
\over {1 - k^2 \;\operatorname{sn}^2 (x)\; \operatorname{sn}^2 (y)}}, \\[8pt]
\operatorname{dn}(x+y) & =
{\operatorname{dn}(x)\;\operatorname{dn}(y)
- k^2 \;\operatorname{sn}(x)\;\operatorname{sn}(y)\;\operatorname{cn}(x)\;\operatorname{cn}(y)
\over {1 - k^2 \;\operatorname{sn}^2 (x)\; \operatorname{sn}^2 (y)}}.
\end{align}

Relations between squares of the functions


-\operatorname{dn}^2(u)+m_1= -m\;\operatorname{cn}^2(u) = m\;\operatorname{sn}^2(u)-m

-m_1\;\operatorname{nd}^2(u)+m_1= -mm_1\;\operatorname{sd}^2(u) = m\;\operatorname{cd}^2(u)-m

m_1\;\operatorname{sc}^2(u)+m_1= m_1\;\operatorname{nc}^2(u) = \operatorname{dc}^2(u)-m

\operatorname{cs}^2(u)+m_1=\operatorname{ds}^2(u)=\operatorname{ns}^2(u)-m

where m + m1 = 1 and m = k2.

Additional relations between squares can be obtained by noting that pq2 · qp2 = 1 and that pq = pr / qr where p, q, r are any of the letters s, c, d, n and ss = cc = dd = nn = 1.

Expansion in terms of the nome

Let the nome be q=\exp(-\pi K'/K) and let the argument be v=\pi u /(2K). Then the functions have expansions as Lambert series

\operatorname{sn}(u)=\frac{2\pi}{K\sqrt{m}}
\sum_{n=0}^\infty \frac{q^{n+1/2}}{1-q^{2n+1}} \sin ((2n+1)v),
\operatorname{cn}(u)=\frac{2\pi}{K\sqrt{m}}
\sum_{n=0}^\infty \frac{q^{n+1/2}}{1+q^{2n+1}} \cos ((2n+1)v),
\operatorname{dn}(u)=\frac{\pi}{2K} + \frac{2\pi}{K}
\sum_{n=1}^\infty \frac{q^{n}}{1+q^{2n}} \cos (2nv).

Jacobi elliptic functions as solutions of nonlinear ordinary differential equations

The derivatives of the three basic Jacobi elliptic functions are:


\frac{\mathrm{d}}{\mathrm{d}z}\, \mathrm{sn}\,(z) = \mathrm{cn}\,(z)\, \mathrm{dn}\,(z),


\frac{\mathrm{d}}{\mathrm{d}z}\, \mathrm{cn}\,(z) = -\mathrm{sn}\,(z)\, \mathrm{dn}\,(z),


\frac{\mathrm{d}}{\mathrm{d}z}\, \mathrm{dn}\,(z) = - k^2 \mathrm{sn}\,(z)\, \mathrm{cn}\,(z).

With the addition theorems above and for a given k with 0 < k < 1 they therefore are solutions to the following nonlinear ordinary differential equations:

\frac{\mathrm{d}^2 y}{\mathrm{d}x^2} + (1+k^2) y - 2 k^2 y^3 = 0
and
 \left(\frac{\mathrm{d} y}{\mathrm{d}x}\right)^2 = (1-y^2) (1-k^2 y^2)
\frac{\mathrm{d}^2 y}{\mathrm{d}x^2} + (1-2k^2) y + 2 k^2 y^3 = 0
and
 \left(\frac{\mathrm{d} y}{\mathrm{d}x}\right)^2 = (1-y^2) (1-k^2 + k^2 y^2)
\frac{\mathrm{d}^2 y}{\mathrm{d}x^2} - (2 - k^2) y + 2 y^3 = 0
and
 \left(\frac{\mathrm{d} y}{\mathrm{d}x}\right)^2 = (y^2 - 1) (1 - k^2 - y^2)

Inverse functions

The inverses of the Jacobi elliptic functions can be defined similarly to the inverse trigonometric functions; if x=\mathrm{sn}(\xi, k), \xi=\mathrm{arcsn}(x, k). They can be represented as elliptic integrals,[2][3] and power series representations have been found.[4]

Map projection

The Peirce quincuncial projection is a map projection based on Jacobian elliptic functions.

See also

Notes

  1. http://nbviewer.ipython.org/github/empet/Math/blob/master/DomainColoring.ipynb
  2. Reinhardt, W. P.; Walker, P. L. (2010), "§22.15 Inverse Functions", in Olver, Frank W. J.; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W., NIST Handbook of Mathematical Functions, Cambridge University Press, ISBN 978-0521192255, MR 2723248
  3. Ehrhardt, Wolfgang. "The AMath and DAMath Special Functions: Reference Manual and Implementation Notes" (PDF). p. 42. Retrieved 17 July 2013.
  4. Carlson, B. C. (2008). "Power series for inverse Jacobian elliptic functions" (PDF). Mathematics of Computation 77: 1615–1621. doi:10.1090/s0025-5718-07-02049-2. Retrieved 17 July 2013.

References

External links