计算光波由光密入射到光疏介质时相位变化
=i$Fl{vH p,
h9D_ clear
C].iCxn )B
T close all
xe]y] .8@$\ZRP n1=1.45,n2=1;
IoxgjUa 5yd MMb theta=0:0.1:90;
Bn^0^J- @ju@WY45$^ a=theta*pi/180;
g,Lq)'N;O ^hG-~z< rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
qz+dmef !02`t4Zc- rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
p#&h=,W} 4;w;'3zq tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
LzG%Z1` sGi"rg# ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
riUwBiVa?2 v:f}XK< arp=angle(rp);
jfP*"uUK zpzK>DH( ars=angle(rs);
],>@";9u" 4qO+_!x{) atp=angle(tp);
"~#3&3HVS SOs:]U-T3 ats=angle(ts);
:nh_k4S@v :yL] ;J figure(1)
}K7#Q 1Lc#m`Jln subplot(1,2,1);
ls_'')yp .?D7dyU l1 plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
w;;BSJ]+[ ~d
o9;8v legend('arg(r_p)','arg(r_s)')
V39g,=`b% !-z'2B*:^ xlabel('\theta_i')
Mz9r5 pf2[,v/ ylabel('\phi')
3]OE}[R J; 3{3 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
z#8~iF1 7yG#Z)VE axis([0 90 -3.5 3.5])
N"pc,Q\xU NC8t)
X7 grid on
-0d0t! OPetj.C/a subplot(1,2,2);
ko5 @qNq P#E &|n7DT plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
,QOG!T4 fsmN)_T legend('arg(t_p)','arg(t_s)')
!xz eM VI <vnHz?71c xlabel('\theta_i')
V8e>l[tH sW Qfr$^A ylabel('\phi')
?# Mr d)B@x` title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
bADnW4N`6; 8&;UO{ axis([0 90 -3.5 3.5])
Eou~P h*t gMv.V{vD grid on
efSM`!%j ZWii)0'PV