计算光波由光密入射到光疏介质时相位变化
@PuJre4!;L 3 D,PbAd clear
|$Dt6{h i`7{q~d= close all
6FG h=~{3, )hK5_]"lmj n1=1.45,n2=1;
A/RHb^ N bO1J#bcZ theta=0:0.1:90;
L
q8}z-? s `xp6\$ a=theta*pi/180;
1B= vr Gq 3;~1rw=$< rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
m8$ 6FN +o(t5O[G rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
W%b <(T;
0z/tceW'F tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Lx,"jA/ q,V JpqQ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
eS fT +UL AuUT 'E@E arp=angle(rp);
k:s}`h _n 9>u2;
'Ls ars=angle(rs);
* c7kB}/ +l@H[r;$ atp=angle(tp);
xB3;%Lc rZ 9bz}K ats=angle(ts);
dsX{ 5 [VI dw 92 figure(1)
meu\j g -"u}lCz> subplot(1,2,1);
^M6v;8EU (~xFd^W9o plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
VYTdK"% LW?] ~| legend('arg(r_p)','arg(r_s)')
W=} l=o!G. znhe]&Fw xlabel('\theta_i')
[U a4{3# u$[
'}z0: ylabel('\phi')
7oA$aJQ ?H eC+=/Z title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
>Mj :' Ks@S5:9sp axis([0 90 -3.5 3.5])
)`F ? {Sg r!V#@Md grid on
l@Ma{*s6=5 ,, mkB6; subplot(1,2,2);
'cy35M nf +8OH7 plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
su j? e6 3ag*dBbs legend('arg(t_p)','arg(t_s)')
\6)l(b; Sd'
uXX@ xlabel('\theta_i')
T
n Ad! \:4WbM:B ylabel('\phi')
v3Tr6[9 gSEj/? title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
+N"A5U sn_]7d+ Q axis([0 90 -3.5 3.5])
|@]J*Kh DF%\ 1C> grid on
8NudY3cU! -q&V V, 
