计算光波由光密入射到光疏介质时相位变化 FG*r'tC~r
qm o9G
clear 46&/gehr
0s2v'A[\
close all s}% M4
7VF LJr t
n1=1.45,n2=1; 'CkIz"Wd
9 j9TPyC/2
theta=0:0.1:90; :A'y+MnK<
7s{GbU\
a=theta*pi/180; ?m? ::R H
.UY^oR=b{
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); nK%LRcAs
ah+iZ}E%
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); BKjS ,2C
^RtIh-Z.9
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); c|@bwat4
*6DB0X_-}
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); >e [i 5
VZmLS 4E
arp=angle(rp); .+A+|yR
JHTSU q
ars=angle(rs); )C] g ld;8
76h ,]xi
atp=angle(tp); J,y[[CdH`
>_"an~Ss
ats=angle(ts); xRLT=.ir
6.nCV 0xA
figure(1) j yUCH*@
<s<n
subplot(1,2,1); iVq' r4S
! \.pq 2
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) EVC]sUT
GH:jH]u!V
legend('arg(r_p)','arg(r_s)') S8j{V5R'
' =8d?aeF
xlabel('\theta_i') nwRc%C``UK
* kDC liL
ylabel('\phi') 7/@TF/V
/zVOK4BqN+
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) P%&0]FCx
qwgPk9l
axis([0 90 -3.5 3.5]) N.{H,oO `
;LSANr&
grid on :hk5 .[
EC!02S
subplot(1,2,2); ?:I* 8Fj
E GU
0)<
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) Q%tX QP .r
ryUQU^v
legend('arg(t_p)','arg(t_s)') c"V"zg22
BdblLUGK#
xlabel('\theta_i') O$j7i:G'5
RF4vtQC=
ylabel('\phi') ' RYIW /a
YdC6k? tzS
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) Mhf5bN|wQ
]e>w }L(gV
axis([0 90 -3.5 3.5]) `l){!rg8IC
VfC <WVYiZ
grid on }DfshZ0QM
*>'V1b4}
