计算光波由光密入射到光疏介质时相位变化
Dk`4bYK +^*5${g;@H clear
ZS XRzH~0 1 /2V.:bg close all
KR=d"t Qw eOehgU5x n1=1.45,n2=1;
K3rBl!7v S\9t4Ki_' theta=0:0.1:90;
`~ _H=l9{ " J
pTE \/ a=theta*pi/180;
A~SL5h i
Y*o;z,~ rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
ypD< 2z^ %{7_E*I@n rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
"Ap$ Jl B 1l+kO,X] tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
n @ &"+ ]3t1=+ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
v7;J%9=0D` )'[x)q arp=angle(rp);
Q(|PZn g S jVsF1d_ ars=angle(rs);
|rHG%VnBH 7m:|u*ij2~ atp=angle(tp);
7uPZuXHxcu a)!![X?\ ats=angle(ts);
.tRr?*V|l YlswSQ figure(1)
)MWUS;O< ME$2P!o subplot(1,2,1);
;)(Sd f[P Xhcn] plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
7/H^<%;y )jyq{ Jb legend('arg(r_p)','arg(r_s)')
8 !4~T,9G K8HIuQ!= xlabel('\theta_i')
w9RF2J d/D,P=j" ylabel('\phi')
;n`R\NO9 lPM3}52Xu title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
I|,pE**T }p~OCW! axis([0 90 -3.5 3.5])
$hkq>i \ .sM <6; grid on
X'fuF2owd Z_H?WGO subplot(1,2,2);
M.>^{n$
z v(DwU! plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
<9P4}`%)3 r&}(9Cq&"y legend('arg(t_p)','arg(t_s)')
,E\h !/X ~xGoJrF\ xlabel('\theta_i')
t,
U)
~wi >+cSPN'i> ylabel('\phi')
jY7=m Ad B:4Ka]{YO title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
!fY7"E{%% [SVhtrx|% axis([0 90 -3.5 3.5])
l)4O . * @;vNX*-J grid on
|%
z ^N* J/2j;,8D 
