计算光波由光密入射到光疏介质时相位变化 ',Y`XP"Q
d^uE 4F}
clear |,YyuCQcL[
C8ss6+k&
close all r 0m A
: jgvg$fd
n1=1.45,n2=1; /G[2
`D
*U@iJ
theta=0:0.1:90; R<\5 q%@G
}ACWSk WK
a=theta*pi/180; GJTKqr|1O
+]?/c>M
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); sA^_I6>M"
{7EnM1]
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
kdmmfw
!+ bLh W`
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); [\i0@
h>0R!Rl8
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); Y9}5 &#
dP[vXhc
arp=angle(rp); Tt ^PiaS!
IZ9L
;"}
ars=angle(rs); }G{ 'Rb
o`tOnwt
atp=angle(tp); :eHD{=
XSfl'Fll D
ats=angle(ts); ;i;2cq
?W Vp,vP
figure(1) wl^7.IR
(w1M\yodV
subplot(1,2,1); v w.rkAGY
Kp]\r-5UD>
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) NY(z 3G
*s=jKV#
legend('arg(r_p)','arg(r_s)') G` ;YB
Wi;w u*
xlabel('\theta_i') -'!K("
3y# U |&]{
ylabel('\phi') a9f!f %9
to2#PXf]y
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) aLo^f= S
CIt%7
\c
axis([0 90 -3.5 3.5]) ?cyBF*o
>]g hme
grid on \zx &5a
#
}q$ 6^y
subplot(1,2,2); K3Sa6"U
m2{DLw".
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) v 0aV>-v
p1BMQ?=($
legend('arg(t_p)','arg(t_s)') ]J '#KT{
a+-X\qN
xlabel('\theta_i') v47S9Vm+
B@+&?%ub:
ylabel('\phi') z~m{'O`
KfPYH\ 0
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) eb#yCDIC
_ArN[]Z
axis([0 90 -3.5 3.5]) Y5/SbQYf1
%&Fsk]T%:
grid on hx.ln6=4
^ ,cwm:B@
