计算光波由光密入射到光疏介质时相位变化
}=pOiILvD 7&t~R}&| clear
~x+Y kq0 T#e4": A&x close all
kbq:U8+k ^?Vq L\V5 n1=1.45,n2=1;
8dV=1O$ / oQBiPN+v.3 theta=0:0.1:90;
!d |8'^gc iQj2UTds3 a=theta*pi/180;
G"f du(.@ B}k'@;G rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
9 d] tjT 1HXlHic rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
z6} P j>1 p[;8 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
P7Z<0Dt\} Z]e4pR6 ! ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Uk"Y/Ddm E~69^ cd arp=angle(rp);
_3 KZME =u[rOU{X"W ars=angle(rs);
^OjvL6 A/p .='3bQ(UZ4 atp=angle(tp);
g(aNyn 2:/u2K ats=angle(ts);
8 +L7E- J{^n=X9M0J figure(1)
~e+0c'n\ -Cml0}.O subplot(1,2,1);
){ jqfkL J,`_,T plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
w2K Wa-BO WkcH 5[ legend('arg(r_p)','arg(r_s)')
Flne=ij6g t3Qm-J}wSB xlabel('\theta_i')
U !.~ XT= 5@CpP-W# ylabel('\phi')
sOjF?bCdO GW:\l~ d title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
]zy~@,\ /H.(d 4C axis([0 90 -3.5 3.5])
gRw.AXR a D@}S t :m} grid on
n8
G F8a A/ h pY a subplot(1,2,2);
+i_'gDy$ okVp\RC plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
L;4[ k;5 tu7+LwF7 legend('arg(t_p)','arg(t_s)')
//6^+ -he %
i %ew4 xlabel('\theta_i')
, uKs>T^ Vx}Yl&*D ylabel('\phi')
D8q3TyCj% zZ<ns+ h title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
?' TA!MR _+}-H'7= axis([0 90 -3.5 3.5])
9 $:QLE+t [..,( grid on
R Co eJ| :QxL 9&" 
