Ecuación de Helmholtz
Propagación libre
Aproximación de Fresnel y Fraunhoffer
David Mas Candela
Sistemas lineales e invariantes
La propagación libre está regulada por la ecuación de onda,
que es una ecuación lineal e invariante a desplazamiento
El paso por una apertura es una condición de contorno que se
impone al campo, de modo que no afecta a la estructura del
mismo
Por tanto, podemos tratar el problema de la propagación
como un sistema LSI
INPUT OUTPUT
Resolviendo la Ec. Helmholtz
El problema de la propagación consiste en encontrar una
solución a la ecuación de Helmholtz:
2π
( ∇ + k ) u(r) = 0
2 2
k=
λ
La función que es solución a este tipo de ecuaciones es de la
forma:
u(x, y, z) = fx(x) × fy(y) × fz(z)
De modo que el Laplaciano sobre la función queda:
2 2 2
2 ∂ u ∂ u ∂ u
∇ u= 2 + 2 + 2
∂x ∂y ∂z
[Link]
[Link]
Resolviendo la Ec. Helmholtz
Con lo que la ecuación de Helmholtz desarrollada
completamente queda:
f′x′(x)fy(y)fz(z) + fx(x)f′y′(y)fz(z) + fx(x)fy(y)f′z′(z) + k 2fx(x)fy(y)fz(z) = 0
O, de modo más compacto:
f′x′(x) f′y′(y) f′z′(z)
+ + + k2 = 0
fx(x) fy(y) fz(z)
Lo que implica que cada sumando es constante ya que si no
fuese así, no podría cancelarse con los otros términos al no
haber dependencias cruzadas.
Por tanto…
[Link]
[Link]
Resolviendo la Ec. Helmholtz
Para cada una de las funciones separadas:
d2 2
fx (x) + kx fx(x) = 0
dx 2
d2 2
fy(y) + k y fy (y) = 0
dy 2
d2 2
fz (z) + kz fz(z) = 0
dz 2
Lo que nos lleva a que la constante k puede escribirse como:
kx2 + ky2 + kz2 = k 2
[Link]
Resolviendo la Ec. Helmholtz
Cada una de las ecuaciones diferenciales separadas que hemos
obtenido tiene como solución una función exponencial, de forma
que la solución completa puede construirse mediante el
producto de las tres soluciones separadas.
u(x, y, z) = e j(kx x+ky y+kz z)
= e j(kx x+ky y)e jkz z
j(kx x+ky y) ±jz k 2
− k 2
− k 2
=e e x y
Esta solución representa una onda propagándose o decayendo
exponencialmente, tal y como veremos a continuación
[Link]
Resolviendo la Ec. Helmholtz
La solución propuesta es una solución particular de la ec. de
Helmholtz. La solución general la construiremos mediante la
combinación lineal de todas las soluciones particulares
posibles, es decir:
+∞ +∞
∫−∞ ∫−∞
j(kx x+ky y) ± jz k 2
− k 2
− k 2
u(x, y, z) = A0(kx, ky)e e x y
dkxdky
Dónde, considerando:
u0(x, y) = u(x, y, z) |z=0
Podemos escribir:
+∞ +∞
∫−∞ ∫−∞
u0(x, y) = A0(kx, ky) e j(kx x+ky y) dkxdky
[Link]
Resolviendo la Ec. Helmholtz
Esta solución mediante superposición de ondas planas es la
base de la llamada Óptica de Fourier ya que las ecuaciones
obtenidas representan una transformada de Fourier.
+∞ +∞
∫−∞ ∫−∞
±jz k 2 − kx2 − ky2 j(k x+k y)
u(x, y, z) = A0(kx, ky)e e x y dkxdky
+∞ +∞
∫−∞ ∫−∞
u0(x, y) = A0(kx, ky) e j(kx x+ky y) dkxdky
+∞ +∞
j(kx x + ky y)
Az (kx, ky, z)
∫−∞ ∫−∞
u(x, y, z) = e dkxdky
(λ λ ) (λ λ ) ( )
α β α β
Az , , z = A0 , ,0 exp jz k 2 − kx2 − ky2
[Link]
Propagación de la luz
Analicemos atentamente la solución encontrada:
+∞ +∞
∫−∞ ∫−∞
jk⋅r
u(x, y, z) = A0(kx, ky)e dkxdky
+∞ +∞
∫−∞ ∫−∞
±jz k 2 − kx2 − ky2 j(k x+k y)
= A0(kx, ky)e e x y dkxdky
Podemos obtener un +∞ sentido
+∞ físico a la solución anterior reescribiendo
en la dirección de∫−∞ ∫−∞vector se expresará en función de los cosenos
k, expresándolo= como su módulo )e j(kx x+ky y)dkxpor
Az(kx, kymultiplicado dky el vector unitario
k. Este
directores de dicho vector unitario:
2π
( z)
k= cos α̂ u + cos β ̂ u + cos γ ̂ u
x y
λ
Propagación de la luz
𝛾!
𝛼! 𝛽$
(cos α)̂ 2 + (cos β)̂ 2 + (cos γ)̂ 2 = 1
Propagación de la luz
Para simplificar la notación, consideraremos el siguiente cambio:
[ ]
cos ̂
α, cos ̂ cos γ ̂ ≡ [α, β, γ]
β,
Con lo que:
2π
k=
λ ( αux + βuy + γuz)
( λ ) (
2
2π
)
2
k = kx2 + ky2 + kz2 = α 2
+ β 2
+ γ 2
(α + β + γ ) = 1
2 2 2
Propagación de la luz
NOTA:
En el desarrollo que se realiza en el libro “Introduction to Fourier
Optics” de J.W. Goodman, se realiza un cambio de variable en todos
los factores de la integral con el fin de hacer explícita la dependencia
angular y mantener la coherencia de las ecuaciones:
α β
kx = 2π ky = 2π
λ λ
En nuestro caso con el fin de facilitar el seguimiento de las ecuaciones
y a costa del rigor matemático, no efectuaremos dicho cambio, salvo
en la raíz cuadrada para simplificar la notación.
No obstante, se marcarán los puntos “sensibles” para evitar
confusiones.
Propagación de la luz
Tal y como habíamos deducido:
+∞ +∞
j(kx x + ky y)
Az (kx, ky, z)
∫−∞ ∫−∞
u(x, y, z) = e dkxdky
( λ )
2π
Az (kx, ky, z) = A0 (kx, ky,0) exp j z 1 − α 2 − β 2
Az es conocido como el espectro angular de la distribución u
Espectro: conjunto de ondas que conforman la señal
Angular: Por la interpretación en forma de ondas planas que viajan en
diferentes direcciones.
Propagación de la luz
Y, en el plano z=0:
+∞ +∞
∫−∞ ∫−∞
u0(x, y) = A0(kx, ky) e j(kx x+ky y) dkxdky
Donde:
+∞ +∞
−j(kx x + ky y)
A0 (kx, ky) =
∫−∞ ∫−∞
u0(x, y) e dxdy
Esto significa que la distribución de entrada también puede
descomponerse en una suma ponderada de ondas planas con
diferentes orientaciones.
Propagación de la luz
Si nos fijamos atentamente en la expresión de 𝐴!
( λ )
2π
Az (kx, ky, z) = A0 (kx, ky,0) exp j z 1 − α 2 − β 2
Este resultado muestra que, cuando (α, β) cumplen:
α2 + β2 <1
el efecto de la propagación en z is un cambio relativo de las fases de
las diferentes componentes del espectro angular.
PERO, si (α, β) cumplen:
α + β >1
2 2
La raíz cuadrada es compleja y se necesita una interpretación
alternativa
Propagación de la luz
α2 + β2 >1
En este caso no podemos hablar de una onda en propagación ya que,
al tener una raíz imaginaria:
1 − α 2 − β 2 = iω
La expresión:
( λ )
2π
Az (kx, ky, z) = A0 (kx, ky,0) exp j z 1 − α 2 − β 2
Se convierte en:
Az (kx, ky, z) = A0 (kx, ky,0) exp (−μz)
Propagación de la luz
α2 + β2 >1
Az (kx, ky, z) = A0 (kx, ky,0) exp (−μz)
Dado que μ se ha definido como un número real y positivo estas
componentes de la onda se atenúan rápidamente. Estas ondas reciben
el nombre de ondas evanescentes.
Recordemos que α y β son los cosenos directores de los ángulos que
forman el vector de propagación k y la dirección de observación, r.
Por tant, la condición de onda evanescente se cumple sólo para
direcciones de propagación en las que α y β son cercanas a 1, es decir,
próximas a la perpendicular de la dirección de propagación.
Propagación de la luz
En el caso del paso de luz a través de una apertura, estas ondas se
darán en las paredes interiores a la apertura, en el lado de la
propagación:
Propagación de la luz
Las ondas evanescentes pueden “capturarse” en un segundo medio
del mismo (o muy similar) índice de refracción que esté lo
suficientemente cerca, de modo que se convierta en una nueva onda
propagada. Este procedimiento es habitual en acoplamientos de
fibras.
Propagación de la luz
Si nos centramos sólo en las ondas que se propagan, recordemos que
+∞ +∞
j(kx x + ky y)
Az (kx, ky, z)
∫−∞ ∫−∞
u(x, y, z) = e dkxdky
Que se convierte en :
+∞ +∞
j(kx x + ky y)
A0 (kx, ky,0) e
∫−∞ ∫−∞
u(x, y, z) =
( ) ( λ )
2π
× circ α + β exp j z 1 − α 2 − β 2 dkx ky
2 2
{0
1 for r < 1
circ(r) =
elsewhere
Hemos introducido la función círculo para enfatizar que hay una
serie de frecuencias que no se propagan y que son cortadas por
dicha función (en realidad no es necesaria)
Propagación de la luz
+∞ +∞
j(kx x + ky y)
A0 (kx, ky,0) e
∫−∞ ∫−∞
u(x, y, z) =
( ) ( λ )
2π
× circ α + β exp j z 1 − α 2 − β 2 dkx ky
2 2
Esta ecuación significa que:
- Una onda propagada puede construirse como una suma conveniente
de ondas con diferentes amplitudes y fases.
- NO TODAS las componentes contribuyen al campo final.
- Parte de la información se pierde debido a la propagación.
- La propagación de la luz (incluyendo la formación de imágenes)
implica pérdida de información.
- Las imágenes perfectas no son estigmáticas sino LIMITADAS POR
DIFRACCIÓN
Propagación de la luz
+∞ +∞
j(kx x + ky y)
A0 (kx, ky)
∫−∞ ∫−∞
u0(x, y) = e dkxdky
+∞ +∞
−j(kx x + ky y)
A0 (kx, ky) =
∫−∞ ∫−∞
u0(x, y) e dxdy
Si consideramos que uo es una función delta, contendrá todas las
frecuencias, con lo que A0=cte para todos los valores de las variables,
según lo visto en el tema anterior.
Propagación de la luz
+∞ +∞
j(kx x + ky y)
A0 (kx, ky,0) e
∫−∞ ∫−∞
u(x, y, z) =
( ) ( λ )
2π
× circ α + β exp j z 1 − α 2 − β 2 dkx ky
2 2
La función circ hace que la constante ya no se extienda en todo el rango
de variables, lo que provocará la desaparición de las frecuencias más
altas.
Propagación de la luz
Dado que estamos trabajando en dominios espaciales y frecuenciales
simultáneamente renombraremos las variables de modo que:
- Usaremos letras latinas para variables espaciales
- Usaremos letras griegas para variables frecuenciales
1 − (λξ) − (λη)
2 2
kx = 2πξ ky = 2πη γ=
Con lo que la ecuación de propagación queda:
+∞ +∞
A0 (ξ, η,0) exp (j2π (ξx + ηy))
∫−∞ ∫−∞
u(x, y, z) =
( ) ( λ )
2π
(λξ) + (λη) exp j z 1 − (λξ) − (λη)
2 2 2 2
× circ dξdη
Propagación de la luz
Reescribamos las ecuaciones con una variable para identificar mejor
los términos:
( λ )
+∞
2π
∫−∞
A0 (ξ,0) exp j z 1 − (λξ) exp (j2πξx) dξ
2
u(x, z) =
+∞
∫−∞
A0 (ξ) = u0(x0) exp (−j2πξx0) dx0
Por lo que:
( λ )
+∞ +∞
2π
∫−∞ ∫−∞
u0 (x0,0) exp (−j2πξx0) exp j z 1 − (λξ) exp (j2πξx) dx0dξ
2
u(x, z) =
O equivalentemente:
( λ )
+∞ +∞
2π
∫−∞ ∫−∞
u0 (x0,0) exp (j2πξ(x − x0)) exp j z 1 − (λξ) dx0dξ
2
u(x, z) =
Propagación de la luz
En el tema anterior definimos la Función de Transferencia como:
ℱ [h (x)] = H (ξ)
Siendo h(x) la respuesta a una función impulsional.
La respuesta impulsional se obtiene como respuesta a una función 𝛿
Si 𝑢" 𝑥" =δ 𝑥"
( λ )
+∞
2π
∫−∞
exp (j2πξx) exp j z 1 − (λξ) dξ
2
h(x, z) =
Que no es más que la transformada de Fourier inversa de la función
𝐻(𝜉)
Propagación de la luz
Con lo que, la función de transferencia de la propagación libre
(recuperamos momentáneamente la notación bidimensional)
( )
1 − (λξ) − (λη)
2π 2 2 2 2 1
exp j z ξ +η <
H (ξ, η) = λ λ
0 resto
Es decir, es una función que solo deja pasar las frecuencias más bajas,
eliminando las más altas.
Aproximación de Fresnel
Con el fin de simplificar la integral, se toman diferentes
aproximaciones. Si retomamos la ecuación de propagación
( λ )
+∞ +∞
2π
∫−∞ ∫−∞
u0 (x0,0) exp (j2πξ(x − x0)) exp j z 1 − (λξ) dx0dξ
2
u(x, z) =
Vemos que, dado que λ para la luz visible es suficientemente
pequeña (<1 micra) la raíz cuadrada puede aproximarse mediante:
(λξ) (λη)
2 2
1 − (λξ) − (λη) ≈ 1 −
2 2
−
2 2
(λξ)
2
1 − (λξ) ≈ 1 −
2
2
Aproximación de Fresnel
Si cambiamos A0(ξ,η,0) por su expresión en términos de u(x0, y0,0) y
reducimos la notación a una sola variable:
[ λ ] ∫−∞ ∫−∞ ( )
+∞ +∞
2π
u(xz, z) = exp i z u0 (x0,0) circ (λξ)2 exp [j2πξ (xz − x0)] exp [−jπλzξ 2] dξdx0
Que puede integrarse en 𝜉 (Ver deducción):
[ λz ]
+∞
π
∫−∞
u0 (x0,0) exp j (xz − x0) dx0
2
u(xz, z) ∝
Que se conoce con el nombre de integral de difracción de Fresnel.
Aproximación de Fresnel
[ λz ]
+∞
π
∫−∞
u0 (x0,0) exp j (xz − x0) dx0
2
u(xz, z) =
[ λz ] [ λz ]
+∞
π 2 2π
∫−∞
= u0 (x0,0) exp j (xz + x0 ) exp −j
2
( xz x0) dx0
[ λz ] [ λz ]
+∞
π 2 2π
∫−∞
∝ u0 (x0,0) exp j (x0 ) exp −j ( xz x0) dx0
La integral de Fresnel:
- Funciona para “campo lejano” (comparado con λ)
- Da la expresión del campo propagado en función de u0
- Puede interpretarse como una convolución
- Puede interpretarse como una Transformada de Fourier
Aproximación de Fresnel
La integral de Fresnel como convolución:
[ λz ]
+∞
π
∫−∞
u0 (x0,0) exp j (xz − x0) dx0
2
u(xz, z) =
Donde, podemos identificar:
+∞
∫−∞
uz (xz) = u0 ⊗ h = u0 (x0)h (x0 − xz) du
Siendo la respuesta impulsional y la función de transferencia
asociadas a la propagación libre:
[ λz ]
π 2
H (ξ, z) ∝ exp [−jπλzξ 2]
TF
h(x) = exp j x
Aproximación de Fresnel
La integral de Fresnel como TF:
( λz )] [ λz ]
+∞
[
π 2 2π
∫−∞
u(xz, z) ∝ u0 (x0,0) exp j x0 exp −j ( xz x0) dx0
Es decir, la TF del objeto multiplicado por una onda esférica.
Si, en el laboratorio tenemos una imagen de la propagación de Fresnel
y tomamos un punto brillante xz=5 cm a una distancia z=10 cm
tenemos que ese punto corresoponde a una frecuencia
xz = 5 ⎫⎪ 5
⎬ξ =
z = 10 ⎪⎭ 10 λ
Conforme propagamos (cambiamos z), el pico representando una
frecuencia determinadase desplazará y cambiará su posición.
Aproximación de Fraunhofer
La expresión de Fresnel:
[ λz ]
+∞
π
∫−∞
u0 (x0,0) exp j (xz − x0) dx0
2
u(xz, z) =
Para aperturas limitadas y distancias muy grandes:
π(xo2 + yo2)max
z≫
2λ
Se convierte en:
( λz ) ∫−∞ [ λz ]
π 2 +∞ 2π
u(xz, z) = exp j xz u0 (x0,0) exp −j xz x0 dx0
Expresión conocida como aproximación de Fraunhofer
Aproximación de Fraunhofer
( λz ) ∫−∞ [ λz ]
+∞
π 2 2π
u(xz, z) = exp j xz u0 (x0,0) exp −j xz x0 dx0
La aproximación de Fraunhofer:
- Es válida para campos muy lejanos
- Puede interpretarse como la transformada de Fourier de la
distribución de entrada con las coordenadas:
(ξ, η) = ( , )
xz yz
λz λz
Aproximación de Fraunhofer
La aproximación de Fraunhofer nos dice que a distancias
suficientemente lejanas, el patrón de luz propagado será equivalente a
la transformada de Fourier de la función de entrada
convenientemente escalada.
En el caso de un haz plano incidiendo sobre una apertura
rectangular, si nos situamos lo suficientemente lejos, veremos la TF de
la apertura, que actuará como función de entrada.
Si miramos una fuente de luz puntual a través de una tela, veremos
una figura en forma de cruz. En realidad estamos observando la TF de
una apertura cuadrada.
En realidad la luz pasa por múltiples aperturas, pero, dado que es
invariante a desplazamiento, veremos la superposición de todas las
TFs en un único punto:
[Link]
[Link]
Aproximación de Fraunhofer
La aproximación de Fraunhofer además informa de un hecho
peculiar.
Supongamos que ponemos un objeto g(x) pegado a la apertura
cuadrada de anchura L. Matemáticamente, eso se expresa mediante
un producto de las dos funciones:
(L)
x
u0(x0) = rect × g(x0)
A una distancia suficientemente alejada obtendremos la TF del
producto, esto es, la convolución de la TF de la apertura con la TF de
la función.
uz(xz) ≈ F(rect) ⊗ F(g)
Es decir, la distribución en campo lejano siempre estará contaminada
por la TF de la apertura.
En el siguiente tema veremos más casos.
