Propagación de ondas electromagnéticas. Ecuaciones de Fresnel.
Campo eléctrico
Cuando un dieléctrico no polar se coloca en un campo eléctrico las cargas positivas y negativas de las moléculas se desplazan unas respecto de las otras, las moléculas se polarizan. El momento dipolar inducido por unidad de volumen , se llama polarización eléctrica
donde σb, C/m2, es la densidad superficial y ρb, C/m3, la densidad (en volumen) de carga ligada
La ley de Gauss se puede expresar
donde la densidad de carga total ρt es la suma de la densidad de carga libre ρf y la densidad de carga ligada ρb
Se define el vector desplazamiento eléctrico
donde εr es la permitividad relativa o constante diléctrica y ε es la permitividad del medio. El vector está relacionado con la densidad de carga libre
Superficie de separación entre dos medios
En la separación de dos medios diferentes, el potencial V debe ser continuo, ya que una discontinuidad implicaría un campo eléctrico infinito. El potencial debe anularse en el infinito si las distribuciones de carga tiene extensión finita y debe ser constante en un conductor si las cargas eléctricas están en reposo.
Componente normal del desplazamiento eléctrico
Componente tangencial de la intensidad del campo eléctrico

Consideremos una superficie cilíndrica de base S y de altura muy pequeña entre los dos medios.
El flujo de solamente se produce a través de las bases, ya que el área lateral es arbitrariamente pequeña.
Si en la superficie de separación existe una densidad superficial de carga libre σf, el flujo del vector a través de la superficie cilíndrica es
Dn1 y Dn2 son las componentes normales de
Si en la separación de dos medios dieléctricos, la densidad superficie la de carga libre es σf es nula, en estas condiciones la componente normal de es continua al cruzar la superficie de separación.
Si la superficie separa un conductor de un dieléctrico, el campo eléctrico es constante en el conductor, y Dn=σf en el dieléctrico, siendo σf la densidad de carga libre en la superficie del conductor.

El campo eléctrico es conservativo. Tomamos un camino cerrado rectangular ABCD, los lados AB y CD tienen longitud l, los lados AD y BC tiene una longitud arbitrariamente pequeña. La circulación
La componente tangencial del campo eléctrico es continua al cruzar la superficie de separación entre dos medios.
Si la separación es la de un conductor y un dieléctrico, entonces en el conductor y Et=0 en el dieléctrico. El campo eléctrico el perpendicular a la superficie del conductor.
De las dos relaciones se deduce que
Dividiendo la segunda entre la primera
El ángulo mayor con la normal ocurre en el medio de permitividad relativa mayor
Campo magnético
La imanación de los materiales magnéticos se corresponde con la polarización de los materiales dieléctricos.
Según el modelo de Ampère, la imanación produce una densidad de corriente equivalente en la superficie y una densidad de corriente equivalente en el volumen donde
La inducción magnética se puede calcular tanto dentro como fuera del material magnético tratando las corrientes equivalentes como si fueran corrientes reales de conducción circulando en el vacío.
La intensidad del campo magnético está relacionado con y del siguiente modo
donde μr es la permeabilidad relativa y μ la permeabilidad
La ley de Ampére se escribe
donde if es la intensidad de la corriente de cargas libres que atraviesa el camino cerrado
La ley de Gauss para el campo magnético es
Superficie de separación entre dos medios
Componente normal del campo magnético
Componente tangencial de

Consideremos una superficie cilíndrica de base S y de altura muy pequeña entre los dos medios.
El flujo de solamente se produce a través de las bases, ya que el área lateral es arbitrariamente pequeña.
La componente normal de es continua al atravesar la superficie de separación

Tomamos un camino cerrado rectangular ABCD, los lados AB y CD tienen longitud l, los lados AD y BC tiene una longitud arbitrariamente pequeña. La circulación
Si no hay corrientes que atraviesen el camino cerrado ABCD
La componente tangencial de es continua al cruzar la superficie de separación entre dos medios.
De las dos relaciones se deduce
Dividiendo la segunda entre la primera
El ángulo mayor con la normal ocurre en el medio de permeabilidad relativa mayor
Ecuaciones de Maxwell. Ondas electromagnéticas
Las ecuaciones de Maxwell en forma diferencial e integral son
Diferencial | Integral |
---|---|
Otra forma alternativa a la cuarta ecuación de Maxwell que se utilizará en la dedución de la ecuación de las ondas electromagnéticas es
Para obtener la segunda, se multiplica la primera por μ y se utiliza las relaciones entre los vectores,
Densidad de carga y corrientes
Nomenclatura
- ρf, es la densidad de carga libre
- , es la densidad de carga ligada
- , es la densidad de corriente de cargas libres
- , es la densidad de corriente de polarización
- , es la densidad de corriente equivalente en la materia imanada
- ε, es la permitividad del material, ε0=8.8542·10-12 F/m es la del vacío
- μ, es la permeabilidad, μ0=4π·10-7 H/m es la del vacío
Ondas electromagnéticas
En este apartado utilizamos la propiedad
Campo eléctrico
Campo magnético
Partimos de la ley de Faraday
Aplicamos a ambos miembros el operador
Obtenemos
En las regiones en las que se anulan ρf y , no hay densidad de carga libre ni corrientes de cargas libres seobtiene la ecuación de ondas homogénea
Partimos de la ley de Ampère-Maxwell
Aplicamos a ambos miembros el operador
Obtenemos
En la región en la que no hay corrientes , se cumple la ecuación de ondas homogénea
Estudiaremos la propagación de ondas electromagnéticas en medios homogéneos, isótropos, lineales y estacionarios. No conductores y no magnéticos
Un medio es homogéneo si sus propiedades no varían con la posición
Un medio es isótropo, si para un punto dado, sus propiedades son las mismas en todas las direcciones
Es lineal si la permitividad ε y la permeabilidad μ son constantes e independientes de y
Un sistema se dice estacionario, cuando las variables que definen su comportamiento no cambian con el tiempo
El caso más sencillo es el de una onda plana que se propaga a lo largo del eje Z. El campo eléctrico es independiente de x e y
Se dice que es una onda electromagnética linealmente polarizada por que el campo eléctrico oscila en el eje X y el campo magnético oscila en el eje Y
Introduciendo en la ecuacuión de onda para el campo eléctrico

v=ω/k se denomina velocidad de fase
Se define índice de refracción n=c/v el cociente entre la velocidad de la luz en el vacío y la velocidad de la luz en un medio
Las ondas electromagnéticas constan de una campo eléctrico y un campo magnético perpendiculares entre sí y a la dirección de propagación , si el campo eléctrico tiene la dirección del eje X, el campo magnético tiene la dirección del eje Y,
como vamos a comprobar utilizando la ley de Faraday
Esta es la relación entre las amplitudes de la onda eléctrica E0 y la onda magnética H0
Representamos la onda eléctrica y la onda magnética con las mismas amplitudes. Las flechas indican el vector campo eléctrico y el vector campo magnético
z=linspace(0,3*pi,200); xx=zeros(1,length(z)); hold on plot3(cos(z), xx,z,'r') plot3(xx, cos(z),z,'b') z=linspace(0,3*pi,20); xx=zeros(1,length(z)); quiver3(xx,xx,z,xx,cos(z),xx,0, 'color','b'); quiver3(xx,xx,z,cos(z),xx,xx,0, 'color','r'); hold off grid on xlabel('x') ylabel('y') zlabel('z') title('Ondas electromagnéticas')
Vector de Poynting
Se denomina vector de Poynting al producto vectorial
Cuando se integra sobre una superficie cerrada nos da el flujo total de energía hacia el exterior por unidad de tiempo.
Teniendo en cuenta que el valor medio <cos2(ωt)>=1/2, el valor medio de ,
Reflexión y refracción de las ondas electromagnéticas
Supongamos dos medios no magnéticos semiinfinitos, lineales, homogéneos e isótropos. Una onda electromagnética que se propaga en el medio 1 incide sobre la superficie de separación de los medios 1 y 2 produciendo una onda reflejada y otra transmitida. Los vectores unitarios , y son normales a los frentes de onda planos (formados por los vectores y ) y tienen el sentido de la propagación
Los ángulos θi, θr, y θt son los ángulos de incidencia, de reflexión y refracción, respectivamente.
Las leyes de la reflexión y la de Snell de la refracción

La onda electromagnética incidente está formada por frentes de onda planos perpendiculares al vector , que señala la dirección de propagación.
Véase en el apartado Vector normal a un plano la ecuación de un plano perpendicular al vector unitario
El campo eléctrico es de la forma
v1 es la velocidad de fase de la onda electromagnética en el medio 1
Tanto la onda relejada como la refractada por la superficie de separación son también planas y polarizadas
v2 es la velocidad de fase de la onda electromagnética en el medio 2
Los tres vectores han de ser funciones del tiempo similares, por tanto, las frecuencias han de ser iguales
ωi=ωr=ωt
En cualquier punto de la superficie de separación se deberá cumplir

La primera ecuación se escribe
La segunda ecuación se escribe
Como el vector está en el plano de separación, el vector es perpendicular a dicho plano, por tanto los ángulos θi=θr, el ángulo de reflexión es igual al de incidencia. El plano formado por los dos vectores se denomina de incidencia.

El vector entre paréntesis debe ser normal a la superficie de separación. Se tiene que cumplir que las componentes (véase la figura de la derecha)
Ecuaciones de Fresnel
Hay dos posibilidades: que el campo eléctrico sea normal al plano de inicidencia y que sea paralelo al mismo
El campo eléctrico es normal al plano de incidencia
En la superficie de separación, las componentes tangenciales de Y son continuas

Teniendo en cuenta la relación entre amplitudes, y la definición de índice de refracción, , la segunda ecuación se transforma en
Para un medio no magnético, la permeabilidad relativa μr≈1. Resolvemos el sistema de dos ecuaciones con dos incógnitas
Despejamos E0r y E0t. El resultado es
Observamos en las fórmulas que , la onda transmitida y la incidente están en fase. Sin embargo, puede ser positivo o negativo.
En fase, n1>n2
De la ley de la refracción n1sinθi=n2sinθt, si n1>n2, entonces θt>θi, cosθi>cosθt, luego , la onda reflejada está en fase con la incidente.
Para el ángulo de incidencia θc tal que el ángulo de refracción es θt=π/2 (90°), se produce la reflexión total, sinθc=n2/n1, para θi≥θc no hay onda refractada
Representamos y para n1=1.5 y n2=1
n1=1.5; n2=1; th_c=asin(n2/n1); %reflexión total hold on r=@(x) ((n1/n2)*cos(x)-cos(asin(n1*sin(x)/n2)))./((n1/n2)*cos(x) +cos(asin(n1*sin(x)/n2))); fplot(r,[0,th_c]) t=@(x) 2*(n1/n2)*cos(x)./((n1/n2)*cos(x)+cos(asin(n1*sin(x)/n2))); fplot(t,[0,th_c]) hold off grid on set(gca,'XTick',0:pi/12:pi/2) set(gca,'XTickLabel',{'0','\pi/12','\pi/6','\pi/4','\pi/3','5\pi/12','\pi/2'}) xlabel('\theta_i') ylabel('E/E_0') legend('E_{0r}/E_{0i}', 'E_{0t}/E_{0i}','Location','best') title('Campo eléctrico, normal')
En oposición de fase, n1<n2
De la ley de la refracción n1sinθi=n2sinθt, si n1<n2, entonces θt<θi, cosθi<cosθt, luego , la onda reflejada está en oposición de fase (desfasada π) con la incidente.
Representamos y para n1=1 y n2=1.5
n1=1; n2=1.5; hold on r=@(x) ((n1/n2)*cos(x)-cos(asin(n1*sin(x)/n2)))./((n1/n2)*cos(x)+ cos(asin(n1*sin(x)/n2))); fplot(r,[0,pi/2]) t=@(x) 2*(n1/n2)*cos(x)./((n1/n2)*cos(x)+cos(asin(n1*sin(x)/n2))); fplot(t,[0,pi/2]) hold off grid on set(gca,'XTick',0:pi/12:pi/2) set(gca,'XTickLabel',{'0','\pi/12','\pi/6','\pi/4','\pi/3','5\pi/12','\pi/2'}) xlabel('\theta_i') ylabel('E/E_0') legend('E_{0r}/E_{0i}', 'E_{0t}/E_{0i}','Location','best') title('Campo eléctrico, normal')
El campo eléctrico es paralelo al plano de incidencia

En la superficie de separación, las componentes tangenciales de Y son continuas
Para un medio no magnético, μr≈1. Resolvemos el sistema de dos ecuaciones con dos incógnitas
Despejamos E0r y E0t. El resultado es
Observamos en las fórmulas que , la onda transmitida y la incidente están en fase. Sin embargo, puede ser positivo o negativo.
En fase si,
En oposición de fase si,
La ley de la refracción, n1/n2=sinθt/sinθi
La desigualdad se transforma en
Teniendo en cuenta la relación trigonométrica sin(A+B)-sin(A-B)=2sinB·cosA
Esta desigualdad se satisface si
Si no se cumplen las dos condiciones anteriores
La fase de la onda reflejada depende de n2/n1 y también de θi y θt. La razón , puede ser positiva o negativa tanto para n2>n1 como para n2<n1
El ángulo de Brewster
El numerador de se hace cero cuando θi+θt=π/2.
Cuando el campo eléctrico es paralelo no hay onda reflejada. Para este caso, el ángulo de incidencia θi se denomina ángulo de Brewster, θB o ángulo polarizante, ya que una onda electromagnética no polarizada se refleja con su vector normal al plano de incidencia
De la ley de la refracción
Comparamos el ángulo crítico a partir del cual se produce la reflexión total, sinθc=n2/n1, con el ángulo de Brewster, tanθB=n2/n1. Representamos θc y θB en función de la razón n1/n2. La onda electromagnética incide en el medio 1
hold on fplot(@(x) atan(1./x),[1,10]); fplot(@(x) asin(1./x),[1,10]); hold off grid on set(gca,'YTick',0:pi/12:pi/2) set(gca,'YTickLabel',{'0','\pi/12','\pi/6','\pi/4','\pi/3','5\pi/12','\pi/2'}) xlabel('n_1/n_2') ylabel('\theta_B, \theta_c') legend('Brewster', 'Reflexión total','Location','best') title('Angulos crítico y Brewster')
El ángulo crítico θm, es mayor que el ángulo de Brewster, θB
Representamos y para n1=1 y n2=1.5. Se señala mediante una línea vertical a trazos el ángulo de Brewster
n1=1; n2=1.5; hold on r=@(x) ((n1/n2)*cos(asin(n1*sin(x)/n2))-cos(x))./((n1/n2)*cos(asin(n1*sin(x)/n2)) +cos(x)); fplot(r,[0,pi/2]) t=@(x) 2*(n1/n2)*cos(x)./((n1/n2)*cos(asin(n1*sin(x)/n2))+cos(x)); fplot(t,[0,pi/2]) ang_B=atan(n2/n1); %ángulo de Brewster line([ang_B, ang_B],[0,1],'lineStyle','--') hold off grid on set(gca,'XTick',0:pi/12:pi/2) set(gca,'XTickLabel',{'0','\pi/12','\pi/6','\pi/4','\pi/3','5\pi/12','\pi/2'}) xlabel('\theta_i') ylabel('E/E_0') legend('E_{0r}/E_{0i}', 'E_{0t}/E_{0i}','Location','best') title('Campo eléctrico, paralelo')
Representamos y para n1=1.5 y n2=1. Se señala mediante una línea vertical a trazos el ángulo de Brewster
n1=1.5; n2=1; th_c=asin(n2/n1); %reflexión total hold on r=@(x) ((n1/n2)*cos(asin(n1*sin(x)/n2))-cos(x))./((n1/n2)* cos(asin(n1*sin(x)/n2))+cos(x)); fplot(r,[0,th_c]) t=@(x) 2*(n1/n2)*cos(x)./((n1/n2)*cos(asin(n1*sin(x)/n2))+cos(x)); fplot(t,[0,th_c]) ang_B=atan(n2/n1); %ángulo de Brewster line([ang_B, ang_B],[0,2.5],'lineStyle','--') hold off grid on set(gca,'XTick',0:pi/12:pi/2) set(gca,'XTickLabel',{'0','\pi/12','\pi/6','\pi/4','\pi/3','5\pi/12','\pi/2'}) xlabel('\theta_i') ylabel('E/E_0') legend('E_{0r}/E_{0i}', 'E_{0t}/E_{0i}','Location','best') title('Campo eléctrico, paralelo')
Coeficientes de reflexión y transmisión
El flujo medio de energía por unidad de área en la onda incidente está dada por el valor medio del vector de Poynting. Para un medio no conductor y no magnético μr≈1
Los coeficientes de reflexión R y transmisión T se definen como los flujos por unidad de tiempo y área en la separación, es el vector unitario normal a la superficie de separación.

El campo eléctrico es perpendicular al plano de inicidencia
Comprobamos que R+T=1
Representamos los coeficientes R y T para n1=1 y n2=1.5
n1=1; n2=1.5; hold on R=@(x) (((n1/n2)*cos(x)-cos(asin(n1*sin(x)/n2)))./((n1/n2)*cos(x)+ cos(asin(n1*sin(x)/n2)))).^2; fplot(R,[0,pi/2]) T=@(x) 4*(n1/n2)*(cos(x).*cos(asin(n1*sin(x)/n2)))./((n1/n2)*cos(x)+ cos(asin(n1*sin(x)/n2))).^2; fplot(T,[0,pi/2]) hold off grid on set(gca,'XTick',0:pi/12:pi/2) set(gca,'XTickLabel',{'0','\pi/12','\pi/6','\pi/4','\pi/3','5\pi/12','\pi/2'}) xlabel('\theta_i') ylabel('R,T') legend('R', 'T','Location','best') title('Campo eléctrico, normal')
Representamos los coeficientes R y T para n1=1.5 y n2=1
n1=1.5; n2=1; th_c=asin(n2/n1); %reflexión total hold on R=@(x) (((n1/n2)*cos(x)-cos(asin(n1*sin(x)/n2)))./((n1/n2)*cos(x)+ cos(asin(n1*sin(x)/n2)))).^2; fplot(R,[0,th_c]) T=@(x) 4*(n1/n2)*(cos(x).*cos(asin(n1*sin(x)/n2)))./((n1/n2)*cos(x)+ cos(asin(n1*sin(x)/n2))).^2; fplot(T,[0,th_c]) hold off grid on set(gca,'XTick',0:pi/12:pi/2) set(gca,'XTickLabel',{'0','\pi/12','\pi/6','\pi/4','\pi/3','5\pi/12','\pi/2'}) xlabel('\theta_i') ylabel('R,T') legend('R', 'T','Location','best') title('Campo eléctrico, normal')
El campo eléctrico es paralelo al plano de inicidencia
Comprobamos que R+T=1
Angulo de Brewster
Teniendo en cuenta estas relaciones, comprobamos que
Representamos los coeficientes R y T para n1=1 y n2=1.5. Se señala mediante una línea vertical el ángulo de Brewster
n1=1; n2=1.5; hold on R=@(x) (((n1/n2)*cos(asin(n1*sin(x)/n2))-cos(x))./((n1/n2)* cos(asin(n1*sin(x)/n2))+cos(x))).^2; fplot(R,[0,pi/2]) T=@(x) 4*(n1/n2)*(cos(x).*cos(asin(n1*sin(x)/n2)))./((n1/n2)* cos(asin(n1*sin(x)/n2))+cos(x)).^2; fplot(T,[0,pi/2]) ang_B=atan(n2/n1); %ángulo de Brewster line([ang_B, ang_B],[0,1],'lineStyle','--') hold off grid on set(gca,'XTick',0:pi/12:pi/2) set(gca,'XTickLabel',{'0','\pi/12','\pi/6','\pi/4','\pi/3','5\pi/12','\pi/2'}) xlabel('\theta_i') ylabel('R,T') legend('R', 'T','Location','best') title('Campo eléctrico, paralelo')
Representamos los coeficientes R y T para n1=1.5 y n2=1. Se señala mediante una línea vertical el ángulo de Brewster
n1=1.5; n2=1; th_c=asin(n2/n1); %reflexión total hold on R=@(x) (((n1/n2)*cos(asin(n1*sin(x)/n2))-cos(x))./((n1/n2)* cos(asin(n1*sin(x)/n2))+cos(x))).^2; fplot(R,[0,th_c]) T=@(x) 4*(n1/n2)*(cos(x).*cos(asin(n1*sin(x)/n2)))./((n1/n2)* cos(asin(n1*sin(x)/n2))+cos(x)).^2; fplot(T,[0,th_c]) ang_B=atan(n2/n1); %ángulo de Brewster line([ang_B, ang_B],[0,1],'lineStyle','--') hold off grid on set(gca,'XTick',0:pi/12:pi/2) set(gca,'XTickLabel',{'0','\pi/12','\pi/6','\pi/4','\pi/3','5\pi/12','\pi/2'}) xlabel('\theta_i') ylabel('R,T') legend('R', 'T','Location','best') title('Campo eléctrico, paralelo')
Incidencia normal
Para la incidencia normal θi=θt=0, el plano de incidencia está indeterminado y las dos pares de fórmulas (paralelo y perpendicular) coinciden.
Representamos los coeficientes R y T para valores de n1/n2 comprendidos entre 0.1 y 10
x=logspace(-1,1); R=@(x) ((x-1)./(x+1)).^2; T=@(x) 4*x./(x+1).^2; semilogx(x,R(x),x, T(x)) grid on xlabel('n_1/n_2') ylabel('R,T') legend('R', 'T','Location','best') title('Incidencia normal')
Onda incidente, reflejada y transmitida

Hay muchas maneras de disponer los ejes. El plano de incidencia es el plano XZ y la perpendicular a este plano es el eje Y. Los vectores unitarios, se expresan en este Sistema de Referencia
A continuación, damos las expresiones del campo eléctrico incidente, reflejado y transmitido
Las ecuaciones de Fresnel establecen las relaciones entre las amplitudes de la onda reflejada E0r y la incidente, E0i; entre las amplitudes de la onda transmitida E0t y la incidente, cuando
- El campo eléctrico es normal ⊥ al plano de incidencia
- El campo eléctrico es paralelo ∥ al plano de incidencia
Reflexión total
n1=1.5; n2=1; fill([0, 4.5,4.5,0],[0,0,-1,-1],'c') th_c=asin(n2/n1); for th=(5:5:35)*pi/180 line([0,tan(th)],[-1,0]) th_t=asin(n1*sin(th)/n2); line([tan(th),tan(th)+tan(th_t)],[0,1]) end line([0,tan(th_c)],[-1,0],'color','r') line([tan(th_c), 4],[0,0],'color','r') for th=(45:5:60)*pi/180 line([0,tan(th)],[-1,0]) line([tan(th),2*tan(th)],[0,-1]) end axis off
Cuando el índice de refracción del medio 1 es mayor que el índice de refracción del medio 2, n1>n2, de acuerdo con la ley de la refracción, n1sinθi=n2sinθt, el ángulo θt>θi. Se denomina ángulo crítico θc al ángulo de la ónda incidente que corresponde a θt=π/2. Si el ángulo de la onda incidente es mayor o igual al crítico θi≥θc, solamente hay onda reflejada, no hay onda transmitida al segundo medio
Para θi>θc, cosθt es un número imaginario y solamente tiene sentido, el signo - delante de la unidad imaginaria
El campo eléctrico es perpendicular al plano de incidencia
Analizamos el caso más simple, el campo eléctrico es perpendicular al plano de incidencia
La onda transmitida se propaga a lo largo del eje X en la superficie de separación
siendo λ1 la longitud de onda en el medio 1
Se amortigua exponencialmente en el segundo medio, la amplitud se reduce en un factor 1/e a una distancia ze que se denomina profundidad de penetración
donde λ2=2π·v2/ω es la longitud de onda en el medio 2
Representamos 2π·ze/λ2 en función del ángulo de incidencia θi para n1=1.5 y n2=1
n1=1.5; n2=1; th_c=asin(n2/n1); %ángulo crítico fplot(@(x) 1./sqrt((n1*sin(x)/n2).^2-1),[th_c,pi/2]); grid on set(gca,'XTick',0:pi/12:pi/2) set(gca,'XTickLabel',{'0','\pi/12','\pi/6','\pi/4','\pi/3','5\pi/12','\pi/2'}) xlabel('\theta_i') ylim([0,6]) ylabel('z_e/\lambda_2') title('Profundidad de penetración')
Diferencia de fase entre la onda incidente y reflejada
Las ecuaciones de Fresnel cuando el campo eléctrico es perpendicular al plano de incidencia son ahora cocientes de números complejos
Recordando la forma polar de un número complejo
Onda reflejada
La primera ecuación de Fresnel se expresan
La amplitud de la onda reflejada es igual a la de la onda incidente, el coeficiente de reflexión es la unidad
La energía se refleja totalmente y no existe flujo neto de energía a través de la separación
Representamos el ángulo α en función del ángulo de incidencia, θi≥θc para n1=1.5 y n2=1
n1=1.5; n2=1; th_c=asin(n2/n1); fplot(@(x) 2*atan(sqrt(sin(x).^2-(n2/n1)^2)./cos(x)),[th_c,pi/2]); grid on set(gca,'XTick',0:pi/12:pi/2) set(gca,'XTickLabel',{'0','\pi/12','\pi/6','\pi/4','\pi/3','5\pi/12','\pi/2'}) set(gca,'YTick',0:pi/6:pi) set(gca,'YTickLabel',{'0','\pi/6','\pi/3','\pi/2','2\pi/3','5\pi/6','\pi'}) xlabel('\theta_i') ylabel('\alpha') title('Reflexión total')
El cambio de fase α en la reflexión varía de cero para el ángulo crítico θc a π (180°) para θi=π/2
Onda transmitida
Para la onda transmitida es evidente que E0t no es cero a pesar de que es nulo el flujo neto de energía a través de la superficie de separación.
Referencias
Paul Lorrain, Dale R. Corson. Campos y ondas electromagnéticas. Selecciones Científicas. 1972