计算光波由光密入射到光疏介质时相位变化
Gt.'_hf Js tQ{/9bN?P clear
C=+9XfP 0 g$s;;V/8e close all
D[yyFo,z 5z_d$.CIc n1=1.45,n2=1;
8)0]cX @N4~|`?U theta=0:0.1:90;
cR3d& /_,U o^/
#i`) a=theta*pi/180;
3Cj) upc elR'e6Q rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
_4N.]jr5 bKpy?5&> rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
~`AB-0t.u P{9:XSa% tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
4Et(3[P71 V8/d27\ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
wGti |7Tu* ZK]qQrIwy arp=angle(rp);
(S!UnBb& _<c"/B ars=angle(rs);
X d+H()nR vR\E;V atp=angle(tp);
=R=V x/O;8 ^b ats=angle(ts);
K} CgFBk 6X@z(EEL figure(1)
aLh(8 ;$ Be|! S_Y P subplot(1,2,1);
zgGysjV r)|~Rs!y, plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
I:jIChT c6[m'cy legend('arg(r_p)','arg(r_s)')
XN<!.RCw b23A& 1X xlabel('\theta_i')
NAD^10 BsFO]F5mmX ylabel('\phi')
dCoP
qK y ?1ey$SSU] title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
uJ2ZHrJ :i!fPN n axis([0 90 -3.5 3.5])
` &A`&-nc= 6!]@ S|vDX grid on
STnM Bz 7 WVeNO,?ytS subplot(1,2,2);
QG* hQh
o:#jvi84F plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
5 &8BO1V. vt5w(}v( legend('arg(t_p)','arg(t_s)')
'^)'q\v 'k G Uu8 N xlabel('\theta_i')
/
<(|4e {Z/iYHv~#c ylabel('\phi')
'f{13-# X@ II Ap-Y~B title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
t#.}0Te7 (n k g axis([0 90 -3.5 3.5])
~S('\h)1 0cG[<\qT grid on
T)e2IXGN <ci( 5M 
