利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
*t{c}Y&@ %:KV2GP 1、光疏射向光密
]t)M}^w rJZs
5g` clear
\.P}`Bpa WWW#s gM% close all
\t# 9zn> 3D{4vMm X n1=1,n2=1.45;
6bKO;^0 Hu8atlpo theta=0:0.1:90;
d[e: }1 no NF;zT a=theta*pi/180;
xeM':hD.o yI.H4Dl< rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
C2rj ]t <n TmZ-; rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
19 wqD IE0 |sq o+E tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
?w37vsN #r.` V!= ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
_$\5Z Ve >xV<nLf/ figure(1)
/~l/_Jct@G (GEi<\16[ subplot(1,2,1);
0t COb 9 Ge:-|*F plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
;%7XU~<a K= Z]#bm legend('r_p','r_s','|r_p|','|r_s|')
?e y&Un" &lPBqw xlabel('\theta_i')
7s8<FyFsjd ;5 Vk01R ylabel('Amplitude')
\As oee F NnLhJPh title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
X!rQ@F3 6+ $d axis([0 90 -1 1])
%rDmW?T frmqBC VJ: grid on
0^y@p&;/. A2|o=mOH subplot(1,2,2);
u@bOEcxK #i +P(xV plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
A_jB|<bjTP 9VUm=Z#` legend('t_p','t_s','|t_p|','|t_s|')
!*?(Q6 @ojg`!, xlabel('\theta_i')
noaN@K[GO 1R2o6`_ ylabel('Amplitude')
q BBYckS. W?^8/1U title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
]~\SR0 pcuMGo-# axis([0 90 0 1])
%< wQ +(<n |~ grid on
_)#=>$k\ i_0 ,BV C Rp=abs(rp).^2;
c3zT(FgO>N K/wiL69 Rs=abs(rs).^2;
@0vC v 6)#- 5m Rn=(Rp+Rs)/2;
g<2lPH
S< EB&P Tp=1-Rp;
fXu~69_ 9B+ zJ Vte Ts=1-Rs;
7O8V1T t q!hy;K`Jd Tn=(Tp+Ts)/2;
n b0 Py>4 D%jD 8 p figure(2)
p iYws<Q kMl @v` subplot(1,2,1);
"# -Nqq vZ[ $H plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
:7$\ X[ .#= j
<& legend('R_p','R_s','R_n')
`z-H]fU t<!+b@l 5 xlabel('\theta_i')
2fkIdy#n@ oh7#cFZZ0 ylabel('Amplitude')
io t.E%G O1x0[sy title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Y!Uu173 O)R 7t3t axis([0 90 0 1])
4V u'r? 'a;ini grid on
+Cn yK(V <qbZG}u subplot(1,2,2);
K k!6B JCw{ ?^F" plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
WFU?o[k-O l zfD)TWb legend('T_p','T_s','T_n')
_ `bS[%CJ {B FT xlabel('\theta_i')
JqI6k6~Q^ v87$NQvwQ ylabel('Amplitude')
M1AZ}b c0] CRZi;7`*1 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
2
) TG C rnB{Z4L axis([0 90 0 1])
g2g`,"T Cz'xG W{ grid on
D^66 p8t N<KKY"?I' 
b@=z rhQ `4VO&lRm 
0s 860Kn
