利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
{vox
x&UX ^
RIWW0 1、光疏射向光密
?cKe~Q?3 ,ePl>m:Z
clear
P&5kO;ia EPd
close all
^2\-zX!bt 3e ?J#; n1=1,n2=1.45;
-s&7zqW 0T.kwZ8 theta=0:0.1:90;
W,bu=2K6 TxvvCV^
a=theta*pi/180;
@r3,|tkrz .kp3<. rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
|oX1J<LM bu,xIT ^ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
M@<r8M]G DsCbMs=Y tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
>TG# 7_CX6: ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
p@[n(?duC. fir#5,*q| figure(1)
I/)*pzt8 )(-aw,iK subplot(1,2,1);
, )PpE& $ 9
k5a plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
^a?g~G 6!Uk c'r legend('r_p','r_s','|r_p|','|r_s|')
Q9{f'B |u_fVQj xlabel('\theta_i')
EnP> f[Xsri ylabel('Amplitude')
0S&C[I
o6 3,)[Q?nKD title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
yM#
%UeZ\ 9\!=i axis([0 90 -1 1])
bA\(oD+: $%.,=~W7 grid on
.Z(Q7j^ &EJ/Rl subplot(1,2,2);
P#-p*4 !112u#V plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
9yWSlbPr] 4P}d/w?'KL legend('t_p','t_s','|t_p|','|t_s|')
b]g#mQ hQwUwfoe@ xlabel('\theta_i')
hp$1c .$a|&P=S ylabel('Amplitude')
9gac7(2`) l._g[qa title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
}gKY_e3 o ]@'R<F(u axis([0 90 0 1])
cv-;fd>' %|B$y;q^3 grid on
5V8`-yO9
+Ou<-EQV Rp=abs(rp).^2;
t@!A1Vr@ -q&7q Rs=abs(rs).^2;
&Xh=bM'/%m ~toR)=Yv Rn=(Rp+Rs)/2;
:
`,#z?Rk \0.!al0 Tp=1-Rp;
wowWq\euY &neB$m3y Ts=1-Rs;
!*PX- ]-jaIvM Tn=(Tp+Ts)/2;
Mo]aB:a [~ !9t9+~ figure(2)
"rHPcp"m c3(0BSv subplot(1,2,1);
\-D[C+1( =yZ6 $ hK plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
{EJ+
F`CDv5 legend('R_p','R_s','R_n')
B$MHn? 'j;i4ie>*x xlabel('\theta_i')
`2`h4[^ [X BO0Y#fs ylabel('Amplitude')
\jr-^n] .F3~eas title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
kH?PEA! \ HG^~7oMf axis([0 90 0 1])
wlpcuz@ .J?RaH{i grid on
7pM&))R 74a k|(! subplot(1,2,2);
7%W@Hr,%F irMBd8WG plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
j-v/;7s/B 2{S*$K[M legend('T_p','T_s','T_n')
5")BCA wsgT`M'J[ xlabel('\theta_i')
'y7<!uo? V o%GO9b; ylabel('Amplitude')
rB_ESNx 3935cxT1U title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
P_(<?0l n1)]. ` axis([0 90 0 1])
6>b'g
~I 9#a/at] grid on
~R7{gCqdr ,irc=0M( W=4|ahk$ :@S=0|:j !'-./LD")