计算光波由光密入射到光疏介质时相位变化
a+\ Gz AP3SOT3I clear
q7 Uu 8JXF SL%4w< close all
/;HytFP +t.T+`
EG n1=1.45,n2=1;
+tfmBZl^ ?{_dW=AQ1 theta=0:0.1:90;
\jq1F9, -Rjn<bTIy a=theta*pi/180;
z RvYN e" ^1- U\ rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
nkN ]z
^j B+zq!+ HJ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Duptles woR((K] #G tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
oH6(L q'q k qwS/s ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
.fo.mC@a .vpQ3m> arp=angle(rp);
B TcxB h /^/'9}7 ars=angle(rs);
4 D\_[(P '| Q=J) atp=angle(tp);
|C>Yd*E,C )E6m}? H5 ats=angle(ts);
r) X?H qV iky=/- figure(1)
9=3V}]^ M ^o?.Rph|i] subplot(1,2,1);
#B+2qD>E / d6mlQS plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
C:4h 5^*
d4[&+ legend('arg(r_p)','arg(r_s)')
db#y ]>^l BvX!n"QIb xlabel('\theta_i')
krz@1[w-j Fzu"&&>0$ ylabel('\phi')
D"A`b{z z)|56
F7' title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
;Kkn7&'F ;2dhue axis([0 90 -3.5 3.5])
a |z{B b |dX#4Mq^, grid on
$)7-wCl</ 1L l@
ocE subplot(1,2,2);
^{~y+1lt' B*D ` KA plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
x]a>Q), ? FM HK\ legend('arg(t_p)','arg(t_s)')
yRvq3>mU W&bh&KzCW xlabel('\theta_i')
~f% gW e j!C^ ylabel('\phi')
.4- ; _|6 {( title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Z#4? /' GhpH7% s axis([0 90 -3.5 3.5])
C'5i >; vZXyc * grid on
A h)7A|0rT {SROg;vA 
