计算光波由光密入射到光疏介质时相位变化
l QxEiDIL *Ag3qnY clear
*ZIX76y<!A Nz:p(X! close all
\H
5t-w= !Wj`U$]; n1=1.45,n2=1;
/#j)GlNp: xl Q]"sm1 theta=0:0.1:90;
L s+zJ1 r{f$n a=theta*pi/180;
# )s
+I2 :lu "14 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
>^SQrB TN<"X :x9 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
sGE %zCB O S1f}< tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
&3SmTg
% z\YLO%M m ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
5Rp mR ErFt5%FN.O arp=angle(rp);
eN`G2eE X9?0`6Li ars=angle(rs);
C3~O6<,Jh FGeKhA 8jT atp=angle(tp);
{REGoe=W% dh1 N/[ ats=angle(ts);
~du U& \ 5Q: %f figure(1)
@'y 8* _ ( B%[NC 6 subplot(1,2,1);
) )t]5Ys%; M !X^2 plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
OGO\u# ?Ss ~!38 legend('arg(r_p)','arg(r_s)')
:ci D!Ly }PM7CZSq xlabel('\theta_i')
q
s: TR x$FcF8 ylabel('\phi')
\jZ)r>US" hZWkw{c title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
F+L%Ho;@P 7Sh1QD YZ axis([0 90 -3.5 3.5])
)`6O SB $GHi9aj_P grid on
8"p rWAN /SyAjZ subplot(1,2,2);
~ _IQ:]k .v-2A );I plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
BkY#wJ' =(~Zm B\ legend('arg(t_p)','arg(t_s)')
EVE"F'Ww,_ :Ys
;)W+R xlabel('\theta_i')
4N=
g l( \rE ] V,,2 ylabel('\phi')
cob??|,\m kq([c r title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
:^`WrcOJ D *IeG>% axis([0 90 -3.5 3.5])
aOQT-C[
O b *3h}n; grid on
Z,"YMUl' -BWWaL 
