利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
SRpPLY{:F a2 `|6M; 1、光疏射向光密
c /ABBvd| @G>Q(a*, clear
!&8HA D= LLm$y
close all
-c'~0g]< e8ZMB$byP n1=1,n2=1.45;
~OQ/ |ws CLX!qw]@ + theta=0:0.1:90;
dd@-9?6M ~xP4}gs1 a=theta*pi/180;
p: 8&&v~I x $
oId{; rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
^SH8*7l7 pz@wbu=($4 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
kc&MO`2 W\ f6-OR]R5 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
` p)$7! ~^pV>>LX| ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
*# 2]`G) pSlosv(6 figure(1)
a jyuk@ W|aFEY subplot(1,2,1);
n%Gk
{h5 Y< drRK! plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
Rr/sxR|0_ zw,=mpf3_ legend('r_p','r_s','|r_p|','|r_s|')
Qt+;b y&$ v@]t1 xlabel('\theta_i')
DU>#eR0 G \ZPmPu9^( ylabel('Amplitude')
*D$[@-7 )cd5i E:FO title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
BLs kUrPF iO_6>&( axis([0 90 -1 1])
hs m%o\ .G|9 :b grid on
#.kDin~ ! Nn%[J+F subplot(1,2,2);
Y^~Dr|5% cK(S{|F plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
I "<ACM Lupug"p0
legend('t_p','t_s','|t_p|','|t_s|')
-l# h^ vUU)zZB ~ xlabel('\theta_i')
ya^zlj\`0e !.nyIA( ylabel('Amplitude')
sF`ELrR \ ClvqI"Rd title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
^k7`:@
z0U FnFJw;:,{ axis([0 90 0 1])
3RyB 0
n `x0GT\O2- grid on
<$]=Vaq e
M T5bn Rp=abs(rp).^2;
Nhnw'9 wgb
e7-{ Rs=abs(rs).^2;
g_l=z`,8 'nO%1BZj+ Rn=(Rp+Rs)/2;
X !vBD E0Y>2HOuL Tp=1-Rp;
lSu\VCG quPNwNy Ts=1-Rs;
&2EimP /d\#|[S Tn=(Tp+Ts)/2;
/Dl{I7W ~RRp5x _ figure(2)
?'dsiA[ 6ofi8( n[ subplot(1,2,1);
NQx`u"= O_u2V'jy9 plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
HoIK^t~VT# ph;ds+b legend('R_p','R_s','R_n')
X_6h8 n}i - 9Ll'fbq xlabel('\theta_i')
l".LtUf- AP8YY8,
ylabel('Amplitude')
P 'dH* }H |H LU5=Y title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
PSM~10l, (")IU{>c6 axis([0 90 0 1])
>*hY 1@N1 GjmPpKIu\ grid on
Y30e7d* qr 2th> +M~A subplot(1,2,2);
Z?7XuELKV p%8v+9+h2 plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
= % O @%v + ~6Nq(kV legend('T_p','T_s','T_n')
3j]P\T oY#62&wk4 xlabel('\theta_i')
Aw38T w _.ny<r:g ylabel('Amplitude')
=%}++7# ]~!jf title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
nbVlP ]%RX\~Q. 4 axis([0 90 0 1])
0gs0[@ ?-y!FD}m& grid on
`pHlGbrW Y ~|C]O 
1Rrl59}5 }3"FQ/6C 
7~2/NU?
