stringtranslate.com

Función de ventana

Una función de ventana popular, la ventana de Hann . Las funciones de ventana más populares son curvas similares en forma de campana.

En procesamiento de señales y estadística , una función de ventana (también conocida como función de apodización o función de reducción gradual [1] ) es una función matemática que tiene valor cero fuera de algún intervalo elegido . Normalmente, las funciones de ventana son simétricas alrededor de la mitad del intervalo, se acercan a un máximo en la mitad y se estrechan desde la mitad. Matemáticamente, cuando otra función o forma de onda/secuencia de datos se "multiplica" por una función de ventana, el producto también tiene valor cero fuera del intervalo: todo lo que queda es la parte donde se superponen, la "vista a través de la ventana". De manera equivalente, y en la práctica real, el segmento de datos dentro de la ventana primero se aísla y luego solo esos datos se multiplican por los valores de la función de la ventana. Por tanto, la reducción gradual, no la segmentación, es el objetivo principal de las funciones de ventana.

Las razones para examinar segmentos de una función más larga incluyen la detección de eventos transitorios y el promedio temporal de espectros de frecuencia. La duración de los segmentos está determinada en cada aplicación por requisitos como resolución de tiempo y frecuencia. Pero ese método también cambia el contenido de frecuencia de la señal mediante un efecto llamado fuga espectral . Las funciones de ventana nos permiten distribuir espectralmente la fuga de diferentes formas, según las necesidades de la aplicación particular. Hay muchas opciones detalladas en este artículo, pero muchas de las diferencias son tan sutiles que resultan insignificantes en la práctica.

En aplicaciones típicas, las funciones de ventana utilizadas son curvas "en forma de campana", suaves y no negativas. [2] También se pueden utilizar rectángulos, triángulos y otras funciones. Una definición más general de funciones de ventana no requiere que sean idénticamente cero fuera de un intervalo, siempre y cuando el producto de la ventana multiplicado por su argumento sea integrable al cuadrado y, más específicamente, que la función vaya lo suficientemente rápido hacia cero. [3]

Aplicaciones

Las funciones de ventana se utilizan en análisis /modificación/ resíntesis espectral , [4] el diseño de filtros de respuesta de impulso finito , fusionando conjuntos de datos multiescala y multidimensionales, [5] [6] así como en la formación de haces y el diseño de antenas .

Figura 2: La ventana de una sinusoide provoca una fuga espectral. Se produce la misma cantidad de fuga ya sea que haya un número de ciclos entero (azul) o no entero (rojo) dentro de la ventana (filas 1 y 2). Cuando se muestrea y se ventana la sinusoide, su transformada de Fourier de tiempo discreto (DTFT) también exhibe el mismo patrón de fuga (filas 3 y 4). Pero cuando el DTFT se muestra escasamente, en un cierto intervalo, es posible (dependiendo de su punto de vista): (1) evitar la fuga, o (2) crear la ilusión de que no hay fuga. Para el caso de la DTFT azul, esas muestras son las salidas de la transformada discreta de Fourier (DFT). El DTFT rojo tiene el mismo intervalo de cruces por cero, pero las muestras de DFT se encuentran entre ellos y se revela la fuga.

Análisis espectral

La transformada de Fourier de la función cos( ωt ) es cero, excepto en la frecuencia ± ω . Sin embargo, muchas otras funciones y formas de onda no tienen transformadas de forma cerrada convenientes. Alternativamente, uno podría estar interesado en su contenido espectral sólo durante un período de tiempo determinado.

En cualquier caso, la transformada de Fourier (o una transformada similar) se puede aplicar en uno o más intervalos finitos de la forma de onda. En general, la transformada se aplica al producto de la forma de onda y una función de ventana. Cualquier ventana (incluida la rectangular) afecta la estimación espectral calculada por este método.

Diseño de filtro

Las ventanas se utilizan a veces en el diseño de filtros digitales , en particular para convertir una respuesta de impulso "ideal" de duración infinita, como una función sinc , en un diseño de filtro de respuesta de impulso finita (FIR). Esto se llama método de ventana . [7] [8] [9]

Estadísticas y ajuste de curvas.

Las funciones de ventana se utilizan a veces en el campo del análisis estadístico para restringir el conjunto de datos que se analizan a un rango cercano a un punto determinado, con un factor de ponderación que disminuye el efecto de los puntos más alejados de la parte de la curva que se ajusta. En el campo del análisis bayesiano y el ajuste de curvas , esto a menudo se denomina núcleo .

Aplicaciones de ventanas rectangulares

Análisis de transitorios

Al analizar una señal transitoria en análisis modal , como un impulso, una respuesta de choque, una ráfaga sinusoidal, una ráfaga de chirrido o una ráfaga de ruido, donde la distribución de energía versus tiempo es extremadamente desigual, la ventana rectangular puede ser la más apropiada. Por ejemplo, cuando la mayor parte de la energía se encuentra al comienzo de la grabación, una ventana no rectangular atenúa la mayor parte de la energía, degradando la relación señal-ruido. [10]

Análisis armónico

Se podría desear medir el contenido armónico de una nota musical de un instrumento particular o la distorsión armónica de un amplificador a una frecuencia determinada. Refiriéndose nuevamente a la Figura 2 , podemos observar que no hay fugas en un conjunto discreto de frecuencias relacionadas armónicamente muestreadas mediante la transformada discreta de Fourier (DFT). (Los nulos espectrales son en realidad cruces por cero, que no se pueden mostrar en una escala logarítmica como esta). Esta propiedad es exclusiva de la ventana rectangular y debe configurarse adecuadamente para la frecuencia de la señal, como se describió anteriormente.

Ventanas superpuestas

Cuando la longitud de un conjunto de datos a transformar es mayor de lo necesario para proporcionar la resolución de frecuencia deseada, una práctica común es subdividirlo en conjuntos más pequeños y ventanarlos individualmente. Para mitigar la "pérdida" en los bordes de la ventana, los conjuntos individuales pueden superponerse en el tiempo. Véase el método de Welch de análisis espectral de potencia y la transformada de coseno discreta modificada .

Ventanas bidimensionales

Las ventanas bidimensionales se utilizan comúnmente en el procesamiento de imágenes para reducir las altas frecuencias no deseadas en la transformada de Fourier de la imagen. [11] Pueden construirse a partir de ventanas unidimensionales en cualquiera de dos formas. [12] La forma separable es trivial de calcular. La forma radial , que involucra el radio , es isotrópica , independiente de la orientación de los ejes de coordenadas. Sólo la función gaussiana es separable e isotrópica. [13] Las formas separables de todas las demás funciones de ventana tienen esquinas que dependen de la elección de los ejes de coordenadas. La isotropía/ anisotropía de una función de ventana bidimensional es compartida por su transformada de Fourier bidimensional. La diferencia entre las formas separable y radial es similar al resultado de la difracción de aberturas rectangulares versus circulares, que se pueden visualizar en términos del producto de dos funciones sinc versus una función de Airy , respectivamente.

Ejemplos de funciones de ventana

Convenciones :

El escaso muestreo de una transformada de Fourier en tiempo discreto (DTFT), como las DFT en la Fig. 2, solo revela la fuga en los contenedores de DFT desde una sinusoide cuya frecuencia también es un contenedor de DFT entero. Los lóbulos laterales invisibles revelan la fuga que se espera de las sinusoides en otras frecuencias. [a] Por lo tanto, al elegir una función de ventana, generalmente es importante muestrear la DTFT más densamente (como lo hacemos a lo largo de esta sección) y elegir una ventana que suprima los lóbulos laterales a un nivel aceptable.

ventana rectangular

ventana rectangular

La ventana rectangular (a veces conocida como furgón o uniforme o ventana de Dirichlet o engañosamente como "sin ventana" en algunos programas [18] ) es la ventana más simple, equivalente a reemplazar todos menos N valores consecutivos de una secuencia de datos por ceros, haciendo que La forma de onda se enciende y apaga repentinamente:

Otras ventanas están diseñadas para moderar estos cambios repentinos, reducir la pérdida de festón y mejorar el rango dinámico (descrito en § Análisis espectral).

La ventana rectangular es la ventana B -spline de primer orden , así como la ventana de potencia de seno de potencia 0 .

La ventana rectangular proporciona la estimación del error cuadrático medio mínimo de la transformada de Fourier en tiempo discreto , a costa de otras cuestiones discutidas.

B-ventanas spline

Las ventanas B -spline se pueden obtener como convoluciones k veces de la ventana rectangular. Incluyen la propia ventana rectangular ( k  = 1), la § ventana triangular ( k  = 2) y la § ventana Parzen ( k  = 4). [19] Las definiciones alternativas muestran las funciones de base B -spline normalizadas apropiadas en lugar de convolucionar ventanas de tiempo discreto. Una función de base B -spline de orden k es una función polinómica por partes de grado k −1 que se obtiene mediante la autoconvolución k veces de la función rectangular .

ventana triangular

Ventana triangular (con L  =  N  + 1)

Las ventanas triangulares vienen dadas por:

donde L puede ser N , [20] N  + 1, [16] [21] [22] o N  + 2. [23] La primera también se conoce como ventana de Bartlett o ventana de Fejér . Las tres definiciones convergen en general  N .

La ventana triangular es la ventana B -spline de segundo orden . La forma L  =  N puede verse como la convolución de dos ventanas rectangulares de N2 de ancho. La transformada de Fourier del resultado son los valores al cuadrado de la transformada de la ventana rectangular de medio ancho.

ventana parzen

ventana parzen

Al definir LN + 1 , la ventana de Parzen, también conocida como ventana de la Vallée Poussin , [16] es la ventana B -spline de cuarto orden dada por:

ventana galesa

Otras ventanas polinomiales

ventana galesa

La ventana de Welch consta de una única sección parabólica :

[23]

El polinomio cuadrático definitorio alcanza un valor de cero en las muestras justo fuera del intervalo de la ventana.

ventana sinusoidal

ventana sinusoidal

La función correspondiente es un coseno sin el desplazamiento de fase π /2. Por eso, la ventana sinusoidal [24] a veces también se denomina ventana coseno . [16] Como representa medio ciclo de una función sinusoidal, también se la conoce variablemente como ventana de media seno [25] o ventana de media coseno . [26]

La autocorrelación de una ventana sinusoidal produce una función conocida como ventana de Bohman. [27]

Ventanas de potencia de seno/coseno

Estas funciones de ventana tienen la forma: [28]

La ventana rectangular ( α  = 0 ), la ventana sinusoidal ( α  = 1 ) y la ventana de Hann ( α  = 2 ) son miembros de esta familia.

Para valores enteros pares de α, estas funciones también se pueden expresar en forma de suma de cosenos:

Ventanas de suma de coseno

Esta familia también se conoce como ventanas de cosenos generalizados .

En la mayoría de los casos, incluidos los ejemplos siguientes, todos los coeficientes a k  ≥ 0. Estas ventanas tienen solo 2 coeficientes DFT de N puntos K  + 1 distintos de cero .

Ventanas de Hann y Hamming

ventana hann
Ventana de Hamming, a 0  = 0,53836 y a 1  = 0,46164. La ventana de Hamming original tendría un 0  = 0,54 y un 1  = 0,46.

Las ventanas habituales de suma de cosenos para el caso K  = 1 tienen la forma:

que se confunde fácilmente (y a menudo) con su versión de fase cero:

La configuración produce una ventana de Hann:

[29]

lleva el nombre de Julius von Hann , y a veces se lo denomina erróneamente Hanning , presumiblemente debido a sus similitudes lingüísticas y formulaicas con la ventana de Hamming. También se conoce como coseno elevado , porque la versión de fase cero es un lóbulo de una función coseno elevado.

Esta función es miembro de las familias de suma de cosenos y potencias de seno. A diferencia de la ventana de Hamming, los puntos finales de la ventana de Hann simplemente tocan cero. Los lóbulos laterales resultantes disminuyen a unos 18 dB por octava. [30]

Establecer aproximadamente 0,54, o más precisamente 25/46, produce la ventana de Hamming , propuesta por Richard W. Hamming . Esa elección coloca un cruce por cero en la frecuencia 5 π /( N  − 1), que cancela el primer lóbulo lateral de la ventana de Hann, dándole una altura de aproximadamente una quinta parte de la de la ventana de Hann. [16] [31] [32] La ventana de Hamming a menudo se denomina señal de Hamming cuando se utiliza para dar forma al pulso . [33] [34] [35]

La aproximación de los coeficientes a dos decimales reduce sustancialmente el nivel de lóbulos laterales, [16] a una condición casi equivalente a la de la ondulación. [32] En el sentido de la equitriple, los valores óptimos para los coeficientes son a 0  = 0,53836 y a 1  = 0,46164. [32] [36]

ventana negra

ventana de Blackman; α  = 0,16

Las ventanas de Blackman se definen como:

Por convención común, el término incondicional ventana de Blackman se refiere a la "propuesta no muy seria" de Blackman de α  = 0,16 ( a 0  = 0,42, a 1  = 0,5, a 2  = 0,08), que se aproxima mucho al Blackman exacto , [37] con a 0  = 7938/18608 ≈ 0,42659, a 1  = 9240/18608 ≈ 0,49656 y a 2  = 1430/18608 ≈ 0,076849. [38] Estos valores exactos colocan ceros en el tercer y cuarto lóbulo lateral, [16] pero dan como resultado una discontinuidad en los bordes y una caída de 6 dB/oct. Los coeficientes truncados no anulan también los lóbulos laterales, pero tienen una caída mejorada de 18 dB/oct. [16] [39]

Ventana Nuttall, primera derivada continua

Ventana Nuttall, primera derivada continua

La forma continua de la ventana de Nuttall y su primera derivada son continuas en todas partes, como la función de Hann . Es decir, la función llega a 0 en x  = ± N /2, a diferencia de las ventanas de Blackman-Nuttall, Blackman-Harris y Hamming. La ventana de Blackman ( α  = 0,16 ) también es continua con derivada continua en el borde, pero la "ventana de Blackman exacta" no lo es.

Ventana de Blackman-Nuttall

Ventana de Blackman-Nuttall

Ventana de Blackman-Harris

Ventana de Blackman-Harris

Una generalización de la familia Hamming, producida añadiendo más funciones sinc desplazadas, destinadas a minimizar los niveles de lóbulos laterales [40] [41]

ventana superior plana

ventana plana

Una ventana de parte superior plana es una ventana de valor parcialmente negativo que tiene una pérdida de festón mínima en el dominio de la frecuencia. Esa propiedad es deseable para la medición de amplitudes de componentes de frecuencia sinusoidales. [17] [42] Sin embargo, su amplio ancho de banda da como resultado un ancho de banda de alto ruido y una selección de frecuencia más amplia, lo que dependiendo de la aplicación podría ser un inconveniente.

Las ventanas de parte superior plana se pueden diseñar utilizando métodos de diseño de filtro de paso bajo, [42] o pueden ser de la variedad habitual de suma de cosenos:

La variante de Matlab tiene estos coeficientes:

Hay otras variaciones disponibles, como lóbulos laterales que caen a costa de valores más altos cerca del lóbulo principal. [17]

Ventanas Rife-Vincent

Las ventanas de Rife-Vincent [43] habitualmente se escalan según el valor promedio unitario, en lugar del valor máximo unitario. Los valores de los coeficientes siguientes, aplicados a la ecuación 1 , reflejan esa costumbre.

Clase I, Orden 1 ( K = 1): Funcionalmente equivalente a la ventana de Hann.

Clase I, Orden 2 ( K = 2):

La clase I se define minimizando la amplitud del lóbulo lateral de alto orden. Se tabulan los coeficientes para pedidos hasta K=4. [44]

La clase II minimiza el ancho del lóbulo principal para un lóbulo lateral máximo dado.

La clase III es un compromiso cuyo orden K  = 2 se asemeja a la ventana de § Blackman. [44] [45]

Ventanas ajustables

ventana gaussiana

Ventana gaussiana, σ  = 0,4

La transformada de Fourier de una gaussiana también es gaussiana. Dado que el soporte de una función gaussiana se extiende hasta el infinito, debe truncarse en los extremos de la ventana o enventanarse con otra ventana terminada en cero. [46]

Dado que el logaritmo de un gaussiano produce una parábola , esto puede usarse para una interpolación cuadrática casi exacta en la estimación de frecuencia . [47] [46] [48]

La desviación estándar de la función gaussiana es σ  ·  N /2 períodos de muestreo.

Ventana gaussiana confinada, σ t  = 0,1

Ventana gaussiana confinada

La ventana gaussiana confinada produce el ancho de frecuencia cuadrático medio más pequeño posible σ ω para un ancho temporal dado ( N + 1) σ t . [49] Estas ventanas optimizan los productos de ancho de banda de tiempo-frecuencia RMS. Se calculan como los vectores propios mínimos de una matriz dependiente de parámetros. La familia de ventanas gaussianas confinadas contiene la ventana § sinusoidal y la ventana § gaussiana en los casos límite de σ t grande y pequeña , respectivamente.

Ventana gaussiana confinada aproximada, σ t  = 0,1

Ventana gaussiana confinada aproximada

Al definir LN + 1 , una ventana gaussiana confinada de ancho temporal L × σ t se aproxima bien mediante: [49]

donde es una función gaussiana:

La desviación estándar de la ventana aproximada es asintóticamente igual (es decir, valores grandes de N ) a L × σ t para σ t < 0,14 . [49]

Ventana normal generalizada

Una versión más generalizada de la ventana gaussiana es la ventana normal generalizada. [50] Manteniendo la notación de la ventana gaussiana anterior, podemos representar esta ventana como

para cualquier incluso . En , esta es una ventana gaussiana y, a medida que se acerca , se aproxima a una ventana rectangular. La transformada de Fourier de esta ventana no existe en forma cerrada para un general . Sin embargo, demuestra los otros beneficios de tener un ancho de banda ajustable y fluido. Al igual que la ventana de § Tukey, esta ventana ofrece naturalmente una "parte superior plana" para controlar la atenuación de amplitud de una serie temporal (sobre la cual no tenemos control con la ventana gaussiana). En esencia, ofrece un buen compromiso (controlable), en términos de fuga espectral, resolución de frecuencia y atenuación de amplitud, entre la ventana gaussiana y la ventana rectangular. Véase también [51] para un estudio sobre la representación tiempo-frecuencia de esta ventana (o función).

ventana tukey

Ventana de Tukey, α  = 0,5

La ventana de Tukey, también conocida como ventana cónica de coseno , puede considerarse como un lóbulo coseno de ancho /2 (que abarca /2 + 1 observaciones) que está convolucionado con una ventana rectangular de ancho N (1 − α /2 ) .

[52] [B] [C]

En α  = 0 se vuelve rectangular y en α  = 1 se convierte en una ventana de Hann.

Ventana cónica de Planck

Ventana cónica de Planck, ε  = 0,25

La llamada ventana "cono de Planck" es una función de tope que se ha utilizado ampliamente [53] en la teoría de particiones de la unidad en variedades . Es suave (una función) en todas partes, pero es exactamente cero fuera de una región compacta, exactamente uno en un intervalo dentro de esa región, y varía suave y monótonamente entre esos límites. Su uso como función de ventana en el procesamiento de señales se sugirió por primera vez en el contexto de la astronomía de ondas gravitacionales , inspirada en la distribución de Planck . [54] Se define como una función por partes :

La cantidad de disminución está controlada por el parámetro ε , y los valores más pequeños dan transiciones más pronunciadas.

Ventana DPSS o Slepian

La DPSS (secuencia esferoidal prolata discreta) o ventana de Slepian maximiza la concentración de energía en el lóbulo principal , [55] y se utiliza en el análisis espectral multicónico , que promedia el ruido en el espectro y reduce la pérdida de información en los bordes de la ventana.

El lóbulo principal termina en un rango de frecuencia dado por el parámetro α . [56]

Las siguientes ventanas de Kaiser se crean mediante una simple aproximación a las ventanas de DPSS:

ventana káiser

La ventana Kaiser, o Kaiser-Bessel, es una aproximación simple de la ventana DPSS que utiliza funciones de Bessel , descubierta por James Kaiser . [57] [58]

[D] [16] : pág. 73 

¿Dónde está la función de Bessel modificada de orden 0 de primer tipo? El parámetro variable determina el equilibrio entre el ancho del lóbulo principal y los niveles de los lóbulos laterales del patrón de fuga espectral. El ancho del lóbulo principal, entre los nulos, viene dado en unidades de contenedores DFT, [65] y un valor típico es 3.

Ventana de Dolph-Chebyshev

Ventana de Dolph-Chebyshev, α  = 5

Minimiza la norma de Chebyshev de los lóbulos laterales para un ancho de lóbulo principal determinado. [66]

La función de ventana de Dolph-Chebyshev de fase cero generalmente se define en términos de su transformada de Fourier discreta de valor real : [67]

T n ( x ) es el n -ésimo polinomio de Chebyshev del primer tipo evaluado en x , que se puede calcular usando

y

es la única solución real positiva para , donde el parámetro α establece la norma de Chebyshev de los lóbulos laterales en −20 α  decibelios. [66]

La función de ventana se puede calcular a partir de W 0 ( k ) mediante una transformada de Fourier discreta inversa (DFT): [66]

La versión retrasada de la ventana se puede obtener mediante:

que para valores pares de N debe calcularse de la siguiente manera:

que es una DFT inversa de

Variaciones:

ventana ultraesférica

El parámetro μ de la ventana ultraesférica determina si las amplitudes de los lóbulos laterales de su transformada de Fourier disminuyen, están niveladas o (como se muestra aquí) aumentan con la frecuencia.

La ventana ultraesférica fue introducida en 1984 por Roy Streit [68] y tiene aplicación en el diseño de conjuntos de antenas, [69] diseño de filtros no recursivos [68] y análisis de espectro. [70]

Like other adjustable windows, the Ultraspherical window has parameters that can be used to control its Fourier transform main-lobe width and relative side-lobe amplitude. Uncommon to other windows, it has an additional parameter which can be used to set the rate at which side-lobes decrease (or increase) in amplitude.[70][71][72]

The window can be expressed in the time-domain as follows:[70]

where is the Ultraspherical polynomial of degree N, and and control the side-lobe patterns.[70]

Certain specific values of yield other well-known windows: and give the Dolph–Chebyshev and Saramäki windows respectively.[68] See here for illustration of Ultraspherical windows with varied parametrization.

Exponential or Poisson window

Exponential window, τ = N/2
Exponential window, τ = (N/2)/(60/8.69)

The Poisson window, or more generically the exponential window increases exponentially towards the center of the window and decreases exponentially in the second half. Since the exponential function never reaches zero, the values of the window at its limits are non-zero (it can be seen as the multiplication of an exponential function by a rectangular window [73]). It is defined by

where τ is the time constant of the function. The exponential function decays as e ≃ 2.71828 or approximately 8.69 dB per time constant.[74]This means that for a targeted decay of D dB over half of the window length, the time constant τ is given by

Hybrid windows

Window functions have also been constructed as multiplicative or additive combinations of other windows.

Bartlett–Hann window

Bartlett–Hann window

Planck–Bessel window

Planck–Bessel window, ε = 0.1, α = 4.45

A § Planck-taper window multiplied by a Kaiser window which is defined in terms of a modified Bessel function. This hybrid window function was introduced to decrease the peak side-lobe level of the Planck-taper window while still exploiting its good asymptotic decay.[75] It has two tunable parameters, ε from the Planck-taper and α from the Kaiser window, so it can be adjusted to fit the requirements of a given signal.

Hann–Poisson window

Hann–Poisson window, α = 2

A Hann window multiplied by a Poisson window. For it has no side-lobes, as its Fourier transform drops off forever away from the main lobe without local minima. It can thus be used in hill climbing algorithms like Newton's method.[76] The Hann–Poisson window is defined by:

where α is a parameter that controls the slope of the exponential.

Other windows

GAP window (GAP optimized Nuttall window)

Generalized adaptive polynomial (GAP) window

The GAP window is a family of adjustable window functions that are based on a symmetrical polynomial expansion of order . It is continuous with continuous derivative everywhere. With the appropriate set of expansion coefficients and expansion order, the GAP window can mimic all the known window functions, reproducing accurately their spectral properties.

[77]

where is the standard deviation of the sequence.

Additionally, starting with a set of expansion coefficients that mimics a certain known window function, the GAP window can be optimized by minimization procedures to get a new set of coefficients that improve one or more spectral properties, such as the main lobe width, side lobe attenuation, and side lobe falloff rate.[78] Therefore, a GAP window function can be developed with designed spectral properties depending on the specific application.

Sinc or Lanczos window

Lanczos window

Asymmetric window functions

The form, according to the convention above, is symmetric around . However, there are window functions that are asymmetric, such as the Gamma distribution used in FIR implementations of Gammatone filters. These asymmetries are used to reduce the delay when using large window sizes, or to emphasize the initial transient of a decaying pulse.[citation needed]

Any bounded function with compact support, including asymmetric ones, can be readily used as a window function. Additionally, there are ways to transform symmetric windows into asymmetric windows by transforming the time coordinate, such as with the below formula

where the window weights more highly the earliest samples when , and conversely weights more highly the latest samples when .[79]

See also

Notes

  1. ^ Some authors limit their attention to this important subset and to even values of N.[16][17] But the window coefficient formulas are still the ones presented here.
  2. ^ This formula can be confirmed by simplifying the cosine function at MATLAB tukeywin and substituting r=α and x=n/N.
  3. ^ Harris 1978 (p 67, eq 38) appears to have two errors: (1) The subtraction operator in the numerator of the cosine function should be addition. (2) The denominator contains a spurious factor of 2. Also, Fig 30 corresponds to α=0.25 using the Wikipedia formula, but to 0.75 using the Harris formula. Fig 32 is similarly mislabeled.
  4. ^ The Kaiser window is often parametrized by β, where β = πα.[59][60][61][62][56][63][7]: p. 474  The alternative use of just α facilitates comparisons to the DPSS windows.[64]

Page citations

  1. ^ Harris 1978, p 57, fig 10.

References

  1. ^ Weisstein, Eric W. (2003). CRC Concise Encyclopedia of Mathematics. CRC Press. ISBN 978-1-58488-347-0.
  2. ^ Roads, Curtis (2002). Microsound. MIT Press. ISBN 978-0-262-18215-7.
  3. ^ Cattani, Carlo; Rushchitsky, Jeremiah (2007). Wavelet and Wave Analysis As Applied to Materials With Micro Or Nanostructure. World Scientific. ISBN 978-981-270-784-0.
  4. ^ "Overlap-Add (OLA) STFT Processing | Spectral Audio Signal Processing". www.dsprelated.com. Retrieved 2016-08-07. The window is applied twice: once before the FFT (the "analysis window") and secondly after the inverse FFT prior to reconstruction by overlap-add (the so-called "synthesis window"). ... More generally, any positive COLA window can be split into an analysis and synthesis window pair by taking its square root.
  5. ^ Ajala, R.; Persaud, P. (2022). "Ground-Motion Evaluation of Hybrid Seismic Velocity Models". The Seismic Record. 2 (3): 186–196. Bibcode:2022SeisR...2..186A. doi:10.1785/0320220022. S2CID 251504921.
  6. ^ Ajala, R.; Persaud, P. (2021). "Effect of Merging Multiscale Models on Seismic Wavefield Predictions Near the Southern San Andreas Fault". Journal of Geophysical Research: Solid Earth. 126 (10). Bibcode:2021JGRB..12621915A. doi:10.1029/2021JB021915. ISSN 2169-9313. S2CID 239654900.
  7. ^ a b Oppenheim, Alan V.; Schafer, Ronald W.; Buck, John R. (1999). "7.2". Discrete-time signal processing (2nd ed.). Upper Saddle River, N.J.: Prentice Hall. pp. 465–478. ISBN 0-13-754920-2.
  8. ^ "FIR Filters by Windowing – The Lab Book Pages". www.labbookpages.co.uk. Retrieved 2016-04-13.
  9. ^ "Mastering Windows" (PDF). www.cg.tuwien.ac.at. Retrieved 2020-02-12.
  10. ^ "The Fundamentals of Signal Analysis Application Note 243" (PDF). hpmemoryproject.org. Retrieved 10 April 2018.
  11. ^ R. Hovden, Y. Jiang, H. Xin, L.F. Kourkoutis (2015). "Periodic Artifact Reduction in Fourier Transforms of Full Field Atomic Resolution Images". Microscopy and Microanalysis. 21 (2): 436–441. arXiv:2210.09024. Bibcode:2015MiMic..21..436H. doi:10.1017/S1431927614014639. PMID 25597865. S2CID 22435248.{{cite journal}}: CS1 maint: multiple names: authors list (link)
  12. ^ Bernstein, Matt A.; King, Kevin Franklin; Zhou, Xiaohong Joe (2004). Handbook of MRI Pulse Sequences. London: Elsevier Academic Press. pp. 495–499. ISBN 0-12-092861-2.
  13. ^ Awad, A.I.; Baba, K. (2011). "An Application for Singular Point Location in Fingerprint Classification". Digital Information Processing and Communications. Communications in Computer and Information Science. Vol. 188. p. 262. doi:10.1007/978-3-642-22389-1_24. ISBN 978-3-642-22388-4.
  14. ^ "Zero Phase Filters". ccrma.stanford.edu. Retrieved 2020-02-12.
  15. ^ Rorabaugh, C.Britton (October 1998). DSP Primer. Primer series. McGraw-Hill Professional. p. 196. ISBN 978-0-07-054004-0.
  16. ^ a b c d e f g h i j Harris, Fredric J. (Jan 1978). "On the use of Windows for Harmonic Analysis with the Discrete Fourier Transform" (PDF). Proceedings of the IEEE. 66 (1): 51–83. Bibcode:1978IEEEP..66...51H. CiteSeerX 10.1.1.649.9880. doi:10.1109/PROC.1978.10837. S2CID 426548. The fundamental 1978 paper on FFT windows by Harris, which specified many windows and introduced key metrics used to compare them.
  17. ^ a b c Heinzel, G.; Rüdiger, A.; Schilling, R. (2002). Spectrum and spectral density estimation by the Discrete Fourier transform (DFT), including a comprehensive list of window functions and some new flat-top windows (Technical report). Max Planck Institute (MPI) für Gravitationsphysik / Laser Interferometry & Gravitational Wave Astronomy. 395068.0. Retrieved 2013-02-10. Also available at https://pure.mpg.de/rest/items/item_152164_1/component/file_152163/content
  18. ^ "Understanding FFTs and Windowing" (PDF). National Instruments. Archived (PDF) from the original on 2024-01-05. Retrieved 2024-02-13.
  19. ^ Toraichi, K.; Kamada, M.; Itahashi, S.; Mori, R. (1989). "Window functions represented by B-spline functions". IEEE Transactions on Acoustics, Speech, and Signal Processing. 37: 145–147. doi:10.1109/29.17517.
  20. ^ "Bartlett Window". ccrma.stanford.edu. Retrieved 2016-04-13.
  21. ^ Tukey, J.W. (1967). "An introduction to the calculations of numerical spectrum analysis". Spectral Analysis of Time Series: 25–46.
  22. ^ "Triangular window – MATLAB triang". www.mathworks.com. Retrieved 2016-04-13.
  23. ^ a b Welch, P. (1967). "The use of fast Fourier transform for the estimation of power spectra: A method based on time averaging over short, modified periodograms". IEEE Transactions on Audio and Electroacoustics. 15 (2): 70–73. Bibcode:1967ITAE...15...70W. doi:10.1109/TAU.1967.1161901. S2CID 13900622.
  24. ^ Bosi, Marina; Goldberg, Richard E. (2003). "Time to Frequency Mapping Part II: The MDCT". Introduction to Digital Audio Coding and Standards. The Springer International Series in Engineering and Computer Science. Vol. 721. Boston, MA: Springer US. p. 106. doi:10.1007/978-1-4615-0327-9. ISBN 978-1-4615-0327-9.
  25. ^ Kido, Ken'iti; Suzuki, Hideo; Ono, Takahiko; Fukushima, Manabu (1998). "Deformation of impulse response estimates by time window in cross spectral technique". Journal of the Acoustical Society of Japan (E). 19 (5): 349–361. doi:10.1250/ast.19.349.
  26. ^ Landisman, M.; Dziewonski, A.; Satô, Y. (1969-05-01). "Recent Improvements in the Analysis of Surface Wave Observations". Geophysical Journal International. 17 (4): 369–403. Bibcode:1969GeoJ...17..369L. doi:10.1111/j.1365-246X.1969.tb00246.x.
  27. ^ "Bohman window – R2019B". www.mathworks.com. Retrieved 2020-02-12.
  28. ^ "Power-of-Cosine Window Family". ccrma.stanford.edu. Retrieved 10 April 2018.
  29. ^ "Hann (Hanning) window - MATLAB hann". www.mathworks.com. Retrieved 2020-02-12.
  30. ^ "Hann or Hanning or Raised Cosine". ccrma.stanford.edu. Retrieved 2016-04-13.
  31. ^ Enochson, Loren D.; Otnes, Robert K. (1968). Programming and Analysis for Digital Time Series Data. U.S. Dept. of Defense, Shock and Vibration Info. Center. p. 142.
  32. ^ a b c "Hamming Window". ccrma.stanford.edu. Retrieved 2016-04-13.
  33. ^ "A digital quadrature amplitude modulation (QAM) Radio: Building a better radio" (PDF). users.wpi.edu. p. 28. Retrieved 2020-02-12.
  34. ^ "Bits to Symbols to Signals and back again" (PDF). users.wpi.edu. p. 7. Retrieved 2020-02-12.
  35. ^ Johnson, C.Richard Jr; Sethares, William A.; Klein, Andrew G. (2011-08-18). "11". Software Receiver Design. Cambridge University Press. ISBN 978-1-139-50145-3. Also https://cnx.org/contents/[email protected]:6R_ztzDY@4/Pulse-Shaping-and-Receive-Filtering
  36. ^ Nuttall, Albert H. (Feb 1981). "Some Windows with Very Good Sidelobe Behavior". IEEE Transactions on Acoustics, Speech, and Signal Processing. 29 (1): 84–91. doi:10.1109/TASSP.1981.1163506. Extends Harris' paper, covering all the window functions known at the time, along with key metric comparisons.
  37. ^ Weisstein, Eric W. "Blackman Function". mathworld.wolfram.com. Retrieved 2016-04-13.
  38. ^ "Characteristics of Different Smoothing Windows - NI LabVIEW 8.6 Help". zone.ni.com. Retrieved 2020-02-13.
  39. ^ Blackman, R.B.; Tukey, J.W. (1959-01-01). The Measurement of Power Spectra from the Point of View of Communications Engineering. Dover Publications. p. 99. ISBN 978-0-486-60507-4.
  40. ^ "Blackman-Harris Window Family". ccrma.stanford.edu. Retrieved 2016-04-13.
  41. ^ "Three-Term Blackman-Harris Window". ccrma.stanford.edu. Retrieved 2016-04-13.
  42. ^ a b Smith, Steven W. (2011). The Scientist and Engineer's Guide to Digital Signal Processing. San Diego, California, USA: California Technical Publishing. Retrieved 2013-02-14.
  43. ^ Rife, David C.; Vincent, G.A. (1970), "Use of the discrete Fourier transform in the measurement of frequencies and levels of tones", Bell Syst. Tech. J., 49 (2): 197–228, doi:10.1002/j.1538-7305.1970.tb01766.x
  44. ^ a b Andria, Gregorio; Savino, Mario; Trotta, Amerigo (1989), "Windows and interpolation algorithms to improve electrical measurement accuracy", IEEE Transactions on Instrumentation and Measurement, 38 (4): 856–863, Bibcode:1989ITIM...38..856A, doi:10.1109/19.31004
  45. ^ Schoukens, Joannes; Pintelon, Rik; Van Hamme, Hugo (1992), "The interpolated fast Fourier transform: a comparative study", IEEE Transactions on Instrumentation and Measurement, 41 (2): 226–232, Bibcode:1992ITIM...41..226S, doi:10.1109/19.137352
  46. ^ a b "Matlab for the Gaussian Window". ccrma.stanford.edu. Retrieved 2016-04-13. Note that, on a dB scale, Gaussians are quadratic. This means that parabolic interpolation of a sampled Gaussian transform is exact. ... quadratic interpolation of spectral peaks may be more accurate on a log-magnitude scale (e.g., dB) than on a linear magnitude scale
  47. ^ "Gaussian Window and Transform". ccrma.stanford.edu. Retrieved 2016-04-13.
  48. ^ "Quadratic Interpolation of Spectral Peaks". ccrma.stanford.edu. Retrieved 2016-04-13.
  49. ^ a b c Starosielec, S.; Hägele, D. (2014). "Discrete-time windows with minimal RMS bandwidth for given RMS temporal width". Signal Processing. 102: 240–246. Bibcode:2014SigPr.102..240S. doi:10.1016/j.sigpro.2014.03.033.
  50. ^ Chakraborty, Debejyo; Kovvali, Narayan (2013). "Generalized normal window for digital signal processing". 2013 IEEE International Conference on Acoustics, Speech and Signal Processing. pp. 6083–6087. doi:10.1109/ICASSP.2013.6638833. ISBN 978-1-4799-0356-6. S2CID 11779529.
  51. ^ Diethorn, E.J. (1994). "The generalized exponential time-frequency distribution". IEEE Transactions on Signal Processing. 42 (5): 1028–1037. Bibcode:1994ITSP...42.1028D. doi:10.1109/78.295214.
  52. ^ Bloomfield, P. (2000). Fourier Analysis of Time Series: An Introduction. New York: Wiley-Interscience.
  53. ^ Tu, Loring W. (2008). "Bump Functions and Partitions of Unity". An Introduction to Manifolds. Universitext. New York: Springer. pp. 127–134. doi:10.1007/978-0-387-48101-2_13. ISBN 978-0-387-48098-5.
  54. ^ McKechan, D.J.A.; Robinson, C.; Sathyaprakash, B.S. (21 April 2010). "A tapering window for time-domain templates and simulated signals in the detection of gravitational waves from coalescing compact binaries". Classical and Quantum Gravity. 27 (8): 084020. arXiv:1003.2939. Bibcode:2010CQGra..27h4020M. doi:10.1088/0264-9381/27/8/084020. S2CID 21488253.
  55. ^ "Slepian or DPSS Window". ccrma.stanford.edu. Retrieved 2016-04-13.
  56. ^ a b Smith, J.O. (2011). "Kaiser and DPSS Windows Compared". ccrma.stanford.edu. Retrieved 2016-04-13.
  57. ^ Kaiser, James F.; Kuo, Franklin F. (1966). System Analysis by Digital Computer. John Wiley and Sons. pp. 232–235. This family of window functions was "discovered" by Kaiser in 1962 following a discussion with B. F. Logan of the Bell Telephone Laboratories. ... Another valuable property of this family ... is that they also approximate closely the prolate spheroidal wave functions of order zero.
  58. ^ Kaiser, James F. (Nov 1964). "A family of window functions having nearly ideal properties". Unpublished Memorandum.
  59. ^ Rabiner, Lawrence R.; Gold, Bernard (1975). "3.11". Theory and application of digital signal processing. Englewood Cliffs, N.J.: Prentice-Hall. p. 94. ISBN 0-13-914101-4.
  60. ^ Crochiere, R.E.; Rabiner, L.R. (1983). "4.3.1". Multirate Digital Signal Processing. Englewood Cliffs, NJ: Prentice-Hall. p. 144. ISBN 0-13-605162-6.
  61. ^ Lin, Yuan-Pei; Vaidyanathan, P.P. (June 1998). "A Kaiser Window Approach for the Design of Prototype Filters of Cosine Modulated Filterbanks" (PDF). IEEE Signal Processing Letters. 5 (6): 132–134. Bibcode:1998ISPL....5..132L. doi:10.1109/97.681427. S2CID 18159105. Retrieved 2017-03-16.
  62. ^ Smith, J.O. (2011). "Kaiser Window". ccrma.stanford.edu. Retrieved 2019-03-20. Sometimes the Kaiser window is parametrized by α, where β = πα.
  63. ^ "Kaiser Window, R2020a". www.mathworks.com. Mathworks. Retrieved 9 April 2020.
  64. ^ "Kaiser Window". www.dsprelated.com. Retrieved 2020-04-08. The following Matlab comparison of the DPSS and Kaiser windows illustrates the interpretation of α as the bin number of the edge of the critically sampled window main lobe.
  65. ^ Kaiser, James F.; Schafer, Ronald W. (1980). "On the use of the I0-sinh window for spectrum analysis". IEEE Transactions on Acoustics, Speech, and Signal Processing. 28: 105–107. doi:10.1109/TASSP.1980.1163349.
  66. ^ a b c "Dolph-Chebyshev Window". ccrma.stanford.edu. Retrieved 2016-04-13.
  67. ^ "Dolph-Chebyshev Window Definition". ccrma.stanford.edu. Retrieved 2019-03-05.
  68. ^ a b c Kabal, Peter (2009). "Time Windows for Linear Prediction of Speech" (PDF). Technical Report, Dept. Elec. & Comp. Eng., McGill University (2a): 31. Retrieved 2 February 2014.
  69. ^ Streit, Roy (1984). "A two-parameter family of weights for nonrecursive digital filters and antennas". Transactions of ASSP. 32: 108–118. doi:10.1109/tassp.1984.1164275.
  70. ^ a b c d Deczky, Andrew (2001). "Unispherical Windows". ISCAS 2001. The 2001 IEEE International Symposium on Circuits and Systems (Cat. No. 01CH37196). Vol. 2. pp. 85–88. doi:10.1109/iscas.2001.921012. ISBN 978-0-7803-6685-5. S2CID 38275201.
  71. ^ Bergen, S.W.A.; Antoniou, A. (2004). "Design of Ultraspherical Window Functions with Prescribed Spectral Characteristics". EURASIP Journal on Applied Signal Processing. 2004 (13): 2053–2065. Bibcode:2004EJASP2004...63B. doi:10.1155/S1110865704403114.
  72. ^ Bergen, Stuart W. A. (2005). "Design of the Ultraspherical Window Function and Its Applications" (PDF). Dissertation, University of Viktoria.
  73. ^ Smith, Julius O. III (2011-04-23). "Poisson Window". ccrma.stanford.edu. Retrieved 2020-02-12.
  74. ^ Gade, Svend; Herlufsen, Henrik (1987). "Technical Review No 3-1987: Windows to FFT analysis (Part I)" (PDF). Brüel & Kjær. Retrieved 2011-11-22.
  75. ^ Berry, C.P.L.; Gair, J.R. (12 December 2012). "Observing the Galaxy's massive black hole with gravitational wave bursts". Monthly Notices of the Royal Astronomical Society. 429 (1): 589–612. arXiv:1210.2778. Bibcode:2013MNRAS.429..589B. doi:10.1093/mnras/sts360. S2CID 118944979.
  76. ^ "Hann-Poisson Window". ccrma.stanford.edu. Retrieved 2016-04-13.
  77. ^ Wesley Beccaro (2020-10-31), "Generalized Adaptive Polynomial Window Function", mathworks.com, retrieved 2020-11-02
  78. ^ "Generalized Adaptive Polynomial Window Function". www.mathworks.com. Retrieved 2020-12-12.
  79. ^ Luo, Jiufel; Xie, Zhijiang; Li, Xinyi (2015-03-02). "Asymmetric Windows and Their Application in Frequency Estimation". Chongqing University. 9 (Algorithms & Computational Technology): 389–412. doi:10.1260/1748-3018.9.4.389. S2CID 124464194.

Further reading

External links