计算光波由光密入射到光疏介质时相位变化
4P-'(4I) Rd@n?qB clear
v"Ud mv " yC&b-y close all
p8}5x 2F ]w&?k:y> n1=1.45,n2=1;
<(6-9(zHa EFT02#F_f theta=0:0.1:90;
R fVV(X hD9' `SQ a=theta*pi/180;
?@,f[ U- jP6oJcZ rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
BSq;R G( 2)j#O rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
eiNk]KXAYX VUg~[ tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
"*H'bz K V+"*A ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
A$9_aqbj ;29X vhS8 arp=angle(rp);
bCac .x#jo *t]&b ;=gE ars=angle(rs);
aaa#/OWQZ +%cr?g atp=angle(tp);
"n2xn%t{ T\b";+!W ats=angle(ts);
4- ?`# (
_ F figure(1)
w|UKMbRMU] lLwQridFXh subplot(1,2,1);
kAe-d /{ #_Um0. plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
7`'fUhB! l4B O@ legend('arg(r_p)','arg(r_s)')
Hh'14n&W pmpn^ZR xlabel('\theta_i')
vM>`CZ pl^"1 Z=* ylabel('\phi')
gm%bxr@X~ S5[RSAbf*t title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
N~K)0RETn AB|VO4-? axis([0 90 -3.5 3.5])
A%2M]];%X fylA 0{ grid on
2
DNzC7}e :tS>D5dz( subplot(1,2,2);
kjKpzdbD jX8 C2}j plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
%dWFg<< | +>9^])K| legend('arg(t_p)','arg(t_s)')
\oZUG =K< I)2
xlabel('\theta_i')
y2hFUq %JH_Nw.P ylabel('\phi')
UFY~D"% / Y`_X@Q title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
:8!3*C-= v]Q _ axis([0 90 -3.5 3.5])
lp%.n= '\ )LMuxj grid on
`]I5WTt*X NCpn^m)Q} 
