利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
nk.Y#+1) 5/x "!J k 1、光疏射向光密
+|}R^x`z 'J^ M`/ clear
7sP;+G LhM{LUi close all
)|;*[S4 / hUuQDJ n1=1,n2=1.45;
&8'.Gw m} ?^Sk1 7G theta=0:0.1:90;
u] C/RDTH >! +. M9 a=theta*pi/180;
Y.I-h l1<r wMy$T<: rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
89@\AjI ~3}Gu^@ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
#Mo`l/Cwp IV{,'+hT tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
|8I #` ;t!n%SnK9! ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
(IX iwu qW]gp7jK4 figure(1)
p:M#F: U3r[ysf subplot(1,2,1);
!:e|M|T'I* >cwyb9;!kK plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
4)D #kP WelB+P2 legend('r_p','r_s','|r_p|','|r_s|')
(H *-b4]/ DO<eBq\O xlabel('\theta_i')
_`lj
3Lm0> H
<CsB ylabel('Amplitude')
QX (x6y>Q KDwz!:ye title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Nf LvK o8 U$a Eby. axis([0 90 -1 1])
1|Y(XB^os( 6cof Zc$ grid on
}*Qd] \fy 'Z2:u!E subplot(1,2,2);
EM/NT/ y7SOz'd plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
jB }O6u[% A7U'>r_. legend('t_p','t_s','|t_p|','|t_s|')
{9TWPB /> Mh C74G xlabel('\theta_i')
Lm+ !/e BqZ^ I eC$ ylabel('Amplitude')
T=cSTS!P;q l n.kEhQ3B title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
GF~ ^-5 *Yv"lB8 axis([0 90 0 1])
3{_A zL }C>Q grid on
dJ\6m!Mp ,2S
<#p! Rp=abs(rp).^2;
-Hzn7L FzmCS@yA Rs=abs(rs).^2;
>(z{1'f{ qG"|,bA
Rn=(Rp+Rs)/2;
iU^ 4a -Nmf}`_ Tp=1-Rp;
7= XQgbY/ Qi Wv Ts=1-Rs;
^&<~6y}U^ `7jm Tn=(Tp+Ts)/2;
O<4i)Lx2 .jMm-vox} figure(2)
s# w+^Mw$ Bp\io$(% subplot(1,2,1);
zF /}s_><* Ff>Y<7CQ
v plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
Gb 61X 6 jIE> t5 fy legend('R_p','R_s','R_n')
K-wjQ|*1 T(+F6 d=1 xlabel('\theta_i')
#%VprcEK o^RdVSkU; ylabel('Amplitude')
n ! qm cb&y8!ci~ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
uuj"Er3 1 x$CpUy{6 axis([0 90 0 1])
oSNB\G< 8 DTk<5mW~ grid on
..FUg"sSO iM2
EEC subplot(1,2,2);
3sGrX"0D MIa].S# plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
L:Y sAv QOuy (GY
legend('T_p','T_s','T_n')
GQqw(2Ub} fLI@;*hL0 xlabel('\theta_i')
/6{`6( p qNHS 1 ylabel('Amplitude')
N(Y9FD;H x+B~ t4A title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
7~\Dzcfk"P :'y axis([0 90 0 1])
oE#HI2X !Go(8`> grid on
Q m`f5-d _zh}%#6L 
:Tuy] ]k u-UUF 
! qrF=a
