计算光波由光密入射到光疏介质时相位变化
nms<6kfzL +' 6ea+$ clear
L"zgBB?K6 H2t pP~!G close all
36Fa9P FCc %RR|QY* n1=1.45,n2=1;
aDJj VD aN); P> theta=0:0.1:90;
d)J] Y=j 9p0HFri[ a=theta*pi/180;
"%@=?X8 C
Wl95 g rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
(H|d 3 QH? 2v rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
zdQu%q `[HoxCV3o tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
t2%bHIG} /3: IE%o ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
K d|l\k! MxvxY,~{0 arp=angle(rp);
! 6i '~x_ ars=angle(rs);
TKs@?Q,J ^eT>R,aB atp=angle(tp);
m_O=X8uj"D 5O;oo@A:[ ats=angle(ts);
{] ^%?]e p 7E{es|J figure(1)
F#C 6.`B U3iyuE subplot(1,2,1);
kQiW 5 /r 2.j3:l plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
s/'h LkxI )k<~}wvQ0 legend('arg(r_p)','arg(r_s)')
xg %]\# TKutO0 xlabel('\theta_i')
KKTfxNxJn ykl=KR ylabel('\phi')
Fm+)mmJP DB= cc title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
HC, 0" W g,]5&C T3v axis([0 90 -3.5 3.5])
H"
g& xM"XNT6b grid on
*:\9 T#h Dve5m= subplot(1,2,2);
xB] v V<I${i$]0 plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
c15^<6]g +OP:"Q_# legend('arg(t_p)','arg(t_s)')
D`@U[ ` Sw 'cc4Y~0s xlabel('\theta_i')
T k=3"y+u[ +s 0Bt ' ylabel('\phi')
<(lSNGv5N 'r < B aL title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
u:kY4T+Z G~L#v AY axis([0 90 -3.5 3.5])
<Q ~7a
hF cr;`0 grid on
H<Taf%J T &\br_ 
