De Rham curve explained

In mathematics, a de Rham curve is a continuous fractal curve obtained as the image of the Cantor space, or, equivalently, from the base-two expansion of the real numbers in the unit interval. Many well-known fractal curves, including the Cantor function, Cesàro–Faber curve (Lévy C curve), Minkowski's question mark function, blancmange curve, and the Koch curve are all examples of de Rham curves. The general form of the curve was first described by Georges de Rham in 1957.[1]

Construction

(M,d)

(generally

R

2 with the usual euclidean distance), and a pair of contracting maps on M:

d0:M\toM

d1:M\toM.

By the Banach fixed-point theorem, these have fixed points

p0

and

p1

respectively. Let x be a real number in the interval

[0,1]

, having binary expansion

x=

infty
\sum
k=1
bk
2k

,

where each

bk

is 0 or 1. Consider the map

cx:M\toM

defined by

cx=

d
b1

\circ

d
b2

\circ\circ

d
bk

\circ,

where

\circ

denotes function composition. It can be shown that each

cx

will map the common basin of attraction of

d0

and

d1

to a single point

px

in

M

. The collection of points

px

, parameterized by a single real parameter x, is known as the de Rham curve.

Continuity condition

The construction in terms of binary digits can be understood in two distinct ways. One way is as a mapping of Cantor space to distinct points in the plane. Cantor space is the set of all infinitely-long strings of binary digits. It is a discrete space, and is disconnected. Cantor space can be mapped onto the unit real interval by treating each string as a binary expansion of a real number. In this map, the dyadic rationals have two distinct representations as strings of binary digits. For example, the real number one-half has two equivalent binary expansions:

h1=0.1000 …

and

h0=0.01111 …

This is analogous to how one has 0.999...=1.000... in decimal expansions. The two points

h0

and

h1

are distinct points in Cantor space, but both are mapped to the real number one-half. In this way, the reals of the unit interval are a continuous image of Cantor space.

The same notion of continuity is applied to the de Rham curve by asking that the fixed points be paired, so that

d0(p1)=d1(p0)

With this pairing, the binary expansions of the dyadic rationals always map to the same point, thus ensuring continuity at that point. Consider the behavior at one-half. For any point p in the plane, one has two distinct sequences:

d0\circd1\circd1\circd1\circ(p)

and

d1\circd0\circd0\circd0\circ(p)

corresponding to the two binary expansions

1/2=0.01111 …

and

1/2=0.1000 …

. Since the two maps are both contracting, the first sequence converges to

d0(p1)

and the second to

d1(p0)

. If these two are equal, then both binary expansions of 1/2 map to the same point. This argument can be repeated at any dyadic rational, thus ensuring continuity at those points. Real numbers that are not dyadic rationals have only one, unique binary representation, and from this it follows that the curve cannot be discontinuous at such points. The resulting de Rham curve

px

is a continuous function of x, at all x.

In general, the de Rham curves are not differentiable.

Properties

De Rham curves are by construction self-similar, since

p(x)=d0(p(2x))

for

x\in[0,1/2]

and

p(x)=d1(p(2x-1))

for

x\in[1/2,1].

The self-symmetries of all of the de Rham curves are given by the monoid that describes the symmetries of the infinite binary tree or Cantor space. This so-called period-doubling monoid is a subset of the modular group.

The image of the curve, i.e. the set of points

\{p(x),x\in[0,1]\}

, can be obtained by an Iterated function system using the set of contraction mappings

\{d0,d1\}

. But the result of an iterated function system with two contraction mappings is a de Rham curve if and only if the contraction mappings satisfy the continuity condition.

Detailed, worked examples of the self-similarities can be found in the articles on the Cantor function and on Minkowski's question-mark function. Precisely the same monoid of self-similarities, the dyadic monoid, apply to every de Rham curve.

Classification and examples

The following systems generate continuous curves.

Cesàro curves

Cesàro curves, also known as Cesàro–Faber curves or Lévy C curves, are De Rham curves generated by affine transformations conserving orientation, with fixed points

p0=0

and

p1=1

.

a

such that

|a|<1

and

|1-a|<1

.

The contraction mappings

d0

and

d1

are then defined as complex functions in the complex plane by:

d0(z)=az

d1(z)=a+(1-a)z.

For the value of

a=(1+i)/2

, the resulting curve is the Lévy C curve.

Koch - Peano curves

In a similar way, we can define the Koch - Peano family of curves as the set of De Rham curves generated by affine transformations reversing orientation, with fixed points

p0=0

and

p1=1

.

These mappings are expressed in the complex plane as a function of

\overline{z}

, the complex conjugate of

z

:

d0(z)=a\overline{z}

d1(z)=a+(1-a)\overline{z}.

The name of the family comes from its two most famous members. The Koch curve is obtained by setting:

a
Koch=1
2

+i

\sqrt{3
},

while the Peano curve corresponds to:

a
Peano=(1+i)
2

.

The de Rham curve for

a=(1+ib)/2

for values of

b

just less than one visually resembles the Osgood curve. These two curves are closely related, but are not the same. The Osgood curve is obtained by repeated set subtraction, and thus is a perfect set, much like the Cantor set itself. The construction of the Osgood set asks that progressively smaller triangles to be subtracted, leaving behind a "fat" set of non-zero measure; the construction is analogous to the fat Cantor set, which has a non-zero measure. By contrast, the de Rham curve is not "fat"; the construction does not offer a way to "fatten up" the "line segments" that run "in between" the dyadic rationals.

General affine maps

The Cesàro–Faber and Peano–Koch curves are both special cases of the general case of a pair of affine linear transformations on the complex plane. By fixing one endpoint of the curve at 0 and the other at 1, the general case is obtained by iterating on the two transforms

d0=\begin{pmatrix} 1&0&0\ 0&\alpha&\delta\ 0&\beta&\varepsilon \end{pmatrix}

and

d1=\begin{pmatrix} 1&0&0\ \alpha&1-\alpha&\zeta\ \beta&-\beta&η \end{pmatrix}.

Being affine transforms, these transforms act on a point

(u,v)

of the 2-D plane by acting on the vector

\begin{pmatrix} 1\\ u\\ v\end{pmatrix}.

The midpoint of the curve can be seen to be located at

(u,v)=(\alpha,\beta)

; the other four parameters may be varied to create a large variety of curves.

The blancmange curve of parameter

w

can be obtained by setting

\alpha=\beta=1/2

,

\delta=\zeta=0

and

\varepsilon=η=w

. That is:

d0=\begin{pmatrix} 1&0&0\ 0&1/2&0\ 0&1/2&w \end{pmatrix}

and

d1=\begin{pmatrix} 1&0&0\ 1/2&1/2&0\ 1/2&-1/2&w \end{pmatrix}.

Since the blancmange curve for parameter

w=1/4

is a parabola of the equation

f(x)=4x(1-x)

, this illustrates the fact that on some occasions, de Rham curves can be smooth.

Minkowski's question mark function

Minkowski's question mark function is generated by the pair of maps

d0(z)=

z
z+1

and

d1(z)=

1
2-z

.

Non-examples

Given any two functions

d0

and

d1

, one can define a mapping from Cantor space, by repeated iteration of the digits, exactly the same way as for the de Rham curves. In general, the result will not be a de Rham curve, when the terms of the continuity condition are not met. Thus, there are many sets that might be in one-to-one correspondence with Cantor space, whose points can be uniquely labelled by points in the Cantor space; however, these are not de Rham curves, when the dyadic rationals do not map to the same point.

Julia set of the Mandelbrot set

The Mandelbrot set is generated by a period-doubling iterated equation

zn+1

2+c.
=z
n
The corresponding Julia set is obtained by iterating the opposite direction. This is done by writing

zn=\pm\sqrt{zn+1-c}

, which gives two distinct roots that the forward iterate

zn+1

"came from". These two roots can be distinguished as

d0(z)=+\sqrt{z-c}

and

d1(z)=-\sqrt{z-c}.

Fixing the complex number

c

, the result is the Julia set for that value of

c

. This curve is continuous when

c

is inside the Mandelbrot set; otherwise, it is a disconnected dust of points. However, the reason for continuity is not due to the de Rham condition, as, in general, the points corresponding to the dyadic rationals are far away from one-another. In fact, this property can be used to define a notion of "polar opposites", of conjugate points in the Julia set.

Generalizations

It is easy to generalize the definition by using more than two contraction mappings. If one uses n mappings, then the n-ary decomposition of x has to be used instead of the binary expansion of real numbers. The continuity condition has to be generalized in:

di(pn-1)=di+1(p0)

, for

i=0\ldotsn-2.

This continuity condition can be understood with the following example. Suppose one is working in base-10. Then one has (famously) that 0.999...= 1.000... which is a continuity equation that must be enforced at every such gap. That is, given the decimal digits

b1,b2,,bk

with

bk\ne9

, one has

b1,b2,,bk,9,9,9,=b1,b2,,bk+1,0,0,0,

Such a generalization allows, for example, to produce the Sierpiński arrowhead curve (whose image is the Sierpiński triangle), by using the contraction mappings of an iterated function system that produces the Sierpiński triangle.

Multifractal curves

Ornstein and others describe a multifractal system, where instead of working in a fixed base, one works in a variable base.

Consider the product space of variable base-

mn

discrete spaces

\Omega=\prodn\inNAn

for

An=Z/mnZ=\{0,1,,mn-1\}

the cyclic group, for

mn\ge2

an integer. Any real number in the unit interval can be expanded in a sequence

(a1,a2,a3,)

such that each

an\inAn

. More precisely, a real number

0\lex\le1

is written as
infty
x=\sum
n=1
an
n
\prodmk
k=1
This expansion is not unique, if all

an=0

past some point

K<n

. In this case, one has that

a1,a2,,aK,0,0,=a1,a2,,aK-1,mK+1-1,mK+2-1,

Such points are analogous to the dyadic rationals in the dyadic expansion, and the continuity equations on the curve must be applied at these points.

For each

An

, one must specify two things: a set of two points
(n)
p
0
and
(n)
p
1
and a set of

mn

functions
(n)
d
j

(z)

(with

j\inAn

). The continuity condition is then as above,
(n)
d
j
(n+1)
(p
1)=d
(n)
j+1
(n+1)
(p
0)
, for

j=0,,mn-2.

Ornstein's original example used

\Omega=\left(Z/2Z\right) x \left(Z/3Z\right) x \left(Z/4Z\right) x

See also

Further reading

Notes and References

  1. Georges de Rham, Sur quelques courbes definies par des equations fonctionnelles. Univ. e Politec. Torino. Rend. Sem. Mat., 1957, 16, 101  - 113