计算光波由光密入射到光疏介质时相位变化
#"{8Z& Z Q d4T?5 vG clear
19u? ^w `~]ReJ!X% close all
-N(MEzAE u" g p"> n1=1.45,n2=1;
6;pREM+ ksaC[G;}: theta=0:0.1:90;
tp^'W7E `x VA]GR4c a=theta*pi/180;
E}= ,"i @#tSx rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
NLGr=*dq n_Y]iAoc` rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
//@=Q!MW ,AM-cwwT:u tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
0cUt "(] xH[yIfHkG@ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
<,4(3 >js ~ J{[]wi arp=angle(rp);
1T,PC?vr{ ><R.z( 4% ars=angle(rs);
rI+w1';C1 ?(8%SPRk atp=angle(tp);
uE`| 0 yKEF ne8^ ats=angle(ts);
_i@eOqoC zN)\2 figure(1)
8 {]Gh 0+ f\U&M,L\ ' subplot(1,2,1);
;;hyjFGq% }k0-?_Z=1 plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
?lN8~Ze 5q kuK F legend('arg(r_p)','arg(r_s)')
_I-VWDCk jZT :-w xlabel('\theta_i')
X/K)kIi PFy;qk ylabel('\phi')
S5u#g`I] {V% O4/ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
)z235}P
OrEuQ-,i@ axis([0 90 -3.5 3.5])
RrdtU7i3 g+ 1=5g grid on
TT ZxkK Df:7P> subplot(1,2,2);
56SS
>b )Q CM2 plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
l()MYuLNV qJXsf M6 legend('arg(t_p)','arg(t_s)')
x3 6 #x 0O,T=z[+> xlabel('\theta_i')
tb?F}MEe A :bPIXb ylabel('\phi')
-jtC>_/ Sja{$zL+W title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
\y H 3Y 4!~
.6cp3 axis([0 90 -3.5 3.5])
3R%UPT0> JV]u(PL grid on
"R2t&X[9 q
X%v Rf0 
