计算光波由光密入射到光疏介质时相位变化
k 1a?yH)= m] yUcj{F clear
Eg&:yF}? ( C}|.z close all
lcLDCt ? U[2;Fkapi n1=1.45,n2=1;
3 l}9'j ACy}w?D< theta=0:0.1:90;
"TP~TjXf q bIlNA )g a=theta*pi/180;
p>i8aN pN)x,<M) rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
d&AO 4^ c&P/v#U_ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
=1k%T {> q7r b3d tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
5} Id[%.x *~YdL 7f)J ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
\#]C ! JQ <Y6 zJ#BD arp=angle(rp);
$P#x>#+[A $f _C~O ars=angle(rs);
4JU 2x 1Jdx#K atp=angle(tp);
~-[!>1!% @/?i |!6 ats=angle(ts);
,(Ol]W} UmvnVmnv figure(1)
gaxM# xJAQ'ANr subplot(1,2,1);
XI |k,Ko< 8V}| (b# plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
$U. | D/Ok legend('arg(r_p)','arg(r_s)')
.%)uCLZr $ @87Y/_l xlabel('\theta_i')
=PQ4S2Q JV@G 9PT ylabel('\phi')
\+nV~Pi"A $./aK J1B title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
/V{UTMSz & !ds#- axis([0 90 -3.5 3.5])
0;@>jo6,! i-w$-2w grid on
[OTn>/W' -)^vO*b 0 subplot(1,2,2);
u8.Tu7~ 8B &EH + plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
( &Jo.
< @8>bp#x/1 legend('arg(t_p)','arg(t_s)')
ZfX$q\7 B
@QWr; xlabel('\theta_i')
^+CWo @. >qOG^{&x ylabel('\phi')
~-6_-Y| SepwMB4@ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
n[gE[kw $pAJ$0=sw axis([0 90 -3.5 3.5])
GC7 WRA A-:k4] {%P grid on
lq "X_M$ ,dTmI{@O 
