计算光波由光密入射到光疏介质时相位变化
vgH3<pDiU6
!sQY&* clear
i@)i$i4 aW)-?(6> close all
IEsEdw]aZE p.v0D:@& n1=1.45,n2=1;
`>f6) C- e, vvzs o theta=0:0.1:90;
Rg\z<wPBG eTI%^d| a=theta*pi/180;
G'Q -An%z AN8`7F1 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
B}Z63|/N zZPuha8 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
DC IxRPw dx5#\"KX=, tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
R _W+Ylob Fq~yL!#! ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
J%v=yBC2 p{amC ;cI$ arp=angle(rp);
3G,Oba[$< `#F{Waww' ars=angle(rs);
+Mo9kC "pi=$/RD9 atp=angle(tp);
fi-WZ C%$edEi ats=angle(ts);
A!s`[2 Z q0i(i.h figure(1)
J*yf2& lI5 u)X=Qm) subplot(1,2,1);
H[s(e5 6z 7 E r23Q
plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
gp };D T*8VDY7 legend('arg(r_p)','arg(r_s)')
\\P jKAsh 1n >X[!
8x xlabel('\theta_i')
3 GaQk- a,@]8 r-" ylabel('\phi')
9_L[w\P|4 6]V4muz#c title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
\f]k CB +X#vVD3" axis([0 90 -3.5 3.5])
W]"zctE oX#Q<2z* grid on
3 R=,1< &oP +$;Y subplot(1,2,2);
n Ot&pq7 oRm L
{UDZ plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
K SbK EA 7L`A{L legend('arg(t_p)','arg(t_s)')
A;;fACF8e %q/62f7 ? xlabel('\theta_i')
"Q}#^ h]F |'.*K]Yp ylabel('\phi')
G"-?&)M#a 6LOnU~l, title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
p#01gB ~ZmN44?R axis([0 90 -3.5 3.5])
:8L8q< U chcbd
y>C grid on
gFeO }otm R+2 +-j4 
