利用菲涅爾公式計(jì)算光波在兩種介質(zhì)表面折反射率及折反射能流密度
j:k}6]p} w2 %u;D% 1、光疏射向光密
iB-h3/ NWL\"xp
`t clear
B8=r^!jEL ={'*C7K)oK close all
^ &UezDTS 8&K1;l } n1=1,n2=1.45;
VT@,RlB0 8nV#\J9 theta=0:0.1:90;
y&A0}>a:d [v0[,K a=theta*pi/180;
u
V[:e|v /i!3Fr" rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Wb)>APL [l`_2{: rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
@$:T]N3m 6(M^`&fl tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
8VWkUsOoI J~jxmh ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
&Hl*Eg
f 4k7
LM] figure(1)
E8gbm&x* fC4#b?Q subplot(1,2,1);
JyiP3whW LA +BH_t& plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
pYxdE|2j s8O.yL legend('r_p','r_s','|r_p|','|r_s|')
9>QGsf.3 PQ0l