Articulo de referencia

Algoritmo adaptativo-aditivo

En los estudios de óptica de Fourier , síntesis de sonido , interferometría estelar , pinzas ópticas y elementos ópticos difractivos (DOE), a menudo es importante conocer la fas...

En los estudios de óptica de Fourier , síntesis de sonido , interferometría estelar , pinzas ópticas y elementos ópticos difractivos (DOE), a menudo es importante conocer la fase de frecuencia espacial de una fuente de onda observada. Para reconstruir esta fase , se puede utilizar el algoritmo adaptativo-aditivo (o algoritmo AA ), que deriva de un grupo de algoritmos adaptativos (de entrada-salida). El algoritmo AA es un algoritmo iterativo que utiliza la transformada de Fourier para calcular una parte desconocida de una onda que se propaga, normalmente la fase de frecuencia espacial (espacio k). Esto se puede hacer cuando se dan las contrapartes conocidas de la fase, generalmente una amplitud observada ( espacio de posición ) y una amplitud inicial supuesta (espacio k). Para encontrar la fase correcta, el algoritmo utiliza la conversión de error, o el error entre las intensidades deseadas y teóricas .

El algoritmo

Historia

El algoritmo adaptativo-aditivo se creó originalmente para reconstruir la fase de frecuencia espacial de la intensidad de la luz en el estudio de la interferometría estelar . Desde entonces, el algoritmo AA se ha adaptado para su uso en los campos de la óptica de Fourier por Soifer y el Dr. Hill, la materia blanda y las pinzas ópticas por el Dr. Grier, y la síntesis de sonido por Röbel.

Algoritmo

  1. Defina la amplitud de entrada y la fase aleatoria.
  2. Transformada de Fourier directa
  3. Separe la amplitud y la fase transformadas.
  4. Comparar la amplitud/intensidad transformada con la amplitud/intensidad de salida deseada.
  5. Comprobar las condiciones de convergencia
  6. Mezcla la amplitud transformada con la amplitud de salida deseada y combínala con la fase transformada.
  7. Transformada inversa de Fourier
  8. Separe la nueva amplitud y la nueva fase.
  9. Combinar la nueva fase con la amplitud de entrada original
  10. Volver a la transformada de Fourier directa

Ejemplo

Para el problema de reconstruir la fase de frecuencia espacial ( espacio k ) para una intensidad deseada en el plano de la imagen ( espacio x ). Supongamos que la amplitud y la fase inicial de la onda en el espacio k sonA0{\displaystyle A_{0}}yϕnortek{\displaystyle \phi _{n}^{k}}respectivamente. Transformar la onda en el espacio k al espacio x mediante la transformada de Fourier.

A0miiϕnortekFFTAnorteFmiiϕnorteF{\displaystyle A_{0}e^{i\phi _{n}^{k}}{\xrightarrow {FFT}}A_{n}^{f}e^{i\phi _{n}^{f}}}

Luego compare la intensidad transformada.InorteF{\displaystyle I_{n}^{f}}con la intensidad deseadaI0F{\displaystyle I_{0}^{f}}, dónde

InorteF=(AnorteF)2,{\displaystyle I_{n}^{f}=\left(A_{n}^{f}\right)^{2},}
ε=(InorteF)2(I0)2.{\displaystyle \varepsilon ={\sqrt {\left(I_{n}^{f}\right)^{2}-\left(I_{0}\right)^{2}}}.}

Controlarε{\displaystyle \varepsilon }contra los requisitos de convergencia. Si no se cumplen los requisitos, entonces mezcle la amplitud transformada.AnorteF{\displaystyle A_{n}^{f}}con la amplitud deseadaAF{\displaystyle A^{f}}.

A¯norteF=[aAF+(1a)AnorteF],{\displaystyle {\bar {A}}_{n}^{f}=\left[aA^{f}+(1-a)A_{n}^{f}\right],}

donde a es la relación de mezcla y

AF=I0{\displaystyle A^{f}={\sqrt {I_{0}}}}.

Tenga en cuenta que a es un porcentaje, definido en el intervalo 0 ≤ a ≤ 1.

Combine la amplitud mixta con la fase en el espacio x y la transformada inversa de Fourier .

A¯FmiiϕnorteFiFFTA¯nortekmiiϕnortek.{\displaystyle {\bar {A}}^{f}e^{i\phi _{n}^{f}}{\xrightarrow {iFFT}}{\bar {A}}_{n}^{k}e^{i\phi _{n}^{k}}.}

SeparadoA¯nortek{\displaystyle {\bar {A}}_{n}^{k}}yϕnortek{\displaystyle \phi _{n}^{k}}y combinarA0{\displaystyle A_{0}}conϕnortek{\displaystyle \phi _{n}^{k}}Incrementar el bucle en uno.nortenorte+1{\displaystyle n\to n+1}y repetir.

Límites

  • Sia=1{\displaystyle a=1}Entonces el algoritmo AA se convierte en el algoritmo de Gerchberg-Saxton .
  • Sia=0{\displaystyle a=0}entoncesA¯nortek=A0{\displaystyle {\bar {A}}_{n}^{k}=A_{0}}.

Véase también

Referencias

  • Dufresne, Eric; Grier, David G; Spalding (diciembre de 2000), "Matrices de pinzas ópticas holográficas generadas por computadora", Review of Scientific Instruments , 72 (3): 1810, arXiv : cond-mat/0008414 , Bibcode : 2001RScI...72.1810D , doi : 10.1063/1.1344176 , S2CID 14064547 .
  • Grier, David G (10 de octubre de 2000), Algoritmo adaptativo-aditivo.
  • Röbel, Axel (2006), "Modelado aditivo adaptativo con trayectorias de parámetros continuos", IEEE Transactions on Audio, Speech, and Language Processing , 14 (4): 1440– 1453, doi : 10.1109/TSA.2005.858529 , S2CID 73476 .
  • Röbel, Axel, Síntesis adaptativa-aditiva del sonido , ICMC 1999, CiteSeerX 10.1.1.27.7602 {{citation}}: CS1 mantenimiento: ubicación ( enlace )
  • Soifer, V. Kotlyar; Doskolovich, L. (1997), Métodos iterativos para el cálculo de elementos ópticos difractivos , Bristol, PA: Taylor & Francis, ISBN 978-0-7484-0634-0
  • Presentación del laboratorio de David Grier sobre pinzas ópticas y fabricación del algoritmo AA.
  • Síntesis aditiva adaptativa para sonido no estacionario. Dr. Axel Röbel.
  • Laboratorios Hill, Universidad de Maryland, College Park .]