利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
0(5qVJ12 l [ m_<1L 1、光疏射向光密
e}gGl<((g U0|wC,7" clear
Cjt].XR@ Gf7r!Ur;g close all
FBi&M Z` </9c=GoJ n1=1,n2=1.45;
$XyDw|z[ v araBFD theta=0:0.1:90;
=585TR;
V CC6]A M(i a=theta*pi/180;
/L`qOr2E ) fuAdG rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
| v> W j`u2\ ; rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
q.hpnE~#lh c8l\1 ce?7 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
nKwOSGPQt }P!:0w3 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
\= 6dF,V 3 cqc< figure(1)
o|KmKC n> dpdp0 subplot(1,2,1);
NljcHe}Qy !':y8(Ou plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
^6UE/4x!y d<_IC7$u> legend('r_p','r_s','|r_p|','|r_s|')
R -Tf9?) %
y` tDR xlabel('\theta_i')
>U:.5Tch'V 2Zq_zvKUt ylabel('Amplitude')
){J ,Z*& V[* <^% title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
LqQ&4I \ 0F
ey9c axis([0 90 -1 1])
=VSUE
Pq ]|Iczg- grid on
|'k7 ;UW ]D >\Z(b subplot(1,2,2);
}1Hy[4B(k\ 9Kc;]2m plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
Xj 6 ?,J AbMf8$$3SH legend('t_p','t_s','|t_p|','|t_s|')
ti ic>j\D !5zj+N xlabel('\theta_i')
&K"qnng/y (ap,3$ hS ylabel('Amplitude')
M;X}v#l|XI " V4ru&a title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Ql nI &o F`f8q\Fc axis([0 90 0 1])
8 }OII\ /:w.Zf>B9 grid on
sc# q03 Kb#Z(C9 Rp=abs(rp).^2;
K"=I,Vr: ?h-:,icR Rs=abs(rs).^2;
*edhJUT uWr vkLGN Rn=(Rp+Rs)/2;
tb'O:/ X./7b{P ax Tp=1-Rp;
V%w]HIhq X|p Ow," Ts=1-Rs;
\ci[<CP :&=`xAX- Tn=(Tp+Ts)/2;
{r[g .@ |("zW7g figure(2)
&n2dL->*# d j:6c@n subplot(1,2,1);
m^YYdyn]M UV)!z gP plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
oAq<ag\qV iJ EKLv legend('R_p','R_s','R_n')
kKNrCv@64d 05[k@f$n xlabel('\theta_i')
/c52w"WW ]n&Eb88 ylabel('Amplitude')
>}!})]Xw9 af;~<o a title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
GPU,.s"&( y/6%'56uF axis([0 90 0 1])
ood,k{ b#{[Pk,w9 grid on
Q$= X
?{ -@G,Ry-\t subplot(1,2,2);
Z/[ww8b. KIHr% plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
5(&'/ U^ <lHVch"(^$ legend('T_p','T_s','T_n')
[<A|\d'x H6%%n
X xlabel('\theta_i')
m\ (crkN
A\ };^Y ylabel('Amplitude')
x`g sD3C )Vnqz
lI5 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
vvcA-k? IO]%AL( .; axis([0 90 0 1])
.Fn7yTQ% VF] ~J=>i grid on
Ny )N ~jn~M_}K 
gKWsmx
T9O3$1eqfo
