计算光波由光密入射到光疏介质时相位变化
D$/*Z5Z)] H/BU2s a clear
t9+ME| h;`]rK;g close all
z\ r|5Z [I4&E > n1=1.45,n2=1;
Ialbz\;F2% Gw0MDV&[ theta=0:0.1:90;
;'xd8Jf BHEZ<K[U
a=theta*pi/180;
8
!+eq5S3 A3c&VT6Q rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
g
xf|L>= 7(< z= F rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
*m 9,_~t Nw|m"VLb tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
xXm:S{I {rWFgn4Li ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
#E #.`/4 @N,I}_ 9- arp=angle(rp);
#Hy fj j _/%,ZoZ2 ars=angle(rs);
'q9='TOk cri.kr9Y atp=angle(tp);
zM W [Xx! GUqhm$6a ats=angle(ts);
N>'|fNx] d:1TSJff%/ figure(1)
_;yi/)-2 dL IZ)16& subplot(1,2,1);
6pi^ rpo E]26a,^L plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
QwL'5ws{q K%/:V legend('arg(r_p)','arg(r_s)')
@Nk]f [r"`r Bw xlabel('\theta_i')
\"^ .>+ T($6L7 j9 ylabel('\phi')
L4C_qb k;: "8|a4Y+F title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
U*7Yi-"/* %.3] F2_Q axis([0 90 -3.5 3.5])
>:yU bo) 0zvA>4cq) grid on
"Ooc;xD3< e$+/;MRq subplot(1,2,2);
eX9Hwq4X44 k yI -nE plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
DHn u F@M _>"f&nb O legend('arg(t_p)','arg(t_s)')
95%,
8t (iw)C)t*u xlabel('\theta_i')
Z 71.* Q+]9Glz9 ylabel('\phi')
X_nxC6[m% -g;cg7O#( title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
"2~%-;c S66. .sa axis([0 90 -3.5 3.5])
!Hp H *MBu5
+u%e grid on
`[*n UdG C0b OPn 
