计算光波由光密入射到光疏介质时相位变化
WRWcB Zk4( clear
o5 fXe}pl@ -5Km 9X8 close all
J6) &b7 U/|H%b n1=1.45,n2=1;
^uy2qO4Yw &Vtgh3I theta=0:0.1:90;
y+C.2 ca N` DLIv8i; a=theta*pi/180;
l5m5H ,` .skR4f,h rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
t!g9,xG<X 3>`CZ]ip} rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Y\qiYra ? 8d7/KZO tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
DM-8azq $ BE;J/ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
XQ{G ) >}43xIRRCq arp=angle(rp);
eDm~B (G$ f4s^$Q{Q ars=angle(rs);
O>E2G]K]\ 6PS #Zydb atp=angle(tp);
e$E~@{[1) 1 D0_k ats=angle(ts);
">n38:?R 7I6bZ;}d figure(1)
yW$ja|^ E R|% R-J] subplot(1,2,1);
z|Hc=AU8y ;
)Vro plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
N
F[v/S r|XNS>V ,$ legend('arg(r_p)','arg(r_s)')
?
[?{X~uq ZD9UE3- xlabel('\theta_i')
IeqJ>t: o@`&
h}
$ ylabel('\phi')
?I` BbT} y@;4F n/ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
h^P>,dy0 OXC7
m axis([0 90 -3.5 3.5])
BX|+"AeF ` L(AvSR grid on
Jlw
oSe:S +w-UK[p subplot(1,2,2);
O[3 AI^2 1NO<K` plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
|H.(?!nTb TqNEU<S/t legend('arg(t_p)','arg(t_s)')
XelY?Ph,, q/,W'lQ\; xlabel('\theta_i')
z.g'8#@ DT # 1*&- ylabel('\phi')
K1Tq7/N "t+VF 4r title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
m:<cLc :. LSQ 2pB2V axis([0 90 -3.5 3.5])
,TD@s$2x tvptaw A. grid on
+XU*NAD,! .J3lo: 
