Articulo de referencia

Curva de blancmange

Un gráfico de la curva del blancmange. En matemáticas , la curva blancmange es una curva fractal autoafín que se puede construir mediante subdivisión por puntos medios. También ...

Un gráfico de la curva del blancmange.

En matemáticas , la curva blancmange es una curva fractal autoafín que se puede construir mediante subdivisión por puntos medios. También se la conoce como curva Takagi , en honor a Teiji Takagi, quien la describió en 1901, o como curva Takagi-Landsberg , una generalización de la curva que lleva el nombre de Takagi y Georg Landsberg . El nombre blancmange proviene de su parecido con un pudín de blancmange . Es un caso particular de la curva de De Rham, más general .

Definición

La función blancmange se define en el intervalo unitario por

blanco(incógnita)=norte=0s(2norteincógnita)2norte,{\displaystyle \operatorname {blanc} (x)=\sum _{n=0}^{\infty }{s(2^{n}x) \over 2^{n}},}

dóndes(incógnita){\displaystyle s(x)}es la onda triangular , definida pors(incógnita)=minnorteZ|incógnitanorte|{\displaystyle s(x)=\min _{n\in {\mathbf {Z} }}|xn|}, eso es,s(incógnita){\displaystyle s(x)}es la distancia desde x al entero más cercano .

La curva de Takagi-Landsberg es una ligera generalización, dada por

Tw(incógnita)=norte=0wnortes(2norteincógnita){\displaystyle T_{w}(x)=\sum _{n=0}^{\infty }w^{n}s(2^{n}x)}

para un parámetrow{\displaystyle w}; por lo tanto, la curva del blancmange es el casow=1/2{\displaystyle w=1/2}. El valorH=registro2w{\displaystyle H=-\log _{2}w}se conoce como el parámetro de Hurst .

La función se puede extender a toda la recta real: al aplicar la definición dada anteriormente, se observa que la función se repite en cada intervalo unitario.

Definición de ecuación funcional

La versión periódica de la curva de Takagi también puede definirse como la única solución acotada.T=Tw:RR{\displaystyle T=T_{w}:\mathbb {R} \to \mathbb {R} }a la ecuación funcional

T(incógnita)=s(incógnita)+wT(2incógnita).{\displaystyle T(x)=s(x)+wT(2x).}

De hecho, la función del blancmangeTw{\displaystyle T_{w}}ciertamente está acotado y resuelve la ecuación funcional, ya que

Tw(incógnita):=norte=0wnortes(2norteincógnita)=s(incógnita)+norte=1wnortes(2norteincógnita){\displaystyle T_{w}(x):=\sum _{n=0}^{\infty }w^{n}s(2^{n}x)=s(x)+\sum _{n=1}^{\infty }w^{n}s(2^{n}x)}=s(incógnita)+wnorte=0wnortes(2norte+1incógnita)=s(incógnita)+wTw(2incógnita).{\displaystyle =s(x)+w\sum _{n=0}^{\infty }w^{n}s(2^{n+1}x)=s(x)+wT_{w}(2x).}

Por el contrario, siT:RR{\displaystyle T:\mathbb {R} \to \mathbb {R} }es una solución acotada de la ecuación funcional, que itera la igualdad que se tiene para cualquier N

T(incógnita)=norte=0nortewnortes(2norteincógnita)+wnorte+1T(2norte+1incógnita)=norte=0nortewnortes(2norteincógnita)+o(1), para norte,{\displaystyle T(x)=\sum _{n=0}^{N}w^{n}s(2^{n}x)+w^{N+1}T(2^{N+1}x)=\sum _{n=0}^{N}w^{n}s(2^{n}x)+o(1),{\text{ para }}N\to \infty ,}

De dóndeT=Tw{\displaystyle T=T_{w}}. Por cierto, las ecuaciones funcionales anteriores poseen infinitas soluciones continuas y no acotadas, por ejemplo Tw(incógnita)+do|incógnita|registro2w.{\displaystyle T_{w}(x)+c|x|^{-\log _{2}w}.}

Construcción gráfica

La curva de blancmange se puede construir visualmente a partir de funciones de onda triangulares si la suma infinita se aproxima mediante sumas finitas de los primeros términos. En las ilustraciones siguientes, se añaden a la curva funciones triangulares progresivamente más finas (mostradas en rojo) en cada etapa.

Propiedades

Convergencia y continuidad

La suma infinita que defineTw(incógnita){\displaystyle T_{w}(x)}converge absolutamente para todosincógnita.{\displaystyle x.}Desde0s(incógnita)1/2{\displaystyle 0\leq s(x)\leq 1/2}a pesar deincógnitaR,{\displaystyle x\in \mathbb {R} ,}

norte=0|wnortes(2norteincógnita)|12norte=0|w|norte=1211|w|{\displaystyle \sum _{n=0}^{\infty }|w^{n}s(2^{n}x)|\leq {\frac {1}{2}}\sum _{n=0}^{\infty }|w|^{n}={\frac {1}{2}}\cdot {\frac {1}{1-|w|}}}

si|w|<1.{\displaystyle |w|<1.}La curva de Takagi del parámetrow{\displaystyle w}se define en el intervalo unitario (oR{\displaystyle \mathbb {R} }) si|w|<1{\displaystyle |w|<1}. La función de Takagi de parámetrow{\displaystyle w}es continuo . Las funcionesTw,norte{\displaystyle T_{w,n}}definido por las sumas parciales

Tw,norte(incógnita)=k=0nortewks(2kincógnita){\displaystyle T_{w,n}(x)=\sum _{k=0}^{n}w^{k}s(2^{k}x)}

son continuas y convergen uniformemente haciaTw:{\displaystyle T_{w}:}

|Tw(incógnita)Tw,norte(incógnita)|=|k=norte+1wks(2kincógnita)|=|wnorte+1k=0wks(2k+norte+1incógnita)||w|norte+1211|w|{\displaystyle {\begin{aligned}\left|T_{w}(x)-T_{w,n}(x)\right|&=\left|\sum _{k=n+1}^{\infty }w^{k}s(2^{k}x)\right|\\&=\left|w^{n+1}\sum _{k=0}^{\infty }w^{k}s(2^{k+n+1}x)\right|\\&\leq {\frac {|w|^{n+1}}{2}}\cdot {\frac {1}{1-|w|}}\end{aligned}}}

para todo x cuando|w|<1.{\displaystyle |w|<1.}Este límite disminuye a medida quenorte.{\displaystyle n\to \infty .}Por el teorema del límite uniforme ,Tw{\displaystyle T_{w}}es continua si | w |  <  1.

Subaditividad

Dado que el valor absoluto es una función subaditiva, también lo es la función.s(incógnita)=minnorteZ|incógnitanorte|{\displaystyle s(x)=\min _{n\in {\mathbf {Z} }}|x-n|}y sus dilatacioness(2kincógnita){\displaystyle s(2^{k}x)}Dado que las combinaciones lineales positivas y los límites puntuales de funciones subaditivas son subaditivos, la función de Takagi es subaditiva para cualquier valor del parámetro.w{\displaystyle w}.

El caso especial de la parábola

Paraw=1/4{\displaystyle w=1/4}, se obtiene la parábola : la construcción de la parábola por subdivisión del punto medio fue descrita por Arquímedes .

Diferenciabilidad

Para valores del parámetro0<w<1/2,{\displaystyle 0<w<1/2,}la función TakagiTw{\displaystyle T_{w}} es diferenciable en el sentido clásico en cualquierincógnitaR{\displaystyle x\in \mathbb {R} }que no es un racional diádico . Por derivación bajo el signo de la serie, para cualquier racional no diádicoincógnitaR,{\displaystyle x\in \mathbb {R} ,}uno encuentra

Tw(incógnita)=norte=0(2w)norte(2bnorte1){\displaystyle T_{w}^{\prime }(x)=\sum _{n=0}^{\infty }(2w)^{n}\,(2b_{n}-1)}

dónde(bnorte)nortenorte{0,1}norte{\displaystyle (b_{n})_{n\in \mathbb {N} }\in \{0,1\}^{\mathbb {N} }}es la secuencia de dígitos binarios en la expansión en base 2 deincógnita{\displaystyle x}:

incógnita=norte=kbnorte2norte1.{\displaystyle x=\sum _{n=-k}^{\infty }b_{n}2^{-n-1}\;.}

De forma equivalente, los bits en la expansión binaria pueden entenderse como una secuencia de ondas cuadradas , las ondículas de Haar , escaladas a ancho2norte.{\displaystyle 2^{-n}.}Esto se deduce del hecho de que la derivada de la onda triangular es simplemente la onda cuadrada:

ddincógnitas(incógnita)=sgn(1/2(incógnitamod1)){\displaystyle {\frac {d}{dx}}s(x)=\operatorname {sgn}(1/2-(x\!\!\!\mod 1))}

y entonces

Tw(incógnita)=norte=0(2w)nortesgn(1/2(2norteincógnitamod1)){\displaystyle T_{w}^{\prime }(x)=\sum _{n=0}^{\infty }(2w)^{n}\operatorname {sgn}(1/2-(2^{n}x\!\!\!\mod 1))}

Para el parámetro0<w<1/2,{\displaystyle 0<w<1/2,}la función Tw{\displaystyle T_{w}}es Lipschitz de constante1/(12w).{\displaystyle 1/(1-2w).}En particular por el valor especialw=1/4{\displaystyle w=1/4}uno encuentra, para cualquier racional no diádicoincógnita[0,1]{\displaystyle x\in [0,1]}T1/4(incógnita)=24incógnita{\displaystyle T_{1/4}'(x)=2-4x}, de acuerdo con lo mencionado T1/4(incógnita)=2incógnita(1incógnita).{\displaystyle T_{1/4}(x)=2x(1-x).}

Paraw=1/2{\displaystyle w=1/2}la función blancmangeTw{\displaystyle T_{w}}es de variación acotada en ningún conjunto abierto no vacío; ni siquiera es localmente Lipschitz, pero es cuasi-Lipschitz, de hecho, admite la función ω(t):=t(|registro2t|+1/2){\displaystyle \omega (t):=t(|\log _{2}t|+1/2)}como módulo de continuidad .

expansión de la serie de Fourier

La función de Takagi-Landsberg admite un desarrollo en serie de Fourier absolutamente convergente:

Tw(incógnita)=metro=0ametroporque(2πmetroincógnita){\displaystyle T_{w}(x)=\sum _{m=0}^{\infty }a_{m}\cos(2\pi mx)}

cona0=1/4(1w){\displaystyle a_{0}=1/4(1-w)}y, parametro1{\displaystyle m\geq 1}

ametro:=2π2metro2(4w)ν(metro),{\displaystyle a_{m}:=-{\frac {2}{\pi ^{2}m^{2}}}(4w)^{\nu (m)},}

dónde2ν(metro){\displaystyle 2^{\nu (m)}}es la potencia máxima de2{\displaystyle 2}que dividemetro{\displaystyle m}De hecho, la onda triangular anteriors(incógnita){\displaystyle s(x)}tiene una expansión en serie de Fourier absolutamente convergente

s(incógnita)=142π2k=01(2k+1)2porque(2π(2k+1)incógnita).{\displaystyle s(x)={\frac {1}{4}}-{\frac {2}{\pi ^{2}}}\sum _{k=0}^{\infty }{\frac {1}{(2k+1)^{2}}}\cos {\big (}2\pi (2k+1)x{\big )}.}

Mediante convergencia absoluta, se puede reordenar la serie doble correspondiente paraTw(incógnita){\displaystyle T_{w}(x)}:

Tw(incógnita):=norte=0wnortes(2norteincógnita)=14norte=0wnorte2π2norte=0k=0wnorte(2k+1)2porque(2π2norte(2k+1)incógnita):{\displaystyle T_{w}(x):=\sum _{n=0}^{\infty }w^{n}s(2^{n}x)={\frac {1}{4}}\sum _{n=0}^{\infty }w^{n}-{\frac {2}{\pi ^{2}}}\sum _{n=0}^{\infty }\sum _{k=0}^{\infty }{\frac {w^{n}}{(2k+1)^{2}}}\cos {\big (}2\pi 2^{n}(2k+1)x{\big )}\,:}

puestametro=2norte(2k+1){\displaystyle m=2^{n}(2k+1)}produce la serie de Fourier anterior paraTw(incógnita).{\displaystyle T_{w}(x).}

Autosimilitud

La definición recursiva permite dar el monoide de auto-simetrías de la curva. Este monoide viene dado por dos generadores, g y r , que actúan sobre la curva (restringida al intervalo unitario) como

[gramoTw](incógnita)=Tw(gramoincógnita)=Tw(incógnita2)=incógnita2+wTw(incógnita){\displaystyle [g\cdot T_{w}](x)=T_{w}\left(g\cdot x\right)=T_{w}\left({\frac {x}{2}}\right)={\frac {x}{2}}+wT_{w}(x)}

y

[rTw](incógnita)=Tw(rincógnita)=Tw(1incógnita)=Tw(incógnita).{\displaystyle [r\cdot T_{w}](x)=T_{w}(r\cdot x)=T_{w}(1-x)=T_{w}(x).}

Un elemento general del monoide tiene entonces la formaγ=gramoa1rgramoa2rrgramoanorte{\displaystyle \gamma =g^{a_{1}}rg^{a_{2}}r\cdots rg^{a_{n}}}para algunos números enterosa1,a2,,anorte{\displaystyle a_{1},a_{2},\cdots ,a_{n}}Esto actúa sobre la curva como una función lineal :γTw=a+bincógnita+doTw{\displaystyle \gamma \cdot T_{w}=a+bx+cT_{w}}para algunas constantes a , b y c . Debido a que la acción es lineal, se puede describir en términos de un espacio vectorial , con la base del espacio vectorial :

1mi1=[100]{\displaystyle 1\mapsto e_{1}={\begin{bmatrix}1\\0\\0\end{bmatrix}}}
incógnitami2=[010]{\displaystyle x\mapsto e_{2}={\begin{bmatrix}0\\1\\0\end{bmatrix}}}
Twmi3=[001]{\displaystyle T_{w}\mapsto e_{3}={\begin{bmatrix}0\\0\\1\end{bmatrix}}}

En esta representación , la acción de g y r viene dada por

gramo=[1000121200w]{\displaystyle g={\begin{bmatrix}1&0&0\\0&{\frac {1}{2}}&{\frac {1}{2}}\\0&0&w\end{bmatrix}}}

y

r=[110010001]{\displaystyle r={\begin{bmatrix}1&1&0\\0&-1&0\\0&0&1\end{bmatrix}}}

Es decir, la acción de un elemento general.γ{\displaystyle \gamma }mapea la curva de blancmange en el intervalo unitario [0,1] a un subintervalo[metro/2pag,norte/2pag]{\displaystyle [m/2^{p},n/2^{p}]}para algunos enteros m , n , p . La aplicación viene dada exactamente por[γTw](incógnita)=a+bincógnita+doTw(incógnita){\displaystyle [\gamma \cdot T_{w}](x)=a+bx+cT_{w}(x)}donde los valores de a , b y c se pueden obtener directamente multiplicando las matrices anteriores. Es decir:

γ=[1metro2paga0nortemetro2pagb00do]{\displaystyle \gamma ={\begin{bmatrix}1&{\frac {m}{2^{p}}}&a\\0&{\frac {n-m}{2^{p}}}&b\\0&0&c\end{bmatrix}}}

Tenga en cuenta quepag=a1+a2++anorte{\displaystyle p=a_{1}+a_{2}+\cdots +a_{n}}es inmediato.

El monoide generado por g y r se denomina a veces monoide diádico ; es un submonoide del grupo modular . Al hablar del grupo modular, la notación más común para g y r es T y S , pero esta notación entra en conflicto con los símbolos utilizados aquí.

La representación tridimensional anterior es solo una de las muchas representaciones posibles; muestra que la curva de blancmange es una posible realización de la acción. Es decir, existen representaciones para cualquier dimensión, no solo para 3; algunas de ellas dan lugar a las curvas de De Rham .

Integración de la curva de Blancmange

Dado que la integral deblanco(incógnita){\displaystyle \operatorname {blanc} (x)}de 0 a 1 es 1/2, la identidadblanco(incógnita)=blanco(2incógnita)/2+s(incógnita){\displaystyle \operatorname {blanc} (x)=\operatorname {blanc} (2x)/2+s(x)}permite calcular la integral sobre cualquier intervalo mediante la siguiente relación. El cálculo es recursivo con un tiempo de cálculo del orden del logaritmo de la precisión requerida. Definiendo

I(incógnita)=0incógnitablanco(y)dy{\displaystyle I(x)=\int _{0}^{x}\operatorname {blanc} (y)\,dy}

uno tiene eso

I(incógnita)={I(2incógnita)/4+incógnita2/2si 0incógnita1/21/2I(1incógnita)si 1/2incógnita1norte/2+I(incógnitanorte)si norteincógnita(norte+1){\displaystyle I(x)={\begin{cases}I(2x)/4+x^{2}/2&{\text{if }}0\leq x\leq 1/2\\1/2-I(1-x)&{\text{if }}1/2\leq x\leq 1\\n/2+I(x-n)&{\text{if }}n\leq x\leq (n+1)\\\end{cases}}}

La integral definida viene dada por:

abblanco(y)dy=I(b)I(a).{\displaystyle \int _{a}^{b}\operatorname {blanc} (y)\,dy=I(b)-I(a).}

Se puede obtener una expresión más general definiendo

S(incógnita)=0incógnitas(y)dy={incógnita2/2,0incógnita12incógnita2/2+incógnita1/4,12incógnita1norte/4+S(incógnitanorte),(norteincógnitanorte+1){\displaystyle S(x)=\int _{0}^{x}s(y)dy={\begin{cases}x^{2}/2,&0\leq x\leq {\frac {1}{2}}\\-x^{2}/2+x-1/4,&{\frac {1}{2}}\leq x\leq 1\\n/4+S(x-n),&(n\leq x\leq n+1)\end{cases}}}

lo cual, combinado con la representación de la serie, da

Iw(incógnita)=0incógnitaTw(y)dy=norte=0(w/2)norteS(2norteincógnita){\displaystyle I_{w}(x)=\int _{0}^{x}T_{w}(y)dy=\sum _{n=0}^{\infty }(w/2)^{n}S(2^{n}x)}

Tenga en cuenta que

Iw(1)=14(1w){\displaystyle I_{w}(1)={\frac {1}{4(1-w)}}}

Esta integral también es autosimilar en el intervalo unitario, bajo una acción del monoide diádico descrito en la sección Autosimilitud . Aquí, la representación es de 4 dimensiones, teniendo la base{mi1,mi2,mi3,mi4}={1,incógnita,incógnita2,Iw(incógnita)}{\displaystyle \{e_{1},e_{2},e_{3},e_{4}\}=\{1,x,x^{2},I_{w}(x)\}}La acción de g sobre el intervalo unitario es el diagrama conmutativo

[gramoIw](incógnita)=Iw(gramoincógnita)=Iw(incógnita2)=incógnita28+w2Iw(incógnita).{\displaystyle [g\cdot I_{w}](x)=I_{w}\left(g\cdot x\right)=I_{w}\left({\frac {x}{2}}\right)={\frac {x^{2}}{8}}+{\frac {w}{2}}I_{w}(x).}

A partir de esto, se pueden leer inmediatamente los generadores de la representación tetradimensional:

gramo=[100001200001418000w2]{\displaystyle g={\begin{bmatrix}1&0&0&0\\0&{\frac {1}{2}}&0&0\\0&0&{\frac {1}{4}}&{\frac {1}{8}}\\0&0&0&{\frac {w}{2}}\end{bmatrix}}}

y

r=[11114(1w)012000100001]{\displaystyle r={\begin{bmatrix}1&1&1&{\frac {1}{4(1-w)}}\\0&-1&-2&0\\0&0&1&0\\0&0&0&-1\end{bmatrix}}}

Las integrales repetidas se transforman bajo una representación de 5, 6, ... dimensiones.

Relación con los complejos simpliciales

Dejar

norte=(nortett)+(nortet1t1)++(nortejj),nortet>nortet1>>nortejj1.{\displaystyle N={\binom {n_{t}}{t}}+{\binom {n_{t-1}}{t-1}}+\ldots +{\binom {n_{j}}{j}},\quad n_{t}>n_{t-1}>\ldots >n_{j}\geq j\geq 1.}

Definir la función Kruskal-Katona

κt(norte)=(nortett+1)+(nortet1t)++(nortejj+1).{\displaystyle \kappa _{t}(N)={n_{t} \choose t+1}+{n_{t-1} \choose t}+\dots +{n_{j} \choose j+1}.}

El teorema de Kruskal-Katona establece que este es el número mínimo de ( t 1)-símplexes que son caras de un conjunto de N t- símplexes.  

A medida que t y N se aproximan al infinito, κt(norte)norte{\displaystyle \kappa _{t}(N)-N}(debidamente normalizada) se aproxima a la curva del blancmange.

Véase también

Referencias

  • Weisstein, Eric W. "Función de Blancmange" . MathWorld .
  • Takagi, Teiji (1901), "Un ejemplo sencillo de función continua sin derivada", Proc. Phys.-Math. Soc. Jpn. , 1 : 176– 177, doi : 10.11429/subutsuhokoku1901.1.F176
  • Benoit Mandelbrot , "Paisajes fractales sin pliegues y con ríos", que aparece en The Science of Fractal Images , editado por Heinz-Otto Peitgen y Dietmar Saupe; Springer-Verlag (1988), págs. 243-260 .
  • Linas Vepstas, Simetrías de mapas de duplicación de período , (2004)
  • Donald Knuth , El arte de la programación informática , volumen 4a. Algoritmos combinatorios, parte 1. ISBN 0-201-03804-8Véanse las páginas 372-375.

Lecturas adicionales

  • Allaart, Pieter C.; Kawamura, Kiko (11 de octubre de 2011), La función Takagi: una revisión , arXiv : 1110.1691 , Bibcode : 2011arXiv1110.1691A
  • Lagarias, Jeffrey C. (17 de diciembre de 2011), La función de Takagi y sus propiedades , arXiv : 1112.4205 , Bibcode : 2011arXiv1112.4205L
  • Explorador Takagi
  • (Algunas propiedades de la función de Takagi)