计算光波由光密入射到光疏介质时相位变化
i5<Va@ru!s uP =_-ZUW clear
`;qv} &/hr-5k close all
A Wh* <H Fl+tbF n1=1.45,n2=1;
mYjiiql~ WUW b 5xA theta=0:0.1:90;
q\b
?o!# _ d)GkXl l1D a=theta*pi/180;
e=]oh$] jbhJ;c : rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Go+xL/f 3Ra\2(bR rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
G2+)R^FSC 8P<UO tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
YV3TxvX MR d! QD vO ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
<#p|z`N .FbZVY c] arp=angle(rp);
gj egzKU mf]( 3ZL ars=angle(rs);
@2~O^5[> aC8,Y$>?E` atp=angle(tp);
}M9DqZ;I qgC-@I ats=angle(ts);
[%iUg\'7d KG-k$glD figure(1)
@Uj _+c
q Z:o
86~su subplot(1,2,1);
;c
m wh< ,F?~'- K plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
ok>P [
&! hKH$AEHEU} legend('arg(r_p)','arg(r_s)')
AeY$.b 4j@kMe;RjZ xlabel('\theta_i')
Dmm r]~ o9T@u Wh+ ylabel('\phi')
H7tQ# @RoRNat title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
.i3lG(
YG <(bCz>o| axis([0 90 -3.5 3.5])
*t?~)o 7 LA;V}%y ? grid on
#lyM+.T ^iV`g?z subplot(1,2,2);
.e FOf V) KT
lP:pB; plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
+J o 3rX'` ez g^5o; legend('arg(t_p)','arg(t_s)')
]w_ i!.I;@ xlabel('\theta_i')
G>[
NZE XAQ\OX# ylabel('\phi')
X@ bn?? _d\u!giy title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
43{_Y] ^EC)~HP@C axis([0 90 -3.5 3.5])
;e~{TkD }1a <{& grid on
,# P,B ;r~ #Cg}!3 8 
