计算光波由光密入射到光疏介质时相位变化
84tuN YXXU Yi~!f clear
AWC zu5ve Dr6"~5~9w close all
K
!8+~[ bT- (lIU n1=1.45,n2=1;
"lLt=s2>L
A t%g^ theta=0:0.1:90;
'dM &~L SQ ^ cN- a=theta*pi/180;
8
W8ahG} A=|LMJMWR rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
q:#,b0|bv ;5# P? rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
*{tn/ ro6a N%6jZmKip tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
!Jb?r SJ.h ?
Ldw\ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
6}R^L (^M >4 GhI65 arp=angle(rp);
3?Y%|ZVM kM]:~b2 ars=angle(rs);
|tz1'YOB ',8]vWsl atp=angle(tp);
Gh0H)
q /MTS>[E ats=angle(ts);
Q@%VJPLv. lT$Vv= M figure(1)
Nt67Ye3; k<=.1cFh subplot(1,2,1);
xuO5|{ h {.SN plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
hU5[k/ q hF+YZU]rT legend('arg(r_p)','arg(r_s)')
tc@v9`^_ jD0^,aiG xlabel('\theta_i')
\A:m<:: S +|aCRS ylabel('\phi')
Yg/ e 8Q2 O(,Ezy x title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
& bh?jW 9cFFQM|o axis([0 90 -3.5 3.5])
~^"
cNv `j[)io k grid on
: <gmgI Rg<y8~|'} subplot(1,2,2);
^SxB b,\ 7u,56V?X plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
(rt DT 82Nh;5T r legend('arg(t_p)','arg(t_s)')
JoKD6Q1D Z4}Yw{=f xlabel('\theta_i')
B9iH+
]W ~ 6!=_" ylabel('\phi')
m^U\l 9LE ]N-K`c] title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Y {Klwn Gw! jYnU axis([0 90 -3.5 3.5])
R[wy{4<y Qz{:m grid on
Y1{6lhxgE &AGV0{NMh] 
