利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
2u>
[[U1: ,%L>TD'48s 1、光疏射向光密
neHozmm| $_5@ NOZ,M clear
, IODV`L +yh-HYo` close all
Vd' KN2Jm UT^-!L
LB] n1=1,n2=1.45;
|.s#m^" K?4/x4p@ theta=0:0.1:90;
Dn}Wsd= )R QX1("O a=theta*pi/180;
N~w4|q!] gm-m_cB< rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
[q MFLY$ 2K 8?S rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
_QhB0/C @k ~_ w# tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
GmA5E L POZA` ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
} - e HfF4BQxm figure(1)
kRyt|ryWh y[ }O( subplot(1,2,1);
Ix"hl0Kh LGOeBEAMV^ plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
.Fz5K&E= 4/Vy@h"A3 legend('r_p','r_s','|r_p|','|r_s|')
w84
] s%y A ko}v"d xlabel('\theta_i')
T@GR Tg 'pUJlPGx ylabel('Amplitude')
WF6'mg^^ ? G isI/ Ir[ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
{GaQV -t S[TJ{ L( axis([0 90 -1 1])
+18)e;
QSn%~o05 grid on
<z>K{:+> xEL nik_L2 subplot(1,2,2);
`q ;79t %Jrdr`< plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
yrDWIU(8;6 $ljgFmR_ legend('t_p','t_s','|t_p|','|t_s|')
U#B,Q6~ DIk\=[{2q xlabel('\theta_i')
-zeodv7 Z66b>.<8 ylabel('Amplitude')
:Rs% (Z xLE+"6;W title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
V/0?0VKG 0I.9m[<Fc axis([0 90 0 1])
a{lDHk `Wf bb6
~H grid on
sP0pw]! xHml" Y1 Rp=abs(rp).^2;
(1bz.N8z N.J;/!%! Rs=abs(rs).^2;
@17hB h AUloP?24 Rn=(Rp+Rs)/2;
CS 7"mE`{ q\q V~G` Tp=1-Rp;
eQk ~YA] K ub "(,k P Ts=1-Rs;
zS Yh ?NB5 ZCQ7xQD Tn=(Tp+Ts)/2;
7'[C+/: HQ%-e5 Q figure(2)
$*| :A (D'Z 4Y subplot(1,2,1);
TQ? D*& )Oq N\ plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
4#5w^ i<g|+}I legend('R_p','R_s','R_n')
`_]Z#X&&h WUid5e2 xlabel('\theta_i')
U*Z P>V v p[(VhbN ylabel('Amplitude')
mMq T-jT \TG!M]D: title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
%Fc, $ = I/bED~Z:a axis([0 90 0 1])
xMsos?5} ;Ef:mr"Nu grid on
r}i< cyL %/dYSC
subplot(1,2,2);
}>JFO :v& D4yJ:ATO& plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
$]CZ]EWts DBl.bgf legend('T_p','T_s','T_n')
Xj^6ZJc L<> ;E xlabel('\theta_i')
'sJ=h0d_[V 1<]g7W ylabel('Amplitude')
Wh%q v V6] y
D.S" title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
v{ .-x\; 09"C &X~ axis([0 90 0 1])
R@``MC0 V *=T o grid on
uROt h_/ oC>QJ(o,8 
>Dpz0v E|D~:M%~ 
GzK{. xf
