计算光波由光密入射到光疏介质时相位变化
6nxX~k IooAXwOF clear
.+ 1.??8:+ V!!E)I close all
*ea%KE": 21Z}Zj n1=1.45,n2=1;
nic7RN ?F< CXqU< a& theta=0:0.1:90;
R~40,$e{ fIOI a=theta*pi/180;
u&)+~X b!W!Vvf^x rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
|Sg
FHuA sQt]Y&_/@ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
[t^Z2a { 1+;Z0$edxz tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
KiXXlaOs
NYwE=b~I ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
%E q} H o!TG8aeb arp=angle(rp);
NABwtx >. '1?b ?nVo ars=angle(rs);
Q v/}WnBk G(7!3a+ atp=angle(tp);
zyNg?_SM _tE`W96J ats=angle(ts);
#7S[Ch}O UEeD Nl$^u figure(1)
O ][R "5d s;-(dQ{O subplot(1,2,1);
qaK9 E@l 2/.E uf plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
}K\ ] M@ _F|oL| legend('arg(r_p)','arg(r_s)')
jOs&E^">&B h
WvQh xlabel('\theta_i')
Obd@#uab # biI=S ylabel('\phi')
c]]OV7;)>
hS) X`M title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Z1j3 F 9pN},F91n: axis([0 90 -3.5 3.5])
]qZs^kQ uW*)B_c grid on
BM+v,hGY O)g\/uRy subplot(1,2,2);
.Y}~2n m
Cvgs plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
-nP
y?>p"| }^).Y7{g[ legend('arg(t_p)','arg(t_s)')
M*|,05> G]Fp}, xlabel('\theta_i')
O gp Zwwk Dh=?Hzw ylabel('\phi')
eh=.Q<N s!* m^zx title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
ay}} v7)GM 'l, ym~R axis([0 90 -3.5 3.5])
msKWb311u $@ t ]0 grid on
}Gmwm|`* nM*-Dy3ou 
