计算光波由光密入射到光疏介质时相位变化
jBZlN Ew 2;Vss<hR4A clear
5.m&93P x93h{K f close all
'NJGez'b , w0qrh\3du n1=1.45,n2=1;
E Q
'L" Y-!~x 0-H theta=0:0.1:90;
7k,pUC-w7c aF*KY<w a=theta*pi/180;
\ #<.&`8B <;Q1u,Mc rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
W>f q 9 !d nCr R rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
er@"4R0 tfB}U. tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
X$*MxMNs
&
-r^Q ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
gJa48 pi %+ln_lgD: arp=angle(rp);
w`BY>Xft0 SeuC7!q{ ars=angle(rs);
xgDd5 `W +85#`{ D atp=angle(tp);
jfmHc(fX4 p7{2/m j ats=angle(ts);
o\/&05rp] h/oRWl0r figure(1)
&10l80vj _{j'` # subplot(1,2,1);
OBb m?`[ -e_91W I plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
gB7kb$J ))T@U?r legend('arg(r_p)','arg(r_s)')
5(wmy-x\ x[?N[>uw xlabel('\theta_i')
0gRj3al( ;x<5F+b ylabel('\phi')
rX7GVg@H *y+N-uq title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
TxJoN]Z. oW}nr<G{< axis([0 90 -3.5 3.5])
m}UcF oaO 6jw 9p+. grid on
|>^5G@e &9Y ^/ W subplot(1,2,2);
6a*?m{ Pj-.oS2dA plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
m+UWvUB) Yg /g9 $' legend('arg(t_p)','arg(t_s)')
45.<eWH$*( e{Q;,jsh xlabel('\theta_i')
2LfiaHO DC samOA~ ylabel('\phi')
Z8Iqgz7|y !hs33@*u~ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
ag V z
&Fiesi!tET axis([0 90 -3.5 3.5])
V`
T l$EF -F. c<@*E grid on
H(hE;|q/ 0g 2? 
