利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
2 Il8f x)SralWb 1、光疏射向光密
{$d <1y^ VWx]1\ clear
f'X9HU{Cz a 7#J2 r close all
mT @ nn, ,:)`+v< n1=1,n2=1.45;
+m)q% I> 5G[x }4U theta=0:0.1:90;
7" 4z+w k?*Knf Vh! a=theta*pi/180;
w##^}nHOR >oGiIYq rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
ST1'\Eo ?A?F.n` rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
#Ex p51 7I9aG.; tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
(YPi&w~S KkcXNjPVS ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
0|c}p([~ os4{0Mxu figure(1)
h !gk s-0 m&*JMA;^ subplot(1,2,1);
I9?Ec6a_ Fh8lmOL;? plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
w(9*7p p E5</h"1 legend('r_p','r_s','|r_p|','|r_s|')
*bd[S0l / 3!fA=+ xlabel('\theta_i')
>yB( lKV jQ\zG J3 ylabel('Amplitude')
j$Ab >}g] zmI] cD@G title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
k^\pU\J i#/]Ks Sp axis([0 90 -1 1])
- +> 1r :| +Qe e grid on
S >yLqPp $q$7^ r@ subplot(1,2,2);
JH8}Ru%Z ]QVNn ?PA8 plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
a 5Y IUVCv ? *>]")[> legend('t_p','t_s','|t_p|','|t_s|')
9 QQyl\ gNYqAUG5 xlabel('\theta_i')
up>c$jJ ;G\8jP '
ylabel('Amplitude')
q1?&Ev^ ZfrVjUB title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
-fS.9+k0/ \C,p
WW axis([0 90 0 1])
D^ Jk @<* ^r6!l. grid on
m_+sR!\H8 F + m }#p Rp=abs(rp).^2;
s EM Q +{<#(} Rs=abs(rs).^2;
Dre2J<QL $+p ?Y)h . Rn=(Rp+Rs)/2;
Fz#X= gmG f| _u7"OX Tp=1-Rp;
t>=fTkB _g%TSumvq< Ts=1-Rs;
^9Qy/Er' 5GA C`}} Tn=(Tp+Ts)/2;
WHXj8*]6 AIMSX]m figure(2)
3?2 FP|G8 H0m|1
7 subplot(1,2,1);
@,Gxk
.p0Clr! plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
h/?l4iR* 7X@mSXis legend('R_p','R_s','R_n')
.*m>\>Gsgw *na?n2Yzt xlabel('\theta_i')
<sK4#!K q 9Y9w( ylabel('Amplitude')
[
ol9|sdu T,aW8| title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
=[b)1FUp ]vwW]O7 axis([0 90 0 1])
n]E?3UGD@W , ]bB9tid grid on
zR 2B-
&]H ]a5 f2l E subplot(1,2,2);
C74a(Bk}H o2<#s)GpY plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
jv&*uYm M#(+c _(r legend('T_p','T_s','T_n')
;4Y%PV z~D Z& ;uh_EC xlabel('\theta_i')
6I@h9uIsze ;L']e"G ylabel('Amplitude')
uj :%#u >=[(^l title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
B95B|tU>. ")d`dj\o axis([0 90 0 1])
0`^&9nR z^T`x_mF grid on
Q~Hy%M%R3 :j + ZI3@ 
#Y$hNQQ$F iof-7{+3_ 
PYGRsrcFd#
