计算光波由光密入射到光疏介质时相位变化
AsO)BeUD D!rD-e clear
gQ$ 0 |0O vRMGNz_P7[ close all
cpu|tK.t ,% *Jm n1=1.45,n2=1;
m2v'zJd}g |\T!,~ theta=0:0.1:90;
R47tg&k6[ S9{&.[O a=theta*pi/180;
{)wl`mw3 6ANA oWg* rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
sNs H l ^ 2LqKo\T rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
8I<LZ{a10 R_gON* 9 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
[[ uZC Ki :LLz$[c8 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
xV.UM8 EfqC_,J*3 arp=angle(rp);
^~W s4[Guo Y@MFH>* ars=angle(rs);
UQO?hZ!y/. S4D~`"4 $/ atp=angle(tp);
Xp~O?2:3l V`xE&BI ats=angle(ts);
!yu-MpeG C
A $R
figure(1)
%h}3}p#4 nmoC(| r subplot(1,2,1);
uO4R5F|tL @!mjjeG+1 plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
-~_;9[uV T;#:Y legend('arg(r_p)','arg(r_s)')
Am=O-;
b'8 0pA> w8 mh xlabel('\theta_i')
H iEQs|""' lFD/hz7lc ylabel('\phi')
[-2Tj)P
C v7mg8' title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
[t"#4[ ]]/p.#oD, axis([0 90 -3.5 3.5])
"tz`@3,5dN YLid2aF grid on
%kL ]-Z 28O 3N;a subplot(1,2,2);
w"OeS;#e: c1 h?aP plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
%$=}ePD .`+N+B(4 legend('arg(t_p)','arg(t_s)')
6B@{X^6y !==C@cH<N xlabel('\theta_i')
5mg] su&# E[ tEW0ub ylabel('\phi')
2qP Q3-' 4H hQzVM{ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
6mC% zXR5 P |;=dX#- axis([0 90 -3.5 3.5])
g4 2f*~l _}Ec [c grid on
@^^,VgW[ 'Ug-64f> 
