计算光波由光密入射到光疏介质时相位变化
@_ tA"E \)6?u_(u clear
*b 7
^s,? ?eg@
7n close all
#ui%=ja[:~ {t IoC;Y n1=1.45,n2=1;
B#/~U`t* y^Xxa' y theta=0:0.1:90;
x:D<Mu# d' eM(4R@ a=theta*pi/180;
*dn-,Q%` ) F9%^a( rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
V 1+o3g{} W} +6L| rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
I$P7%} eC1c`@C: tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
d T -O8 ?[|4 QzR ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
CEJG =*3 5'}!v arp=angle(rp);
tGy%n[ \ u/{_0-+P ars=angle(rs);
]H@uuPT! S
g_?.XZc[ atp=angle(tp);
]r{ #268 f>JzG,- ats=angle(ts);
w})&[d xN~<<PIZ figure(1)
[}g5Z=l 0eu$ oel- subplot(1,2,1);
&T 4Cn@ kO\&mL&
qD plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
K
{N;k- L SOwa legend('arg(r_p)','arg(r_s)')
jC }u>AB WlJRK M2 xlabel('\theta_i')
Y&Nv>o_}5 mwiPvwHrg ylabel('\phi')
0~I)
/T hCx#H eh title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
IaZAP jI pcMN< axis([0 90 -3.5 3.5])
R5YtCw]i= ' k) P(H grid on
D2mAyU - 53#5p;k
subplot(1,2,2);
X=7vUb,\gB Wh&Z *J plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
="*C&wB^ e @MCumc~+ legend('arg(t_p)','arg(t_s)')
zCN;LpbEJY S6Xw+W02 xlabel('\theta_i')
S %%q n W;j)ux7jMY ylabel('\phi')
bJu,R-f A}+r;Y8[h title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
T%b^|="@ lIT2 AFX+ axis([0 90 -3.5 3.5])
s`#j8>`M
WxI]Fcb< grid on
/]UN N~( vTa23YDW 
