计算光波由光密入射到光疏介质时相位变化
F]o&m::/K )y>o;^5' clear
7|Z=#3INw &j =Fx F9o close all
?pSb ,kN}' kS_oj n1=1.45,n2=1;
/g$cQ=c U9^o"vT theta=0:0.1:90;
~*"]XE?M pT3p!/pl3 a=theta*pi/180;
]^aOYtKX #9{ N[t rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
`;KU^dH F<FNZQ@<U rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Mn$w_Z? ZqT8G tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
jw63sn .quui\I3 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
DD 8uG`< w7Fz(`\ arp=angle(rp);
)@lZ~01~d y[QQopy4: ars=angle(rs);
st~
1[in `| nC r atp=angle(tp);
;Q YUi R P`z#tDT^" ats=angle(ts);
( 4:&tm/; Kf=6l#J7 figure(1)
Y:o\qr!Y U|tUX)9O subplot(1,2,1);
]M^ k~Xa 4)- ?1?) plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
o8<0#W@S BY*{j&^ legend('arg(r_p)','arg(r_s)')
kH;DAphk 6( #fGH&[ xlabel('\theta_i')
Q= B>Q k OYF]^uJ ylabel('\phi')
K<k!sh B=<>OYH title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
2 pr #qh8 u.\FNa axis([0 90 -3.5 3.5])
LWH(b s9U "gt-bo., grid on
WG~|sLg pr~%%fCh subplot(1,2,2);
Tt;h? MH wjJ plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
x}^ :Bs+j ?=u/ &3Cw legend('arg(t_p)','arg(t_s)')
Go^a~Sf$ j 3/ I = xlabel('\theta_i')
1gK< d g gu1:%raXd ylabel('\phi')
N@qP}/}8 +,;"?j6<p title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
6[.#B !;9 T [? wbYfW axis([0 90 -3.5 3.5])
a:q>7V|%$ cj[a^ ZH grid on
g3 V
bP S['rfD>9 
