计算光波由光密入射到光疏介质时相位变化
.F$AmVTN V# %spW clear
%Yn )t3d )h!l%72 close all
1X}Tp\e wLp
t2 b8S n1=1.45,n2=1;
9@a;1Wr/f ,^Srd 20 theta=0:0.1:90;
3jVm[c5%] N~KRwsD H a=theta*pi/180;
^ "#rDP"v ) CTM rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
>43yty\
~F6gF7] z rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
?B!ZqJ # ^W05Z!} tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
JX<W[P>M e %Xf*64 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
;I@L d iL l>z arp=angle(rp);
k0uwG'(z9 TRok4uc ars=angle(rs);
:P 1c>:j[ #iDFGkK / atp=angle(tp);
>U/g*[> OYy !4Fp ats=angle(ts);
Czxrn2p/ D*Q#G/TF3 figure(1)
4N0nU _ 7.y4zQJ subplot(1,2,1);
0l@+xS; i)3\jO0&GU plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
oA%[x i?=.;
0[| legend('arg(r_p)','arg(r_s)')
`\0a5UFR "uC*B4` xlabel('\theta_i')
B9-[wg#0G 04d$_1:}a ylabel('\phi')
cN >i3}fq e .^9&Fk"N title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
X5<.%@Z xB+H7Ya axis([0 90 -3.5 3.5])
iCKwd 9?) 1hS~!r'qqv grid on
c/U6K
yiK ,4,c-
subplot(1,2,2);
I !O S&8:u !l^A Kn| plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
kR(=VM JU w5jZI|
legend('arg(t_p)','arg(t_s)')
+}/!yQtH mkA|gM[g7 xlabel('\theta_i')
O+j:L J,Ap9HJt ylabel('\phi')
LL
[>Uu?Y VC7F#a*V title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
J@iN':l- +cH(nZ* f axis([0 90 -3.5 3.5])
ii9/ UtIQ `p|vutk)U grid on
s+- aHn xrnH= >.;m 
