利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
NvjJ b-u 0:KE@ = 1、光疏射向光密
6I_Hd> 4 {5*|C-WWtG clear
0^'B3$> \bC X=E- close all
Kt>X[o3m, 9@ $,o M= n1=1,n2=1.45;
52zGJ I*
t<:D@J]a theta=0:0.1:90;
+a=
0\lpOy URcR a=theta*pi/180;
!k$}Kj)I N1 2:{U rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Hg}I]!B p<WFqLe(": rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
VVWM9x P,gdnV
^ tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
"f QRk P4
ul[zZ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
LE?u`i,e=+ Mu_i$j$vvP figure(1)
vd#,DU=p! ?3sT" r_d@ subplot(1,2,1);
5T4"j;_.BL {gf>* plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
Kn+=lCk QL"gWr`R legend('r_p','r_s','|r_p|','|r_s|')
(U.**9b; 5IMSNGS xlabel('\theta_i')
o}8{Bh^ a
ib}`l ylabel('Amplitude')
h\ ,5/ )Y b 5ie <s title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
i% (yk#=V ~
}KzJiL axis([0 90 -1 1])
AvRcS]@= l%"`{ grid on
7{BnXN[ y{rn-?`{ subplot(1,2,2);
m% bE-# ;R1B9-, plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
^JhFI* QJ4AL3
^6 legend('t_p','t_s','|t_p|','|t_s|')
jW5iqU"{* | 1V2tx xlabel('\theta_i')
.K9l*-e [= dt-K ylabel('Amplitude')
nlfPg-78B+ ;e;\q;GP title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
]xq::a{Oy n85r^W axis([0 90 0 1])
g5R,% 6 z}5<$K_U grid on
yKDZ+3xK] \y*j4 0 Rp=abs(rp).^2;
;/ Q6 i
M'u=H Rs=abs(rs).^2;
_ @prm Sc _om[VKJ d Rn=(Rp+Rs)/2;
{'IFWD. 5 W(k:Pl# Tp=1-Rp;
*( *z| 2 : HZ;Po Ts=1-Rs;
+Jw{qQ R/* 8IWw jyRr Tn=(Tp+Ts)/2;
* -z4 <LAa $r"A@69^RS figure(2)
XM!M%.0WS 5^F]tRz- subplot(1,2,1);
~PYFYjHC F50 JJZ plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
Yq0# #__ 4$i} Xk# 3 legend('R_p','R_s','R_n')
5Bd(>'ig_ 7)PJ:4IqS xlabel('\theta_i')
5DSuUEvWcL Q [:<S/w ylabel('Amplitude')
/| f[u s-w )OUU ]MUH title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
.Cfp'u%\; U?EG6t axis([0 90 0 1])
WY.5K
=} 7sguGwg) _ grid on
N?^_= KE@ CXa Ld7nMX subplot(1,2,2);
96)v#B?p RO$*G
jQd plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
@H4wHlb <{ # <5 8 legend('T_p','T_s','T_n')
vrO %XvXW C M^r|4 K xlabel('\theta_i')
h&$P y f\=6I3z ylabel('Amplitude')
z"3 H { A `OLB';D title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
fr]Hc+7 q06@SD$
axis([0 90 0 1])
D{3fhPNU<b *cQz[S@F grid on
:1NYpsd.i [[$Mh_MD 
D zE E:&*= i2U/RXu 
) =Y-f?o!
