Optica: optica fisica

coherencia temporal y espacial

Figura 1. El experimento de Young, en la figura (a) tenemos la fuente luminosa que ilumina el plano paralelo que tiene dos rendijas que difractan la luz al pasar por ellas. En la figura (b) la fuente es puntual y pasa por los orificios S1 Y S2.

Si en la figura 1, la fuente primaria que se encuentra en el primer plano de la imagen (a), que llamaremos S, se encoge hasta convertirse en una fuente puntual sobre el eje central como en la imagen (b) y tiene un ancho de banda de frecuencia finito, los efectos de coherencia temporal predominarán.

Las perturbaciones ópticas en ${S}_{1}\$ y ${S}_{2}\$ son entonces idénticas, donde ${S}_{1}\$ y ${S}_{2}\$ son las ventanas que muestra la figura (b) en el plano intermedio. La coherencia mutua ${\tilde {\Gamma }}_{12}(\tau )=\left\langle \mathbf {V} _{1}(t)\mathbf {V} _{2}^{*}(t+\tau )\right\rangle$ entre los dos puntos donde las perturbaciones son idénticas es en este caso la autocoherencia del campo.

De aquí que podamos escribirla como

${\tilde {\Gamma }}(S_{1},S_{2},\tau )={\tilde {\Gamma }}_{12}(\tau )={\tilde {\Gamma }}_{11}(\tau )$ o´ ${\tilde {\gamma }}_{12}(\tau )={\tilde {\gamma }}_{11}(\tau ).$

lo mismo obtenemos cuando ${S}_{1}\$ y ${S}_{2}\$ se combinan y ${\tilde {\gamma }}_{11}(\tau )$ se denomina grado de coherencia temporal complejo donde $\tau \$ es el intervalo de tiempo que se separan.

Esta interpretación es un poco diferente para un interferómetro de división de amplitud (interferometro de Michelson[1]), donde $\tau \$ es igual a la diferencia de recorrido dividida por c. La expresión para la irradiancia I $I=I_{1}+I_{2}+2(I_{1}I_{2})^{\frac {1}{2}}Re{\tilde {\gamma }}_{12}(\tau ),$ contendra entonces a ${\tilde {\gamma }}_{11}(\tau )$ en lugar de ${\tilde {\gamma }}_{12}(\tau )$.

Ahora si una onda luminosa con el interferómetro se divide en dos perturbaciones del campo idénticas de la forma ${\tilde {\mathbf {V} }}(t)=\mathbf {V} _{0}e^{i\phi (t)}$ donde $\phi \ (t)$ se refiere a la fase de la onda y la tilde expresa en V que puede tener una representacion compleja, si el interferometro después las vuelve a combinar para generar una distribución de franjas, entonces

${\tilde {\gamma }}_{11}(\tau )={\frac {\left\langle {\tilde {\mathbf {V} }}(t+\tau ){\tilde {\mathbf {V} }}^{*}(t)\right\rangle }{{\tilde {|\mathbf {V} |}}^{2}}}={\frac {\left\langle e^{i\phi (t+\tau )}e^{-i\phi (t)}\right\rangle }{{\tilde {|\mathbf {V} |}}^{2}}}{\tilde {|\mathbf {V} |}}^{2}=\left\langle e^{i\phi (t+\tau )}e^{-i\phi (t)}\right\rangle$

Para saber que tipo de resultados nos generará la ecuación anterior y cómo varía veamos primero que pasa con una onda plana monocromática de longitud de coherencia infinita Error al representar (error de sintaxis): \phi(t)= \mathbf{k • r}-\omega t

donde Error al representar (error de sintaxis): \Delta\phi=\mathbf{k} {•} \mathbf{r}-\omega(t+\tau)-\mathbf {k} •\mathbf{r}+\omega t=-\omega\tau que nos lleva a ${\tilde {\gamma }}_{11}(\tau )=e^{i\omega \tau }=\cos \omega \tau -i\sin \omega \tau$

de tal modo que $|{\tilde {\gamma }}_{11}|=1$, llamada coherencia completa.

Pasemos al caso opuesto, con una onda monocromatica con $\tau >\$ longitud de coherencia de tal modo que $\Delta \phi \$ sera al azar, variando entre 0 y $2\pi \$ de tal manera que la integral se promedia a cero, $|{\tilde {\gamma }}_{11}(\tau )|=0$, llamada incoherencia completa.

Es importante entender que la transformada de Fourier de la función de autocoherencia es el espectro energético que describe la distribución de la energía espectral de la luz, esto lo veremos más adelante cuando trabajemos con las aplicaciones del teorema que nos muestra la relación entre el grado de coherencia complejo en una región del espacio y la distribución de la irradiancia correspondiente en la fuente extensa que da lugar a los campos luminosos, este teorema es el Teorema de Van Cittert Zernike,

Teorema de Cittert-Zernike

Archivo:Ls3 unit6A76.gif

Figura 2.Fuente arbitraria S emisora de luz, en nuestras ecuaciones $P_{1}S\$ es $R_{1}\$ y $P_{2}S\$ es $R_{2}\$

Se presenta una fuente incoherente extensa cuasimonocromática, S, situada en el plano $\sigma \$ (representado por el segmento cuadriculado que se ve en la figura 2) y con una irradiancia proporcionada por I(\epsilon,\eta). También se muestra una pantalla de observación en la que se hallan dos puntos, $P_{1}\$ y $P_{2}\$ a una distancia $R_{1}\$ y $R_{2}\$ respectivamente, desde un elemento diminuto de S.

Vamos a determinar la intensidad mutua $J_{12}\$ y el grado de coherencia complejo $j_{12}\ (\mu _{12})$ para los puntos de $P_{1}\$ y $P_{2}\$ en la pantalla de observación iluminada por un haz de luz proviniente de $\sigma \$.

Por simplicidad $\sigma \$ sera tomada como un plano paralelo al plano de observación y vamos a suponer que el medio entre la fuente y la pantalla es homogéneo. También que las dimensiones lineales de $\sigma \$ son pequeñas en comparación con la distancia $0^{\prime }0\$ entre los planos de la fuente y la pantalla, y que el ángulo entre las líneas que unen un punto de la fuente S a $P_{1}\$ y $P_{2}\$ es pequeño.

Imaginemos entonces que la fuente se divide en elementos $d\sigma _{1}\ ,d\sigma _{2},...$ centrado en los puntos $S_{1}\ ,S_{2},...$ de dimensiones lineales pequeñas comparadas con la longitud de onda ${\tilde {\lambda }}\$ donde la tilde en este caso representa el promedio de $\lambda \$. Es decir, si $V_{m1}(t)\$ y $V_{m2}(t)\$ son las perturbaciones del campo complejo en $P_{1}\$ y $P_{2}\$, debido al elemento $d\sigma _{1}\$, el total de las perturbaciones en estos puntos son ahora las vibraciones de luz que surgen de los diferentes elementos de la fuente, puede ser asumida como estadísticamente independientes (incoherentes entre sí), de modo que $V_{1}\ =\sum _{m}V_{m1}(t)$ y $V_{2}\ =\sum _{m}V_{m2}(t)$,

siendo la intensidad mutua

$J(P_{1},P_{2})=\left\langle {V}_{1}(t){V}_{2}^{*}(t)\right\rangle$

Podemos escribir la ecuación anterior como una suma, separandola en dos partes al no ser necesariamente los subindices iguales. Es decir una cuando los subindices son iguales y otra cuando no lo son.

$J(P_{1},P_{2})=\sum _{m}\left\langle {V}_{m1}(t){V}_{m2}^{*}(t)\right\rangle +\sum \sum _{nm}\left\langle {V}_{m1}(t){V}_{n2}^{*}(t)\right\rangle$

Al ser una fuente incoherente, cuando m es distinto de n la suma vale cero.

$J(P_{1},P_{2})=\sum _{m}\left\langle {V}_{m1}(t){V}_{m2}^{*}(t)\right\rangle +\sum \sum _{nm}0$

Si $R_{m1}\$ y $R_{m2}\$ son las distancias de $P_{1}\$ y $P_{2}\$ del elemento de origen $d\sigma _{m}\$ las perturbaciones del campo las podemos escribir como

$V_{m1}(t)=A_{m}\left(t-{\frac {R_{m1}}{v}}\right){\frac {e^{-2i{\tilde {\nu }}\pi (t-{\frac {R_{m1}}{v}})}}{R_{m1}}}$

;

$V_{m2}(t)=A_{m}\left(t-{\frac {R_{m2}}{v}}\right){\frac {e^{-2i{\tilde {\nu }}\pi (t-{\frac {R_{m2}}{v}})}}{R_{m2}}}$

Donde $|A_{m}|\$ caracteriza a la intensidad I y el argumento $A_{m}\$ la fase de la de la radiación del elemento m-esimo, v es la velocidad de la luz en el medio entre la fuente y la pantalla. Por lo tanto

$\left\langle \mathbf {v} _{m1}(t)\mathbf {v} _{m2}^{*}(t)\right\rangle =\left\langle A_{m}\left(t-{\frac {R_{m1}}{v}}\right)A_{m}^{*}\left(t-{\frac {R_{m1}}{v}}\right)\right\rangle {\frac {e^{2i\pi {\tilde {\nu }}{\frac {(R_{m1}-R_{m2})}{v}}}}{R_{m1}R_{m2}}}=\left\langle A_{m}\left(t\right)A_{m}^{*}\left(t-{\frac {R_{m2}-R_{m1}}{v}}\right)\right\rangle {\frac {e^{2i\pi {\tilde {\nu }}{\frac {(R_{m1}-R_{m2})}{v}}}}{R_{m1}R_{m2}}}$

si la diferencia de caminos opticos $(R_{m1}-R_{m2})$ es pequeña comparada con la longitud de coherencia de la Luz podemos despreciar el termino ${\frac {(R_{m1}-R_{m2})}{v}}$ en el argumento de $A_{m}^{*}\$ y obtener

$J(P_{1},P_{2})=\sum _{m}\left\langle A_{m}\left(t\right)A_{m}^{*}\left(t\right)\right\rangle {\frac {e^{2i\pi {\tilde {\nu }}{\frac {(R_{m1}-R_{m2})}{v}}}}{R_{m1}R_{m2}}}$

Se sigue que $\left\langle A_{m}\left(t\right)A_{m}^{*}\left(t\right)\right\rangle$ caracteriza la intensidad de radiación de los elementos de la fuente $d\sigma \$ asi que podemos sustituirla por I(S) que es la intensidad por unidad de área de la fuente, dejándonos así

$J(P_{1},P_{2})=\int _{\sigma }I(S){\frac {e^{i{\tilde {k}}(R_{1}-R_{2})}}{R_{1}R_{2}}}dS$

Regresando a la notación de $R_{1}\$ y $R_{2}\$ que ya definimos y el numero de onda en el medio como ${\tilde {k}}={\frac {2\pi {\tilde {\nu }}}{v}}={\frac {2\pi }{\tilde {\lambda }}}$, siendo las cantidades con tilde los promedios.

Obtenemos el grado complejo de coherencia

$j(P_{1},P_{2})={\frac {J(P_{1},P_{2})}{\sqrt {I(P_{1})I(P_{2})}}}={\frac {1}{\sqrt {I(P_{1})I(P_{2})}}}\int _{\sigma }I(S){\frac {e^{i{\tilde {k}}(R_{1}-R_{2})}}{R_{1}R_{2}}}dS...(1)$

Con las intensidades de $P_{1}\$ y $P_{2}\$

$I(P_{1})=J(P_{1},P_{1})=\int _{\sigma }{\frac {I(S)}{{R_{1}}^{2}}}dS$

$I(P_{2})=J(P_{2},P_{2})=\int _{\sigma }{\frac {I(S)}{{R_{2}}^{2}}}dS$

El cálculo da como resultado una integral con una estructura conocida cuya forma y cuyos resultados son iguales a los de una integral de difracción, con tal de que cada término se vuelva a interpretar correctamente. Notemos que esta integral (1), se puede calcular de una manera diferente, sobre la base del principio de Huygens-Fresnel, de la alteración compleja en el patrón de difracción derivado de la difracción de una onda esférica en una apertura en una pantalla opaca. De manera más precisa, la integral (1) implica que la igualdad de grado de coherencia complejo, que describe la correlación de las vibraciones en un punto fijo $P_{2}\$ y el punto variable $P_{1}\$ en un plano iluminado por una fuente primaria extensa cuasimonocromática (incoherente), es igual a la normalización de la amplitud compleja en el punto correspondiente $P_{1}\$ en un cierto patrón de difracción, centrado en $P_{2}\$. Este patrón se obtendría en la sustitución de la fuente por una abertura de difracción de la misma forma y tamaño de la fuente, y llenado con una onda esférica convergente para $P_{2}\$, la apertura es proporcional a la distribución de intensidad a través de la fuente.

Aplicaciones.

En la mayoría de las aplicaciones de la intensidad I(S) que puede suponerse como independiente de la posición de S en la superficie (la intensidad es uniforme). El problema de difracción correspondiente es entonces el de la difracción de una onda esférica de amplitud de uniforme por una abertura de la misma forma y tamaño de la fuente.

Sean $(\epsilon ,\eta )\$ las coordenadas de un punto de fuente típico S, en el plano sigma con el origen en 0, y $(X_{1},Y_{1})\$ y $(X_{2},Y_{2})\$ las coordenadas de $P_{1}\$ y $P_{2}\$, en el plano paralelo A y con origenel 0´. Entonces, si R representa la distancia de 0 a 0´

${R_{1}}^{2}=({X_{1}-\epsilon })^{2}+({Y_{1}-\eta })^{2}+R^{2}$

$\psi ={\frac {{\tilde {k}}[(X_{1}^{2}+Y_{1}^{2})-({X_{2}}^{2}-{Y_{2}}^{2})]}{R2}}$

Entonces, si las dimensiones lineales de una fuente y la distancia entre $P_{1}\$ y $P_{2}\$ son pequeñas comparadas con la distancia de esos puntos a la fuente, el grado de coherencia $j_{12}\$ es igual al valor absoluto de la transformada de Fourier normalizada de la función intensidad de la fuente.

$j(P_{1},P_{2})={\frac {e^{i\psi }\int \int _{\sigma }I(\epsilon ,\eta )e^{i{\tilde {k}}({\frac {(X_{1}-X_{2})\epsilon }{r}}+{\frac {(Y_{1}-Y_{2})\eta }{r}}}d\epsilon d\eta }{\int \int _{\sigma }I(\epsilon ,\eta )d\epsilon d\eta }}...(2)$

Para una fuente circular uniforme de radio Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): \rho\ con centro en O, la expresión anterior la podemos reescribir como

$j(P_{1},P_{2})=({\frac {2J_{1}(v)}{v}})e^{i\psi }$

donde

$v={\frac {{\tilde {k}}\rho }{R}}{\sqrt {(X_{1}-X_{2})^{2}+(Y_{1}-Y_{2})^{2}}}$

Figura 4. Fuente circular donde $P_{0}\$ es $P_{2}\$ y $P\$ es $P_{1}\$ .

Figura 3. Distribución de la irradiancia para una apertura circular

$J_{1}\$ es la función de Bessel[2] de primer orden. Donde $|{\frac {2J_{1}(v)}{v}}|$ decrece constantemente desde el valor 1 cuando $v=0\$, al valor 0 cuando $v=3.83\$; así como los puntos $P_{1}\$ y $P_{2}\$ son separados mas y mas, el grado de coherencia decrece constantemente y hay una completa incoherencia cuando $P_{1}\$ y $P_{2}\$ son separados por una distancia

$P_{1}P_{2}={\sqrt {(X_{1}-X_{2})^{2}+(Y_{1}-Y_{2})^{2}}}={\frac {0.61R{\tilde {\lambda }}}{\rho }}$

Un incremento adicional en v reintroduce una cantidad pequeña de coherencia, pero el grado de coherencia permanece menor a 0.14, y hay una completa incoherencia para $v=7.02\$. Puesto que $J_{1}(v)\$ cambia de signo cada vez que v pasa a través de cada 0 de $J_{1}(v)\$ la fase $B_{12}\$ es equivalente al argumento de $j_{12}\$ cambia por $\pi \$, en consecuencia la posición de las franjas oscuras y brillantes son intercambiadas después de cada desaparición de las franjas.

$P_{1}P_{2}={\frac {0.16R{\tilde {\lambda }}}{\rho }}$

Entonces ${\tilde {\Gamma }}_{12}(0)$ es sencillamente una función sinc cuando la fuente es una rendija y una función Bessel cuando es circular.

Figura 6.Fuente rectangular donde $P_{0}\$ es $P_{1}\$ y $P\$ es $P_{2}\$ en nuestras ecuaciones

Figura 5. Distribución de la irradiancia para una apertura cuadrada

$\beta ={\frac {kb}{2}}sen\theta$

y

$\theta \thickapprox sen\theta .$

$\beta ={\frac {kb\theta }{2}}$

y

$\theta ={\frac {y}{l}}$

por lo tanto

$\left|{\tilde {\gamma }}_{12}(0)\right|=\left|sinc({\frac {\pi by}{l\lambda }})\right|.$

Otro ejemplo es el de la luz solar.

Si consideramos al sol como un disco observado sobre la superficie de la tierra, el diámetro angular 2α que se genera con la proyección del disco solar es aproximadamente de 0 grados 32 min (0.0093 Rad).

Entonces, si despreciamos la variación de la luminosidad del disco solar, esto es si consideramos que no se pueden detectar los cambios en la luminosidad de la fuente, el diámetro del área de coherencia es aproximadamente $0.16{\tilde {\lambda }}\$ (tomando la media de la longitud de onda) que es el resultado que obtuvimos para la fuente circular, es decir los conceptos que utilizamos y el razonamiento serán aplicados de la misma manera para el disco solar. Tomando la longitud de onda promedio como 0.000055 cm, esto proporciona un diámetro de aproximadamente 0.019 mm.

PRINCIPIO DE FERMAT

Cualquiera pensaría que el camino mas corto que seguiría la luz de un punto $P_{1}$ a un punto $P2$ en una superficie reflectora seria una trayectoria geométrica, pero no es así la trayectoria que sigue la luz es un camino que en óptica se conoce como un camino óptico. Son puntos de vista que llega mediante el principio de Fermat.

Fermat reformuló la afirmación de Hero de la siguiente manera: la trayectoria real que adopta la luz entre dos puntos es aquella recorrida en el tiempo mas mínimo.

El principio de Fermat puede enunciarse como sigue: La trayectoria real (rayo) que sigue la luz para viajar desde un punto $P_{1}$ hasta otro punto $P_{2}$ es aquella por la cual el camino óptico D entre ambos puntos es estacionario es decir

$\delta C=\delta \int _{P_{1}}^{P_{2}}nds=0$...........(1)

donde el símbolo $\delta$ significa variacion y s es el parametro del arco sobre la curva considerada. El camino óptico entre dos puntos a lo largo de un rayo puede ser maximo o minimo o mantenerse constante con respecto a la trayectoria proximas a él.

llamaremos $\mathbf {r}$ al vector de posicion de un cierto punto de la curva que representa el rayo que une los puntos $P_{1}$ y $P_{2}$ y sea $\mathbf {R} _{\epsilon }$ el vector de posición asociada a una familia de curvas próximas al rayo

$\mathbf {R} _{\varepsilon }=\mathbf {r} +\mathbf {\varepsilon } \mathbf {f}$...........(2)

donde ${\varepsilon }$ es un numero real y f un vector (arbitraio pero fijo) funcion del punto tal que

$\mathbf {f} (P_{1})=\mathbf {f} (P_{2})=0$............(3)

se tiene por lo tanto, $R_{\varepsilon }(P_{i})=\mathbf {r} (P_{i}),i=1,2$ El camino optico entre $P_{1}$ y $P_{2}$ a lo largo de una cualquiera de las curvas de la familia anterior viene dado por: ${[P_{1}P_{2}]}_{R_{\varepsilon }}=\int _{P_{1}}^{P_{2}}n(\mathbf {R} _{\varepsilon })ds=D({\varepsilon })=0$............(4)

$\delta C=\delta \int _{P_{1}}^{P_{2}}nds=0$

donde la función D(ε) indica explícitamente que $[P_{1}P_{2}]$ depende únicamente de ε (fijados f y el rayo). La condicion 1.1 para D(ε) es equivalente, por tanto, a:

${\frac {dD}{d{\varepsilon }}}=0$.............5

Demostraremos en lo que sigue que esta condicion se satisface para ${\varepsilon }=0$, es decir cuando el camino óptico se calcula sobre el rayo. Si, en notación paramétrica, ${\mathbf {R} _{\varepsilon }}(\alpha )=(x(\alpha ),y(\alpha ),z(\alpha )$, se tiene entonces:

$||{\dot {\mathbf {R} _{\varepsilon }}}||^{2}=({\frac {dx}{d\alpha }})^{2}+({\frac {dy}{d\alpha }})+({\frac {dz}{d\alpha }})^{2}$.............(6)

y sobre la curva, se cumple además

$ds=||{\dot {\mathbf {R} _{\varepsilon }}}||d\alpha$.............(7)

donde, de nuevo, el punto indica derivación con respecto a α (se esta suponiendo que $\mathbf {f}$ admite derivada primera continua). La ecuacion (4) se convierte, por tanto.

$D(\varepsilon )=\int _{\alpha _{1}}^{\alpha _{2}}n(\mathbf {R} _{\varepsilon })||{\dot {\mathbf {R} _{\varepsilon }}}||d\alpha$.............(8)

donde $\alpha _{1}$ y $\alpha _{2}$ son los valores que el parametro α alcanza $P_{1}$ y $P_{2}$ respectivamente. La derivada del camino óptico con respecto a $\varepsilon$ es, por tanto:

${\frac {dD}{d\varepsilon }}=\int _{\alpha _{1}}^{\alpha _{2}}{\frac {d}{d\varepsilon }}[n(\mathbf {R} _{\varepsilon })||{\dot {\mathbf {R} _{\varepsilon }}}||]d\alpha$.............(9)

Ahora bien, sobre la familia de curvas dada (2) se tiene

${\frac {d}{d\varepsilon }}(n||{\dot {\mathbf {R} _{\varepsilon }}}||)=||{\dot {\mathbf {R} _{\varepsilon }}}||\nabla n\bullet f+n{\frac {{\dot {\mathbf {R} }}_{\varepsilon }\bullet {\dot {\mathbf {f} }}}{||{\dot {\mathbf {R} _{\varepsilon }}}||}}$.............(10)

de donde

${\frac {dD}{d\varepsilon }}=\int _{\alpha _{1}}^{\alpha _{2}}||{\dot {\mathbf {R} _{\varepsilon }}}||(\nabla n\bullet f)d\alpha +\int _{\alpha _{1}}^{\alpha _{2}}n{\frac {{\dot {\mathbf {R} }}_{\varepsilon }\bullet {\dot {\mathbf {f} }}}{||{\dot {\mathbf {R} _{\varepsilon }}}||}}d\alpha$.............(11)

intengrando por partes la segunda integral de la ecuación anterior y teniendo en cuenta (3) se llega a:

${\frac {dD}{d\varepsilon }}=\int _{\alpha _{1}}^{\alpha _{2}}[||{\dot {\mathbf {R} _{\varepsilon }}}||\nabla n-{\frac {d}{d\alpha }}[n{\frac {{\dot {\mathbf {R} }}_{\varepsilon }}{||{\dot {\mathbf {R} _{\varepsilon }}}||}}]\bullet \mathbf {f} d\alpha$.............(12)

en particular, si se hace $\varepsilon =0$, es decir, se considera estrictamente el rayo que une a $P_{1}$ y $P_{2}$, la ecuacion (12) toma la forma

${\frac {dD}{d\varepsilon }}|_{\epsilon \to 0}=\int _{\alpha _{1}}^{\alpha _{2}}[||{\dot {\mathbf {r} }}||\nabla n-{\frac {d}{d\alpha }}(n{\frac {\dot {\mathbf {r} }}{||{\dot {\mathbf {r} }}||}})]\bullet \mathbf {f} d\alpha$.............(13)

donde $\mathbf {r} (\alpha )$ satisface la ecuación de los rayos. Recuerdese que, en notación paramétrica, y esta dada por:

${\frac {d}{d\alpha }}(n{\frac {\dot {\mathbf {r} }}{||{\dot {\mathbf {r} }}||}})=||{\dot {\mathbf {r} }}||\nabla n$.............(14)

sustituyendo (14) en (13) se obtiene finalmete

formula${\frac {dD}{d\epsilon }}=0$.............(15)

como queriamos demostrar. Observece que este resultado es independiente del vector f elegido. Si se considera ahora la demostracion en sentido inverso, es ecir partiendo del principio de Fermat, puede definirse un rayo luminos como aquella curva que hace estacionario el camino optico entre dos punto. Imponiendo entoces la condicion (15) a D(ε) definido en (4) para una familia arbitraria de curvas, la demostracion anterio conduce a la igualdad.

$\int _{\alpha _{1}}^{\alpha _{2}}[||{\dot {\mathbf {r} }}||\nabla n-{\frac {d}{d\alpha }}(n{\frac {\dot {\mathbf {r} }}{||{\dot {\mathbf {r} }}||}})]\bullet \mathbf {f} d\alpha =0$.............(16)

--Yonatan Aldana Saldañas 14:38 29 mar 2012 (UTC)

Mas sobre el teorema: [[3]] [[4]] [[5]]

Sobre Fourier Optica: Optica de Fourier

REFERENCIAS.

1. Óptica, autor: Eugene Hecht, Tercera edición, Addison Wesley Iberoamericana, Madrid, 2000.

2. Principles of optica, autores: Max Born y Emil Wolf, Quinta edición 1975, Pergamon Press.