利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
n%Oq"`w4
15OzO.Ud
1、光疏射向光密 =sF4H_B
>=;hnLu
clear 9FPqd8(]*V
6C\WX(@4
close all ,aIkiT
w|ahb
n1=1,n2=1.45; `ZT/lB`
wN^^_
theta=0:0.1:90; I'[;E.KU
i)$ySlEh
a=theta*pi/180; HE>V\+
AL
_9q byhS7
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); 2/\I/QkTs
sE
^YOT<
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); 0.2stBw
xzOn[.Fi
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); AGKT* l.-
dI>cPqQ
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); 5K-,k^T}
xxwbX6^d
figure(1) GMB3`&qh
|*M07Hc x
subplot(1,2,1); fzOh3FO+
#Nad1C/]
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2) <$d2m6 J
_>;{+XRX[
legend('r_p','r_s','|r_p|','|r_s|') 'K01"`#
<PM.4B@
xlabel('\theta_i') z?F`)}
I
zVc
ylabel('Amplitude') +SW|/oIU
S!b?pl
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) ,_s.amL3O{
%~kE,^
axis([0 90 -1 1]) !6a;/ys
J8ScKMUN2
grid on k'+y
$&1D l
subplot(1,2,2); ;=k{[g 'gv
'St?nW3
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2) At
!:d3
.}T- R?
legend('t_p','t_s','|t_p|','|t_s|') W;os4'h$
?M6)O?[
xlabel('\theta_i') )1gT&sU