计算光波由光密入射到光疏介质时相位变化
'U/X<LCl XRx^4]c clear
1 LUvs~Qu N* NGC!p`N close all
8^2E77s4U >9i%Yuy]( n1=1.45,n2=1;
5`QN<4?% e X3|<Bf theta=0:0.1:90;
^EF VjGM (D~NW*, 9 a=theta*pi/180;
E~K5n2CI z 5q( rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
<C{uodFll l`(pV ;{W rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
> uy(N >'g>CD! tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
R^+,D D wtvtglqV ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
gWLhO|y 5JggU arp=angle(rp);
:9Mqwgk,;3 v~`'!N8 ars=angle(rs);
m!LJK`g A m$!Ex}2 atp=angle(tp);
*DS>#x@3*i OkfnxknZ| ats=angle(ts);
w3E#v&"=Y ] GH_; figure(1)
rcU*6`IWA wG}Rh, subplot(1,2,1);
] 3&BLq 8h'*[-]70u plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
.z} *!
SsfHp legend('arg(r_p)','arg(r_s)')
)7;E,m<:tO m>&:)K}m xlabel('\theta_i')
Gq0Q}[53 F$nc9x[S ylabel('\phi')
:/~TV s^zX9IVnp title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
i=AQ1X\s uB>OS 1= axis([0 90 -3.5 3.5])
7L !$hk >))K%\p
grid on
JSu+ /rI1 l0nm>ps'D subplot(1,2,2);
rJ w
Ws bW?cb5C plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
PCs`aVZ ;q &2$Mb legend('arg(t_p)','arg(t_s)')
%Gc)$z/Wd : @]%n~x xlabel('\theta_i')
i&Xjbcbp 5Zy%Nam'gN ylabel('\phi')
~tB#Q6`nB hzV= 7 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
+Vw]DLWR )m<CmYr2 axis([0 90 -3.5 3.5])
wS @-EcCB :O/QgGZN$ grid on
H}PZJf_E N"-U)d-. 
