计算光波由光密入射到光疏介质时相位变化
c. >oe*+ 0l-Ef 1 clear
px}|Mu7z~ mg*qiScfW close all
/f|X(docI .6`9H 1 n1=1.45,n2=1;
joiL{ d`
jjGEj theta=0:0.1:90;
0@H|n^Md# MLRK74D a=theta*pi/180;
">y%iE T>R0T{A rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
wtH?
[>S;) J6L K rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
;#oie<
Vit >g+?Oebgw tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
9983aFam QlO0qbG[ y ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
}j*KcB_ f] J M / arp=angle(rp);
% l ,,_:7{ hvDNz"ec{ ars=angle(rs);
CS==A57I C4~;y hz atp=angle(tp);
} ^/9G17 n&-qaoNl ats=angle(ts);
<;#~l* <iH`rP# figure(1)
-3;*K4z$/ h $L/<3oP6 subplot(1,2,1);
pO ml8SQf L"{JR bh[ plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
UG'9*( * X&0 uI*r legend('arg(r_p)','arg(r_s)')
^!9~Nwn -5I2ga xlabel('\theta_i')
}T%}wdj 'J6
M*vO ylabel('\phi')
++V=s\d7
\9/ b!A title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
~Uxsn@nLr m3mp/g.> axis([0 90 -3.5 3.5])
?)`L$Vr= d"6&AJ5a grid on
/4 vG3 TJ1+g
\ subplot(1,2,2);
^sVX) % _c, '>aH= plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
o9F/y=.r= [BT/~6ovrZ legend('arg(t_p)','arg(t_s)')
/ew
Ukc8, XI9js{p xlabel('\theta_i')
_ B 5gR *{y K
8 ylabel('\phi')
Dg}
Ka7H p~9vP)74u title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
4Rvf C@bm axis([0 90 -3.5 3.5])
T`Mf]s)* k&. Jk
B" grid on
'*?WU_L(g N9vP7 
