利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
wXZ-%,R -D [ cH/Y2[ 1、光疏射向光密
6 Sl [ Dt\E4 clear
@%TQ/L^| m&%N4Q~X> close all
2cDC6rul 49#-\=<gt n1=1,n2=1.45;
mrbIoN==` K)14v ;@ theta=0:0.1:90;
|/s.P NP2 ~W#f,mf a=theta*pi/180;
MVj@0W33m ?y
'.sQ rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Q\r qG | - R::gm rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
F+6ZD 5/ E`s_Dr}K tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
42C:cl} ." PQ[TTLG\& ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
PY2`RZ/ @ y#MLxm figure(1)
z_H2 L"Z Q,4F =b subplot(1,2,1);
SB;Wa% .d fTv/n plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
#[si.rv-> h1L p:@:| legend('r_p','r_s','|r_p|','|r_s|')
Dxe|4"%^ ]@f6O *&= xlabel('\theta_i')
:.~a[\C@V< V.F 's(o ylabel('Amplitude')
dOhV`8l SN+Bmdu p title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Vy9n3W"FB1 d!d
3r W;A axis([0 90 -1 1])
)FP|}DCxQ G^R;~J*TDE grid on
^4[\-L8Lpq S ~_ % subplot(1,2,2);
\w:u&6,0O j\vK`.z plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
8x{vgx @M J.&q[ legend('t_p','t_s','|t_p|','|t_s|')
OBl8kH(b> B -K Of xlabel('\theta_i')
=j{jyl C e\dT~) c ylabel('Amplitude')
Z^|C~lp;n *"5a5.`%, title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
{$t*Mb0 ^Pf&C0xXv axis([0 90 0 1])
zbR.Lb x o%iL grid on
=oTYwU 7 L]?)2= Rp=abs(rp).^2;
^M9 oTNk2 ~
/[Cgh0 Rs=abs(rs).^2;
mx[^LaR>v So^`L s;S Rn=(Rp+Rs)/2;
)L!R~F
C g-^m \>B Tp=1-Rp;
Jv<)/Km` HLk"a -+' Ts=1-Rs;
</Id' ;|v 8. ^U6xA Tn=(Tp+Ts)/2;
z]YhQIU4n8 \,EPsQV0? figure(2)
>(rB[ZJ =tNiIU subplot(1,2,1);
5{a(
+ ' l! y
_P plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
;40m goN i6CYD legend('R_p','R_s','R_n')
` Y;gMrp jq]"6/xxb xlabel('\theta_i')
Yel(}Ny Mh|`XO.5I ylabel('Amplitude')
O)|4>J*B )r i3ds title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
l{U 3; O-5U|wA axis([0 90 0 1])
@>@Nu g2 gk 1S"H grid on
`Cf
en 8 5)7mjyo% subplot(1,2,2);
=)YYx8gR N)*e^ Nfb plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
mv.I.EL I0vn d7 legend('T_p','T_s','T_n')
X@&uu0JJ yTt (fn:; xlabel('\theta_i')
}hitU( 5t0 $tHwJ!<$& ylabel('Amplitude')
.4W>9
8 a.IF%hP0xo title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
X9A[
25n (&NV axis([0 90 0 1])
|Du13i4].& $ cK
B+} grid on
{U '&9_y
z~e~K`S 
]MLLr'6? K_|~3g 
0o$HC86w
