计算光波由光密入射到光疏介质时相位变化
J)66\h= ~v2_vEu}JX clear
8K{
TRPy JGJQ5zt close all
^oj)#(3C S&9{kt|BI n1=1.45,n2=1;
9Y~A2C C:W}hA! theta=0:0.1:90;
#;H+Kb5O T-eeYw?Yf a=theta*pi/180;
7kHEY5s
" i9_ZK/* rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
nx=Zl:Q} w $pBACX rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
1S <V,9( T0v;8E e tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
JhIgq W2 $ TWt[ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
F.K7w G=(F-U;* arp=angle(rp);
R9~%ORI#; _a^%V9t ars=angle(rs);
nRJcYl~
Y RD<l<+C^~ atp=angle(tp);
0d~>zKho <S%M*j ats=angle(ts);
hQvSh\p 6 I"Q9( figure(1)
,i e84o Kr!8H/Z subplot(1,2,1);
g!i\ AMG? \5cAOBja plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
,vl][MhM 6*le(^y` legend('arg(r_p)','arg(r_s)')
1 f).J Yu`b []W xlabel('\theta_i')
nJNdq`y2 J[ du>1D ylabel('\phi')
+2 !F6"hP AN ZD7v6a title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
'/G.^Zl9 P 'xq+ Q axis([0 90 -3.5 3.5])
weYP^>gH' G BV]7. grid on
ggIz) </ I MpEp}7 subplot(1,2,2);
8lh { R |]w0ytL>(2 plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
.A//Q|ot! I6ffp!^}Y legend('arg(t_p)','arg(t_s)')
/{#1w\ "WV]|
TS"] xlabel('\theta_i')
a`|& rggN bzz=8n ylabel('\phi')
: $lx] 9X/c%:)\= title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
< {$zOF} 1-!u=]JDE axis([0 90 -3.5 3.5])
3@f@4t@5V Zu951+&` grid on
LS}dt?78`V 6lpfk& 
