利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
.)eJ L KDt@Xi 6|| 1、光疏射向光密
4JOw@/nE D4';QCwo clear
. W[[Z;D nMz~.^Q- close all
Kr;7~`$[ >9?BJv2 n1=1,n2=1.45;
[ij8h,[~] zYL</!6a[ theta=0:0.1:90;
_PI w""ssr ( C1@f!Z a=theta*pi/180;
?/EyfTex 6Vq]AQx rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
$ U~3$*R O( P
,! rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
^N{L au gWqO5C~h tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
x+mf QcSD& R78=im7 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
9!aQ@ J^ 4^? J BpBZ figure(1)
GQU 9UXe HOEjLwH subplot(1,2,1);
ch^tq",1> pONBF3H8 plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
m{~p(sQL #<^ngoOj legend('r_p','r_s','|r_p|','|r_s|')
o*d+W7l b[mAkm?9+1 xlabel('\theta_i')
g{]C@,W %`o3YR ylabel('Amplitude')
ao@"j}c &n5Lc` title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
CB7 6 z j[/~ I axis([0 90 -1 1])
&7gL&AY8 N[(ovr grid on
vCPiT2 G upD 2vtU subplot(1,2,2);
9} \{0;9 `2 vv8cg^ plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
t1y
hU"(J /1h
0 l; legend('t_p','t_s','|t_p|','|t_s|')
0Q2P"1>KT/ R0 g- xlabel('\theta_i')
)$h<9 e ;bC163[ ylabel('Amplitude')
ev bqBb21b 6NvdFss'A{ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
m'
LRP:9v L uLnmnmB axis([0 90 0 1])
%*>ee[^L , `ViF Y
grid on
R V*7? y%3 K#O8P+n5[ Rp=abs(rp).^2;
@0XqUcV xz @/^Cj Rs=abs(rs).^2;
h{dR)#)GF< Bt8 Rn=(Rp+Rs)/2;
Z+z x*(X i2PPVT Tp=1-Rp;
S\qYw(G H*l2, 0&W Ts=1-Rs;
rUb`_ W@ E7XFt#P. Tn=(Tp+Ts)/2;
[Rqv49n*V 3#d5.Ut figure(2)
z}Qt6na]- ;NyX9&@ subplot(1,2,1);
{V> >a *?k~n 9n5U plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
Lyx \ s; :/Zy=F9: legend('R_p','R_s','R_n')
t+R8{9L- Ol }^'7H xlabel('\theta_i')
u P'x{Pr) zi^T ?<t ylabel('Amplitude')
6[-N}) H#/}FoBiS title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
C
*\
=Q Qx9lcO_ axis([0 90 0 1])
S nMHk3 (\ D6 2xC5 grid on
g;bfi{8s_ R:=
%gl! subplot(1,2,2);
t%<@k)hd~G RhJ{#G~:% plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
LuM[*_8 w8| 38m legend('T_p','T_s','T_n')
B( ]=I@L=W :R/szE*Ak xlabel('\theta_i')
+Om(&\c(6 '.n0[2> ylabel('Amplitude')
bt=%DMTn =Q % F~ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
@` qhQ |Rh%w J axis([0 90 0 1])
mk)F3[ ke __}j
{Buk grid on
AE:IX P|c (lDbArqy 
VEolyPcsg& &zl=}xeA 
/V66P@[>
