计算光波由光密入射到光疏介质时相位变化
|n8^Xsx4w l:HuG! clear
xH uyfQLk rD g l@B3 close all
3E-&8x7uYR 'z/hj>B< n1=1.45,n2=1;
/FY_LM 8wOPp dc theta=0:0.1:90;
^gImb`<6- VXp
X#O a=theta*pi/180;
u#<]>EtbB ]LjW,b" rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
fsoS!6h0k qS>el3G rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
=UN: IzT '1<QK tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
; V8 =B8w gC.T5,tn ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
4j5plm= ="4 )! arp=angle(rp);
N T0q!r/! "-
AiC6u ars=angle(rs);
7JbrIdDl| t [ f] atp=angle(tp);
&I8ZVtg ~Q5HM ats=angle(ts);
QMP:} %?ad.F+7 figure(1)
p6p_B ! WN r09` subplot(1,2,1);
E@p9vf-> u56cT/J 1 plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
G6FknY j |yr }g-m legend('arg(r_p)','arg(r_s)')
>K3Lww)Ln =x> KA*O1 xlabel('\theta_i')
kq+L63fZ xQ4Q '9 ylabel('\phi')
6Y=)12T Q =Z-vTD+ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
3$_wAt4w 6;Bqu5_Cj axis([0 90 -3.5 3.5])
mSs%g L]g gwyz)CUkL grid on
9#+X?|p+0 h8Gp>b subplot(1,2,2);
N]3XDd|q ^V D14V 3 plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
;TYkJH" 8WMC ~ legend('arg(t_p)','arg(t_s)')
s&4Y+dk93 yfj<P/aA+ xlabel('\theta_i')
cz afBO6 3LG)s :p$/ ylabel('\phi')
qbjRw!2?w 9kcAMk1K title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
&W1c#]q@r = Z 2sQQVS axis([0 90 -3.5 3.5])
u?!p[y6 Gmc0yRN grid on
z'
@F@k6 =73wngw 
