利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
'iMI&?8u xJ>hN@5}i 1、光疏射向光密
lEZ[0oa Tr/wG clear
^lRXc .c z 2;[75(l6|} close all
WU#bA|Cf UH%?{>oRh n1=1,n2=1.45;
in # qV PM =I theta=0:0.1:90;
[JMz~~ F w:'$Uf8] a=theta*pi/180;
x gaN0! '4dnC2a] rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
S][: b }[D~#Z!k rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
8xg:ItJaA0 _*bXVJ
] tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
c
!P9`l~MQ e
d4T_O; ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
f:"es: Fb L V33vy figure(1)
>\} 2("bv \RMYaI^+; subplot(1,2,1);
kt X(\Hf! E;V W6[M plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
wzo-V^+q i0p"q p legend('r_p','r_s','|r_p|','|r_s|')
G0> ' H1 Z tPC8/ntP8 xlabel('\theta_i')
jW2z3.w
6=A ++H @ ylabel('Amplitude')
:$/lGIz +Z)||MR" title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
oc\rQ? /kx:BoV axis([0 90 -1 1])
I7B fA,mZ7 zMtK_c cQ grid on
vWe)c J Ik;~u8j1e subplot(1,2,2);
Zn #ri 8S i1 SP plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
"St, 4 b A3'i
- legend('t_p','t_s','|t_p|','|t_s|')
K.l?R#G`,F e
*;"$7o9 xlabel('\theta_i')
^x 4,}' ( x3wyIio* ylabel('Amplitude')
K Qub%`n ZW+{<XTof4 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
z?j~ 2K<4 7:M%w'oR axis([0 90 0 1])
{*jkx ,| x-y=Jor grid on
nV:.-JR lmr {Ib2a Rp=abs(rp).^2;
H[]j 6D jn)~@~c Rs=abs(rs).^2;
eSa ]6 xOTm-Cm9L Rn=(Rp+Rs)/2;
?>RJ8\Sj [uAfE3 Tp=1-Rp;
iKp4@6a n Sw'DS Ts=1-Rs;
2!]':(8mR F P
mLost Tn=(Tp+Ts)/2;
c/T]=S[ 8 {Zgvqbb figure(2)
f*oL8"?u& + ` Em& subplot(1,2,1);
G _42ckLq >K50 h plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
[(#ncR8B J, {sRb% legend('R_p','R_s','R_n')
<]CO }r
$ZQP f xlabel('\theta_i')
R['q BHQ? l6l )M ylabel('Amplitude')
] MUuz'< zpg*hlv title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
}p 8a'3@Z KS(s<ip| axis([0 90 0 1])
g<UjB Z[0xqGYLB grid on
_5.7HEw>/ Q4c>gds` subplot(1,2,2);
"'6KQnpZ - I4@` V plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
^P&y9dC. q'K=Ly+ legend('T_p','T_s','T_n')
lv$t p,+ z4{|? 0=C xlabel('\theta_i')
f>p;Jh{2fn tOp>O oD ylabel('Amplitude')
RE *UIh*O 2,NQ(c_c$ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
0%'&s)# g w }t.3} axis([0 90 0 1])
K-p 1v!IC Sfh \ 4h$H grid on
AbLOq@lrK LRB#|PW 
sbk WJy ZHw)N &Qn 
5&+
qX
2b
