计算光波由光密入射到光疏介质时相位变化
l%5%oN`4 =v=! x clear
0!\q x#VUEu]8 close all
cW^) $>A O ;t?@!_ n1=1.45,n2=1;
\k-juF80 =0yJ2[R7Do theta=0:0.1:90;
yC*B OJS
I_s(yO4pw a=theta*pi/180;
ho$%7mc la 7QN QW rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Rw8 m 5U |^{" 2l"j rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
TN_$E&69I V
*y tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
r @UY$z ^#nWgo7{7 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
~y^lNgujO $s!meg@s arp=angle(rp);
5 FKd{V' ZU'^%)6~o~ ars=angle(rs);
C>VZf,JE1 4x=Y9w0?8 atp=angle(tp);
0J</`/g H *lO+^\HXD ats=angle(ts);
hkI);M+@6 0 d]G figure(1)
6^u(PzlA|~ Wy )g449 subplot(1,2,1);
}"k(kH M X\-)e# plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
X0TGJ,yW( H2cc).8 " legend('arg(r_p)','arg(r_s)')
>AFpO*q" 'A2"&6m)28 xlabel('\theta_i')
.w)t<7 y K a|\gl;V ylabel('\phi')
X#< Sv>c^ ^uzVz1%mM title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
ji=po;g=E k@U`? 7X axis([0 90 -3.5 3.5])
_S1uJ~j;E nJg2O@mRJ grid on
Xy&#}S}9 6<nO2 GW subplot(1,2,2);
:h~!#;w_ vn0*KIrX plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
W7"sWaOhW 5Oq ;V: 7 legend('arg(t_p)','arg(t_s)')
h DsORh!i 0yN lf-O xlabel('\theta_i')
$g;xw ?~# ro@BmRMW ylabel('\phi')
k0?6.[ku &nP rozC title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
SiT & p &r%3)Z8Et axis([0 90 -3.5 3.5])
DBDfB b T7'$A!c grid on
^ + <uHd> 9Z
4 R!Q 
