计算光波由光密入射到光疏介质时相位变化
\K{ 0L a}BY ov clear
&0OG*}gi hOu3 bA close all
.9 on@S uk< 4+x,2) n1=1.45,n2=1;
+,T RfP
Fb -aPg#ub theta=0:0.1:90;
|mdVdD~go h5{'Q$Erl a=theta*pi/180;
<;eW=HT+uq j^j1 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
hYT 0l$Ng nA-.mWD_C rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
H1pO!>M [fya)} tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Xtq_y'I c)TPM/>(p ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
dUeN*Nq&(, <ktrPlNuM arp=angle(rp);
B 4c]}r+ N=T< _`$5 ars=angle(rs);
JIEK*ui =r?hg GWe atp=angle(tp);
b,l$1{ 0U(@ = 7V ats=angle(ts);
G\/zkrxmv ^EtMxF@D figure(1)
VuhGx :Xl =\&;Fi] subplot(1,2,1);
6"L cJ%o -j# 2}[J7 plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
1y4|{ 7bb )0. kv2o. legend('arg(r_p)','arg(r_s)')
b$d;Q x !=*g@mg F xlabel('\theta_i')
r_)' Ps xBT hq?N ? ylabel('\phi')
0r QMLx :KSV4>X[%a title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
<$ A !1jBC. G1 axis([0 90 -3.5 3.5])
Q Uwd [ #px+;k 5 grid on
1E[J%Rh\ l -/wtI subplot(1,2,2);
[N-Di" }Sm(]y plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
s [RAHU e/KDw legend('arg(t_p)','arg(t_s)')
2RVN\?s: {g'(~ qv xlabel('\theta_i')
IA fc T!{ g+8OekzB5 ylabel('\phi')
[SjqOTon{ +xh`Q=A title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
i^/ T MD}w Y><C axis([0 90 -3.5 3.5])
}kw#7m54 EKYY6S2 grid on
9a[9i}_ yJ[0WY8<kC 
