利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
.{y
uo{u a7 N!B' y 1、光疏射向光密
T)r9-wOq |2{wG 4 clear
M1KqY: 9E >j D[X5Y close all
/cU<hApK 0}_[DAd6 n1=1,n2=1.45;
5mb]Q)f9- "
Hd|7F'u= theta=0:0.1:90;
+\v?d&.f0 f bUr`~Y" a=theta*pi/180;
B^g ?=|{ j~*L~7 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
b0oMs=uBn XQw>EZdj_N rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
N<~ku<nAU 0jEL<TgC tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
)ZN|t?| JQ"U4GVp ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
>,td(= : 720D V +o figure(1)
*QH28%^ K'GBMnjD subplot(1,2,1);
w K}T`*k AvW:<}a, plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
MdW]MW{ )= =Jfn y legend('r_p','r_s','|r_p|','|r_s|')
-,U3fts >3Eo@J,?d xlabel('\theta_i')
R6(oZp h H-
$)3"K ylabel('Amplitude')
/FRm2m83 UeRj< \"Q title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
~;N^g4s OynXkH]0T+ axis([0 90 -1 1])
R~PA 1wDZ >G%o,9i grid on
s=;uc] 9g qw^uPs7Uw subplot(1,2,2);
[C'JH//q*t _WRFsDZ' plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
5rU[ T ir r4SXE\
G legend('t_p','t_s','|t_p|','|t_s|')
~Yy>zUH^X h-[VH% xlabel('\theta_i')
J|qZ+A[z Rx }$0c0 ylabel('Amplitude')
;'cN<x)% | 9J2NH|]c title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
rp ;b" q z62;cv axis([0 90 0 1])
0E3[N:s vw xXgk grid on
$adbCY \ tg<EY!WY Rp=abs(rp).^2;
lIW
}EM oY933i@l)P Rs=abs(rs).^2;
$}!p+$ =\H!GT Rn=(Rp+Rs)/2;
;6>2"{NW !1$])VQWI Tp=1-Rp;
e4 >_v( ' ;)D];u|_ Ts=1-Rs;
-;^j:L{ hpO`] Tn=(Tp+Ts)/2;
"eB$k40- YoBDvV":@ figure(2)
AP'*Nh@Ik( R#%(5-Zu#R subplot(1,2,1);
7/I, HxXp! i OW#>66d plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
5kCUaPu E87 Ww,z8 legend('R_p','R_s','R_n')
e4? >- lh7jux xlabel('\theta_i')
leTf &W % C| n9* ylabel('Amplitude')
;D FSzbF` #h`
V>; title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
`p2+&&]S ;:\<gVi: axis([0 90 0 1])
VY
| _d k E &2OD [iX grid on
UQ?XqgUM *U&0<{|T subplot(1,2,2);
=8AO: I|gB@|_~ plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
PQ5QA61 xhV O3LW' legend('T_p','T_s','T_n')
PgtJ3oq [} -GhP9; d xlabel('\theta_i')
K|zZS%?$ :XZU&Sr" ylabel('Amplitude')
1OCeN%4]Qk 9g'Lk P title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
g{OwuAC_ l;R%= P?'F axis([0 90 0 1])
<D<4BnZ( I*{4rDt grid on
CZud&
< e!N:,`R
5 
,?8qpEG~#+ *W,]>v0%T 
_x%7@ .TB
