计算光波由光密入射到光疏介质时相位变化
$ Yz &x%Lb D,eJR(5I clear
b,a\`% m} B>?Y(" E close all
p _e-u- 4Uf+t?U9 n1=1.45,n2=1;
g]xZ^M+ l @*/1O)v theta=0:0.1:90;
#|XEBOmsQ V 8&/O)} o a=theta*pi/180;
K23_1-mbe Kc udWW] rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
4Sg!NPuu7& A2&&iL=j/ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
`Ik}X w savz>E & tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
7IJb$af:; &SN$D5U' ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
/&j4I lT %,,h )9 arp=angle(rp);
l3MbCBX2 ;oFaDTX] ars=angle(rs);
@@|E1 'c7 *gz {:}NX atp=angle(tp);
SA>;]6)` ( ap%o\&T; ats=angle(ts);
G^mk<pH ,.`^Wx6F figure(1)
3}:pD]`h w78Ius, subplot(1,2,1);
jPU# {Wo# 1]>KuXd
r plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
cn\& ;5 5v ^/n[5@6H legend('arg(r_p)','arg(r_s)')
ex!^&7Q( 8v^i%Gg xlabel('\theta_i')
b "FsT T (2,iG8 ylabel('\phi')
}BogE$tc "}HQ)54& title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
c?R.SBr,' 8e\v5K9 axis([0 90 -3.5 3.5])
pn $50c 1M;)$m: grid on
\[oHt:$do l.) N subplot(1,2,2);
dY-a,ch"8p zRJy3/> plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
hE6tu' 2mlE;.}8 legend('arg(t_p)','arg(t_s)')
#P9VX5Tg fBLd5 xlabel('\theta_i')
8&UuwZ6i- b<( W}$x ylabel('\phi')
Ez$5wY^J `p ?E{k.N title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
-.*\J|S@g Ys_YjlMIbl axis([0 90 -3.5 3.5])
;Z`)*TRp4 |TUpv*pq grid on
JZXc1R| 9 wwAT@=X*} 
