利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
%/UV_@x& ~hS3*\^~M 1、光疏射向光密
B3 y?. = E&b= clear
F;IP3tD =LlLE<X"%x close all
CTl (_g qq
OxTG] n1=1,n2=1.45;
ooTc/QEYi `+roQX.p theta=0:0.1:90;
/3J z3 Yuwc$Qp) a=theta*pi/180;
O4m(Er@a @)o0GHNP rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
uzHT.iBn J6*f Uh rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
;c1relR2 Wt>J` tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
j'3j}G%\T JT4wb]kdV ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
HwB {8S?sm t:\l&R& figure(1)
r)gCTV(kb <v]z6B@9! subplot(1,2,1);
uCw>}3 z<a$q3!# plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
i*X{^A73" /r2 76 Q legend('r_p','r_s','|r_p|','|r_s|')
wz|DT3"Xs '9 'l=Sh xlabel('\theta_i')
;xfO16fNk 3I|O^ ylabel('Amplitude')
w'A tf - ^yXLa;D title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
<x0)7xX BYS lKTh axis([0 90 -1 1])
$ K})Q3FNi C`)_i3
^ grid on
O#vn)+Y,* 4AQ[igTDP subplot(1,2,2);
G%S6$@: >W
r$Y{ plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
b6|Z"{TI
_ ~F; CE"3A legend('t_p','t_s','|t_p|','|t_s|')
cQX:%Ix= :V-k'hm
& xlabel('\theta_i')
W@^J6sH .l \r9I( ylabel('Amplitude')
k=?^){[We t)} \9 ^Uo title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
#EJP(wXa 6a*83 G,k axis([0 90 0 1])
H8&p<= \QGa 4_# grid on
wZ jlHe #1[z;Mk0 Rp=abs(rp).^2;
{SJsA)9:# :N2E}hxk Rs=abs(rs).^2;
^2E hlK^) /Pk:4, Rn=(Rp+Rs)/2;
3htq[Ren \]8VwsP Tp=1-Rp;
'd/*BjNp) Q ]"jD#F Ts=1-Rs;
4d`+CD C G6V/S aD Tn=(Tp+Ts)/2;
9OyN i ]v\^&7pW figure(2)
X FQNr ` H -,TS^W subplot(1,2,1);
)t%h[0{{ F.?01,J=1 plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
F+uk AT
s]Gd-j legend('R_p','R_s','R_n')
zb :kanb- =IZ[_ /@ xlabel('\theta_i')
@}DFp`~5| k][ h9' ylabel('Amplitude')
PNpu*# Z` uE`r /=4 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
v L}T~_=3 &yU>2=/T axis([0 90 0 1])
\GL*0NJ ^ H3m\!h grid on
`\N]wlB2/b E!!
alc{ subplot(1,2,2);
,q#2:b<E !n<o)DsZR plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
CxDcY (.?ZKL legend('T_p','T_s','T_n')
m5*[t7@% SkHYXe"] xlabel('\theta_i')
. I==-| aGK@)&h$ ylabel('Amplitude')
-Sz_mr Wp[9beI*M title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
o=_c2m
()\jCNLT axis([0 90 0 1])
! E0! -UpY aSL6zye
, grid on
>zng J$ [T^?Q%h 
:( , mL2[ ? ?[ g}> 
8(lR!!=q
