利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度 $'CS/U`E}
}f;TG:6
1、光疏射向光密 DW,ERQ^
lGahwn:
clear =4+2y '
zfDfy!\2_
close all u `w w
i;~. kgtq4
n1=1,n2=1.45; O D5qPovsd
T0fm6
J
theta=0:0.1:90; 5xKod0bA
^vh!1"T
a=theta*pi/180; n\ Ls m
:s+?"'DP
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); Zt41f PQ
,^
,R .T
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); T*B`8P
V G7#C@>Z
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); &b:y#gvJ:
rgXX,+cO
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); I3.JAoB>!
4#W$5_Ny
figure(1) lrkgsv6
/AX)n:,
subplot(1,2,1); "MzBy)4Q
g'k m*EV
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2) =_E$* }
0 7qjW o/t
legend('r_p','r_s','|r_p|','|r_s|') 9`qw,X&AK_
BJHWx,v
xlabel('\theta_i') GZ5 DI+3
) X *_oH=
ylabel('Amplitude') (oCpQDab@
Yt_tAm
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) !j #8zN
MsIaMW _
axis([0 90 -1 1]) k=d _{2 ~
fZap\
grid on &< &eK q
zGd[sjL
subplot(1,2,2); GRj [2I7:
DV?c%z `YO
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2) lM#/F\
;O>zA]Z8r
legend('t_p','t_s','|t_p|','|t_s|') YJwI@E (l$
9^sz,au B
xlabel('\theta_i') eGKvzu
2sqH
> fen
ylabel('Amplitude') M?sTz@tqq
\
D>!&
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) |'}r-}
,Cm1~ExJ
axis([0 90 0 1]) X6 !KFc
B|^=2 >8s
grid on C@ XnV=J
aY, '^S
Rp=abs(rp).^2; UY *Z`$
;x|? N*
Rs=abs(rs).^2; Hjs#p{t[
o<s~455m/
Rn=(Rp+Rs)/2; n[a%*i6x
Xa'b @*o&
Tp=1-Rp; um". Z4S
^=-W8aVi>
Ts=1-Rs; L Do~
g_Y$5ft`
Tn=(Tp+Ts)/2; oO
&%&;[/A
'./ qB J
figure(2) Z_jV0[\v0P
9,>c;7s X
subplot(1,2,1); }*P;kV
s<C66z
plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2) v~dUH0P<>e
ey>tUmt6?
legend('R_p','R_s','R_n') Cw.DLg
HZ{DlH;&
xlabel('\theta_i') "Q.C1#W}.
+q432ZG
ylabel('Amplitude') iqd7
,0, &
L
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) /T&+vzCF
)q{e L $
axis([0 90 0 1]) \rbvlO?}
)Y1 +F,C
grid on cR6 #$-a
bHs},i6
subplot(1,2,2); 'A/{7*,
m0Uk*~Gz
plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2) XPi 5E"
?,Wm|xY
legend('T_p','T_s','T_n') riglEA[^
I' R| B \
xlabel('\theta_i') srU*1jD)
:7qJ[k{g
ylabel('Amplitude') hZO=$Mm4p
mL s>RR#b
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) F
B&l|#e
"?
5@j/
e`
axis([0 90 0 1]) }Mc&yjhMrg
{`% q 0Nr
grid on Y-"7R>^I
3)dtl!VMW[
+227SPLd
&nn+X% m9g
[k,FJ5X
