Esta página, es una aplicación de las integrales de Fresnel y de la espiral de Cornu.

Las funciones fresnelc y fresnels de MATLAB, calculan las integrales de Fresnel C(x) y S(x), respectivamente

En los tres ejemplos, se indica el procedimiento geométrico para el cálculo de la intensidad debida a la difracción producida por

• Una abertura rectangular
• El borde de una placa semiinfinita
• Un obstáculo en forma de placa rectangular

Difracción producida por una rendija

Se observa la evolución desde la difracción Fresnel a la difracción Fraunhofer a medida que se disminuye la anchura de la rendija.

dx=2; %ancho rendija
x=0:0.02:5;
s1=fresnels(x+dx/2);
s2=fresnels(x-dx/2);
c1=fresnelc(x+dx/2);
c2=fresnelc(x-dx/2);
y=((s2-s1).^2+(c2-c1).^2)/2;
plot(x,y)
grid on
xlabel('x')
ylabel('Intensidad')
title('difracción rendija')

Actividades

• Se introduce la anchura de la rendija en el control titulado Ancho

Se pulsa el botón titulado Nuevo

Se mide la intensidad de la luz en una posición dada, que se puede cambiar en el control titulado Posición

Difracción producida por un borde de una placa

En el programa interactivo observamos que la intensidad no es nula en la zona de sombra de la placa y en las proximidades del borde.

Definimos una intensidad umbral como aquella I=0.01·I₀. Moviendo el pequeño cuadrado de color rojo hacia la izquierda (véase la figura más abajo) vemos que corresponde a u=2.2. (en la parte superior izquierda, se proporcionan los datos de u y de la intensidad  I/I₀.

Por ejemplo, en el caso del sonido si I₀=10^-4 W/m² es decir, 80 decibeles, I=0.01·I₀=10^-6 W/m², es 60 dB.

El parámetro u vale

${\displaystyle u=x\sqrt{\frac{2}{\lambda y}}}$

donde x e y son las coordenadas del punto de observación P.

En la figura, se muestra el lugar geométrico de los puntos tales que u=2.2 o los puntos (x, y) para los cuales la intensidad es el 1% de la intensidad del movimiento ondulatorio que se mediría sin la presencia del obstáculo.

• La curva en color rojo, se ha trazado para un movimiento ondulatorio de longitud de onda λ=0.5 m
• La curva en color azul, para λ=0.1 m.

En el caso del sonido, la velocidad de propagación es 340 m/s,

• Un sonido de frecuencia de 680 Hz le corresponde una longitud de onda λ=0.5 m
• Un sonido de 3400 Hz le corresponde una una longitud de onda λ=0.1 m

Como apreciamos en la figura, oiremos sonido detrás de un edificio, alrededor de la esquina.

En el caso de la luz visible, su longitud de onda es del orden de λ=5·10^-7 m. La curva que se dibujaría sería indistinguible de la recta horizontal.

El sonido por tanto se difracta dentro de la región de sombra de la placa semiinfinita. En una situación real, el sonido escuchado detrás de un obstáculo puede provenir no solamente del efecto de la difracción, sino también de la  reflexión por otros obstáculos que puedan estar presentes.

x=-3:0.02:5;
s=fresnels(-x);
c=fresnelc(-x);
y=((0.5-s).^2+(0.5-c).^2)/2;
plot(x,y)
grid on
xlabel('x')
ylabel('Intensidad')
title('difracción borde')

Actividades

Se mide la intensidad de la luz en una posición dada, que se establece en el control titulado Posición

Difracción producida por un obstáculo en forma de placa rectangular

Observamos que siempre hay una región iluminada en el eje central del obstáculo

dx=2; %ancho obstáculo
x=0:0.02:5;
s1=fresnels(x+dx);
s2=fresnels(x-dx);
c1=fresnelc(x+dx);
c2=fresnelc(x-dx);
y=((1-s1+s2).^2+(1-c1+c2).^2)/4;
plot(x,y)
grid on
xlabel('x')
ylabel('Intensidad')
title('difracción obstáculo')

Actividades

• Se introduce la anchura de la rendija en el control titulado Ancho

Se pulsa el botón titulado Nuevo

Se mide la intensidad de la luz en una posición dada, que se puede cambiar en el control titulado Posición

Referencias

Hecht E., Zajac A. Óptica. Editorial Addison-Wesley Iberoamericana (1986), págs. 388-418

Ferguson J. L., Why can we hear but not see around a corner?. Am. J. Phys. 54 (7) July 1986, pp. 661-662

Press W. H., Teukolsky S. A., Vetterling W. T., Flannery B. P. Numerical Recipes in C, Second edition, Special functions. Fresnel Integrals, Cosine and Sine Integrals Chapter 6º. Cambridge University Press. Código en C adaptado por el autor al lenguaje Javascript