利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
$&"V ^@ NTS# sgP 1、光疏射向光密
!<UJ6t} hFk3[zTy clear
ov wQ2TuK f)g7
3= close all
F e.t/amS/ MB%Q WU n1=1,n2=1.45;
[tg^GOf ' rz "txN theta=0:0.1:90;
xGU(n _ Y {
.*y a=theta*pi/180;
E;d7 ch ?%VI{[y#> rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
M;0]u.D*= @x eAc0.^ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Y!WG)u5 Fbu5PWhlc tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
PG8 ^.)]M ?-tVSRKQ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
MwfOy@|N >7roe []-| figure(1)
Ja SI^go Yp0/Ab(v subplot(1,2,1);
d/}SAvtt u7xDau(c plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
<BSc* 9Q ]uspx [UIc legend('r_p','r_s','|r_p|','|r_s|')
gt YAH i n39t}`WIl xlabel('\theta_i')
ltkI}h,e ;SwM u@tg ylabel('Amplitude')
X*Z8CM_ 4 O:W#bx title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
~V/?H!r'{} t6BHGX {o axis([0 90 -1 1])
<" @zn H"Klj_<dH0 grid on
l c_E!"1 W5^.-B,(K subplot(1,2,2);
fK&e7j`qO {Qbv R*gv plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
Iyd?|f" '+
xu#R legend('t_p','t_s','|t_p|','|t_s|')
x!_<z'' NxVw!TsR xlabel('\theta_i')
7 S a1;%R N$x tHtz8" ylabel('Amplitude')
'TTUN=y }<PxWZ`,\ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
^1Xt]T`e R.LL#u}; axis([0 90 0 1])
s/:Fwr4q#a kD) $2I? grid on
K^rIG6 ^ FZ^6* Rp=abs(rp).^2;
3x~{QG5Gn n6(i`{i Rs=abs(rs).^2;
JlKM+UE :
@HBEt^! Rn=(Rp+Rs)/2;
M0| 'f' 5j0{p$ '9 Tp=1-Rp;
uY<
H#k O)kg B rB Ts=1-Rs;
f'q 28lVf xyH/e*a Tn=(Tp+Ts)/2;
b dJ+@r \<vNVz7.D figure(2)
v(l eide YpUp@/" subplot(1,2,1);
"+?Cz !i hh~n#7w~IR plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
O+=vEp( I
_i6-<c.Q legend('R_p','R_s','R_n')
`o79g"kxe AM} brO xlabel('\theta_i')
1Uf*^WW4 CK_(b" ylabel('Amplitude')
|3K)$.6~ !3oKmL5 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
sA9 &/p/ Gg\ G'QU axis([0 90 0 1])
7ou46v|m5 bXYA5wG grid on
E3a_8@ZB7 .bf<<+'o subplot(1,2,2);
Gjz [1d P6Bl
*@G plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
qQ7w&9r.M j%<}jw[2 legend('T_p','T_s','T_n')
4R>zPEo Je4Z(kj 0 xlabel('\theta_i')
xx*2? i 5&CDHc7Oj ylabel('Amplitude')
{dPgf Io2mWvu?5 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
oK FT? "[X lH|LdlX axis([0 90 0 1])
OMihXt[ RV-h IdAU grid on
# $\fh;!W r^HA a GpC 
N0n^ L|(R 7.<^j[? 
g"Q h]:
