计算光波由光密入射到光疏介质时相位变化 ^k7`:@
z0U
vdLBf+Zi
clear ,J3s1 ]~^
<$]=Vaq
close all A=S_5y
nr t3wqJ
n1=1.45,n2=1; dq&N;kk
|
IY='tw
theta=0:0.1:90; }c$@0x;YQ
bm+
#OI
a=theta*pi/180; @{|vW
dO{a!Ca
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); np#RBy
"DniDA
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); SQ_w~'(
d/fg
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); cn~M:LW23
M'
d ,TV[
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); XDWR]
cy6lsJ"?
arp=angle(rp); ph;ds+b
B?cn5
ars=angle(rs); ah(lH5r
X4"D Lt"
atp=angle(tp); /kLG/ry8l:
m1mA:R\zM
ats=angle(ts); 2Bf]#l{z
{O4y Y=G
figure(1) rk$$gXg9/
ZT\=:X*e
subplot(1,2,1); aOj(=s
W.fsW<{4j
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) Vlk]
-f+U:/'.>v
legend('arg(r_p)','arg(r_s)') BO3#*J5S\
jgfl|;I?pg
xlabel('\theta_i') a=m7pe^
zuq7 x7
ylabel('\phi') wstH&^
Vh WF(*
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) )9.i'{{ 0
m3cO{
1I
axis([0 90 -3.5 3.5]) `H_.<``>
|)7dh B
grid on B-@ ]+W
b*7:{FXg
subplot(1,2,2); QL<uQ`>(
w _*|u
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) XpFoSW#K
-27uh
legend('arg(t_p)','arg(t_s)') "`aNNIG&
f^tCD'Vmi
xlabel('\theta_i') &_