利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
Q}cti / jT/P+2hMW 1、光疏射向光密
6=qC/1,l fP6 \Ur clear
U ^nv) q<Qjc close all
e"*1l>g ]')y(_{ n1=1,n2=1.45;
r)Vpt
fg; A@L r(L theta=0:0.1:90;
Qv
g_|~n U#bmM H a=theta*pi/180;
R NQq"c\ ; <mcvm rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
q [+ h ~) _CTg")0o rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
c`&g.s@N\ .C &kWM&j tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
fUfd5W1" O} (sn ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
cM CM> *X ';!-a] N figure(1)
:_]0 8 `.d w G3R subplot(1,2,1);
kr@!j@j$ +v'2s@e`
# plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
F FcIOn 8aC=k @YE legend('r_p','r_s','|r_p|','|r_s|')
V#|/\-@ >I<}:= xlabel('\theta_i')
IOF!Ra:w 8 R7w$3pp\ ylabel('Amplitude')
_ker,;{9C ` AD} 6O+x title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
'rS \9T 7+}WU 4 axis([0 90 -1 1])
G mE`YW mlLqQ< grid on
$CJf 0[| " FhC"}N subplot(1,2,2);
z@o6[g/*Q *M*WjEOA plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
F6{/iF ,grx'to(X legend('t_p','t_s','|t_p|','|t_s|')
Q+wO\TtE J] w3iYK xlabel('\theta_i')
T8)X ?>CIW mdQe)> ylabel('Amplitude')
a7uL {*ZR `IJ)'$pn title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
ya5HAs Yk )fBPHr axis([0 90 0 1])
MxUbx+_N yPe9KN_ grid on
2{Dnfl'k BOR$R}q Rp=abs(rp).^2;
;DhAw 1 B0A y Rs=abs(rs).^2;
fAz4>_4 E.sZjo1 Rn=(Rp+Rs)/2;
);y ZyWDV ` s Im&.d Tp=1-Rp;
.qKfhHJ W`c$2KS?DO Ts=1-Rs;
u"%D ; CB,2BTtRE Tn=(Tp+Ts)/2;
I<,~>'cq. I Z*) figure(2)
?Q+*[YEJ5 "L `BuAB subplot(1,2,1);
B2*>7 kc_s r..f$FF)\ plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
f2pA+j5[ *JZ9'|v_H legend('R_p','R_s','R_n')
tS5J{j>T *GY8#Az xlabel('\theta_i')
5Ma."?rW
}EHL
}Q ylabel('Amplitude')
lq]8zm<\)] -P-8D6 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
'vu]b #l3 zBwqI JfM axis([0 90 0 1])
&OE-+z m\CU,9;;( grid on
,quUG S ^c9_ F9N subplot(1,2,2);
%Ie,J5g5 R
>SZ E" plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
s]`6u yW" k ka{u[ruA legend('T_p','T_s','T_n')
qmGHuQVe \ZhkOl xlabel('\theta_i')
~ ;pP@DA i92Z`jiR ylabel('Amplitude')
,3eN& WlY\R>x# title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
v$x)$ /]n |
2<zYY axis([0 90 0 1])
[Mc H l1a "X._:||8
grid on
d2US~.;>l J#4pA{0 1w 
\L$]2"/v- [ + CFQf> 
Z,\(bW
qF
