利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
_FxeZ4\ m/3,;P.6 1、光疏射向光密
_:]g:F[
# > f'aW clear
c~+KrWbZ~ t:"3M iM=c close all
IG OEqUw* !/e*v>3u& n1=1,n2=1.45;
sC A ,KCxNdg^#- theta=0:0.1:90;
i5aY{ 3! Zs4NN 2 ~ a=theta*pi/180;
hHZ'*,9 y 5z&>NI rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
7,zE?KG / _2Mpzv rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
qOy0QZ#0 /0o#V-E) tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Sm{> 8e}UE _F5*\tQ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
/;k Sa}"Q ]!j%Ad figure(1)
Cu!]-c{ O7p=N8 V subplot(1,2,1);
/+1(,S 3W%j^nM plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
mV^~ BOW BD@y legend('r_p','r_s','|r_p|','|r_s|')
7pou(U
SsPZva xlabel('\theta_i')
*% Fu/ Sy' ]fGvx ylabel('Amplitude')
Nv^b yWqu ia (`3r title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
jtUqrJFlQ 4,>9N 9.?9 axis([0 90 -1 1])
-
AU {Y`j SF`(`h0e grid on
+~~2OU L }*C*!?pcd subplot(1,2,2);
)*`h)`\y \2]_NU5. plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
ITg<u?z_ 0 ?}n( f!S legend('t_p','t_s','|t_p|','|t_s|')
]T'7+5w a{@}vZx>3 xlabel('\theta_i')
T];dFv-GT )XHn.>]nc ylabel('Amplitude')
LM+d3|gSV P8Wv&5A title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
[Ky3WppR Rw|'LaW axis([0 90 0 1])
l*6Zh "o: Le:(;:eL>t grid on
h[&"KA YM};85 K Rp=abs(rp).^2;
* k<@ #=VYq4B = Rs=abs(rs).^2;
G[`1Yw$ J/O{x Rn=(Rp+Rs)/2;
dt^h9 I2O Zazff@O * Tp=1-Rp;
loO"[8i.k Bp3E)l Ts=1-Rs;
&!OEd ] DzQ Tn=(Tp+Ts)/2;
DY9]$h*y ]8}51 y8 figure(2)
?C#E_ N0.|Mb"?t subplot(1,2,1);
DU0/if9. Xl$, f`f~ plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
jj1\oyQ8 A4'5cR9T! legend('R_p','R_s','R_n')
jn:_2g[ /("7*W 2 xlabel('\theta_i')
K)BQ0v.:[ *8WB($T} ylabel('Amplitude')
,wwO0,"y7 ?oFd%|I title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
ATl?./T u Y}1c>5{bE axis([0 90 0 1])
xEp?|Q$ \ l#eW
x grid on
X!p`|i PO`p.("h subplot(1,2,2);
2lCFE) ~/]]H;;^u plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
o`,~#P| 0z8?6~M;< legend('T_p','T_s','T_n')
=9X1 +x lI 4tW= xlabel('\theta_i')
8HQ.MXKP 5$:
toL ylabel('Amplitude')
viG,z4Zf @Fb1D"! title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
%'yrI R ?B1Zfu0 axis([0 90 0 1])
iCE!Tm DT u3C_Xz grid on
l! F$V;R Lr:Qc#2 
y
%Q. ( ch8a 
:<t=??4m
