Volume of an n-ball

From testwiki
Revision as of 01:46, 28 October 2024 by imported>Citation bot (Altered date. | Use this bot. Report bugs. | Suggested by Susmuffin | Category:CS1 maint: date format | #UCB_Category 163/423)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Template:Short description

Volumes of balls in dimensions 0 through 25; unit ball in red.

In geometry, a ball is a region in a space comprising all points within a fixed distance, called the radius, from a given point; that is, it is the region enclosed by a sphere or hypersphere. An Template:Math-ball is a ball in an Template:Math-dimensional Euclidean space. The volume of a Template:Math-ball is the Lebesgue measure of this ball, which generalizes to any dimension the usual volume of a ball in 3-dimensional space. The volume of a Template:Mvar-ball of radius Template:Mvar is RnVn, where Vn is the volume of the [[unit ball|unit Template:Mvar-ball]], the Template:Mvar-ball of radius Template:Math.

The real number Vn can be expressed via a two-dimension recurrence relation. Closed-form expressions involve the gamma, factorial, or double factorial function. The volume can also be expressed in terms of An, the area of the [[unit n-sphere|unit Template:Mvar-sphere]].

Formulas

Template:See also

The first volumes are as follows:

Dimension Volume of a ball of radius Template:Math Radius of a ball of volume Template:Math
0 1 (all 0-balls have volume 1)
1 2R V2=0.5×V
2 πR23.142×R2 V1/2π0.564×V12
3 4π3R34.189×R3 (3V4π)1/30.620×V1/3
4 π22R44.935×R4 (2V)1/4π0.671×V1/4
5 8π215R55.264×R5 (15V8π2)1/50.717×V1/5
6 π36R65.168×R6 (6V)1/6π0.761×V1/6
7 16π3105R74.725×R7 (105V16π3)1/70.801×V1/7
8 π424R84.059×R8 (24V)1/8π0.839×V1/8
9 32π4945R93.299×R9 (945V32π4)1/90.876×V1/9
10 π5120R102.550×R10 (120V)1/10π0.911×V1/10
11 64π510395R111.884×R11 (10395V64π5)1/110.944×V1/11
12 π6720R121.335×R12 (720V)1/12π0.976×V1/12
13 128π6135135R130.911×R13 (135135V128π6)1/131.007×V1/13
14 π75040R140.599×R14 (5040V)1/14π1.037×V1/14
15 256π72027025R150.381×R15 (2027025V256π7)1/151.066×V1/15
n Vn(R) Rn(V)

Closed form

The Template:Math-dimensional volume of a Euclidean ball of radius Template:Math in Template:Math-dimensional Euclidean space is:[1]

Vn(R)=πn/2Γ(n2+1)Rn,

where Template:Math is Euler's gamma function. The gamma function is offset from but otherwise extends the factorial function to non-integer arguments. It satisfies Template:Math if Template:Math is a positive integer and Template:Math if Template:Math is a non-negative integer.

Two-dimension recurrence relation

The volume can be computed without use of the Gamma function. As is proved below using a vector-calculus double integral in polar coordinates, the volume Template:Math of an Template:Math-ball of radius Template:Math can be expressed recursively in terms of the volume of an Template:Math-ball, via the interleaved recurrence relation:

Vn(R)={1if n=0,2Rif n=1,2πnR2×Vn2(R)otherwise.

This allows computation of Template:Math in approximately Template:Math steps.

Alternative forms

The volume can also be expressed in terms of an Template:Math-ball using the one-dimension recurrence relation:

V0(R)=1,Vn(R)=Γ(n2+12)πΓ(n2+1)RVn1(R).

Inverting the above, the radius of an Template:Math-ball of volume Template:Math can be expressed recursively in terms of the radius of an Template:Math- or Template:Math-ball:

Rn(V)=(12n)1/n(Γ(n2)V)2/(n(n2))Rn2(V),Rn(V)=Γ(n2+1)1/nΓ(n2+12)1/(n1)V1/(n(n1))Rn1(V).

Using explicit formulas for particular values of the gamma function at the integers and half-integers gives formulas for the volume of a Euclidean ball in terms of factorials. For non-negative integer Template:Mvar, these are:

V2k(R)=πkk!R2k,V2k+1(R)=2(k!)(4π)k(2k+1)!R2k+1.

The volume can also be expressed in terms of double factorials. For a positive odd integer Template:Math, the double factorial is defined by

(2k+1)!!=(2k+1)(2k1)531.

The volume of an odd-dimensional ball is

V2k+1(R)=2(2π)k(2k+1)!!R2k+1.

There are multiple conventions for double factorials of even integers. Under the convention in which the double factorial satisfies

(2k)!!=(2k)(2k2)422/π=2kk!2/π,

the volume of an Template:Math-dimensional ball is, regardless of whether Template:Math is even or odd,

Vn(R)=2(2π)(n1)/2n!!Rn.

Instead of expressing the volume Template:Math of the ball in terms of its radius Template:Math, the formulas can be inverted to express the radius as a function of the volume:

Rn(V)=Γ(n2+1)1/nπV1/n=(n!!V2(2π)(n1)/2)1/nR2k(V)=(k!V)1/(2k)π,R2k+1(V)=((2k+1)!V2(k!)(4π)k)1/(2k+1).

Approximation for high dimensions

Stirling's approximation for the gamma function can be used to approximate the volume when the number of dimensions is high.

Vn(R)1nπ(2πen)n/2Rn.
Rn(V)(πn)1/(2n)n2πeV1/n.

In particular, for any fixed value of Template:Math the volume tends to a limiting value of 0 as Template:Math goes to infinity. Which value of Template:Mvar maximizes Template:Math depends upon the value of Template:Mvar; for example, the volume Template:Math is increasing for Template:Math, achieves its maximum when Template:Math, and is decreasing for Template:Math.[2]

Also, there is an asymptotic formula for the surface area[3]limn1nlnAn1(n)=12(ln(2π)+1)

Relation with surface area

Surface areas of hyperspheres in dimensions 0 through 25

Let Template:Math denote the hypervolume of the [[n-sphere|Template:Math-sphere]] of radius Template:Math. The Template:Math-sphere is the Template:Math-dimensional boundary (surface) of the Template:Math-dimensional ball of radius Template:Math, and the sphere's hypervolume and the ball's hypervolume are related by:

An1(R)=ddRVn(R)=nRVn(R).

Thus, Template:Math inherits formulas and recursion relationships from Template:Math, such as

An1(R)=2πn/2Γ(n2)Rn1.

There are also formulas in terms of factorials and double factorials.

Proofs

There are many proofs of the above formulas.

The volume is proportional to the Template:Mathth power of the radius

An important step in several proofs about volumes of Template:Math-balls, and a generally useful fact besides, is that the volume of the Template:Math-ball of radius Template:Math is proportional to Template:Math:

Vn(R)Rn.

The proportionality constant is the volume of the unit ball.

This is a special case of a general fact about volumes in Template:Math-dimensional space: If Template:Math is a body (measurable set) in that space and Template:Math is the body obtained by stretching in all directions by the factor Template:Math then the volume of Template:Math equals Template:Math times the volume of Template:Math. This is a direct consequence of the change of variables formula:

V(RK)=RKdx=KRndy=RnV(K)

where Template:Math and the substitution Template:Math was made.

Another proof of the above relation, which avoids multi-dimensional integration, uses induction: The base case is Template:Math, where the proportionality is obvious. For the inductive step, assume that proportionality is true in dimension Template:Math. Note that the intersection of an n-ball with a hyperplane is an Template:Math-ball. When the volume of the Template:Math-ball is written as an integral of volumes of Template:Math-balls:

Vn(R)=RRVn1(R2x2)dx,

it is possible by the inductive hypothesis to remove a factor of Template:Math from the radius of the Template:Math-ball to get:

Vn(R)=Rn1RRVn1(1(xR)2)dx.

Making the change of variables Template:Math leads to:

Vn(R)=Rn11Vn1(1t2)dt=RnVn(1),

which demonstrates the proportionality relation in dimension Template:Math. By induction, the proportionality relation is true in all dimensions.

The two-dimension recursion formula

A proof of the recursion formula relating the volume of the Template:Math-ball and an Template:Math-ball can be given using the proportionality formula above and integration in cylindrical coordinates. Fix a plane through the center of the ball. Let Template:Math denote the distance between a point in the plane and the center of the sphere, and let Template:Math denote the azimuth. Intersecting the Template:Math-ball with the Template:Math-dimensional plane defined by fixing a radius and an azimuth gives an Template:Math-ball of radius Template:Math. The volume of the ball can therefore be written as an iterated integral of the volumes of the Template:Math-balls over the possible radii and azimuths:

Vn(R)=02π0RVn2(R2r2)rdrdθ,

The azimuthal coordinate can be immediately integrated out. Applying the proportionality relation shows that the volume equals

Vn(R)=2πVn2(R)0R(1(rR)2)(n2)/2rdr.

The integral can be evaluated by making the substitution Template:Math to get

Vn(R)=2πVn2(R)[R2n(1(rR)2)n/2]r=0r=R=2πR2nVn2(R),

which is the two-dimension recursion formula.

The same technique can be used to give an inductive proof of the volume formula. The base cases of the induction are the 0-ball and the 1-ball, which can be checked directly using the facts Template:Math and Template:Math. The inductive step is similar to the above, but instead of applying proportionality to the volumes of the Template:Math-balls, the inductive hypothesis is applied instead.

The one-dimension recursion formula

The proportionality relation can also be used to prove the recursion formula relating the volumes of an Template:Math-ball and an Template:Math-ball. As in the proof of the proportionality formula, the volume of an Template:Math-ball can be written as an integral over the volumes of Template:Math-balls. Instead of making a substitution, however, the proportionality relation can be applied to the volumes of the Template:Math-balls in the integrand:

Vn(R)=Vn1(R)RR(1(xR)2)(n1)/2dx.

The integrand is an even function, so by symmetry the interval of integration can be restricted to Template:Math. On the interval Template:Math, it is possible to apply the substitution Template:Math. This transforms the expression into

Vn1(R)R01(1u)(n1)/2u12du

The integral is a value of a well-known special function called the beta function Template:Math, and the volume in terms of the beta function is

Vn(R)=Vn1(R)RB(n+12,12).

The beta function can be expressed in terms of the gamma function in much the same way that factorials are related to binomial coefficients. Applying this relationship gives

Vn(R)=Vn1(R)RΓ(n2+12)Γ(12)Γ(n2+1).

Using the value Template:Math gives the one-dimension recursion formula:

Vn(R)=RπΓ(n2+12)Γ(n2+1)Vn1(R).

As with the two-dimension recursive formula, the same technique can be used to give an inductive proof of the volume formula.

Direct integration in spherical coordinates

The volume of the n-ball Vn(R) can be computed by integrating the volume element in spherical coordinates. The spherical coordinate system has a radial coordinate Template:Math and angular coordinates Template:Math, where the domain of each Template:Math except Template:Math is Template:Math, and the domain of Template:Math is Template:Math. The spherical volume element is:

dV=rn1sinn2(φ1)sinn3(φ2)sin(φn2)drdφ1dφ2dφn1,

and the volume is the integral of this quantity over Template:Math between 0 and Template:Math and all possible angles:

Vn(R)=0R0π02πrn1sinn2(φ1)sin(φn2)dφn1dφ1dr.

Each of the factors in the integrand depends on only a single variable, and therefore the iterated integral can be written as a product of integrals:

Vn(R)=(0Rrn1dr)(0πsinn2(φ1)dφ1)(02πdφn1).

The integral over the radius is Template:Math. The intervals of integration on the angular coordinates can, by the symmetry of the sine about Template:Sfrac, be changed to Template:Math:

Vn(R)=Rnn(20π/2sinn2(φ1)dφ1)(40π/2dφn1).

Each of the remaining integrals is now a particular value of the beta function:

Vn(R)=RnnB(n12,12)B(n22,12)B(1,12)2B(12,12).

The beta functions can be rewritten in terms of gamma functions:

Vn(R)=RnnΓ(n212)Γ(12)Γ(n2)Γ(n21)Γ(12)Γ(n212)Γ(1)Γ(12)Γ(32)2Γ(12)Γ(12)Γ(1).

This product telescopes. Combining this with the values Template:Math and Template:Math and the functional equation Template:Math leads to

Vn(R)=2πn/2RnnΓ(n2)=πn/2RnΓ(n2+1).

Gaussian integrals

The volume formula can be proven directly using Gaussian integrals. Consider the function:

f(x1,,xn)=exp(12i=1nxi2).

This function is both rotationally invariant and a product of functions of one variable each. Using the fact that it is a product and the formula for the Gaussian integral gives:

𝐑nfdV=i=1n(exp(12xi2)dxi)=(2π)n/2,

where Template:Math is the Template:Math-dimensional volume element. Using rotational invariance, the same integral can be computed in spherical coordinates:

𝐑nfdV=0Sn1(r)exp(12r2)dAdr,

where Template:Math is an Template:Math-sphere of radius Template:Math (being the surface of an Template:Math-ball of radius Template:Math) and Template:Math is the area element (equivalently, the Template:Math-dimensional volume element). The surface area of the sphere satisfies a proportionality equation similar to the one for the volume of a ball: If Template:Math is the surface area of an Template:Math-sphere of radius Template:Math, then:

An1(r)=rn1An1(1).

Applying this to the above integral gives the expression

(2π)n/2=0Sn1(r)exp(12r2)dAdr=An1(1)0exp(12r2)rn1dr.

Substituting Template:Math:

0exp(12r2)rn1dr=2(n2)/20ett(n2)/2dt.

The integral on the right is the gamma function evaluated at Template:Math.

Combining the two results shows that

An1(1)=2πn/2Γ(n2).

To derive the volume of an Template:Math-ball of radius Template:Math from this formula, integrate the surface area of a sphere of radius Template:Math for Template:Math and apply the functional equation Template:Math:

Vn(R)=0R2πn/2Γ(n2)rn1dr=2πn/2nΓ(n2)Rn=πn/2Γ(n2+1)Rn.

Geometric proof

The relations Vn+1(R)=Rn+1An(R) and An+1(R)=(2πR)Vn(R) and thus the volumes of n-balls and areas of n-spheres can also be derived geometrically. As noted above, because a ball of radius R is obtained from a unit ball Bn by rescaling all directions in R times, Vn(R) is proportional to Rn, which implies dVn(R)dR=nRVn(R). Also, An1(R)=dVn(R)dR because a ball is a union of concentric spheres and increasing radius by ε corresponds to a shell of thickness ε. Thus, Vn(R)=RnAn1(R); equivalently, Vn+1(R)=Rn+1An(R).

An+1(R)=(2πR)Vn(R) follows from existence of a volume-preserving bijection between the unit sphere Sn+1 and S1×Bn:

(x,y,z)(xx2+y2,yx2+y2,z)

(z is an n-tuple; |(x,y,z)|=1; we are ignoring sets of measure 0). Volume is preserved because at each point, the difference from isometry is a stretching in the xy plane (in 1/x2+y2 times in the direction of constant x2+y2) that exactly matches the compression in the direction of the gradient of |z| on Sn (the relevant angles being equal). For S2, a similar argument was originally made by Archimedes in On the Sphere and Cylinder.

Balls in Template:Math norms

There are also explicit expressions for the volumes of balls in [[Lp space|Template:Math norm]]s. The Template:Math norm of the vector Template:Math in Template:Math is

xp=(i=1n|xi|p)1/p,

and an Template:Math ball is the set of all vectors whose Template:Math norm is less than or equal to a fixed number called the radius of the ball. The case Template:Math is the standard Euclidean distance function, but other values of Template:Math occur in diverse contexts such as information theory, coding theory, and dimensional regularization.

The volume of an Template:Math ball of radius Template:Math is

Vnp(R)=(2Γ(1p+1))nΓ(np+1)Rn.

These volumes satisfy recurrence relations similar to those for Template:Math:

Vnp(R)=(2Γ(1p+1))ppnRpVnpp(R)

and

Vnp(R)=2Γ(1p+1)Γ(n1p+1)Γ(np+1)RVn1p(R),

which can be written more concisely using a generalized binomial coefficient,

Vnp(R)=2(n/p1/p)RVn1p(R).

For Template:Math, one recovers the recurrence for the volume of a Euclidean ball because Template:Math.

For example, in the cases Template:Math (taxicab norm) and Template:Math (max norm), the volumes are:

Vn1(R)=2nn!Rn,Vn(R)=2nRn.

These agree with elementary calculations of the volumes of cross-polytopes and hypercubes.

Relation with surface area

For most values of Template:Math, the surface area An1p(R) of an Template:Math sphere of radius Template:Math (the boundary of an Template:Math Template:Math-ball of radius Template:Math) cannot be calculated by differentiating the volume of an Template:Math ball with respect to its radius. While the volume can be expressed as an integral over the surface areas using the coarea formula, the coarea formula contains a correction factor that accounts for how the Template:Math-norm varies from point to point. For Template:Math and Template:Math, this factor is one. However, if Template:Math then the correction factor is Template:Math: the surface area of an Template:Math sphere of radius Template:Math in Template:Math is Template:Math times the derivative of the volume of an Template:Math ball. This can be seen most simply by applying the divergence theorem to the vector field Template:Math to get

nVn1(R)=Template:Oiint = Template:Oiint =RnTemplate:Oiint =RnAn11(R).

For other values of Template:Math, the constant is a complicated integral.

Generalizations

The volume formula can be generalized even further. For positive real numbers Template:Math, define the Template:Math ball with limit Template:Math to be

Bp1,,pn(L)={x=(x1,,xn)𝐑n:|x1|p1++|xn|pnL}.

The volume of this ball has been known since the time of Dirichlet:[4]

V(Bp1,,pn(L))=2nΓ(1p1+1)Γ(1pn+1)Γ(1p1++1pn+1)L1p1++1pn.

Comparison to Template:Math norm

Using the harmonic mean p=n1p1+1pn and defining R=Lp, the similarity to the volume formula for the Template:Math ball becomes clear.

V({x𝐑n:|x1|p1++|xn|pnpR})=2nΓ(1p1+1)Γ(1pn+1)Γ(np+1)Rn.

See also

References

Template:Reflist

Further reading

  1. Equation 5.19.4, NIST Digital Library of Mathematical Functions. http://dlmf.nist.gov/5.19#E4, Release 1.0.6 of 2013-05-06.
  2. Smith, David J. and Vamanamurthy, Mavina K., "How Small Is a Unit Ball?", Mathematics Magazine, Volume 62, Issue 2, 1989, pp. 101–107, https://doi.org/10.1080/0025570X.1989.11977419.
  3. Template:Cite web
  4. Template:Cite journal