计算光波由光密入射到光疏介质时相位变化
!/i{l MH9q ;?.J clear
zFff`]^` c>:wd@w close all
3>`mI8 $t .Una+Z n1=1.45,n2=1;
RF53J yt 9BBmw(M} theta=0:0.1:90;
Z!zF\<r cZU=o\ a=theta*pi/180;
F {4bo$~ > :L@?2), rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
/,&<6c-Q@W qCpp6~]Um rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
9YQb & 1.{z3_S21: tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
O6a<`]F < ?}-$ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
<~'"<HwtK qqr?!vem6 arp=angle(rp);
Pz |>"' s6v ; ars=angle(rs);
cyv`B3} {Y=WW7:Qx atp=angle(tp);
1&evG-#<: u9GQU ats=angle(ts);
A _
N;
2E/"hQ w figure(1)
)E@.!Ut4o 0s3%Kqi[ subplot(1,2,1);
= eXU@B dIa+K?INX plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
(\hx` Yh=> [{<`o5qR legend('arg(r_p)','arg(r_s)')
# 9}D4i.`} b v r^zH,C xlabel('\theta_i')
FR4QUk ukfQ e }I ylabel('\phi')
E +R1 !. AFDq}*2Qb title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
7>Ouqxh21 A gg<tM{yB axis([0 90 -3.5 3.5])
(*nT(Adk 6YLj^w] % grid on
QP^Cx= 3kIN~/<R+7 subplot(1,2,2);
1 s2>C!\ ZaDyg"Tw+ plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
{ MSkHf= B^Nf #XN( legend('arg(t_p)','arg(t_s)')
eJVj uG qL&[K>2z xlabel('\theta_i')
8#
>op6^ H* QIB_ ylabel('\phi')
U; V7 u/{ @+M
/& title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
*:YiimOY" ?'#`
nx(! axis([0 90 -3.5 3.5])
3yY}04[9< D},>mfzF grid on
D>@I+4{p +0%w ;'9z 
