计算光波由光密入射到光疏介质时相位变化
WCa>~dF> YH_mWN\Wu clear
kzMul<>sl ixK& E#
close all
h$F.(N IYe RQaB _bg7 n1=1.45,n2=1;
jO` b&]0 8KxBN)fO; theta=0:0.1:90;
,'FH[2 ^,P#
<,D, a=theta*pi/180;
uc.dtq! uc0 1{t0, rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
HR.^
y$IE :5Y
yI.T rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
7(ni _|$| E5^P*6c ( tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
)@vhq Vv? z%lu% ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
QNH-b9u>8 DC&3=Nd arp=angle(rp);
Z&f@)j *=sMJY9#jE ars=angle(rs);
}01c7/DRP< A* 1-2 atp=angle(tp);
]l T8Z-h@ we4k VAn ats=angle(ts);
K;:_UJ>t ^M:Y$9r_s figure(1)
D d: T FZo iy<|<*s2D subplot(1,2,1);
y4 @zi "G Y/%(4q*' plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
q ocN:Of1 q
< , b legend('arg(r_p)','arg(r_s)')
(D.B'V#> cO8':P5Q xlabel('\theta_i')
e;|:W A $#ve^.VHv ylabel('\phi')
9<e%('@[ ~CtLSyB title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
#[
H4`hZ (6y[,lYH axis([0 90 -3.5 3.5])
} @)r\t4m (G>S` B grid on
I pp#{'Do '-,$@l# subplot(1,2,2);
[NxC7p:Lo KQf=t0Z=Ce plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
J]!&E~Y s6DmZ ^Y% legend('arg(t_p)','arg(t_s)')
Kl' u jk Q%b.a xlabel('\theta_i')
7J[DD5 7 R4t%^F ylabel('\phi')
^v j} azs lNL title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
!J6;F}Pd/ {
R`"Nk axis([0 90 -3.5 3.5])
\&6^c= 2= 3)6TnY/u6{ grid on
=O1py_m y6hb-:
#1 
