Plücker coordinates: Difference between revisions

Sintel (talk | contribs)
No edit summary
No edit summary
 
(13 intermediate revisions by 3 users not shown)
Line 1: Line 1:
{{Expert}}
{{Expert|Wedgie}}
[[File:Plucker_embedding.png|thumb|600px|right|Schematic illustration of the Plücker embedding. Linear subspaces of <math>\mathbb{R}^n</math> (here lines) get mapped to points on a quadric surface in projective space.]]
[[File:Plucker_embedding.png|thumb|600px|right|Schematic illustration of the Plücker embedding. Linear subspaces of <math>\mathbb{R}^n</math> (here lines) get mapped to points on a quadric surface in projective space.]]
{{Wikipedia|Plücker embedding}}
{{Wikipedia|Plücker embedding}}


In [[exterior algebra]] applied to [[regular temperament theory]], '''Plücker coordinates''' (also known as the [[wedgie]]) are a way to assign coordinates to abstract temperaments, by viewing them as elements of some projective space.
In [[exterior algebra]] applied to [[regular temperament theory]], '''Plücker coordinates''' (also known as the '''wedgie''') are a way to assign coordinates to abstract temperaments, by viewing them as elements of some projective space.


The usual way to write down an abstract temperament is via its mapping matrix, but Plücker coordinates give us a unique description that is useful for some calculations.
The usual way to write down an abstract temperament is via its mapping matrix, but Plücker coordinates give us a unique description that is useful for some calculations.
The definition here is given in terms of temperament matrices, but by duality, we can also embed interval spaces in the same way.
More specifically, the interval subspace spanned by the commas of some temperament can also be used to give unique coordinates to that temperament.
These two representations are related via the [[Hodge dual]].


== Definition ==
== Definition ==
Line 16: Line 19:


The Plücker map takes a temperament and embeds it into a projective space by taking the wedge product of the basis vectors:
The Plücker map takes a temperament and embeds it into a projective space by taking the wedge product of the basis vectors:
$$
:<math>
\begin{align}
\begin{align}
\iota: \mathrm{Gr} (k, n)  
\iota: \mathrm{Gr} (k, n)  
& \to \mathbf{P}\left(\Lambda^{k} \mathbb{R}^n \right) \\
& \to \mathbf{P}\left(\Lambda^{k} \, \mathbb{R}^n \right) \\
\operatorname {span} (m_1, \ldots, m_k)  
\operatorname {span} (m_1, \ldots, m_k)  
& \mapsto \left[ m_1 \wedge \ldots \wedge m_k \right] \, .
& \mapsto \left[ m_1 \wedge \ldots \wedge m_k \right] \, .
\end{align}
\end{align}
$$
</math>


Here, <math>\Lambda^{k} \mathbb{R}^n</math> is the k-th exterior power (the subspace containing all k-vectors). This construction is independent of the basis we choose.
Here, <math>\Lambda^{k} \, \mathbb{R}^n</math> is the k-th exterior power (the subspace containing all k-vectors). This construction is independent of the basis we choose.
While the original space of temperaments has dimension <math>k(n-k)</math>, the space of Plücker coordinates is typically larger, with dimension <math>\binom{n}{k} - 1</math>.
While the original space of temperaments has dimension <math>k(n-k)</math>, the space of Plücker coordinates is typically larger, with dimension <math>\binom{n}{k} - 1</math>.


Line 34: Line 37:
The simplest non-trivial case is <math>\mathrm{Gr} (2, 4)</math>.
The simplest non-trivial case is <math>\mathrm{Gr} (2, 4)</math>.
An element <math>M</math> spanned by two lines <math>x, y</math>, can be represented as the matrix
An element <math>M</math> spanned by two lines <math>x, y</math>, can be represented as the matrix
$$
:<math>
\begin{equation}
\begin{equation}
\begin{bmatrix}
\begin{bmatrix}
Line 41: Line 44:
\end{bmatrix} \, .
\end{bmatrix} \, .
\end{equation}  
\end{equation}  
$$
</math>


These are not 'proper' coordinates, as doing row operations on this matrix preserves the row-span.
These are not 'proper' coordinates, as doing row operations on this matrix preserves the row-span.


The projective coordinates can be calculated by taking the determinants of all <math>2 \times 2</math> sub-matrices
The projective coordinates can be calculated by taking the determinants of all <math>2 \times 2</math> sub-matrices
$$
:<math>
p_{ij} =  
p_{ij} =  
\begin{vmatrix}
\begin{vmatrix}
Line 52: Line 55:
y_i & y_j
y_i & y_j
\end{vmatrix} \, ,
\end{vmatrix} \, ,
$$
</math>


which finally gives us
which finally gives us


$$
:<math>
\begin{equation}
\begin{equation}
\iota (M) = \left[ x \wedge y \right] = \left[ p_{12} : p_{13} : p_{14} : p_{23} : p_{24} : p_{34} \right] \, .
\iota (M) = \left[ x \wedge y \right] = \left[ p_{12} : p_{13} : p_{14} : p_{23} : p_{24} : p_{34} \right] \, .
\end{equation}
\end{equation}
$$
</math>


Note the use of colons to signify that these coordinates are homogeneous.
Note the use of colons to signify that these coordinates are homogeneous.
Line 70: Line 73:
For the example above on <math>\mathrm{Gr} (2, 4)</math>, the Plücker relation is
For the example above on <math>\mathrm{Gr} (2, 4)</math>, the Plücker relation is


$$
:<math>
p_{12} p_{34} - p_{13} p_{24} + p_{14} p_{23} = 0 \, .
p_{12} p_{34} - p_{13} p_{24} + p_{14} p_{23} = 0 \, .
$$
</math>


Note that in this case, there is only one such relation, but in higher dimensions there will be many.
Note that in this case, there is only one such relation, but in higher dimensions there will be many.


== Rational points ==
== Rational points ==
A rational point <math>P</math> on <math>\mathrm{Gr}(k, n)</math> is a k-dimensional subspace such that <math>P \cap \mathbb{Z}^n</math> is a rank k sublattice of <math>\mathbb{Z}^n</math>. Abstract temperaments correspond exactly to these rational points (although the vast majority of them will be terrible temperaments).
A '''rational point''' <math>P</math> on <math>\mathrm{Gr}(k, n)</math> is a k-dimensional subspace such that <math>\mathcal{L} = P \cap \mathbb{Z}^n</math> is a rank k sublattice of <math>\mathbb{Z}^n</math>. Abstract temperaments correspond exactly to these rational points, although most have no practical musical use.


The same relations as above can be derived, where we represent P as integer matrix <math>M \in \mathbb{Z} ^ {k \times n}</math> and the projective coordinates similarly have entries in <math>\mathbb{Z}</math> instead.
The same relations as above can be derived, where we represent P as integer matrix <math>M \in \mathbb{Z} ^ {k \times n}</math>, whose rows span <math>\mathcal{L}</math>.
Because the Plücker coordinates are homogeneous, we can always put them in a 'canonical' form by dividing all entries by their GCD and ensuring the first element is non-negative.
The projective coordinates similarly have integer entries.
Because the Plücker coordinates are homogeneous, we can always put them in a canonical form by dividing all entries by their greatest common divisor (GCD) and ensuring the first element is non-negative.
 
An advantage of studying rational points is that we do not have to worry about [[torsion]].
The quotient group <math>\mathbb{Z}^n / \mathcal{L}</math> is a finitely generated abelian group. 
When the Plücker coordinates are normalized (GCD = 1), we ensure that
<math>
\mathbb{Z}^n / \mathcal{L} \cong \mathbb{Z}^{n-k},
</math>
which is torsion-free.


== Height ==
== Height ==
Line 86: Line 98:


We can define the height of a rational point simply as the Euclidean norm on its Plücker coordinates <math>X = \iota (P)</math>.
We can define the height of a rational point simply as the Euclidean norm on its Plücker coordinates <math>X = \iota (P)</math>.
$$
:<math>
H(P) = \left\| X \right\| = \left\| p_1 \wedge \ldots \wedge p_n \right\| \\
H(P) = \left\| X \right\| = \left\| m_1 \wedge \ldots \wedge m_n \right\| \\
$$
</math>


This is equivalent to the covolume of the lattice defined by P (also know as the lattice determinant), which can be easily computed using the Gram matrix.
In terms of the lattice defined by P, this definition is equivalent to the volume of the {{w|fundamental domain}}, also known as the lattice determinant.
$$
It is easy to show that this does not depend on the basis we choose.
 
The height can be easily computed using the {{w|Gram matrix}}:
:<math>
\begin{align}
\begin{align}
\mathrm{G}_{ij} &= \left\langle p_i, p_j \right\rangle \\
\mathrm{G}_{ij} &= \left\langle m_i, m_j \right\rangle \\
\sqrt{\det(\mathrm{G})} &= \left\| p_1 \wedge \ldots \wedge p_n \right\| = \left\| X \right\| \, .
\sqrt{\det(\mathrm{G})} &= \left\| m_1 \wedge \ldots \wedge m_n \right\| = \left\| X \right\| \, .
\end{align}
\end{align}
$$
</math>


In regular temperament theory, this height is usually known as simply the [[Tenney-Euclidean temperament_measures #TE complexity|complexity]].
In regular temperament theory, this height is usually known as simply the [[Tenney-Euclidean temperament_measures #TE complexity|complexity]].
Line 106: Line 121:
Given vectors <math>a, b \in \mathbb{R^n}</math>, we famously have
Given vectors <math>a, b \in \mathbb{R^n}</math>, we famously have


$$
:<math>
\frac{a \cdot b}{\left\| a \right\|  \left\| b \right\| } = \cos (\theta) \, .
\frac{a \cdot b}{\left\| a \right\|  \left\| b \right\| } = \cos (\theta) \, .
$$
</math>


In projective space, there is an analogous formula, using the wedge product instead.
In projective space, there is an analogous formula, using the wedge product instead.
Given some real point <math>j \in \mathbb{R^n}</math> with homogeneous coordinates <math>y</math>, and a linear subspace <math>P \in \mathrm{Gr} (k, n)</math> with Plücker coordinates <math>X</math>, we define the projective distance as
Given some real point <math>j \in \mathbb{R^n}</math> with homogeneous coordinates <math>y</math>, and a linear subspace <math>P \in \mathrm{Gr} (k, n)</math> with Plücker coordinates <math>X</math>, we define the projective distance as


$$
:<math>
d(P, j) = \frac{ \left\| X \wedge y \right\| }{\left\| X \right\|  \left\| y \right\| } = \sin (\theta) \, .
d(P, j) = \frac{ \left\| X \wedge y \right\| }{\left\| X \right\|  \left\| y \right\| } = \sin (\theta) \, .
$$
</math>


Where we can take <math>j</math> to be the usual n-limit vector of log primes, so that <math> y = \left[ 1 : \log_2 (3) : \ldots : \log_2 (p_n) \right] </math>.
Where we can take <math>j</math> to be the usual n-limit vector of log primes, so that <math> y = \left[ 1 : \log_2 (3) : \ldots : \log_2 (p_n) \right] </math>.
Line 121: Line 136:


Since for any decent temperament this angle will be extremely small, we can take <math>\sin (\theta) \approx \theta</math>.
Since for any decent temperament this angle will be extremely small, we can take <math>\sin (\theta) \approx \theta</math>.
== See also ==
* [[Wedgie supplement]] - Supplementary page going over additional information on wedgies
* [[Exterior algebra]] - exterior product, which produces wedgies
* [[Interior product]] - interior product, dual of the exterior product
* [[Hodge dual]] - acts on wedgies
[[Category:Exterior algebra]]