计算光波由光密入射到光疏介质时相位变化
u?H@C)P i`R (7Z clear
9lKRL'QR 5
>'66gZ close all
)U]:9) }iMXXXBOT n1=1.45,n2=1;
?THa5%8f O/(3 87= U theta=0:0.1:90;
LNaeB(z" dV)Y,Yx0${ a=theta*pi/180;
z}iSq$ (X*'y*: rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
:MK:T JV hN[X 1* rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
g
,`F<CF9 6={IMkmA tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
aXX,Zu^ y?"$(%3| ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
J_|7$
l/ F|6
nwvgq arp=angle(rp);
J`4Z<b53 Kq$:\B)<c ars=angle(rs);
cT\O v
P*_ MST:.x ; atp=angle(tp);
:D DO= qI(W$ ats=angle(ts);
oN_ S}o
=sxkr ih figure(1)
7F"3 <U@J "v@);\-V subplot(1,2,1);
TFNUv<>X xT:qe plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
4tbw*H5!5 C(f$!~M4b legend('arg(r_p)','arg(r_s)')
]U_5\$ T:be 9 5!, xlabel('\theta_i')
3Wjq >\ TViBCed40 ylabel('\phi')
4s[`yV
"(Mvl1^BT title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
o^8*aH)I>Y Jw2B&)k/ axis([0 90 -3.5 3.5])
k&WUv 0 :8](&B68gE grid on
?$UH9T9) =s`XZkh subplot(1,2,2);
&,^mM'
C c$@`P plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
wYZy e^7 p^T&jE8])# legend('arg(t_p)','arg(t_s)')
/7c2OI=\ >_rzT9gX& xlabel('\theta_i')
sGXp}{E9 fx]\) 0n ylabel('\phi')
-0{T 7h(HG ?2Y title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
x*NqA( r KK@
&q axis([0 90 -3.5 3.5])
> .a+: M' &J _g grid on
,O$Z,J4VL _cW_u?0X: 
