计算光波由光密入射到光疏介质时相位变化
[0D .K}7 | n#_$\
p>Yd clear
%sQ^.` 2 Rr|VD@% close all
XACm[NY_ qPfQy
n1=1.45,n2=1;
85xR2 <: 'Ne@e)s9 theta=0:0.1:90;
a PbE;"
f q.`NtsW!\+ a=theta*pi/180;
y:qUn!3 $(>+VH`l rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
_/5H l` 1 Nd2{( rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
9j:"J` ' c. F6~IHu7 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
cJ
@Wt>YI M} v/tRI ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
W#WV fr Y#$%iF arp=angle(rp);
n3
r3"~i k-OPU , ars=angle(rs);
*>qp:;,DKP A` o8'+`C atp=angle(tp);
Jcm&RI"{ ;)^`3` ats=angle(ts);
0$fpIz |02gup qqi figure(1)
Qtv&ijFC CDR@
`1- subplot(1,2,1);
tQYM &6g j@uOOhy plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
o:P}Wg/NK \WB<86+z legend('arg(r_p)','arg(r_s)')
Z+. '> )6,=f. % xlabel('\theta_i')
ktBj|-'> @x'"~"%7b ylabel('\phi')
@wo(tf=@P J!pygn O title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
9}F*P669f LKOwxF#TKT axis([0 90 -3.5 3.5])
=sJ7=39 v~V!ayn)wQ grid on
8sM|%<$=j '2nhv,|.U subplot(1,2,2);
4UvZ)^r m\r@@! plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
W R y aKM h"_;IUZ! legend('arg(t_p)','arg(t_s)')
p,>5\Zre~ Be=u&T:~ xlabel('\theta_i')
2Mvrey) ^l=!JP=M= ylabel('\phi')
QaS7z#/?. vAh6+K.e title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
(# Gw 1 XJ` ]ga axis([0 90 -3.5 3.5])
dx &'fe*? Raxrb=7 grid on
aM2 l2 `V]5 sE]G 
