AM

Capı́tulo 3

UN
Optica con Matrices.

SS
Durante el curso hemos visto parte importante de la óptica geométrica
principalmente dentro de la aproximación paraxial. Con las leyes de reexión
y de refracción es posible calcular la trayectoria que recorre la luz al pasar
-D
por componentes ópticos como lentes gruesas, lentes delgadas y supercies
esféricas. En la gura 3.1 mostramos el ejemplo sencillo de un arreglo de lentes
que puede ser utilizado para formar imágenes de objetos sobre el sensor de
una cámara, por ejemplo. Con la óptica de matrices es posible predecir como
se formará esta imagen si uno cuenta con las caracterı́sticas de cada lente.
M

Una vez hecho este cálculo uno puede tratar al arreglo óptico como una caja
negra que arroja un resultado según las caracterı́sticas de la luz que entra.
La gran utilidad de esto es el poder construir óptica útil, para lo cual es
JJ

necesario jar un objetivo y calcular las caracterı́sticas necesarias de cada

elemento para lograr un n determinado.
Es indispensable poder predecir el efecto que tiene cada material sobre la
trayectoria de la luz puesto que cambia al cruzar cada interfaz. El método
ica

que hemos estudiado es suciente pero poco práctico para arreglos de dos
o más elementos ópticos. Por ejemplo, si la posición del objeto cambia, hay
que rehacer una serie de cálculos para calcular su imagen, lo que requiere
pt

cuidado y concentración. En este capı́tulo estudiaremos la óptica Gaussiana

que permite hacer estos cálculos de una manera más directa y eciente.
Como muestra la gura 3.2 cada vez que la luz atraviesa la supercie de
O

separación entre dos medios tiene un cambio de trayectoria. Caracterizare-

mos este cambio por dos variables usando como referencia el eje óptico del
sistema: el ángulo ↵ que hace el rayo, y la altura y del punto en el que cruza
la supercie, ambos medidos respecto del eje óptico. Esencialmente todos los

81
 82 CAPÍTULO 3. OPTICA CON MATRICES.

AM
Figura 3.1: El objetivo de la óptica con matrices es simplicar el tratamiento
matemático de cualquier arreglo de manera tal que con la información de
entrada sea suciente para saber la trayectoria que seguirá la luz a la salida,
como si fuese una caja negra.

UN
SS
-D
M

Figura 3.2: Cambio de trayectoria de un rayo al atravesar la supercie de

separación de dos medios.
JJ

componentes ópticos cuentan con una supercie plana o esférica. Ası́ que a
estas alturas del curso ya podemos planear sistemas ópticos relativamente
complejos. Estos, pueden ser imaginados como una serie de materiales subse-
ica

cuentes, cada uno con ı́ndice de refracción constante como muestra la Figura
??. Hemos visto que se trata de una tarea sencilla. Sin embargo calcular a
mano la trayectoria de la luz considerando varios elementos ópticos puede
pt

llevar largo tiempo y existe una probabilidad alta de error.

Nos hace falta un forma eciente de computar sistemas compuestos por
un número arbitrario de supercies. Un método estándar utilizado tanto en
O

investigación e industria es conocido como óptica de matrices. Fue inventado

en los años 30 por el Inglés Thomas Smith [1883-1969]. Consiste en reducir
cada cambio de dirección causado por un cambio de ı́ndice de refracción a un
sistema de dos ecuaciones con dos incógnitas. Una vez encontrado el sistema
 83

AM
UN
Figura 3.3: Sistema óptico formado por varios medios separados por super-
cies. SS
-D
de ecuaciones es trivial encontrar la matriz que nos permite resolverlo:
✓ ◆ ✓ ◆
y y
=M . (3.1)
n2 ↵ 2 n1 ↵ 1
Ası́, aunque tengamos un número arbitrario de componentes ópticos, siempre
M

podemos encontrar la matriz equivalente al sistema completo. Por ejemplo,

si tenemos 5 materiales distintos, como muestra la gura 3.3, la matriz equi-
valente está dada por el producto de todas las matrices empezando por el
JJ

último componente, i. e.
M = M5 M 4 M 3 M 2 M 1 . (3.2)
para el sistema que muestra la gura 3.3. Con este método es trivial calcular
ica

sistemas ópticos complejos llevando acabo solamente operaciones de aritméti-

ca y de álgebra lineal. Lo cual es trivial para cualquier computadora casera
en nuestros dı́as. Son pocas las matrices necesarias para modelar arreglos
pt

ópticos complejos. En este capı́tulo nos concentraremos en describir lo que

sucede a través de un medio homogéneo, una interfaz plana y una esférica.
Con ellas podremos construir la matriz equivalente para una lente gruesa y
O

después deducir la matriz que describe a una lente delgada. Veremos que,
bajo ciertas condiciones, es posible encontrar la lente delgada equivalente a
un arreglo óptico que podrı́a ser más complejo. Lo cual nos permite diseñar
artefactos ópticos más complejos.
 84 CAPÍTULO 3. OPTICA CON MATRICES.

AM
UN
Figura 3.4: Un rayo luz cruzando un medio homogéneo no cambia su trayec-
SS
toria rectilı́nea original. En este caso lo consideramos viajando con un ángulo
↵ con respecto a la horizontal. Lo cual hace tener alturas distintas en cada
momento denotadas por y1 y y2 .
-D
3.1. Matriz para un medio homogéneo y una
supercie esférica
Dentro de un medio con el mismo ı́ndice de refracción y sin obstáculos la
M

luz no tiene razón para desviarse. Por lo que mantendrá una trayectoria en
lı́nea recta. Esto ocurre en el interior de casi todos los elementos ópticos fabri-
JJ

cados actualmente. La gura 3.4 ilustra este caso con un rayo siguiendo una
trayectoria recta caracterizada por un ángulo ↵ que causa un desplazamien-
to vertical de y1 a y2 después de viajar una distancia t sobre la coordenada
horizontal t. Para construir la matriz correspondiente siguiendo la prescrip-
ica

ción general anteriormente mencionada consideraremos explı́citamente que el

ángulo es constante, es decir, ↵1 = ↵2 = ↵. Invitamos al lector a comprobar
que la matriz que describe a esta transformación, siguiendo la forma general
pt

descrita por la ecuación (3.1), es

✓ ◆
1 nt
Mt = . (3.3)
O

0 1

Con la ecuación (3.3) en combinación con la matriz correspondiente a la

transformación causada por una supercie esférica podemos construir a la
matriz equivalente de cualquier lente. Para construir esta última recordemos
 3.1. MATRIZ PARA UN MEDIO HOMOGÉNEO Y UNA SUPERFICIE ESFÉRICA85

AM
UN
SS
Figura 3.5: Cambio de trayectoria de un rayo al atravesar una supercie
-D
esférica.

la geometrı́a vista previamente, donde consideramos a un rayo que viaja

de izquierda a derecha incidiendo sobre una supercie esférica con radio de
curvatura +R (centro localizado a la derecha del vértice), como se muestra
M

en la gura 3.5. Aplicamos la ley de Snell paraxial en el punto donde el rayo

atraviesa la supercie:
n 1 ✓i = n 2 ✓r . (3.4)
JJ

De la gura sabemos que

✓i = ↵ 1 +  (3.5)
ica

✓r = ↵2 + . (3.6)
Por lo tanto
n1
↵ 2 = ✓r   = ( + ↵1 )  , (3.7)
pt

n2
pero  = y/R por lo que
✓ ◆
O

n1 y n2  n1
↵2 = ↵1  . (3.8)
n2 R n2
Si tomamos en cuenta que el punto en el que el rayo toca a la supercie es
común a las dos trayectorias (y2 = y1 = y) y reescribimos la Ecuación 3.8,
 86 CAPÍTULO 3. OPTICA CON MATRICES.

AM
UN
SS
Figura 3.6: Una lente gruesa está compuesta en general por un material
óptico con dos radios de curvatura R1 y R2 distintos. Aunque normalmente
-D
está inmerso en un medio, el dentro del caso general debemos considerar un
ı́ndice de refracción distinto, n1 y n3 , para cada lado. Por lo que, describir
completo este pequeño sistema óptico, necesitamos construir tres matrices
equivalentes: Mn1,n2,R1 , Mn2,t y Mn2,n1,R2 .
M

recordando que PR = (n2  n1 )/R es la dioptrı́a o potencia de una supercie

esférica con radio R,
JJ

y n1
↵2 =  PR + ↵ 1 .
n2 n2
ica

o bien, en forma matricial,

✓ ◆ ✓ ◆✓ ◆
y2 1 0 y
= . (3.9)
pt

n2 ↵ 2 PR 1 n1 ↵ 1

Es importante notar que en este caso la matriz equivalente depende de dos

O

ı́ndices de refracción, a diferencia de la matriz que corresponde a un medio

homogéneo. Como veremos más adelante las matrices consideradas en las
transformaciones (3.3) y (3.9) son el bloque de contracción de prácticamente
cualquier elemento óptico.
 3.2. MATRIZ PARA UNA LENTE GRUESA Y SU LENTE DELGADA EQUIVALENTE87

3.2. Matriz para una lente gruesa y su lente

delgada equivalente

AM
La geometrı́a de una lente gruesa está descrita en la Figura 3.6. Revi-
sando la gura izquierda a derecha observamos que la luz atraviesa por dos
supercies. El cambio a través de su primer interfaz, de n1 a n2 , lo acabamos
de calcular. Este es seguido por una trayectoria recta únicamente a lo largo
del material que compone la lente es decir, una trayectoria a lo largo de un

UN
medio homogéneo, bien descrita por la ecuación 3.3. Aunque normalmente
los elementos ópticos están inmersos en un solo medio, lo cual jarı́a n1 = n3 ,
en la gura 3.6 consideramos el caso general en el que atravesar la segunda
supercie implica también un cambio de medio de refracción. La matriz para
describir esta etapa nal del proceso óptico deberá tomar en cuenta esto y el
cambio de signo en el radio de curvatura. SS
Para construir la matriz completa a una lente gruesa debemos juntar los
resultados resultados para una supercie esférica
-D
✓ ◆
1 0
Mn1 ,n2 ,R = . (3.10)
 n2 n
R
1
1
y para un medio homogéneo
✓ ◆
M

1 t/n
Mn,t = (3.11)
0 1
Dejamos como ejercicio vericar que en los dos casos el determinante es
JJ

igual a 1; el producto de matrices unitarias tendrá determinante unitario

también. Este es un requisito indispensable para que el álgebra de matrices
funcione y para que pueda describir a un elemento óptico sin ambigüedades.
ica

Para terminar de obtener la matriz equivalente a la lente gruesa en cuestión

únicamente nos falta precisar la diferencia entre las matrices correspondientes
a la primer y segunda supercie esférica. Adoptemos la notación de la Figura
3.6. Usaremos R1 y R2 para denotar a los radios de curvatura izquierdo y
pt

derecho respectivamente. Desde ahora notamos que R1 > 0 y R2 < 0. Sin

perder generalidad sustituimos directamente en la Ecuación 3.9 para obtener,
O

por ejemplo, la matriz correspondiente a la supercie con radio de curvatura

R2
✓ 0 ◆ ✓ ◆✓ ◆
y 1 0 y
= , (3.12)
n2 ↵ 0 PR2 1 n1 ↵
 88 CAPÍTULO 3. OPTICA CON MATRICES.

donde
n1  n2
PR2 = . (3.13)
R2
Haciendo el producto de los tres segmentos que componen a la lente gruesa

AM
obtenemos su matriz equivalente

M` = Mn2 ,n3 ,R2 Mn2 ,t Mn1 ,n2 ,R1

✓ ◆✓ ◆✓ ◆
1 0 1 nt2 1 0
=

UN
n3 n2 n2 n1
 R2 1 0 1  R1 1
✓ ◆ !
1 0 1  t(nn22n
R1
1) t
n2
= n3 n2 n n
 R2 1  R1 2 1
1
!
⇣ 1  t(n⌘n22n
R1
1) t
n2
=
 n2 n1
R1
 n3Rn
2
2 SS+ t(n2 n 3 )(n1 n2 )
n 2 R1 R2

Este caso ilustra la forma general en la que se computa cada segmento

1 t(n3 n2 )
n 2 R2
(3.14)
-D
en la propagación de un haz. Dado un arreglo arbitrario de materiales y
fronteras es posible encontrar su matriz equivalente abcd con la que será
posible predecir la trayectoria del rayo trasmitido,
✓ 0 ◆ ✓ ◆ ✓ ◆✓ ◆
y y a b y
= M = . (3.15)
M

E
ni ↵ 0 nt ↵ c d nt ↵

El objetivo es predecir que sucederá con la trayectoria de la luz antes y

JJ

después de un sistema óptico. Igualmente a la inversa: dada una necesidad

en el laboratorio, ¿qué conjunto de componentes debo seleccionar para lograr
la manipulación deseada? De tal manera que podamos reducir el problema a
un bulto a donde metemos luz y la escupe con las caracterı́sticas necesarias.
ica

Para esto es útil identicar formas matriciales recurrentes.

En el Capı́tulo 2 dedujimos la ecuación 2.10 de Gauss para una lente
delgada a partir de hacer tender a cero el grosor t de una lente gruesa. En
pt

esta sección mostraremos que la misma aproximación es válida dentro del

tratamiento matricial. Es decir, dada la matriz de una lente gruesa encon-
traremos una forma general de deducir la matriz que corresponde a la lente
O

delgada equivalente. Posteriormente mostraremos que es posible encontrar

esta matriz incluso en el caso de sistemas ópticos más complejos.
En su lı́mite t ! 0 la matriz (3.14), que es válida para cualquier lente
gruesa, tiende a
 3.3. LENTE DELGADA EQUIVALENTE A CUALQUIER FORMADOR DE IMAGEN89

!
⇣ 1 ⌘ 0
Mdelg. = n2 n1
 R1
 n3Rn 2
1

AM
2
✓ ◆
1 0
= . (3.16)
p 1
En este caso es necesario denir la potencia de la lente en términos de los
dos radios de curvatura y las diferencias entre los ı́ndices de refracción,

UN
✓ ◆
n2  n1 n2  n3
P =  (3.17)
R1 R2
.
Solo en casos muy especializados es necesario considerar n3 distinto de
n1 . Normalmente los lentes son colocados inmersos en un medio, dando pie
SS
a simplicar las ecuaciones (3.14) y (3.16) estableciendo n3 = n1 . Tenemos
que para la lente gruesa
!
-D
1⇣ t(nn22n 1)
R1 ⌘ n2
t
Mt,R1,R2 = 2 n1 )
2 (3.18)
(n2  n1 ) R11  R12  t(n n 2 R1 R2
1 + t(nn22n
R2
1)
.
Para obtener la matriz de la lente equivalente al lente grueso dentro del
mismo medio exterior debemos tomar el lı́mite de la matriz (3.18) cuando
M

t ! 0, como hicimos en el caso general. Con lo cual obtenemos

!
1⇣ ⌘ 0
JJ

Mdelg. =
(n2  n1 ) R11  R12 1
✓ ◆ ✓ ◆
1 0 1 0
= = . (3.19)
 nf1 1 P 1
ica

En la matriz (3.19) hemos denido a las distancia focal f de acuerdo a la

ecuación (2.11), que describe el funcionamiento de una lente delgada y dene
su distancia focal f .
pt

3.3. Lente delgada equivalente a cualquier for-

O

mador de imagen
Hemos visto que para la lente delgada y la supercie esférica es posible
encontrar la trayectoria del rayo de salida si conocemos sus caracterı́sticas
 90 CAPÍTULO 3. OPTICA CON MATRICES.

H H’

AM
D D’

UN
Figura 3.7: Para encontrar la lente delgada equivalente a un sistema óptico
cualquiera (negro) lo extendemos usando dos rebanadas con grosor D y D0
(lila). La equivalencia existirá siempre y cuando las ecuaciones que relacionan
SS
a los dos sistemas ópticos tenga solución para un par de valores D y D0 .
-D
geométricas de entrada. En ambos casos la forma de pensar es la misma:
obtenemos a los nuevos parámetros mediante una transformación dada por
la matriz (3.19). En otras palabras, justo después de atravesar la lente en-
contramos que el vector
✓ 0 ◆ ✓ ◆✓ ◆
y 1 0 y
M

=
n0 ↵ 0 P 1 n↵
✓ ◆
y
JJ

= (3.20)
P y + n↵

describe su geometrı́a de propagación. Como la lente es delgada la altura y

no cambia inmediatamente; el ángulo ↵ si lo hace, de acuerdo a
ica

n P
↵0 = ↵  0 y. (3.21)
n0 n
pt

A continuación veremos que existen ciertas condiciones bajo las cuales es

posible generalizar la relación para encontrar la lente delgada (o supercie
esférica) equivalente de una lente gruesa a un sistema óptico arbitrario. De
O

inicio solo supondremos que este sistema, ilustrado como una caja negra en
la gura 3.7, está descrito por una matriz abcd cuyo determinante es igual a
uno. Comencemos por extender al sistema arbitrario con dos rebanadas de
propagación homogénea. Estas extensiones tienen grosores D y D0 , aparecen
 3.3. LENTE DELGADA EQUIVALENTE A CUALQUIER FORMADOR DE IMAGEN91

en color lila la gura 3.7. Traducido al lenguaje de matrices esto es,

✓ D0 ◆ ✓ ◆✓ ◆
1 n1 a b 1  nD1
MEquiv. =
0 1 c d 0 1

AM
✓ ◆ ✓ aD ◆
1 t0 /n a  n1 + b
=
0 1 c  cDn1
+d
0
⇣ ⌘!
D0
a + cDn
 aD
n1
+ b + n1
 cD
n1
+ d
= . (3.22)

UN
c  cD
n1
+ d

Para encontrar la matriz de la lente delgada equivalente a nuestra caja negra

debemos igualar entrada por entrada las matrices (3.16) y (3.22). Con lo cual
arribamos al siguiente sistema de ecuaciones

✓
cD0
a+
n1
◆
=1
SS
aD D0 cD
-D
 +b+  +d =0
n1 n1 n
1
c = P (= )
f
cD
d = 1. (3.23)
M

n1
En primer lugar notamos que la potencia del sistema está dada por c = 1/f .
JJ

Ahora a, b, c y d están dadas por las caracterı́sticas del sistema óptico pero
podemos escoger D y D0 de tal manera que el sistema de ecuaciones anterior
sea soluble. De la primera y última ecuación del sistema (3.23) encontramos
las condiciones
ica

1a a1
D0 = n1 = n1
c P
y
pt

1d d1
D= n1 =  n1 . (3.24)
c P
O

Las ecuaciones (3.24) nos dicen el grosor que debe tener cada rebanada de
extensión al sistema. Los planos perpendiculares que cruzan al eje óptico
en los puntos D y D0 se llaman planos principales, denotados por H y H 0
en la gura 3.7. En las ecuaciones (3.24) encontramos que c 6= 0 es una
 92 CAPÍTULO 3. OPTICA CON MATRICES.

AM
Figura 3.8: Cualquier sistema óptico formador de imágenes puede ser visto
como una lente delgada equivalente. Esta debe cumplir con la ley de Gauss

UN
que relaciona su distancia focal f con la distancia del respectivo plano prin-
cipal al objeto So y a la imagen Si . Como es de esperarse la imagen es sujeta
a una amplicación MT que relaciona su tamaño y orientación con las pro-
piedades geométricas del objeto.

SS
condición necesaria para que exista la equivalencia entre nuestra caja negra
óptica y una lente delgada. Esto tiene sentido puesto que el sistema arbitrario
debe ser un formador de imagen; los planos principales son la referencia que
-D
usaremos para ubicar al objeto y a la imagen en el marco de la lente delgada
equivalente.
Ahora, si sustituimos D y D0 en la tercer ecuación del sistema (3.23)
encontramos que
✓ ◆
M

aD D0 cD a(1  d) d(1  a) (1  a)(1  d)

 +b+  +d = +b+ +
n1 n1 n1 c c c
JJ

a(1  d) + bc + d(1  a) + (1  a)(1  d)

=
c
a  da + bc + d  da + 1  d  a + da
=
c
ica

1  (da  bc)
= = 0. (3.25)
c
Con la ecuación (3.25) notamos que nuestra hipótesis de inicio, det(MEquiv. ) =
pt

bc  da = 1, es la segunda condición necesaria para que sea posible encontrar

una lente delgada equivalente a cualquier sistema óptico descrito por una
matriz abcd.
O

Recordemos que estamos hablando de un sistema óptico creador de imáge-

nes, visualizado como una la caja negra en la gura 3.7. Por medio del pro-
cedimiento anterior demostramos que todo elemento que pertenezca a esta
categorı́a tiene una lente delgada equivalente. De tal manera que el sistema
 3.3. LENTE DELGADA EQUIVALENTE A CUALQUIER FORMADOR DE IMAGEN93

AM
UN
SS
Figura 3.9: La matriz de transferencia cubre la trayectoria de la luz entre
-D
dos puntos determinados. En el caso de un sistema generador de imagen lo
útil es que abarque desde la ubicación So1 del objeto hasta el punto Si5 de la
imagen. Los subı́ndices indican que en este caso las distancias son medidas
con respecto al vértice del primer objeto y con respecto al vértice del último
elemento óptico, el cinco en este caso.
M

mostrado en las guras 3.7 y 3.8 son equivalentes. Sin embargo, debemos ser
JJ

precavidos al calcular las propiedades de la imagen resultante. Normalmente

sabremos la distancia So1 del objeto al vértice anterior del primer elemen-
to óptico del sistema. Para conocer So debemos considerar la distancia D
entre este vértice y el plano principal correspondiente i. e., So = So1  D.
ica

Una vez conocida So podemos aplicar directamente la ecuación de Gauss

para encontrar Si . De esta manera se vuelve trivial calcular la amplica-
ción MT = Si /So . De igual manera, en términos de la distancia del último
vértice (el quinto en la gura 3.3) a la imagen, Si = Si5  D0 .
pt
O

3.3.1. La matriz de transferencia

Podemos corroborar el comportamiento predicho por la lente delgada
equivalente calculando su matriz de transferencia. Esto signica construir
paso a paso a la matriz que cubre la trayectoria de la luz desde el objeto hasta
 94 CAPÍTULO 3. OPTICA CON MATRICES.

la imagen pasando por cada elemento del sistema óptico sin abstracciones. La
gura 3.9 ilustra esta situación. Recordemos que la caja negra considerada
tiene diversos objetos ópticos, cinco en este caso. Es importante notar que

AM
en la construcción de esta matriz la distancia So1 es medida desde el vértice
del primer elemento óptico y que la distancia Si5 tiene al vértice posterior
del último objeto como referencia.
Traducido a la lengua matricial lo anterior es que la matriz de transfe-
rencia T nos regresa el vector de salida

UN
✓ ◆ ✓ ◆
yi yo
=T (3.26)
ni ↵ i no ↵ o

que contiene a los parámetros de la sección cuyo lı́mite es la imagen, yi y ↵i ,

en función de los parámetros yo y ↵o que corresponden a la rebanada inicial.
Explı́citamente la matriz T es
✓ Si5 ◆ ✓
1 n a b
◆✓
SS
1  Sno1
◆ ✓
A B
◆
T= = . (3.27)
-D
0 1 c d 0 1 C D

Si existe la lente equivalente al sistema que imaginamos como una caja negra
su matriz de transferencia entre objeto e imagen debe ser la misma. Esto es,

✓ ◆✓ ◆✓ ◆
M

1 Si /n 1 0 1 So /n
T=
0 1 P 1 0 1
✓ ◆✓ ◆
JJ

1 Si /n 1 So /n
=
0 1 P P So /n + 1
✓ (So +Si )n+P Si So
◆ ✓ ◆
1  Si P 2 T11 T12
= n = . (3.28)
P P So /n + 1 T11 T22
ica

Al factorizar (Si So )/n de entrada T12 obtenemos

 
So Si P 1 1
pt

   + = 0, (3.29)
n n So S i
que debe ser igual a cero por consistencia con la ley de Gauss para una lente
O

delgada. Ambas matrices, ABCD y Tij representan al mismo sistema con

solución única. De manera tal que B = T12 = 0 es el requisito indispensable
para que el sistema óptico representado la matriz T sea un formador de ima-
gen y tenga sentido buscar la lente delgada que le sea equivalente. Haciendo
 3.4. EJEMPLOS 95

(a)

AM
H H’
(b)

UN
Figura 3.10: Geometrı́a correspondiente al ejemplo 1 de este capı́tulo. La
SS
echa roja en (a) representa a un objeto colocado a una distancia de 1.2
m del centro de una esfera con 0.2 m de diámetro. La esfera está hecha con
vidrio con n = 1.5. En (b) hacemos un acercamiento a la esfera representando
-D
con la lı́nea punteada azul a los planos principales.

cero al segundo elemento de T encontramos que en general el vector de salida

al atravesar un sistema creador de imagen es de la forma
✓ ◆ ✓ ◆✓ ◆
M

yi A 0 yo
= . (3.30)
ni ↵ i C D no ↵ o
JJ

La ecuación asociada al primer renglón del vector (3.30) nos provee de una
forma muy práctica para determinar la amplicación transversal, puesto que
yi
MT = = A. (3.31)
yo
ica

3.4. Ejemplos
pt

Dentro de este capı́tulo hemos introducido una nueva herramienta para

estudiar sistemas ópticos. Con matrices podemos describir desde los siste-
mas ópticos más sencillos hasta complicarlos arbitrariamente. Para asentar
O

el aprendizaje de tan útil método recomendamos seguir los siguientes ejem-

plos.
Ejemplo 1
Una esfera de vidrio con n = 1.5 y diámetro d = 0.2 m se coloca a de forma
 96 CAPÍTULO 3. OPTICA CON MATRICES.

tal que su centro queda a 1.2 m del centro de un objeto, como muestra la
gura 3.10 (a). Debemos encontrar la posición de su imagen.
Comencemos por notar que el sistema está compuesto por dos supercies

AM
esféricas de radio R = d/2 = 0.1 m con signo opuesto, separadas por una
sección de transición con longitud t = d. La matriz que conjunta a estos
elementos es
✓ ◆✓ ◆✓ ◆
1 0 1 d/n2 1 0
Mesf = (3.32)
 n1 n 2
1 0 1  n2 n1 1

UN
R
✓ ◆✓ ◆✓ R ◆
1 0 1 0.2/1.5 1 0
= (3.33)
 11.5 1 0 1  1.51
1
✓ 0.1 ◆ ✓ 2
◆ ✓ ◆ 0.1
1 0 1 15 1 0
= (3.34)
5 1 0 1 5 1
✓ ◆ ✓ ◆
=
1/3 2/15
20/3 1/3
SS =
a b
c d
, (3.35)

en donde hemos dado por hecho que la esfera está en el aire por lo que
-D
n1 = n3 = 1. Directamente de nuestra relación general de la matriz abcd con
la distancia focal de la lente delgada equivalente y la posición de los planos
principales obtenemos
1 3
f = = = 0.15 m
M

c 20
1d 3 2 1
D= = = m
c 20 3 10
JJ

1a 1
D0 = =  m. (3.36)
c 10
De las ecuaciones (3.36) notamos que las distancias D y D0 tienen signo con-
ica

trario y tienen una longitud igual al radio de la esfera. Es decir, en este caso
lo dos planos principales H y H 0 parten la esfera por la mitad, ver gura 3.10
de la esfera. Por otro lado, conocemos So que podemos sustituir directamente
en la ecuación de Gauss para encontrar Si puesto que también hemos calcu-
pt

lado la distancia focal para la lente delgada equivalente equivalente dentro

de las ecuaciones (3.36):
O

1 1 1
+ = ) Si = 0.1714 m
1.2 m Si 0.15 m
El cálculo está hecho con los parámetros de la lente delgada equivalente; la
notación Si es adecuada puesto la distancia está medida desde el segundo
 3.4. EJEMPLOS 97

plano principal, 0.10 m a la izquierda del segundo vértice. Por lo tanto la

distancia de este vértice a la imagen es Si2 = 0.1714  0.1 = 0.0714 m. Ya
contamos con la información necesaria para calcular a amplicación trans-

AM
versal
Si 0.1714
MT = = = 0.1452.
So 1.2
Encontraremos el mismo resultado a partir de la matriz de transferencia de
la esfera. Recordemos que esta se obtiene expandiendo su matriz equivalente

UN
una distancia So hacia la izquierda y una distancia Si a la derecha a partir
de los vértices correspondientes. Es decir,
✓ Si2 ◆ ✓ ◆✓ ◆
1 n 1/3 2/15 1  Sno1
Tesfera = Msalida Mesfera Mentrada = .
0 1 20/3 1/3 0 1
(3.37)
SS
Los datos del problema, So1 = 1.1 mm y n = 1, son suciente para calcular
Si2 a partir de las entradas de la matriz Tesfera , solo debemos desarrollarla:
-D
✓ ◆✓ ◆✓ ◆
1 Si2 1/3 2/15 1 1.10
Tesfera = (3.38)
0 1 20/3 1/3 0 1
✓ ◆
0.333  7.333Si2 0.5  7Si2
= (3.39)
7.333 7
M

En la sección 3.3.1 aprendimos que en cualquier matriz de transferencia co-

rrespondiendo a un objeto óptico formador de imagen el segundo elemento
JJ

B = 0 Por lo tanto Si2 = 1/14 = 0.0714 m. Sustituyendo este valor en el

primer elemento A que, también aprendimos, directamente nos arroja la am-
plicación concluimos MT = 0.333  7.333Si2 = 0.142858.
ica

Ejemplo 2
Se encuentran dos lentes delgadas inmersas en aire, colocadas con una dis-
tancia t entre ellas como muestra la gura 3.11. La lente de la izquierda tiene
pt

distancia focal f1 mientras que la lente de la derecha tiene distancia focal

f2 2. ¿Cual es la matriz de la lente delgada equivalente al sistema bifocal?
O

En la sección 3.2 aprendimos que la matriz de una lente con distancia

focal f es
✓ ◆
1 0
Mf = ; (3.40)
 f1 1
 98 CAPÍTULO 3. OPTICA CON MATRICES.

AM
UN
Figura 3.11: Geometrı́a correspondiente al ejemplo 2. Dos lentes con distan-
cias focales f1 y f2 están colocadas con una distancia de t entre ellas. Se
busca la matriz de la lente delgada equivalente tomando en cuenta que el
medio rodeando las lentes es aire.

través de la distancia entre los lentes

✓
SS
el espacio entre las dos lentes está considerado en la matriz de transición a
◆
1 nt
-D
MT = (3.41)
0 1
con n = 1. Por lo tanto la matriz que considera al sistema completo es el
producto Mf 2 Mt Mf 1 ,
✓ ◆✓ ◆✓ ◆
M

1 0 1 nt 1 0
Mf1 ,t,f2 =
 f12 1 0 1  f11 1
✓ ◆✓ ◆
JJ

1 0 1  ft1 t
=
 f12 1  f11 1
✓ ◆
1  ft1 t
= . (3.42)
ica

1 t 1 t
 f2 + f2 f1  f1  f2 + 1
Una vez calculada la matriz del sistema completo podemos identicar a cada
entrada con la matriz general abcd de una lente delgada. Es decir
pt

a = 1  t/f1
b=t
O

1 t 1
c= + 
f2 f2 f1 f1
t
d=1 . (3.43)
f2
 3.4. EJEMPLOS 99

Con lo cual tenemos conocimiento total de la trayectoria de cualquier rayo

de luz paraxial antes y después del travesar el sistema óptico. La distancia
focal está dada por
1 1 1 t

AM
= c = +  , (3.44)
f f2 f1 f2 f1
mientras que las distancias entre los vértices y los respectivos planos princi-
pales son

UN
✓ ◆
1d 1 t 1 t f
D= = 11+ = = t
c c f2 c f2 f2
y
✓ ◆
0 1a 1 t 1 t f
D = = 11+ = = t (3.45)
c c f1 SS c f1 f1

usando la ecuación del fabricante de lentes obtenemos la distancia del segundo

-D
plano principal a la imagen Si en términos de la distancia So del primer plano
principal y de la distancia focal
So f
Si = . (3.46)
f + So
M

Metamos números para jar ideas. Designemos f1 = - 30 cm y f2 = 40 cm y

t = 25 cm. Además colocamos un objeto con 4 cm de altura a 10 cm de la
primer lente. ¿Cuál será la posición y el tamaño de la imagen nal? Iniciamos
JJ

calculando el la distancia focal directamente de c,

1 1 1 25 1
=  = ,
f 30 40 (30)(40) 80
ica

obteniendo f = 80 cm. Ahora sustituimos en las ecuaciones (3.45) para

encontrar las distancias a los planos principales,
pt

(25)(80)
D= = 50 cm
40
O

(25)(80) 200
D0 =  = = 66.666 cm
30 3
y ası́, poder calcular las propiedades de la imagen. Nuevamente calculamos
la distancia del objeto al primer plano principal, So = 10  50 = 60 cm.
 100 CAPÍTULO 3. OPTICA CON MATRICES.

Lente 1 Lente 2 H H´

AM
Figura 3.12: Esquemático a escala del ejemplo 2 incluyendo los planos prin-
cipales H en rojo y H 0 en azul.

UN
Con lo que es directo encontrar la distancia del segundo plano principal a
la imagen usando la ley de Gauss, Si = - 240 cm. Por lo tanto la distancia
del segundo vértice a la imagen es Si2 = - 240 + 66.666 = -173.333 cm.
Finalmente obtenemos la amplicación del sistema calculando el cociente
MT = Si /So = 240/60 = 4. La gura 3.12 es una ilustración a escala del
SS
sistema incluyendo a sus planos principales H y H 0 .
-D
3.5. Suplemento. Aperturas y diafragmas.
3.5.1. Diafragma de apertura y diafragma de campo.
Consideremos un par objeto-imagen en un sistema óptico con eje, por
M

ejemplo, una lente convergente como se muestra en la gura 3.13. Solo un

conjunto de los rayos luminosos que salen de un punto en el objeto pasará
por el sistema óptico y servirá para formar su imagen. En este caso el borde
JJ

de la lente limita el paso de rayos luminosos y dene el cono de rayos que

servirán para la formación de la imagen.
Cuando se tiene un sistema con más elementos, habrá uno de ellos que
limite el ingreso de rayos luminosos desde un punto del objeto. En ocasiones
ica

también se agregan diafragmas circulares, centrados en el eje óptico, cuyo

objetivo es justamente limitar el paso de rayos luminosos y ası́ mejorar la
calidad de la imagen. A este tipo de diafragmas, que limitan el cono de
pt

luz que entra a un sistema óptico, se les llamma diafragma de apertura. Por
ejemplo, en la gura 3.14 se muestra el ejemplo de un diafragma que se coloca
entre dos lentes que forman un sistema óptico.
O

Por otra parte, el sistema óptico también limita la extensión de objetos en

el plano del objeto que aparecerán en la imagen. Por ejemplo, en el sistema
formado por una lente simple y un sensor de luz rectangular como el que se
muestra en la gura 3.15. El tamaño nito del sensor hace que algunos objetos
 3.5. SUPLEMENTO. APERTURAS Y DIAFRAGMAS. 101

AM
UN
SS
-D

Figura 3.13: Limitación del paso de rayos luminosos por una lente.
M

en el plano del objeto se salgan de la imagen. Ası́ en nuestro ejemplo, la echa

JJ

roja sı́ aparece en la imagen, pero la echa azul no sale completa. Decimos
entonces que la imagen de la echa azul se sale de campo. En ocasiones
también se agregan diafragmas que denen el campo de un sistema óptico.
Llamaremos diafragma de campo al elemento que determina el área máxima
ica

en el plano de la imagen.
El objetivo de este suplemento es establecer las principales deniciones
de los diafragmas de apertura, ası́ como estudiar su efecto en la formación
pt

de la imagen en sistemas ópticos.

O

3.5.2. Pupilas.
Analicemos el comportamiento del par objeto-imagen en un sistema ópti-
co. Es claro que, para este par, uno de los elementos del sistema es el que
limita la apertura del cono de rayos luminosos que salen de un punto del ob-
 102 CAPÍTULO 3. OPTICA CON MATRICES.

AM
UN
Figura 3.14: Un diafragma colocado entre dos lentes.

SS
jeto. A la apertura circular de este elemento es a lo que llamamos el difragma
de apertura. Ahora nos preguntamos cómo se ve el diafragma de apertura
-D
visto desde el objeto. Para responder a esta pregunta lo que tenemos que
hacer es obtener la imagen que del diafragma de apertura forman todos los
elementos que se encuentran a su izquierda, esto es, la imagen del diafragma
por los elementos del sistema entre el objeto y el diafragma. La pupila de
entrada es la imagen circular del diafragma de apertura vista desde el ob-
M

jeto. Por construcción, el cono de rayos luminosos que salen del objeto en
dirección de la pupila de salida serán los rayos que pasarán por el sistema
y formarán la imagen. De igual manera, si construimos la imagen que del
JJ

diafragma de apertura forman todos los elementos que están a su derecha

(entre el diafragma y el plano de la imagen), habremos obtenido la pupila
de salida del sistema. Esto quiere decir que el cono de rayos que sirvió para
formar un punto de la imagen pareció haber salido de la pupila de salida.
ica

El caso de una lente delgada es muy simple. El diafragma de apertura,

la pupila de entrada y la pupila de salida coinciden con el borde de la lente.
Cuando se tienen más de dos elementos hay que hacer un análisis más cuida-
pt

doso, para determinar el diafragma de apertura y las pupilas. Por ejemplo,

en la gura 3.16 se muestra la posición de la pupila de entrada y la pupila
de salida para el caso del sistema de la gura 3.14 formado por dos lentes
O

y un diafragma entre ellas. . El objeto está colocado en el de la gura 3.14

punto A, de tal manera que el diafragma de apertura (D.A.) corresponde al
diafragma entre las lentes. La pupila de entrada (P.E.) es la imagen del dia-
fragma que produce la lente L1 . Es importante subrayar que para localizar
 3.5. SUPLEMENTO. APERTURAS Y DIAFRAGMAS. 103

AM
UN
SS
-D

Figura 3.15: El área del sensor determina cuáles de los objetos aparecen en
M

la imagen.
JJ

esta imagen estamos propagando la luz desde el diafragma hacia la lente y

luego al objeto (de derecha a izquierda), en sentido contrario a como hemos
acostumbrado a obtener la imagen en sistemas ópticos. La lente L1 produce
una imagen virtual del diafragma, que está indicada como la base del cono
ica

de rayos de entrada al sistema (cono amarillo). La pupila de salida (P.S.) es

la imagen del diafragma producida por la lente L2 . Se trata también de una
imagen virtual, que es la base del cono de rayos de salida que convergen en
pt

el punto B (cono naranja).

O

3.5.3. Rayo guı́a. Rayos marginales.

Antes de continuar con ejemplos de determinación de diafragmas y pupilas
conviene denir un conjunto de rayos que son muy útiles en la construcción
de la imagen. Llamamos rayo guı́a al rayo que sale de un punto del objeto
 104 CAPÍTULO 3. OPTICA CON MATRICES.

AM
UN
Figura 3.16: Posición de pupila de entrada y pupila de salida de la g. 3.14.

SS
-D

Figura 3.17: Rayo guı́a y un rayo marginal para el sistema de la g. 3.14.

en dirección del centro de la pupila de entrada. Por construcción este rayo

M

pasará por el centro del diafragma de apertura y por el centro de la pupila de

salida. De manera semejante, se denen los rayos marginales como aquéllos
JJ

que salen de un punto en el objeto en dirección al borde de la pupila de

entrada. Por tanto, estos rayos marginales pasan por el borde del diafragma
de apertura y salen del sistema pasando por el borde de la pupila de salida.
En la gura 3.17 se muestra el rayo guı́a y un rayo marginal que parten del
ica

extremo del objeto (A) en el sistema de la gura 3.14.

3.5.4. Uso de matrices para identicar el diafragma de

pt

apertura.
Las matrices sirven para identicar el diafragma de apertura en sistemas
O

ópticos más complicados, que tienen varios elementos refractores y varios

diafragmas. Hemos aprendido a calcular la matriz de transferencia hasta uno
de los elementos, digamos el número j. Esa matriz es el producto de todas las
matrices de transferencia anteriores, que incluyen las de los elementos ópticos
 3.5. SUPLEMENTO. APERTURAS Y DIAFRAGMAS. 105

y las de espacio libre. Colocamos el objeto en la posición x0 = 0, quedando

el elemento j en la posición xj . Escribimos la ecuación de transferencia desde
el objeto hasta el elemento j como

AM
✓ ◆ ✓ ◆✓ ◆ ✓ ◆
yj aj b j 0 bj ↵ 0
= = , (3.47)
↵j cj d j ↵0 dj ↵0

Recordemos que yj = bj ↵o es la altura con la que el rayo llega al elemento

j habiendo salido desde un punto en el objeto con y0 = 0 y que forma un

UN
ángulo ↵0 con el eje. Esta expresión nos permite entonces calcular el ángulo
↵0j que tendrı́a que tener un rayo que sale del objeto para tocar el borde del
elemento j, que suponemos tiene una altura rj . Este ángulo está dado por
 
 
1  rj 
↵0j = tan  
bj SS
De esta manera podemos obtener los ángulos de salida ↵0j para tocar los
-D
bordes de cada uno de los elementos del sistema óptico. Por denición, el
diafragma de apertura corresponde al elemento que limita el cono de rayos
luminosos que entran al sistema, y por tanto corresponde al mı́nimo

mı́n {↵01 , ↵02 , ..., ↵0n }

M

Una vez que se ha identicado el diafragma de apertura se pueden emplear

las matrices antes y después de este elemento para calcular la posición y el
JJ

tamaño de las pupilas de entrada y de salida.

3.5.5. Ejemplos de análisis de diafragmas de apertura

ica

y pupilas.
Concluimos este suplemento con varios ejemplos de cómo identicar el
diafragma de apertura de sistemas ópticos. Una vez que se identica el dia-
pt

fragma de apertura se mostrará la manera de obtener las pupilas de entrada

y de salida. Es importante subrayar que el diafragma de apertura de un sis-
tema depende de la posición del objeto. Si se mueve el objeto es posible que
O

cambie el elemento del sistema que dene el diafragma de apertura.

Ejemplo 1. Se coloca una apertura circular de 3.0 cm de diámetro frente

a una lente delgada, a una distancia de 1.5 cm. El diámetro de la lente
 106 CAPÍTULO 3. OPTICA CON MATRICES.

es 4.80 cm y su distancia focal es 3.5 cm. Se coloca un objeto a 8 cm

de la lente. Localizar la posición y el tamaño de las pupilas de entrada
y de salida.

AM
El sistema es lo sucientemente sencillo como para poder estudiarlo
con la ecuación de Gauss para lente delgada. Calculamos los ángulos
que tendrı́a que tener un rayo que parte del objeto para llegar al borde
del diafragma o al borde de la lente. E diafragma está a 8.0  1.5 = 6.5
cm del objeto. El primer ángulo es entonces

UN
1.5
↵01 = tan1 = 12.99
6.5
El segundo ángulo es
2.4
SS
↵02 = tan1
8
= 16.7

Como el menor de estos ángulos es el primero, el diafragma de apertura

-D
del sistema es el diafragma frente a la lente. Como el diafragma de
apertura queda justo frente al objeto, también será la pupila de entrada.
Entonces, la pupila de entrada está 1.5 cm delante de la lente y tiene
un diámetro de 3.0 cm. Para obtener la pupila de salida calculamos la
imagen que la lente produce del diafragma. Empleando la ecuación de
M

Gauss obtenemos la posición de la imagen

1 1 1
 + =
JJ

1.5 si 3.5
obteniendo si = 2.625 cm. Se trata de una imagen virtual (queda
antes de la lente). Con este resultado podemos calcular el diámetro de
ica

la imagen del diafragma.

si 2.625
Di = MT D = = ⇥ 3.0 = 5.25 cm
so 1.5
pt

Entonces, la pupila de salida se encuentra a 2.625 cm delante de la lente

y tiene un diámetro de 5.25 cm. En la gura 3.18 se muestra la gura
O

que corresponde a este ejemplo.

Ejemplo 2. En el ejemplo anterior, determinar las posiciones del ob-
jeto para las que la apertura circular funciona como el diafragma de
apertura del sistema.
 3.5. SUPLEMENTO. APERTURAS Y DIAFRAGMAS. 107

AM
UN
SS
-D
M

Figura 3.18: Diafragma y una lente del ejemplo 1.

Aquı́ conviene preguntarnos si hay una posición del objeto (a la izquier-

JJ

da de la lente) para la que los dos ángulos son iguales. Si x es la posición

de la lente respecto al objeto, el primer ángulo, hasta el diafragma, es
1.5
↵01 = tan1
ica

x  1.5
mientras que el segundo ángulo, hasta la lente, es
2.4
↵02 = tan1
pt

x
Los dos ángulos serán iguales cuando
O

1.5 2.4
=
x  1.5 x
Resolviendo se obtiene x = 1.11 cm. Se puede demostrar que para
mayores valores de x el primer ángulo es siempre menor, y para x < 1.11
 108 CAPÍTULO 3. OPTICA CON MATRICES.

cm el segundo ángulo es siempre menor. En el primer caso el diafragma

es el diafragma de apertura y en el segundo lo es la lente.

Ejemplo 3. Se tiene un sistema de dos lentes delgadas (f1 = 8 cm,D1 =

AM
11 cm, f2 = 10 cm, D2 = 6 cm) separadas una distancia de 6 cm.
Se coloca un diafragma cuyo diámetro es 8 cm entre las lentes, a 2 cm
de la primera lente. Se coloca un objeto a 22 cm de la primera lente.
Demostrar que para esta posición del objeto el diafragma actúa como

UN
diafragma de apertura. Calcular la posición y el tamaño de la pupila
de entrada y de la pupila de salida del sistema. Estos son los datos con
los que se generó la gura 3.14.
Este sistema sirve para ilustrar las ventajas del empleo de matrices.
Con el objeto situado en x = 0 escribimos las matrices de propagación.

M1 =
SS
La primera es propagación en medio homogéneo (n = 1) por 22 cm:
✓
1 22
◆

0 1
-D
Le sigue la matriz de lente delgada con una distancia focal f1 = 8 cm:
✓ ◆
1 0
M2 =
 18 1
M

Después tenemos la matriz de propagación hasta el diafragma

✓ ◆
JJ

1 2
M3 =
0 1

seguida de la matriz de propagación hasta la segunda lente

ica

✓ ◆
1 4
M4 =
0 1
pt

Aquı́ conviene notar que el producto de estas dos últimas matrices

es igual a la matriz de propagación de la primera a la segunda lente.
O

Aunque no la necesitamos para el cálculo del diafragma de apertura,

la última matriz del sistema es la de la segunda lente (f2 = 10 cm.
✓ ◆
1 0
M5 = 1
10
1
 3.5. SUPLEMENTO. APERTURAS Y DIAFRAGMAS. 109

Calculamos el primer ángulo, que es el que se forma con el borde de la

primera lente. Está dado por
5.5 1

AM
↵1 = tan1 = tan1
M1 (1, 2) 4

Para el segundo ángulo necesitamos propagar hasta el diafragma con

el producto ✓ 3 ◆
37

UN
M3 M2 M1 = 4 2
 18  74
El segundo ángulo está dado por el cociente del radio del diafragma
entre el elemento (2, 1) de este producto:

8
↵2 = tan1 SS
37
Por último necesitamos el ángulo que se forma con el borde de la se-
-D
gunda lente. Necesitamos propagar hasta la tercera lente con la matriz
✓ 1 23
◆
M 4 M 3 M2 M1 = 4 2 (3.48)
1 7
8 4
M

resultando un ángulo
3⇥2
↵3 = tan1
23
JJ

El diafragma de apertura está determinado por el elemento que da el

menor de estos ángulos, que claramente es ↵2 . Por tanto hemos demos-
trado que el diafragma entre las lentes es el diafragma de apertura.
ica

Ahora calculamos la posición y el tamaño de la pupila de entrada.

Propagamos desde el diafragma hacia la izquierda para localizar la
imagen que produce la primera lente:
pt

✓ ◆✓ ◆✓ ◆ ✓ ◆
1 si 1 0 1 2 1  s8i 2 + 3s4 i
=
0 1  18 1 0 1  18 3
4
O

Para la posición de la pupila de salida obtenemos el valor de si que

hace al elemento (1, 2) de este producto igual a cero, lo que da

xP E = si = 2.333
 110 CAPÍTULO 3. OPTICA CON MATRICES.

respecto de la lente. Esto nos dice que está la pupila de entrada está
del mismo lado que el diafragma. El tamaño (diámetro) de la pupila
de salida se obtiene al sustituir este valor de si en el elemento (1, 1)

AM
del producto de matrices, lo que da la apmplicación lateral, que luego
hay que multiplicar por el diámetro del diafragma:
✓ ◆
2.333
1+ ⇥ 8 = 10.667
8

UN
De manera similar, si se calcula la imagen del diafragma producida por
la segunda lente encontramos que la pupila de salida se encuentra a
si = 2, 857 cm de la segunda lente (del mismo lado que el diafragma)
y que su diámetro es 5.714 cm.
SS
-D
M
JJ
ica
pt
O