计算光波由光密入射到光疏介质时相位变化
WoB'B|% 0F_hXy@K clear
qa@;S,lp kb"Fw :0
close all
c /^:vTF PZ'|) n1=1.45,n2=1;
FJ!`[.t1AU L;vglS=l; theta=0:0.1:90;
p$x>I3C(\ No[9m_ a=theta*pi/180;
NKB["+S< T|h!06 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
xZ%3e
sp X@Z t4)2# rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
%x@bP6d[ e:4,rfF1 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
*\}$,/m[' ht6}v<x.eA ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
/g9^g( ?6ss SjR} arp=angle(rp);
_99 +Vjy t.RDS2N| ars=angle(rs);
aQY. 96yo >$CNR*}@ atp=angle(tp);
a;U)#*(5|v a_ [+id ats=angle(ts);
bf1$:09 `
-SC,qHw figure(1)
x5smJ__/ hS AI G subplot(1,2,1);
Z[Iej:o5 N.]~%)K:{ plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
aL;zN%Tw Ge?DD,a c legend('arg(r_p)','arg(r_s)')
9fTl 6?x fU%Ys9:wU xlabel('\theta_i')
Gf.xr%mUZr jVGAgR=[ G ylabel('\phi')
%RF$Y=c'C ;QCGl$8A title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Lp }V 94xT Mg8ciV}\xY axis([0 90 -3.5 3.5])
Er$&}9G+- +"WNG grid on
"W4|}plnu S{]3e-? subplot(1,2,2);
L r\(7r pBBKfv plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
n 4XkhY| |pMP- legend('arg(t_p)','arg(t_s)')
P@5-3]m= Y Kp@ n8A xlabel('\theta_i')
G\k & s F 3^q9ll7Op ylabel('\phi')
.) ,9a , 'h~IbP title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
eW3?3l`fvt F8tMZ,: axis([0 90 -3.5 3.5])
aWLA6A+C& 9<P%?Q grid on
D?* du#6 ]0% { IgB 
