Articulo de referencia

función de Hann

Función de Hann (izquierda) y su respuesta en frecuencia (derecha). La función de Hann recibe su nombre del meteorólogo austriaco Julius von Hann . Es una función de ventana que...

Función de Hann (izquierda) y su respuesta en frecuencia (derecha).

La función de Hann recibe su nombre del meteorólogo austriaco Julius von Hann . Es una función de ventana que se utiliza para realizar el suavizado de Hann o hanning . [ 1 ] [ 2 ] La función, con longitudL{\displaystyle L}y amplitud1/L,{\displaystyle 1/L,}está dado por:

w0(incógnita){1L(12+12porque(2πincógnitaL))=1Lporque2(πincógnitaL),|incógnita|L/20,|incógnita|>L/2}.{\displaystyle w_{0}(x)\triangleq \left\{{\begin{array}{ccl}{\tfrac {1}{L}}\left({\tfrac {1}{2}}+{\tfrac {1}{2}}\cos \left({\frac {2\pi x}{L}}\right)\right)={\tfrac {1}{L}}\cos ^{2}\left({\frac {\pi x}{L}}\right),\quad &\left|x\right|\leq L/2\\0,\quad &\left|x\right|>L/2\end{array}}\right\}.} [ a ]

Para el procesamiento de señales digitales , la función se muestrea simétricamente (con espaciadoL/norte{\displaystyle L/N}y amplitud1{\displaystyle 1}) :

w[norte]=Lw0(Lnorte(nortenorte/2))=12[1porque(2πnortenorte)]=pecado2(πnortenorte)},0nortenorte,{\displaystyle \left.{\begin{aligned}w[n]=L\cdot w_{0}\left({\tfrac {L}{N}}(nN/2)\right)&={\tfrac {1}{2}}\left[1-\cos \left({\tfrac {2\pi n}{N}}\right)\right]\\&=\sin ^{2}\left({\tfrac {\pi n}{N}}\right)\end{aligned}}\right\},\quad 0\leq n\leq N,}

que es una secuencia denorte+1{\displaystyle N+1}muestras ynorte{\displaystyle N}Puede ser par o impar. También se conoce como ventana de coseno elevado , filtro de Hann , ventana de von Hann , ventana de Hanning , etc. [ 2 ] [ 3 ] [ 4 ]

transformada de Fourier

Arriba: Ventana de Hann par de 16 muestras de la transformada discreta de Fourier (DFT ). Abajo: Su transformada discreta de Fourier (DTFT) y los 3 valores distintos de cero de su transformada discreta de Fourier (DFT).

La transformada de Fourier dew0(incógnita){\displaystyle w_{0}(x)}está dado por:

W0(F)=12desde(LF)(1L2F2)=pecado(πLF)2πLF(1L2F2){\displaystyle W_{0}(f)={\frac {1}{2}}{\frac {\operatorname {sinc} (Lf)}{(1-L^{2}f^{2})}}={\frac {\sin(\pi Lf)}{2\pi Lf(1-L^{2}f^{2})}}} [ b ]
Derivación

Utilizando la fórmula de Euler para expandir el término coseno enw0(incógnita),{\displaystyle w_{0}(x),}podemos escribir:

w0(incógnita)=1L(12recto(incógnita/L)+14mii2πincógnita/Lrecto(incógnita/L)+14mii2πincógnita/Lrecto(incógnita/L)),{\displaystyle w_{0}(x)={\tfrac {1}{L}}\left({\tfrac {1}{2}}\operatorname {rect} (x/L)+{\tfrac {1}{4}}e^{i2\pi x/L}\operatorname {rect} (x/L)+{\tfrac {1}{4}}e^{-i2\pi x/L}\operatorname {rect} (x/L)\right),}

que es una combinación lineal de ventanas rectangulares moduladas :

1Lrecto(incógnita/L)transformada de Fourierdesde(LF)pecado(πLF)πLF.{\displaystyle {\tfrac {1}{L}}\operatorname {rect} (x/L)\quad {\stackrel {\text{Fourier transform}}{\longleftrightarrow }}\quad \operatorname {sinc} (Lf)\triangleq {\frac {\sin(\pi Lf)}{\pi Lf}}.}

Transformando cada término:

W0(F)=12desde(LF)+14desde(L(F1/L))+14desde(L(F+1/L))=12pecado(πLF)πLF+14pecado(π(LF1))π(LF1)+14pecado(π(LF+1))π(LF+1)=12π(pecado(πLF)LF12pecado(πLF)LF112pecado(πLF)LF+1)=pecado(πLF)2π(1LF+1211LF1211+LF)=pecado(πLF)2π1LF(1LF)(1+LF)=12desde(LF)(1L2F2).{\displaystyle {\begin{aligned}W_{0}(f)&={\tfrac {1}{2}}\operatorname {sinc} (Lf)+{\tfrac {1}{4}}\operatorname {sinc} (L(f-1/L))+{\tfrac {1}{4}}\operatorname {sinc} (L(f+1/L))\\&={\tfrac {1}{2}}{\frac {\sin(\pi Lf)}{\pi Lf}}+{\tfrac {1}{4}}{\frac {\sin(\pi (Lf-1))}{\pi (Lf-1)}}+{\tfrac {1}{4}}{\frac {\sin(\pi (Lf+1))}{\pi (Lf+1)}}\\&={\frac {1}{2\pi }}\left({\frac {\sin(\pi Lf)}{Lf}}-{\tfrac {1}{2}}{\frac {\sin(\pi Lf)}{Lf-1}}-{\tfrac {1}{2}}{\frac {\sin(\pi Lf)}{Lf+1}}\right)\\&={\frac {\sin(\pi Lf)}{2\pi }}\left({\frac {1}{Lf}}+{\tfrac {1}{2}}{\frac {1}{1-Lf}}-{\tfrac {1}{2}}{\frac {1}{1+Lf}}\right)\\&={\frac {\sin(\pi Lf)}{2\pi }}\cdot {\frac {1}{Lf(1-Lf)(1+Lf)}}={\frac {1}{2}}{\frac {\operatorname {sinc} (Lf)}{(1-L^{2}f^{2})}}.\end{aligned}}}

Transformaciones discretas

La transformada discreta de Fourier (DTFT) de lanorte+1{\displaystyle N+1}La secuencia de longitud y desplazamiento temporal se define mediante una serie de Fourier , que también tiene un equivalente de 3 términos que se deriva de forma similar a la derivación de la transformada de Fourier:

F{w[norte]}norte=0nortew[norte]mii2πFnorte=miiπFnorte[12pecado(π(norte+1)F)pecado(πF)+14pecado(π(norte+1)(F1norte))pecado(π(F1norte))+14pecado(π(norte+1)(F+1norte))pecado(π(F+1norte))].{\displaystyle {\begin{aligned}{\mathcal {F}}\{w[n]\}&\triangleq \sum _{n=0}^{N}w[n]\cdot e^{-i2\pi fn}\\&=e^{-i\pi fN}\left[{\tfrac {1}{2}}{\frac {\sin(\pi (N+1)f)}{\sin(\pi f)}}+{\tfrac {1}{4}}{\frac {\sin(\pi (N+1)(f-{\tfrac {1}{N}}))}{\sin(\pi (f-{\tfrac {1}{N}}))}}+{\tfrac {1}{4}}{\frac {\sin(\pi (N+1)(f+{\tfrac {1}{N}}))}{\sin(\pi (f+{\tfrac {1}{N}}))}}\right].\end{aligned}}}

La secuencia truncada{w[norte], 0nortenorte1}{\displaystyle \{w[n],\ 0\leq n\leq N-1\}}es una ventana de Hann par (también conocida como periódica ) de DFT . Dado que la muestra truncada tiene valor cero, es evidente, a partir de la definición de la serie de Fourier, que las DTFT son equivalentes. Sin embargo, el enfoque seguido anteriormente da como resultado una expresión de 3 términos de aspecto significativamente diferente, pero equivalente:

F{w[norte]}=miiπF(norte1)[12pecado(πnorteF)pecado(πF)+14miiπ/nortepecado(πnorte(F1norte))pecado(π(F1norte))+14miiπ/nortepecado(πnorte(F+1norte))pecado(π(F+1norte))].{\displaystyle {\mathcal {F}}\{w[n]\}=e^{-i\pi f(N-1)}\left[{\tfrac {1}{2}}{\frac {\sin(\pi Nf)}{\sin(\pi f)}}+{\tfrac {1}{4}}e^{-i\pi /N}{\frac {\sin(\pi N(f-{\tfrac {1}{N}}))}{\sin(\pi (f-{\tfrac {1}{N}}))}}+{\tfrac {1}{4}}e^{i\pi /N}{\frac {\sin(\pi N(f+{\tfrac {1}{N}}))}{\sin(\pi (f+{\tfrac {1}{N}}))}}\right].}

Una DFT de longitud N de la función ventana muestrea la DTFT en frecuenciasF=k/norte,{\displaystyle f=k/N,}para valores enteros dek.{\displaystyle k.} De la expresión inmediatamente anterior, es fácil ver que solo 3 de los N coeficientes DFT son distintos de cero. Y de la otra expresión, es evidente que todos son de valor real. Estas propiedades son atractivas para aplicaciones en tiempo real que requieren transformadas con y sin ventana (con ventana rectangular), porque las transformadas con ventana se pueden derivar eficientemente de las transformadas sin ventana mediante convolución . [ 5 ] [ c ] [ d ]

Nombre

La función recibe su nombre en honor a von Hann, quien utilizó la técnica de suavizado de promedio ponderado de tres términos en datos meteorológicos. [ 6 ] [ 2 ] Sin embargo, también se utiliza convencionalmente el término función de Hanning , [ 7 ] derivado del artículo en el que el término hanning una señal se utilizó para referirse a la aplicación de la ventana de Hann a la misma. [ 4 ] [ 8 ] Es distinta de la función de Hamming , de nombre similar , que recibe su nombre de Richard Hamming .

Véase también

Citas de página

  1. Nuttall 1981 , pág. 84 (3)
  2. Nuttall 1981 , pág. 86 (17)
  3. Nuttall 1981 , pág. 85
  4. Harris 1978 , pág. 62

Referencias

  1. Essenwanger, OM (Oskar M.) (1986). Elementos de análisis estadístico . Elsevier. ISBN 0444424261OCLC 152410575 
  2. 1 2 3 Kahlig, Peter (1993), "Algunos aspectos de la contribución de Julius von Hann a la climatología moderna" , en McBean, GA; Hantel, M. (eds.), Interacciones entre subsistemas climáticos globales: El legado de Hann , Geophysical Monograph Series, vol. 75, American Geophysical Union, pp. 1–7 , doi : 10.1029/gm075p0001 , ISBN   9780875904665, recuperado el 01/07/2019 , Hann parece ser el inventor de un cierto procedimiento de suavizado de datos, ahora llamado "hanning" ... o "suavizado de Hann" ... Esencialmente, es una media móvil de tres términos (media móvil) con pesos desiguales (1/4, 1/2, 1/4).
  3. Smith, Julius O. (Julius Orion) (2011). Procesamiento espectral de señales de audio . Universidad de Stanford. Centro de Investigación Informática en Música y Acústica, Universidad de Stanford. Departamento de Música. [Stanford, California?]: W3K. ISBN 9780974560731OCLC 776892709 
  4. 1 2 Blackman, RB ; Tukey, JW (1958). "La medición de espectros de potencia desde el punto de vista de la ingeniería de comunicaciones — Parte I". The Bell System Technical Journal . 37 (1): 273. doi : 10.1002/j.1538-7305.1958.tb03874.x . ISSN 0005-8580 . 
  5. ↑ Patente estadounidense 6898235 , Carlin, Joe; Collins, Terry y Hays, Peter et al., "Dispositivo de interceptación de comunicaciones de banda ancha y localización de direcciones mediante hipercanalización", publicada el 10 de diciembre de 1999, emitida el 24 de mayo de 2005 , también disponible en https://patentimages.storage.googleapis.com/4d/39/2a/cec2ae6f33c1e7/US6898235.pdf 
  6. von Hann, Julius (1903). Manual de Climatología . Macmillan. pág. 199. Las cifras bajo b se determinan teniendo en cuenta los paralelos a 5° de distancia a cada lado. Así, por ejemplo, para una latitud de 60° tenemos ½[60 + (65 + 55)÷2].     
  7. Harris, Fredric J. (enero de 1978). "Sobre el uso de ventanas para el análisis armónico con la transformada discreta de Fourier" (PDF) . Actas del IEEE . 66 (1): 51–83 . CiteSeerX 10.1.1.649.9880 . doi : 10.1109/PROC.1978.10837 . El nombre correcto de esta ventana es 'Hann'. El término 'Hanning' se utiliza en este informe para reflejar el uso convencional. El término derivado 'Hann'd' también se utiliza ampliamente. 
  8. Blackman, RB (Ralph Beebe) ; Tukey, John W. (John Wilder) (1959). La medición de espectros de potencia desde el punto de vista de la ingeniería de comunicaciones . Nueva York : Dover Publications. págs . 98. LCCN 59-10185 .   {{cite book}}: CS1 mantenimiento: ubicación del editor ( enlace )
  1. Nuttall, Albert H. (febrero de 1981). "Algunas ventanas con un comportamiento de lóbulos laterales muy bueno" . IEEE Transactions on Acoustics, Speech, and Signal Processing . 29 (1): 84– 91. doi : 10.1109/TASSP.1981.1163506 .