计算光波由光密入射到光疏介质时相位变化
Oy U [( %}t<,ex(yO clear
'o%IA)sF os= Pr{ close all
BQ&G7V #b FJ6;g=V n1=1.45,n2=1;
O ' 5xPJ A6?+$ Hr theta=0:0.1:90;
B/P E{ / FlWg Tn> a=theta*pi/180;
RbexsBq 5C03)Go3Z rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
H;#3 S< LyEM^d] rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
q7itznQSKc zF+NS]XK tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
]p,sve vo l4 @ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
/}%$fB `# N j8 arp=angle(rp);
K^H{B& b8 _;o)MTw|' ars=angle(rs);
;< aT| 4 )qg cz<p?W atp=angle(tp);
sTn}:A 6 [2]Ti_
>D ats=angle(ts);
O'<5PwhG ( MI8Kkb1d figure(1)
'N- nFc^ B%~hVpm,eM subplot(1,2,1);
'}\#bMeObg Z *9Qeu-N: plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
"OIra 2O 3LxhQVx2 legend('arg(r_p)','arg(r_s)')
X/=*o;" : >";I3S-t xlabel('\theta_i')
.?16w`Y {F Ir|R& ylabel('\phi')
K>!+5A$6i v/W\k.?q/ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
75u 5zD I n[!g axis([0 90 -3.5 3.5])
\Vc-W|e
: h"Bf@3 grid on
x WZ87 {P/5cw subplot(1,2,2);
Hr?_ `: 0j =xWC plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
Gr1WBYK K,ccM[hu| legend('arg(t_p)','arg(t_s)')
j_3X
1 w)k y:C=Ni&," xlabel('\theta_i')
gpIq4Q<
l ~b ylabel('\phi')
NuL.l__W 3RwDI k?>% title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
%K9 9_ Cl3 X9|={ng)g# axis([0 90 -3.5 3.5])
#@ quuiYq yvQRr75 grid on
vz\^Aa
#fv XL10W ^ 
