计算光波由光密入射到光疏介质时相位变化
TE
Z%|5(] b G:\*1T clear
D;R~!3f./b 3F;C{P! close all
91]|4k93 u$d[&|`>_ n1=1.45,n2=1;
aj5HtP- +?*;#=q theta=0:0.1:90;
kSz+UMC-7: n6oOk nCna a=theta*pi/180;
d0>U-. HP8J\` rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
t W UI?\ =dAAb\: rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
- qP)L;n &Gt{9# tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
j.& dHtp nqy*>X` ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Q4cCg7|0 {xZ Y4b2 arp=angle(rp);
&Sd5]r@+ U6&`s%mIa ars=angle(rs);
q=X <QhK $}&a*c> atp=angle(tp);
uz!8=,DFw lAN&d;NU6Z ats=angle(ts);
@ [;'b$T$ L A Crg figure(1)
$\u\ 4 n +@rFbsyJ. subplot(1,2,1);
E*YmHJ:k ,?8a3% plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
MN<LZC% $ dEk#"cvg legend('arg(r_p)','arg(r_s)')
jF?0,g : TT q
xlabel('\theta_i')
8!2)=8|f
fJ*^4 ylabel('\phi')
":/c|! ` Q*`\-8J title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
E>_Rsw * b\"F6TF: axis([0 90 -3.5 3.5])
RW&o3_Ua 6u"wgX]H grid on
^Gd1 T LaJvPOQ subplot(1,2,2);
$ +WXM$N @}q, ';H7 plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
Pl B3"{}0Q ZjxF@`H legend('arg(t_p)','arg(t_s)')
L gF?1? Nw. )O xlabel('\theta_i')
AZ!/{1 Az i*|HN"! ylabel('\phi')
3 P0z$jh"H I"> "> title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
WXx nOLJr "e-Y?_S7R8 axis([0 90 -3.5 3.5])
4 ?BQ&d g"/n95k< grid on
7g*!6-W[ FsqH:I4O 
