计算光波由光密入射到光疏介质时相位变化
\(a9 rZ9 JTw'ecFev clear
=?`5n|A* SWjOJjn close all
gY],
(*v 3z[yKua\ n1=1.45,n2=1;
8c m,G Ve, _;<F]S theta=0:0.1:90;
s D] W/ *v'&i) J a=theta*pi/180;
H!Y`?Rc rO%
|PRP rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
rl^_RI ?-,v0# rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
P-L<D!25 5|ih>? C/( tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
b;& J2:` =: U63 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
h:[PO6GdX D-.XSIEMu arp=angle(rp);
b|T }mn A 6'G%of
ars=angle(rs);
?op6_a-wm G 3P3 atp=angle(tp);
"Hmo`E B0 w$H=GF?" ats=angle(ts);
^f{+p*i}: }-m/
'Q figure(1)
@3bQ2jn M:OY8=V subplot(1,2,1);
[{_JO+)+n S%b7NK plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
3OV#H% ^3QHB1I legend('arg(r_p)','arg(r_s)')
H~+D2A hq/k}Y xlabel('\theta_i')
+ )?1F VAc -RaA ylabel('\phi')
+nRO< A'qJke= title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
/Cr0jWu
_ JoA^9AYhR axis([0 90 -3.5 3.5])
5t R<aIf #rzq9}9tB grid on
mFOuE 5 i /*&; subplot(1,2,2);
P y6c=&* !z X`M1J plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
Yt2_*K@rC H E+y1f] legend('arg(t_p)','arg(t_s)')
38 Lc|w n H*U xlabel('\theta_i')
Mw7 ~:O`
m2%uGqz ylabel('\phi')
.@.,D% 7< <FmrYwt title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
d c_^ t#oY|G3O} axis([0 90 -3.5 3.5])
!l@IG C { ) 4D1 grid on
O'U0Y8HN {hK$6bD 3^ 
