此为 Mathematica 7 文档,内容基于更早版本的 Wolfram 语言
查看最新文档(版本11.2)

Elliptic Integrals and Elliptic Functions

Even more so than for other special functions, you need to be very careful about the arguments you give to elliptic integrals and elliptic functions. There are several incompatible conventions in common use, and often these conventions are distinguished only by the specific names given to arguments or by the presence of separators other than commas between arguments.
• Amplitude Phi (used by Mathematica, in radians)
• Argument u (used by Mathematica): related to amplitude by Phi=am (u)
• Delta amplitude CapitalDelta (Phi):
• Coordinate x: x=sin (Phi)
• Characteristic n (used by Mathematica in elliptic integrals of the third kind)
• Parameter m (used by Mathematica): preceded by VerticalSeparator, as in I (PhiVerticalSeparatorm)
• Complementary parameter m1: m1=1-m
• Modulus k: preceded by comma, as in I (Phi, k); m=k2
• Modular angle Alpha: preceded by \ , as in I (Phi\Alpha); m=sin2 (Alpha)
• Nome q: preceded by comma in Theta functions; q=exp[-Pi K (1-m)/K{m)]=exp (i Pi OmegaPrime/Omega)
• Invariants g2, g3 (used by Mathematica)
• Half-periods Omega, OmegaPrime: g2=60w-4, g3=140w-6, where w=2rOmega+2sOmegaPrime
• Ratio of periods Tau: Tau=OmegaPrime/Omega
• Discriminant CapitalDelta:
• Parameters of curve a, b (used by Mathematica)
• Coordinate y (used by Mathematica): related by y2=x3+ax2+b x

Common argument conventions for elliptic integrals and elliptic functions.

JacobiAmplitude[u,m]give the amplitude Phi corresponding to argument u and parameter m
EllipticNomeQ[m]give the nome q corresponding to parameter m
InverseEllipticNomeQ[q]give the parameter m corresponding to nome q
WeierstrassInvariants[{Omega,OmegaPrime}]
give the invariants {g2, g3} corresponding to the half-periods {Omega, OmegaPrime}
WeierstrassHalfPeriods[{g2,g3}]
give the half-periods {Omega, OmegaPrime} corresponding to the invariants {g2, g3}

Converting between different argument conventions.

Elliptic Integrals

EllipticK[m]complete elliptic integral of the first kind K (m)
EllipticF[Phi,m]elliptic integral of the first kind F (PhiVerticalSeparatorm)
EllipticE[m]complete elliptic integral of the second kind E (m)
EllipticE[Phi,m]elliptic integral of the second kind E (PhiVerticalSeparatorm)
EllipticPi[n,m]complete elliptic integral of the third kind CapitalPi (nVerticalSeparatorm)
EllipticPi[n,Phi,m]elliptic integral of the third kind CapitalPi (n;PhiVerticalSeparatorm)
JacobiZeta[Phi,m]Jacobi zeta function Z (PhiVerticalSeparatorm)

Elliptic integrals.

Integrals of the form IntegralR (x, y)DifferentialDx, where R is a rational function, and y2 is a cubic or quartic polynomial in x, are known as elliptic integrals. Any elliptic integral can be expressed in terms of the three standard kinds of Legendre-Jacobi elliptic integrals.
The elliptic integral of the first kind EllipticF[Phi, m] is given for -Pi/2<Phi<Pi/2 by . This elliptic integral arises in solving the equations of motion for a simple pendulum. It is sometimes known as an incomplete elliptic integral of the first kind.
Note that the arguments of the elliptic integrals are sometimes given in the opposite order from what is used in Mathematica.
The complete elliptic integral of the first kind EllipticK[m] is given by . Note that K is used to denote the complete elliptic integral of the first kind, while F is used for its incomplete form. In many applications, the parameter m is not given explicitly, and K (m) is denoted simply by K. The complementary complete elliptic integral of the first kind KPrime (m) is given by K (1-m). It is often denoted KPrime. K and iKPrime give the "real" and "imaginary" quarter-periods of the corresponding Jacobi elliptic functions discussed in "Elliptic Functions".
The elliptic integral of the second kind EllipticE[Phi, m] is given for -Pi/2<Phi<Pi/2 by .
The complete elliptic integral of the second kind EllipticE[m] is given by . It is often denoted E. The complementary form is EPrime (m)=E (1-m).
The Jacobi zeta function JacobiZeta[Phi, m] is given by Z (PhiVerticalSeparatorm)=E (PhiVerticalSeparatorm)-E (m)F (PhiVerticalSeparatorm)/K (m).
The Heuman lambda function is given by .
The elliptic integral of the third kind EllipticPi[n, Phi, m] is given by .
The complete elliptic integral of the third kind EllipticPi[n, m] is given by .
Here is a plot of the complete elliptic integral of the second kind E (m).
In[1]:=
Click for copyable input
Out[1]=
Here is K (Alpha) with Alpha=30EmptySmallCircle.
In[2]:=
Click for copyable input
Out[2]=
The elliptic integrals have a complicated structure in the complex plane.
In[83]:=
Click for copyable input
Out[1]=

Elliptic Functions

JacobiAmplitude[u,m]amplitude function am (uVerticalSeparatorm)
JacobiSN[u,m], JacobiCN[u,m], etc.
Jacobi elliptic functions sn (uVerticalSeparatorm), etc.
InverseJacobiSN[v,m], InverseJacobiCN[v,m], etc.
inverse Jacobi elliptic functions sn-1 (vVerticalSeparatorm), etc.
EllipticTheta[a,u,q]theta functions CurlyThetaa (u, q) (a=1, ..., 4)
EllipticThetaPrime[a,u,q]derivatives of theta functions (a=1, ..., 4)
SiegelTheta[Tau,s]Siegel theta function CapitalTheta (Tau, s)
SiegelTheta[v,Tau,s]Siegel theta function CapitalTheta[v] (Tau, s)
WeierstrassP[u,{g2,g3}]Weierstrass elliptic function WeierstrassP (u;g2, g3)
WeierstrassPPrime[u,{g2,g3}]
derivative of Weierstrass elliptic function WeierstrassPPrime (u;g2, g3)
InverseWeierstrassP[p,{g2,g3}]
inverse Weierstrass elliptic function
WeierstrassSigma[u,{g2,g3}]Weierstrass sigma function Sigma (u;g2, g3)
WeierstrassZeta[u,{g2,g3}]Weierstrass zeta function Zeta (u;g2, g3)

Elliptic and related functions.

Rational functions involving square roots of quadratic forms can be integrated in terms of inverse trigonometric functions. The trigonometric functions can thus be defined as inverses of the functions obtained from these integrals.
By analogy, elliptic functions are defined as inverses of the functions obtained from elliptic integrals.
The amplitude for Jacobi elliptic functions JacobiAmplitude[u, m] is the inverse of the elliptic integral of the first kind. If u=F (PhiVerticalSeparatorm), then Phi=am (uVerticalSeparatorm). In working with Jacobi elliptic functions, the argument m is often dropped, so am (uVerticalSeparatorm) is written as am (u).
The Jacobi elliptic functions JacobiSN[u, m] and JacobiCN[u, m] are given respectively by sn (u)=sin (Phi) and cn (u)=cos (Phi), where Phi=am (uVerticalSeparatorm). In addition, JacobiDN[u, m] is given by .
There are a total of twelve Jacobi elliptic functions JacobiPQ[u, m], with the letters P and Q chosen from the set S, C, D and N. Each Jacobi elliptic function JacobiPQ[u, m] satisfies the relation pq (u)=pn (u)/qn (u), where for these purposes nn (u)=1.
There are many relations between the Jacobi elliptic functions, somewhat analogous to those between trigonometric functions. In limiting cases, in fact, the Jacobi elliptic functions reduce to trigonometric functions. So, for example, sn (uVerticalSeparator0)=sin (u), sn (uVerticalSeparator1)=tanh (u), cn (uVerticalSeparator0)=cos (u), cn (uVerticalSeparator1)=sech (u), dn (uVerticalSeparator0)=1 and dn (uVerticalSeparator1)=sech (u).
The notation Pq (u) is often used for the integrals . These integrals can be expressed in terms of the Jacobi zeta function defined in "Elliptic Integrals".
One of the most important properties of elliptic functions is that they are doubly periodic in the complex values of their arguments. Ordinary trigonometric functions are singly periodic, in the sense that f (z+sOmega)=f (z) for any integer s. The elliptic functions are doubly periodic, so that f (z+rOmega+sOmegaPrime)=f (z) for any pair of integers r and s.
The Jacobi elliptic functions sn (uVerticalSeparatorm), etc. are doubly periodic in the complex u plane. Their periods include Omega=4K (m) and OmegaPrime=4iK (1-m), where K is the complete elliptic integral of the first kind.
The choice of p and q in the notation pq (uVerticalSeparatorm) for Jacobi elliptic functions can be understood in terms of the values of the functions at the quarter periods K and iKPrime.
This shows two complete periods in each direction of the absolute value of the Jacobi elliptic function .
In[3]:=
Click for copyable input
Out[3]=
Also built into Mathematica are the inverse Jacobi elliptic functions InverseJacobiSN[v, m], InverseJacobiCN[v, m], etc. The inverse function sn-1 (vVerticalSeparatorm), for example, gives the value of u for which v=sn (uVerticalSeparatorm). The inverse Jacobi elliptic functions are related to elliptic integrals.
The four theta functions CurlyThetaa (u, q) are obtained from EllipticTheta[a, u, q] by taking a to be 1, 2, 3 or 4. The functions are defined by: CurlyTheta1 (u, q)=2q1/4 (-1)nqn (n+1)sin[ (2n+1)u], CurlyTheta2 (u, q)=2q1/4qn (n+1)cos[ (2n+1)u], CurlyTheta3 (u, q)=1+2qn2cos (2nu), CurlyTheta4 (u, q)=1+2 (-1)nqn2cos (2nu). The theta functions are often written as CurlyThetaa (u) with the parameter q not explicitly given. The theta functions are sometimes written in the form CurlyTheta (uVerticalSeparatorm), where m is related to q by q=exp[-PiK (1-m)/K (m)]. In addition, q is sometimes replaced by Tau, given by q=eiPiTau. All the theta functions satisfy a diffusion-like differential equation PartialD2CurlyTheta (u, Tau)/PartialDu2=4Pii PartialDCurlyTheta (u, Tau)/PartialDTau.
The Siegel theta function SiegelTheta[Tau, s] with Riemann square modular matrix Tau of dimension p and vector s generalizes the elliptic theta functions to complex dimension p. It is defined by , where n runs over all p-dimensional integer vectors. The Siegel theta function with characteristic SiegelTheta[Nu, Tau, s] is defined by , where the characteristic Nu is a pair of p-dimensional vectors {Alpha, Beta}.
The Jacobi elliptic functions can be expressed as ratios of the theta functions.
An alternative notation for theta functions is CapitalTheta (uVerticalSeparatorm)=CurlyTheta4 (vVerticalSeparatorm), CapitalTheta1 (uVerticalSeparatorm)=CurlyTheta3 (vVerticalSeparatorm), H (uVerticalSeparatorm)=CurlyTheta1 (v), H1 (uVerticalSeparatorm)=CurlyTheta2 (v), where v=Piu/2K (m).
The Neville theta functions can be defined in terms of the theta functions as CurlyThetas (u)=2K (m)CurlyTheta1 (vVerticalSeparatorm)/ , CurlyThetac (u)=CurlyTheta2 (vVerticalSeparatorm)/CurlyTheta2 (0VerticalSeparatorm), CurlyThetad (u)=CurlyTheta3 (vVerticalSeparatorm)/CurlyTheta3 (0VerticalSeparatorm), CurlyThetan (u)=CurlyTheta4 (vVerticalSeparatorm)/CurlyTheta4 (0VerticalSeparatorm), where v=Piu/2K (m). The Jacobi elliptic functions can be represented as ratios of the Neville theta functions.
The Weierstrass elliptic function WeierstrassP[u, {g2, g3}] can be considered as the inverse of an elliptic integral. The Weierstrass function WeierstrassP (u;g2, g3) gives the value of x for which . The function WeierstrassPPrime[u, {g2, g3}] is given by .
The Weierstrass functions are also sometimes written in terms of their fundamental half-periods Omega and OmegaPrime, obtained from the invariants g2 and g3 using WeierstrassHalfPeriods[{u, {g2, g3}].
The function InverseWeierstrassP[p, {g2, g3}] finds one of the two values of u for which p=WeierstrassP (u;g2, g3). This value always lies in the parallelogram defined by the complex number half-periods Omega and OmegaPrime.
InverseWeierstrassP[{p, q}, {g2, g3}] finds the unique value of u for which p=WeierstrassP (u;g2, g3) and q=WeierstrassPPrime (u;g2, g3). In order for any such value of u to exist, p and q must be related by q2=4p3-g2p-g3.
The Weierstrass zeta function WeierstrassZeta[u, {g2, g3}] and Weierstrass sigma function WeierstrassSigma[u, {g2, g3}] are related to the Weierstrass elliptic functions by ZetaPrime (z;g2, g3)=-WeierstrassP (z;g2, g3) and SigmaPrime (z;g2, g3)/Sigma (z;g2, g3)=Zeta (z;g2, g3).
The Weierstrass zeta and sigma functions are not strictly elliptic functions since they are not periodic.

Elliptic Modular Functions

DedekindEta[Tau]Dedekind eta function Eta (Tau)
KleinInvariantJ[Tau]Klein invariant modular function J (Tau)
ModularLambda[Tau]modular lambda function Lambda (Tau)

Elliptic modular functions.

The modular lambda function ModularLambda[Tau] relates the ratio of half-periods Tau=OmegaPrime/Omega to the parameter according to m=Lambda (Tau).
The Klein invariant modular function KleinInvariantJ[Tau] and the Dedekind eta function DedekindEta[Tau] satisfy the relations .
Modular elliptic functions are defined to be invariant under certain fractional linear transformations of their arguments. Thus for example Lambda (Tau) is invariant under any combination of the transformations TauTau+2 and TauTau/ (1-2Tau).

Generalized Elliptic Integrals and Functions

ArithmeticGeometricMean[a,b]the arithmetic-geometric mean of a and b
EllipticExp[u,{a,b}]generalized exponential associated with the elliptic curve y2=x3+ax2+bx
EllipticLog[{x,y},{a,b}]generalized logarithm associated with the elliptic curve y2=x3+ax2+bx

Generalized elliptic integrals and functions.

The definitions for elliptic integrals and functions given above are based on traditional usage. For modern algebraic geometry, it is convenient to use slightly more general definitions.
The function EllipticLog[{x, y}, {a, b}] is defined as the value of the integral , where the sign of the square root is specified by giving the value of y such that . Integrals of the form can be expressed in terms of the ordinary logarithm (and inverse trigonometric functions). You can think of EllipticLog as giving a generalization of this, where the polynomial under the square root is now of degree three.
The function EllipticExp[u, {a, b}] is the inverse of EllipticLog. It returns the list {x, y} that appears in EllipticLog. EllipticExp is an elliptic function, doubly periodic in the complex u plane.
ArithmeticGeometricMean[a, b] gives the arithmetic-geometric mean (AGM) of two numbers a and b. This quantity is central to many numerical algorithms for computing elliptic integrals and other functions. For positive reals a and b the AGM is obtained by starting with a0=a, b0=b, then iterating the transformation , until an=bn to the precision required.