计算光波由光密入射到光疏介质时相位变化
3U$fMLx]k Qf]!K6eR clear
,Swa DWNO Wh_c<E}& close all
}!Lr!eALr Qj U"|$ n1=1.45,n2=1;
>C3 9`1
N&.p\T&t theta=0:0.1:90;
e90z(EF?0 \ y",Qq? a=theta*pi/180;
_Z2)e*( ,[#f}|s_ rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
iNSJOS Mv =;+?z! rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
jQ}| ]pj+ c'R|Wyf tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
xII!2. tH(#nx8 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
'~J6 mojE Su #1yw> arp=angle(rp);
rzLl M nQ~L.V ars=angle(rs);
U$ bM:d :tG5~sK atp=angle(tp);
4*X$Jle| S~Q";C[& ats=angle(ts);
"O
"@HVF@ h %S#+t(Bf figure(1)
'L/TaP/3 74#@F{ w subplot(1,2,1);
k<H&4Z)d9 d o7{ plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
$*R9LPpk+ @oNrR$7 legend('arg(r_p)','arg(r_s)')
oZt z"B &KWh5S@w xlabel('\theta_i')
0=&]!WRT w2L)f,X ylabel('\phi')
> w2u w"|c;E1;_ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
d M$S|, H ZT#G:a axis([0 90 -3.5 3.5])
Y ~!@ r_m&Jl@4 grid on
3RUB2c4 PV290 4 subplot(1,2,2);
ntejFy9_ m<4Lo0?nS plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
FC#Q tu~J l ,.;dw legend('arg(t_p)','arg(t_s)')
."O(Ig [ ~fe0 Ba4 xlabel('\theta_i')
f9 $q.a* J:a^'' ylabel('\phi')
}s[/b"%y [>86i title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
m/A N*` V CiSG=obw axis([0 90 -3.5 3.5])
@ 2_&ti >smaR ^m grid on
zqkmsFH{ MNz q}(p 
