计算光波由光密入射到光疏介质时相位变化
zDB" r UgR : qjI clear
L x.jrF|& ?mMW*ico close all
L? 8^a G Cc!J1) n1=1.45,n2=1;
8.yC A Tr%F U i theta=0:0.1:90;
)\>r-g$ JdiP>KXV a=theta*pi/180;
*hF^fxLb l =p q:m rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
b,Ke>.m xdZ<|
vMR rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
oSY7IIf%L y@3Q; ~l, tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
n3$g x,KL \, R!S /R# ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
[/cIUQ T`\]!>eb arp=angle(rp);
o\4 CoeG zT7"VbP ars=angle(rs);
UW6VHA> >
H BJk: atp=angle(tp);
fagM7)x 2 z7}+lH ats=angle(ts);
k][ h9' PNpu*# Z` figure(1)
2 n)gpLIJ NEV p8) w subplot(1,2,1);
=(( yWn+t [jrqzB plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
]g,lRG xQU"A 2{}> legend('arg(r_p)','arg(r_s)')
MDHb'<o?y ^=+e?F`:{ xlabel('\theta_i')
uv]{1S{tb 8z}^j TM ylabel('\phi')
w"8V0z BvK QlT title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
8sg|MWS U c -B/~& axis([0 90 -3.5 3.5])
fu R2S70d A nMV < grid on
m#\[m<F g
j xS subplot(1,2,2);
"3}<8 c 3-|3`( plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
Ur!~<4GO ]Blf9h7 legend('arg(t_p)','arg(t_s)')
J,5+47b1}R fu4!t3 1 xlabel('\theta_i')
z%sy$^v@vD >c~RI7uu ylabel('\phi')
{@.Vh] "0Q1qZ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
U~H?4Izl= 66\jV6eH7L axis([0 90 -3.5 3.5])
V%NeZ1{ e H}Z Q?uK; grid on
}x :f%Z5h ~L{l+jK$p 
