利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
*T.V5FB0S Rhw+~gd*F 1、光疏射向光密
biU
?>R
*9`k$' clear
3@1$y`SN &
sXMB close all
MFipX E! hb> uHUb& n1=1,n2=1.45;
c4bv Jy 8 a; &0u> theta=0:0.1:90;
w ,X J8+B yvnvI y a=theta*pi/180;
g3Ul'QJ nk ;+L rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
OJ.oHf=K! V8Z@y&ny rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
h|<;:o?yh ly@%1 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
rAD4}A_w Yfy";C7X ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Ij9=J1c4 E_{P^7Z|Jg figure(1)
$-\%%n0>6 Gbv bGEG subplot(1,2,1);
BLm}mb#/{ @scSW 5+ plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
>VG*La' c BT
f legend('r_p','r_s','|r_p|','|r_s|')
I,r0K] `*i:z' xlabel('\theta_i')
!.O;SG }$ a *XY1 ylabel('Amplitude')
Q*8 x Bi1 Iomx"y]9 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
6 I43a1[s f3<253 1/} axis([0 90 -1 1])
VuH }@ +KIBbXF7 grid on
<W*6=HZ' m=w #l>! subplot(1,2,2);
zJOyr"B'8 ^xr & E plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
,,?X Gx &C#?&AQ legend('t_p','t_s','|t_p|','|t_s|')
tnq Zl S ifmX<'(9 A xlabel('\theta_i')
{H
3wL p# >,{ ylabel('Amplitude')
6f+@@=Xc 8C I\NR{x8 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Z,>owoP4 )?( _vrc< axis([0 90 0 1])
}G" bD8+ ]@l;;Sp grid on
I1p{(f J G}nj
71=H Rp=abs(rp).^2;
WqHp 23 D$`$4mX@hP Rs=abs(rs).^2;
Io('kCOR; On,z # A Rn=(Rp+Rs)/2;
6g)CpZU w4RtIDW: Tp=1-Rp;
Z0M|B v9_ -8S Z}J Ts=1-Rs;
3RI %O CGF )0=H)k0 Tn=(Tp+Ts)/2;
<V|\yH 9 -r[O _[g w figure(2)
R-Y 7I ) Lo hB,? subplot(1,2,1);
^j1i CL! :S+Bu*OyH plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
NH'QMjL) ?VyiR40-Cx legend('R_p','R_s','R_n')
[6nN]U~ Y !]&+g'aC3 xlabel('\theta_i')
/@", 5U# 98CS|NE e ylabel('Amplitude')
P[H 4Yp Bey9P)_Of title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
[MeFj! ( ~Vc`AcW P axis([0 90 0 1])
0 R>!jw 6Zv-kG grid on
= @ 1{LF; t;t;+M|W subplot(1,2,2);
Iz!] LW Z
jXn,W]~ plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
xBE}/F$ 45 VzWH9%w legend('T_p','T_s','T_n')
W'v
o? k
jx<;##R8 xlabel('\theta_i')
7<LCX{Uw kU{a!ca4 ylabel('Amplitude')
}?9 A:& i8 =+ <d title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
.F0Q< s9 Q|7 m9 ~ axis([0 90 0 1])
w[u >*I (?[%u0%_ grid on
z :A_ :\%ZTBLL 
_DYe<f. xQNGlVipZ@ 
)vG xF}I3
