计算光波由光密入射到光疏介质时相位变化
3rv~ r0 $-]PD`wmY clear
771r(X?Fa '~Gk{'Nx" close all
cNRe > q}7(w$& n1=1.45,n2=1;
6~(iLtd# jowR!rqf theta=0:0.1:90;
[IuF0$w=dj |Q~5TL>b a=theta*pi/180;
}2_ i<4,L 2k""/xMF' rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
PxZMH= rexy*Xv`2p rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
9RN! <`H Xc@%_ 6 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
`wLa .Gzj /H[ !v:U ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
&l1CE1 9< sr,8Qd 0M arp=angle(rp);
W(UrG]J*l aruT eJF ars=angle(rs);
oQ8If$a} +LAj h)m atp=angle(tp);
<</
Le% qw%wy j7 ats=angle(ts);
FiJU
* f0lK ,U@P figure(1)
z}m) u 6P* 2Kg` subplot(1,2,1);
q\6ZmKGnT p= U5qM.O plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
` ]:&h' v/lQ5R1 legend('arg(r_p)','arg(r_s)')
(|<.7K N $3MYr5
xlabel('\theta_i')
UK xeN[fv J d eGQ ylabel('\phi')
s\dF7/ b c 7ury L title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
F45UO%/P u0& dDZ axis([0 90 -3.5 3.5])
BAS3&f A (*2kM| grid on
Fps.F hm ?rn#S8nNx< subplot(1,2,2);
2 r=A' l6EDl0~r plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
Hh1OD?N ) <+c6CM$#}V legend('arg(t_p)','arg(t_s)')
:X6A9jmd e7.!=R{6 xlabel('\theta_i')
; eeu 9_$ =-qv[;%& 6 ylabel('\phi')
P7GF"/ (7g1eEK % title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Q@lJ| qrmJJSJ axis([0 90 -3.5 3.5])
U0:tE>3` +wwK#ocw grid on
7BhRt8FSD+ IuQY~! 
