利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
/kVy#sT| <S0!$.Kg*< 1、光疏射向光密
-zz9k= q i\rI j0+ clear
M42D5|tZc i4!n Oyk close all
{\ziy4<II m= fmf( n1=1,n2=1.45;
S-yd-MtQp ld[] f*RuW theta=0:0.1:90;
!HDk ] tQJ@//C\z a=theta*pi/180;
<z )MV
oa 4[.-
a&!} rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
&TWO/F+Y :yk Z 7X& rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
6aQ{EO-]'= ?-\K V ha tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Q:VD 2<2 N M ~e ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
K{>O. 5 +Mm0bqNN figure(1)
|_[mb(<| 7X}_yMxc subplot(1,2,1);
Punbw\9!d, o
ohf)) plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
^Z:x poz, a9jY^E '|n legend('r_p','r_s','|r_p|','|r_s|')
E4y "$U%. n7<<}wcV xlabel('\theta_i')
'A8T.BU NK#Dq&W+& ylabel('Amplitude')
sQihyq6U; nNilT J
title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
`O}.
.N]g F/EHU?_EI axis([0 90 -1 1])
Z^GriL ! >S' eXt grid on
[ V jd )% 'OX6e Y5 subplot(1,2,2);
Rb#Z\e}e- )s= z i" plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
B4|% E$1+ ^ACp_RM legend('t_p','t_s','|t_p|','|t_s|')
BTd'bD~EA V">Uh@[J_ xlabel('\theta_i')
Gh{k ~/B GF/p|I D ylabel('Amplitude')
4 a2&kIn (J^
Tss title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
!& 'xkw ` $yFur[97C axis([0 90 0 1])
/{kyjf[o&* ?ST}0F00} grid on
vCpi|a_eCu dNJK[1e6 Rp=abs(rp).^2;
p6HZ2Q:a VJR'B={ h Rs=abs(rs).^2;
hCxL4LrF y6PAXvv'{ Rn=(Rp+Rs)/2;
1
yzxA( B EI/OGp Tp=1-Rp;
ReK@~#hLY ZXU e4@qfl Ts=1-Rs;
%ri4nKGS 3mm`8!R Tn=(Tp+Ts)/2;
vy"Lsr3 Q^p@ 1I figure(2)
PO%]Jme EhD|\WLx! subplot(1,2,1);
/t9w%Y quKD\hL$ plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
pKH 4?F W#lvH=y legend('R_p','R_s','R_n')
FQ-( #[ e:,.-Kvzp` xlabel('\theta_i')
l.34h Czci6 Lz ylabel('Amplitude')
Kq S2 q!H 3JL title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
#VO2O0GR <7jb4n< axis([0 90 0 1])
Bs1-UI}+ <HzAh<_@F grid on
"FXS;Jf 0}^-, Q, subplot(1,2,2);
9nG] .@ H U1"t |KW8 plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
ROjjN W`W zgRP!q<9tt legend('T_p','T_s','T_n')
KqSa"76R 8uGPyH xlabel('\theta_i')
#P4dx'vm qA~ D*= ylabel('Amplitude')
f C^l9CRY 3w9
]@kU title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
c0qv11,: t JyqFFZ& axis([0 90 0 1])
'0O[ d N : te xl grid on
{fb~`=? \=@r1[d 
F8 89JSZ% V%0.%/<#5 
(O-)uC
