利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
w/L ` #h@J=Ki 1、光疏射向光密
6_4D 9 W BAO| )~1Pd clear
c_"
~ n| X*L;.@xA close all
F!.E5<&7= ^#0U ?9 n1=1,n2=1.45;
HQ{JwW!m pZU9^Z?~6 theta=0:0.1:90;
~l2aNVv; UswZG^Wh a=theta*pi/180;
/ >j';6vi Q%-di= rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
{%#)5l) 'L{8@gq i rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Rml2"9"` !'Q -yoHKD tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
4Yl ; cr;g5C
V ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
b^[F ""!e oc^Br~ Th figure(1)
dZ _zg< 1 n%?l[o subplot(1,2,1);
*kIc9 } !yhh8p3 plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
\R VW +<prgP`v legend('r_p','r_s','|r_p|','|r_s|')
1xc~` ~ &x > B xlabel('\theta_i')
s#uJ
;G _ {|D ylabel('Amplitude')
sm m]6 _f6HAGDN title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
b$eXFi/ 4H+Ked&Oq axis([0 90 -1 1])
*SO{ \bu mR;qMX)0h grid on
Ssk}e=] R=T qj,6 subplot(1,2,2);
j:w{;(1=W qp}Ma8 + plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
7r{ 83_B CB&iI' legend('t_p','t_s','|t_p|','|t_s|')
srV.)Ur nhN);R~o"1 xlabel('\theta_i')
-rKO
)} Gy[O)PEEh ylabel('Amplitude')
we#w H- (Y^X0yA/ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
I L_d:HF|1 jLJ1u/l>; axis([0 90 0 1])
3 v")J*t c/ 5W4_J grid on
G HQ;hN: u0`% +:]0 Rp=abs(rp).^2;
hva2o` S!{t6'8K Rs=abs(rs).^2;
_sy'.Fo oy<W Ub9W Rn=(Rp+Rs)/2;
KFZm`,+69 BA=,7 y&;j Tp=1-Rp;
6'W [{gzl ,&4qgp{) Ts=1-Rs;
<[GYLN[0Q S60IPya Tn=(Tp+Ts)/2;
'$ef+ @y ~}ET?Q7t figure(2)
>*$; >h+349 subplot(1,2,1);
f+.T^es OMk5{-8B plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
`m.eM ;#MB7A
legend('R_p','R_s','R_n')
7a=S i*eAdIi xlabel('\theta_i')
*6BThvg|&X 1oKfy>i e ylabel('Amplitude')
0hZ1rqq8C IcI OC8WC title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
5(
_6+'0 C>Is1i^9 axis([0 90 0 1])
",>H(wJ8 ~m *,mz grid on
-|Kzo_"
v5 u
=gt<1U subplot(1,2,2);
=e?$ M TEsnN i
1 plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
dC }`IR !AJ]j|@VBd legend('T_p','T_s','T_n')
,YhwpkL Iak0 [6Ey xlabel('\theta_i')
gK|R =J 2l5KJlfj>k ylabel('Amplitude')
]S 7^ITn k
n8N,,+
title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
I?Q+9Rmm`J *^ \FI Ud axis([0 90 0 1])
5. tvB <Q<+4Y{R grid on
$\AEWFB kcB+ _ 
djT5 X Eaad,VBtU 
bl/,*Wx:4.
