计算光波由光密入射到光疏介质时相位变化
5?I]\Tb =o HJ_ clear
"3>#[o p,1 RRbyc close all
tr'95'5W. dm 2_Fj n1=1.45,n2=1;
RE7[bM3a uV\=EDno theta=0:0.1:90;
v~xG *e iFDQnt
[t a=theta*pi/180;
(>Yii_Cd k1cBMDSokO rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
X F40;urm <T&$1 m{ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
~ \kRW6 SB F3\ tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
:g\qj? o aY?}4Bx ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
;x7S Y;0* LS_QoS arp=angle(rp);
']rh0? R i3m438 ars=angle(rs);
v
EX <9 x
D r^&rC atp=angle(tp);
~/Ry=8 hQ\]vp7V ats=angle(ts);
jjbw.n+1 JBg>E 3*N figure(1)
$i2gOz ZcQm(my subplot(1,2,1);
npe *A CkflEm fe plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
)a6i8b3 8<G@s`* legend('arg(r_p)','arg(r_s)')
yZ P+ kjmF-\ xlabel('\theta_i')
QJTC @o "8rP?B( ylabel('\phi')
-y AIrvO1q Ue\& title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
7oWv' C/w;g3 axis([0 90 -3.5 3.5])
/bb4nM_E/ LRI_s>7 grid on
ML]?`qv ' 0O:TKgb&C. subplot(1,2,2);
b9v Kux W%wS+3Q/ plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
X >**M z/ i3 legend('arg(t_p)','arg(t_s)')
2<O
hO
^ j C@^/rMh xlabel('\theta_i')
b6i0_fOO *oPSkEA{ ylabel('\phi')
vxm`[s |QC C$hsR& title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
'=ZE*nGC X}Csl~W8in axis([0 90 -3.5 3.5])
O{4 m-; UFl*^j_)] grid on
"K@os< hb_Ia]b 
