计算光波由光密入射到光疏介质时相位变化
jdP )y]c .$%p0Yx+ clear
HU i?\4 ?)xIn)#l s close all
QO/0VB42 _\}'5nmw\
n1=1.45,n2=1;
Xv8-<Ks B7C3r9wj theta=0:0.1:90;
V&M*,#(? pk/#RUfT+ a=theta*pi/180;
G<;~nAo?f0 /\ I6j;$z rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Je1d|1!3 )/H;5 cn rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
`@~e<s`j {ez $kz tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
0bg"Q4 t'eqk#rq ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
n~N>;m P y_Y(Xx3 arp=angle(rp);
Rh!UbEPjC HN5m %R&` ars=angle(rs);
{ "y/;x/ Mt[yY|Ec| atp=angle(tp);
> H!sD\b D yo v}y ats=angle(ts);
[{znwK@ &`"DG$N( figure(1)
YG "Ta|@5 ~ZHjP_5Q subplot(1,2,1);
BQ@7^E[ t&r.Kf9Z\ plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
l& :EKh ~P4C`Q1PT# legend('arg(r_p)','arg(r_s)')
zTz}H* U j#NyNv(jE1 xlabel('\theta_i')
@KN+)q P VF%QM;I[Rc ylabel('\phi')
o`U}u qrO #EzBB*kP
title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
{15j'Qwm <I+k B^ Er axis([0 90 -3.5 3.5])
( "$/sT L\:|95Yq grid on
L*@`i ]jl %xt9k9=vZ subplot(1,2,2);
-JW~_Q[ ]a $6QS plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
?OD$`{1 bLqy7S9x legend('arg(t_p)','arg(t_s)')
Y
},E3< |/Nh# xlabel('\theta_i')
?
SFBUX(p _`$Q6!Z)l ylabel('\phi')
B8& q$QV qmS9*me
{ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
u3T-U_:jSV TdhfX {nk axis([0 90 -3.5 3.5])
3efOgP=L #bOv}1,s grid on
_l
d.Xmvd Kb~nC6yJc 
