计算光波由光密入射到光疏介质时相位变化
~ Frk LP :g :h 0'G clear
87 E3pe NCpn^m)Q} close all
rZCAj C-b% PgA n1=1.45,n2=1;
48NXj\L[y ua>~$`@gX theta=0:0.1:90;
N~<}\0 ?)QBJ9F a=theta*pi/180;
b 0x0CMf W%Nu]9T rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
6rC P]YnF *SG2k .$ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
b2-|e_x v2X0Px_ tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
8!`.%)- 4 eeTaF!W ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
fB"gM2' l7 y`$8Co arp=angle(rp);
Y52xrI vl\ 7tO$'q*h ars=angle(rs);
~W2&z]xD E/-Kd!|" atp=angle(tp);
6uE1&-:L Uht:wEr ats=angle(ts);
"X}F%:HL QgM_SY|Rj figure(1)
w*kFtNBfU ={vtfgxl subplot(1,2,1);
72.IhBNtT )KQv4\0y< plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
>w# 3fTJ dnc!=Z89 legend('arg(r_p)','arg(r_s)')
_ll aH 2s ,n!u
Fd xlabel('\theta_i')
k3h53QTmC b'"% ylabel('\phi')
0c6AQP"=V YH tI% title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
QjbPBk Q ##ea-"m8 axis([0 90 -3.5 3.5])
/4BXF4ksi, ghR]$SG grid on
m:<3d]L b<cM[GaV~ subplot(1,2,2);
)L("t t[oc p;Q plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
w:@M|O4` fT~<C
{ legend('arg(t_p)','arg(t_s)')
~?+m=\ }H^^v[4 xlabel('\theta_i')
BV:,b S 5i&V ~G ylabel('\phi')
KA2B3 \ ?kefRev<#h title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
v@SrEmg 5E(P,!-. axis([0 90 -3.5 3.5])
Gnq~1p5^ ?q&mI* j! grid on
/Loe y
. M $D 
