利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
^{zwIH2I] ~i=/@;wRp 1、光疏射向光密
{(7Dz*0 eWhv X9
< clear
T=A 7f6` :nd
}e close all
P zzX Ds6 EN+WEMro n1=1,n2=1.45;
t'Nu^_# Qo0okir theta=0:0.1:90;
VX[{X8PkS @lc1Ipfk" a=theta*pi/180;
_( 0!bUs> OV;Ho rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
3`5?Zgp '_4u,
\SG rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
q F%wl a'
.o tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
s^AQJ{X oB&s 2~ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
@ J!)o d +O*S>0 figure(1)
`g6ZhG:W <Lyz7R6 subplot(1,2,1);
RCND|X WvAl !^{` plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
3S_KycE{ &+J5GHt@ legend('r_p','r_s','|r_p|','|r_s|')
H
\.E K Z Ipb 4{A&"\ xlabel('\theta_i')
*O$kF.3q [A
yq%MA ylabel('Amplitude')
%[B &JhT _ 0oZgt) title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
dG Qy=T: 4T G g`$e; axis([0 90 -1 1])
V. 1sb
pI
wAVO%8u grid on
pE^L Qi 5u~Ik c~ subplot(1,2,2);
t1n'Ecm ( "P&|e|7 plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
x1|5q/I x*}(l%[ legend('t_p','t_s','|t_p|','|t_s|')
[77]0V7 .^,fw=T|1 xlabel('\theta_i')
j8hb A^3cP , L ylabel('Amplitude')
\%m R* J+ :>F:G%(DK title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
R)nhgp(~ [LjYLm%< axis([0 90 0 1])
yJ/m21f ky#<\K1}' grid on
/+JCi6{sHS U1 1rj,7 Rp=abs(rp).^2;
%vxd($Ti" v $+G_ @ Rs=abs(rs).^2;
I]~xs0$4# Rn=(Rp+Rs)/2;
M GN*i9CE )WwysGkqol Tp=1-Rp;
3P^eD:)
w U 87VaUr Ts=1-Rs;
Qk Gr{ c ?V,a`6 Tn=(Tp+Ts)/2;
^ }U{O A /n@_Ihx figure(2)
v;JY;Uh|
$Z.c9 rY1 subplot(1,2,1);
0=t_ a]+ :}+U?8/"7 plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
E36<Wog 9oKRu6]D- legend('R_p','R_s','R_n')
K4,VSy1byI qf4|!UR{ xlabel('\theta_i')
iMk`t:!;#" zw\"!=r^ ylabel('Amplitude')
]9R?2{"K s^L\ h r title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
03$Ay _ 2 dWI/X axis([0 90 0 1])
$v-lG( &X}9D)\UJ grid on
XL EA |# ]L}<Y9)t subplot(1,2,2);
j(va# f# 0:v7X)St plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
P@k
;Lg" .zBSjh_=H legend('T_p','T_s','T_n')
1?E\2t&K 7QQ3IepP xlabel('\theta_i')
Nf<([8v;t } ..}]J;To ylabel('Amplitude')
JeWW~y`e?{ <ywxz1 i title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
/5U?4l(6[f 8eWb{n uJ> axis([0 90 0 1])
{3_M&$jN sJ3HH0e grid on
kIU"-;5tP YHvmo@ 
G]QD6b9~ B:^5W{ 
T( J&v|FK
