计算光波由光密入射到光疏介质时相位变化
G/y5H;<9M 0I-9nuw,^; clear
g>E LGG |Q L81ZbNU?$ close all
Di ,^% U#7#aeI n1=1.45,n2=1;
/wGM#sFH n K1Slg#U theta=0:0.1:90;
D=A&+6B@- F/,NDZN a=theta*pi/180;
V@.Ior}w zH 72'"w rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
F$y$'Rzu_B ozyX$tp rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
( U D nsF H5/6TX72N tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
rK 8lBy:< Fk&c=V;SU ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
ueogaifvB "@^k)d$ arp=angle(rp);
`z}?"BW| +qN>.y!Y ars=angle(rs);
nUaJz Pl xWH.^o," atp=angle(tp);
c8 )DuJ#U zF`0J ats=angle(ts);
q^@Q"J =v :^lI`9'*R figure(1)
etQCzYIhn B9 _ X;c subplot(1,2,1);
i'<[DjMDlm dM.f]-g plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
wA ,6bj RT4x\ &q legend('arg(r_p)','arg(r_s)')
B&M%I:i 1 &jc/*Z" xlabel('\theta_i')
RXMISt3+{y s2p\]|5 ylabel('\phi')
{f_={k G{~J|{t\yz title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
tn\yI! a LG9+GszX 2 axis([0 90 -3.5 3.5])
G6Axs1a @Rze|
T. grid on
*}qWj_RT b<[Or^X
] subplot(1,2,2);
e-/&$Qq )th<,Lo3# plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
2 0h}
[Q( 4/ ~E4"8 legend('arg(t_p)','arg(t_s)')
AEI>\Y H064BM xlabel('\theta_i')
'T;P;:!\ 79rD7D&g ylabel('\phi')
IxN9&xa kOrZv,qFG[ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
wYXQlxd y .&iawz axis([0 90 -3.5 3.5])
23eX;gL J 9nX"Sb grid on
IJp-BTO{V #4NaL 
