Volume of an n-ball

From testwiki
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 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 πR2≈3.142×R2 V1/2π≈0.564×V12
3 4π3R3≈4.189×R3 (3V4π)1/3≈0.620×V1/3
4 π22R4≈4.935×R4 (2V)1/4π≈0.671×V1/4
5 8π215R5≈5.264×R5 (15V8π2)1/5≈0.717×V1/5
6 π36R6≈5.168×R6 (6V)1/6π≈0.761×V1/6
7 16π3105R7≈4.725×R7 (105V16π3)1/7≈0.801×V1/7
8 π424R8≈4.059×R8 (24V)1/8π≈0.839×V1/8
9 32π4945R9≈3.299×R9 (945V32π4)1/9≈0.876×V1/9
10 π5120R10≈2.550×R10 (120V)1/10π≈0.911×V1/10
11 64π510395R11≈1.884×R11 (10395V64π5)1/11≈0.944×V1/11
12 π6720R12≈1.335×R12 (720V)1/12π≈0.976×V1/12
13 128π6135135R13≈0.911×R13 (135135V128π6)1/13≈1.007×V1/13
14 π75040R14≈0.599×R14 (5040V)1/14π≈1.037×V1/14
15 256π72027025R15≈0.381×R15 (2027025V256π7)1/15≈1.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×Vn−2(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)RVn−1(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(n−2))Rn−2(V),Rn(V)=Γ(n2+1)1/nΓ(n2+12)1/(n−1)V−1/(n(n−1))Rn−1(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)⋅(2k−1)⋯5⋅3⋅1. 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)⋅(2k−2)⋯4⋅2⋅2/π=2k⋅k!⋅2/π, the volume of an Template:Math-dimensional ball is, regardless of whether Template:Math is even or odd, Vn(R)=2(2π)(n−1)/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π)(n−1)/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]limn1nln⁡An−1(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: An−1(R)=ddRVn(R)=nRVn(R).

Thus, Template:Math inherits formulas and recursion relationships from Template:Math, such as An−1(R)=2πn/2Γ(n2)Rn−1. 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 the integral of volumes of Template:Math-balls: Vn(R)=∫−RRVn−1(R2−x2)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)=Rn−1∫−RRVn−1(1−(xR)2)dx. Making the change of variables Template:Math leads to: Vn(R)=Rn∫−11Vn−1(1−t2)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π∫0RVn−2(R2−r2)rdrdθ, The azimuthal coordinate can be immediately integrated out. Applying the proportionality relation shows that the volume equals Vn(R)=2πVn−2(R)∫0R(1−(rR)2)(n−2)/2rdr. The integral can be evaluated by making the substitution Template:Math to get Vn(R)=2πVn−2(R)⋅[−R2n(1−(rR)2)n/2]r=0r=R=2πR2nVn−2(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)=Vn−1(R)∫−RR(1−(xR)2)(n−1)/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 Vn−1(R)⋅R⋅∫01(1−u)(n−1)/2u−12du 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)=Vn−1(R)⋅R⋅B(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)=Vn−1(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)Vn−1(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=rn−1sinn−2(φ1)sinn−3(φ2)⋯sin⁡(φn−2)drdφ1dφ2⋯dφn−1, and the volume is the integral of this quantity over Template:Math between 0 and Template:Math and all possible angles: Vn(R)=∫0R∫0π⋯∫02πrn−1sinn−2(φ1)⋯sin⁡(φn−2)dφn−1⋯dφ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)=(∫0Rrn−1dr)(∫0πsinn−2(φ1)dφ1)⋯(∫02πdφn−1). 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(2∫0π/2sinn−2(φ1)dφ1)⋯(4∫0π/2dφn−1). Each of the remaining integrals is now a particular value of the beta function: Vn(R)=RnnB(n−12,12)B(n−22,12)⋯B(1,12)⋅2B(12,12).

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

Vn(R)=Rnn⋅Γ(n2−12)Γ(12)Γ(n2)⋅Γ(n2−1)Γ(12)Γ(n2−12)⋯Γ(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⁡(−12∑i=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=∫0∞∫Sn−1(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: An−1(r)=rn−1An−1(1). Applying this to the above integral gives the expression (2π)n/2=∫0∞∫Sn−1(r)exp⁡(−12r2)dAdr=An−1(1)∫0∞exp⁡(−12r2)rn−1dr. Substituting Template:Math: ∫0∞exp⁡(−12r2)rn−1dr=2(n−2)/2∫0∞e−tt(n−2)/2dt. The integral on the right is the gamma function evaluated at Template:Math.

Combining the two results shows that An−1(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)rn−1dr=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 that dVn(R)dR=nRVn(R). Also, An−1(R)=dVn(R)dR, since a ball is a union of concentric spheres and increasing the radius by ε sweeps out a shell of thickness ε. Thus, Vn(R)=RnAn−1(R); equivalently, Vn+1(R)=Rn+1An(R).

An+1(R)=(2πR)Vn(R), since a volume-preserving bijection between the unit sphere (Sn+1) and S1×Bn exists: (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 (by 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 (with 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 ‖x‖p=(∑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))ppnRpVn−pp(R)

and

Vnp(R)=2Γ(1p+1)Γ(n−1p+1)Γ(np+1)RVn−1p(R),

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

Vnp(R)=2(n/p1/p)RVn−1p(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 An−1p(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 =RnAn−11(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|pn≤L}. 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|pnp≤R})=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