计算光波由光密入射到光疏介质时相位变化
S(_DR 8 ~Sx\>wBlc clear
&a'LOq+r' c9+yU~( close all
*C.Kdf3w HP:[aR!2P n1=1.45,n2=1;
rGay~\ #j"GS/y" theta=0:0.1:90;
Nwk^r75l q \og2\Oh&gH a=theta*pi/180;
593!;2/@ 0+A MN- rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
*TPWLR ^ T[2<_ nn= rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
FhQb9\g A| y U'k tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
m,ur{B8 : 'a>D+A: ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
=)J<R ; uYijzHQyD arp=angle(rp);
~w}=Oby'y wcGv#J], ars=angle(rs);
uFFC. w GZm=>!T atp=angle(tp);
Jv!f6*&< cvYKZ B ats=angle(ts);
}$s QmR R !wKiMgLS figure(1)
nIqF:6/ @ @=,bO subplot(1,2,1);
(
geV(zT 'h[7AZ&)# plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
ZHCrKp CiTWjE?|7 legend('arg(r_p)','arg(r_s)')
(N5"'`NZA !p_l(@f xlabel('\theta_i')
^!x qOp! |7!B k$(vA ylabel('\phi')
"tu*(>'~5 5[~ C!t; title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
!$<Kp6 ::@ JL axis([0 90 -3.5 3.5])
#z}0]GJKj !e('T@^u6u grid on
/9| 2uw` AG"iS<u subplot(1,2,2);
-?Cu-' R S] N%`] plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
kRH
D{6mol qipS`:TER legend('arg(t_p)','arg(t_s)')
!. :b}t 3M>y.M S xlabel('\theta_i')
7<k@{xI/ -o57 "r^x ylabel('\phi')
(A -Uo
1 K] title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
;|f]e /El mu O; g& axis([0 90 -3.5 3.5])
K] &GSro ,? Q1JZPy@ grid on
{f Mo#`9= |WW'qg]Uu 
