计算光波由光密入射到光疏介质时相位变化
1X9s\JKQ .0~uM!3y clear
HPJ HA , mZjpPlJ close all
;K%/s IIke Z&P\}mm n1=1.45,n2=1;
,|.}6\zl*{ bX` Gv + theta=0:0.1:90;
='m%Iq7X -|>T?
t'K a=theta*pi/180;
#N'9
w . %O<8H7e)V rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
K>b4(^lf .`V$j.a rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
l BAu @M
a6 0rJ#GD tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Hrph>v N-Z=p)] ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
$ ]#WC\Hv 4RYH^9;>K arp=angle(rp);
`(=Kp=b BMu Ef a^ ars=angle(rs);
>xq. bG {<Gp5j atp=angle(tp);
*f`P7q* +oq<}CNr{ ats=angle(ts);
Im@OAR4,R PLMC<4$s figure(1)
n|X heG7: G -+!h4p subplot(1,2,1);
7 ]ysvSM )|^8`f plot(theta,arp,'-',theta,ars,'--','LineWidth',2)
K7K/P{@9[9 o[fg:/5)A legend('arg(r_p)','arg(r_s)')
0 yq w^$C\bCbh xlabel('\theta_i')
`[ U.BVP' isQOt *
i ylabel('\phi')
"42/P4: #jW=K&; title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
\=bKuP(it s&-MJ05y axis([0 90 -3.5 3.5])
q CYu@H o 0<NS1y grid on
MGt>:&s(] E G\;l9T subplot(1,2,2);
oW 6Hufu+o G8klWZAJ plot(theta,atp,'-',theta,ats,'--','LineWidth',2)
f17E2^(I(} UU MB"3e legend('arg(t_p)','arg(t_s)')
*Q/^ib9= n$r`s`} xlabel('\theta_i')
#?jsC) z+ {q Q! ylabel('\phi')
,_B n{T=U L\:m)g,F. title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
U i`{U tm^joK[{|J axis([0 90 -3.5 3.5])
vk 4 8&8 PQ[?zNr SV grid on
,) #rD9ZnC 7Y(Dg`8 G 
