计算光波由光密入射到光疏介质时相位变化
5K'EuI) d!8q+FI clear
%fP^Fh CW\o >yh close all
>&^w\"' NJqALm!( n1=1.45,n2=1;
H[ @}ri < k6JB%m\E theta=0:0.1:90;
HKOJkbVZ2^ BT>*xZLpS a=theta*pi/180;
Cc z: NpK+ lNsP wyCoj rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
?P`wLS^; ^%_B'X9 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
q,nj|9z V 'IfM~9'D tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
L*FmJ{Yf ?Tuh22J{Q ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
>qtB27jV ItM?nyA arp=angle(rp);
2IjqT L 5mS Xf"R^ ars=angle(rs);
w2{k0MW 1< \cMY6 atp=angle(tp);
9}}D -&Mc {h9#JMIA ats=angle(ts);
LV8,nTYv E o\|dm. "f figure(1)
nt;A7pI` 0?p_|X'_ subplot(1,2,1);
,6t0w|@-k Fg#*rzA plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
}$qy_Esl gFDP:I/ ` legend('arg(r_p)','arg(r_s)')
?\vJ8H[bD =Rb, `% xlabel('\theta_i')
\XRViG,|5 fiw~"2U ylabel('\phi')
"#4PU5. Nq*\{rb title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
a*SJHBB gEi" m5po axis([0 90 -3.5 3.5])
!]kn=7 4 2m`7uQ grid on
' J-(v < S\S @3 subplot(1,2,2);
m@zxjIwT <[~x]- plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
3N]pN<3@ t4gD*j6J3 legend('arg(t_p)','arg(t_s)')
#mllVQ 4 uNcp0 xlabel('\theta_i')
<EBp X H[>_LYZ8 ylabel('\phi')
}1 _gemlf
.mok.f<G_m title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
c&0IJ7fZG u |QfCwQ axis([0 90 -3.5 3.5])
;OdUH l >qCT grid on
q1?2
U< 2iH ,U 
