计算光波由光密入射到光疏介质时相位变化 zPt<b!q
r T2Njy1
clear $1f2'_`8~
=2\2S p
close all c^}y9% 4c
C`5'5/-.
n1=1.45,n2=1; 3H2 ~?CaJ
fU>l:BzJ K
theta=0:0.1:90; !XP jRd q
M2Q,&>M
a=theta*pi/180; HP# SR';E
Af 3 |l
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); omE- c
$*q|}Tvl#
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); #&Hi0..y
3 h7RQ:lUi
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); )/RG-L
gAAC>{Wh
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); }g bLWx'iG
fyQOF ItM
arp=angle(rp); F0$w 9p
JFT$1^n
ars=angle(rs); B;k'J:-"
bP18w0>,
atp=angle(tp); * %D_\0;
G6p gG+w
ats=angle(ts); #Y7jNrxE
"Zp&7hI
figure(1) ] Lv3XMa
ddQ +EY@!
subplot(1,2,1); u)oAQ<w
DU^.5f
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) ^{M$S0g|N
X1# Ar)
legend('arg(r_p)','arg(r_s)') c
p"K ?)
~?\U];l
xlabel('\theta_i') ]8A*uyi
$ nt&'Xnv
ylabel('\phi') j
iKHx_9P
Xm#W}Y'
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) "OLg2O^
n xZz{ &
axis([0 90 -3.5 3.5]) Kd:l8%+
3x~7N
grid on
S kj G}
r8 9o
subplot(1,2,2); qPDNDkjDD
N@d~gE&^
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) 5wue2/gl
Skr (C5 T
legend('arg(t_p)','arg(t_s)')
p9"dm{
Ix bQ6
xlabel('\theta_i') }#ink4dK:
x?S86,RW
ylabel('\phi') mG[jR*JW
jg]KE8(
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) ;W-
A2g
[*
<x)
axis([0 90 -3.5 3.5]) ,p|Q/M^
"rB B&l
grid on coXg]bUKo
ZID- ~
6
