计算光波由光密入射到光疏介质时相位变化
x|b52<dLL& Wbq0K6X clear
:F |ll? ` !udU,|N close all
Y>/T+ub =bBV
A0y n1=1.45,n2=1;
Dru iiA !*?|*\B^I theta=0:0.1:90;
|erG cKk F@tfbDO? a=theta*pi/180;
o7r7HmA@ &KYPi'C9!z rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
%eE0a4^". fN_qJm#:$y rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
vg-Ah6BC{ RoFOjCc>D. tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
2q
NA\-0i> i V$TvD+ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
kl~)<,/@ w;{= arp=angle(rp);
6 MD9DqD %lnVzGP ars=angle(rs);
8nu@6 )# 7D KTd^^M atp=angle(tp);
Ud_7>P$a j* ZU}Ss ats=angle(ts);
%/4_|.8u r)
u@,P figure(1)
mDXG~*1 >T#" Im- subplot(1,2,1);
M
ZAz= )- _f1;Hhoa plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
h+ms%tNT uq2C|=M-x\ legend('arg(r_p)','arg(r_s)')
5SMV3~*P n *EGOS xlabel('\theta_i')
B1V+CP3t w:pPd;nz0Y ylabel('\phi')
:jf/$]p LVNA`|> title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
<t&Qa~mA ]6
HR axis([0 90 -3.5 3.5])
`Frr?.3&- |vw], r6 grid on
DDq? 4 S,s#D9NU subplot(1,2,2);
#tIeI6 Qw ? *v*fs0 plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
v:P=t2q /^L <q legend('arg(t_p)','arg(t_s)')
1=h5Z3/fj !Ln 'Mi_B xlabel('\theta_i')
PfsUe,* %>y`VN
D ylabel('\phi')
zsLMRO o3 )<f4 F!?,A title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
sd"eu s f.z(o axis([0 90 -3.5 3.5])
;Dbx5-t [1Aoj| grid on
I)kc[/^j$ [C/{ ru&E 
