计算光波由光密入射到光疏介质时相位变化
mXJ`t5v^l ^"GDaMF clear
fb_q2p}
G N@3&e;y close all
::k
cV'* Rwe!xY^d8 n1=1.45,n2=1;
`6FH@" |I {Z_?7J&z theta=0:0.1:90;
.gs:.X)TG9 4#&w -W a=theta*pi/180;
Z11 I1)%s /"
6Gh' rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
HzZ.q2Zz% xW+ XN`77 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Z,:}H6Mj9 ot;
]?M tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Xp?W oC N &.chqP(| ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
VmP5`):?b ^wW{7Uq> arp=angle(rp);
=(P k7{ r [4dGt ars=angle(rs);
Y=(%t:#_
XM< atp=angle(tp);
h{.x:pPXy M2piJ'T4u ats=angle(ts);
<ql:n ZL+46fj figure(1)
3fq'<5 ^ M <c cfU ! subplot(1,2,1);
4R28S]Gb QB6.
o6 plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
4 mwLlYZ K?gO ]T{6 legend('arg(r_p)','arg(r_s)')
[9 ,34/i l;vA"b=] xlabel('\theta_i')
.v[!_bk8C F0tx.]uS ylabel('\phi')
o>MB8 [r _GsHT\ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
8c' 5P =fcRH:B: axis([0 90 -3.5 3.5])
^#&P Tq> ~ 't+X grid on
17S<6j#H5 +W#["%kw subplot(1,2,2);
g] m}@b6(h
S)W(@R+@4 plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
#p/'5lA&j ]Ag{#GJ5D legend('arg(t_p)','arg(t_s)')
:Nl.< 6+ {IT;g9x xlabel('\theta_i')
O+PRP"$g" jGFDj"Y ylabel('\phi')
V
u1|5 4S*7*ak{ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
h45RwQ5Z 8rM1kOCf axis([0 90 -3.5 3.5])
'OvyQ/T
v0,&wdi grid on
Q vm[2mb L ]QBh\ 
