利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
8:Jc2K *| YR 8f 1、光疏射向光密
a)M#O\i` pH!8vnoA clear
'sAs# ^sZ, (sc{G close all
gl"1;C 8` f=E h n1=1,n2=1.45;
(Z};(Hn JdA3O{mT) theta=0:0.1:90;
\}}) U# Ru^j ~Cj5 a=theta*pi/180;
}^3ICwzm ^U@E rc#d rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
hX sH9R
7S]akcT/ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
`Ot;KDz T,Zfz9{n tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
x4bj ?=+ r081.< ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
i_qR&X 095Z Z20 figure(1)
1W2hd!J7C "G
@(AE( subplot(1,2,1);
TY h_uox6 B[6y2+6$0 plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
aJ}Cq k H$n{|YO ` legend('r_p','r_s','|r_p|','|r_s|')
JRl`evTS 3XomnL { xlabel('\theta_i')
h\qM5Qx+Q
MfNguh ylabel('Amplitude')
} h0
) ~Uw<E:?v title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
#}{1>g{sXt jZvQMW axis([0 90 -1 1])
Yy:Q/zw o %?[H=v(b grid on
x_TtS| L[FNr& subplot(1,2,2);
kdH P
v=/U e^ygQ<6 % plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
#4<Rs|K F( Iq8DV legend('t_p','t_s','|t_p|','|t_s|')
B Ji }vLK-V v xlabel('\theta_i')
<rs"$JJV BVpRkUC" ylabel('Amplitude')
;^j 2>Azn ?u]%T]W title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
0au\X$)Q H7"m/Bia axis([0 90 0 1])
=5sUpP V( g [~"c} grid on
3Vj,O?(Z $'2yPoR Rp=abs(rp).^2;
pK~K>8\ $ zvqjT:> Rs=abs(rs).^2;
ekj@;6
d] ^SouA[ Rn=(Rp+Rs)/2;
.(8sa8{N m
pWmExQ Tp=1-Rp;
p kR+H| w Kq-|yf, Ts=1-Rs;
'h&"xXv4| @Gj|X> 0 Tn=(Tp+Ts)/2;
3K57xJzK HjO- 6F#s figure(2)
l}$Pv?T,2 ag$mc8-p[ subplot(1,2,1);
$0MP*TFWa ? }2]G'7? plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
D`b H_1X n^ Qt !~ legend('R_p','R_s','R_n')
X Y;cz b uRK\ C xlabel('\theta_i')
{ T]^C 6^ ]Y]) ylabel('Amplitude')
EfR3$sp ;n Pjyu'g title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
'o#ve72z1 QJaF6>m axis([0 90 0 1])
:Q\{LB c $Y!$I .+ grid on
CmNd0S4v xsq+RBJi subplot(1,2,2);
os]P6TFFX? Rmrv@ .dr! plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
Olg@ Ri { L(Q|bB legend('T_p','T_s','T_n')
J0e ~s ,rB"ag ! xlabel('\theta_i')
M 8
E8r
"`Q .z~ ylabel('Amplitude')
<MlRy%3Z 2sJj -3J title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
IQFt4{aK3 S?bG U8R5 axis([0 90 0 1])
CV~\xYY 0{/P1 grid on
l ;I)$=={= /N`l
z>^~ 
2\xv Yf- G &'eP 
)19As8rL/o
