利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
J[]YG+r 2NI3&;{4 1、光疏射向光密
#*?5 `2Ff2D^ ? clear
aBol9`6 %mh
K1, close all
6g( 2O[n. Q9?/)&3Bu n1=1,n2=1.45;
@S<=Okrlj ~#@sZ0/< theta=0:0.1:90;
1R1J/Z*V/ Z2hRTJJ[A a=theta*pi/180;
O`~T:N|D xZbm,.v rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
G^j/8e Vcc/ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
wrw~J [%z~0\lu8 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
B9Y "J 9LzQp`In ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Hi<{c 7BdvJ" figure(1)
ysDfp'C, M}3>5*!= subplot(1,2,1);
MQw{^6Z>1 ="<+^$7:k plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
gmy_ZVU' V+cHL legend('r_p','r_s','|r_p|','|r_s|')
~h tV*R Yb6(KT xlabel('\theta_i')
pH'#v]" Y7]N.G3,] ylabel('Amplitude')
Bk~WHg>@G Ah)_mxK title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
)m
\}ITf X=mzo\Aos axis([0 90 -1 1])
xgnt)&7T Z/6B[,V grid on
O0';j!?X &*/8Ojv)9 subplot(1,2,2);
dX,2cK[aG 7bCTR2e\@w plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
``$At ,m ko$bCG% legend('t_p','t_s','|t_p|','|t_s|')
a~DR$^m N:\I]M xlabel('\theta_i')
gO1`zP!9Z j*rra ylabel('Amplitude')
Tg)Fr) )9{?C4NQ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
<Y9((QSM4 f[!N]* axis([0 90 0 1])
%}x/fq $ZH$x3; grid on
l7
j3;Ly _{TGO
jZr Rp=abs(rp).^2;
rhHX0+ q --NLm@; Rs=abs(rs).^2;
2.Th29] srw5&s(3X Rn=(Rp+Rs)/2;
M
g1E1kXe mc~d4<$`! Tp=1-Rp;
g\OPidY ?SkYFa`u* Ts=1-Rs;
)g(2xUk-y $DJp|(8 Tn=(Tp+Ts)/2;
FzXVNUMP =YR/X@& figure(2)
2_
< b6'%nR*f subplot(1,2,1);
A d=NJhzl 4?jXbC k~x plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
(|Y[5O) JGHQ_AI legend('R_p','R_s','R_n')
M%/ML=eLi ?r"m*fY% xlabel('\theta_i')
6,ylkf3 %1 9TJn%J$ ylabel('Amplitude')
#(?EL@5 j$4Tot title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
.zAB)rNc
| .fk!~8b[Q+ axis([0 90 0 1])
3&/5!zOg) <2HI. @^ grid on
G sm5L<rx w[qWr@
subplot(1,2,2);
sV@kQ:
!E'jd72O plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
o6P)IZ1 d/k&f5 legend('T_p','T_s','T_n')
Ie`kzssM J0~Ha u xlabel('\theta_i')
'3 xvQFg _S7GkpoK ylabel('Amplitude')
s_yY,Z: T_lexX[\ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
{*bXO8vi(( KA#4iu{ axis([0 90 0 1])
DI'wZySS^ Vf`n> grid on
-5l74f!i ?_3K]i1IS
f>g>7OsD] >j{phZ
Y .Na9&-(