Tenney–Euclidean tuning: Difference between revisions
→Otherwise normed tunings: correction (benedetti height isn't a norm, recte wilson height) |
Unify nomenclature (k -> r, A -> V, B -> M); also the tuning map shouldn't be weighted |
||
| Line 1: | Line 1: | ||
'''Tenney-Euclidean tuning''' ('''TE tuning'''), also known as '''TOP-RMS tuning''', is a tuning technique for regular temperaments which leads to the least sum of squared errors of the Tenney-weighted basis. | '''Tenney-Euclidean tuning''' ('''TE tuning'''), also known as '''TOP-RMS tuning''', is a tuning technique for regular temperaments which leads to the least sum of squared errors of the Tenney-weighted basis. | ||
If we have '' | If we have ''r'' linearly independent [[Vals and tuning space|vals]] of dimension ''n'', they will span a subspace of [[Vals and tuning space|tuning space]]. This subspace defines a regular temperament of rank ''r'' in the prime limit ''p'', where ''p'' is the ''n''-th prime. Similarly, starting from ''n'' - ''r'' independent commas for the same regular temperament, the corresponding monzos span an ''n'' - ''r'' dimensional subspace of [[Monzos and interval space|interval space]]. Both the subspace of tuning space and the subspace of interval space characterize the temperament completely. A question then arises as to how to choose a specific tuning for this temperament, which is the same as asking how to choose a point (vector) in this subspace of tuning space which provides a good tuning. One answer to this is the weighted RMS ([[Wikipedia: Root mean square|root mean squared]]) tuning discussed right here. | ||
TE tuning can be viewed as a variant of [[TOP tuning]] since it employs the [[Tenney-Euclidean metrics #TE norm|TE norm]] in place of the [[Tenney height]] as in TOP tuning. Just as TOP tuning minimizes the maximum Tenney-weighted ''L''<sub>1</sub> error of any interval, TE tuning minimizes the maximum Tenney-weighted ''L''<sub>2</sub> error of any interval. | TE tuning can be viewed as a variant of [[TOP tuning]] since it employs the [[Tenney-Euclidean metrics #TE norm|TE norm]] in place of the [[Tenney height]] as in TOP tuning. Just as TOP tuning minimizes the maximum Tenney-weighted ''L''<sub>1</sub> error of any interval, TE tuning minimizes the maximum Tenney-weighted ''L''<sub>2</sub> error of any interval. | ||
| Line 12: | Line 12: | ||
<math>W = \operatorname {diag} (1/\log_2 (Q))</math> | <math>W = \operatorname {diag} (1/\log_2 (Q))</math> | ||
If | If V is the mapping of the [[Regular temperament|abstract temperament]] whose rows are (not necessarily independent) vals, then V<sub>W</sub> = VW is the mapping in the weighted space. If J is the row vector of targeted JI intervals (i.e. the [[JIP]]), then J<sub>W</sub> = JW is the JI intervals in the weighted space, in the case of Tenney-weighting it is {{val| 1 1 … 1 }}. Let us also denote the row vector of TE generators G. TE tuning then defines a [[Wikipedia: Least squares|least squares problem]] of the following overdetermined linear equation system: | ||
<math> | <math> \displaystyle | ||
GV_W = J_W | |||
</math> | |||
The system simply says that the sum of ''v''<sub>''kl''</sub> steps of generator ''g''<sub>''k''</sub> for all ''k'''s should equal the ''l''-th targeted JI interval ''j''<sub>''l''</sub>. | The system simply says that the sum of (''v''<sub>''w''</sub>)<sub>''kl''</sub> steps of generator ''g''<sub>''k''</sub> for all ''k'''s should equal the ''l''-th targeted JI interval (''j''<sub>''w''</sub>)<sub>''l''</sub>. | ||
There are a number of methods to solve least squares problems. One common way is to use the [[Wikipedia: Moore–Penrose pseudoinverse|Moore–Penrose pseudoinverse]]. | There are a number of methods to solve least squares problems. One common way is to use the [[Wikipedia: Moore–Penrose pseudoinverse|Moore–Penrose pseudoinverse]]. | ||
| Line 35: | Line 37: | ||
In the pseudoinverse method, the (not necessarily independent) TE generators which correspond to the rows of V are given by | In the pseudoinverse method, the (not necessarily independent) TE generators which correspond to the rows of V are given by | ||
<math>G = | <math>\displaystyle | ||
G = J_W V_W^+ | |||
</math> | |||
Applying the weighted val list to the generators, The TE tuning map is given by | Applying the weighted val list to the generators, The TE tuning map is given by | ||
<math>T = GV = | <math> | ||
\displaystyle T = GV = J_W V_W^+ V | |||
</math> | |||
We may also obtain the TE tuning from a projection matrix. P = V<sup>+</sup>V is the orthogonal projection matrix that maps onto the space spanned by the rows of V. This space corresponds to the temperament, and so does P. However, P is independent of how the temperament is defined; it does not depend on whether the vals are linearly independent, how many of them there are, or whether [[contorsion]] has been removed. The tuning map giving the tuning of each prime number is found by multiplying by the | We may also obtain the TE tuning from a projection matrix. P = V<sub>W</sub><sup>+</sup>V<sub>W</sub> is the orthogonal projection matrix that maps onto the space spanned by the rows of V. This space corresponds to the temperament, and so does P. However, P is independent of how the temperament is defined; it does not depend on whether the vals are linearly independent, how many of them there are, or whether [[contorsion]] has been removed. The tuning map giving the tuning of each prime number is found by multiplying by the JIP: T = JP where J is the JIP, which is the nearest point in the subspace corresponding to the temperament to J. | ||
We may find the same projection matrix starting from a list of weighted monzos rather than vals. If M is a rank ''n'' matrix whose columns are weighted monzos, and I is the ''n''×''n'' identity matrix, then P = I - | We may find the same projection matrix starting from a list of weighted monzos rather than vals. If M<sub>W</sub> is a rank-''n'' matrix whose columns are weighted monzos, and I is the ''n''×''n'' identity matrix, then P = I - M<sub>W</sub>M<sub>W</sub><sup>+</sup> is the same projection matrix as V<sub>W</sub><sup>+</sup>V<sub>W</sub> so long as the temperament defined by the vals is the same as the temperament defined by the monzos. Again, it is irrelevant if the monzos are independent or how many of them there are. | ||
== Enforcement == | == Enforcement == | ||
=== Pure-octaves TE tuning === | === Pure-octaves TE tuning === | ||
{{ | {{Main| POTE tuning }} | ||
We may call pure-octaves Tenney-Euclidean tuning the | We may call pure-octaves Tenney-Euclidean tuning the ''POTE tuning''. If T = JP = GV is the TE tuning map, then a corresponding pure-octaves map can be found by [[wikipedia: Scalar multiplication|scalar multiplication]], T/T<sub>1</sub>, where T<sub>1</sub>, the first entry, is the tuning of 2. The justification for this is that T does not only define a point, but a line through the origin lying in the subspace defining the temperament, or in other words, a point in the linear subspace of projective space corresponding to the temperament, and hence is a projective object. Another way to say this is that T defines not only the closest point to J, but the closest direction in terms of angular measure between the line through T and the line through J. | ||
=== Constrained TE tuning === | === Constrained TE tuning === | ||
{{Main| | {{Main| Constrained tuning }} | ||
Another way to enforce the pure octave is by adding the constraint before the optimization process. This is the | Another way to enforce the pure octave is by adding the constraint before the optimization process. This is the ''CTE tuning''. The result, under the constraint of pure octaves, remains TE optimal. | ||
== Otherwise normed tunings == | == Otherwise normed tunings == | ||
| Line 60: | Line 66: | ||
We may also do the same things starting from nonweighted vals. This leads to a different tuning, the '''Frobenius tuning''', which is perfectly functional but has less theoretical justification than TE tuning. However, if greater weight needs to be attached to the larger primes than TE tuning attaches, Frobenius tuning may be preferred; people who feel that larger primes require more tuning care than smaller ones may well prefer it. | We may also do the same things starting from nonweighted vals. This leads to a different tuning, the '''Frobenius tuning''', which is perfectly functional but has less theoretical justification than TE tuning. However, if greater weight needs to be attached to the larger primes than TE tuning attaches, Frobenius tuning may be preferred; people who feel that larger primes require more tuning care than smaller ones may well prefer it. | ||
The list of Frobenius generators, G | The list of Frobenius generators, G, is given by: | ||
<math> | <math>G = JV^+</math> | ||
where J | where J is the nonweighted JIP and V is the nonweighted mapping. | ||
The Frobenius tuning map, T | The Frobenius tuning map, T, is given by: | ||
<math> | <math>T = GV = JV^+V</math> | ||
However, the main value of unweighted vals is that the pseudoinverse and projection matrix have rational entries, so that the rows of the matrix are [[fractional monzos]]. The Frobenius projection matrix therefore, like the [[wedgie]], defines a completely canonical object not depending on any arbitrary definition (e.g. how Hermite normal form or LLL reduction is specifically defined) which corresponds one-to-one with temperaments, and which does not depend on whether the monzos or vals from which it is computed are [[Mathematical theory of saturation|saturated]]. It may be found starting either from a set of vals or a set of commas, since if Q is the projection matrix found by treating monzos in the same way as vals, P = I - Q is the same projection matrix as would be found if starting from a set of vals defining the same temperament. | However, the main value of unweighted vals is that the pseudoinverse and projection matrix have rational entries, so that the rows of the matrix are [[fractional monzos]]. The Frobenius projection matrix therefore, like the [[wedgie]], defines a completely canonical object not depending on any arbitrary definition (e.g. how Hermite normal form or LLL reduction is specifically defined) which corresponds one-to-one with temperaments, and which does not depend on whether the monzos or vals from which it is computed are [[Mathematical theory of saturation|saturated]]. It may be found starting either from a set of vals or a set of commas, since if Q is the projection matrix found by treating monzos in the same way as vals, P = I - Q is the same projection matrix as would be found if starting from a set of vals defining the same temperament. | ||
Spelling this out, if | Spelling this out, if V is a matrix whose rows are vals, then P = V<sup>+</sup>V is a [[wikipedia: Positive-definite matrix|positive-semidefinite]] [[wikipedia: Symmetric matrix|symmetric matrix]] with rational matrix entries, which exactly specifies the regular temperament defined by the vals of V. If M is a matrix with columns of monzos which spans the subspace of interval space containing the commas, then this same matrix P is given by I - MM<sup>+</sup>. | ||
If the vals defining | If the vals defining V are linearly independent, then P = V<sup>T</sup>(VV<sup>T</sup>)<sup>-1</sup>V. If the columns of M are independent, then we likewise have P = I - M(M<sup>T</sup>M)<sup>-1</sup>M<sup>T</sup>. | ||
=== Benedetti-Euclidean tuning === | === Benedetti-Euclidean tuning === | ||