計算光波由光密入射到光疏介質時相位變化
e4FR)d0x U$WGe >, clear
;Irn{O [Pl''[ close all
Y+eDE:4 ]U&<y8Q_6 n1=1.45,n2=1;
q+
`QiPj \46*4?pP theta=0:0.1:90;
K^IB1U$ Bh7hF?c Sj a=theta*pi/180;
Q]<6voyy tBVtIOm9 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
G +AP."M? _<'?s>(U' rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
P1R[M|Fx `A<2wd; tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
TLq^5,qG 7CwWf ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
r}}2Kl "q
KVGd arp=angle(rp);
JvkTfTE7 LtPaTe ars=angle(rs);
TW}].A_- ZC`VuCg2O atp=angle(tp);
S~)_=4Z 9CAu0N5< ats=angle(ts);
6{I6'+K~ ! F<::fN figure(1)
"\u<\CL Awr(}){ subplot(1,2,1);
s1t