A Compendium of Musical Mathematics
A Compendium of Musical Mathematics
Musical
Mathematics
This page intentionally left blank
A Compendium of
Musical
Mathematics
Franck Jedrzejewski
CEA, Paris-Saclay University, France
World Scientific
NEW JERSEY • LONDON • SINGAPORE • BEIJING • SHANGHAI • HONG KONG • TAIPEI • CHENNAI • TOKYO
Published by
World Scientific Publishing Co. Pte. Ltd.
5 Toh Tuck Link, Singapore 596224
USA office: 27 Warren Street, Suite 401-402, Hackensack, NJ 07601
UK office: 57 Shelton Street, Covent Garden, London WC2H 9HE
Printed in Singapore
About the Author
v
This page intentionally left blank
Contents
vii
viii A Compendium of Musical Mathematics
xi
xii A Compendium of Musical Mathematics
References
1
2 A Compendium of Musical Mathematics
∃ n ∈ ZN \{0}, Tn (A) = A
HA = {n ∈ ZN , Tn (A) = A}
is not reduce to {0}. Let L be the set of the sets of limited transpo-
sition.
For enumeration purpose, we consider class of sets of limited
transposition.
Definition 1.12. Let L be the quotient set L = L/∼ of the set class
of limited transposition by the cyclic group, up to transpositions.
A generator set is a set of elements of the quotient set L contain-
ing {0}.
This choice is well defined because if A is a set of limited trans-
position (A ∈ L), then all its transpositions belong to L:
∀m ∈ ZN , Tmn (a) ∈ A
Tritone (0, 6) 6
Aug. Fifth (0, 4, 8) 4
Dim. Seventh (0, 3, 6, 9) 3
Alt. Seventh (0, 2, 6, 8) 6
Double Fifths (0, 1, 6, 7) 6
Mode 1 (0, 2, 4, 6, 8, 10) 2
Mode 2 (0, 1, 3, 4, 6, 7, 9, 10) 3
Mode 3 (0, 1, 2, 4, 5, 6, 8, 9, 10) 4
Mode 4 (0, 1, 2, 3, 6, 7, 8, 9) 6
Mode 5 (0, 1, 2, 6, 7, 8) 6
Mode 6 (0, 1, 2, 4, 6, 7, 8, 10) 6
Mode 7 (0, 1, 2, 3, 4, 6, 7, 8, 9, 10) 6
Mode A (0, 1, 4, 5, 8, 9) 4
Mode B (0, 1, 3, 6, 7, 9) 6
Mode C (0, 1, 4, 6, 7, 10) 6
Chromatic Scale (0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11) 1
8 A Compendium of Musical Mathematics
The set
The set
The sets
The set E8 ∪T1 (E8 )∪T4 (E8 ) = E14 = {0, 1, 4, 6, 7, 10} is the mode C.
The following sets are transposition of the preceding set:
11
12 A Compendium of Musical Mathematics
1 r−k+2
r r−k+1 r N
KN = (−1)k+1 ··· 2 pi1 pi2 ···pik
2
k=1 i1 =1 i2 =i1 +1 ik =ik−1 +1
Musical Set Theories 13
N
N
PN = PN,k and MN = MN,k
k=1 k=1
1
PN,k − LN,k = (MN,k − KN,k )
k
As the number of pcset classes of k notes and the number of pcset
classes of N − k notes are equal, the sum over k leads to
2
PN − LN = (MN − KN )
N
From this relation, we get the expression of LN . On the other hand,
the number of modes of k notes is
N −1 (N − 1)!
MN,k = =
k−1 (N − k)!(k − 1)!
N
N −1
N −1
MN = MN,k = = 2N −1
k
k=1 k=0
fm : MN/m → KN
r
r r−k+2
r−k+1
r
Aj = k+1
(−1) ··· Ai1 ∩ Ai2 ∩ · · · ∩ Aik
j=1 k=1 i1 =1 i2 =i1 +1 ik =ik−1 +1
we get
r r−k+2
r−k+1
KN = |KN | = (−1)k+1 ···
k=1 i1 =1 i2 =i1 +1
r
fpi1 ··· pik (MN/pi1 ··· pi )
k
ik =ik−1 +1
and then
N
fpi1 ..pik (MN/pi1 ..pi ) = MN/pi1 ..pi = 2 p i1 p i2 · · · p ik
k k
Theorem 1.7. Let N be the number of notes per octave in the equal
temperament and k an integer less than N . Let’s denote m the num-
ber of common prime factors of N and k, and q1 , q2 , . . . , qm these
factors. The number of modes MN,k is given by
m m−k+2
m−k+1
m N
− 1
qi1 qi2 ···qik
MN,k = (−1)k+1
···
k=1 i1 =1 i2 =i1 +1 ik = ik−1 +1
k
qi1 qi2 ···qik −1
The number of pcset classes (Pcsc) in the N -EDO and the number
of limited transposition set classes (Ltsc) are given in the following
table for different values of N (in the cyclic classification).
For divisions of m-th of tone (m > 5, m = N/6), we get large
values.
G/H = {gh | h ∈ H}
G(h) = {ghg−1 | g ∈ G}
Gh = {g ∈ G, gh = hg}
|G/H| = |H\G|
Proof.
1 1 |G| 1
|G\\X| = = = |Gx |
|G(x)| |G| |G(x)| |G|
x∈X x∈X x∈X
1 1 1
= 1= 1= |Xg |
|G| |G| |G|
x∈X g∈Gx g∈G x∈Xg g∈G
In order to introduce Pólya’s enumeration theorem, we need a
generalized version of this lemma. Let R be a commutative ring such
that Q is a subring of R and G be a finite group acting on a finite
set X. A weight function is a function w : X → R constant on each
G-orbit, i.e.
w(gx) = w(x) ∀g ∈ G, ∀x ∈ X
Proof.
w(u) w(u)
w(u) = =
|u| |G(x)|
u∈G\\X u∈G\\X x∈u u∈G\\X x∈u
1 |G| 1
= w(x) = |Gx | w(x)
|G| |G(x)| |G|
x∈X x∈X
1 1
= w(x) = w(x)
|G| |G|
x∈X g∈Gx g∈G x∈Xg
In the following, X, Y are two finite sets, R is a commutative
ring containing Q as a subring and G is a group acting on X.
The set of configurations Y X is the set of functions from X to Y.
Musical Set Theories 21
(g, f ) → f ◦ g −1
1
|X|
j (g)
P(G,X) (t1 , . . . , t|X| ) = tkk
|G|
g∈G k=1
⎛ ⎞jk (g)
1
|X|
ω(u) = ⎝ h(y)k ⎠
X
|G|
u∈G\\Y g∈G k=1 y∈Y
1
|X|
j (g c )
P(G,X) (t1 , . . . , t|X| ) = |c| tkk
|G|
c∈C k=1
where C is the set of all conjugacy classes. If for the weights of the
elements of Y we take power of an independent variable z, the power
series in Pólya’s formula are called the configuration counting series
and the figures counting series. The coefficient of order k of the first
series is the number of classes in Y X of weight z k , and the coefficients
of the second series are the numbers of elements of Y of weight z k .
j k
1 tk
P(SN ,X) =
jk ! k
j k
N
kjk = N
k=1
tk = 1 + z k
1
P(A2 ,Z2 ) (t1 , t2 ) = (t21 + t2 )
2
1
P(A4 ,Z4 ) (t1 , t2 , t3 , t4 ) = (t41 + 2t21 t2 + 3t22 + 2t4 )
8
and for a ≥ 3
⎛
a−1
1 ⎝22(a−1)−1 t2a + a−i
P(A2a ,Z2a ) = 22(i−1) + ϕ(2i−1 )2a−1 t22i
22a−1
i=1
2a−i−1⎞
a−2 a−i a−i−1 −1
i
+ ϕ(2i ) 2i t21 + 2a−1 t21 t22 t2k ⎠
i=0 k=1
N m
P (t1 , t2 , . . . , tN ) = a(i) tikk , Q(t1 , t2 , . . . , tm ) = b(j) tjkk
(i) k=1 (j) k=1
N m
i j gcd(k,l)
P ∗Q= a(i) b(j) k l
tlcm(k,l)
(i) (j) k=1l=1
r
P(AN ,ZN ) = ∗ P
j=1 A aj ,Z aj
p p
j j
1 3
P(A3 ,Z3 ) (t1 , t2 , t3 ) = (t + 3t1 t2 + 2t3 )
6 1
1.5 Exercises
Exercise 1.1
Let us study the sets of limited transposition for N = 6. By imitating
the proof of Theorem 1.4, find the four sets of limited transposition.
Exercise 1.2
In the N = 6 universe, using the results of the previous exercise,
explain the enumeration of chords as it is done in Section 1.4. Find
these results by computing the generator polynomial (1.1). Establish
the table of the number of k-chords up to transposition and the
number of sets of limited transposition.
Musical Set Theories 29
Exercise 1.3
Let us study the sets of limited transposition for N = 8. By imitating
the proof of Theorem 1.4, find the five sets of limited transposition.
Exercise 1.4
In the N = 8 universe, using the results of the previous exercise,
explain the enumeration of chords as it is done on Section 1.4. Find
these results by computing the generator polynomial (1.1). Establish
the table of the number of k-chords up to transposition and the
number of sets of limited transposition.
Exercise 1.5
Write down the generator polynomial (1.2) for the action of the dihe-
dral group for N = 6.
Exercise 1.6
Write down the generator polynomial (1.2) for the action of the dihe-
dral group for N = 8.
Exercise 1.7
In the N = 8 universe, we consider the group T, I, M3 generated by
the transpositions, the inversions and the multiplication by 3 mod 8.
Show that the inverse of M3 is M5 the multiplication by 5 mod 8 and
write down the associated permutation. Compute the order of this
group.
References
Loquin, A. (1874). Tableau de tous les effets harmoniques d’une à cinq notes
inclusivement, Féret, Bordeaux.
Loquin, A. (1884). Une révolution dans la science des accords, algèbre de
l’harmonie, traité complet d’harmonie moderne, sans notes de musiques ni
signes équivalents, avec des nombres pour représenter les effets harmoniques
et des lettres pour exprimer les mouvements de basse, Richault, Paris.
Mazzola, G. (1985). Gruppen und Kategorien in der Musik. Entwurf einer math-
ematischen Musiktheorie, Heldermann Verlag, Berlin.
Messiaen, O. (1944). Techniques de mon langage musical, Leduc, Paris.
Messiaen, O. (1996–2002). Traité du rythme, de couleur et d’ornithologie, Vol. 7,
Alphonse Leduc, Paris.
Morris, R. (1991). Class Notes for Atonal Music Theory, Frog Peak Music,
Hanover.
Parzysz, B. (1984). Musique et mathématique, Publication de l’Association des
professeurs de mathématiques de l’enseignement public, Paris.
Piňos, A. (1971). Tónové skupiny, Edition Supraphon, Prague. English transla-
tion, Tone Groups, Janáčkova akademie múzických uměnı́, Brno, 2001.
Starr, D. (1978). Sets, invariance and partitions, Journal of Music Theory 22(1),
1–42.
Wei, W.-D. and Xu, J.-Y. (1993). Cycle index of direct product of permutation
groups and number of equivalence classes of subsets of Zv , Discrete Math-
ematics 123, 179–188.
Chapter 2
31
32 A Compendium of Musical Mathematics
Example 2.1.
(1) The triple (Z, G, int) where G is the additive group of integers
(Z, +) is a GIS for the function int(x, y) = y − x.
(2) The triple (Z12 , G, int) where G is the additive group of integers
modulo 12 (Z12 , +) is a GIS for the function int(x, y) = y − x
mod 12.
(3) The triple (G, G, int) where S is a multiplicative group G is a
GIS for the function int(x, y) = yx−1 .
(4) If S is the diatonic scale isomorphic to Z7 , the triple
(Z7 , (Z, +), int) is a GIS for the function int(x, y) = y −x mod 7.
(5) In an Equal Division of the Octave with N degrees, the triple
(ZN , G, int) where G is the additive group of integers modulo
N (ZN , +) is a GIS for the function int(x, y) = y − x mod N.
(6) If G = D12 is the dihedral group of transpositions and inversions,
the triple (Z12 , D12 , int) is the space of Allen Forte’s set theory.
(7) If G is the permutation group on S, the triple (S, G, int) is a GIS
for the function int(x, y) = σ if and only if σ(x) = y.
(8) In order to modelize rhythms and durations, one can consider
the set of real numbers and its additive group (R, (R, +), int).
The function int(x, y) defines the number of time units between
x and y.
(9) If G = Q is the group of rational numbers and S is a set of
frequencies of the pitches, the triple (S, Q, int) is a GIS for the
function int : S × S → Q int(x, y) = x/y. This is the GIS of Just
intonation.
Definition 2.2. Let (S, G, int) be a GIS. A function f : S → S
preserves intervals if for all x, y ∈ S
int(x, Tn (x)) = n
implies
int(a, x) = int(a, y)
This implies that there exists a point x ∈ S such that Iab (x) = y.
Thus the function Iab is onto, and bijective since Iab is one-to-one and
onto. The previous equality,
shows that Iab (b) = a. In the same way, we show that Iab (a) = b.
Definition 2.5. Let X, Y be two finite sets of the GIS (S, G, int).
The interval function ifunc(X, Y ) of X and Y is defined by, for all
n ∈ G,
ifunc(X, Y )(n) = card{(x, y) ∈ X × Y, int(x, y) = n}
= ifunc(X, Y )(n + i)
(2) In the same way, we have
ifunc(X, Ti (Y ))(n) = 1X (k)1Ti (Y ) (n + k)
k
= 1X (k)1Y (n + k − i)
k
= ifunc(X, Y )(n − i)
Generalized Interval Systems 37
Definition 2.6. In the common set Z12 , the interval vector of a set
of pitch class set X is the vector whose coordinates are the values of
the interval function
for n = 1, 2, . . . , 6.
Example 2.3. The interval vector of set X = {0, 4, 7} is iv(X) =
(3, 0, 0, 1, 1, 1, 0, 1, 1, 1, 0, 0), and its interval-class vector is icv(X) =
(3, 0, 0, 1, 1, 1, 0). The interval vector of set X = {0, 4, 6} is
iv(X) = (3, 0, 1, 0, 1, 0, 2, 0, 1, 0, 1, 0) and its interval-class vector is
icv(X) = (3, 0, 1, 0, 1, 0, 1).
Theorem 2.6. Let X be a subset S = Z12 , and let iv(X) = (n0 ,
n1 , . . . , n11 ) be the interval vector of X and icv(X) = (m0 , m1 , . . . ,
m6 ) the interval-class vector, then
38 A Compendium of Musical Mathematics
(1) nj = mj for j = 0, 1, . . . , 5,
(2) n6 = 2m6
(3) nj = m12−j for j = 7, 8, . . . , 11.
In symbolic form, the interval function is made up of the number
of elements of X, the interval-class vector and its retrogradation.
−→ ←−
ifunc(X, X) = iv(X) = (card(X), icv(X), icv(X))
This result is generalized in the N -EDO.
Theorem 2.7. Let X be a subset S = ZN , and let iv(X) =
(n0 , n1 , . . . , nN −1 ) be the interval vector of X and icv(X) =
(m0 , m1 , . . . , mN/2 ) the interval-class vector, then
(1) nj = mj for j = 0, 1, . . . , N/2 − 1,
(2) nN/2 = 2mN/2
(3) nj = mN −j for j = N/2 + 1, . . . , N − 1.
where x is the ceiling function. In other words, if N = 2k + 1 is
odd
icv(X) = (n0 , n1 , . . . , nk−1 , nk )
and if N = 2k is even
icv(X) = (n0 , n1 , . . . , nk−1 , nk /2)
Proof. The proof follows from the definition of the injection num-
ber and the bijectivity of f ,
inj(X, Y )(f ) = 1(f (x)∈Y ) = card (f (X) ∩ Y )
x∈X
Furthermore,
card f −1 (Y ) ∩ X = card (Y ∩ f (X))
it follows that
= inj(Ac , Ac )(f )
Theorem 2.13. Let (S, G, int) be a GIS, and X a subset of S. For
all bijection f , the injection number satisfies
inj(X, X)(f ) + inj(X, X c )(f ) = card(X)
In particular, in Z12 , if A is a hexachord (|A| = |Ac | = 6), then
ifunc(A, Ac ) = 6 − ifunc(A, A)
Proof. By definition, we get
inj(X, X)(f ) + inj(X, X c )(f ) = 1f (x)∈X + 1f (x)∈X c
x∈X x∈X
= 1 =card(X)
x∈X
More generally, we have the following result:
Theorem 2.14. Let (S, G, int) be a GIS, and f a bijection on S.
We have the following properties for any subset X and Y of S:
(1) inj(X, Y c )(f ) = card(X) − inj(X, Y )(f )
(2) inj(X c , Y )(f ) = card(Y ) − inj(X, Y )(f )
(3) inj(X c , Y c )(f ) = card(Y c ) − card(X) + inj(X, Y )(f )
42 A Compendium of Musical Mathematics
Proof.
(1) For any subset X and Y, we have
inj(X, Y c )(f ) + inj(X, Y )(f ) = 1f (x)∈Y c + 1f (x)∈Y
x∈X x∈X
= 1 =card(X)
x∈X
icv(A) = icv(Ac )
it induces that the interval function are the same for all n
ifunc(A, A) = ifunc(Ac , Ac )
iv(A) = iv(Ac )
icv(A, A) = icv(Ac , Ac )
The following corollary generalizes Babbitt’s hexachord theorem
for a pitch-class set X of any cardinality and X ⊂ ZN .
Corollary 2.2. Let (S, G, int) be a GIS, and X a subset of S = ZN .
Let iv(X) = (n0 , n1 , . . . , nN −1 ) be the interval vector of X a subset
of any cardinality and set
m = card(X c ) − card(X)
thus
iv(X c ) = iv(X) + m
m = card(X c ) − card(X)
and
Definition 2.9. Let (S, G, int) be a GIS. The function sum : S×S →
G is defined for all x, y, z ∈ S, by
sum(y, x) = sum(x, y)
Proof. By definition
and
Definition 2.10. Let (S, G, int) be a GIS, and X, Y two finite sub-
sets of S. The function interval sum sfunc(X, Y ) is defined by, for all
n∈G
= card(In (X) ∩ Y )
Proof.
(1) From the previous result,
sfunc(Ti (X), Y )(n) = 1∗Ti (X) (j)1Y (n + j)
j
= 1Ti (X) (−j)1Y (n + j)
j
= 1X (−j − i)1Y (n + j)
j
= 1X (−k)1Y (n − i + k)
k
= sfunc(X, Y )(n − i)
Theorem 2.20. Let (S, G, int) be a commutative GIS. Set Sab the
sum sum(a, b) and Iab inversion function for a and b in S. For two
finite subsets X, Y ⊂ S, we have
Remark 2.3. In the common EDO ZN , the inversion Iab is the usual
inversion In (x) for n = a + b = sum(a, b). In this particular case, we
find the equality
3
2iπ3kx
− 4 if x ∈ {0, 4, 8}
FX (x) = e 12 =
0 otherwise
k=0
We will not detail here the results on the Discrete Fourier Trans-
form which are too numerous. We refer the reader to the book of
Emmanuel Amiot (2016) and to the special issue of the Journal of
Mathematics and Music (Fiore, 2017).
The elementary properties of the discrete Fourier transform are
as follows:
∀n ∈ ZN , −n
1A (x) = A(ωN )
N −1
A ∗ A∗ (x) = ci xi = iv(A)(i)xi
ci ∈ifunc(A,A) i=0
50 A Compendium of Musical Mathematics
n
= xai −aj = ΔA(x)
i,j=1
ΔA(x) = A(x)A(x−1 )
fact that QA and QB are relatively prime. Since RA (x) and RA (x−1 )
are relatively prime to both RB (x) and RB (x−1 ), RA (x) = ±xα and
RB (x) = ±xβ . Therefore, A(x) = ±xα P (x)V (x) and B(x) = ±
xβ P (x)V (x−1 ). Denote U (x) = ±xα P (x) and set a = β − α, we get
A(x) = U (x)V (x) and B(x) = ±xa U (x)V (x−1 ).
For difference sets (see Chapter 8), we have the following
characterization.
Theorem 2.26. Let A be a k-subset of ZN . Then A is a (N, k, λ)-
difference set if and only if its Patterson function is in Z [ZN ]:
A(x)A(x−1 ) = k + λ(T (x) − 1) mod xN − 1
where T is the polynomial
T (x) = 1 + x + x2 + · · · + xN −1
Another way to approach homometry is to introduce the notion
of spectral unit as it was done by Rosenblatt (1984).
Definition 2.16. A unit S(x) of the algebra group Q[ZN ] is called
a spectral unit if S(x)S(x−1 ) = 1 mod (xN − 1).
Generalized Interval Systems 53
{0, 1, 3, 4}) and three 5-notes homometric pairs in Z10 . In the Z13
universe studied by Erv Wilson and W. Sethares, there is a unique
4-notes pair ({0, 1, 3, 9}, {0, 1, 4, 6}) and only two 5-notes pairs
({0, 1, 2, 3, 6, 10}, {0, 1, 2, 4, 5, 8}) and ({0, 1, 2, 4, 7, 9}, {0, 1, 3,
5, 7, 8}).
2.8 Exercises
Exercise 2.1
Let’s consider the quarter-tone space, which we present in relation
to the usual notes Z24 = Z12 × Z2 . Each note of Z24 is coded (x, n),
where x is a point of Z12 (a note of our classical tempered space) and
n = 0 or 1. If n = 0, the note x is not altered and if n = 1, the note x
is raised by a quarter tone. In this system, (0,0) is C, (0,1) is C+1/4
tone, (1,0) is C , (1,1) is C sharp raised by a quarter tone (C+3/4
tone), (2,0) is D, etc. We also consider the additive group Z24 . Define
an interval function int(x, y) on Z12 ×Z2 such that (Z12 ×Z2 , Z24 , int)
is a GIS.
Exercise 2.2
In the N = 8 universe, we introduce M3 the multiplication by 3
mod 8. Let us consider the injection number inj(X, Y ) of two pcsets
Generalized Interval Systems 55
X and Y
Show that
and
Exercise 2.3
In the N = 12 universe, we introduce M5 the multiplication by
5 mod 12. Let us consider the injection number inj(X, Y ) of two
pcsets X and Y
Show that
and
Exercise 2.4
In the N = 12 universe, we introduce the function ibang defined by
Exercise 2.5
Show that the discrete Fourier transform of Tn (X) where X is a pcset
verifies
Exercise 2.6
Show that X = {0, 1, 2, 5, 7} and Y = {0, 1, 3, 5, 6} are homomet-
ric sets in Z10 by calculating their interval vector. Compute their
Patterson polynomial and show using Theorem 2.26 that they are
(N, k, λ)-difference sets.
Exercise 2.7
Verify that the Z-relation is stable by multiplication. Let X be a set
in ZN , m an integer such that gcd(m, N ) = 1 and m = 1, N − 1.
Denote by Mm the multiplication modulo N , and if X and Y are
homometric sets, then Mm X and Mm Y are also homometric. For
instance, let N = 12, show that the sets X = {0, 1, 2, 3, 5, 6} and
Y = {0, 1, 2, 3, 4, 7} are homometric and compute their interval con-
tent. Verify that the sets M5 X and M5 Y are also homometric and
compute their interval content.
References
57
58 A Compendium of Musical Mathematics
where x is the floor function (the greatest integer less than or equal
to x) and Douthett’s J-set is
m
JN,k = {JN,k
m
(x), x = 0, 1, . . . , k − 1}
Definition 3.9. A scale S has the Myhill property if each scale inter-
val comes in two chromatic sizes. In other terms, S has exactly two
c-distances for every d-distance.
Example 3.6. The major scale {0, 2, 4, 5, 7, 9, 11} satisfies the
Myhill property.
Definition 3.10 (Clough et al., 1993). A P-cycle is a cycle of
three or more pairwise distinct sets (scales or chords) from the same
set class (sets equivalent under transposition and inversion) such that
there is a map between each pair of adjacent sets that leaves all but
one of its notes fixed. If the sets in the cycle exhaust the set class,
then the cycle is said to be unidirectional, otherwise it is a toggling
P-cycle.
All P-cycles are unidirectional for our purposes.
Example 3.7. The cycle of fifths is unidirectional
7T 7 T 7 T 7 T 7 T
{0, 7} −→ {7, 2} −→ {2, 9} −→ · · · −→ {5, 0} −→ {0, 7}
For each choice of index m, there exists only one class of set
of progressive transposition. More precisely, we have the following
result.
X = {x, Tm (x), . . . , Tm
k−1
(x)}
u = max(−r, r) (3.1)
Proof. The three consecutive notes are (m, 3m, 5m) with m =
(N + 1)/2. Since 3(N + 1)/2 = (N + 1)/2 + N + 1 = (N + 1)/2 + 1
mod N , 5(N + 1)/2 = (N + 1)/2 + N + 2 = (N + 1)/2 + 2 mod N .
A = {3j, j = 0, 1, 2, . . . , (N + 1)/4}
and
B = {3j + m, j = 0, 1, 2, . . . , (N − 3)/4}
X = {x, Tm (x), . . . , Tm
k−1
(x)}
68 A Compendium of Musical Mathematics
u≤k≤v (3.3)
Dia11
13 = {0, 1, 3, 5, 7, 9, 11}
From these scales, we can choose a unique diatonic class scale (up
to inversion and transposition) by choosing the smallest m.
which verifies the properties of the diatonic scale. But only Dia11
20
verifies the condition on m (the smallest integer such that 2m >
N + 1).
Proposition 3.5. The complement of X does not contain two con-
secutive notes if and only if X c ⊂ T1 (X).
Proof. If (x, x+1) ∈ X c , then T−1 (x, x+1) = (x−1, x) ∈ / X since x
is in X. Thus, (x, x+1) ∈ / T1 (X). Conversely, since the transpositions
are bijections, if (x, x+1) ∈ X, then T1 (x, x+1) = (x+1, x+2) ∈ / X c.
Ld = {dj, j = 0, 1, . . . , k}
Ld = {dj, j = 0, 1, . . . , q − 1}
Generalized Diatonic Scales 71
3 0 11 8 5 2 13 10 7
has nine elements. For k = 8, the first eight multiples define the
major scale with a unique tritone
• Case N even
If N is even, we have a graphical representation of these scales.
The set {N − m, (N/2 − 1)m} is a generalized tritone. The major
scale with tonic 0 in the N universe is
In this case, the diatonic scale is equal to the major scale. Its graphic
representation is in the frame (N − m, N/2) which has the multiples
of N − m along the x-axis and the multiples of N/2 along the y-axis.
A linear representation can be obtain if the two first notes {0, N −m}
are added after m on the right of the graph. The two ends of the
graph are then separated by the unique generalized tritone of the
major scale:
0 — N −m
|
(N/2 − 1)m — · · · — 4m — 3m — 2m — m
The relative minor scale has one more generalized tritone {2m, 2m +
N/2}
The note m is replaced by 2m + N/2. The two scales are in the same
neighborhood (differ by only one note):
0 — N −m 2m + N/2
| |
(N/2 − 1)m — · · · — 4m — 3m — 2m
Wk = {0, 1, 3, 5, . . . , 2k + 1, 2k + 2, 2k + 4, . . . , 4k + 2}
N = m + |Wk | = 4k + 4
3.5 Exercises
Exercise 3.1
Compute the hierarchy of the well-formed scales for N = 12, m = 7
and k from 2 to 12 starting at x = 5. Compute Douthett’s J-sets
7
J12,k for k = 2, . . . , 7. Compare with the well-formed scales hierarchy.
7
Show that J12,7 is the major scale.
Exercise 3.2
Consider the C harmonic minor scale X = {0, 2, 3, 5, 7, 8, 11}, the
ascending C melodic minor scale Y = {0, 2, 3, 5, 7, 9, 11} and the
descending melodic minor scale Z = {0, 2, 3, 5, 7, 8, 10}. Are they
maximally even? Which scales have Myhill’s property?
Exercise 3.3
Consider the scales X = {0, 2, 3, 5, 6, 8, 9, 11} and Y = {0, 3, 6, 10}.
Are they maximally even?
Exercise 3.4
Consider the subset X = {0, 2, 4} of the whole tone scale
S = {0, 2, 4, 6, 8, 10}. Prove that cardinality equals variety and struc-
ture implies multiplicity are false.
Exercise 3.5
Consider the X = {0, 2, 4} of the pentatonic scale S = {0, 2, 4, 7, 9}.
Prove that cardinality equals variety and structure implies multiplic-
ity are true.
Generalized Diatonic Scales 75
Exercise 3.6
For N = 13 and for the generator m = 8, compute the generated scale
X of k = 8 notes. Show that X is a set of progressive transposition.
8
Compute the Douthett’s scale J13,8 and connect it to X. Deduce that
X is maximally even. Compute the value of u and v using 3.1 and
3.2. Deduce that X is a diatonic scale of the 13-EDO.
Exercise 3.7
For N = 16 and for generator m = 9, compute the generated scale
X of 9 notes. Show that X is a PT set and a diatonic scale. Deduce
that X is also the generalized major scale.
References
77
78 A Compendium of Musical Mathematics
restricted and does not correspond to the set of all the actions of G
on V , for musical purposes.
(1) they must have the same configuration of vertices and edges,
(2) there must be some isomorphism F that maps the
transformation-system used to label the edges of one network
into the transformation-system used to label the edges of the
other,
(3) the transformation X labels an edge of the one network, then the
transformation F (X) labels the corresponding edge of the other.
We will see that for the T /I group there are four types of automor-
phisms which define four types of isographies. Two types of isogra-
phies are said to be positive and negative. In this section and in
the next one, we follow David Lewin (1990) where one can find the
demonstrations.
Voice Leading and Neo-Riemannian Transformations 79
Lemma 4.2.
(1) Tn M5 has order 2 if n is even and order 4 if n is odd.
(2) Tn M7 has order 2 if n is a multiple of 3 and of order 6 otherwise.
Proof. From (Tn M5 )(Tn M5 ) = Tn T5n = T6n , we see that if n is
even, T6n = T0 and so Tn M5 has order 2. If n is odd, T6n = T6 and
so Tn M5 has order 4. From (Tn M7 )(Tn M7 ) = Tn+7n M1 = T8n , we
see that if n is a multiple of 3, T8n = T0 and so Tn M7 has order 2.
Otherwise, T8n = T8 if n = 1 mod 3, and T8n = T4 if n = 2 mod 3.
So, Tn M7 is of order 7.
Corollary 4.1. The PLR group acts simply transitively on the set
of consonant triads.
In this section, we use the term JQZ group synonymously to the term
T /I group. This is in analogy to the usage of the term PLR group
synonymously to the term S/W group (Schritt/Wechsel group). This
particular system of generators J, Q and Z has very similar properties
to the generators P , L and R of the Schritt–Wechsel group. They
are not contextually defined and can be applied pointwise. Their
definition is unique up to conjugation. Our concrete choice depends
upon the choice of the C-major triad {0, 4, 7} as a distinctive chord of
reference. The choice of any other consonant triad f ({0, 4, 7}) in this
Voice Leading and Neo-Riemannian Transformations 87
Definition 4.11. The JQZ group is the group of order 24, generated
by the JQZ transformations, with presentation:
J, Q, Z = J, Q, Z | J 2 = Q2 = Z 2 = (JQ)3
= (JZQ)2 = (JZ)4 = 1
Bm D Bm Db
R QZQ
ZJZQ
E / C ZJQJ / D ZJQJ/ E ZJQJ / E
O RP RL P RP L P RP L P RP L
RL QZ
JQJZ QJ QZ
E o Eo Co F o JZJZ B
P LP R PL RL P RP R
ZQ
B
J
Dm / Gm JQZ C ZQZ Dm J
B
ZQZQ
/ C ZQ / F
L RL P RL RLR L RLRL LR
As in the case of the P LR group, two subgroups of the JQZ group are
important: the hexatonic group generated by transformations J and
Q and the octatonic group generated by transformations J and Z.
The group acts on the major and minor chords, but likewise it acts
on any other set class. The 3-chord {0, 6, 11} has the 6-element orbit
{{0, 6, 11}, {1, 7, 8}, {3, 4, 10}, {3, 4, 9}, {2, 7, 8}, {0, 5, 11}}
Since the groups T/I, PLR and JQZ are isomorphic, they have
the same center. This centre is the group of order 2 consisting of
{1, T6 }.
Introduced by Hook (2002), the group of uniform triadic trans-
formations (U T T ) acts on the set of 24 major and minor triads Σ.
Its elements are denoted by (p, q, σ) with p, q ∈ {0, 1, . . . , 11} and
σ = ±1. Let (x)M be a major triad rooted on x and (x)m a minor
triad rooted on x, the (p, q, +1) U T T elements send the major triad
(x)M to (x+p)M and the minor triad (x)m to (x+q)m . The (p, q, −1)
U T T elements send the major triad (x)M to the minor triad (x+p)m
and the minor triad (x)m to the major triad (x + q)M .
Proposition 4.4. The group U of uniform triadic transformation
(UTT) of order 288 is isomorphic to the wreath product of the cyclic
group C12 of order 12 by the cyclic group C2 of order 2:
U = C12 C2
The orbits in the first pair are triads in proper hexatonic collections,
i.e. they are also orbits of the hexatonic PL-group P, L. The orbits
in the second pair form orbits under a conjugate of P, L, namely
R P, L R. This leads to the following result.
Theorem 4.12. With respect to the first two orbits, the group J, Q
is subdual to P, L. With respect to the second two orbits, J, Q is
subdual to R P, L R.
The octatonic JZ group J, Z has three orbits of eight triads
each. The first one {F m, F , E m, A, Cm, C, Am, E } is also the
orbit of the octatonic P R group. The two others {Bm, C , Gm, E,
Fm, G, Dm, B} and {Bm, D, Gm, F, Em, G, C m, B } are orbits
under conjugation (RL) P, R (RL)−1 . This lead to the following
result.
Theorem 4.13. With respect to the first orbit, the group J, Z
is subdual to P, R. With respect to the second two orbits, J, Z is
subdual to (RL) P, R (RL)−1 .
Some other groups have been introduced to generalized neo-
Riemannian transformations. Alexander Popoff (2013) defined the
Z3 Q8 groups as follows, where Q8 is the quaternion group of
presentation:
Q8 = i, j, k | i2 = j 2 = k 2 = ijk
It has also a presentation with two generators:
Q8 = u, v | u4 = 1, u2 = v 2 , uvu = v
This group acts simply transitively on the set of major and minor
triads. It is generated by the transformations u and v such that
u(xM ) = (x + 1)M u(xm ) = (x − 1)m
v(xM ) = xm v(xm ) = (x + 6)M
Verification of relationships is obvious. For instance, for the last
one, we have
uvu(xM ) = uv(x + 1)M = u(x + 1)m = xm = v(xM )
uvu(xm ) = uv(x − 1)m = u(x + 5)M = (x + 6)M = v(xm )
96 A Compendium of Musical Mathematics
The study of the interplay of the PLR group and the JQZ group is
helpful in atonal analysis. We will illustrate this for the the atonal
triad 0, 1, 6.
The main advantage of JQZ transformations is to be able to con-
sider all types of chords and not only consonant chords. For instance,
if the JQZ group acts on the atonal triad 016 {c, c, f }, the action
leads to a new lattice. In the other hand, if we replace the set Σ of
consonant triads by the set Γ of atonal triads rooted on x of the form
x, x + 1, x + 6 and x, x − 1, x − 6, we get new relations
P x, x + 1, x + 6 = x + 6, x + 5, x ,
P x, x + 11, x + 6 = x + 6, x + 7, x
R x, x + 1, x + 6 = x + 1, x, x + 7 ,
R x, x + 1, x + 6 = x + 11, x, x + 5
L x, x + 1, x + 6 = x + 7, x + 6, x + 1 ,
L x, x + 11, x + 6 = x + 5, x + 6, x + 11 .
2
107
JZ / 43A QZ / BA5 JZ / 218 QJ / 650 (ZJ) / 0B6
(LP ) 3 RP 3(LP ) 4 (LP ) LR
The right hand evolves in the same way as the left hand: the
atonal triads are linked by the transformations (JZ, QJ). The same
transformations are used in both hands, except in the last triads.
But if we consider a crossing of the hands, 781 goes to 0B6 by the
transformation J, just as 650 goes to 127 by the same transformation.
Another interpretation is possible. As we saw for the PLR group,
we can consider two kinds of atonal triads: one of structure (1, 5) and
the other of structure (5, 1). From the relations of P LR seen above,
we can introduce the transformations for atonal triads {x, x+1, x+6},
JQ JQ JQ JQ JQ
650 / 54B / 43A / 329 / 218 / 107 / 107
PL PL PL PL PL
4
RP L ZJ Q JZ RP ZQ RL ZJ RP L (LP ) ZQJ (LP )3 ZJ J Q J Q (RP )2
JQ
A94 / A94 / A94 / A94 / A94 / A94 / 0B6
LP
T2 R123 T2 T1 R123
EM 7 / F m7 / Gm7 / AM 7 T11 R123/ Gm7 T10
/ F m7 T10 R123/ EM 7
RC3 = E8 × C12
R1 , R2 , R3 , T | R12 = R22 = R32 = T 12 = (R1 R2 )2
= = (R2 R3 )2 = (R1 R3 )2 = T −1 R1 T R1
= T −1 R2 T R2 = T −1 R3 T R3
T5 T5 T5 T9 R1 T5 R13
D7 /F / C7/E / F 7/E / BΔ /D / Gm7 5 / C7
T5 R3 T3 R3 T5 R3 T3 R3 T5 R3 T6 R123 T5 R12
D7 / GΔ / B7 / EΔ / F 7 / BΔ / F m7 / B7
T5 R3
B7 / EΔ T6 R123/ Am7 T5 R12 / D7 T5 R3
/ GΔ T6 R123
/ Cm7 T5 R12 / F 7
T5 R3
F 7 / BΔ T6 R3 123/ F m T5 R12
/ B7 T5 R3
/ EΔ T10 R123/ C7 T5 R12
/ F 7
R1 , R2 , T, I | R12 = R22 = I 2 = T 12 = (T I)2
RD3 = E4 × D12 = = (R1 R2 )2 = R1 T R1 T −1 = R2 T R2 T −1
= (R1 I)2 = (R2 I)2 = 1
Voice Leading and Neo-Riemannian Transformations 103
T9 R12 I1
Gm7 / E7 / Bm75 T5 R2 / EmM 7 I0
/ AΔ T3 R1 / Cm7
Fig. 4.5.
Action of the group RD3 on the set U.
104 A Compendium of Musical Mathematics
R1 , . . . , Rn−1 , T, I | (Ri Rj )2 = I 2 = T 12
RDn = E2n−1 × D12 = = Ri T Ri T −1 = (Ri I)2
= (T I)2 , i, j = 1, . . . , n − 1
The group has order |RDn | = 2n−1 × 24 = 2n × 12 = |RCn |. The
groups are included in each other:
RD1 ⊂ RD2 ⊂ RD3 ⊂ RD4 ⊂ · · ·
Voice Leading and Neo-Riemannian Transformations 105
4.10 Exercises
Exercise 4.1
In 1880, Hugo Riemann presented two kinds of operations in the
tonnetz, Schritt and Wechsel. The Schritt Sn transposes a major
triad up by n semitones or a minor triad down by the same amount.
S0 is the identity. By indexing the major triad by capital M and
minor triad by m, Sn is defined by its values on the triads: S1 (CM ) =
CM , S1 (Cm ) = Bm , S2 (CM ) = DM , . . . , S11 (CM ) = BM , S11 (Cm ) =
Cm . The Wechsel Wn first transposes a major triad up and a minor
triad down by n semitones and then converts it to a triad of opposite
mode on the same root. Thus, Wn = P Sn , with P the usual neo-
Riemannian transformation: W0 (CM ) = Cm , W1 (CM ) = Cm ,. . . ,
W11 (CM ) = Bm . (1) Verify that W0 = P , W4 = L and W9 = R. (2)
For all triad X, show that (indices are considered mod 12)
Si Sj (X) = Si+j (X)
Si Wj (X) = Wj−i (X)
Wi Sj (X) = Wj+i (X)
Wi Wj (X) = Sj−i (X)
Exercise 4.2
Compare the action of the T /I subgroup {T0 , T4 , T8 , I3 , I7 , I11 } on
the set {CM , EM , AbM , Cm , Em , Abm } with the action of the S/W
subgroup {S0 , S4 , S8 , W0 , W4 , W8 } on the same set.
Exercise 4.3
Assign order numbers to the 24 triads such that C major is 1, C
major is 2, D major is 3, . . . , B major is 12, and C minor is 13,
C minor is 14, . . . , B minor is 24. Consider the permutation of the
Schritt
S = (1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12)
(24, 23, 22, 21, 20, 19, 18, 17, 16, 15, 14, 13)
and the Wechsel
W = (1, 13)(2, 14)(3, 15)(4, 16)(5, 17)(6, 18)(7, 19)
(8, 20)(9, 21)(10, 22)(11, 23)(12, 24)
106 A Compendium of Musical Mathematics
F m / Dm
References
Rec(M ) = Rat(M )
which is a rotation of ω.
Example 5.12. The converse is false. For instance, the bossa nova
rhythm 33433 which is written on the circle {0, 3, 6, 10, 13}16 has the
following distances
0 3 6 10 13
0 0 3 6 6 3
3 3 0 3 7 6
6 6 3 0 4 7
10 6 7 4 0 3
13 3 6 7 3 0
and occurrences
d 1 2 3 4 5 6 7
i 0 0 4 1 0 3 2
Example 5.13. The bossa nova rhythm {0, 3, 6, 10, 13}16 is E-deep;
it is the set of multiples for m = 13, starting at 6, {6, 3, 0, 13, 10}.
The rhythm {0, 4, 5, 9, 10, 14, 15}16 is W-deep and thus E-deep. It is
the set of multiples of m = 5, {0, 5, 10, 15, 4, 9, 14}.
or
a if (n + 1)α + ρ − nα + ρ = 0
ωn =
b otherwise
C(8/5) = aabaababaabab
The upper Christoffel path is the mirror image of the lower Christoffel
path. Looking at squares crossed by the line D, the discrete path
joining the center of each square encodes the billiard word:
B(8/5) = xyxxyxyxxyxy
m Nx + m
JN,k (x) = mod N
k
David Clampitt and Thomas Noll (see Clampitt and Noll, 2011;
Noll, 2009) have developed a theory of diatonic modes based on the
Sturmian monoid that we now define.
Definition 5.40. A morphism f is Sturmian if for every Sturmian
word ω, f (ω) is a Sturmian word.
Example 5.25. The morphisms identity 1 : a → a, b → b and swap
S : a → b, b → a are Sturmian. The morphisms
La : a → a, b → ab, Ra : a → a, b → ba
and
Lb : a → ba, b → b, Rb : a → ab, b → b
are Sturmian.
La L−1 −1 −1
b La = Lb La Lb , La Rb−1 La = Rb−1 La Rb−1
Ra L−1 −1 −1
b Ra = Lb Ra Lb , Ra Rb−1 Ra = Rb−1 Ra Rb−1
(La L−1 4
b Ra ) = 1, (L−1 −1 4
b Ra Rb ) = 1
σi σj = σj σi if |j − i| ≥ 2
and
π π =
1 / 2Z(B3 ) / B3 / GL2 (Z) det / Z2 /1
|u|a |v|a
∈ SL2 (Z)
|u|b |v|b
ω = C(5/8) = aabaababaabab
= (u, v) = (aabaabab, aabab)
and
|u|a |v|a 53
det = det = +1
|u|b |v|b 32
and
ω = C(4/7) = xxyxxyxxyxy
ω ∗ = C(3/8) = xxxyxxxyxxy
In 1961, Schensted stated that the length of the first row of P (ω)
is equal to the maximal length of nondecreasing subwords (in our
examples, the three subwords 145, 125, 123 have length 3) and the
length of the first column of P (ω) is equal to the maximal length
of the decreasing subwords. Some years later, Greene introduced the
equivalence relations on two words of A∗ if they have the same Schen-
sted tableau:
ω ∼ ω ⇐⇒ P (ω) = P (ω )
Moreover, in 1970, Knuth (1970) proved that the equivalence ∼
between two words coincides with the plactic congruence.
Another tool for studying the plactic monoid is the RSK
(Robinson–Schensted–Knuth) correspondence.
Definition 5.47. The RSK map is associate with each word ω ∈
A∗ over the linearly ordered alphabet A = {1, 2, . . .}, the pair
(P (ω), Q(ω)), where P (ω) is the Schensted semistandard Young
tableau and Q(ω) is a standard Young tableau.
Proposition 5.17. The RSK map ω → (P (ω), Q(ω)) is one-to-one.
A standard tableau is a semistandard tableau where each k ∈ A
appears exactly once. By definition, Q(ω) is the standard tableau
encoding the chain of shapes in the construction of P (ω). In the
previous example, Q(142531) is the tableau:
1 2 4
3 5
6
Combinatorics on Musical Words 135
std(142531) = 153642
6 1 3 bbaaab 9
6 1 4 baaaba 9
7 1 2 bababbb 14
8 1 2 babababa 14
8 1 2 aabababb 20
8 1 2 abababab 28
8 1 3 aaaababa 20
9 1 2 aaaababab 27
9 1 2 aaabababa 48
10 1 2 aaaaaababa 35
Combinatorics on Musical Words 137
(2) If ev(ω) = (3, 2). Each word ω has three letters a and two
letters b. In the 12-tone equal temperament, the modes correspond
to a = 2, b = 3 semitones. Since there are three possible shapes,
there are three plactic classes. The first class has only one element
a a a b b
The second class has four elements given as follows. In this class, the
pentatonic minor mode (22323) is in relation with the Blues major
mode (23223):
aabba ∼ aabab ∼ abaab ∼ baaab
22332 22323 23223 32223
All elements of this class have the same Schensted tableau P (ω):
a a a b
b
Finally, the third class has five elements. In this class, the Egyptian
mode (23232) is linked to the pentatonic major mode (32232) which
is related by plactic congruence to the Blues minor mode (32322):
abbaa ∼ ababa ∼ baaba ∼ babaa ∼ bbaaa
23322 23232 32232 32322 33222
The Schensted tableau has the form
a a a
b b
baabbbb
|
ababbbb
|
abbabbb − abbbabb − abbbbab − abbbbba
locrian phrygian natakapriya kikilapriya
140 A Compendium of Musical Mathematics
The second class is the plactic class of the major mode (bbabbba
with a = 1, b = 2). These modes appear in John Foulds’s classifica-
tion under the Vth class, and three of them are used in the Essays
in the modes for piano (1920–1927). The plactic heptatonic class has
14 elements and its graph is as follows (see Fig. 5.3).
1 r s−1
hi−2 = Nar(r, s, i) =
r i i−1
144 A Compendium of Musical Mathematics
For the noncrossing partitions point of view, the number of (r, s)-
Dyck paths with αj vertices runs of length j is the Kreweras number:
(s − 1)!
Krew(r, s, α) =
α0 !α1 ! · · · αr !
The f -vector is the integer sequence (f−1 , f0 , . . . , fr−2 ) of Ass(r, s)
with f−1 = 1, and fi is the number of i-dimensional faces 0 ≤ i ≤
r − 2, also given by the Kirkman numbers:
1 r s+i−1
fi−2 = Kir(r, s, i) =
r i i−1
These numbers verify the relation
r−2
r−2
fi (t − 1)r−2−i = hi tr−2−i
i=−1 i=−1
Example 5.34. For instance, for the Ass(3, 5) of Fig. 5.5, h-vector
= (1, 4, 2), f -vector = (1, 6, 7) and the reduced Euler characteristic
χ = −2 which is equal to (−1)3 Cat (3/2) = −Cat(2) = −2.
Dyck paths are related by plactic relations. Dyck path of Ass(3, 5)
is a class of this monoid with three letters x and five letters y, namely
000011122236
134534534547
268787676858
At each vertex, there are three blocks which form a parallel class
(a partition of Z9 ), except when the ribbon wraps.
5.8 Exercises
Exercise 5.1
Write down the Euclidean string E(4, 7). (It is the Bulgarian rhythm
called Ruchentiza).
Exercise 5.2
Is the Arab rhythm Agsag-Samai 12222 a Euclidean string?
Exercise 5.3
For integer n ≥ 2, write down the Euclidean string E(n, 2n + 1),
which is the hierarchy of some Aksak rhythms.
Exercise 5.4
Give two different proofs that Euclidean rhythm E(4, 9) is E-deep.
Proof that the Aksak rhythms hierarchy is E-deep.
Exercise 5.5
Fume-fume is an African bell pattern of 12 pulses {0, 2, 4, 7, 9}12 .
Show that fume-fume is well-formed, maximally even and deep.
Exercise 5.6
Show that the Bulgarian rhythm E(7, 18) = 3232323 has the rhyth-
mic oddity property (1) by using the definition (2) by using Theorem
5.13 and (3) by using Theorem 5.16.
Combinatorics on Musical Words 147
Exercise 5.7
Using Theorem 5.16, show that the word ω = 332232322 has the
rhythmic oddity property.
Exercise 5.8
Let us consider the Christoffel word ω = C(8/3) on the alphabet
{a, b}. By drawing the Dick path, write down the standard factor-
ization of ω = (ω1 , ω2 ). Verify that
|ω1 |a |ω2 |a
∈ SL2 (Z)
|ω1 |b |ω2 |b
References
Rhythmic Canons
The first canons appeared in the 12th century. The french chace
and the italian caccia are primitive forms of canons. Despite the
fact that canon is an old procedure, it has been used at any time,
from Renaissance composers (G. Ockeghem, M. Forestier, etc.) to
21st century composers. Famous canons have been composed by
Jean-Sébastien Bach (Musical Offering, Goldberg Variations, Canon-
ical Variations Vom Himmel hoch BWV 769 ), Gabriel Fauré
(Romances sans paroles, Thèmes et variations, Prélude 6 ), Arnold
Schoenberg (Opus 25, 28, 40, 41), Alban Berg (Kammerkonzert),
Anton Webern (Opus 2, 16, 28, 30, 31), Igor Stravinsky (Octuor,
Sextuor, In memoriam Dylan Thomas), Paul Hindemith (Fuga
undecima in B in Ludus Tonalis) and many others. For example,
Olivier Messiaen uses a rhythmic canon in the seventh song of Harawi
called Adieu. The rhythms are counted relatively to the double qua-
ver (e.g. 3 means three double quavers). The upper voice juxtaposes
three nonretrogradable rhythms: 35853, 43534 and 2235322. The
same rhythms appear in the two other voices shifted in time. From
a mathematical point of view, canons are linked to the decomposi-
tion of finite Abelian groups. Some results have been set by Amiot
(2005, 2011). Enumeration of rhythmic canons have been studied by
Fripertinger (2001).
151
152 A Compendium of Musical Mathematics
6.1 Tilings
Let’s start with the most general situation. Let G be a locally com-
pact Abelian group and G its dual group. Let Ω be a subset of G of
finite nonzero Haar measure.
Definition 6.1. The set A is a tile if there exists a set B in G, called
the tiling set, such that
G=A⊕B
In other words, if every element g of G has unique decomposition,
g = a + b with a ∈ A and b ∈ B.
The pair (A, B) is called a tiling pair, or a factorization of the
group G. B is sometimes called a complementary tile to A. Usual
cases are G = Rd , R2 , R, Z and ZN = Z/N Z. In the first approach,
only the cyclical group ZN interests the musicians.
Example 6.1. The tile A = {0, 2} tiles the set Z4 = {0, 1, 2, 3} with
translations B = {0, 1}, i.e. with the sets A and A + 1 = {1, 3}.
This tiling corresponds to the decomposition of Z4 in the direct sum
Z4 = A ⊕ B = {0, 2} ⊕ {0, 1} = 2 · I2 ⊕ I2
where IN is the least residue system modulo N:
IN = {0, 1, . . . , N − 1}
In the musical domain, a factorization (A, B) is called a canon or
a tiling canon. The set A represents the set of events of the ground
voice (dux ) and the set B represents the time attacks: A + b with
b ∈ B are the different voices (comes) of the canon. The exchange
between the inner and the outer voices leads to the same results:
if (A, B) is a canon, then (B, A) is a canon and conversely.
Remark 6.1. If (A, B) is a canon of ZN , then |A| × |B| = N . In par-
ticular, |A| and |B| are divisors of N and thus the product of disjoint
sets of prime factors of N .
Rhythmic Canons 153
d|n Φd (x) = x − 1 =
Proof. The proposition follows from n
Φp (x) = 1 + x + x2 + · · · + xp−1
and
α
Φpα+1 (x) = Φp (xp )
Proof. As d|p Φd (x) = xp − 1 = Φ1 (x)Φp (x), so
xp − 1 xp − 1
Φp (x) = = = 1 + x + x2 + · · · + xp−1
Φ1 (x) x−1
α+1 α+1 α+1
Φ1 (x) Φpk (x) = Φs (x) = xp − 1 = (x − 1)(xp + · · · + 1)
k+1 s|pα+1
It follows
α+1
1 + x + x2 + · · · + xp −1
Φpα+1 (x) =
Φp (x)Φp2 (x) · · · Φpα (x)
α α α α
1 + xp + x2p + · · · + x(p−1)p 1 + x + x2 + · · · + xp −1
=
1 + x + x2 + · · · + xpα −1
α α α α
= 1 + xp + x2p + · · · + x(p−1)p = Φp (xp )
156 A Compendium of Musical Mathematics
Φs (xt ) = Φrs (x)
r|t
Proof.
n
rj
n
rj
= Φs (1) Φpl (1) = Φs (1) plj
j
j=1 l=1 j=1 l=1
n
r
= Φs (1) pj j = Φs (1)s
j=1
thus s = 1.
Rhythmic Canons 157
implies that A ⊕ B = ZN .
158 A Compendium of Musical Mathematics
(3) ⇒ (4). If (3), then there exists a polynomial Q(x) such that
for some integer polynomials Q(x), R(x) with deg(R) < N . Since
the product A(x)B(x) is by hypothesis divisible by the product of
all Φs (x), it is divisible by 1 + x + · · · + xN −1 . It follows that R(x) =
C(1 + x + · · · + xN −1 ), where C is a constant. Taking x = 1 leads to
A(1)B(1) = R(1) = N = CN . Therefore, C = 1, and
A(x)B(x) = Q(x)(xN − 1) + (1 + x + · · · + xN −1 )
≡ (1 + x + · · · + xN −1 ) mod (xN − 1)
A − A = {a − a ; a, a ∈ A}
and let Z(A) the set of the zeroes of A
Z(A) = {k ∈ ZN , A(k) = 0}
given by the following table. The set Z15 is tiled by five augmentations
of A = {0, 1, 2}:
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14
A+2 0 0 1 1 1 0 0 0 0 0 0 0 0 0 0
2A + 4 0 0 0 0 0 0 0 0 1 0 1 0 1 0 0
4A + 5 0 0 0 0 0 1 0 0 0 1 0 0 0 1 0
5A + 1 0 1 0 0 0 0 1 0 0 0 0 1 0 0 0
7A 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1
Rhythmic Canons 161
751114572425247
Therefore,
G=A⊕B
is a Vuza canon of ZN .
Proof. Since kIN can be seen as a subgroup of ZkN , we can apply
the previous theorem. Thus, (A, B ∩ kIN ) is a Vuza canon of kZN .
Taking the quotient by k, since A , B are in ZN , (A , B ) is a Vuza
canon of ZN .
Example 6.7. For instance, consider the Vuza canon A = {0, 16,
32, 36, 52, 68} and
B = {0, 7, 12, 15, 24, 33, 34, 45, 46, 55, 57, 58, 63, 72, 84, 96, 103,
105, 106, 111, 117, 118, 129, 130}
K1 = n2 n3 (Ip2 ⊕ p2 n1 Ip1 )
K2 = n1 n3 (Ip1 ⊕ p1 n2 Ip2 )
is a Vuza canon.
170 A Compendium of Musical Mathematics
Ia ⊕ aIb = Iab
cIa ⊕ aIbc = Iabc mod abc gcd(a, c) = 1
B1 = K1 ∪ T1 (K2 ) = {0, 1, 5, 6, 12, 25, 29, 36, 42, 48, 49, 53}
and
gives the Vuza canon (A1 , B1 ) for the non-Hajós group Z72 . There
are three nonisomorphic solutions for the inner rhythm, namely A1 =
{0, 18}⊕{0, 8, 16}, A2 = {0, 18}⊕{0, 8, 40}, A3 = {0, 18}⊕{0, 8, 64},
but there is only one nonisomorphic solution since A2 = 5A1 and
A3 = 7A1 .
xN − 1
| A(x)
xk − 1
Proof. We have
xN − 1
A(x)
xk − 1
in other words iff for all divisors of N , different of N , there exist a
divisor d of N which is not a divisor of k such that the cyclotomic
polynomial Φd (x) A(x).
Theorem 6.14. Let N = n1 n2 n3 p1 p2 ∈ N where n1 , n2 , n3 , p1 ,
p2 > 1, and gcd(n1 p1 , n2 p2 ) = 1. Let us denote
A1 = n1 p1 n3 In2 A2 = n2 p2 n3 In1
U1 = n1 n2 n3 p1 Ip2 U2 = n1 n2 n3 p2 Ip1
V1 = n2 n3 Ip2 V2 = n1 n3 Ip1
K1 = {0} K2 = {1, 2, . . . , n3 − 1}
A = A1 ⊕ A2
B = (U1 ⊕ V2 ⊕ K1 ) (U2 ⊕ V1 ⊕ K2 )
A1 ⊕ U1 ⊕ V2 = A1 ⊕ U1 ⊕ V1
We have
A1 ⊕ U1 ⊕ V2 = n1 p1 n3 In2 p2 ⊕ n1 n3 Ip1
= n1 n3 (Ip1 ⊕ p1 In2 p2 )
= n1 n3 In2 p1 p2
A1 ⊕ U1 ⊕ V1 = n1 p1 n3 In2 p2 ⊕ n1 n2 n3 p2 Ip1
= n1 n3 (p1 In2 p2 ⊕ n2 p2 Ip1 )
= n1 n3 In2 p1 p2
Rhythmic Canons 173
6.7 Exercises
Exercise 6.1
Let N = 16 and A = {0, 5, 8, 13}. Write down the mask polynomial of
A(x) and its decomposition into cyclotomic polynomials. Find a set
B such that A⊕B = ZN . Using Coven–Meyerowitz (CM) conditions,
prove that A tiles Z16 .
Exercise 6.2
Let us consider the rhythmic canon A ⊕ B = ZN , k ∈ N× and the
set
(1) Prove that C ⊕ B = ZkN or, in other words, that the concatena-
tion k times of A tiles ZkN with the same tile B. (2) Prove that the
CM condition (T2 ) is stable by concatenation, that is, if A tiles ZN
and satisfies (T2 ), then the set C verifies (T2 ).
Exercise 6.3
Prove that the CM condition (T1 ) is not sufficient for A to tiles. Hint:
Consider the set A = {0, 1, 3, 4, 6, 7}.
Rhythmic Canons 175
Exercise 6.4
Let A ⊂ N be finite and n ∈ N. Let us consider B = A + n. (1) Prove
that A satisfies (T1 ) iff B satisfies (T1 ). (2) Prove that A satisfies
(T2 ) iff B satisfies (T2 ).
Exercise 6.5
Let A ⊂ N be finite and n ∈ N. Let us consider B = nA. (1) Prove
that A satisfies (T1 ) iff B satisfies (T1 ). (2) Prove that if n = p is
prime, A satisfies (T2 ) iff B satisfies (T2 ).
References
Redéi, L. (1950). Ein Beitrag zum Problem der Faktorisation von endlichen
Abselschen Gruppen, Acta Mathematica Academiae Scientiarum
Hungarica 1, 197–207.
Sands, A.D. (1957). On the factorisation of finite Abelian groups. I, Acta
Mathematica Hungarica 8, 65–86.
Sands, A.D. (1962). On the factorisation of finite Abelian groups. II, Acta
Mathematica Hungarica 13, 153–169.
Sands, A.D. (2004). Factoring finite Abelian groups, Journal of Algebra 275,
540–549.
Szabo, S. and Sands, A.D. (2009). Factoring Groups into Subsets, CRC Press,
Boca Raton.
Vuza, D.T. (1991). Supplementary sets and regular complementary unending
canons, Perspectives of New Music 29(2), 22–49; 30(1), 184–207; 30(2),
102–125; 31(1), 270–305.
Chapter 7
Serial Knots
177
178 A Compendium of Musical Mathematics
forms use the same chord diagram, viewed in the space, possibly
mirroring or rotating. That is why chord diagram can be used to
classify 12-tone rows.
So they are also 26 · 6! series, for each n odd. Finally, there are
12! + 7 · 26 · 6!
fixed tone rows. Thus, we get |S| = (12! + 7t)/48 = 9985920 tone row
classes.
More generally, for the N -EDO, we have the following result,
which is a direct application of the de Bruijn extension of the Pólya
enumeration theorem.
Theorem 7.2. In the N -EDO (N ≥ 3), there are exactly N !
tone rows of N notes. Under the equivalence of the derived forms
(i.e. under the action of the dihedral group), there are
1 (N −4)/2 1 N + 1 ! if N is even
4 (N − 1)! + 2 N 2
1
(N −5)/2 N −1 !
4 (N − 1)! + 2 2 if N is odd
where ϕ(i) is Euler’s function and νn is defined for all the divisors
of 2n by
⎧ n/i
⎨ i (2n/i − 1)!!
⎪ if i is odd
νn (i) = ni 2n/i
⎪
⎩ ik (2k − 1)!! if i is even
k=0 2k
Under the action of the dihedral group, the number of Gauss diagrams
in the 2n-EDO is given by
1 1
dn = cn + (κn−1 + κn )
2 2
Serial Knots 181
with
2
n
n!
κn =
k!(n − 2k)!
k=0
With the Gauss order, the 554 serial chord diagrams can be
classified. The Gauss word can be written in a form slightly dif-
ferent by locating the chords of different types: the chord a connects
two consecutive points of the Gauss diagram, b connects the two
extreme points of three consecutive points of the circle, and so on,
until the type f which connects the two extreme points of seven con-
secutive points of the circle. Moreover, one uses the opposite (e.g.
e−1 is a chord linking the two extreme points of eight consecutive
points). For example, the series of Jean Barraqué Au-delà du hasard
087142t3e569 is represented by the word
X = abd−1 dab
The structural vector of a chord diagram is made of six compo-
nents which count the number of chords of each type a, b, c, d, e, f .
For instance, the structural vector of the previous chord diagram is
220200. This means that the diagram includes two chords of type a,
two chords of type b, none of type c, two chords of type d and none
of types e and f .
The intertwining properties of chords in chord diagrams, found in
the Gauss or structural word, are very important for the analysis of
musical serial structures, as well as indexing under the same chord
diagram.
Example 7.1. The chord diagram used by Karel Goeyvaerts
in his Sonata for two pianos (123679t540e8) has a Gauss word
121324546356, a structural word bcf bdc and a structural vector
022101.
A statistical study of 1132 12-tone rows (see Jedrzejewski, 2021)
has shown that the diagram most used by composers is the diagram
of the Symphony, op. 21 from Webern (5876t934012e).
Many tone rows do not belong to the four previous types. For
instance, the tone rows of Pierre Boulez’s Sonatine for flute and
piano (0e7184392t56), the Cycle de l’artisan furieux of Le Marteau
sans Maı̂tre (314596t78e20), Karel Goeyvaerts’s Opus 2 for 13
instruments (e65478329t10), Jacques Lenot’s Third Sonate for piano
(58726039t14e), Peter Schat’s Polonaise opus 2 (689e2t013457) and
Jan van Vlijmen’s Costruzione per due pianoforti (015473298t6e) do
not belong to the set T ∪ T ∪ R ∪ R .
7.5 Combinatoriality
d(X, Y ) = d(F X, F Y )
Consider the 12-tone row of Jean Barraqué’s Piano Sonata, writ-
ten in compact form X = 46e5102398t7 and the same slightly dis-
turbed tone row Y = 46e8102395t7. The tone row X is also written
in the form of a permutation
0 1 2 3 4 5 6 7 8 9 t e
X=
4 6 e 5 1 0 2 3 9 8 t 7
that we write in cycle notation
X = (0, 4, 1, 6, 2, e, 7, 3, 5 )(6, 9)
are different:
der(X, Y ) = card{n, xn = yn }
it is also the number of points where the n-th element sn of the tone
row XY −1 is different from n:
der(X, Y ) = card{n, sn = n}
The analytical study of the groups associated with all the string
diagrams leads to the following result.
7.8 Exercises
Exercise 7.1
Build all chord diagrams in Z6 . List their Gaussian words and their
intervallic structure. Show that the symmetry of the diagrams is
such that there is no difference between the action of the cyclic
group and the action of the dihedral group. Find this result using
194 A Compendium of Musical Mathematics
0123456789te
σ=
97823e t10645
= (096t4328)(17)(5e)
Prove that σ 8 (S) = S and draw the chord diagrams for each
σ k (S),for k = 1, . . . , 8. Draw the graph of the permutational evolu-
tion of the series used by Barraqué.
Serial Knots 195
Exercise 7.6
Let us consider two series of Elisabeth Lutyens’s Requiem: serie
A = 32e78014569t and serie B = 370e8465t912. Compute the Cayley
distance of the two series.
References
Combinatorial Designs
Δ = {d2 − d1 | d1 , d2 ∈ D, d1 = d2 }
k(k − 1) = v(v − 1)
197
198 A Compendium of Musical Mathematics
D1 ∪ D2 ∪ · · · ∪ Dr = G\{0}
Mp (x) = px mod v
M3 (D) = {6, 9, 2, 7} = D + 4
0 0 0 1 1 2 3
1 2 4 2 5 3 4
3 6 5 4 6 5 6
0 0 0 1 1 2 2
1 3 5 3 4 3 4
2 4 6 5 6 6 5
The dual of the Fano plane exchanges points and lines (see
Fig. 8.1).
z 2 = (k − λ)x2 + (−1)(v−1)/2 λy 2
and k = r,
det(N ) = k2 (k − λ)v−1
Example 8.9. The symmetric design (49, 16, 5) leads to the dio-
phantine equation z 2 = 11x2 + 5y 2 which (1, 1, 4) is a nonzero
solution.
Combinatorial Designs 203
Remark 8.2. It has been proved by Lam in 1989 that the converse
of the Bruck–Ryser–Chowla (BRC) theorem is false. The projec-
tive planes are symmetric designs with parameters (v, k, λ) =
(n2 + n + 1, n + 1, 1), where n is the order of the projective plane.
Lam shows that there are no projective plane of order 10. If there
were a projective plane of order n = 10, it would be a symmetric
design (111, 11, 1). Since 111 is odd, and the diophantine equation
z 2 = 10x2 − y 2 has the solution (1, 3, 1), the design (111, 11, 1) passes
the BRC test. But since there are no projective plane of order 10,
the BRC conditions are not sufficient to guarantee the existence of
symmetric designs.
this is not enough and he will often consider other relations, such as
the absence of common points or the movement of only one point of
each block.
The easy case is when blocks are constructed from a set of gen-
erators B = B | T1v (B) ≡ 1 with action of the cyclic group. These
designs are represented by cyclic representations. A special case are
the designs
m
p − 1 pm−1 − 1 pm−1 − 1
, ,
p−1 p−1 p−1
208 A Compendium of Musical Mathematics
Tom Johnson is the composer who has worked the most on block
designs. In 2009, he composed Septet II for two flutes, oboe, clarinet,
two violins and viola based on the Fano plane (7, 3, 1). In Vermont
Rhythms (2008), a 42 × 11 rhythm, he used the design (11, 6, 3) of
11 blocks.
Fifty-Five Chords (2008) for organ, 23 minutes of organ music
all derived from a (11, 4, 6) block design of 55 blocks. This last
piece raises the question of the graphical representation of a block
design. Johnson gives several graphs of the same design: Cosmo-
logical view in which every single chord has no notes in common
with exactly four chords; Pentagonal view where each chord has one
pair of notes in common with one chord, the other pair in common
with one other chord, and no notes in common with the adjacent
chords; Spider web view linking chords with three notes in com-
mon; and Starfish view where three pairs of notes combine to form
three chords (see Jedrzejewski and Johnson, 2014). On the wheel of
the graph drawn by the composer (see Fig. 8.6), the first block is
on the left opening (11, 10, 3, 2). Blocks are connected if they have
three notes in common, and move clockwise until one reach the
block (2, 9, 3, 10), which has also three notes in common with the
first one.
Block Design for piano (2005) is based on design 4-(12, 6, 10) of
330 blocks. It could be built using 30 base blocks and one automor-
phism that adds 1 to each point except 11 which remains always
the same. In the commentary he gives of his work, the composer
shows how the block design technique is for him as important as the
210 A Compendium of Musical Mathematics
8.7 Exercises
Exercise 8.1
Verify that D = {1, a, a3 , b, a2 b2 } is a difference set in
G = a, b | a7 = b3 = 1, ba = a2 b
Exercise 8.2
Verify that D = {a, a2 , a4 , b, b2 } is a difference set in
G = a, b | a7 = b3 = 1, ba = ab
Exercise 8.3
Verify that D = {1, 3, 13, 16, 17} is a (additive) difference set (21,5,1).
Find a multiplier.
Exercise 8.4
In the additive group Z21 , find a multiplier of the difference set
D = {1, 4, 5, 10, 12}.
Exercise 8.5
Let G = Z31 . Find a multiplier of the difference set D = {1, 5, 11,
24, 25, 27}.
Exercise 8.6
Prove that it cannot exist a symmetric design (67, 12, 2).
Exercise 8.7
Consider the block design (7,3,2) of 14 blocks:
1 1 1 1 1 1 2 2 2 2 3 3 3 3
2 2 3 4 5 6 3 4 5 6 4 4 5 5
3 4 5 6 7 7 7 5 6 7 6 7 6 7
Find a graphic representation.
References
Tuning Systems
213
214 A Compendium of Musical Mathematics
81/64 34 81
K= = 4 = ≈ 22 cents
5/4 2 ·5 80
P 531441/524288 32805 38 · 5
S= = = = 15 ≈ 2 cents
K 81/80 32768 2
CH = 21/53 ≈ 23 cents
216 A Compendium of Musical Mathematics
CF = 21/31 ≈ 39 cents
16/9 26 64
C7 = = 2 = ≈ 27 cents
7/4 3 ·7 63
The leimma is the difference between three octaves and five just
fifths:
23 28 256
L= 5
= 5
= ≈ 90 cents
(3/2) 3 243
The apotome is the difference between seven just fifths and four
octaves:
(3/2)7 37 2187
A= = = ≈ 114 cents
24 211 2048
ΔFi = ai K + bi S
In the same way, the difference between the major third of the tem-
perament and the pure major third (5/4) is denoted by ΔTM and
Tuning Systems 217
the difference between the minor third of the temperament and the
pure minor third (6/5) is denoted by ΔTm :
4 −6 s4 r4
4 3 2 k s a 5 K + b5 S (1 + s85 )K + r58 S −(1 + s42 )K − r24 S
11 −17 s11 r11
5 3 2 k s a12 K + b12 S −s4 K − (1 + r4 )S −(1 + s11
11 11 11
9 )K − r9 S
6 −9 s6 r6
6 3 2 k s a7 K + b7 S (1 + s7 )K + r7 S −(1 + s4 )K − r46 S
10 10 6
−1 s1 r1
7 3.2 k s a 2 K + b2 S (1 + s52 )K + r25 S s10 10
2 K + (1 + r2 )S
8 −12 s8 r8
8 3 2 k s a9 K + b9 S −s8 K − (1 + r8 )S −(1 + s6 )K − r68 S
8
3 −4 s3 r3
9 3 2 k s a 4 K + b4 S (1 + s74 )K + r47 S −(1 + s3 )K − r3 S
10 −15 s10 r10
10 3 2 k s a11 K + b11 S −s3 K − (1 + r3 )S −(1 + s10
10 10 10
8 )K − r8 S
5 −7 s5 r5 9 9
11 3 2 k s a 6 K + b6 S (1 + s6 )K + r6 S −(1 + s3 )K − r35 S
5
The coefficients verify the relations (in what follows, all sums are
considered modulo 12)
12
12
ai = bi = −1
i=1 i=1
and
b
n
rab = bi rn = r1n = bi
i=a i=1
a
12
12 j+6
md = 1 + 7 sj
12
j=1
218 A Compendium of Musical Mathematics
where a is defined by
37
a = 1200 log 2
213/12
about 14 cents.
Over the years, there have been several solution in pre-20th cen-
tury Europe to temper the fifths. Perhaps, the most successful tuning
was the meantone tuning which appeared sometime around the late
15th century. Pietro Aaron (ca 1485–1545) proposed to diminish the
fifths by a fourth of a syntonic comma except one.
In equal temperament, all fifths are diminished by the twelfth of
a Pythagorean comma.
Continued fractions appear, apparently for the first time, in the work
of Hindu mathematician Aryabhata in the late fifth century. A simple
continued fraction is a representation of a real number α given by
1
α = a0 +
1
a1 +
1
a2 +
a3 + · · ·
or in a compact abbreviated notation
α = [a0 ; a1 , a2 , . . .]
with a semicolon to separate the integer part. If the integer part is
null (a0 = 0), the notation α = [a1 , a2 , . . .] is also used. For typo-
graphic convenience, the continued fraction is also written as
1 1 1
a0 + ···
a1 + a2 + a3 +
It is well known that every real number α has a unique continued
fraction expansion. This expansion is finite if α is rational (α ∈ Q)
and infinite if α is irrational (α ∈ R\Q). Let α be the floor function
of α, the greatest integer less or equal α. Define the sequence xn by
the relations x0 = α and
1
xn+1 =
xn − xn
The numbers an are then the integer part of xn
an = xn
Define inductively the numbers pn and qn by the relations
p0 = a0, p1 = a0 a1 + 1, pn+1 = pn an+1 + pn−1 , n≥1
q0 = 1, q 1 = a1 , qn+1 = qn an+1 + qn−1 , n≥1
Definition 9.4. The convergents are the rational numbers cn
defined by
pn 1
cn = = [a0 ; a1 , . . . , an ] = a0 +
qn 1
a1 +
..
. 1
+
an
Tuning Systems 221
Proof.
(1) It is easy to see by induction that (qn ) is a strictly increasing
sequence, since q0 = 1 < q1 = a1 , and if qn−1 < qn , we have
qn+1 qn−1
= an+1 + > an+1 > 1
qn qn
(2) Eliminating an between the equations
pn+1 = pn an+1 + pn−1
qn+1 = qn an+1 + qn−1
leads to
pn+1 qn − pn qn+1 = −(pn qn−1 − pn−1 qn−1 )
and since
p1 q 0 − p0 q 1 = 1
we get
pn+1 qn − pn qn+1 = (−1)n
(3) Dividing the previous equation by qn qn+1 , we get
pn+1 pn (−1)n+1
− =
qn+1 qn qn qn+1
or
(−1)n+1
cn+1 − cn =
qn qn+1
Thus, the distance between cn and cn−2 is
pn pn−2
cn − cn−2 = −
qn qn−2
pn−1 pn−2 (−1)n−1
= − +
qn−1 qn−2 qn−1 qn
(−1)n−2 (−1)n−1
= +
qn−2 qn−1 qn−1 qn
(−1)n−1 1 1
= −
qn−1 qn qn−2
(−1)n an
=
qn−2 qn−1
Tuning Systems 223
Given a real number α, there are two ways to define a best approx-
imation of α. These two kinds of approximation are associated with
two distinct musical measurements and lead to two methods of scale
construction. All results of this section are demonstrated in Khinchin
and Ya (1964).
The best approximation of the second kind is also the best approx-
imation of the first kind, but the converse is false.
Example 9.4. To find q just fifths (3/2) lying in p octaves, one has
to solve the Diophantine equation with p, q = 0:
q
p 3
2 =
2
p0 1 (−1)n
= + + ··· + + ···
q0 q0 q1 qn qn+1
we get
2q = Qq = 5q/4
3a · 5b · 7c · · · n
{1, 8/7, 7/6, 6/5, 5/4, 4/3, 7/5, 10/7, 3/2, 8/5, 5/3, 12/7, 7/4}
{1, 10/9, 9/8, 8/7, 7/6, 6/5, 5/4, 9/7, 4/3, 7/5, 10/7, 3/2,
14/9, 8/5, 5/3, 12/7, 7/4, 16/9, 9/5}
{1, 12/11, 11/10, 10/9, 9/8, 8/7, 7/6, 6/5, 11/9, 5/4, 14/11, 9/7,
4/3, 11/8, 7/5, 10/7, 16/11, 3/2, 14/9, 11/7, 8/5, 18/11, 5/3,
12/7, 7/4, 16/9, 9/5, 20/11, 11/16}
9.10 Exercises
Exercise 9.1
Compute the first convergents of ln(6/5)/ ln(2) and show that
6 19
5 ≈ 25 . Evaluate the difference in cents.
Exercise 9.2
Compute the first convergents of ln(5/4)/ ln(2) and show that
5 28
4 ≈ 29 . Evaluate the difference in cents.
Exercise 9.3
Compute the tuning induced by the Farey sequence F6 .
Exercise 9.4
Compute the Euler–Fokker tuning 34 · 5.
Exercise 9.5
Study the cyclic tuning for the major third w = 5/4.
Exercise 9.6
Study the cyclic tuning for the minor third w = 6/5.
Exercise 9.7
Let G be the free group G = 2, 3, 5 and define the degree by
deg(2x 3y 5z ) = 7x + 11y + 16z
Compute the Kassel–Hellgouarch tuning. Find two kets |x, y, z solu-
tions of
7, 11, 16 |x, y, z = 0
and compute the associated interval 2x 3y 5z .
Exercise 9.8
Let us consider the degree
deg(2x 3y 7z ) = 7x + 11y + 20z
Compute the Kassel–Hellgouarch tuning. Find two kets |x, y, z
solutions of
7, 11, 20 |x, y, z = 0
and compute the associated interval 2x 3y 7z .
Tuning Systems 243
References
Solution 1.1
For N = 6, the divisors are {1, 2, 3, 6}. For |HA | = 6, we get the
set E1 = {0, 1, 2, 3, 4, 5}. For |HA | = 3, E2 = {0, 2, 4}. And for
|HA | = 2, E3 = {0, 3} and E4 = E3 ∪ T1 (E3 ) = {0, 1, 3, 4}. All others
combinations are transpositions of these sets.
Solution 1.2
For intervals (k = 2), 62 = 15 = 2 × 6 + 1 × 3. There are 1 interval
of limited transposition (the set {0, 3} of the previous exercise) with
3 transpositions, and 2 others
intervals with 6 transpositions.
For 3-chords (k = 3), 63 = 20 = 3 × 6 + 1 × 2. There are
four chords up to transpositions, one of which is a set of limited
transposition (the set {0, 1, 3, 4}).
Using the polynomial (1.1), we get
245
246 A Compendium of Musical Mathematics
Solution 1.3
For N = 8, the divisors are {1, 2, 4, 8}. For |HA | = 8, we get the
set E1 = {0, 1, 2, 3, 4, 5, 6, 7}. For |HA | = 4, E2 = {0, 2, 4, 6}. And
for |HA | = 2, E3 = {0, 4}, E4 = E3 ∪ T1 (E3 ) = {0, 1, 4, 5}, and
E5 = E3 ∪ T1 (E3 )∪ T2 (E3 ) = {0, 1, 2, 4, 5, 6}. All others combinations
are transpositions of these sets.
Solution 1.4
For intervals (k = 2), 82 = 28 = 3 × 8 + 1 × 4. There are 1 interval
of limited transposition (the set {0, 4} of the previous exercise) with
4 transpositions, and 3 others
8 intervals with 8 transpositions.
For 3-chords (k = 3), 3 = 56 = 7 × 8. There are 7 chords up to
transpositions and no set of limited transposition.
For 4-chords (k = 4), 84 = 70 = 8 × 8 + 1 × 2 + 1 × 4. There
are 10 chords up to transpositions, two of which are sets of limited
transposition (the set {0, 2, 4, 6} has 2 transpositions and the set
{0, 1, 4, 5} has 4 transpositions).
Using the polynomial (1.1), we get
k 1 2 3 4 5 6 7 8
k-chords 1 4 7 10 7 4 1 1
SLT 0 1 0 1 0 1 0 1
Solution 1.5
Solution 1.6
Solution 1.7
The inverse of M3 is M3 .
Solutions to Exercises 247
The M3 permutation is
0 1 2 3 4 5 6 7
M3 =
0 3 6 1 4 7 2 5
= (1, 3)(2, 6)(5, 7)
The order of the group is 32.
Solution 2.1
For two points (x, n) and (y, m) of Z12 ×Z2 , we introduce the interval
function
int ((x, n), (y, m)) = (2y + m) − (2x + n)
It is easy to verify that (Z12 × Z2 , Z24 , int) is a GIS.
Solution 2.2
Using the properties of the injection number, for all bijective func-
tion f
inj(X, Y )(f ) = inj(Y, X)(f −1 )
we already get the result since M3 is the inverse of M3 and M5 is the
inverse of M5 . Moreover,
inj(M3 (X), Y )(M3 ) = inj(X, Y )(M3 M3 ) = inj(X, Y )(id) = card(X ∩ Y )
and
inj(X, M5 (Y ))(M3 ) = inj(X, Y )(M5 M5 ) = inj(X, Y )(id) = card(X ∩ Y )
Solution 2.3
The proofs are the same as in the previous exercise. Using the prop-
erties of the injection number, for all bijective function f
inj(X, Y )(f ) = inj(Y, X)(f −1 )
we already get the result since M5 is the inverse of M7 . Moreover,
inj(M5 (X), Y )(M7 ) = inj(X, Y )(M7 M5 ) = inj(X, Y )(id) = card(X ∩ Y )
and
inj(X, M7 (Y ))(M7 ) = inj(X, Y )(M5 M7 ) = inj(X, Y )(id) = card(X ∩ Y )
248 A Compendium of Musical Mathematics
Solution 2.4
By definition,
ibang(X, Y )(n) = 1X (k)1Y (n + 5k)
k
Solution 2.5
We have
FTn (X) (x) = e−2iπkx/N = e−2iπ(j+n)x/N
k∈Tn (X) j∈X
= e−2iπnx/N e−2iπjx/N = e−2iπnx/N FX (x)
j∈X
Solution 2.6
By calculating the table of differences, we get
Since their interval vectors are equal, the two sets are homometric.
The Patterson polynomial of X is
Solution 2.7
Since X and Y are two complementary sets, they are homometric.
Their interval content is ic(X) = 433221. By multiplying by
5 mod 12, the sets M5 X = {0, 1, 3, 5, 6, 10} and M5 Y = {0, 3, 5,
8, 10, 11} are also homometric. Their interval content is equal to
ic(M5 X) = 233231.
Solution 3.1
By starting the cycle of fifths on F (x = 5) and retaining the first k
fifths, we obtain the desired scales. For some values of k, the scale
is not well formed. For instance, for k = 4, the scale 0, 2, 5, 7 is not
WF because they are two notes between 5 and 0 and three notes
between 0 and 7. The scales for k = 6, 8, 9, 10, 11 are also not WF.
250 A Compendium of Musical Mathematics
Solution 3.2
By computing the generic and specific intervals, we get
X, generic Specific Y, generic Specific Z, generic Specific
1 1, 2, 3 1 1, 2 1 1, 2
2 3, 4 2 3, 4 2 3, 4
3 4, 5, 6 3 4, 5, 6 3 5, 6
4 6, 7, 8 4 6, 7, 8 4 6, 7
5 8, 9 5 8, 9 5 8, 9
6 9, 10, 11 6 10, 11 6 10, 11
Solution 3.4
Transpositions of X lead to only one pattern M2-M2. Thus, the car-
dinality (of X = 3) is not equal to variety (= 1). Moreover, on the
circle of fifths (0, 2, 4, 6, 8, 10, 0), the structure 1-1-4 is different of the
multiplicity.
Solution 3.5
The transpositions of X are {0, 2, 4}, {2, 4, 7}, {4, 7, 9}, {7, 9, 0} and
{9, 0, 2}: one has structure 2-2-8, two have structure 2-3-7 and the
two others have structure 3-2-7. They are 3 patterns for 3 elements,
thus cardinality equals variety. On the circle of fifths (0, 7, 2, 9, 4, 0),
the structure 2-2-1 is equal to the multiplicity.
Solution 3.6
Using the generator m, we get 0, 8, 3, 11, 6, 1, 9, 4, 12, 7, 2, 10, 5. The
scale X is the 8 first values
X = {0, 1, 3, 4, 6, 8, 9, 11}
Solution 4.1
(1) Denote xM the major chord rooted on x, and xm the minor chord
rooted on x, we have
P (xM ) = xm = W0 (xM )
L(xM ) = (x + 4)m = W4 (xM )
R(xM ) = (x + 9)m = W9 (xM )
0M 4M 8M 0m 4m 8m
S0 0M 4M 8M 0m 4m 8m
S4 4M 8M 0M 8m 0m 4m
S8 8M 0M 4M 4m 8m 0m
W0 0m 4m 8m 0M 4M 8M
W4 4m 8m 0m 8M 0M 4M
W8 8m 0m 4m 4M 8M 0M
Solutions to Exercises 253
Solution 4.3
S is clearly the permutation associated with S1 , and W the permu-
tation of W0 . Thus, Wn = P Sn = W0 Sn = W S n . (2) P = W0 = W ,
L = W4 = W0 S4 = W S 4 , R = W9 = W S 9 . Consequently, since
Wi Wj (X) = Sj−i (X), we have
(RL)3 R = (W9 W4 )3 W9 = S73 W9 = S21 W9 = S9 W9 = W0 = P
using Si Sj (X) = Si+j (X) we see that Sn3 = S3n (mod 12 on the
indices).
Solution 4.4
DM
LP L / Bbm DM
J / Bbm
PL PL JQJ JQ
F m / Dm F m / Dm
LP L QJ
Solution 4.5
We have for the major triad
U V (xM ) = U (xm ) = (x + 1)m
V U (xM ) = V (x + 1)M = (x + 1)m
and for the minor triad
U V (xm ) = U (x + 6)M = (x + 7)M
V U (xm ) = V (x + 1)m = (x + 7)M
Therefore, U V = V U for all triads.
The permutations are
U = (1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12)
(13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24)
V = (1, 13, 7, 19)(2, 14, 8, 20)(3, 15, 9, 21)
(4, 6, 10, 22)(5, 17, 11, 23)(6, 18, 12, 24)
The computation leads to
U V = V U = (1, 14, 9, 22, 5, 18)(2, 15, 10, 23, 6, 19)
(3, 16, 11, 24, 7, 20)(4, 17, 12, 13, 8, 21)
254 A Compendium of Musical Mathematics
Solution 4.6
Verifying the relations, we get u4 (xM ) = xM and u4 (xm ) = xm .
Furthermore,
thus u2 (xM ) = v 2 (xM ). The same holds for minor triads u2 (xm ) =
v 2 (xm ). The last relation is also verified
Solution 5.2
We have
Solution 5.3
We have
Solution 5.4
By computing the distances and occurrences, we see that the rhythm
E(4, 9) = 2223 which is {0, 2, 4, 6, }9 is E-deep:
0 2 4 6
0 0 2 4 3 d 1 2 3 4
2 2 0 2 4
i 0 3 1 2
4 4 2 0 2
6 3 4 2 0
Solution 5.5
The fume-fume rhythm is isomorphic to the pentatonic scale
0-7-2-9-4. It is well formed and maximally even since they are two
c-distances for each generic interval:
0 2 4 7 9
d 1 2 3 4 5 6
0 0 2 4 5 3
i 0 3 2 1 4 0
2 2 0 2 5 5
4 4 2 0 3 5
c 1 2 3 4
7 5 5 3 0 2
d 2, 3 4, 5 7, 8 9, 10
9 3 5 5 2 0
Solution 5.6
(1) The height of ω = 3232323 is even and no cyclic shift of ω can
be factorized in two words uv such that h(u) = h(v). (2) The length
of ω is 2 + 1 with = 3. The height of the prefixes of length 3 (323,
332, 233, 232) is 8 or 7 for 232. (3) Let u = 323 and v = 2323, the
256 A Compendium of Musical Mathematics
Thus, the matrix belongs to SL2 (Z). The standard factorization with
ω1 = C(3, 1) = aaab = aub with u = aa and ω2 = C(5, 2) =
aaabaab = avb with v = aabaa can be rewritten as the palin-
dromic factorization ω = (ava, bab). Considering the multiples of
3 (or 8 reading from left to right) (3, 6, 9, 1, 4, 7, 10, 2, 5, 8, 0), we
get the black keys ω1 = 3, 6, 9, 1 = {1, 3, 6, 9} and the white keys
ω2 = 4, 7, 10, 2, 5, 8, 0 which is the diatonic scale:
Solution 6.1
We have
Since A(1) = 4 = A and Φ16 | A(x), conditions (T1 ) and (T2 ) are
satisfied. A tiles Z16 by CM theorem.
Solutions to Exercises 257
Solution 6.2
(1) Let us set
xN − 1
QN (x) = = 1 + x + x2 + · · · + xN −1
x−1
We have A(x)B(x) = QN (x), and C(x) = A(x)Qk (xN ). Thus,
C(x)B(x) = A(x)B(x)Qk (xN ) = QN (x)Qk (xN ) = QkN (x)
(2) This means that the cyclotomic factors of C are those of A, and
the cyclotomic polynomial Φd is such that d | kN and d N . Let
pα ∈ SA and q β ∈ SQkN , then since A tiles ZN p | N implies that
pα q β | kN . Thus, Φpα qβ | QN | C.
Solution 6.3
Clearly, the shape of A: 11 11 11 shows that A does not tile. More-
over,
A(x) = 1 + x + x3 + x4 + x6 + x7 = (1 + x)(1 + x3 + x6 ) = Φ2 (x)Φ9 (x)
A verifies (T1 ) since SA = {2, 32 } and A(1) = 6 = 2 × 3 = |A|. But
Φ2×9 = Φ18 does not divide A(x).
Solution 6.4
(1) Since B(x) = xn A(x), A(1) = B(1) = |A| = |B|. (2) For
every cyclotomic polynomial Φd (x), we have Φd (x) | A(x) ⇔
Φd (x) | B(x).
Solution 6.5
(1) We have B(x) = A(xn ), and (T1 ) follows since A(1) = B(1).
(2) Let SA = {q1α1 , . . . , qm
αm } be powers of distinct primes. From the
property,
p Φps (x) if p is a factor of s
Φs (x ) = (9.4)
Φs (x)Φps (x) if p is not a factor of s
we see that if s ∈ SB = SpA , then s = pqiαi if p is a factor of qiαi or
s = qiαi if p is not factor of qiαi . Consequently,
SB = {pα+1 , pα ∈ SA } ∪ {q β ∈ SA , q prime = p}
and A satisfies (T2 ) iff B satisfies (T2 ). The result is also true for
any n. The proof consists in iterating over each prime resulting from
the decomposition of n.
258 A Compendium of Musical Mathematics
Solution 7.1
Let’s label the circle clockwise from 0 to 5. The generalized tritone
is (0, 3). We write the chord diagrams by connecting the points of
3 elementary intervals. There exists 5 chord diagrams whose Gauss
words are as follows:
Solution 7.2
The group G has 24 elements
1 b b2 b3 b4 b5
b6 b−1 b−2 b−3 b−4 b−5
a ab ab2 ab3 ab4 ab5
c cb cb2 cb3 cb4 cb5
Solution 7.3
The series S = 03e4128t5976 is all-interval series. Moreover, the set
Ac = {5, 6, 7, 8, 9, t} of its second hexachord is the transposition a
tritone up (or down) T6 (Ac ) = A = {0, 1, 2, 3, 4, e} of its first hexa-
chord. Since T6 (A) = Ac , S ∈ T is a transpose tone row, and since
T9 (−A) = Ac , S ∈ R is a reverse series.
Solution 7.4
The series S is not an all-interval series. Furthermore, by listing the
transpositions and inversions of the first hexachord, we see that S ∈ /
T ∪ R ∪ T ∪ R .
Solutions to Exercises 259
Solution 7.5
We have
σ(S) = 901738425BA6 D111 σ 5 (S) = 34179A065B82 D10
σ 2 (S) = 69712038B54A D36 σ 6 (S) = 2371649AB508 D105
σ 3 (S) = A61789205B34 D105 σ 7 (S) = 8217A3645B90 D36
σ 4 (S) = 4A710689B523 D10 σ 8 (S) = 087142A3B569 D111
2 3 4
& & &
1 5
+ D111e D36 e D105e D10 n
8 7 6
Solution 7.6
We have
−1 0123456789t e
B =
2te057614983
and
−1 0123456789t e
AB =
e9t304128657
= (0e72t54)(196)
Solution 8.1
Compute all differences (written multiplicatively):
1 a a3 b a2 b2
1 1 a6 a4 b2 a3 b
a a 1 a5 ab2 a4 b
a3 a3 a2 1 a3 b2 a6 b
b b ba6 ba4 1 ba3 b
a2 b2 a2 b2 a2 b2 a6 a2 b2 a4 a2 b4 a2 b2 a3 b
260 A Compendium of Musical Mathematics
Solution 8.6
Applying the BRC theorem, we consider the equation z 2 = 10x2 −
2y 2 . It has no solution, since −2 = 8 mod 10. The design (67, 12, 2)
does not exist.
Solution 8.7
For instance, Tom Johnson finds the graphic representations (see
Fig. 9.1).
Solution 9.1
The first convergents are 1/2, 1/4, 5/19, 111/422, 449/1707, etc.
6 19
5 ≈ 31.948 and 25 = 32. This means that after 19 minor thirds,
we are in the vicinity of 5 octaves. The difference is
25 519
c= = , and 1200 log 2 (c) = 2.8 cents
(6/5)19 214 319
Solution 9.2
The first convergents are 1/3, 9/28, 19/59, 47/146, 207/643, etc.
5 28
4 ≈ 516.988 and 29 = 512. This means that after 28 major
thirds, we are in the vicinity of 9 octaves. The difference is
(5/4)28 528
c= = , and 1200 log 2 (c) = 16.8 cents
29 265
The next approximation 19/59 is better, and the comma is
219 2137
c= 59
= 59 , 1200 log 2 (c) = 7.5 cents
(5/4) 5
Solution 9.3
Solution 9.4
From the graph
5 3×5 32 × 5 33 × 5 34 × 5
1 3 32 33 34
Solutions to Exercises 263
Solution 9.5
#A Structure a b a (cents) b (cents)
3 aba 5/4 32/25 386 427
7 abababa 128/125 625/512 41 345
13 a ba3 ba3 ba2
2 128/125 510 /223 41 263
19 a3 ba5 ba5 ba3 128/125 516 /237 41 181
25 a4 ba7 ba7 ba4 128/125 522 /251 41 99
31 a5 ba9 ba9 ba5 128/125 528 /265 41 17
Solution 9.6
#A Structure a b a (cents) b (cents)
3 aba 6/5 25/18 316 569
7 abababa 125/108 648/625 253 63
11 aba2 baba2 ba 648/625 57 /25 37 63 191
19 a2 ba4 ba3 ba4 ba2 648/625 515 /211 315 63 65
23 a2 ba5 ba5 ba5 ba2 648/625 519 /214 319 63 3
Solution 9.7
The tuning is [1, 9/8, 5/4, 4/3, 3/2, 5/3, 9/5, 2]. The simplest solu-
tions of equation 7x + 11y + 16z = 0 are |−3, −1, 2 associated with
25/24 and |−7, 3, 1 associated with 135/128.
Solution 9.8
The tuning is [1, 8/7, 7/6, 4/3, 3/2, 12/7, 7/4, 2]. The kets |6, −2, −1
with comma 64/63 (27 cents) and |1, 3, −2 with comma 54/49
(168 cents) are solutions of 7x + 11y + 20z = 0.
This page intentionally left blank
Index
265
266 A Compendium of Musical Mathematics
N R
n-limit, 233 Rameau groups, 104
Narayana numbers, 143 Rameau–Schillinger operators, 99
negative isographies, 78 rational languages, 112
Netto’s theorem, 208 recognizable languages, 112
Noll, Thomas, xiii, 94, 109, 127 Reiner, David, 179
Nono, Luigi, 184, 190 resolvable designs, 204
rhythmic oddity, 119
O Riemann, Hugo, xii, 77
octatonic groups, 93 Riotte, André, 182
octatonic modes, 140 Rosenblatt, Joseph, 51
orbits, 17 Roux, Sébastien, 167
overtones, 213 RSK correspondence, 134
Russel, Georges, 1
P
Papadopoulos, Athanase, xiv, 213 S
parallel class, 204 S/W group, 86
Partch, Harry, 234 Sands, Arthur, xiii, 163
partials, 213 Savart, Félix, 214
Parzysz, Bernard, 10, 183 Saygun, Ahmed, 119
Patterson function, 49 scale, 3
pentatonic modes, 137 Schat, Peter, 185
perfect tilings, 160 Schensted, Craige, 133
periodic words, 123 Schönberg, Arnold, xiii, 151, 184, 187,
Pernazza, Ludovico, 171 192, 210
Piňos, Aloı̈s, xi, 16 Schützenberger, Marcel-Paul, 133
pitch class set, 2 semi-combinatorial, 185
pitch classes, 2 semiconvergents, 221
plactic modal class, 136 serial groups, 192
plactic monoid, 132
set class, 3
PLR group, 86
similarity measures, 187
Pólya’s theorem, 12, 21
Skalkottas, Nikos, 183
Popoff, Alexandre, 95
Slonimsky, Nicolas, 1
Posadas, Alberto, 114
Sós, Vera T., 234
positive isographies, 78
special Sturmian monoid, 129
primitive word, 110
spectral sets, 165
progressive transposition sets, 58
spectral unit, 52
Pythagoras, 229
stabilizer, 17
Pythagorean temperaments, 235, 239
star-free languages, 112
Pythagorean tunings, 228
Starr, Daniel, 16
Steedman, Marc, 113
Q
Steiner systems, 199
quarter tones, 10 Steiner triple systems, 199, 206
Index 269