Fractal Geometry - Notes
Fractal Geometry - Notes
Introduction
There is no hard and fast definition of a fractal set, but just a list of characteristic properties.
Suppose that 𝐹 is a fractal set, then 𝐹 has some (or all) of the following properties:
(i) 𝐹 has fine structure, that is, detail on arbitrary small scales;
(ii) 𝐹 is too irregular to be described in traditional geometrical language, both locally and
globally;
(iii) often 𝐹 has some form of self-similarity, perhaps approximate or statistical;
(iv) usually the fractal dimension of 𝐹 (defined in some way) is greater than its topological
dimension;
(v) in most cases of interest, 𝐹 is defined in a very simple way, perhaps recursively.
Fractal Geometry - Notes
Mathematical Background
Basic Set Theory
Theorem A set 𝐴 is closed if, and only if:
(i) its complement is open;
(ii) it is equal to its closure, i.e. 𝐴 = 𝐴;
(iii) it contains all of its limit points.
Theorem A set 𝐴 is open if, and only if, it is equal to its interior, i.e. 𝐴 = int 𝐴.
Theorem The union of any collection of open sets is open, as is the intersection of any finite
number of open sets.
Theorem The intersection of any collection of closed sets is closed, as is the union of any
finite number of closed sets.
Theorem A point 𝑥 ∈ ℝ𝑛 is a limit point of a set 𝐴 if, and only if, for every 𝑟 > 0 the open
ball 𝐵 o (𝑥, 𝑟) contains points of 𝐴 other than 𝑥.
Theorem A point 𝑥 ∈ ℝ𝑛 is a boundary point of a set 𝐴 if, and only if, the ball 𝐵(𝑥, 𝑟)
intersects both 𝐴 and its complement for all 𝑟 > 0.
Theorem The diameter of a set 𝐴 is equal to the diameter of its closure, i.e. |𝐴| = |𝐴|.
Theorem Let 𝐴 ⊆ ℝ𝑛 , then 𝐴 is compact if, and only if, 𝐴 is closed and bounded.
Theorem The intersection of any collection of compact sets is compact, as is the union of
any finite number of compact sets.
Theorem Suppose that 𝐴1 ⊃ 𝐴2 ⊃ ⋯ is a decreasing sequence of non-empty compact sets.
Then the intersection 𝐴 = ⋂∞ 𝑖=1 𝐴𝑖 is non-empty. Moreover, if 𝐴 is contained in some open set
𝑉, then the finite intersection ⋂𝑘𝑖=1 𝐴𝑖 is contained in 𝑉 for some 𝑘.
Theorem Suppose that 𝐴 and 𝐵 are compact subsets of ℝ𝑛 . Then 𝐴 and 𝐵 are a positive
distance apart, that is
inf{|𝑥 − 𝑦|: 𝑥 ∈ 𝐴, 𝑦 ∈ 𝐵} > 0.
Theorem Let 𝑓, 𝑔: ℝ𝑛 → ℝ be functions. If there exists 𝛿 > 0 such that 𝑓(𝑥) ≤ 𝑔(𝑥) for 𝑥 ∈
(𝑎, 𝛿), then
lim 𝑓(𝑥) ≤ lim 𝑔(𝑥) and lim 𝑓(𝑥) ≤ lim 𝑔(𝑥).
𝑥→𝑎 𝑥→𝑎 𝑥→𝑎 𝑥→𝑎
Note that, even if 𝑓(𝑥) < 𝑔(𝑥), the upper and lower limits could still be equal.
Holder and Lipschitz functions are continuous functions.
Congruences, similarities and affine transformations on ℝ𝑛 are examples of
homeomorphisms.
Theorem Let 𝑋 ⊆ ℝ𝑛 and 𝑌 ⊆ ℝ𝑚 . If 𝑓: 𝑋 → 𝑌 is a continuous function and 𝐴 is an open
subset of 𝑌, then 𝑓 −1 [𝐴] is an open subset of 𝑋.
Theorem Let 𝑓: 𝑋 → 𝑌 be continuous on 𝑋 and let 𝐴 be a compact subset of 𝑋. Then 𝑓[𝐴] is
also compact.
Mean Value Theorem Let the function 𝑓: ℝ → ℝ be continuous on the closed, finite interval
[𝑎, 𝑏] and differentiable on the open interval (𝑎, 𝑏). Then there exists a point 𝑐 in the open
interval (𝑎, 𝑏) such that
𝑓(𝑏) − 𝑓(𝑎)
𝑓 ′ (𝑐) = .
𝑏−𝑎
Geometrically this states the there exists a point 𝑐 ∈ (𝑎, 𝑏) for which the tangent to the graph
of the function 𝑓 at (𝑐, 𝑓(𝑐)) is parallel to the chord determined by the points (𝑎, 𝑓(𝑎)) and
(𝑏, 𝑓(𝑏)).
Theorem Let {𝑓𝑛 } be a sequence of functions defined on a common domain 𝐷 ⊂ ℝ. If {𝑓𝑛 }
converges pointwise to 𝑓 on 𝐷, then {𝑓𝑛 } converges uniformly to 𝑓 on 𝐷.
Theorem Let {𝑓𝑛 } be a sequence of functions defined on a common domain 𝐷 ⊂ ℝ, and let
𝑥0 ∈ 𝐷. If the sequence {𝑓𝑛 } converges uniformly to some function 𝑓 on 𝐷, and if each of the
functions 𝑓𝑛 is continuous at 𝑥0 , then the function 𝑓 is continuous at 𝑥0 . In particular, if each of
the functions 𝑓𝑛 is continuous on 𝐷, then 𝑓 is continuous on 𝐷.
𝜇 (⋃ 𝐴𝑛 ) ≤ ∑ 𝜇(𝐴𝑛 )
𝑛=1 𝑛=1
with equality
∞ ∞
𝜇 (⋃ 𝐴𝑛 ) = ∑ 𝜇(𝐴𝑛 )
𝑛=1 𝑛=1
𝜇 (⋃ 𝐴𝑛 ) = lim 𝜇(𝐴𝑛 ).
𝑛→∞
𝑛
A simple extension of this is that if, for 𝛿 > 0, 𝐴δ are Borel sets that are increasing as 𝛿
decreases, that is 𝐴𝛿′ ⊂ 𝐴δ for 0 < 𝛿 < 𝛿 ′ , then
𝜇 (⋃ 𝐴𝑛 ) = lim 𝜇(𝐴𝛿 ).
𝛿→0
𝛿>0
𝐸∞ = ⋂ 𝐸𝑘 .
𝑘=1
Fractal Geometry - Notes
Box-counting Dimension
A dimension indicates, in some way, how much space a set occupies near to each of its points.
Fundamental to most definitions of dimension is the idea of ‘measurement of a set at scale 𝛿’.
For each 𝛿 we measure a set in a way that detects irregularities of size 𝛿, and we see how
these measurements behave as 𝛿 → 0.
The following are desirable properties of a dimension;
Monotonicity - if 𝐸 ⊂ 𝐹, then dim 𝐸 ≤ dim 𝐹;
Range of values - if 𝐹 ⊂ ℝ𝑛 , then 0 ≤ dim 𝐹 ≤ 𝑛;
Stability - dim(𝐸 ∪ 𝐹) = max(dim 𝐸 , dim 𝐹);
Countable stability - dim(⋃∞
𝑖=1 𝐹𝑖 ) = sup𝑖 (dim 𝐹𝑖 );
Divider Dimension
Divider dimension Let 𝐶 be a curve and let 𝛿 > 0. Furthermore, let 𝑀𝛿 (𝐶) be the maximum
number of points 𝑥0 , 𝑥1 , … , 𝑥𝑚 on the curve, in that order, such that |𝑥𝑘 − 𝑥𝑘−1 | = 𝛿 for 𝑘 =
1,2, … , 𝑚. Then, the divider dimension is defined as
log 𝑀𝛿 (𝐶)
lim
𝛿→0 − log 𝛿
assuming that the limit exists (otherwise we may define upper and lower divider dimensions
using upper and lower limits).
The quantity (𝑀𝛿 (𝐶) − 1)𝛿 may be thought of as the ‘length’ of the curve 𝐶 measured using a
pair of dividers with points set at a distance 𝛿 apart. Hence the name ‘divider dimension’.
The divider dimension of a line is one.
The divider dimension can be estimated as the gradient of a log-log graph plotted over a
suitable range of values of 𝛿.
Fractal Geometry - Notes
Box-counting Dimension
Given a subset 𝐹 of ℝ𝑛 , for each 𝛿 > 0 we find the smallest number of sets of diameter at most
𝛿 that can cover the set 𝐹 and call this number 𝑁𝛿 (𝐹), indicating the number of ‘clumps’ of
size about 𝛿 into which 𝐹 may be divided. The dimension of 𝐹 reflects the way in which 𝑁𝛿 (𝐹)
grows as 𝛿 → 0. If 𝑁𝛿 (𝐹) obeys, at least approximately, a power law
𝑁𝛿 (𝐹) ≃ 𝑐𝛿 −𝑠
for positive constants 𝑐 and 𝑠, we say that 𝐹 has box dimension 𝑠.
Taking logarithms gives
log 𝑁𝛿 (𝐹) log 𝑐
log 𝑁𝛿 (𝐹) ≃ log 𝑐 − log 𝛿 ⇒ 𝑠≃ + .
− log 𝛿 log 𝛿
We hope to obtain 𝑠 as
log 𝑁𝛿 (𝐹)
𝑠 = lim
𝛿→0 − log 𝛿
and
log 𝑁𝛿 (𝐹)
dimB 𝐹 = lim
𝛿→0 − log 𝛿
and
log 𝑁𝛿𝑘 (𝐹)
dimB 𝐹 = lim .
𝑘→∞ − log 𝛿𝑘
These features are not desirable for a dimension, and severely limit the usefulness of box
dimension.
Fractal Geometry - Notes
Hausdorff Measure and Dimension
Hausdorff measure Suppose that 𝐹 is a subset of ℝ𝑛 and 𝑠 ≥ 0 is a real number. For 𝛿 > 0
define
∞
for all 𝛿-covers {𝑈𝑖 } of 𝐹. So, given an arbitrary 𝜀 > 0, we can find a 𝛿-cover of 𝐹 such that
∞
ℋ (⋃ 𝐹𝑖 ) ≤ ∑ ℋ 𝑠 (𝐹𝑖 )
𝑠
𝑖=1 𝑖=1
with equality
∞ ∞
ℋ 𝑠 (⋃ 𝐹𝑖 ) = ∑ ℋ 𝑠 (𝐹𝑖 )
𝑖=1 𝑛=1
increasing with 𝑠, so ℋ 𝑠 (𝐹) is also non-increasing. Furthermore, if 𝑡 > 𝑠 and {𝑈𝑖 } is a 𝛿-cover
of 𝐹, then
∞ ∞ ∞
0 𝑠
0 dimH 𝐹 𝑛
If 𝑡 > 𝑠, then
𝑁𝛿 (𝐹)𝛿 𝑡 = 𝛿 𝑡−𝑠 𝑁𝛿 (𝐹)𝛿 𝑠 .
Taking limits as 𝛿 → 0 gives
Fractal Geometry - Notes
lim 𝑁𝛿 (𝐹)𝛿 𝑡 = (lim 𝛿 𝑡−𝑠 ) (lim 𝑁𝛿 (𝐹)𝛿 𝑠 ).
𝛿→0 𝛿→0 𝛿→0
The value of 𝑠 at which lim 𝑁𝛿 (𝐹)𝛿 𝑠 jumps from ∞ to 0 is equal to the box dimension of 𝐹.
𝛿→0
It follows that there exists a 𝛿0 < 1 such that, for 0 < 𝛿 < 𝛿0
log(𝑁𝛿 (𝐹)) + 𝑠 log 𝛿 log(𝑁𝛿 (𝐹)) 𝑠 log 𝛿 log(𝑁𝛿 (𝐹))
= + >0 ⇒ > 𝑠.
− log 𝛿 − log 𝛿 − log 𝛿 − log 𝛿
Hence
log(𝑁𝛿 (𝐹))
dimB 𝐹 = lim ≥ 𝑠.
𝛿→0 − log 𝛿
Furthermore, if lim 𝑁𝛿 (𝐹)𝛿 𝑠 = 0, then taking logs gives
𝛿→0
It follows that there exists a 𝛿1 < 1 such that, for 0 < 𝛿 < 𝛿1
log(𝑁𝛿 (𝐹)) + 𝑠 log 𝛿 log(𝑁𝛿 (𝐹)) 𝑠 log 𝛿 log(𝑁𝛿 (𝐹))
= + <0 ⇒ < 𝑠.
− log 𝛿 − log 𝛿 − log 𝛿 − log 𝛿
Hence
log(𝑁𝛿 (𝐹))
dimB 𝐹 = lim ≤ 𝑠.
𝛿→0 − log 𝛿
This gives a definition of box dimension that looks very similar to the definition of Hausdorff
dimension. For sets that are 'reasonably regular' these two dimensions turn out to be equal. If,
however, a set can be covered more efficiently by sets of very different sizes (rather than by
sets of exactly the same size), then the expressions ℋ𝛿𝑡 (𝐹) and 𝑁𝛿 (𝐹)𝛿 𝑠 can be very different,
which leads to differing values for the Hausdorff and box dimensions.
Proposition 3.4 For every non-empty bounded 𝐹 ⊂ ℝ𝑛
dimH 𝐹 ≤ dimB 𝐹 ≤ dimB 𝐹.
A fundamental property of dimension is that Hausdorff, lower box and upper box dimensions
are invariant under bi-Lipschitz transformations. Thus, if two sets have different dimensions
there cannot be a Lipschitz mapping from one onto the other. Therefore, one approach to
fractal geometry is to regard two sets as equivalent if there is a bi-Lipschitz mapping between
them.
Proposition 3.5 Every set 𝐹 ⊂ ℝ𝑛 with dimH 𝐹 < 1 is totally disconnected.
Fractal Geometry - Notes
Calculation of Hausdorff Dimension
Typically, Hausdorff dimension calculations involve finding an upper, and a lower estimate for
dimH 𝐹, which will hopefully give the same values. Each of these estimates usually involves a
geometric observation followed by a calculation.
Calculation of Hausdorff measures and dimension can be rather complicated, even for fairly
simple fractal sets. For the upper estimate, there are often natural coverings that can be used.
It is the lower estimate that is usually more difficult since, according to the definition of
Hausdorff measure, all possible 𝛿-covers have to be considered to obtain the infimum value.
There is a simple 'heuristic' method that often gives the correct result when considering the
Hausdorff dimension of a self-similar fractal set. If 𝐹 = ⋃∞ 𝑖=1 𝐹𝑖 , where each 𝐹𝑖 is geometrically
similar to 𝐹, but scaled by a factor 𝑟𝑖 , then, provided that the 𝐹𝑖 are disjoint Borel sets, we can
write
∞ ∞ ∞
∑ 𝑟𝑖𝑠 = 1,
𝑖=1
Equivalent Definitions
The definition of Hausdorff measure involves coverings by arbitrary sets of diameter at most
𝛿. However, there are some restricted classes of covering sets that lead to the same values of
Hausdorff dimension.
The following covering sets give the same values of Hausdorff dimension:
arbitrary sets of diameter at most 𝛿;
balls of diameter at most 𝛿;
open, or closed sets of diameter at most 𝛿.
Binary interval A binary interval is an interval of the form [𝑟2−𝑘 , (𝑟 + 1)2−𝑘 ] where 𝑘 =
0,1,2, … and 𝑟 = 0,1, … , 2𝑘 − 1.
Net measures Let 𝐹 be a subset of the interval [0,1]. Define
∞
Then
ℳ 𝑠 (𝐹) = lim ℳ𝛿𝑠 (𝐹)
𝛿→0
Basic Methods
It is often relatively simple to obtain an upper bound for the dimension of a set. However, it is
usually much more difficult to find a lower bound for the dimension of a set.
To obtain an upper bound for the Hausdorff dimension of a set 𝐹 it is enough to evaluate sums
of the form ∑|𝑈𝑖 |𝑠 for specific coverings of 𝐹.
Proposition 4.1 Suppose that a set 𝐹can be covered by 𝑛𝑘 sets of diameter at most 𝛿𝑘 with
𝛿𝑘 → 0 as 𝑘 → ∞. Then
log 𝑛𝑘
dimH 𝐹 ≤ dimB 𝐹 ≤ lim .
𝑘→∞ − log 𝛿𝑘
Moreover, if 𝑛𝑘 𝛿𝑘𝑠 remains bounded as 𝑘 → ∞, then ℋ 𝑠 (𝐹) < ∞. If 𝛿𝑘 → 0 but 𝛿𝑘+1 ≥ 𝑐𝛿𝑘 for
some 0 < 𝑐 < 1, then
log 𝑛𝑘
dimB 𝐹 ≤ lim .
𝑘→∞ − log 𝛿𝑘
To obtain a lower bound for the Hausdorff dimension of a set 𝐹 we must show that sums of
the form ∑|𝑈𝑖 |𝑠 are greater than some positive constant for all 𝛿-coverings of 𝐹.
Proposition 4.2 (mass distribution principle) Let 𝜇 be a mass distribution on a set 𝐹, and
suppose that for some 𝑠 > 0 there exist numbers 𝑐 > 0 and 𝜀 > 0 such that
𝜇(𝑈) ≤ 𝑐|𝑈|𝑠
for all sets 𝑈 with |𝑈| ≤ 𝜀. Then
𝜇(𝐹)
ℋ 𝑠 (𝐹) ≥
𝑐
and
𝑠 ≤ dimH 𝐹 ≤ dimB 𝐹 ≤ dimB 𝐹.
The following theorem is useful when working with Proposition 4.1 and the mass distribution
principle.
Theorem Let 𝑎 and 𝑏 be positive real numbers, then
1
𝑎log 𝑏⁄log 𝑎 = 𝑏 and 𝑎log 𝑏⁄− log 𝑎 = .
𝑏
𝐹 = ⋂ 𝐸𝑘
𝑘=0
Fractal Geometry - Notes
is a totally disconnected subset of [0,1] that is generally a fractal.
The above general construction of fractal subsets of ℝ may be thought of as a generalisation of
the middle third Cantor set.
If the intervals in 𝐸𝑘 may have different lengths, it is not possible to use Proposition 4.1 to
obtain an upper bound for dimH 𝐹. The following generalisation of proposition 4.1 may be
used instead.
Generalisation of Proposition 4.1 Suppose that 𝐸𝑘 is a cover of 𝐹 by sets 𝑈𝑖 of diameter at
most 𝛿𝑘 with 𝛿𝑘 → 0 as 𝑘 → ∞. If there exists 0 < 𝑀 < ∞ such that
∑ |𝑈𝑖 |𝑠 ≤ 𝑀
𝑈𝑖 ∈𝐸𝑘
𝐹 = ⋂ 𝐻𝑘 ,
𝑘=0
then dimH 𝐹 = 𝑠.
Jarnik's theorem Suppose that 𝛼 > 2. Let 𝐹 be the set of real numbers 𝑥 ∈ [0,1] for which
the inequality
‖𝑞𝑥‖ ≤ 𝑞1−𝛼
is satisfied by infinitely many positive integers 𝑞. Then
2
dimH 𝐹 = .
𝛼
It should be noted that ‖𝑞𝑥‖ ≤ 𝑞1−𝛼 if, and only if
𝑛 𝑛
𝑥 ∈ [ − 𝑞 −𝛼 , + 𝑞 −𝛼 ]
𝑞 𝑞
for some integer 𝑛.
Fractal Geometry - Notes
Fractals Defined by Transformations
Many fractals are made up of parts that are, in some way, similar to the whole. These self-
similarities are not only properties of the fractals: they may actually be used to define them.
Iterated function systems do this in a unified way and, moreover, often lead to a simple way of
finding dimensions.
Iterated function system A finite family of contractions {𝑆1 , 𝑆2 , … , 𝑆𝑚 } with 𝑚 ≥ 2 is called
an iterated function system, or IFS (technically, this is the definition of a hyperbolic, or
contracting IFS).
Invariant set (attractor) Let 𝐷 be a closed subset of ℝ𝑛 . A nonempty compact subset 𝐹 of
𝐷 is an invariant set, or attractor, for the iterated function system {𝑆1 , 𝑆2 , … , 𝑆𝑚 } if, and only if
𝑚
𝐹 = ⋃ 𝑆𝑖 (𝐹).
𝑖=1
Theorem 9.1 Consider the iterated function system given by the contractions {𝑆1 , 𝑆2 , … , 𝑆𝑚 }
on 𝐷, which is a closed subset of ℝ𝑛 , so that
|𝑆𝑖 (𝑥) − 𝑆𝑖 (𝑦)| ≤ 𝑐𝑖 |𝑥 − 𝑦| (𝑥, 𝑦 ∈ 𝐷)
with 𝑐𝑖 < 1 for each 𝑖. Then there is a unique attractor 𝐹, i.e. a non-empty compact subset 𝐹 of
𝐷 such that
𝑚
𝐹 = ⋃ 𝑆𝑖 (𝐹).
𝑖=1
𝑆(𝐸) = ⋃ 𝑆𝑖 (𝐸)
𝑖=1
for 𝐸 ∈ 𝒮, and write 𝑆 𝑘 for the 𝑘th iterate of 𝑆 (so 𝑆 0 = 𝐸 and 𝑆 𝑘 (𝐸) = 𝑆(𝑆 𝑘−1 (𝐸)) for 𝑘 ≥ 1),
then
∞
𝐹 = ⋂ 𝑆 𝑘 (𝐸).
𝑘=0
𝐹 = ⋃ 𝑆𝑖 (𝐹)
𝑖=1
However, there may be other sets satisfying this equation which are not compact.
There are two main problems that arise in connection with iterated function systems:
(i) to determine an IFS for a given unique attractor;
(ii) to determine a unique attractor for a given IFS.
Fractal Geometry - Notes
Finding an IFS that has a given 𝐹 as its unique attractor can often be done by inspection, at
least if 𝐹 is self-similar or self-affine. For example, the Cantor set is easily seen to be the
attractor of the two similarities which give the basic self-similarities of the set
1 1 2
𝑆1 (𝑥) = 𝑥 and 𝑆1 (𝑥) = 𝑥 + .
3 3 3
However, there is no general method that can be applied to any 𝐹.
The transformation 𝑆 introduced in Theorem 9.1 is the key to computing the attractor of an
IFS.
The sequence of iterates 𝑆 𝑘 (𝐸) converges to the attractor 𝐹 for any initial set 𝐸 ∈ 𝒮 in the
sense that the Hausdorff metric 𝑑(𝑆 𝑘 (𝐸), 𝐹) → 0. Thus the 𝑆 𝑘 (𝐸) provide increasingly good
approximations to 𝐹. If 𝐹 is a fractal, these approximations may be called pre-fractals for 𝐹.
For each 𝑘, we have
where the union is over the set 𝐼𝑘 of all 𝑘-term sequences (𝑖1 , 𝑖2 , … , 𝑖𝑘 ). If 𝑆𝑖 (𝐸) is contained in
𝐸 for each 𝑖, and 𝑥 is a point of 𝐹, it follows from the fact that
∞
𝐹 = ⋂ 𝑆 𝑘 (𝐸)
𝑘=0
that there is a (not necessarily unique - see point 10 in the course notes) sequence (𝑖1 , 𝑖2 , … 𝑖𝑘 )
such that
𝑥 ∈ 𝑆𝑖1 ∘ 𝑆𝑖2 ∘ ⋯ ∘ 𝑆𝑖𝑘 (𝐸)
for all 𝑘. This sequence provides a natural coding for 𝑥, with
∞
so that 𝐹 = ⋃{𝑥𝑖1 ,𝑖2 ,… }. This expression for 𝑥𝑖1 ,𝑖2 ,… is independent of 𝐸 provided that 𝑆𝑖 (𝐸) is
contained in 𝐸 for all 𝑖.
For example, the Cantor set has the following IFS
1 1 2
𝑆1 (𝑥) = 𝑥 and 𝑆1 (𝑥) = 𝑥 + .
3 3 3
Now, if 𝑥𝑖1 ,𝑖2 ,… is a point in the Cantor set with base-3 expansion 0. 𝑎1 𝑎2 ⋯, then
𝑎𝑘 = 0 if 𝑖𝑘 = 1
and
𝑎𝑘 = 2 if 𝑖𝑘 = 2.
Note that, if the union
𝑚
𝐹 = ⋃ 𝑆𝑖 (𝐹)
𝑖=1
is disjoint, then 𝐹 must be totally disconnected (provided that the 𝑆𝑖 are injections).
Fractal Geometry - Notes
The pre-fractals 𝑆 𝑘 (𝐸) provide the usual construction of many fractals for a suitably chosen
initial set 𝐸, in which case, the 𝑆𝑖1 ∘ 𝑆𝑖2 ∘ ⋯ ∘ 𝑆𝑖𝑘 (𝐸) are called the 𝑘-level sets of the
construction.
𝑉 ⊃ ⋃ 𝑆𝑖 (𝑉)
𝑖=1
𝐹 = ⋃ 𝑆𝑖 (𝐹)
𝑖=1
∑ 𝑐𝑖𝑠 = 1.
𝑖=1
Some Variations
The calculations underlying Theorem 9.3 may be adapted to estimate the dimension of the
attractor of an IFS consisting of contractions that are not similarities.
Proposition 9.6 Let 𝐹 be the attractor of an IFS consisting of contractions 𝑆1 , 𝑆2 , … , 𝑆𝑚 on a
closed subset 𝐷 of ℝ𝑛 such that
|𝑆𝑖 (𝑥) − 𝑆𝑖 (𝑦)| ≤ 𝑟𝑖 |𝑥 − 𝑦| (𝑥, 𝑦 ∈ 𝐷)
with 0 < 𝑟𝑖 < 1 for each 𝑖. Then dimH 𝐹 ≤ dimB 𝐹 ≤ 𝑠, where 𝑠 is given by
Fractal Geometry - Notes
𝑚
∑ 𝑟𝑖𝑠 = 1.
𝑖=1
𝐹 = ⋃ 𝑆𝑖 (𝐹)
𝑖=1
with this union disjoint. Then 𝐹 is totally disconnected, and dimH 𝐹 ≥ 𝑠, where 𝑠 is given by
𝑚
∑ 𝑟𝑖𝑠 = 1.
𝑖=1
Note that the union ⋃𝑚 𝑖=1 𝑆𝑖 (𝐹) will certainly be disjoint if there is some non-empty compact
set 𝐸 with 𝑆𝑖 (𝐸) ⊂ 𝐸 for all 𝑖 and with the 𝑆𝑖 (𝐸) disjoint.
Dimensions of Graphs
Proposition 11.1 Let 𝑓: [0,1] → ℝ be a continuous function. Let 0 < 𝛿 < 1, and let 𝑚 be the
least integer greater or equal to 1⁄𝛿 . Then, if 𝑁𝛿 is the number of squares of the 𝛿-mesh that
intersect graph 𝑓
𝑚−1 𝑚−1
1 1
∑ 𝑅𝑓 [𝑖𝛿, (𝑖 + 1)𝛿] ≤ 𝑁𝛿 ≤ 2𝑚 + ∑ 𝑅𝑓 [𝑖𝛿, (𝑖 + 1)𝛿].
𝛿 𝛿
𝑖=0 𝑖=0
Self-affine Curves
Affine transformation A transformation defined by 𝑆(𝑥) = 𝑇(𝑥) + 𝑎 where 𝑇 is a non-
singular linear transformation and 𝑎 is a vector in ℝ𝑛 is called an affine transformation, or
affinity.
An affinity may be thought of as a shearing transformation; its contracting or expanding effect
need not be the same in every direction.
Let {𝑆1 , 𝑆2 , … , 𝑆𝑚 } be affine transformations represented in matrix notation with respect to
(𝑡, 𝑥) coordinates by
𝑡 1⁄𝑚 0 𝑡 (𝑖 − 1)⁄𝑚
𝑆𝑖 [ ] = [ ][ ] + [ ]
𝑥 𝑎𝑖 𝑟𝑖 𝑥 𝑏𝑖
that is
𝑡+𝑖−1
𝑆𝑖 (𝑡, 𝑥) = ( , 𝑎𝑖 𝑡 + 𝑟𝑖 𝑥 + 𝑏𝑖 ).
𝑚
Then, for each 𝑖, 𝑆𝑖 transforms vertical lines {(𝑡0 , 𝑥): 𝑥 ∈ ℝ} to vertical lines {(𝑡0 +
(𝑖 − 1)⁄𝑚 , 𝑥 ′ ): 𝑥 ′ ∈ ℝ}, with the vertical strip 0 ≤ 𝑡 ≤ 1 mapped onto the vertical strip
(𝑖 − 1)⁄𝑚 ≤ 𝑡 ≤ 𝑖 ⁄𝑚. We suppose that
1
≤ 𝑟𝑖 < 1
𝑚
for each 𝑖 so that contractions in the 𝑡 direction are stronger than in the 𝑥 direction.
Now
0 1⁄𝑚 0 0 0 0
𝑆1 [ ⁄(1 ] = [ ] [ ⁄(1 ] + [ ] = [ ]
𝑏1 − 𝑟1 ) 𝑎1 𝑟1 𝑏1 − 𝑟1 ) 𝑏1 𝑏1 ⁄(1 − 𝑟1 )
and
1 1⁄𝑚 0 1 (𝑚 − 1)⁄𝑚
𝑆𝑚 [(𝑎 ) ⁄(1 ) ]=[ ][ ]+[ ]
𝑚 + 𝑏𝑚 − 𝑟1 𝑎𝑚 𝑟𝑚 (𝑎𝑚 + 𝑏𝑚 )⁄(1 − 𝑟1 ) 𝑏𝑚
1
= [(𝑎 ].
𝑚 + 𝑏𝑚 )⁄(1 − 𝑟1 )
Therefore, 𝑝1 = (0, 𝑏1 ⁄(1 − 𝑟1 )) and 𝑝𝑚 = (1, (𝑎𝑚 + 𝑏𝑚 )⁄(1 − 𝑟1 )) are fixed points of 𝑆1 and
𝑆𝑚 , respectively.
The transformation 𝑆𝑖 maps the straight line segment from 𝑝1 to 𝑝𝑚 to the straight line
segment from 𝑆𝑖 (𝑝1 ) to 𝑆𝑖 (𝑝𝑚 ), which we denote [𝑆𝑖 (𝑝1 ), 𝑆𝑖 (𝑝𝑚 )]. We let 𝐸0 = [𝑝1 , 𝑝𝑚 ] and
choose the matrix entries so that
Fractal Geometry - Notes
𝑆𝑖 (𝑝𝑚 ) = 𝑆𝑖+1 (𝑝1 )
to ensure that the line segments [𝑆𝑖 (𝑝1), 𝑆𝑖 (𝑝𝑚 )] join up to form a polygonal curve 𝐸1 . To avoid
trivial cases, we assume that the points 𝑆1 (𝑝1 ) = 𝑝1 , 𝑆2 (𝑝1), … , 𝑆𝑚 (𝑝1 ), 𝑝𝑚 are not all collinear.
The attractor 𝐹 of the iterated function system {𝑆1 , 𝑆2 , … , 𝑆𝑚 } may be constructed by
repeatedly replacing line segments by affine images of the ‘generator’ set 𝐸1 . The condition
𝑆𝑖 (𝑝𝑚 ) = 𝑆𝑖+1 (𝑝1 ) ensures that the segments join up with the result that 𝐹 is the graph of
some continuous function 𝑓: [0,1] → ℝ.
It should be noted that the above conditions do not necessarily imply that the transformations
𝑆𝑖 are contractions with respect to Euclidean distance. However, it is possible to redefine
distance in the (𝑡, 𝑥) plane in such a way that the 𝑆𝑖 become contractions, and, in this context,
the IFS theory guarantees a unique attractor.
Let 𝐹 = graph 𝑓 be the self-affine curve described above, then
log(𝑟1 + ⋯ + 𝑟𝑚 )
dimB (graph 𝑓) = 1 + .
log 𝑚
Self-affine functions are useful for fractal interpolation. Suppose we wish to find a fractal
curve of a given dimension passing through the points with coordinates (𝑖 ⁄𝑚 , 𝑥𝑖 ) for 𝑖 =
0,1, … , 𝑚. By choosing transformations in such a way that 𝑆𝑖 maps the segment [𝑝1 , 𝑝𝑚 ] onto
the segment [((𝑖 − 1)⁄𝑚 , 𝑥𝑖−1 ), (𝑖 ⁄𝑚 , 𝑥𝑖 )] for each 𝑖, the construction described above gives a
self-affine function with graph passing through the given points. By adjusting the values of the
matrix entries, we can ensure that the curve has the required box dimension.
Self-affine curves can be generalised so that the 𝑆𝑖 do not have the same contraction in the 𝑡
direction. This leads to fractal interpolation between points at unequally spaced intervals of 𝑡.
The calculation of Example 11.4 may be extended to give the box dimension of such curves.
Fractal Geometry - Notes
Dynamical Systems
Discrete dynamical system Let 𝐷 ⊆ ℝ𝑛 , and let 𝑓: 𝐷 → 𝐷 be a continuous mapping. A
discrete dynamical system is an iterative scheme {𝑓 𝑘 } where 𝑓 𝑘 denotes the 𝑘th iterate of 𝑓.
Orbit The orbit of a discrete dynamical system about a point 𝑥 is the sequence of iterates
𝑘=1 .
{𝑓 𝑘 (𝑥)}∞
Dense orbit Suppose that 𝐹 is an attractor, or repeller of a function 𝑓 and that 𝑦 ∈ 𝐹. Then
the orbit {𝑓 𝑘 (𝑦)} is said to be dense in 𝐹 if, given any point 𝑥 ∈ 𝐹 and any 𝜀 > 0, there exists a
𝑘 ∈ ℕ such that the distance between 𝑓 𝑘 (𝑦) and 𝑥 is at most 𝜀.
The behaviour 𝑓 𝑘 (𝑥)as 𝑘 → ∞ may be one of the following:
(i) convergence to a fixed point 𝑤 ∈ 𝐷, i.e. 𝑓(𝑤) = 𝑤;
(ii) convergence to a period-𝑝 orbit, i.e. {𝑤, 𝑓(𝑤), … , 𝑓 𝑝−1 (𝑤)} where 𝑝 is the least positive
integer with 𝑓 𝑝 (𝑤) = 𝑤 in the sense that |𝑓 𝑘 (𝑥) − 𝑓 𝑘 (𝑤)| → 0 as 𝑘 → ∞;
(iii) random sequence of iterates that remain bounded, these often remain close to a certain
set which may be a fractal;
(iv) divergence to infinity.
Attractor Let 𝐷 ⊆ ℝ𝑛 and let 𝑓: 𝐷 → 𝐷 be a continuous mapping. A subset 𝐹 of 𝐷 is called
an attractor for 𝑓 if 𝐹 is a closed set that is invariant under 𝑓 (i.e. with 𝑓(𝐹) = 𝐹) such that the
distance from 𝑓 𝑘 (𝑥) to 𝐹 converges to zero as 𝑘 tends to infinity for all 𝑥 in an open set 𝑉
containing 𝐹, and 𝐹 has no proper subset satisfying these conditions. The largest such open
set 𝑉 is called the basin of attraction of 𝐹.
Repeller Let 𝐷 ⊆ ℝ𝑛 and let 𝑓: 𝐷 → 𝐷 be a continuous mapping. A subset 𝐹 of 𝐷 is called a
repeller for 𝑓 if 𝐹 is a closed set that is invariant under 𝑓 (i.e. with 𝑓(𝐹) = 𝐹) such that there
exists some 𝑘 for which 𝑓 𝑘 (𝑥) ∉ 𝑉 for all 𝑥 ∈ 𝑉 ∖ 𝐹 in an open set 𝑉 containing 𝐹, and 𝐹 has no
proper subset satisfying these conditions.
Roughly speaking, an attractor is a set to which all nearby orbits converge, whereas a repeller
is a set from which all nearby orbits diverge.
Since 𝑓(𝐷) ⊂ 𝐷, repeatedly applying 𝑓 gives
𝑓 𝑘 (𝐷) ⊂ 𝑓 𝑘−1 (𝐷) ⊂ ⋯ ⊂ 𝑓(𝐷) ⊂ 𝐷
so that
𝑘
⋂ 𝑓 𝑖 (𝐷) = 𝑓 𝑘 (𝐷).
𝑖=1
𝑘
Thus, the set 𝐹 = ⋂∞ 𝑖=1 𝑓 (𝐷) is invariant under 𝑓. Since 𝑓 (𝑥) ∈ ⋂𝑖=1 𝑓 (𝐷) for all 𝑥 ∈ 𝐷, the
𝑖 𝑘 𝑖
𝑊 ≡ ⋃ 𝑓 𝑘 (𝑈)
𝑘=𝑗
is the whole of ℂ, except possibly for a single point. Any such exceptional point is not in 𝐽(𝑓),
and is independent of 𝑤 and 𝑈.
Corollary 14.8 Let 𝑓(𝑧) be a polynomial.
(a) The following holds for all 𝑧 ∈ ℂ with at most one exception: if 𝑈 is an open set
intersecting 𝐽(𝑓), then 𝑓 −𝑘 (𝑧) intersects 𝑈 for infinitely many values of 𝑘.
(b) If 𝑧 ∈ 𝐽(𝑓), then 𝐽(𝑓) is the closure of ⋃∞
𝑘=1 𝑓
−𝑘 (𝑧)
.
Proposition 14.9 Let 𝑓(𝑧) be a polynomial, then 𝐽(𝑓) is a perfect set (i.e. closed and with no
isolated points) and is therefore uncountable.
Theorem 14.10 Let 𝑓(𝑧) be a polynomial, then 𝐽(𝑓) is the closure of the repelling periodic
points of 𝑓.
Basin of attraction If 𝛼 is an attracting fixed point of an analytic function 𝑓, then the basin
of attraction 𝐴(𝛼) of 𝛼 under 𝑓 is the set
{𝑧: 𝑓 𝑘 (𝑧) → 𝛼 as 𝑘 → ∞}.
Lemma 14.11 Let 𝑤 be an attractive fixed point of the polynomial 𝑓. Then
Fractal Geometry - Notes
𝜕𝐴(𝑤) = 𝐽(𝑓).
The same is true if 𝑤 = ∞.
In summary then, the Julia set 𝐽(𝑓) of a polynomial 𝑓 is the boundary of the set of points 𝑧 ∈ ℂ
such that 𝑓 𝑘 (𝑧) → ∞. It is an uncountable non-empty compact set containing no isolated
points and is invariant under 𝑓 and 𝑓 −1 . Furthermore, 𝐽(𝑓) = 𝐽(𝑓 𝑃 ) for every positive integer
𝑝. If 𝑧 ∈ 𝐽(𝑓), then 𝐽(𝑓) is the closure of ⋃∞
𝑘=1 𝑓
−𝑘 (𝑧)
. The Julia set is the boundary of the basin
of attraction of each attractive fixed point of 𝑓, including ∞, and is the closure of the repelling
periodic points of 𝑓.