Sheng Li (gmachine1729) wrote,
Sheng Li

How to compute the volume of an arbitrary parallelotope embedded in R^n

Originally published at 狗和留美者不得入内. You can comment here or there.

[公式] linearly independent vectors in [公式] result in a [公式] -dimensional parallelotope (a generalization of parallelogram or parallelpiped in higher dimensions). We wish to determine its generalized volume with orientation in [公式] dimensions. The determinant of an [公式] matrix gives the (signed) volume of the parallelotope induced by the column vectors, the order of which affects the sign. (The proof of this is rather straightforward and will be left to the reader.) Independent of coordinates, a linear transformation [公式] between two [公式] -dimensional vector spaces [公式] the positively oriented orthogonal bases of which are [公式] respectively assigns each element in [公式] to a linear combination of the elements of [公式] , and the determinant is the volume of the parallelotope generated by elements of [公式] [公式] .

Let [公式] be the vectors for our [公式] -dimensional parallelotope in [公式] . Together they define a map [公式] from [公式] to [公式] that by the rank-nullity theorem is onto, with [公式] . Let [公式] be a positively oriented orthogonal basis of [公式] , the subspace of [公式] perpendicular to the kernel of [公式] with [公式] . We then define the same map on a restricted domain as

[公式]We then have


The coordinate free definition adjoint operator (or transpose) and its determinant

If [公式] is a linear map from [公式] and inner products are defined on [公式] -dimensional vector space [公式] with respect to their bases such that

[公式] Then the adjoint of [公式] , which we denote as [公式] is the map from [公式] such that

[公式] In the language of transposes, we have

[公式] Note how the left hand side of [公式] gives for every element in an ordered set of vectors [公式] in [公式] with respect to an orthonormal basis of [公式] , and then prescribes an ordered set of vectors [公式] in [公式] the coordinates of the [公式] th of which, with respect to an orthonormal basis of [公式] , is the [公式] th coordinate of the elements of [公式] , and vice versa.

We’ve essentially defined an [公式] matrix the elements of which are [公式] , and then the elements of transpose matrix [公式] . Applying the permutation based determinant formula gives us


Computing the volume of the parallelotope

The matrix formed by linearly independent [公式] column vectors in [公式] , [公式] , corresponds directly to a linear isomorphism from [公式] , the codomain of which is of course a [公式] dimensional subspace of [公式] . We showed in the previous section that the adjoint of [公式] , namely [公式] , has the same determinant. Since determinant is a multiplicative function, [公式] thus gives us the square of the volume of the parallelotope.

Per the rule of matrix multiplication, given a linear isomorphism [公式] which takes the elements of the orthonormal basis of its domain to vectors [公式] , regardless of basis or coordinates in the range, its composition with its adjoint is represented with respect to aforementioned basis by

[公式] wherein the invoked inner product is, of course, per the properties of inner product invariant with respect to coordinate transformations, which means [公式] is well-defined. This corresponds to the matrix formed from our [公式] vectors in [公式] ,


mutiplied by its transpose on the left, which gives us the result prescribed by [公式] , wherein the inner product is the Euclidean inner product in [公式] restricted to a [公式] dimensional subspace within it. The resulting Gramian matrix


is, as we’ve already explained, such that

[公式] or in words the square of the volume of the parallelotope generated by [公式] .

Decomposing a [公式]volume element into orthogonal components in [公式] -form space

There is also intimate connection here with differential forms, exterior products, and Hodge dual. In [公式] of [1], we defined an inner product on the [公式] th exterior product space with the Gramian determinant invoked in [公式] such that

[公式] Let [公式] be an orthonormal positively oriented basis for our vector space, which by definition results in the equivalence relation

[公式] as far as [公式] -dimensional volume is concerned, with which each of our [公式] s decomposes to [公式] . Substitution of this decomposition into [公式] yields for the coefficient of [公式] the determinant of the matrix assembled from the [公式] rows of our [公式] column vectors, with an appropriate sign adjustment. This coefficient is necessarily also an anti-symmetric [公式] th rank tensor.

Geometrically this is the oriented [公式] -dimensional volume obtained if only the components corresponding to indices [公式] are considered. The result in [公式] of [1], which was calculated via the Hodge star in a way that equates to the Gramian matrix definition of the inner product, yields for the value of the inner product given in [公式] the sum of squares of the coefficients with respect to our [公式] basis elements [公式] . Essentially there is a basis of [公式] of dimension [公式] volume elements of equivalent volume in [公式] space represented by [公式] -forms, which are mutually orthogonal with respect to the inner product we defined on the space of [公式] -forms, and we have projected our arbitrary [公式] dimensional volume element onto each of them. In this sense it is natural that the square of its norm or size would be the sum of the squares of its components.

I am dedicating this article to Seki Takakazu (1642-1708), who based on no more than 13th century Chinese algebra and arithmetic obtained results regarding determinants and resultants decades before the West and who discovered Bernoulli numbers (or Takakazu numbers) in connection to the closed formula for sum of the first [公式] th powers around the same time as did Jacob Bernoulli.


Tags: uncategorized

  • Post a new comment


    Anonymous comments are disabled in this journal

    default userpic

    Your reply will be screened