利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度 lh w]?\
l+Dl~o }
1、光疏射向光密 dB +x,+%u+
%_0,z`f
clear
O<Qa1Ow7f
v_)a=I%o&2
close all JZ Qkr
S(9Xbw)T
n1=1,n2=1.45; R $HI JM
"D}PbT[V
theta=0:0.1:90; f u?5gzT+b
b Gq 0k&
a=theta*pi/180; )f-u x5
%Fig`qX
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); *vc=>AEc
F`S OF O
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); 79U
Th@r}
AFcsbw
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); iDt^4=`
@qmONQ eb
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); P*oKcq1R
s
`HSTq2
figure(1) F(>']D9$.
CB<1 ] Z
subplot(1,2,1); &sR=N60n
0@d )DLM?
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2) #qXE[%
Xt~`EN
legend('r_p','r_s','|r_p|','|r_s|') zvf:*N a")
@P#uH5U
xlabel('\theta_i') qIcQPJn!}
P1;T- .X~&
ylabel('Amplitude') -FytkM^]6
$C0Nv Jf
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) , C2qP3yg
mt3j- Mw
axis([0 90 -1 1]) b/Y9fQ n
I'2:>44>I6
grid on @/*{8UBP
6NH.!}"G9
subplot(1,2,2); lS]<~
<8Ek-aNNt
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2) WLW' .
/AV
[g^x2
legend('t_p','t_s','|t_p|','|t_s|') x7K
s:lar4>kM
xlabel('\theta_i') %^[45 e
(__=*e w
ylabel('Amplitude') 3dfG_a61y
$T{,3;kt
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) *cx mQ
3":ef|w]
axis([0 90 0 1]) {Md xIp[
# eqt{
grid on rs {e6
`fh^[Q|4n0
Rp=abs(rp).^2; *vv <@+ gA
6eE%x?#
Rs=abs(rs).^2; $ r|R`n =
Vd3'dq8/?
Rn=(Rp+Rs)/2; p
2x OjS1
B3K!>lz
Tp=1-Rp; #q`[(` Bx
0}Rxe
Ts=1-Rs; n={} ='
VTk6.5!8
Tn=(Tp+Ts)/2; H+vONg
+jYO?uaT
figure(2) E@TX>M-&
4O_z|K_ k|
subplot(1,2,1); _F>1b16:/P
vF"<r,pg
plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2) `?LQd2p
7IW:,=Zk8+
legend('R_p','R_s','R_n') JPfNf3<@My
B04%4N.g"X
xlabel('\theta_i') J9S9r ir &
QEL^0c8 ~
ylabel('Amplitude') ! utgo/n
I2kqA5>)j
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) m 7 6**X
O$ u;]cg
axis([0 90 0 1]) 1rh\X[@
0/hX3h
grid on %y.9S=,v,
/tRzb8`
subplot(1,2,2); _?>!Bz
m
mN +~fu h
plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2) 9+@"DuYc6
#S2LQ5U
legend('T_p','T_s','T_n') kwNXKn/
h
Vz%{R"
xlabel('\theta_i') =o;QvOS;
X<@yt HBv
ylabel('Amplitude') ,Hh7' `
nL!h hseH
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) nR4L4tdS
XSt5s06TM
axis([0 90 0 1]) >SpXB:wx
dVc;Tt
grid on zOYkkQE3mJ
PbxuD*LQ.
.O @bX)
zMs]9o
1<A+.W
