# Principio de Fermat

El principio de Fermat, en óptica, es un principio de tipo extremal y que establece:
Este enunciado no es completo y no cubre todos los casos, por lo que existe una forma moderna del principio de Fermat. Esta dice que:
Esto quiere decir que, si se expresa el trayecto recorrido por la luz entre dos puntos $O_{1}$ y $O_{2}$ por medio de una funcional llamada camino óptico definida como ${\mathcal {L}}_{O_{1}O_{2}}[n({\vec {r}})]$ la trayectoria real de la luz seguirá un camino extremal respecto de esta funcional:
   $\delta {\mathcal {L}}_{O_{1}O_{2}}[n({\vec {r}})]=\delta \int _{O_{1}}^{O_{2}}{n({\vec {r}})ds}=0.$
La característica importante, como dice el enunciado, es que los trayectos próximos al verdadero requieren tiempos aproximadamente iguales. En esta forma, el principio de Fermat recuerda al principio de Hamilton o a las ecuaciones de Euler-Lagrange.
El principio en su forma moderna fue declarado por Pierre de Fermat en una carta de 1662, de ahí que lleve su nombre.
Siguen ahora algunos ejemplos de la aplicación del principio para deducir las leyes de la óptica geométrica.

## Ecuación de la trayectoria de un rayo luminoso

La ecuación de la trayectoria de un rayo luminoso real en un sistema óptico es:
${\vec {\nabla }}n({\textbf {r}})-{d \over ds}\left[n({\textbf {r}}){{\textbf {dr}} \over {ds}}\right]=0$
y se deduce a partir del principio de Fermat.

## Deducción de la ecuación

Tenemos un rayo de luz que se desplaza por un medio con índice de refracción continuo desde un punto $A$ a un punto $B$.
Sea ${\textbf {r}}$ el vector posición y ${\textbf {u}}$ el vector unitario tangente a la trayectoria.
Tenemos que $L=\int _{A}^{B}n({\textbf {r}})ds=\int dL$. De esto se deduce que $nds=dL$. Además se tiene que $ds={\textbf {u}}\cdot {\textbf {dr}}$ y que $dL={\vec {\nabla }}L\cdot {\textbf {dr}}$.
Por lo tanto $n{\textbf {u}}\cdot {\textbf {dr}}={\vec {\nabla }}L\cdot {\textbf {dr}}$. Al ser ${\textbf {dr}}$ de un vector arbitrario tenemos que $n{\textbf {u}}={\vec {\nabla }}L$ y por tanto que $\mid {{\vec {\nabla }}L}\mid =n$.
De esto se obtiene que ${{d(n{\textbf {u}})} \over {ds}}={d \over ds}{\vec {\nabla }}L$. Notamos que${{d(n{\textbf {u}})} \over {ds}}={d \over ds}\left[n({\textbf {r}}){{\textbf {dr}} \over ds}\right]$.
Por otro lado tenemos que (se empleará cartesianas pero sirve para el resto de bases ortonormales):
${d \over ds}\left({{\partial L} \over {\partial x}}\right)={\vec {\nabla }}\left({{\partial L} \over {\partial x}}\right)\cdot {{\textbf {dr}} \over ds}={\vec {\nabla }}\left({{\partial L} \over {\partial x}}\right)\cdot {\textbf {u}}$; porque ${d}\left({{\partial L} \over {\partial x}}\right)={\vec {\nabla }}\left({{\partial L} \over {\partial x}}\right)\cdot {\textbf {dr}}$.
${d \over ds}\left({{\partial L} \over {\partial x}}\right)={\vec {\nabla }}\left({{\partial L} \over {\partial x}}\right)\cdot {1 \over n}{\vec {\nabla }}L$; porque $n{\textbf {u}}={\vec {\nabla }}L$.
${d \over ds}\left({{\partial L} \over {\partial x}}\right)={1 \over n}\left({{\partial } \over {\partial x}}{{\partial L} \over {\partial x}},{{\partial } \over {\partial y}}{{\partial L} \over {\partial x}},{{\partial } \over {\partial z}}{{\partial L} \over {\partial x}}\right)\cdot \left({{\partial L} \over {\partial x}},{{\partial L} \over {\partial y}},{{\partial L} \over {\partial z}}\right)={1 \over n}\left[{{\partial } \over {\partial x}}{{\partial L} \over {\partial x}}{{\partial L} \over {\partial x}}+{{\partial } \over {\partial y}}{{\partial L} \over {\partial x}}{{\partial L} \over {\partial y}}+{{\partial } \over {\partial z}}{{\partial L} \over {\partial x}}{{\partial L} \over {\partial z}}\right]$.
${d \over ds}\left({{\partial L} \over {\partial x}}\right)={1 \over n}\left[{{\partial } \over {\partial x}}{{\partial L} \over {\partial x}}{{\partial L} \over {\partial x}}+{{\partial } \over {\partial x}}{{\partial L} \over {\partial y}}{{\partial L} \over {\partial y}}+{{\partial } \over {\partial x}}{{\partial L} \over {\partial z}}{{\partial L} \over {\partial z}}\right]={1 \over {2n}}{{\partial } \over {\partial x}}\left[\left({{\partial L} \over {\partial x}}\right)^{2}+\left({{\partial L} \over {\partial y}}\right)^{2}+\left({{\partial L} \over {\partial z}}\right)^{2}\right]$; porque ${{\partial } \over {\partial \beta }}\left({{\partial L} \over {\partial \alpha }}\right)^{2}=2{{\partial } \over {\partial \beta }}{{\partial L} \over {\partial \alpha }}{{\partial L} \over {\partial \alpha }}$.
${d \over ds}\left({{\partial L} \over {\partial x}}\right)={1 \over {2n}}{{\partial } \over {\partial x}}n^{2}={1 \over {2n}}2n{{\partial n} \over {\partial x}}={{\partial n} \over {\partial x}}$; porque$\left[\left({{\partial L} \over {\partial x}}\right)^{2}+\left({{\partial L} \over {\partial y}}\right)^{2}+\left({{\partial L} \over {\partial z}}\right)^{2}\right]=\mid {\vec {\nabla }}L\mid ^{2}=n^{2}$.
Realizando las mismas componentes para $y$ e $z$ se obtiene que:
${d \over ds}\left({\vec {\nabla }}L\right)={d \over ds}\left({{\partial L} \over {\partial x}}\right){\textbf {i}}+{d \over ds}\left({{\partial L} \over {\partial y}}\right){\textbf {j}}+{d \over ds}\left({{\partial L} \over {\partial z}}\right){\textbf {k}}={{\partial n} \over {\partial x}}{\textbf {i}}+{{\partial n} \over {\partial y}}{\textbf {j}}+{{\partial n} \over {\partial z}}{\textbf {k}}={\vec {\nabla }}n$.
Reemplazando se obtiene que:
${{d(n{\textbf {u}})} \over {ds}}={\vec {\nabla }}n$ o ${d \over ds}\left[n({\textbf {r}}){{\textbf {dr}} \over ds}\right]={\vec {\nabla }}n$.

## Analogía con la mecánica lagrangiana

El camino óptico se puede equiparar a la acción en la mecánica lagrangiana. Se puede tratar el índice de refracción como un lagrangiano compuesto por un potencial. De este modo el problema se puede resolver con las ecuaciones de Lagrange. El rayo de luz se dirige hacia la zona de mayor índice de refracción, de lo que el potencial equivalente sería el opuesto del índice de refracción.
Tenemos: $\int _{A}^{B}n(\lbrace q_{i}\rbrace )ds=\int _{A}^{B}n(\lbrace q_{i}\rbrace )Rd\lambda \equiv \int _{A}^{B}L(\lbrace q_{i}\rbrace )d\lambda$ .
Con: $L(\lbrace q_{i}\rbrace )=n(\lbrace q_{i}\rbrace )R$ , además $R={\sqrt {\sum \left({{dq_{i}} \over {d\lambda }}\right)^{2}}}={\sqrt {\sum {{\dot {q}}_{i}}^{2}}}$ y $d\lambda ={ds \over R}$ .
De modo que si definimos $q_{i}(\lambda )$ obtenemos que: ${d \over {d\lambda }}\left({{\partial (nR)} \over {\partial {\dot {q}}_{i}}}\right)-{{\partial (nR)} \over {\partial q_{i}}}=0$.
Se obtiene: $R{d \over {ds}}\left({{\partial (n)} \over {\partial {\dot {q}}_{i}}}R+{{\partial (R)} \over {\partial {\dot {q}}_{i}}}n\right)-R{{\partial (n)} \over {\partial q_{i}}}=0$ con ${d \over {d\lambda }}={1 \over {1/R}}{d \over {ds}}=R{d \over ds}$
Además se tiene que ${{\partial (n)} \over {\partial {\dot {q}}_{i}}}=0$ y que ${{\partial R} \over {\partial {\dot {q}}_{i}}}={{2{\dot {q}}_{i}} \over {2R}}={{{\dot {q}}_{i}} \over {R}}={1 \over R}{dq_{i} \over d\lambda }={1 \over R}R{dq_{i} \over ds}={dq_{i} \over ds}$.
De donde: $R{d \over {ds}}\left(n{dq_{i} \over ds}\right)-R{{\partial (n)} \over {\partial q_{i}}}=0$ o ${d \over {ds}}\left(n{dq_{i} \over ds}\right)-{{\partial (n)} \over {\partial q_{i}}}=0$.
Creando un sistema con las tres coordenas se obtiene que: ${d \over {ds}}\left(n{{\textbf {dr}} \over ds}\right)-{\vec {\nabla }}n=0$.

## Ecuación de la trayectoria de un rayo luminoso en un medio homogéneo

Partimos de que:
${d \over ds}\left[n({\textbf {r}}){{\textbf {dr}} \over ds}\right]={\vec {\nabla }}n$.
Con $n$constante, de lo que:
${d \over ds}\left[n({\textbf {r}}){{\textbf {dr}} \over ds}\right]=n{d \over ds}{{\textbf {dr}} \over ds}$y ${\vec {\nabla }}n=0$.
Por tanto:
${d \over ds}{{\textbf {dr}} \over ds}=0$ de lo que ${{\textbf {dr}} \over ds}={\textbf {l}}$ siendo ${\textbf {l}}$ una constante, de lo que ${\textbf {r}}={\textbf {l}}s+C$ siendo ${\textbf {C}}$ una constante.
El resultado es una recta de punto inicial ${\textbf {C}}$ y vector director ${\textbf {l}}$.

## Teorema de Malus-Dupin

Si sobre cada rayo emitido por un foco recorremos caminos ópticos iguales, entonces los puntos que los delimitan forman una superficie normal a todos los rayos. Denominamos a dicha superficie frente de ondas. Coincide con el frente de onda dado por la teoría oscilatoria. Al deducirse del principio de Fermat es válido a pesar del número de reflexiones o refracciones que pueda sufrir el rayo antes de llegar a su destino.

## Ley de la reflexión

Si se supone que un rayo de luz sale del punto A en dirección a la superficie plana, que suponemos reflectora, y viaja hasta el punto B ¿Cuál será la trayectoria seguida por la luz? En este caso la luz viaja durante todo el camino por el mismo medio, con el mismo índice de refracción y, por tanto, a la misma velocidad. Así, el tiempo necesario para recorrer el camino entre A y B (pasando por la superficie P) será la distancia APB dividida por la velocidad de la luz en ese medio. Como la velocidad es una constante, la trayectoria real, según el principio de Fermat, será la más corta.
Es fácil ver que la distancia APB es la misma que la distancia A'PB, donde A' es la imagen de A. A' está sobre la recta perpendicular al espejo que pasa por A, a la misma distancia del espejo que A y al otro lado del mismo. La distancia mínima A'PB es, obviamente, la línea recta A'P2B, con lo que la trayectoria real es AP2B. El análisis completo de la situación muestra que P2 es tal que los ángulos de incidencia y de reflexión en el punto son iguales, de lo que se deduce la fórmula de la ley de la reflexión: $\theta _{i}=\theta _{r}$

## Ley de la refracción

Sea un medio de propagación con índice de refracción $n_{1}\$ y un segundo medio de propagación con índice de refracción $n_{2}\$ tales que situamos la superficie que separa los dos medios de modo que coincida con el eje de las abcisas.
Sean $A=(x_{A},\;y_{A})$ y $B=(x_{B},\;y_{B})$ dos puntos fijos situados del plano, de modo que A está situado en el primer medio, y B en el segundo medio.
Sea un rayo de luz que se propaga de A a B atravesando la superficie que separa los dos medios en el punto $P=(x,\;0)$.
El siguiente paso es deducir el tiempo que tarda el rayo en recorrer ${\overline {AP}}$ y ${\overline {PB}}$.
Sean $v_{1}$ y $v_{2}$ la velocidad de propagación de la luz en el primer y segundo medio respectivamente.
   $t_{1}={\frac {\overline {AP}}{v_{1}}}={\frac {\sqrt {(x_{A}\ -x)^{2}\ +{y_{A}}^{2}}}{v_{1}}}$; $t_{2}={\frac {\overline {PB}}{v_{2}}}={\frac {\sqrt {(x-x_{B})^{2}\ +{y_{B}}^{2}}}{v_{2}}}$
   $t={\frac {\sqrt {(x_{A}\ -x)^{2}\ +{y_{A}}^{2}}}{v_{1}}}+{\frac {\sqrt {(x-x_{B})^{2}\ +{y_{B}}^{2}}}{v_{2}}}$
Si se busca el valor de $x\$ cuando $t\$ es mínimo, es equivalente si encontramos el valor de $x\$ para el cual la función derivada de $t\$ toma el valor 0.
   ${\frac {dt}{dx}}=-{\frac {x_{A}\ -x}{v_{1}\ {\sqrt {(x_{A}\ -x)^{2}\ +{y_{A}}^{2}}}}}+{\frac {x-x_{B}\ }{v_{2}\ {\sqrt {(x-x_{B}\ )^{2}\ +{y_{B}}^{2}}}}}=0$
   ${\frac {x_{A}\ -x}{v_{1}\ {\sqrt {(x_{A}\ -x)^{2}\ +{y_{A}}^{2}}}}}={\frac {x_{B}-x\ }{v_{2}\ {\sqrt {(x-x_{B}\ )^{2}\ +{y_{B}}^{2}}}}}$
   ${\frac {x_{A}\ -x}{v_{1}\ {\overline {AP}}}}={\frac {x_{B}-x\ }{v_{2}\ {\overline {PB}}}}$
   ${\frac {1}{v_{1}\ }}\sin {\alpha _{1}}={\frac {1}{v_{2}\ }}\sin {\alpha _{2}}$
   ${\frac {c}{v_{1}\ }}\sin {\alpha _{1}}={\frac {c}{v_{2}\ }}\sin {\alpha _{2}}$
   $n_{1}\ \sin {\alpha _{1}}=n_{2}\ \sin {\alpha _{2}}$

## Historia

Herón de Alejandría (Heron) (c. 60) describió un principio de reflexión, que declaraba que un rayo de luz que va desde el punto A al punto B, sufriendo cualquier número de reflexiones en espejos planos, en el mismo medio, tiene un trayecto con menor longitud que cualquier trayecto cercano.
Ibn al-Haytham (Alhazen), en su Libro de Óptica (1021), amplió el principio tanto a la reflexión como a la refracción, y expresó una temprana versión del principio del menor tiempo. Sus experimentos se basaron en trabajos anteriores sobre la refracción realizados por el científico griego Claudio Ptolomeo.
El principio generalizado del menor tiempo en su forma moderna fue declarado por Pierre de Fermat en una carta fechada el 1 de enero de 1662 enviada a Cureau de la Chambre. Se encontró con las objeciones efectuadas en mayo de 1662 por Claude Clerselier, un experto en óptica y líder portavoz de los cartesianos en ese momento. Entre sus objeciones, Clerselier establecía:
El original, en francés, de Mahoney, es el siguiente:
De hecho el principio de Fermat no se sostiene por sí solo, y ahora se sabe que se puede derivar de principios anteriores, como el principio de Huygens. Históricamente, el principio de Fermat ha servido como principio rector en la formulación de las leyes de la Física con el uso del cálculo variacional (véase el principio de mínima acción).