stringtranslate.com

Ondícula

Una wavelet es una oscilación similar a una onda con una amplitud que comienza en cero, aumenta o disminuye y luego regresa a cero una o más veces. Las ondas se denominan "oscilación breve". Se ha establecido una taxonomía de wavelets, basada en el número y dirección de sus pulsos. Las wavelets están dotadas de propiedades específicas que las hacen útiles para el procesamiento de señales .

Onda sísmica

Por ejemplo, se podría crear una wavelet para que tenga una frecuencia de C medio y una duración corta de aproximadamente una décima de segundo. Si esta wavelet se convolucionara con una señal creada a partir de la grabación de una melodía, entonces la señal resultante sería útil para determinar cuándo apareció la nota C media en la canción. Matemáticamente, una wavelet se correlaciona con una señal si una parte de la señal es similar. La correlación es el núcleo de muchas aplicaciones prácticas de wavelets.

Como herramienta matemática, las wavelets se pueden utilizar para extraer información de muchos tipos de datos, incluidas señales de audio e imágenes. Se necesitan conjuntos de wavelets para analizar los datos por completo. Las ondas "complementarias" descomponen una señal sin espacios ni superposiciones, de modo que el proceso de descomposición es matemáticamente reversible. Por lo tanto, los conjuntos de wavelets complementarios son útiles en algoritmos de compresión/descompresión basados ​​en wavelets , donde es deseable recuperar la información original con una pérdida mínima.

En términos formales, esta representación es una representación en serie de ondas de una función integrable al cuadrado con respecto a un conjunto ortonormal completo de funciones básicas , o un conjunto o marco sobrecompleto de un espacio vectorial , para el espacio de Hilbert de funciones integrables al cuadrado . . Esto se logra a través de estados coherentes .

En física clásica , el fenómeno de difracción se describe mediante el principio de Huygens-Fresnel que trata cada punto en un frente de onda en propagación como una colección de ondas esféricas individuales. [1] El patrón de flexión característico es más pronunciado cuando una onda de una fuente coherente (como un láser) encuentra una rendija/apertura que es comparable en tamaño a su longitud de onda . Esto se debe a la adición, o interferencia , de diferentes puntos en el frente de onda (o, de manera equivalente, cada ondícula) que viajan por caminos de diferentes longitudes hasta la superficie de registro. Múltiples aberturas estrechamente espaciadas (por ejemplo, una rejilla de difracción ) pueden dar como resultado un patrón complejo de intensidad variable.

Etimología

La palabra wavelet se ha utilizado durante décadas en el procesamiento de señales digitales y en la geofísica de exploración. [2] Morlet y Grossmann utilizaron la palabra francesa equivalente ondelette, que significa "ola pequeña" , a principios de los años 1980.

Teoría de las ondículas

La teoría wavelet es aplicable a varios temas. Todas las transformadas wavelet pueden considerarse formas de representación de tiempo-frecuencia para señales de tiempo continuo (analógicas) y, por lo tanto, están relacionadas con el análisis armónico . La transformada de ondas discretas (continua en el tiempo) de una señal de tiempo discreto (muestreada) mediante el uso de bancos de filtros de tiempo discreto de configuración diádica (banda de octava) es una aproximación de ondas a esa señal. Los coeficientes de dicho banco de filtros se denominan coeficientes de desplazamiento y escala en la nomenclatura de wavelets. Estos bancos de filtros pueden contener filtros de respuesta de impulso finita (FIR) o de respuesta de impulso infinita (IIR). Las ondas que forman una transformada de ondas continua (CWT) están sujetas al principio de incertidumbre de la teoría de muestreo respectiva del análisis de Fourier: dada una señal con algún evento, no se puede asignar simultáneamente una escala exacta de respuesta de tiempo y frecuencia a ese evento. El producto de las incertidumbres del tiempo y la escala de respuesta en frecuencia tiene un límite inferior. Por lo tanto, en el escalaograma de una transformada wavelet continua de esta señal, tal evento marca una región completa en el plano de la escala de tiempo, en lugar de solo un punto. Además, las bases de ondículas discretas pueden considerarse en el contexto de otras formas del principio de incertidumbre. [3] [4] [5] [6]

Las transformadas Wavelet se dividen en términos generales en tres clases: continuas, discretas y basadas en multiresolución.

Transformadas wavelet continuas (parámetros de escala y cambio continuo)

En las transformadas wavelet continuas , una señal dada de energía finita se proyecta en una familia continua de bandas de frecuencia (o subespacios similares del espacio funcional L p L 2 ( R ) ). Por ejemplo, la señal puede representarse en cada banda de frecuencia de la forma [ f , 2 f ] para todas las frecuencias positivas f > 0. Luego, la señal original puede reconstruirse mediante una integración adecuada de todos los componentes de frecuencia resultantes.

Las bandas de frecuencia o subespacios (subbandas) son versiones escaladas de un subespacio a escala 1. Este subespacio, a su vez, se genera en la mayoría de las situaciones por los cambios de una función generadora ψ en L 2 ( R ), la wavelet madre . Para el ejemplo de la escala de una banda de frecuencia [1, 2] esta función es

función sinc

El subespacio de escala a o banda de frecuencia [1/ a , 2/ a ] es generado por las funciones (a veces llamadas ondas secundarias )

ababR +R.

La proyección de una función x en el subespacio de escala a tiene entonces la forma

coeficientes wavelet

Para el análisis de la señal x , se pueden ensamblar los coeficientes wavelet en un escalaograma de la señal.

Vea una lista de algunas wavelets continuas .

Transformadas wavelet discretas (parámetros de escala y desplazamiento discretos, continuos en el tiempo)

Es computacionalmente imposible analizar una señal utilizando todos los coeficientes wavelet, por lo que uno puede preguntarse si es suficiente elegir un subconjunto discreto del semiplano superior para poder reconstruir una señal a partir de los coeficientes wavelet correspondientes. Uno de esos sistemas es el sistema afín para algunos parámetros reales a > 1, b > 0. El subconjunto discreto correspondiente del semiplano consta de todos los puntos ( a m , nb a m ) con m , n en Z . Las ondas secundarias correspondientes ahora se dan como

Una condición suficiente para la reconstrucción de cualquier señal x de energía finita mediante la fórmula

base ortonormalL 2R

Transformadas wavelet discretas basadas en resolución múltiple (continuas en el tiempo)

Ondícula D4

En cualquier transformada wavelet discretizada, sólo hay un número finito de coeficientes wavelet para cada región rectangular acotada en el semiplano superior. Aún así, cada coeficiente requiere la evaluación de una integral. En situaciones especiales, esta complejidad numérica se puede evitar si las ondas escaladas y desplazadas forman un análisis multiresolución . Esto significa que tiene que existir una función auxiliar , la wavelet padre φ en L 2 ( R ), y que a es un número entero. Una elección típica es a = 2 y b = 1. El par de ondas padre y madre más famoso es la ondícula de 4 pulsaciones de Daubechies . Tenga en cuenta que no todas las ondas discretas ortonormales pueden asociarse a un análisis multiresolución; por ejemplo, la wavelet Journe no admite análisis multiresolución. [7]

A partir de las ondas madre y padre se construyen los subespacios.

De estos se requiere que la secuencia

análisis multiresoluciónL 2W mV mV m −1

Por analogía con el teorema de muestreo se puede concluir que el espacio V m con una distancia de muestreo de 2 m cubre más o menos la banda base de frecuencia de 0 a 1/2 m -1 . Como complemento ortogonal, W m cubre aproximadamente la banda [1/2 m −1 , 1/2 m ].

De aquellas inclusiones y relaciones de ortogonalidad, especialmente , se desprende la existencia de sucesiones y que satisfacen las identidades

ecuación de refinamientotransformada wavelet rápida

Del análisis multiresolución se deriva la descomposición ortogonal del espacio L 2 como

Ondas causales del tiempo

Para procesar señales temporales en tiempo real, es esencial que los filtros wavelet no accedan a valores de señal del futuro y que se puedan obtener latencias temporales mínimas. Szu et al [8] y Lindeberg, [9] han desarrollado representaciones de wavelets causales en el tiempo, y este último método también implica una implementación recursiva en el tiempo eficiente en memoria.

Ola madre

Para aplicaciones prácticas, y por razones de eficiencia, se prefieren funciones continuamente diferenciables con soporte compacto como wavelet (funciones) madre (prototipo). Sin embargo, para satisfacer los requisitos analíticos (en el WT continuo) y en general por razones teóricas, se eligen las funciones wavelet de un subespacio del espacio. Este es el espacio de funciones medibles de Lebesgue que son absolutamente integrables y integrables al cuadrado en el sentido de que

Estar en este espacio asegura que uno pueda formular las condiciones de media cero y norma cuadrática uno:

Para que ψ sea una wavelet para la transformada wavelet continua (ver allí para la declaración exacta), la wavelet madre debe satisfacer un criterio de admisibilidad (en términos generales, una especie de media diferenciabilidad) para obtener una transformada establemente invertible.

Para la transformada wavelet discreta , se necesita al menos la condición de que la serie wavelet sea una representación de la identidad en el espacio L 2 ( R ). La mayoría de las construcciones de WT discreta hacen uso del análisis multiresolución , que define la wavelet mediante una función de escala. Esta función de escala en sí misma es una solución a una ecuación funcional.

En la mayoría de las situaciones es útil restringir ψ para que sea una función continua con un mayor número M de momentos de fuga, es decir, para todo número entero m < M

La wavelet madre se escala (o dilata) por un factor de a y se traslada (o desplaza) por un factor de b para dar (según la formulación original de Morlet):

Para el WT continuo, el par ( a , b ) varía en todo el semiplano R + × R ; para el WT discreto este par varía en un subconjunto discreto del mismo, que también se denomina grupo afín .

Estas funciones a menudo se denominan incorrectamente funciones básicas de la transformada (continua). De hecho, al igual que en la transformada continua de Fourier, no existe ninguna base en la transformada wavelet continua. La interpretación tiempo-frecuencia utiliza una formulación sutilmente diferente (según Delprat).

Restricción:

  1. cuando a 1 = a y b 1 = b ,
  2. tiene un intervalo de tiempo finito

Comparaciones con la transformada de Fourier (tiempo continuo)

La transformada wavelet a menudo se compara con la transformada de Fourier , en la que las señales se representan como una suma de sinusoides. De hecho, la transformada de Fourier puede verse como un caso especial de la transformada wavelet continua con la elección de la wavelet madre . La principal diferencia en general es que las wavelets se localizan tanto en tiempo como en frecuencia, mientras que la transformada estándar de Fourier sólo se localiza en frecuencia . La transformada de Fourier de corto tiempo (STFT) es similar a la transformada wavelet, en el sentido de que también está localizada en el tiempo y la frecuencia, pero existen problemas con el compromiso de resolución entre frecuencia y tiempo.

En particular, suponiendo una región de ventana rectangular, se puede pensar en STFT como una transformación con un núcleo ligeramente diferente.

uel teorema de Parseval
u

y el cuadrado del soporte espectral de la ventana actuando sobre una frecuencia

La multiplicación con una ventana rectangular en el dominio del tiempo corresponde a la convolución con una función en el dominio de la frecuencia, lo que da como resultado artefactos de timbre espurios para ventanas temporales cortas/localizadas. Con la Transformada de Fourier de tiempo continuo, y esta convolución es con una función delta en el espacio de Fourier, lo que da como resultado la verdadera transformada de Fourier de la señal . La función de ventana puede ser algún otro filtro apodizante , como un gaussiano . La elección de la función de ventana afectará el error de aproximación relativo a la verdadera transformada de Fourier.

El producto tiempo-ancho de banda de una determinada celda de resolución no puede excederse con el STFT. Todos los elementos básicos de STFT mantienen un soporte espectral y temporal uniforme para todos los cambios o compensaciones temporales, logrando así una resolución igual en el tiempo para frecuencias más bajas y más altas. La resolución está determinada puramente por el ancho de muestreo.

Por el contrario, las propiedades multiresolución de la transformada wavelet permiten grandes soportes temporales para frecuencias más bajas mientras mantienen anchos temporales cortos para frecuencias más altas mediante las propiedades de escala de la transformada wavelet. Esta propiedad extiende el análisis de tiempo-frecuencia convencional al análisis de escala de tiempo. [10]

Átomos de tiempo-frecuencia STFT (izquierda) y átomos de escala de tiempo DWT (derecha). Los átomos de tiempo-frecuencia son cuatro funciones básicas diferentes utilizadas para el STFT (es decir, se requieren cuatro transformadas de Fourier separadas ). Los átomos de escala de tiempo del DWT logran anchos temporales pequeños para frecuencias altas y buenos anchos temporales para frecuencias bajas con un conjunto de bases de transformación única .

La transformada wavelet discreta es menos compleja desde el punto de vista computacional y requiere un tiempo O( N ) en comparación con O( N  log  N ) para la transformada rápida de Fourier . Esta ventaja computacional no es inherente a la transformada, pero refleja la elección de una división logarítmica de frecuencia, en contraste con las divisiones de frecuencia equiespaciadas de la FFT (transformada rápida de Fourier), que utiliza las mismas funciones básicas que la DFT (transformada discreta de Fourier). . [11] También es importante tener en cuenta que esta complejidad solo se aplica cuando el tamaño del filtro no tiene relación con el tamaño de la señal. Una wavelet sin soporte compacto como la wavelet de Shannon requeriría O( N 2 ). (Por ejemplo, también existe una transformada logarítmica de Fourier con complejidad O( N ), pero la señal original debe muestrearse logarítmicamente en el tiempo, lo que sólo es útil para ciertos tipos de señales. [12] )

Definición de wavelet

Una wavelet (o una familia de wavelets) se puede definir de varias maneras:

Filtro de escala

Una wavelet ortogonal está completamente definida por el filtro de escala: un filtro de respuesta de impulso finito (FIR) de paso bajo de longitud 2 N y suma 1. En las wavelets biortogonales , se definen filtros separados de descomposición y reconstrucción.

Para el análisis con wavelets ortogonales, el filtro de paso alto se calcula como el filtro de espejo en cuadratura del paso bajo, y los filtros de reconstrucción son el tiempo inverso de los filtros de descomposición.

Las ondas Daubechies y Symlet se pueden definir mediante el filtro de escala.

Función de escala

Las wavelets se definen por la función wavelet ψ( t ) (es decir, la wavelet madre) y la función de escala φ( t ) (también llamada wavelet padre) en el dominio del tiempo.

La función wavelet es, en efecto, un filtro de paso de banda y un escalamiento que para cada nivel reduce a la mitad su ancho de banda. Esto crea el problema de que para cubrir todo el espectro se necesitaría un número infinito de niveles. La función de escala filtra el nivel más bajo de la transformación y garantiza que se cubra todo el espectro. Véase [13] para una explicación detallada.

Para una wavelet con soporte compacto, φ( t ) puede considerarse de longitud finita y es equivalente al filtro de escala g .

Las wavelets de Meyer se pueden definir mediante funciones de escala

Función wavelet

La wavelet solo tiene una representación en el dominio del tiempo como la función wavelet ψ( t ).

Por ejemplo, las wavelets del sombrero mexicano se pueden definir mediante una función wavelet. Vea una lista de algunas ondas continuas .

Historia

El desarrollo de las wavelets puede vincularse a varias líneas de pensamiento distintas, comenzando con el trabajo de Haar a principios del siglo XX. Un trabajo posterior de Dennis Gabor produjo átomos de Gabor (1946), que se construyen de manera similar a las ondas y se aplican con propósitos similares.

Las contribuciones notables a la teoría de las wavelets desde entonces pueden atribuirse al descubrimiento de Zweig de la transformada wavelet continua (CWT) en 1975 (originalmente llamada transformada coclear y descubierta mientras estudiaba la reacción del oído al sonido), [14] Pierre Goupillaud, La formulación de Grossmann y Morlet de lo que ahora se conoce como CWT (1982), los primeros trabajos de Jan-Olov Strömberg sobre wavelets discretas (1983), el banco de filtros no ortogonales de 5/3 grifos de Le Gall-Tabatabai (LGT) con fase lineal (1988), [15] [16] [17] Wavelets ortogonales con soporte compacto de Ingrid Daubechies (1988), marco multiresolución no ortogonal de Mallat (1989), Binomial QMF de Ali Akansu (1990), Nathalie La interpretación tiempo-frecuencia de Delprat del CWT (1991), la transformada wavelet armónica de Newland (1993) y la partición de conjuntos en árboles jerárquicos (SPIHT) desarrollada por Amir Said con William A. Pearlman en 1996. [18]

El estándar JPEG 2000 fue desarrollado entre 1997 y 2000 por un comité del Grupo Conjunto de Expertos en Fotografía (JPEG) presidido por Touradj Ebrahimi (más tarde presidente de JPEG). [19] A diferencia del algoritmo DCT utilizado por el formato JPEG original , JPEG 2000 utiliza algoritmos de transformada wavelet discreta (DWT). Utiliza la transformada wavelet CDF 9/7 (desarrollada por Ingrid Daubechies en 1992) para su algoritmo de compresión con pérdida , y el banco de filtros de tiempo discreto Le Gall-Tabatabai (LGT) 5/3 (desarrollado por Didier Le Gall y Ali J. Tabatabai en 1988) por su algoritmo de compresión sin pérdidas . [20] La tecnología JPEG 2000 , que incluye la extensión Motion JPEG 2000 , fue seleccionada como estándar de codificación de vídeo para cine digital en 2004. [21]

Línea de tiempo

transformaciones wavelet

Una wavelet es una función matemática que se utiliza para dividir una función determinada o una señal de tiempo continuo en diferentes componentes de escala. Generalmente se puede asignar un rango de frecuencia a cada componente de la escala. Luego, cada componente de escala se puede estudiar con una resolución que coincida con su escala. Una transformada wavelet es la representación de una función mediante wavelets. Las ondas son copias escaladas y traducidas (conocidas como "ondas hijas") de una forma de onda oscilante de longitud finita o de rápida decadencia (conocida como "ondas madre"). Las transformadas Wavelet tienen ventajas sobre las transformadas de Fourier tradicionales para representar funciones que tienen discontinuidades y picos agudos, y para deconstruir y reconstruir con precisión señales finitas, no periódicas y/o no estacionarias .

Las transformadas wavelet se clasifican en transformadas wavelet discretas (DWT) y transformadas wavelet continuas (CWT). Tenga en cuenta que tanto DWT como CWT son transformaciones de tiempo continuo (analógicas). Se pueden utilizar para representar señales de tiempo continuo (analógicas). Los CWT operan en todas las escalas y traducciones posibles, mientras que los DWT utilizan un subconjunto específico de valores de escala y traducción o cuadrícula de representación.

Hay una gran cantidad de transformadas wavelet, cada una de las cuales es adecuada para diferentes aplicaciones. Para obtener una lista completa, consulte la lista de transformaciones relacionadas con wavelets, pero las más comunes se enumeran a continuación:

Transformaciones generalizadas

Hay una serie de transformadas generalizadas de las cuales la transformada wavelet es un caso especial. Por ejemplo, Yosef Joseph Segman introdujo la escala en el grupo de Heisenberg , dando lugar a un espacio de transformación continuo que es función del tiempo, la escala y la frecuencia. El CWT es un corte bidimensional a través del volumen de frecuencia-escala de tiempo 3D resultante.

Otro ejemplo de una transformada generalizada es la transformada chirplet en la que el CWT también es un corte bidimensional a través de la transformada chirplet.

Un área de aplicación importante para las transformadas generalizadas involucra sistemas en los que la resolución de alta frecuencia es crucial. Por ejemplo, las transformadas ópticas de electrones de campo oscuro intermedias entre el espacio directo y recíproco se han utilizado ampliamente en el análisis armónico de la agrupación de átomos, es decir, en el estudio de cristales y defectos cristalinos . [22] Ahora que los microscopios electrónicos de transmisión son capaces de proporcionar imágenes digitales con información a escala picométrica sobre la periodicidad atómica en nanoestructuras de todo tipo, se amplía el rango de aplicaciones de reconocimiento de patrones [23] y deformación [24] / metrología [25] para transformadas intermedias. con resolución de alta frecuencia (como brushlets [26] y ridgelets [27] ) está creciendo rápidamente.

La transformada wavelet fraccional (FRWT) es una generalización de la transformada wavelet clásica en los dominios de la transformada fraccional de Fourier. Esta transformada es capaz de proporcionar información en el dominio de tiempo y fraccional simultáneamente y representar señales en el plano de frecuencia fraccional de tiempo. [28]

Aplicaciones

Generalmente, se utiliza una aproximación a DWT para la compresión de datos si ya se ha muestreado una señal, y el CWT para el análisis de la señal . [29] Por lo tanto, la aproximación DWT se usa comúnmente en ingeniería e informática, [30] y la CWT en investigación científica. [31]

Al igual que otras transformaciones, las transformaciones wavelet se pueden utilizar para transformar datos y luego codificar los datos transformados, lo que da como resultado una compresión efectiva. Por ejemplo, JPEG 2000 es un estándar de compresión de imágenes que utiliza wavelets biortogonales. Esto significa que aunque el marco está sobrecompleto, es un marco ajustado (ver tipos de marcos de un espacio vectorial ), y las mismas funciones de marco (excepto la conjugación en el caso de wavelets complejas) se utilizan tanto para el análisis como para la síntesis, es decir , tanto en la transformada directa como en la inversa. Para obtener más información, consulte compresión wavelet .

Un uso relacionado es el de suavizar/eliminar ruido de datos basándose en el umbral del coeficiente de wavelet, también llamado contracción de wavelet. Umbralizando de forma adaptativa los coeficientes wavelet que corresponden a componentes de frecuencia no deseados, se pueden realizar operaciones de suavizado y/o eliminación de ruido.

Las transformadas Wavelet también se están empezando a utilizar para aplicaciones de comunicación. Wavelet OFDM es el esquema de modulación básico utilizado en HD-PLC (una tecnología de comunicaciones por línea eléctrica desarrollada por Panasonic ), y en uno de los modos opcionales incluidos en el estándar IEEE 1901 . Wavelet OFDM puede lograr muescas más profundas que FFT OFDM tradicional, y wavelet OFDM no requiere un intervalo de guarda (lo que generalmente representa una sobrecarga significativa en los sistemas FFT OFDM). [32]

Como representación de una señal.

A menudo, las señales se pueden representar bien como una suma de sinusoides. Sin embargo, considere una señal no continua con una discontinuidad abrupta; Esta señal todavía se puede representar como una suma de sinusoides, pero requiere un número infinito, lo cual es una observación conocida como fenómeno de Gibbs . Esto, entonces, requiere un número infinito de coeficientes de Fourier, lo que no es práctico para muchas aplicaciones, como la compresión. Las wavelets son más útiles para describir estas señales con discontinuidades debido a su comportamiento localizado en el tiempo (tanto la transformada de Fourier como la wavelet están localizadas en frecuencia, pero las wavelets tienen una propiedad adicional de localización en el tiempo). Debido a esto, muchos tipos de señales en la práctica pueden no ser dispersas en el dominio de Fourier, pero muy dispersas en el dominio wavelet. Esto es particularmente útil en la reconstrucción de señales, especialmente en el campo recientemente popular de la detección comprimida . (Tenga en cuenta que la transformada de Fourier de tiempo corto (STFT) también está localizada en el tiempo y la frecuencia, pero a menudo hay problemas con el equilibrio entre la resolución frecuencia-tiempo. Las wavelets son mejores representaciones de señales debido al análisis de resolución múltiple ).

Esto motiva por qué las transformadas wavelet se están adoptando ahora para una gran cantidad de aplicaciones, reemplazando a menudo a la transformada de Fourier convencional . Muchas áreas de la física han visto este cambio de paradigma, incluida la dinámica molecular , la teoría del caos , [33] cálculos ab initio , la astrofísica , el análisis de datos transitorios de ondas gravitacionales , [34] [35] la localización de la matriz de densidad , la sismología , la óptica , la turbulencia y la cuántica. mecánica . Este cambio también se ha producido en el procesamiento de imágenes , EEG , EMG , [36] análisis de ECG , ritmos cerebrales , análisis de ADN , análisis de proteínas , climatología , análisis de la respuesta sexual humana, [37] procesamiento general de señales , reconocimiento de voz , acústica, señales de vibración, [38] gráficos por computadora , análisis multifractal y codificación dispersa . En visión por computadora y procesamiento de imágenes , la noción de representación espacial de escala y operadores derivados gaussianos se considera una representación canónica de múltiples escalas.

Eliminación de ruido de ondas

Eliminación de ruido de señal mediante umbralización de transformada wavelet

Supongamos que medimos una señal ruidosa , donde representa la señal y representa el ruido. Supongamos que tiene una representación escasa en una determinada base de ondículas, y

Sea la transformada wavelet de , donde es la transformada wavelet del componente de señal y es la transformada wavelet del componente de ruido.

La mayoría de los elementos son 0 o cercanos a 0, y

Como es ortogonal, el problema de estimación equivale a la recuperación de una señal en ruido gaussiano iid . Como es escaso, un método consiste en aplicar un modelo de mezcla gaussiana para .

Supongamos un prior , donde es la varianza de los coeficientes "significativos" y es la varianza de los coeficientes "insignificantes".

Entonces , se llama factor de contracción, que depende de las varianzas anteriores y . Al establecer en cero los coeficientes que caen por debajo de un umbral de contracción, una vez que se aplica la transformación inversa, se pierde una pequeña cantidad esperada de señal debido al supuesto de escasez. Se espera que los coeficientes más grandes representen principalmente la señal debido a la escasez, y estadísticamente se espera que muy poca señal, aunque la mayor parte del ruido, esté representada en coeficientes de magnitud tan baja... por lo tanto, se espera que la operación de puesta a cero Elimina la mayor parte del ruido y poca señal. Normalmente, los coeficientes superiores al umbral no se modifican durante este proceso. Algunos algoritmos para la eliminación de ruido basados ​​en wavelets también pueden atenuar coeficientes más grandes, basándose en una estimación estadística de la cantidad de ruido que se espera que se elimine mediante dicha atenuación.

Por último, aplique la transformada wavelet inversa para obtener

Red climática multiescala

Agarwal et al. propuso métodos avanzados lineales [39] y no lineales [40] basados ​​en wavelets para construir e investigar el clima como redes complejas en diferentes escalas de tiempo. Las redes climáticas construidas utilizando conjuntos de datos de TSM en diferentes escalas de tiempo afirmaron que el análisis de procesos climáticos a múltiples escalas basado en ondas promete una mejor comprensión de la dinámica del sistema que puede pasarse por alto cuando los procesos se analizan en una sola escala de tiempo [41]

Lista de ondas

Ondas discretas

Ondas continuas

Valor real

Valor complejo

Ver también

Referencias

  1. ^ Comunicaciones inalámbricas: principios y práctica, serie de tecnologías emergentes e ingeniería de comunicaciones de Prentice Hall, T. S. Rappaport, Prentice Hall, 2002, p. 126.
  2. ^ Ricker, normando (1953). "Contracción de ondas, expansión de ondas y control de la resolución sísmica". Geofísica . 18 (4): 769–792. Código bibliográfico : 1953Geop...18..769R. doi :10.1190/1.1437927.
  3. ^ Meyer, Yves (1992), Wavelets y operadores, Cambridge, Reino Unido: Cambridge University Press, ISBN 0-521-42000-8 
  4. ^ Chui, Charles K. (1992), Introducción a las Wavelets, San Diego, CA: Academic Press, ISBN 0-12-174584-8 
  5. ^ Daubechies, Ingrid. (1992), Diez conferencias sobre Wavelets, SIAM, ISBN 978-0-89871-274-2 
  6. ^ Akansu, Ali N.; Haddad, Richard A. (1992), Descomposición de señales multiresolución: transformadas, subbandas y ondículas, Boston, MA: Academic Press, ISBN 978-0-12-047141-6 
  7. ^ Larson, David R. (2007), Análisis y aplicaciones de wavelets (ver: sistemas unitarios y conjuntos de wavelets) , Appl. Número. Armon. Anal., Birkhäuser, págs. 143-171
  8. ^ Szu, Harold H.; Telfer, Brian A.; Lohmann, Adolf W. (1992). "Transformada wavelet analítica causal". Ingeniería Óptica . 31 (9): 1825. Bibcode : 1992OptEn..31.1825S. doi :10.1117/12.59911.
  9. ^ Lindeberg, T. (23 de enero de 2023). "Una representación de espacio de escala covariante de escala causal y recursiva en el tiempo de señales temporales y tiempo pasado". Cibernética biológica . 117 (1–2): 21–59. doi : 10.1007/s00422-022-00953-6 . PMC 10160219 . PMID  36689001. 
  10. ^ Mallat, Stéphane. "Un recorrido por ondas sobre el procesamiento de señales. 1998". 250-252.
  11. ^ La guía para científicos e ingenieros sobre el procesamiento de señales digitales por Steven W. Smith, Ph.D. capítulo 8 ecuación 8-1: http://www.dspguide.com/ch8/4.htm
  12. ^ Haines, VG. V.; Jones, Alan G. (1988). "Transformada logarítmica de Fourier" (PDF) . Revista geofísica (92): 171–178. doi :10.1111/j.1365-246X.1988.tb01131.x. S2CID  9720759.
  13. ^ "Una guía realmente sencilla sobre wavelets: PolyValens". www.polyvalens.com .
  14. ^ Weisstein, Eric W. "Zweig, George - del mundo de la biografía científica de Eric Weisstein". scienceworld.wolfram.com . Consultado el 20 de octubre de 2021 .
  15. ^ Sullivan, Gary (8 a 12 de diciembre de 2003). "Características generales y consideraciones de diseño para la codificación de vídeo de subbanda temporal". UIT-T . Grupo de expertos en codificación de vídeo . Consultado el 13 de septiembre de 2019 .
  16. ^ Bovik, Alan C. (2009). La guía esencial para el procesamiento de vídeos. Prensa académica . pag. 355.ISBN 9780080922508.
  17. ^ Gall, Didier Le; Tabatabai, Ali J. (1988). "Codificación de subbanda de imágenes digitales mediante filtros de núcleo cortos simétricos y técnicas de codificación aritmética". ICASSP-88., Conferencia Internacional sobre Acústica, Habla y Procesamiento de Señales . págs. 761–764 vol.2. doi :10.1109/ICASSP.1988.196696. S2CID  109186495.
  18. ^ Dijo, Amir; Pearlman, William A. (junio de 1996). "Un nuevo códec de imágenes rápido y eficiente basado en la partición de conjuntos en árboles jerárquicos". Transacciones IEEE sobre circuitos y sistemas para tecnología de vídeo . 6 (3): 243–250. doi : 10.1109/76.499834. ISSN  1051-8215.
  19. ^ Taubman, David; Marcelino, Miguel (2012). Fundamentos, estándares y práctica de la compresión de imágenes JPEG2000: Fundamentos, estándares y práctica de la compresión de imágenes. Medios de ciencia y negocios de Springer . ISBN 9781461507994.
  20. ^ Unser, M.; Blu, T. (2003). "Propiedades matemáticas de los filtros wavelet JPEG2000" (PDF) . Transacciones IEEE sobre procesamiento de imágenes . 12 (9): 1080–1090. Código Bib : 2003ITIP...12.1080U. doi :10.1109/TIP.2003.812329. PMID  18237979. S2CID  2765169. Archivado desde el original (PDF) el 13 de octubre de 2019.
  21. ^ Swartz, Charles S. (2005). Comprensión del cine digital: un manual profesional. Taylor y Francisco . pag. 147.ISBN 9780240806174.
  22. ^ P. Hirsch, A. Howie, R. Nicholson, DW Pashley y MJ Whelan (1965/1977) Microscopía electrónica de cristales finos (Butterworths, Londres/Krieger, Malabar FLA) ISBN 0-88275-376-2 
  23. ^ P. Fraundorf, J. Wang, E. Mandell y M. Rose (2006) Cuadros digitales de campo oscuro, Microscopía y microanálisis 12 : S2, 1010-1011 (cf. arXiv:cond-mat/0403017)
  24. ^ Hÿtch, MJ; Snoeck, E.; Kilaas, R. (1998). "Medición cuantitativa de campos de desplazamiento y deformación a partir de micrografías HRTEM". Ultramicroscopía . 74 (3): 131-146. doi :10.1016/s0304-3991(98)00035-7.
  25. ^ Martin Rose (2006) Mediciones de espaciado de franjas de red en imágenes HRTEM mediante descomposición digital en campo oscuro (Tesis de maestría en Física, U. Missouri - St. Louis)
  26. ^ FG Meyer y RR Coifman (1997) Análisis armónico computacional y aplicado 4 :147.
  27. ^ AG Flesia, H. Hel-Or , A. Averbuch, EJ Candes , RR Coifman y DL Donoho (2001) Implementación digital de paquetes ridgelet (Academic Press, Nueva York).
  28. ^ Shi, J.; Zhang, N.-T.; Liu, XP.-P. (2011). "Una novedosa transformada wavelet fraccionaria y sus aplicaciones". Ciencia. Inf. China. Ciencia . 55 (6): 1270-1279. doi :10.1007/s11432-011-4320-x. S2CID  255201598.
  29. ^ AN Akansu, WA Serdijn e IW Selesnick, Aplicaciones emergentes de wavelets: una revisión, Comunicación física, Elsevier, vol. 3, número 1, págs. 1-18, marzo de 2010.
  30. ^ Lyakhov, Pavel; Semyonova, Nataliya; Nagornov, Nikolay; Bergerman, Maxim; Abdulsalyamova, Albina (14 de noviembre de 2023). "Procesamiento de imágenes Wavelet de alta velocidad utilizando el método Winograd con reducción de resolución". Matemáticas . 11 (22): 4644. doi : 10.3390/math11224644 . ISSN  2227-7390. Las wavelets se utilizan activamente para resolver una amplia gama de problemas de procesamiento de imágenes en diversos campos de la ciencia y la tecnología, por ejemplo, eliminación de ruido, reconstrucción, análisis y análisis y procesamiento de vídeo. Los métodos de procesamiento de wavelets se basan en la transformada de wavelets discreta utilizando filtrado digital 1D.
  31. ^ Dong, Liang; Zhang, Shaohua; Gan, Tiansiyu; Qiu, Yan; Canción, Qinfeng; Zhao, Yongtao (1 de diciembre de 2023). "Análisis de las características de frecuencia del potencial tubería-suelo bajo interferencia de corriente parásita del metro utilizando el método de transformada wavelet continua". Materiales de Construcción y Construcción . 407 : 133453. doi : 10.1016/j.conbuildmat.2023.133453. ISSN  0950-0618. S2CID  263317973.
  32. ^ Stefano Galli; O. Logvinov (julio de 2008). "Desarrollos recientes en la estandarización de las comunicaciones por líneas eléctricas dentro del IEEE". Revista de comunicaciones IEEE . 46 (7): 64–71. doi :10.1109/MCOM.2008.4557044. S2CID  2650873.Una descripción general de la propuesta P1901 PHY/MAC.
  33. ^ Wotherspoon, T.; et., al. (2009). "Adaptación al borde del caos con retroalimentación de ondas aleatorias". J. Física. química . 113 (1): 19-22. Código Bib : 2009JPCA..113...19W. doi :10.1021/jp804420g. PMID  19072712.
  34. ^ Abbott, Benjamín P.; et al. (Colaboración científica LIGO y Colaboración Virgo) (2016). "Observación del transitorio de ondas gravitacionales GW150914 con suposiciones mínimas". Física. Rev. D. 93 (12): 122004. arXiv : 1602.03843 . Código Bib : 2016PhRvD..93l2004A. doi : 10.1103/PhysRevD.93.122004. S2CID  119313566.
  35. ^ V Necula, S Klimenko y G Mitselmakher (2012). "Análisis transitorio con rápida transformada tiempo-frecuencia de Wilson-Daubechies". Revista de Física: Serie de conferencias . 363 (1): 012032. Código bibliográfico : 2012JPhCS.363a2032N. doi : 10.1088/1742-6596/363/1/012032 .
  36. ^ J. Rafiee y col. Extracción de características de señales EMG de antebrazo para prótesis, Sistemas expertos con aplicaciones 38 (2011) 4058–67.
  37. ^ J. Rafiee y col. Respuestas sexuales femeninas utilizando técnicas de procesamiento de señales, The Journal of Sexual Medicine 6 (2009) 3086–96. (pdf)
  38. ^ Rafiee, J.; Tse, Peter W. (2009). "Uso de autocorrelación en coeficientes wavelet para diagnóstico de fallas". Sistemas Mecánicos y Procesamiento de Señales . 23 (5): 1554–72. Código Bib : 2009MSSP...23.1554R. doi :10.1016/j.ymssp.2009.02.008.
  39. ^ Agarwal, Ankit; Maheswaran, Rathinasamy; Marwan, Norberto; César, Levke; Kurths, Jürgen (noviembre de 2018). "Medida de similitud multiescala basada en wavelets para redes complejas" (PDF) . La revista física europea B. 91 (11): 296. Código bibliográfico : 2018EPJB...91..296A. doi :10.1140/epjb/e2018-90460-6. eISSN  1434-6036. ISSN  1434-6028. S2CID  125557123.
  40. ^ Agarwal, Ankit; Marwan, Norberto; Rathinasamy, Maheswaran; Merz, Bruno; Kurths, Jürgen (13 de octubre de 2017). "Análisis de sincronización de eventos a múltiples escalas para desentrañar procesos climáticos: un enfoque basado en ondas". Procesos no lineales en geofísica . 24 (4): 599–611. Código Bib : 2017NPGeo..24..599A. doi : 10.5194/npg-24-599-2017 . eISSN  1607-7946. S2CID  28114574.
  41. ^ Agarwal, Ankit; César, Levke; Marwan, Norberto; Maheswaran, Rathinasamy; Merz, Bruno; Kurths, Jürgen (19 de junio de 2019). "Identificación y caracterización de teleconexiones en red a diferentes escalas". Informes científicos . 9 (1): 8808. Código bibliográfico : 2019NatSR...9.8808A. doi :10.1038/s41598-019-45423-5. eISSN  2045-2322. PMC 6584743 . PMID  31217490. 
  42. ^ Caja de herramientas de Matlab - URL: http://matlab.izmiran.ru/help/toolbox/wavelet/ch06_a32.html
  43. ^ Erik Hjelmås (21 de enero de 1999) URL de Gabor Wavelets : http://www.ansatt.hig.no/erikh/papers/scia99/node6.html

Otras lecturas

enlaces externos