Articulo de referencia

Función refinable

En matemáticas , en el área del análisis de ondículas , una función refinable es una función que cumple algún tipo de autosimilitud . Una función φ {\displaystyle \varphi } se d...

En matemáticas , en el área del análisis de ondículas , una función refinable es una función que cumple algún tipo de autosimilitud . Una funciónφ{\displaystyle \varphi }se denomina refinable con respecto a la máscarah{\displaystyle h}si

φ(incógnita)=2k=0norte1hkφ(2incógnitak){\displaystyle \varphi (x)=2\cdot \sum _{k=0}^{N-1}h_{k}\cdot \varphi (2\cdot xk)}

Esta condición se denomina ecuación de refinamiento , ecuación de dilatación o ecuación de dos escalas .

Utilizando la convolución (denotada por una estrella, *) de una función con una máscara discreta y el operador de dilatación.D{\displaystyle D}Se puede escribir de forma más concisa:

φ=2D1/2(hφ){\displaystyle \varphi =2\cdot D_{1/2}(h*\varphi)}

Esto significa que se obtiene la función, de nuevo, si se convoluciona la función con una máscara discreta y luego se reduce su escala. Existe una similitud con los sistemas de funciones iteradas y las curvas de De Rham .

El operadorφ2D1/2(hφ){\displaystyle \varphi \mapsto 2\cdot D_{1/2}(h*\varphi)}es lineal. Una función refinable es una autofunción de ese operador. Su valor absoluto no está definido de forma única. Es decir, siφ{\displaystyle \varphi }es una función refinable, entonces para cadado{\displaystyle c}la funcióndoφ{\displaystyle c\cdot \varphi }También es refinable.

Estas funciones desempeñan un papel fundamental en la teoría de ondículas como funciones de escala .

Propiedades

Valores en puntos enteros

Una función refinable se define solo implícitamente. También puede ser que existan varias funciones que sean refinables con respecto a la misma máscara. Siφ{\displaystyle \varphi }si tiene soporte finito y se desean los valores de la función en argumentos enteros, entonces la ecuación de dos escalas se convierte en un sistema de ecuaciones lineales simultáneas .

Dejara{\displaystyle a}sea ​​el índice mínimo yb{\displaystyle b}sea ​​el índice máximo de elementos no nulos deh{\displaystyle h}, entonces se obtiene (φ(a)φ(a+1)φ(b))=(haha+2ha+1haha+4ha+3ha+2ha+1hahbhb1hb2hb3hb4hbhb1hb2hb)(φ(a)φ(a+1)φ(b)).{\displaystyle {\begin{pmatrix}\varphi (a)\\\varphi (a+1)\\\vdots \\\varphi (b)\end{pmatrix}}={\begin{pmatrix}h_{a}&&&&&\\h_{a+2}&h_{a+1}&h_{a}&&&\\h_{a+4}&h_{a+3}&h_{a+2}&h_{a+1}&h_{a}&\\\ddots &\ddots &\ddots &\ddots &\ddots &\ddots \\&h_{b}&h_{b-1}&h_{b-2}&h_{b-3}&h_{b-4}\\&&&h_{b}&h_{b-1}&h_{b-2}\\&&&&&h_{b}\end{pmatrix}}{\begin{pmatrix}\varphi (a)\\\varphi (a+1)\\\vdots \\\varphi (b)\end{pmatrix}}.}

Utilizando el operador de discretización , llámeloQ{\displaystyle Q}aquí, y la matriz de transferencia deh{\displaystyle h}, llamadoTh{\displaystyle T_{h}}, esto se puede escribir de forma concisa como Qφ=ThQφ.{\displaystyle Q\varphi =T_{h}Q\varphi .}

Esta es de nuevo una ecuación de punto fijo . Pero esta ahora puede considerarse como un problema de autovectores y autovalores . Es decir, una función refinable con soporte finito existe solo (pero no necesariamente), siTh{\displaystyle T_{h}}tiene el valor propio  1.

Valores en puntos diádicos

A partir de los valores en los puntos enteros se pueden derivar los valores en los puntos diádicos, es decir, puntos de la formak2j{\displaystyle k\cdot 2^{-j}}, conkZ{\displaystyle k\in \mathbb {Z} }yjnorte{\displaystyle j\in \mathbb {N} }.

φ=D1/2(2(hφ)){\displaystyle \varphi =D_{1/2}(2\cdot (h*\varphi ))}
D2φ=2(hφ){\displaystyle D_{2}\varphi =2\cdot (h*\varphi)}
Q(D2φ)=Q(2(hφ))=2(hQφ){\displaystyle Q(D_{2}\varphi )=Q(2\cdot (h*\varphi ))=2\cdot (h*Q\varphi )}

La estrella denota la convolución de un filtro discreto con una función. Con este paso se pueden calcular los valores en puntos de la formak2{\displaystyle {\frac {k}{2}}}. Reemplazando iterativamenteφ{\displaystyle \varphi }porD2φ{\displaystyle D_{2}\varphi }Obtendrás los valores en todas las escalas más finas.

Q(D2j+1φ)=2(hQ(D2jφ)){\displaystyle Q(D_{2^{j+1}}\varphi )=2\cdot (h*Q(D_{2^{j}}\varphi ))}

Circunvolución

Siφ{\displaystyle \varphi }es refinable con respecto ah{\displaystyle h}, yψ{\displaystyle \psi }es refinable con respecto agramo{\displaystyle g}, entoncesφψ{\displaystyle \varphi *\psi }es refinable con respecto ahgramo{\displaystyle h*g}.

Diferenciación

Siφ{\displaystyle \varphi }es refinable con respecto ah{\displaystyle h}y el derivadoφ{\displaystyle \varphi '}entonces existeφ{\displaystyle \varphi '}es refinable con respecto a2h{\displaystyle 2\cdot h}Esto puede interpretarse como un caso especial de la propiedad de convolución, donde uno de los operandos de convolución es una derivada del impulso de Dirac .

Integración

Siφ{\displaystyle \varphi }es refinable con respecto ah{\displaystyle h}y hay un antiderivadoΦ{\displaystyle \Phi }conΦ(t)=0tφ(τ)dτ{\textstyle \Phi (t)=\int _ {0}^{t}\varphi (\tau )\,\mathrm {d} \tau }, luego el antiderivadotΦ(t)+do{\displaystyle t\mapsto \Phi (t)+c}es refinable con respecto a la máscara12h{\textstyle {\frac {1}{2}}\cdot h}donde la constantedo{\displaystyle c}debe cumplirdo(1jhj)=jhjΦ(j){\textstyle c\cdot \left(1-\sum _{j}h_{j}\right)=\sum _{j}h_{j}\cdot \Phi (-j)}.

Siφ{\displaystyle \varphi }tiene soporte acotado , entonces podemos interpretar la integración como una convolución con la función de Heaviside y aplicar la ley de convolución.

productos escalares

El cálculo de los productos escalares de dos funciones refinables y sus traslaciones se puede desglosar en las dos propiedades anteriores. SeaT{\displaystyle T}ser el operador de traducción. Contiene φ,Tkψ=φψ,Tkδ=(φψ)(k){\displaystyle \langle \varphi ,T_{k}\psi \rangle =\langle \varphi *\psi ^{*},T_{k}\delta \rangle =(\varphi *\psi ^{*})(k)} dóndeψ{\displaystyle \psi ^{*}}es el adjunto deψ{\displaystyle \psi }con respecto a la convolución , es decir,ψ{\displaystyle \psi ^{*}}es la versión conjugada compleja y volteada deψ{\displaystyle \psi }, es decir,ψ(t)=ψ(t)¯{\displaystyle \psi ^{*}(t)={\overline {\psi (-t)}}}.

Debido a la propiedad anterior,φψ{\displaystyle \varphi *\psi ^{*}}es refinable con respecto ahgramo{\displaystyle h*g^{*}}y sus valores en argumentos enteros pueden calcularse como autovectores de la matriz de transferencia. Esta idea puede generalizarse fácilmente a integrales de productos de más de dos funciones refinables. [ 1 ]

Suavidad

Una función refinable suele tener forma fractal. El diseño de funciones refinables continuas o suaves no es obvio. Antes de abordar la imposición de suavidad, es necesario medir la suavidad de las funciones refinables. Utilizando la máquina de Villemoes [ 2 ], se puede calcular la suavidad de las funciones refinables en términos de exponentes de Sobolev .

En un primer paso, la máscara de refinamientoh{\displaystyle h}está dividido en un filtrob{\displaystyle b}, que es una potencia del factor de suavidad(1,1){\displaystyle (1,1)}(esta es una máscara binomial) y un restoq{\displaystyle q}En términos generales, la máscara binomialb{\displaystyle b}hace suavidad y q{\displaystyle q}representa un componente fractal, que reduce la suavidad nuevamente. Ahora el exponente de Sobolev es aproximadamente del orden deb{\displaystyle b}menos el logaritmo del radio espectral deTqq{\displaystyle T_{q*q^{*}}}.

Generalización

El concepto de funciones refinables puede generalizarse a funciones de más de una variable, es decir, funciones deRdR{\displaystyle \mathbb {R} ^{d}\to \mathbb {R} }. La generalización más simple se refiere a productos tensoriales . Siφ{\displaystyle \varphi }yψ{\displaystyle \psi }son refinables con respecto ah{\displaystyle h}ygramo{\displaystyle g}, respectivamente, entoncesφψ{\displaystyle \varphi \otimes \psi } es refinable con respecto ahgramo{\displaystyle h\otimes g}.

El esquema se puede generalizar aún más a diferentes factores de escala con respecto a diferentes dimensiones o incluso a la mezcla de datos entre dimensiones. [ 3 ] En lugar de escalar por un factor escalar como 2, la señal, las coordenadas se transforman por una matrizMETRO{\displaystyle M}de enteros. Para que el esquema funcione, los valores absolutos de todos los valores propios deMETRO{\displaystyle M}debe ser mayor que uno. (Tal vez también sea suficiente que|detMETRO|>1{\displaystyle \left|\det M\right|>1}.)

Formalmente, la ecuación de dos escalas no cambia mucho: φ(incógnita)=|detMETRO|kZdhkφ(METROincógnitak){\displaystyle \varphi (x)=\left|\det M\right|\cdot \sum _{k\in \mathbb {Z} ^{d}}h_{k}\cdot \varphi (M\cdot x-k)}φ=|detMETRO|DMETRO1(hφ){\displaystyle \varphi =\left|\det M\right|\cdot D_{M^{-1}}(h*\varphi )}

Ejemplos

  • Si la definición se extiende a distribuciones , entonces el impulso de Dirac se puede refinar con respecto al vector unitario.δ{\displaystyle \delta }, que se conoce como delta de Kronecker .norte{\displaystyle n}La derivada -ésima de la distribución de Dirac es refinable con respecto a2norteδ{\displaystyle 2^{n}\cdot \delta }.
  • La función de Heaviside es refinable con respecto a12δ{\displaystyle {\frac {1}{2}}\cdot \delta }.
  • Las funciones de potencia truncadas con exponentenorte{\displaystyle n}son refinables con respecto a12norte+1δ{\displaystyle {\frac {1}{2^{n+1}}}\cdot \delta }.
  • La función triangular es una función refinable. [ 4 ] Las funciones B-spline con nodos integrales sucesivos son refinables, debido al teorema de convolución y a la refinabilidad de la función característica para el intervalo[0,1){\displaystyle [0,1)}(una función de vagón de carga ).
  • Todas las funciones polinómicas son refinables. Para cada máscara de refinamiento existe un polinomio que está definido de forma única salvo un factor constante. Para cada polinomio de gradonorte{\displaystyle n}Hay muchas máscaras de refinamiento que difieren entre sí por un tipo de máscara.v(1,1)norte+1{\displaystyle v*(1,-1)^{n+1}}para cualquier máscarav{\displaystyle v}y el poder de convolución(1,1)norte+1{\displaystyle (1,-1)^{n+1}}. [ 5 ]
  • Una función racionalφ{\displaystyle \varphi }es refinable si y solo si puede representarse utilizando fracciones parciales comoφ(incógnita)=iZsi(incógnitai)k{\displaystyle \varphi (x)=\sum _{i\in \mathbb {Z} }{\frac {s_{i}}{(x-i)^{k}}}}, dóndek{\displaystyle k}es un número natural positivo ys{\displaystyle s}es una sucesión real con un número finito de elementos distintos de cero (un polinomio de Laurent ) tal ques|(s2){\displaystyle s|(s\uparrow 2)}(leer:h(z)R[z,z1] h(z)s(z)=s(z2){\displaystyle \exists h(z)\in \mathbb {R} [z,z^{-1}]\ h(z)\cdot s(z)=s(z^{2})}). El polinomio de Laurent2k1h{\displaystyle 2^{k-1}\cdot h}es la máscara de refinamiento asociada. [ 6 ]

Referencias

  1. Dahmen, Wolfgang; Micchelli, Charles A. (1993). "Uso de la ecuación de refinamiento para evaluar integrales de ondículas". SIAM Journal on Numerical Analysis . 30 (2): 507– 537. doi : 10.1137/0730024 .
  2. Villemoes, Lars. "Regularidad de Sobolev de ondículas y estabilidad de bancos de filtros iterados" . Archivado del original (PostScript) el 11 de mayo de 2002.
  3. Berger, Marc A.; Wang, Yang (1992), "Ecuaciones de dilatación de dos escalas multidimensionales (capítulo IV)", en Chui, Charles K. (ed.), Wavelets: A Tutorial in Theory and Applications , Wavelet Analysis and its Applications, vol. 2, Academic Press, Inc., pp . 295–323  
  4. Nathanael, Berglund. "Reconstrucción de funciones refinables" . Archivado del original el 4 de abril de 2009. Consultado el 24 de diciembre de 2010 .
  5. Thielemann, Henning (2012-01-29). "Cómo refinar funciones polinomiales". arXiv : 1012.2453 [ math.FA ].
  6. Gustafson, Paul; Savir, Nathan; Spears, Ely (14 de noviembre de 2006), "Una caracterización de funciones racionales refinables" (PDF) , American Journal of Undergraduate Research , 5 (3): 11–20 , doi : 10.33697/ajur.2006.021

Véase también

Obtenido de " https://en.wikipedia.org/w/index.php?title=Refinable_function&oldid=1310660740 "