Imagen de un pez dentro de una pecera

Representamos la pecera por una esfera de plástico transparente o de vidrio de radio R, llena de agua. Supondremos que el recipiente está hecho de un material de pequeño espesor que no afecta a la marcha de los rayso de luz. El índice de refracción de la esfera es n1 que es el del agua y el índice de refracción del aire es n2

Situamos el origen en el centro C de la pecera y los ejes X (horizontal) y Z (vertical)

Vamos a determinar la posición x1 y tamaño z1 de la imagen (flecha de color rojo) del pez representado por una flecha vertical azul de altura z0, situada en la posición x0

Trazamos un rayo de luz que pasa por el origen C, que no se desvía al pasar del agua al aire.

Trazamos otro paralelo al eje X, que se refracta en el punto P. El ángulo del rayo incidente es θ1 y el refractado es θ2. La ley de Snell de la refracción es

n 1 sin θ 1 = n 2 sin θ 2

El ángulo incidente θ1 tiene un valor límite que es aquél para el cual θ2=π/2, reflexión total

sin θ c = n 2 n 1

Dado el valor del radio R de la pecera, la altura del objeto z0<Rsinθc

La posición y tamaño de la imagen virtual se obtiene de la intersección del rayo que pasa por el centro de la esfera y de la prolongación del rayo refractado en el punto P

La ecuación de la recta z=mx+b que pasa por el centro C de la esfera y por el punto (x0, z0) es

b=0 z=mx z 0 =m x 0 ,m= z 0 x 0 z= z 0 x 0 x

Determinamos la ecuación de la recta que pasa por el punto P y forman un ángulo π-α con el eje X. La pendiente m es

m=tan( πα )=tan( θ 2 θ 1 )

La ordenada en el origen b se determina, sabiendo que pasa por el punto P ( R 2 z 0 2 , z 0 )

z 0 =tan( θ 2 θ 1 ) R 2 z 0 2 +b

La ecuación de la recta es

z=tan( θ 2 θ 1 )x+ z 0 +tan( θ 2 θ 1 ) R 2 z 0 2 z= z 0 tan( θ 2 θ 1 )( x R 2 z 0 2 )

La intersección de ambas rectas nos da la posición x1 y la altura z1 de la imagen

{ z= z 0 tan( θ 2 θ 1 )( x R 2 z 0 2 ) z= z 0 x 0 x x 1 = z 0 +tan( θ 2 θ 1 ) R 2 z 0 2 z 0 x 0 +tan( θ 2 θ 1 ) z 1 = z 0 x 0 x 1

Representamos la recta que pasa por el centro y la prolongación del rayo refractado. El objeto (en azul) y la imagen (en rojo)

n1=4/3; %agua
n2=1; %aire
z0=2/15; %tamaño del objeto
R=1; %radio de la esfera
x0=-0.5; %posición del obejto
hold on
th=linspace(0,360)*pi/180;
fill(cos(th),sin(th),'c')

fplot(@(t) cos(t), @(t) sin(t),[0,2*pi])
m=z0/x0;
line([-2,1],[-2*m,m])
th_1=asin(z0/R);
th_2=asin(n1*sin(th_1)/n2);
m=tan(th_2-th_1);
line([-2,2],[z0-m*(-2-sqrt(R^2-z0^2)),z0-m*(2-sqrt(R^2-z0^2))])
x1=(z0+m*sqrt(R^2-z0^2))/(z0/x0+m);
z1=z0*x1/x0;
line([x0,x0],[0,z0],'lineWidth',1.5,'color','b') %objeto
line([x1,x1],[0,z1],'lineWidth',1.5,'color','r') %imagen
hold off
grid on
axis equal
xlabel('x_0/R')
ylabel('M')
title('Objeto e imagen')

Posición de la imagen

Expresamos la posición x1 y tamaño z1 de la imagen en términos de la posición x0/R y tamaño z0/R del objeto y del cociente de los índices de refracción n2/n1 (aire y agua)

sin θ 1 = z 0 R ,tan θ 1 = sin θ 1 1 sin 2 θ 1 = z 0 R 1 ( z 0 R ) 2 = z 0 R 2 z 0 2 tan θ 2 = sin θ 2 1 sin 2 θ 2 = n 1 n 2 sin θ 1 1 n 1 2 n 2 2 sin 2 θ 1 = n 1 n 2 z 0 R 1 n 1 2 n 2 2 z 0 2 R 2 = n 1 z 0 n 2 2 R 2 n 1 2 z 0 2 tan( θ 2 θ 1 )= tan θ 2 tan θ 1 1tan θ 2 tan θ 1 = n 1 z 0 n 2 2 R 2 n 1 2 z 0 2 z 0 R 2 z 0 2 1 n 1 z 0 n 2 2 R 2 n 1 2 z 0 2 z 0 R 2 z 0 2 = z 0 n 1 R 2 z 0 2 n 2 2 R 2 n 1 2 z 0 2 n 2 2 R 2 n 1 2 z 0 2 R 2 z 0 2 n 1 z 0 2 x 1 = x 0 z 0 +tan( θ 2 θ 1 ) R 2 z 0 2 z 0 + x 0 tan( θ 2 θ 1 ) = x 0 z 0 + z 0 n 1 R 2 z 0 2 n 2 2 R 2 n 1 2 z 0 2 n 2 2 R 2 n 1 2 z 0 2 R 2 z 0 2 n 1 z 0 2 R 2 z 0 2 z 0 + x 0 z 0 n 1 R 2 z 0 2 n 2 2 R 2 n 1 2 z 0 2 n 2 2 R 2 n 1 2 z 0 2 R 2 z 0 2 n 1 z 0 2 x 1 = x 0 n 2 2 R 2 n 1 2 z 0 2 R 2 z 0 2 n 1 z 0 2 +( n 1 R 2 z 0 2 n 2 2 R 2 n 1 2 z 0 2 ) R 2 z 0 2 n 2 2 R 2 n 1 2 z 0 2 R 2 z 0 2 n 1 z 0 2 + x 0 ( n 1 R 2 z 0 2 n 2 2 R 2 n 1 2 z 0 2 ) = x 1 R = x 0 R ( n 2 n 1 ) 2 ( z 0 R ) 2 1 ( z 0 R ) 2 ( z 0 R ) 2 +( 1 ( z 0 R ) 2 ( n 2 n 1 ) 2 ( z 0 R ) 2 ) 1 ( z 0 R ) 2 ( n 2 n 1 ) 2 ( z 0 R ) 2 1 ( z 0 R ) 2 ( z 0 R ) 2 + x 0 R ( 1 ( z 0 R ) 2 ( n 2 n 1 ) 2 ( z 0 R ) 2 )

Representamos la posición x1/R de la imagen de un objeto de altura z0/R=2/15, en función de la posición del objeto en el intervalo 1 z 0 2 < x 0 < 1 z 0 2 . El cociente de los índices de refracción, n2/n1=3/4

n=3/4; %n2/n1=aire/agua
z0=2/15; %altura objeto
x1=@(x) x*(sqrt(n^2-z0^2)*sqrt(1-z0^2)-z0^2+(sqrt(1-z0^2)-sqrt(n^2-z0^2))*
sqrt(1-z0^2))./(sqrt(n^2-z0^2)*sqrt(1-z0^2)-z0^2+x*(sqrt(1-z0^2)-sqrt(n^2-z0^2)));
fplot(x1,[-sqrt(1-z0^2),sqrt(1-z0^2)])
line([-sqrt(1-z0^2),sqrt(1-z0^2)],[-sqrt(1-z0^2),sqrt(1-z0^2)],'color','k')
grid on
xlabel('x_0/R')
ylabel('x_1/R')
title('Posición x/R')

Comparando la representación de x1/R en color azul, con la recta x1/R=x0/R, en color negro, apreciamos

Aumento

La altura z1 de la imagen es z 1 = z 0 x 0 x 1

El aumento M

M= z 1 z 0 = x 1 x 0

Representamos el aumento M en función de la posición x0 del objeto

n=3/4; %n2/n1
z0=2/15;
M=@(x) (sqrt(n^2-z0^2)*sqrt(1-z0^2)-z0^2+(sqrt(1-z0^2)-sqrt(n^2-z0^2))*
sqrt(1-z0^2))./(sqrt(n^2-z0^2)*sqrt(1-z0^2)-z0^2+x*(sqrt(1-z0^2)-sqrt(n^2-z0^2)));
fplot(M,[-sqrt(1-z0^2),sqrt(1-z0^2)])
grid on
xlabel('x_0/R')
ylabel('M')
title('Aumento')

El tamaño de la imagen se incrementa con x0, a medida que el pez se aleja del observador A

Referencias

R De Luca. Watching a goldfish swimming in a fishbowl. Phys. Educ. 60 (2025) 055005