利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
n% Oq"`w4 15OzO.Ud 1、光疏射向光密
=sF4H_B >=;hnLu clear
9FPqd8(]*V 6C\WX(@4 close all
,a IkiT 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 byh S7 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);
fzO h3FO+ #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
zV c 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 0 ] SLeWs ylabel('Amplitude')
y u &mu CA `^
a:1^ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
# p] V? )>U"WZ'< axis([0 90 0 1])
{47Uu%XT ~wOMT grid on
qg|ark*1u S)W?W}*R\ Rp=abs(rp).^2;
h9!4\{V;h
+U%epq Rs=abs(rs).^2;
94|ZY }8|f d$xvM Rn=(Rp+Rs)/2;
Bjj =UtI k\9kOZW Tp=1-Rp;
[>\e @ = 4\k{E-x $ Ts=1-Rs;
aQf2}kD !%DE(E*'(
Tn=(Tp+Ts)/2;
<#199`R k sOc,4A figure(2)
?L H[,8z r!w*y3 subplot(1,2,1);
cLvnLaA} *5? Qam3 plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
5|>ms)[RQ uEG4^ legend('R_p','R_s','R_n')
QMmZvz\^ }_}
xlabel('\theta_i')
L/,g D.h^ A f!`7l- ylabel('Amplitude')
#<U@SMv [O|c3; title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
*uP;rUY fe"w--v axis([0 90 0 1])
Da!vGr 1zw,;m n grid on
0pl'*r*9 .j"heYF) subplot(1,2,2);
/u`Opv&I ( ]0F3@k#s plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
' V*}d w5rtYT I legend('T_p','T_s','T_n')
Ey%[t lbw+!{Ch xlabel('\theta_i')
u$aN~6 HG gB+CM?
LKq ylabel('Amplitude')
i_+e&Bjd4j Z= ;=9<vA title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
qW|h"9sr 2']0c
z axis([0 90 0 1])
c,ct=m.|6A
fCnwDT grid on
:_2:Fh.}3~ HU|qeSyel 
zQt)>Qx_ H!,V7R 
1=E} X5
