计算光波由光密入射到光疏介质时相位变化
}opMf6`w k&f/f clear
Fi;OZ>;a vZ$E
[EG} close all
5|Z8UzL cwtlOg n1=1.45,n2=1;
}fo_"bs@ /4;A.r`; theta=0:0.1:90;
.;ofRx< 98?O[= a=theta*pi/180;
v.>K
)%`# \ )> #`X rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
YN<vOv W$;qh B rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
^P]5@d v A<TYt
M tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
1ZYo-a;) h# Z ,ud_ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
+( af O ~9 J6 5:MaS arp=angle(rp);
QUvSeNSp PN<Vqt W ars=angle(rs);
z\8s |! Pi9?l> atp=angle(tp);
/cUu]#h
iD])E/ ats=angle(ts);
|_7k*:#q: >#l: ]T figure(1)
`"yxmo*0 W+U0Y,N 6 subplot(1,2,1);
XE2rx2k v#<{Y' K plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
yJ0 %6],^g ^ACrWk~UY legend('arg(r_p)','arg(r_s)')
*_#2|96) V mQ' xlabel('\theta_i')
oz,.gP% |mbD q\U ylabel('\phi')
7ZyP BOcD?rrZ0 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
%l a1-r~ r@v t.t0# axis([0 90 -3.5 3.5])
K\8zhY yqL" YD grid on
PUZcb+%]h !.t D.(XP subplot(1,2,2);
`zw % ZnzO] plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
g)TZ/,NQ{ K{`R`SXD legend('arg(t_p)','arg(t_s)')
_`^AgRE 'kY/=*=Q xlabel('\theta_i')
yE,qLiH w3sU& |N ylabel('\phi')
UA2 KY}pz5 -\; 0gnf{J title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
]$9y 7Bhj. n=F
r v*"Z axis([0 90 -3.5 3.5])
IpSWg :b %2qBv grid on
|mHf 7gCX -Q
JP J. 
