计算光波由光密入射到光疏介质时相位变化
-|$* l
Q ;HT0w_, clear
hF9B?@n?B o8mo=V4j close all
cnC_#kp `lvh\[3^ n1=1.45,n2=1;
\cFAxL( &F86SrsI theta=0:0.1:90;
qY# m*R j@_nI~7f} a=theta*pi/180;
zAu}hVcW O<Jwaap rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
B_b8r7Vn` 8QPT\~ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
@<VG8{ cyDiA(ot& tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Za34/ro/T ^]KIgGv\ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
M'b:B*>6 JPHUmv6 arp=angle(rp);
_y|[Z; M2a}x+5' ars=angle(rs);
-.^@9
a> d!w1t=2H atp=angle(tp);
;;D%
l^m+ ,7QBJ_-;QJ ats=angle(ts);
Bfu/9ad T27:"LVw figure(1)
rrW! X q jw%fN!? subplot(1,2,1);
2f s9JP{^0 WvIK=fdZ$ plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
fYv ;TV>73 32TP Mk legend('arg(r_p)','arg(r_s)')
iM4mkCdOO aS3P(s L xlabel('\theta_i')
[>Ikitow p<<6}3~ ylabel('\phi')
YNuewD 3AX?B~s title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
f|;HS!$ Y}85J:q] axis([0 90 -3.5 3.5])
Qf/j: e-&L\M grid on
B/kcb(5v :zRB)hd subplot(1,2,2);
k?!TjBKm ")fOup@ ^a plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
,^UcRZ8.H 3 (Gygq# legend('arg(t_p)','arg(t_s)')
x(exx
)w 1uK)1%vK xlabel('\theta_i')
JDIz28 Ww { mK pD ylabel('\phi')
cL-6M^!a s0x@
u title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
OJh+[bf" 7{?lEQ&UE axis([0 90 -3.5 3.5])
Y}eZPG.h 8^&fZL', grid on
\ x>NB eh8<?(eK