计算光波由光密入射到光疏介质时相位变化
Pp8S\%z~h NfgXOLthM clear
QHk\ Z #aP#r4$ close all
}\"EI<$s 7* 5 B n1=1.45,n2=1;
jdxHWkQ q#K{~: theta=0:0.1:90;
_\WR3Q!V A WR : ~{ a=theta*pi/180;
>f]/VaMH{ AjVC{\Ik rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
CY1WT E=s h^Q(A rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
%6m/ve HGQ</5Z tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
}% < ?] sl' 4AK~\ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
oB:7 R^a 11H`WOTQF arp=angle(rp);
-+".ut:R C(N' =-;Kl ars=angle(rs);
V"/.An| `a83RX_\ atp=angle(tp);
yZleots1 |a( KVo ats=angle(ts);
]>n{~4a 02J/=AC5 figure(1)
-$d?e%}# ~MG6evm & subplot(1,2,1);
_{*} )&!M Y)rK'OY' plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
W{6QvQD8 /JD}b[J$ legend('arg(r_p)','arg(r_s)')
|L<JOQ t>bzo6cj xlabel('\theta_i')
PucNu8 _}l(i1o,/ ylabel('\phi')
+b|F_ 4) 8k?iC* title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Vk6 c^/v km %r{ axis([0 90 -3.5 3.5])
*6s_7{; e={O&9Z grid on
2I.FSR_G? ` H'G"V subplot(1,2,2);
W{,fpm pS? D ~0Nb plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
Ia2WBs = VfnL-bDGV legend('arg(t_p)','arg(t_s)')
odP<S. $fb%?n{ xlabel('\theta_i')
P;5)Net1X R`j"iC2 ylabel('\phi')
t^ #1=nK M\1CDU+*Ns title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
xdSMYH{2A N8 sT? axis([0 90 -3.5 3.5])
{v"f){ %['NPs%B grid on
a"( Ws]K 1g;2e##) 
