計(jì)算光波由光密入射到光疏介質(zhì)時(shí)相位變化 }H2#H7!H  
H0+:XF\M  
clear xU:4Y0y8  
;AJ< LC  
close all rei5{PC  
#<]Iz'\`  
n1=1.45,n2=1; :$2Yg[Zc3  
\j:AR4  
theta=0:0.1:90; 7*MU2gb  
vzcz<i )  
a=theta*pi/180; o X@nP?\  
<5k&)EoT  
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); Z)5klg$c  
l#ZyB|  
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); SBz/VQ  
*l{4lu  
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); d EIa=e|  
3A d*
,>!  
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); >Fp&8p`am  
1\IZcJ {  
arp=angle(rp); 7jT#BWt  
zsQF,7/}B  
ars=angle(rs); kU1 %f o  

Yw1Y-M  
atp=angle(tp); F^]aC98]1  
"*t6t4/Q  
ats=angle(ts); ?uig04@3  
V(Dj
F=8  
figure(1) 80HEAv,O  
) G|"jFP  
subplot(1,2,1); hR0]8l|  
PF`:1;PU  
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) hwp/jO:7\  
61kO1,Uz*  
legend('arg(r_p)','arg(r_s)') ?;fv!'?%  
@BXaA0F4  
xlabel('\theta_i') s<,"Hsh^CR  
>pnz_MQ  
ylabel('\phi') ]a=l^Pc(xN  
Rd0?zEKV  
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) T%w(P ^qk  
t9]r  
axis([0 90 -3.5 3.5]) 0OEtU5lf`y  
F=VoFmF@  
grid on ONNW.xHp  
g{]e