Long-range correlation and multifractality in Bach’s Inventions
pitches
G. R. Jafari, P. Pedram
Department of physics, Shahid Beheshti University, Evin, Tehran 19839, Iran
Department of nano-science, IPM, P. O. Box 19395-5531, Tehran, Iran
L. Hedayatifar
Department of physics, Yazd University, Yazd, Iran
(Dated: October 28, 2018)
Abstract
We show that it can be considered some of Bach pitches series as a stochastic process with
scaling behavior. Using multifractal deterend fluctuation analysis (MF-DFA) method, frequency
series of Bach pitches have been analyzed. In this view we find same second moment exponents
(after double profiling) in ranges (1.7 − 1.8) in his works. Comparing MF-DFA results of original
series to those for shuffled and surrogate series we can distinguish multifractality due to long-range
correlations and a broad probability density function. Finally we determine the scaling exponents
and singularity spectrum. We conclude fat tail has more effect in its multifractality nature than
long-range correlations.
PACS numbers: 02.50.Fz, 05.45.Tp
http://arxiv.org/abs/0704.0726v2
Pythagoras knew it, but Bach demonstrated it: without mathematics there is no music [1].
I. INTRODUCTION
Many people think that, mathematics and music have some vague sort of affinity, but
most often supposed relationship between two fields turns out to be in details that are
not central to either. The mathematical proportions in musical works are hidden from
listener from old days. Thus it is forced to make use of interpretative techniques in order
to search for them, which is problematic from a methodological point of view. Some music
historians have very little time for numerology. In addition to mathematics being seen as
numerical symbolism, music is closely linked to absolute physical entities, such as frequency
and relation between intervals (an interval is space between two notes). Already in antiquity
this was seen as natural or cosmic premise on which music relied. It is illustrated fact that
not just musical notation, but also relationship between music and time has something to
do with mathematics and with one of the most significant transformations in music history.
Complexity of music can be especially attracting scientific interest. Among great variety
of complex and disordered systems most of music parameters such as frequency and pitch
(Pitch is the sound frequency of any given note.) [2, 3, 4, 5, 6], Amplitude or Dynamics
(Dynamics are the changes in volume during a musical piece.) [7, 8], Intervals (Intervals are
the distances between notes in the musical scale.), Rhythm (Rhythm is the structure of the
placement of notes in musical time.) can be consider as stochastic processes. Also, some
authors try to clustering the music [9, 10].
In all technicalities, music can be composed of notes. A note is a sign used in music
to represent the relative duration and pitch of sound. In traditional music theory pitch
classes are represented by the first seven letters of the Latin alphabet (A, B, C, D, E, F,
and G) or (Do - Re - Mi - Fa - Sol - La - Si), in Italian notation. Each note is assigned
a specific vertical position on a staff. Since the physical causes of music are vibrations of
mechanical systems, their frequencies are often measured in hertz (Hz). These frequencies
are mathematically related to each other, and are defined around the central note. The
current ”standard pitch” for this note is 440 Hz. Any note is exactly an integer number
of half-steps away from central note. Let this distance be denoted n. Then, the desired
FIG. 1: Typical sheet music, Invention no.1 by Bach.
frequency is given by
ν = 440× 2n/12 Hz. (1)
In this paper we would like to characterize complex behavior of frequency of note signal of
Bach’s Inventions and Sinfonias through computation of signal parameters scaling exponents,
which quantifies correlation exponents and multifractality of signal. Inventions and Sinfonias
are collection of short pieces which Bach wrote for musical education of Bach young pupils.
These are among finest examples of artistic gems ever written for this purpose, and probably
because of this, they became very popular among Bach’s pupils and others ever since they
were written. Inventions and Sinfonias contain two and three music voices respectively. The
number of the data in frequency series is dependent to the piece and will obtain from one
voice which contains 1000 notes in average. Because of non-stationary nature of frequency of
music series, and due to finiteness of available data sample, we should apply a methods which
are insensitive to non-stationarities, like trends. In order to separate trends from correlations
we need to eliminate trends in our frequency data. Several methods are used effectively for
this purpose: detrended fluctuation analysis (DFA) [11], rescaled range analysis (R/S) [12]
and wavelet techniques (WT) [13].
We use MF-DFA method for analysis and eliminating trend from data set. This method
are the modified version of DFA method to detect multifractal properties of time series.
DFA method introduced by Peng et al. [11] has became a widely used technique for the
determination of (mono-) fractal scaling properties and the detection of long-range correla-
tions in noisy, non-stationary time series [14, 15, 16, 17]. It has successfully been applied to
diverse fields such as DNA sequences [11, 18], heart rate dynamics [19, 20, 21], neuron spik-
ing [22], human gait [23], wind speed, [24], long-time weather records [25], cloud structure
[26], geology [27], ethnology [28], economical time series [29], and solid state physics [30].
One reason to employ DFA method is to avoid spurious detection of correlations that are
artefacts of non-stationarity in time series.
The focus of present paper is on intriguing statistical properties and multifractal nature
of frequency series. In Particular, Figures 1 and 2 show the score of Invention no.1 and the
frequency fluctuation of Invention no. 1 and Sinfoina no. 1, respectively. In general, two
different types of multifractality in frequency series can be distinguished: (i) Multifractality
due to a fatness of probability density function (PDF) of time series. In this case mul-
tifractality cannot be removed by shuffling the series. (ii) Multifractality due to different
long-range correlations in small and large scale fluctuations. In this case data may have a
PDF with finite moments, e.g. a Gaussian distribution. Thus corresponding shuffled time
series will exhibit mono-fractal scaling, since all long-range correlations are destroyed by
shuffling procedure. If both kinds of multifractality are present, shuffled series will show
weaker multifractality than original series.
The paper is organized as follows: In section II we describe MF-DFA methods in detail
and show that, scaling exponents determined by MF-DFA method are identical to those
obtained by standard multifractal formalism based on partition functions. In section III,
we analysis of frequency series of Bach’s Inventions also examine source of multifractality
in frequency data by comparison MF-DFA results for remaining data set to those obtained
by MF-DFA for shuffled and surrogate series. section IV closes with a conclusion.
II. MULTIFRACTAL DETRENDED FLUCTUATION ANALYSIS
The simplest type of multifractal analysis is based upon standard partition function multi-
fractal formalism, which has been developed for multifractal characterization of normalized,
stationary measurements [31, 32, 33, 34]. Unfortunately, this standard formalism does not
give us correct results for non-stationary time series that are affected by trends or those which
cannot be normalized. In the early 1990s an improved multifractal formalism wavelet trans-
form modulus maxima (WTMM) method [13], has been developed. This method is based
on wavelet analysis and involves tracing the maxima lines in continuous wavelet transform
100 200
100 200
FIG. 2: Up to down: typical frequency series of Invention and Sinfonias no.1 by Bach.
over all scales. Other method like, multifractal detrended fluctuation analysis (MF-DFA), is
based on identification of scaling of qth-order moment depending on signal length, and it is
generalization of standard DFA method in which q = 2. In contrast to WTMM method MF-
DFA does not require modulus maxima procedure, and hence does not require more effort in
programming and computing time than conventional DFA. On the other hand, one should
find correct scaling behavior of fluctuations, from experimental data which often affected by
non-stationary like trends. This have to be well distinguished from intrinsic fluctuations of
system. In addition often in collected data we do not know reasons, or even worse scales, for
underlying trends, and also, usually available record data is small. So for reliable detection
of correlations, it is essential to distinguish trends for intrinsic fluctuations from collected
data. Hurst rescaled-range analysis [12] and other non-detrending methods work well when
records are long and do not involve trends, otherwise it might give wrong results. DFA is a
well established method for determining scaling behavior of noisy data, where data presence
of trends and we don’t know their origin and shape [11, 20, 35, 36, 37].
A. Description of MF-DFA method
Modified multifractal DFA (MF-DFA) procedure consists of five steps. The first three
steps are essentially identical to conventional DFA procedure (see e.g. [11, 14, 15, 16, 17]).
Suppose that xk is a series of length N , and it is of compact support, i.e. xk = 0 for an
insignificant fraction of the values only.
• Step 1: Determine the “profile”
Y (i) ≡
[xk − 〈x〉] , i = 1, . . . , N. (2)
Subtraction of the mean 〈x〉 from xk is not compulsory, since it would be eliminated by later
detrending in third step.
• Step 2: Divide profile Y (i) into Ns ≡ int(N/s) non overlapping segments of equal lengths
s. Since length N of series is often not a multiple of considered time scale s, a short part at
the end of profile may remain. In order not to disregard this part of series, same procedure
should be repeated starting from the opposite end. Thereby, 2Ns segments are obtained
altogether.
• Step 3: Calculate local trend for each of 2Ns segments by a least-square fit of series. Then
determine the variance
F 2(s, ν) ≡
{Y [(ν − 1)s+ i]− yν(i)}
, (3)
for each segment ν = 1, . . . , Ns and
F 2(s, ν) ≡
{Y [N − (ν −Ns)s+ i]− yν(i)}
, (4)
for ν = Ns + 1, . . . , 2Ns. Where, yν(i) is fitted polynomial in segment ν. Linear, quadratic,
cubic, or higher order polynomials can be used in the fitting procedure (conventionally
called DFA1, DFA2, DFA3, . . ., DFAm) [11, 21]. Since detrending of time series is done by
subtraction of fited polynomial from profile, different order DFA differ in their capability of
eliminating trends in series. In (MF-)DFAm trend of order m in profile (and equivalently,
order m − 1 in original series) are eliminated. Thus a comparison of results for different
orders of DFA allows one to estimate type of the polynomial trend in time series [15, 16].
• Step 4: Average over all segments to obtain q-th order fluctuation function, defined by:
Fq(s) ≡
F 2(s, ν)
, (5)
where, in general, variable q can take any real value except zero. q = 2, standard DFA
procedure is retrieved. Generally we are interested to know how generalized q dependent
fluctuation functions Fq(s) depend on time scale s for different values of q. Hence, we must
repeat steps 2, 3 and 4 for several time scales s. It is apparent that Fq(s) will increase with
increasing s. Of course, Fq(s) depends on DFA order m. By construction, Fq(s) is only
defined for s ≥ m+ 2.
• Step 5: Determine scaling behavior of fluctuation functions by analyzing log-log plots of
Fq(s) versus s for each value of q. If series xi are long-range power law correlated, then Fq(s)
, for large values of s, increases as a power-law i.e.,
Fq(s) ∼ s
h(q). (6)
In general, exponent h(q) may depend on q. For stationary time series such as, fractional
Gaussian noise (fGn), Y (i) in eq. 2, will have a fractional Brownian motion (fBm) signal,
so, 0 < h(q = 2) < 1.0. The exponent h(2) is identical to well known Hurst exponent H
[11, 14, 31]. Also for non-stationary signal, such as fBm noise, Y (i) in eq. 2, will be a sum
of fBm signal, so corresponding scaling exponent of Fq(s) is identified by h(q = 2) > 1.0
[11, 38]. For monofractal time series, h(q) is independent of q, since scaling behavior of
variance F 2(s, ν) is identical for all segments ν, and averaging procedure in eq.(5) will just
give this identical scaling behavior for all values of q. If we consider positive values of q, the
segments ν with large variance F 2(s, ν) (i.e. large deviation from the corresponding fit) will
dominate average Fq(s). Thus, for positive values of q, h(q) describes scaling behavior of
segments with large fluctuations. For negative values of q, segments ν with small variance
F 2(s, ν) will dominate average Fq(s). Hence, for negative values of q, h(q) describes scaling
behavior of segments with small fluctuations.
B. Relation to standard multifractal analysis
For a stationary, normalized series multifractal scaling exponent h(q) defined in eq.(6)
are directly related to scaling exponent τ(q) defined by standard partition function based
on multifractal formalism as shown below. Suppose that series xk of length N is a sta-
tionary, normalized sequence, then detrending procedure in step 3 of MF-DFA method is
not required, since no trend has to be eliminated. Thus, DFA can be replaced by standard
fluctuation analysis (FA), which is identical to DFA except definition of variance, which is
simplified for each segment ν = 1, . . . , Ns. In step 3 eq.(3) now becomes :
F 2FA(s, ν) ≡ [Y (νs)− Y ((ν − 1)s)]
2. (7)
Inserting this simplified definition into eq.(5) and using eq.(6), we obtain
|Y (νs)− Y ((ν − 1)s)|q
∼ sh(q). (8)
For simplicity we can assume that length N of series is an integer multiple of scale s,
obtaining Ns = N/s and therefore
|Y (νs)− Y ((ν − 1)s)|q ∼ sqh(q)−1. (9)
This corresponds to the multifractal formalism used e.g. in [32, 34]. In fact, a hierarchy of
exponents Hq similar to our h(q) has been introduced based on eq.(9) in [32]. In order to
relate this to standard textbook box counting formalism [31, 33], we employ definition of
profile in eq.(2). It is evident that term Y (νs)− Y ((ν − 1)s) in eq.(9) is identical to sum of
numbers xk within each segment ν of size s. This sum is known as box probability ps(ν) in
standard multifractal formalism for normalized series xk,
ps(ν) ≡
k=(ν−1)s+1
xk = Y (νs)− Y ((ν − 1)s). (10)
The scaling exponent τ(q) is usually defined via partition function Zq(s),
Zq(s) ≡
|ps(ν)|
q ∼ sτ(q), (11)
where q is a real as in MF-DFA method, discussed above. Using eq.(10) we see that eq.(11)
is identical to eq.(9), and one can obtain analytical relation between two sets of multifractal
scaling exponents,
τ(q) = qh(q)− 1. (12)
Thus, we see that h(q) defined in eq.(6) for MF-DFA is directly related to classical mul-
tifractal scaling exponents τ(q). Note that, h(q) is different from generalized multifractal
dimension
D(q) ≡
q − 1
qh(q)− 1
q − 1
, (13)
50 100
FIG. 3: The log-log plot F (s) versus s for q = 2.0 after double profiling for Invention no. 1.
that are used instead of τ(q) in some papers. In this case, while h(q) is independent of q for a
monofractal time series, D(q) depends on q. Another way to characterize a multifractal series
is looking to singularity spectrum f(α), which is related to τ(q) via a Legendre transform
[31, 33],
α = τ ′(q) and f(α) = qα− τ(q). (14)
Here, α is the singularity strength or Hölder exponent, while f(α) denotes dimension of
subset of series that is characterized by α. Using eq.(12), we can directly relate α and f(α)
to h(q),
α = h(q) + qh′(q) and f(α) = q[α− h(q)] + 1. (15)
Hölder exponent denotes monofractality, while in multifractal case, different parts of struc-
ture are characterized by different values of α, leading to existence of spectrum f(α).
III. ANALYSIS OF MUSIC FREQUENCY SERIES
As mentioned in section II, a spurious of correlations may be detected if time series is non-
stationarity, so direct calculation of correlation behavior, spectral density exponent, fractal
dimensions etc., don’t give us a reliable results. It can be checked out that, frequency series
is non-stationary. One can verified non-stationarity property experimentally by measuring
stability of average and variance in a moving window for example using scale s. According
to MF-DFA1 method, generalized Hurst exponents h(q) in eq.(6) can be found by analyzing
-20 -10 0 10 20
-20 -10 0 10 20
2 2.5
FIG. 4: The q dependence of the exponents h∗(q), τ∗(q) and singularity spectrum f∗(α∗), after
double profiling, are shown in the upper to lower panels respectively for Invention no. 1 .
log-log plots of Fq(s) versus s for each q. Hurst exponent is between 0 < H < 1. However,
the MF-DFA method can only determine positive generalized Hurst exponents, in order to
refine the analysis near the fGn/fBm boundary or strongly anti-correlated signals when is
close to zero. The simplest way to analyze such data is to integrate the time series before the
MF-DFA procedure. Hence, we replace the single summation in Eq. (1), which is describing
10 20 30
FIG. 5: The second moment exponents, h∗(q = 2), after double profiling, of various Inventions
and Sinfonias in their number ordering.
the profile from the original data, by a double summation by using Classification by the
signal summation conversion method (SSC). After using SSC method, fGn switch to fBm
and fBm switch to sum-fBm. In this case the relation between the new exponent, h∗(q = 2),
and h(q = 2) is h(q = 2) = h∗(q = 2)− 1 [38, 39, 40] (recently Movahed et al. have proven
the relation between derived exponent from double profile of series in DFA method and
h(q = 2) exponent in Appendix of [39]). We find h∗(q = 2) = 1.74 ± 0.03 of Invention no.
1 by using SSC method. Note that, the exponents of the new series, after double profiling
(SSC), are different from before SSC. Therefore we name these new exponents as h∗(q),
τ ∗(q) and f ∗(α∗).
Results using MF-DFA1 method for frequency signal are shown in Figure 4, which shows
that frequency series is a multifractal process as indicated by strong q dependence of the
exponents h∗(q = 2) and τ ∗(q)[40]. q dependence of the multifractal scaling exponent τ ∗(q)
has linear behaviors for q > 0 and q < 0 and slopes of τ ∗(q) are 1.53± 0.01 and 2.61± 0.01,
respectively. The values which derived for quantities of MF-DFA1 method for Invention no.
1 are given in Table I. We have calculated exponent h∗(q) for other Inventions and Sinfonias
as well, and all of them are in the 1.7− 1.8 range (Fig. 5). Figure . 4(c) shows the width of
singularity spectrum, f ∗(α∗), i.e. ∆α = α∗(qmax)− α
∗(qmin) for the series is approximately,
1.70. It value shows the power of multifractality of interevents is very strong [41]. The values
of derived quantities from MF-DFA1 method, are given in Table I. There are some specific
lengths in music that have important effects in music’s statistical parameters. a) Rhythm,
b) Scale is a set of tones each having a definite pitch (perceived fundamental frequency of a
sound) and each having a specific frequency ratio compared to the others, that is, each having
specific interval relative to all the other pitches, arranged in a sequence from low to high,
or alternatively, from high to low. Composers often transform musical patterns by moving
every note in the pattern by a constant number of scale intervals. Since the intervals of a
scale can have various sizes, this process introduces subtle melodic and harmonic variation
into the music. This variation is what gives scalar music much of its complexity.), c) Bar or
measure is a segment of time defined as a given number of beats (the basic time unit of a
piece) of a given duration. More over in the Persian music (Eastern music), there are others
characteristic lengths such as Gushe which is based on the change of the tonic or stop notes.
Existence of these characteristic lengths can effect on music complexity. The music style
and more importantly composers can use various the characteristic lengths in the music. We
check a few pieces of other genders of music like jazz and Persian traditional music (Eastern
music) that obtained h∗(2) in ranges, (1.75−1.9) and (1.3−1.7), respectively. We find more
characteristic lengths in the Persian traditional music which leads to decrease the exponents.
Usually, in MF-DFA method, deviation from a straight line in log-log plot of eq.(6) occurs
for small scales s. This deviation limits capability of DFA to determine correct correlation
behavior for very short scales and in regime of small s. The modified MF-DFA is defined as
follows [15]:
Fmodq (s) = Fq(s)
〈[F shufq (s
′)]2〉1/2 s1/2
〈[F shufq (s)]
2〉1/2 s′1/2
(for s′ ≫ 1),
where 〈[F shufq (s)]
2〉1/2 denotes the usual MF-DFA fluctuation function, defined in eq.(5),
averaged over several configurations of shuffled data taken from original time series, and
s′ ≈ N/40. The values of h∗(q) exponent obtained by modified MF-DFA1 methods for
frequency series time series is 1.70±0.03. The relative deviation of the new exponent which
is obtained by modified MF-DFA1 in comparison to MF-DFA1 for original data is less than
Now, we are interested in to determine source of multifractality. In general, two different
-20 -10 0 10 20
Origenal
Shuffeld
Surrogate
FIG. 6: The new exponent, h∗(q), after double profiling, as a function of q for the original,
surrogate and shuffled data.
types of multifractality in time series can be distinguished: (i) multifractality due to a
fatness of probability density function (PDF) of the time series. In this case multifractality
cannot be removed by shuffling series. (ii) multifractality due to different correlations in
small and large scale fluctuations. In this case data may have a PDF with finite moments,
e.g. a Gaussian distribution. Thus corresponding shuffled time series will exhibit mono-
fractal scaling, since all long-range correlations are destroyed by shuffling procedure. If
both kinds of multifractality are present, shuffled series will show weaker multifractality
than original series. The easiest way to clarify the type of multifractality, is by analyzing
corresponding shuffled and surrogate time series. Shuffling of time series destroys long-range
correlation, therefore if multifractality belongs only to long-range correlation, we should find
h∗shuf(q) = 1.5. The multifractality nature due to fatness of PDF signals is not affected by
shuffling procedure. On the other hand, to determine multifractality due to broadness of
PDF, phase of discrete fourier transform (DFT) coefficients of frequency series time series are
replaced with a set of pseudo independent distributed uniform (−π, π) quantities in surrogate
method. The correlations in surrogate series do not change, but probability function changes
to Gaussian distribution. If multifractality in time series is due to a broad PDF, h∗(q)
obtained by surrogate method will be independent of q. If both kinds of multifractality
are present in frequency series time series, shuffled and surrogate series will show weaker
multifractality than original one.
To check nature of multifractality, we compare fluctuation function Fq(s), for origi-
nal series (after cancelation of sinusoidal trend) with result from corresponding shuffled,
F shufq (s) and surrogate series F
q (s). Differences between these two fluctuation functions
with original one, directly indicate presence of long-range correlations or broadness of prob-
ability density function in original series. These differences can be observed in ratio plot
of Fq(s)/F
q (s) and Fq(s)/F
q (s) with respect s [40]. Since anomalous scaling due to a
broad probability density affects both of Fq(s) and F
q (s) in same way, only multifractality
due to correlations will be observed in Fq(s)/F
q (s). Scaling behavior of these ratios are
Fq(s)/F
q (s) ∼ s
h∗(q)−h∗
(q) = sh
cor(q), (17)
Fq(s)/F
q (s) ∼ s
h∗(q)−h∗
(q) = sh
(q). (18)
If only fatness of PDF is responsible for multifractality, one should have h∗(q) = h∗shuf(q)
and h∗cor(q) = 0. On the other hand, deviations from h
cor(q) = 0 indicates existence of
correlations, and q dependence of h∗cor(q) indicates that multifractality is due to long rage
correlation. If only correlation multifractality is present, one finds h∗shuf(q) = 1.5. If both
distribution and correlation multifractality are present, both of h∗shuf(q) and h
sur(q) will
depend on q. The q dependence of exponent h∗(q) for original, surrogate and shuffled time
series are shown in Figures 6. q dependence of hcor and h
PDF shows that multifractality
nature of frequency series time series is due to both broadness of PDF and long-range
correlation. Absolute value of h∗PDF(q) is greater than h
cor(q), so multifractality due to
correlation is weaker than mulifractality due to fatness.
Deviation of h∗sur(q) and h
shuf(q) from h
∗(q) can be determined by using χ2 test as follows:
[h∗(qi)− h
(qi)]
σ(qi)2 + σ⋄(qi)2
, (19)
symbol “ ⋄ ” can be replaced by “sur” and “shuf”, to determine the confidence level of h∗sur
and h∗shuf to the new exponents, h
∗(q), of original series, respectively. Reduced chi-square
χ2ν⋄ =
(N is number of degree of freedom) for shuffled and surrogate time series are
16.07, 21.30, respectively. Width of singularity spectrum, f ∗(α∗), i.e. ∆α = α∗(qmin) −
α∗(qmax) for original, surrogate and shuffled time series are approximately, 1.04, 0.65 and
0.92 respectively. These values conclude that multifractality due to fat tail is dominant[41].
TABLE I: Values of h∗(q = 2), τ∗(q = 2) exponents and Width of singularity spectrum, f∗(α∗),
∆α for q = 2.0 of Invention no.1 obtained by MF-DFA1.
h∗(2) τ∗(2) ∆α
Original 1.74 ± 0.03 2.48 ± 0.03 1.70
Surrogate 1.71 ± 0.03 2.44 ± 0.03 0.80
Shuffled 0.50 ± 0.03 0.02 ± 0.03 0.50
Values of the exponents, h∗(q = 2.0), τ ∗(q = 2) and width of singularity spectrum,
f ∗(α∗), ∆α for the original, shuffled and surrogate of frequency series obtained with MF-
DFA1 method are reported in Table I.
IV. CONCLUSION
MF-DFA method allows us to determine multifractal characterization of non-stationary
and stationary time series. We have shown that MF-DFA1 result of frequency series of
Bach Inventions. Applying MF-DFA1 method is demonstrated that frequency series have
long term correlation. We calculated the second moment exponent after the new profile for
other Inventions and Sinfonias and they are in the 1.7 − 1.8 range. q dependence of h∗(q)
and τ ∗(q), shows that frequency series has multifractal behavior. By comparing the second
moment exponent of original time series with shuffled and surrogate one’s, we have found
that multifractality due to broadness of probability density function has more contribution
than correlation in Inventions.
V. ACKNOWLEDGMENT
We would like to thank M Sadegh Movahed for reading the manuscript and useful com-
ments. GRJ would like to acknowledge the hospitality extended during his visits at the
IPAM, UCLA, where this work was started.
[1] http://www.ntnu.no/
http://www.ntnu.no/
[2] Julyan H. E. Cartwright, Diego L. Gonzalez, and Oreste Piro, 1999 Phys. Rev. Lett. 82, 5389
[3] Heather D. Jennings, Plamen Ch. Ivanov, Allan de M. Martins, P.C. da Silva, Viswanathan
G M, Physica A 336 (2004) 585 594.
[4] Jafari G R, Pedram P, Ghafori K, AIP Proceeding (2006).
[5] Gonzalez D L, Morettini L, Sportolari F, Rosso O, Cartwright J H E and Piro O,
arXiv:chao-dyn/9505001, (1995).
[6] YU SHI, arXiv:adap-org/9509001, (1995).
[7] Diodati p and Piazz P, Eur. Phys. J. B 17, 143-145.
[8] Jean Pierre Boon and Olivier Decroly, 1995, Chaos 5(3) 501-508.
[9] Rudi Cilibrasi, Paul Vitanyi and Ronald de Wolf, Computer Music Journal, 28:4, pp. 4967,
Winter 2004.
[10] Rudi Cilibrasi and Paul Vitanyi, 2005, IEEE TRANSACTIONS ON INFORMATION THE-
ORY, 51(4), 15231545.
[11] Peng C K, Buldyrev S V, Havlin S, Simons M, Stanley H E, and Goldberger A L, 1994 Phys.
Rev. E 49, 1685 ; Ossadnik S M, Buldyrev S B, Goldberger A L, Havlin S, Mantegna R N,
Peng C K, Simons M and Stanley H E, 1994 Biophys. J. 67, 64
[12] Hurst H E, Black R P and Simaika Y M, 1965 Long-term storage. An experimental study
(Constable, London)
[13] Muzy J F, Bacry E and Arneodo A, 1991 Phys. Rev. Lett. 67, 3515
[14] Taqqu M S, Teverovsky V and Willinger W, 1995 Fractals 3, 785
[15] Kantelhardt J W, Koscielny-Bunde E, Rego H H A, Havlin S and Bunde A, 2001 Physica A
295, 441
[16] Hu K, Ivanov P Ch, Chen Z, Carpena P and Stanley H E, 2001 Phys. Rev. E 64, 011114
[17] Chen Z, Ivanov P Ch, Hu K and Stanley H E, 2002 Phys. Rev. E 65, preprint physics/041107.
[18] Buldyrev S V, Goldberger A L, Havlin S, Mantegna R N, Matsa M E, Peng C K, Simons M
and Stanley H E, 1995 Phys. Rev. E 51, 5084; Buldyrev S V, Dokholyan N V, Goldberger A
L, Havlin S, Peng C K, Stanley H E and Viswanathan G M, 1998 Physica A 249, 430
[19] Ivanov P Ch, Bunde A, Amaral L A N, Havlin S, Fritsch-Yelle J, Baevsky R M, Stanley H E
and Goldberger A L, 1999 Europhys. Lett. 48, 594; Ashkenazy Y, Lewkowicz M, Levitan J,
Havlin S, Saermark K, Moelgaard H, Thomsen P E B, Moller M, Hintze U and Huikuri H V,
2001 Europhys. Lett. 53, 709; Ashkenazy Y, Ivanov P Ch, Havlin S, Peng C K, Goldberger
http://arxiv.org/abs/chao-dyn/9505001
http://arxiv.org/abs/adap-org/9509001
A L and Stanley H E, 2001 Phys. Rev. Lett. 86, 1900
[20] Peng C K, Havlin S, Stanley H E and Goldberger A L, 1995 Chaos 5 82
[21] Bunde A, Havlin S, Kantelhardt J W, Penzel T, Peter J H and Voigt K, 2000 Phys. Rev. Lett.
85, 3736
[22] Blesic S, Milosevic S, Stratimirovic D and Ljubisavljevic M, 1999 Physica A 268, 275; Bahar
S, Kantelhardt J W, Neiman A, Rego H H A, Russell D F, Wilkens L, Bunde A and Moss F,
2001 Europhys. Lett. 56, 454
[23] Hausdorff J M, Mitchell S L, Firtion R, Peng C K, Cudkowicz M E, Wei J Y and Goldberger
A L, 1997 J. Appl. Physiology 82, 262
[24] Govindon R B and Kantz H, Europhysics letters, 68 (2), pp: 184-190 (2004)
[25] Koscielny-Bunde E, Bunde A, Havlin S, Roman H E, Goldreich Y and Schellnhuber H J, 1998
Phys. Rev. Lett. 81, 729; Ivanova K and Ausloos M, 1999 Physica A 274, 349; Talkner P and
Weber R O, 2000 Phys. Rev. E 62, 150
[26] Ivanova K, Ausloos M, Clothiaux E E and Ackerman T P, 2000 Europhys. Lett. 52, 40
[27] Malamud B D and Turcotte D L, 1999 J. Stat. Plan. Infer. 80, 173
[28] Alados C L and Huffman M A, 2000 Ethnology 106, 105
[29] Mantegna R N and Stanley H E, 2000 An Introduction to Econophysics (Cambridge University
Press, Cambridge); Liu Y, Gopikrishnan P, Cizeau P, Meyer M, Peng C K and Stanley H E,
1999 Phys. Rev. E 60, 1390; Vandewalle N, Ausloos M and Boveroux P, 1999 Physica A 269,
[30] Kantelhardt J W, Berkovits R, Havlin S and Bunde A, 1999 Physica A 266, 461; Vandewalle
N, Ausloos M, Houssa M, Mertens P W and Heyns M M, 1999 Appl. Phys. Lett. 74, 1579
[31] Feder J, 1988 Fractals (Plenum Press, New York)
[32] Barabási A L and Vicsek T, 1991 Phys. Rev. A 44, 2730
[33] Peitgen H O, Jürgens H and Saupe D, 1992 Chaos and Fractals (Springer-Verlag, New York),
Appendix B
[34] Bacry E, Delour J and Muzy J F, 2001 Phys. Rev. E 64, 026103
[35] Fano U, 1947 Phys. Rev. 72 26
[36] Barmes J A and Allan D W, 1996 Proc. IEEE 54 176
[37] Buldyrev S V, Goldberger A L, Havlin S, Mantegna R N, Matsa M E, Peng C K, Simons M,
Stanley H E, 1995 Phys. Rev. E 51 5084
[38] Eke A, Herman P, Kocsis L and Kozak L R, 2002 Physiol. Meas. 23, R1-R38
[39] M Sadegh Movahed, G R Jafari, F Ghasemi , Sohrab Rahvar and M Reza Rahimi Tabar, J.
Stat. Mech. (2006) P02003.
[40] Kantelhardt J W, Zschiegner S A, Kosciliny-Bunde E, Bunde A, Pavlin S and Stanley H E,
2002 Physica A 316, 78-114.
[41] Oświȩcimka P and et. al, [arXive:cond-mat/0504608]
http://arxiv.org/abs/cond-mat/0504608
	Introduction
	Multifractal Detrended Fluctuation Analysis
	Description of MF-DFA method
	Relation to standard multifractal analysis
	Analysis of music frequency series
	Conclusion
	Acknowledgment
	References

  We show that it can be considered some of Bach pitches series as a stochastic
process with scaling behavior. Using multifractal deterend fluctuation analysis
(MF-DFA) method, frequency series of Bach pitches have been analyzed. In this
view we find same second moment exponents (after double profiling) in ranges
(1.7-1.8) in his works. Comparing MF-DFA results of original series to those
for shuffled and surrogate series we can distinguish multifractality due to
long-range correlations and a broad probability density function. Finally we
determine the scaling exponents and singularity spectrum. We conclude fat tail
has more effect in its multifractality nature than long-range correlations.

Introduction to Econophysics (Cambridge University
Press, Cambridge); Liu Y, Gopikrishnan P, Cizeau P, Meyer M, Peng C K and Stanley H E,
1999 Phys. Rev. E 60, 1390; Vandewalle N, Ausloos M and Boveroux P, 1999 Physica A 269,
[30] Kantelhardt J W, Berkovits R, Havlin S and Bunde A, 1999 Physica A 266, 461; Vandewalle
N, Ausloos M, Houssa M, Mertens P W and Heyns M M, 1999 Appl. Phys. Lett. 74, 1579
[31] Feder J, 1988 Fractals (Plenum Press, New York)
[32] Barabási A L and Vicsek T, 1991 Phys. Rev. A 44, 2730
[33] Peitgen H O, Jürgens H and Saupe D, 1992 Chaos and Fractals (Springer-Verlag, New York),
Appendix B
[34] Bacry E, Delour J and Muzy J F, 2001 Phys. Rev. E 64, 026103
[35] Fano U, 1947 Phys. Rev. 72 26
[36] Barmes J A and Allan D W, 1996 Proc. IEEE 54 176
[37] Buldyrev S V, Goldberger A L, Havlin S, Mantegna R N, Matsa M E, Peng C K, Simons M,
Stanley H E, 1995 Phys. Rev. E 51 5084
[38] Eke A, Herman P, Kocsis L and Kozak L R, 2002 Physiol. Meas. 23, R1-R38
[39] M Sadegh Movahed, G R Jafari, F Ghasemi , Sohrab Rahvar and M Reza Rahimi Tabar, J.
Stat. Mech. (2006) P02003.
[40] Kantelhardt J W, Zschiegner S A, Kosciliny-Bunde E, Bunde A, Pavlin S and Stanley H E,
2002 Physica A 316, 78-114.
[41] Oświȩcimka P and et. al, [arXive:cond-mat/0504608]
http://arxiv.org/abs/cond-mat/0504608
	Introduction
	Multifractal Detrended Fluctuation Analysis
	Description of MF-DFA method
	Relation to standard multifractal analysis
	Analysis of music frequency series
	Conclusion
	Acknowledgment
	References
