计算光波由光密入射到光疏介质时相位变化
n
`&/ D d1AioQ9 clear
`6lr4Kk @R bwD,YC close all
="Ho%*@6 K {|p~B n1=1.45,n2=1;
eJxw) zd7 u%}n w :> theta=0:0.1:90;
D^l%{IG
g!lW u[d a=theta*pi/180;
Cs9.&Y W+UfGk}A rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
i BF|&h(\ 9
Vkb>yFX' rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
#gz
M| n>ULRgiT:o tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
o] y l ;I F Sw\_[^CQ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
piPR=B+ .
$BUw arp=angle(rp);
:~2vJzp@? gp>3I!bo[K ars=angle(rs);
`UD/}j@ sH{4Y-J atp=angle(tp);
-w9pwB &dM.
d! ats=angle(ts);
YC++& Nk c3 jx+Q
figure(1)
OGK}EI |bTPtrT8 subplot(1,2,1);
sDPs
G5q< w,#>G07D plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
zHA!%> %' ED=V8';D legend('arg(r_p)','arg(r_s)')
/]K^
rw[ ^2Fei.?T. xlabel('\theta_i')
oI
}VV6vO #|L8tuWW ylabel('\phi')
z@!`:'ak U{2BVqM title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
ULs \+U */sS`/Lx axis([0 90 -3.5 3.5])
b$N 2z X{5vXT\/y grid on
eD,.~Y#?= GeyvId 03H subplot(1,2,2);
]{3)^axW; }AB, 8n` plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
C 3^JAP !% CWZZ 6u legend('arg(t_p)','arg(t_s)')
jyZ (RB `=4r+ xlabel('\theta_i')
RIlw dt
E{XH? _xo ylabel('\phi')
t h Q_"]+i]s@ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
uG wm
r my(2;IJ#{ axis([0 90 -3.5 3.5])
mWoAO@}Y /)YNs7gR grid on
Sdp&jZY )#?"Gjf~ 
