Midpoint displacement

In mathematics, the blancmange curve is a self-affine curve constructible by midpoint subdivision. It is also known as the Takagi curve, after Teiji Takagi who described it in 1901, or as the Takagi–Landsberg curve, a generalization of the curve named after Takagi and Georg Landsberg. The name ''blancmange'' comes from its resemblance to a Blancmange pudding. It is a special case of the more general de Rham curve; see also fractal curve.

Definition

The blancmange function is defined on the unit interval by

: $(x) = \sum_^\infty ,$

where $s(x)$ is the triangle wave, defined by $s(x)=\min_, x-n,$,
that is, $s(x)$ is the distance from ''x'' to the nearest integer.

The Takagi–Landsberg curve is a slight generalization, given by

: $T_w(x) = \sum_^\infty w^n s(2^x)$

for a parameter $w$; thus the blancmange curve is the case $w=1/2$. The value $H=-\log_2 w$ is known as the Hurst parameter.

The function can be extended to all of the real line: applying the definition given above shows that the function repeats on each unit interval.

The function could also be defined by the series in the section Fourier series expansion.

Functional equation definition

The periodic version of the Takagi curve can also be defined as the ''unique bounded solution $T=T_w:\R\to\R$ to the functional equation''

: $T(x) = s(x) + w T(2x)$.

Indeed, the blancmange function $T_w$ is certainly bounded, and solves the functional equation, since
:$T_w(x) := \sum_^\infty w^n s(2^x)= s(x) + \sum_^\infty w^n s(2^x)$$= s(x) + w\sum_^\infty w^n s(2^x)= s(x) + wT_w(2x)$.
Conversely, if $T:\R\to\R$ is a bounded solution of the functional equation, iterating the equality one has for any ''N''
: $T(x) =\sum_^N w^n s(2^x) + w^ T(2^x) =\sum_^N w^n s(2^x) + o(1)$, for $N\to\infty,$

whence $T=T_w$. Incidentally, the above functional equations possesses infinitely many continuous, non-bounded solutions, e.g.
$T_w(x)+ c , x, ^.$

Graphical construction

The blancmange curve can be visually built up out of triangle wave functions if the infinite sum is approximated by finite sums of the first few terms. In the illustrations below, progressively finer triangle functions (shown in red) are added to the curve at each stage.

Properties

Convergence and continuity

The infinite sum defining $T_w(x)$ converges absolutely for all $x$: since $0\le s(x) \le 1/2$ for all $x\in \mathbb$, we have:

: $\sum_^\infty , w^n s(2^n x), \le \frac \sum_^\infty , w, ^n = \frac \cdot \frac$ if $, w, <1$.

Therefore, the Takagi curve of parameter $w$ is defined on the unit interval (or $\mathbb$) if $, w, <1$.

The Takagi function of parameter $w$ is continuous. Indeed, the functions $T_$ defined by the partial sums $T_(x) = \sum_^n w^k s(2^k x)$ are continuous and converge uniformly toward $T_w$, since:

: $\left, T_w(x) - T_(x)\ = \left, \sum_^\infty w^k s(2^k x)\ = \left, w^ \sum_^\infty w^k s(2^ x)\ \le \frac \cdot \frac$ for all x when $, w, < 1$.

This value can be made as small as we want by selecting a big enough value of ''n''. Therefore, by the uniform limit theorem, $T_w$ is continuous if , ''w'', <1.

File:Blancmange k1.5.png, parameter w=2/3
File:Blancmange k2.png, parameter w=1/2
File:Blancmange k3.png, parameter w=1/3
File:Blancmange k4.png, parameter w=1/4
File:Blancmange k8.png, parameter w=1/8

Subadditivity

Since the absolute value is a subadditive function so is the function $s(x)=\min_, x-n,$, and its dilations $s(2^kx)$; since positive linear combinations and point-wise limits of subadditive functions are subadditive, the Takagi function is subadditive for any value of the parameter $w$.

The special case of the parabola

For $w=1/4$, one obtains the parabola: the construction of the parabola by midpoint subdivision was described by Archimedes.

Differentiability

For values of the parameter $0< w < 1/2$ the Takagi function $T_w$ is differentiable in classical sense at any $x\in\R$ which is not a dyadic rational. Precisely,
by derivation under the sign of series, for any non dyadic rational $x\in\R$ one finds

:$T_w'(x) = \sum_^\infty (2w)^n \,(-1)^$

where $(x_n)_\in\^\Z$ is the sequence of binary digits in the base 2 expansion of $x$, that is, $x=\sum_2^n x_n$. Moreover, for these values of $w$ the function $T_w$ is Lipschitz of constant $1\over 1-2w$. In particular for the special value $w=1/4$ one finds, for any non dyadic rational x\in ,1/math> $T_'(x) = 2 - 4x$, according with the mentioned $T_(x) = 2x(1 - x).$

For $w=1/2$ the blancmange function $T_w$ it is of bounded variation on no non-empty open set; it is not even locally Lipschitz, but it is quasi-Lipschitz, indeed, it admits the function $\omega(t):=t(, \log_2 t, +1/2)$ as a modulus of continuity .

Fourier series expansion

The Takagi-Landsberg function admits an absolutely convergent Fourier series expansion:
:$T_w(x) =\sum_^\infty a_m\cos(2\pi m x)$
with $a_0=1/4(1-w)$ and, for $m\ge 1$
:$a_m:=-\frac(4w)^,$
where $2^$ is the maximum power of $2$ that divides $m$.
Indeed, the above triangle wave $s(x)$ has an absolutely convergent Fourier series expansion
:$s(x)=\frac-\frac\sum_^\infty\frac\cos\big(2\pi (2k+1)x\big).$
By absolute convergence, one can reorder the corresponding double series for $T_w(x)$:

:$T_w(x):=\sum_^\infty w^n s(2^nx)= \frac\sum_^\infty w^n -\frac\sum_^\infty\sum_^\infty \frac\cos\big(2\pi 2^n(2k+1)x\big)\, :$

putting $m=2^n(2k+1)$ yields the above Fourier series for $T_w(x).$

Self similarity

The recursive definition allows the monoid of self-symmetries of the curve to be given. This monoid is given by two generators, ''g'' and ''r'', which act on the curve (restricted to the unit interval) as

: \cdot T_wx) = T_w\left(\frac\right) = \frac + w T_w(x)
and

: \cdot T_wx) = T_w(1-x) = T_w(x).

A general element of the monoid then has the form $\gamma=g^ r g^ r \cdots r g^$ for some integers $a_1, a_2, \cdots, a_n$ This acts on the curve as a linear function: $\gamma \cdot T_w = a + bx + cT_w$ for some constants ''a'', ''b'' and ''c''. Because the action is linear, it can be described in terms of a vector space, with the vector space basis:

:$1 \mapsto e_1 = \begin 1 \\ 0 \\ 0 \end$
:$x \mapsto e_2 = \begin 0 \\ 1 \\ 0 \end$
:$T_w \mapsto e_3 = \begin 0 \\ 0 \\ 1 \end$

In this representation, the action of ''g'' and ''r'' are given by

:$g=\begin 1 & 0 & 0 \\ 0 & \frac & \frac \\ 0 & 0 & w \end$
and
:$r=\begin 1 & 1 & 0 \\ 0 & -1 & 0 \\ 0 & 0 & 1 \end$

That is, the action of a general element $\gamma$ maps the blancmange curve on the unit interval ,1to a sub-interval /2^p, n/2^p/math> for some integers ''m'', ''n'', ''p''. The mapping is given exactly by gamma \cdot T_wx) = a + bx + cT_w(x) where the values of ''a'', ''b'' and ''c'' can be obtained directly by multiplying out the above matrices. That is:

:$\gamma=\begin 1 & \frac & a \\ 0 & \frac & b \\ 0 & 0 & c \end$

Note that $p=a_1+a_2+\cdots +a_n$ is immediate.

The monoid generated by ''g'' and ''r'' is sometimes called the dyadic monoid; it is a sub-monoid of the modular group. When discussing the modular group, the more common notation for ''g'' and ''r'' is ''T'' and ''S'', but that notation conflicts with the symbols used here.

The above three-dimensional representation is just one of many representations it can have; it shows that the blancmange curve is one possible realization of the action. That is, there are representations for any dimension, not just 3; some of these give the de Rham curves.

Integrating the Blancmange curve

Given that the integral of $(x)$ from 0 to 1 is 1/2, the identity $(x)= (2x)/2+s(x)$ allows the integral over any interval to be computed by the following relation. The computation is recursive with computing time on the order of log of the accuracy required. Defining

:$I(x) = \int_0^x(y)\,dy$
one has that
:$I(x) =\begin I(2x)/4+x^2/2 & \text 0 \leq x \leq 1/2 \\ 1/2-I(1-x) & \text 1/2 \le x \le 1 \\ n/2+I(x-n) & \text n \le x \le (n+1) \\ \end$
The definite integral is given by:
:$\int_a^b(y)\,dy = I(b) - I(a).$

A more general expression can be obtained by defining
:$S(x)=\int_0^x s(y)dy = \begin x^2/ 2, & 0 \le x \le \frac \\ - x^2 / 2 +x - 1/4, & \frac \le x \le 1 \\ n/4 + S(x-n), & (n \le x \le n+1) \end$
which, combined with the series representation, gives

:$I_w(x)= \int_0^x T_w(y) dy = \sum_^\infty (w/2)^n S(2^n x)$

Note that
:$I_w(1)=\frac$

This integral is also self-similar on the unit interval, under an action of the dyadic monoid described in the section '' Self similarity''. Here, the representation is 4-dimensional, having the basis $\$. Re-writing the above to make the action of ''g'' more clear: on the unit interval, one has

: \cdot I_wx) = I_w\left(\frac\right) = \frac + \fracI_w(x).

From this, one can then immediately read off the generators of the four-dimensional representation:

:$g=\begin 1 & 0 & 0 & 0\\ 0 & \frac & 0 & 0 \\ 0 & 0 & \frac & \frac \\ 0 & 0 & 0 & \frac \end$
and
:$r=\begin 1 & 1 & 1 & \frac \\ 0 & -1 & -2 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & -1 \end$

Repeated integrals transform under a 5,6,... dimensional representation.

Relation to simplicial complexes

Let
: $N=\binom+\binom+\ldots+\binom,\quad n_t > n_ > \ldots > n_j \geq j\geq 1.$

Define the Kruskal–Katona function

:$\kappa_t(N)= + + \dots + .$

The Kruskal–Katona theorem states that this is the minimum number of (''t'' − 1)-simplexes that are faces of a set of ''N'' ''t''-simplexes.

As ''t'' and ''N'' approach infinity,
$\kappa_t(N)-N$ (suitably normalized) approaches the blancmange curve.

See also

* Cantor function (also known as the Devil's staircase)
* Minkowski's question mark function
* Weierstrass function
* Dyadic transformation

References

*
*
* Benoit Mandelbrot, "Fractal Landscapes without creases and with rivers", appearing in ''The Science of Fractal Images'', ed. Heinz-Otto Peitgen, Dietmar Saupe; Springer-Verlag (1988) pp 243–260.
* Linas Vepstas,
Symmetries of Period-Doubling Maps
', (2004)
* Donald Knuth, The Art of Computer Programming, volume 4a. Combinatorial algorithms, part 1. . See pages 372–375.

Further reading

*
*

External links

Takagi Explorer

(Some properties of the Takagi function)

{{DEFAULTSORT:Blancmange Curve
De Rham curves
Theory of continuous functions