利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
6( CDNMzj ^(HUGl_ 1、光疏射向光密
<r#eL39I hrG M|_BE clear
~Wo)?q8UY, ngohtB^] close all
?L&'- e@ },+wJ 1 n1=1,n2=1.45;
="wzq+ U R\@/U=iqR theta=0:0.1:90;
&|3
$!S nw>8GivO a=theta*pi/180;
u_(VE fs4 li ~d?> rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
InXn%9]p] J_ x13EaV0 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
nD5 gP $6OkIP. tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
aT>'.*\ ] l&iq5}[n& ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
}by;F9&B 5[0
O'%$ figure(1)
s<qe,' Y X/iT)R]b subplot(1,2,1);
o(DOQ Gl S;$-''o?9 plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
At7>V-f} F a'k0/_j legend('r_p','r_s','|r_p|','|r_s|')
8"sb; z!l.:F xlabel('\theta_i')
Vn*tp bz y W$0\E6<r ylabel('Amplitude')
?c [*:N( j@YU|-\qh title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
yE}}c{hSn GB$;n? axis([0 90 -1 1])
$ s9Vrw0Z dZi "$ g grid on
s2Mb[#:a" :<}=e@/~| subplot(1,2,2);
*:ZDd I 'V4D[H5 plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
N5a*7EJv+ ;W>k@L legend('t_p','t_s','|t_p|','|t_s|')
-$\+ '
\ ,%uo6% xlabel('\theta_i')
zuUW|r W[Ls|<Q ylabel('Amplitude')
&YF^ j2 Ney/[3 A title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
j'A_'g'^ mV3cp rRqv axis([0 90 0 1])
Pi]19boM. :]\([Q+a grid on
9I&xfvD, ;M)QwF1 Rp=abs(rp).^2;
;7} VBkH ,6-:VIHQ Rs=abs(rs).^2;
Tj:B!>> 0*f)=Q' Rn=(Rp+Rs)/2;
*MKO
I' "*In+ !K Tp=1-Rp;
&J+CSv,39 <
jJ Ts=1-Rs;
#ZB~ x6i6 kqFP)!37 Tn=(Tp+Ts)/2;
wB.&}p9p 9[<)WQe6M figure(2)
}H ^+A77v P)P*Xq r#: subplot(1,2,1);
&litXIvT> ?l9XAW t\ plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
4 o Fel.o o]4*|ARPs legend('R_p','R_s','R_n')
5 > [u ` Ff)8Q.m xlabel('\theta_i')
[ CQ+p!QZ 'Gj3:-xqL ylabel('Amplitude')
PvPOU" 4K\G16'$v title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
~E17L]ete fU/>z]K axis([0 90 0 1])
\NPmym_ 6J oKuI0-*mR grid on
'=b/6@& V<GHpFi0 subplot(1,2,2);
R!}H;[c dYJ(!V& plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
EJMM9(DQ7 0XE4<U legend('T_p','T_s','T_n')
Te"ioU?. g6j?,c|y xlabel('\theta_i')
Kp%2 k^U > qnko9 V ylabel('Amplitude')
0X6YdW _2X ua3~iQj- title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
LSL/ZvSP (C\]-E> axis([0 90 0 1])
"@V Y h4fJvOk|! grid on
E(>=rD /+ cr7 }^ s 
6' k<+IR f%][}NN)Xr 
XP!S$Q]D
