计算光波由光密入射到光疏介质时相位变化
Wf<L R3 >MK98(F clear
a> )f=uS i&k7-< close all
a6H%5N 6,uX,X5 n1=1.45,n2=1;
qVPeB,kIz {| \.i theta=0:0.1:90;
Mq8L0%j Ha ]YJ} a=theta*pi/180;
+O5hH8<&b ,
dp0;nkr rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
iCyf Oh Ewm9\qmg rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
zT[!o
j7 <ih[TtZ tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
<1%$Vq dTtSUA|V7" ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
b6 M 8V(pugJ arp=angle(rp);
MW{8VH6+ `W-Fssu ars=angle(rs);
.m AjfP* o Rzi> rr atp=angle(tp);
.~}1+ \~5 j7c3(*Pl ats=angle(ts);
y? :.;%!E JCaOK2XT; figure(1)
:Yks|VJ1 j=J/x:w_e subplot(1,2,1);
;>YzEo ,(4K4pN plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
2 m[<]$ gNhQD*+>{ legend('arg(r_p)','arg(r_s)')
m 6djeOl 4a Y|TN/| xlabel('\theta_i')
l +R+&b^ XrPfotj1 ylabel('\phi')
q]M0md _ y>~
yZx title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
8*fv' k'"%.7 $U! axis([0 90 -3.5 3.5])
9S -9.mvop O<\@~U grid on
&M '*6A 2Gdd*=4z subplot(1,2,2);
c=+!>Z&i$G ^VACf|0 plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
rt| 7h>RQ F/A|(AH' legend('arg(t_p)','arg(t_s)')
ow#1="G,= ,=:D xlabel('\theta_i')
Bnd [X ~= -RK$= ylabel('\phi')
+@:x!q|^ Ki~1q u: title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
VQ{f ne< ,u
g@f-T axis([0 90 -3.5 3.5])
2>H24F : \}(&
> grid on
9$m|'$p3sG z"4~P3>{g 
