Plücker coordinates: Difference between revisions

Sintel (talk | contribs)
make notation a bit more consistent
No edit summary
 
(17 intermediate revisions by 3 users not shown)
Line 1: Line 1:
{{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.]]
{{Wikipedia|Plücker embedding}}
{{Wikipedia|Plücker embedding}}
'''Plücker coordinates''' (also known as the '''[[wedgie]]'''), are a way to assign coordinates to 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 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 ==
We have a Grassmannian variety <math>\mathrm{Gr} (k, n)</math> consisting of the k-dimensional subspaces of <math>\mathbb{R}^n</math>.
A temperament can be viewed as a point in what is called a Grassmannian variety, written as <math>\mathrm{Gr} (k, n)</math>.
The rational points on this variety can be identified with rank-k temperaments on a JI space with n primes.
This variety contains all possible k-dimensional subspaces of <math>\mathbb{R}^n</math>.
In musical terms, k represents the rank of the temperament (how many independent generators it has), and n is the number of primes we're considering in our [[just intonation subgroup]].


Let <math>M</math> be an element of <math>\mathrm{Gr} (k, n)</math>, spanned by basis vectors <math>m_1, \ldots, m_k</math>.
Let <math>M</math> be an element of <math>\mathrm{Gr} (k, n)</math>, spanned by basis vectors <math>m_1, \ldots, m_k</math>.
We can embed the Grassmannian into a projective space using the Plücker map:
These basis vectors are the rows of the temperament mapping matrix.
$$
 
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 of our original space <math>\mathbb{R}^n</math>. The dimension of <math>\mathrm{Gr} (k, n)</math> is <math>k(n-k)</math>, while the dimension of <math>\Lambda^{k} \mathbb{R}^n</math> is <math>\binom{n}{k}</math>, which is typically much larger.
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>.


== Examples ==
== Examples ==
Line 25: 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 32: 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.
Put another way, we can always multiply by some <math>g \in GL_k (\mathbb{R})</math>.


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 45: 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 63: 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>.
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 ==
A height function is a way to measure the 'arithmetic complexity' of a rational point. For example, the rational numbers <math>\frac{3}{2}</math> and <math>\frac{3001}{2001}</math> are close to eachother, but intuitively the second is much more complicated.
A height function is a way to measure the 'arithmetic complexity' of a rational point. For example, the rational numbers <math>\frac{3}{2}</math> and <math>\frac{3001}{2001}</math> are close to eachother, but intuitively the second is much more complicated.


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]].


== Projective distance ==
== Projective distance ==
Given a temperament, we want to have some notion of distance, so that we can measure how well the temperament approximates JI. Since we are talking about linear subspaces (which all intersect at the origin), the only thing that is sensible to measure is the angle between them.
Given a temperament, we want to have some notion of distance, so that we can measure how well the temperament approximates JI. Since we are talking about linear subspaces (which all intersect at the origin), the only thing that is sensible to measure is the angle between them.


Line 101: 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 116: 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]]