计算光波由光密入射到光疏介质时相位变化
7"S|GEs: 0u0Hl% nl clear
E*'O)) BS{">lPmx close all
aH >.o 1; \K)q$E <! n1=1.45,n2=1;
@cXY"hP` b0&dpMgh: theta=0:0.1:90;
D)!k '~a!~F~> a=theta*pi/180;
xAoozDj ] #J ]f rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
S ":-5S6 h.8J6;36 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
BE m%x 0y f^]2q oN tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
&XZ>}^lD^ AU7c =
H:? ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
fNlU c FB</~
g arp=angle(rp);
k r/[|.bq F4:ssy^ ars=angle(rs);
N,;Bl&EU DLz ~$TF^ atp=angle(tp);
0_j! t g;*~ xo ats=angle(ts);
c5]1 aFKz WRNO) f< figure(1)
)>?K:y8I~ <2R=!n@b\ subplot(1,2,1);
z?K+LTf8 iKdC2m plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
M9i u#6P PgxU;N7Y legend('arg(r_p)','arg(r_s)')
Lu <'A4Q1 ~+|p.(I xlabel('\theta_i')
3 =S.- $A-b-`X ylabel('\phi')
|M+ !O93 |m ?0h.O, title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
bS0LjvY9g rv\<Q-uQ8 axis([0 90 -3.5 3.5])
Uyv FR @ _@HMk"A grid on
Q#vur o Z3Vi il: subplot(1,2,2);
J NsK _El= M0 plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
qUVV374N d_OHQ pfK legend('arg(t_p)','arg(t_s)')
hd '!f ni3^J5X W xlabel('\theta_i')
T3Tk:r PDtaL ylabel('\phi')
Ldig/: ?6a:!^eL title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
2. CI^.5& 2"2b\b}my axis([0 90 -3.5 3.5])
5Rc
5/ m xr o grid on
TMq\}k-I5 i)eub`uMy 
