计算光波由光密入射到光疏介质时相位变化
k \|Hd"T fDf:Jec`[ clear
'F1NBL UW!!! close all
\xw E4K '~2;WF0h n1=1.45,n2=1;
w[t!? (![> p "u5wJ_ theta=0:0.1:90;
Z':}ZXy] .xS}/^8iD a=theta*pi/180;
K[O'@v npZ=x-ce rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
b k 30d ULj'DzlfH rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
ex1b j M7 GHfsq|*j,Z tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Qn_*(CSp d
H]'&&M ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
"*a^_tsT?i Mz_*`lR N arp=angle(rp);
ZBi|B D j' g': U ars=angle(rs);
N^H~VG&D( 4fgA3% atp=angle(tp);
se2ay_<F+ a!vF;J-Zqa ats=angle(ts);
q3n(Z {dTtYL$ '" figure(1)
>8\EdN59{ Q 0s!]Dk subplot(1,2,1);
|p}qK
Fdi 0x*L"HD plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
0P_qtS 2Ck'A0d legend('arg(r_p)','arg(r_s)')
1--_E,Su> \L Gj]mb1 xlabel('\theta_i')
: _X9x{ gmd-$%" ylabel('\phi')
g P:mZ7 ~RgO9p(dY title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
wGr5V! T*e>_\Tx axis([0 90 -3.5 3.5])
2|8e7q: +* n$&xVaF| grid on
D-J G0.@ ~H`~& ? subplot(1,2,2);
)%FRBO] 6qp2C]9= plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
?B7n,!&~ kwF4I )6 legend('arg(t_p)','arg(t_s)')
9M8 n ,e<(8@BBL xlabel('\theta_i')
qn=~4rg]R X]
cI ? ylabel('\phi')
'@t} 8J [Vzp D 4 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
n o ).70K 4'9yM XR axis([0 90 -3.5 3.5])
v 9\2/B XqX6UEVR4 grid on
" -bsWC y(!J8( yA 
