计算光波由光密入射到光疏介质时相位变化 WT d})
s
85X^T]zo
clear uD5i5,q1Hs
me9R nPe:
close all %Cr- cR0
q!7\`>.2:{
n1=1.45,n2=1; WAQv4&xGM
"5e]- u'
theta=0:0.1:90; A+dY~@*a
\myc n/e
a=theta*pi/180; C= Zuy^
& v`kyc
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); : Z.mM5
y"]> Rr
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); n^A=ar.
Pgo5&SQb
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); kBT cN D|
H11 Wb(6Wu
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); Kzmgy14o
x0 #+yP
arp=angle(rp); ?@u
&3/&
mH hm~u
ars=angle(rs); v|e\o~2D`
!P;qc
atp=angle(tp); I?:V EN:
xq]&XlA:ug
ats=angle(ts); &fTCY-W[
zZy>XHR
H
figure(1) FX'W%_f,
*R:nB)(6<
subplot(1,2,1); [R~HhM
B[*i}k%i
plot(theta,arp,'-',theta,ars,'--','LineWidth',2) >s>5k
O
MRa>@Jn??A
legend('arg(r_p)','arg(r_s)') 9hq 7:
> SzTZ3!E
xlabel('\theta_i') A76=^ iw
Oll\T GXP!
ylabel('\phi') v14[G@V~\
bv ] ZUF0
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) 9DT}sCLz:B
q}8R>`Z{
axis([0 90 -3.5 3.5]) Rbl(oj#
9*x9sfCv9
grid on =*'` \}];"
abK/!m[q
subplot(1,2,2); ePF9Vzq
-8TLn l~[
plot(theta,atp,'-',theta,ats,'--','LineWidth',2) )oNomsn
|ST& ,a$(
legend('arg(t_p)','arg(t_s)') M5q7`
}>G
+uH1rF_&@
xlabel('\theta_i') g,1\Gj%y
<;Xj4
J
ylabel('\phi') W[vak F
1]qhQd-u
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) >sP;B5S
O [m+5+
axis([0 90 -3.5 3.5]) x o*a9H?@
rVO+
vhih
grid on I7}[%(~Sf/
T|=8 jt,
