利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
WxW7q t ,O.3&Nz,c 1、光疏射向光密
d1T,eJ} /4t j3B, clear
7bqBk,`9 %F\?R[^5 close all
VK}fsOnj0 | B. 0TdF n1=1,n2=1.45;
GRGzP&}@ e@c0Wl Wa theta=0:0.1:90;
Kpb#K[(]& 4?0vso*X<: a=theta*pi/180;
H:!7 : Z&ZP"P4 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
iDMJicW!+F pV.Av rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
UH]l9 Aq$P dArDP[w tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
A{UULVp :/y1yM ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
e&7JpT UjaK&K+M? figure(1)
'#s05h r 9v?N+Rb subplot(1,2,1);
P9=?zh 6G. =jlt5 z plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
&Ts-a$Z7?S 8dV=[+ legend('r_p','r_s','|r_p|','|r_s|')
7#@c z5Su W4[V}s5u xlabel('\theta_i')
~vs}.kb 5Ycco,x ylabel('Amplitude')
u1t% (_h T;@;R % title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
K/A*<<r
~ =g/{% ; axis([0 90 -1 1])
A6GE,FhsG u @~Ji iC% grid on
eAX
)^q )\sc83L subplot(1,2,2);
"J+ 3w K#*reJ}K plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
w:s]$:MA8 .7i` (F) legend('t_p','t_s','|t_p|','|t_s|')
.{ 44a$) 2rmSo&3@s xlabel('\theta_i')
+6UVn\9Q +
` s@ ylabel('Amplitude')
m_=$0m J$ ^\\Tx*#i title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
~\=1'D^6CK d@ Y}SW TB axis([0 90 0 1])
{LX.iH9}l yyv8gH grid on
M7+nW ; e% _ cy 2z Rp=abs(rp).^2;
Z ZcEt &0 7]LF$] Rs=abs(rs).^2;
0 GB:GBhZ Xv< B1 Rn=(Rp+Rs)/2;
GytXFL3`: -:30 :oq Tp=1-Rp;
43={Xy r`AuvwHPs[ Ts=1-Rs;
q6_u@:3u 5WC+guK7 Tn=(Tp+Ts)/2;
rF aF
Bd Eq$&qV-?( figure(2)
=
QQ5f5\l `!D s6 subplot(1,2,1);
LEW'G"+ _@
*+~9%8p plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
;3\3q1oX u}!@ ,/) legend('R_p','R_s','R_n')
si&S% 4( Zm=(+
f xlabel('\theta_i')
OR}c)|1 )\6& 12rj ylabel('Amplitude')
#{k|I$ cgl*t+o& title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
A81ls#is #<e\QE '! axis([0 90 0 1])
i /~1F_ `}BF${vF grid on
oI}kH=<, U
f|>
(C subplot(1,2,2);
q4}PM[K?=\ s lg ]#Dy plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
S+'rG+NJ usEwm,b) legend('T_p','T_s','T_n')
t5'V6nv EI_ xlabel('\theta_i')
deM7fN4lTi ? [)}l9 ylabel('Amplitude')
%g1,N k Tj HwjRa title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
UW9?p} F ~zS Cg|"r axis([0 90 0 1])
}0u 8r` 0
;b[QRmy grid on
Pu"R,a p/U+0f 
mp9{m`Jb* dbOdq 
9lGOWRxR)
