利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
G0&'B6I> &L`yX/N2 1、光疏射向光密
'e5,%"5(c v7@O ,% clear
Sxg&73;ZV %y _AT2A close all
+ 3%i7 $ Ov#^wfA n1=1,n2=1.45;
->Bx>Y ruK, Z,3Q theta=0:0.1:90;
VaLl$w 3P[u>xE a=theta*pi/180;
+U= !svE > J>|+W rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
z1`z
k 0 F,>-+~L= rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
]n$&|@ pSh$#]mZ` tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
l5fF.A7TT F}dq~QCzw ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
r,-9 ]?i vB;$AFh{ figure(1)
rN5;W 6> Ca O subplot(1,2,1);
9o|#R& 0 ? =a, plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
+KKx\m* 9a;8^?Ld%S legend('r_p','r_s','|r_p|','|r_s|')
dpI9DzA; ~s UWXw7~ xlabel('\theta_i')
;_<K>r* ZFO*D79:K ylabel('Amplitude')
kFsq23Ne X"p p l 7o title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
"}'Sk( W99Hq1W;r axis([0 90 -1 1])
6A \Z221E lM'yj}:~ grid on
cAAyyc"yJ y.m;4(( subplot(1,2,2);
h<ULp &g |nTZ/MXbw plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
Q1(6U6L %*`yd.L0W legend('t_p','t_s','|t_p|','|t_s|')
GtNGrJU Q($aN- xlabel('\theta_i')
vhaUV#V" XOi[[G} ylabel('Amplitude')
<{yQNXf[ y$^.HI02jP title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
RS`~i8e' Y%iimbBY| axis([0 90 0 1])
SuU %x2 C;_*vi2u grid on
ezR !ngt d[Lr`=L; Rp=abs(rp).^2;
WCPl}7> A}cGag+sp Rs=abs(rs).^2;
Rf .b_Y@O L4,Ke Rn=(Rp+Rs)/2;
J)xc mK 3x{2Dh i Tp=1-Rp;
b!ea(D!: P Zc{wbjp& Ts=1-Rs;
F&I ;E i &QQ8 ut, ; Tn=(Tp+Ts)/2;
(`& `vf Oor &1 figure(2)
}|7y.* (~<9\ZJs subplot(1,2,1);
ugI9rxT]Kv 30Z RKrW"~ plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
&@MiR 8 3h|:ew[ legend('R_p','R_s','R_n')
$SmmrM +Eh^j3W xlabel('\theta_i')
7~V,=WEe *4tJ|m6"Y6 ylabel('Amplitude')
)KLsa`RV: `=h`:` title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
<
<vE . =(~UK9` axis([0 90 0 1])
uM^eoh_ -b4#/q+bb+ grid on
ZAG ia q dw
v(8 subplot(1,2,2);
~KufSt * 2AAZZx +$ plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
?T(>! m T&] J3T FJ legend('T_p','T_s','T_n')
=a<};X hF!t{ Lf3 xlabel('\theta_i')
TFWx(}1 #i Ooi9( ylabel('Amplitude')
?1?zma S h o7L@NR title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
>*PZ&" }M x}W,B,q axis([0 90 0 1])
xnW3,:0 )-_]y|/D:r grid on
G' ~Z' ]mN'Qoc 
[`~E)B1Y +%'S>g0W= 
q}PeXXH
