利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
_N&]w*ce -H;p +XAY 1、光疏射向光密
+r
2\v auP6\kp Me clear
h% T$m_ t/9,JG close all
#`9D,+2iB% #aV2+ `d n1=1,n2=1.45;
a;$'A[hq t!"XQ$g' theta=0:0.1:90;
umD[4aP~; R0!qweGi@ a=theta*pi/180;
l;gj] ,* Pxn;]!Z # rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
JZ .`ppp!:a4 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
5%E.UjC #r 4S% tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
M~1 n# !z.C}n5F ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
D9|?1+Kc + ^9;<>P figure(1)
}m6j6uAR6) "/-T{p;. subplot(1,2,1);
@Yy:MdREA uJP9J U
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
K{)YnY_E; -gP4| r8& legend('r_p','r_s','|r_p|','|r_s|')
D|1pBn.b]' #' =rv xlabel('\theta_i')
.RNr^ *AQ ;uC +5g` ylabel('Amplitude')
= yH#Iil "c S?t title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
h*qoe(+ZD G
39 axis([0 90 -1 1])
P|S'MS';: mQOYjy3 grid on
, DT =( 2@ (Qd3N( subplot(1,2,2);
J6H3X;vxQw UJ'}p&E plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
`euk&]/^.) r"MKkS EM legend('t_p','t_s','|t_p|','|t_s|')
:W/,V^x} T F !Lp: xlabel('\theta_i')
Ij_VO{]G'l I\0mmdi73 ylabel('Amplitude')
[]{g9CO dpBG)Xzoyv title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
A@wRP8<GKj %0PZZl5b axis([0 90 0 1])
$cLtAo^W 'ErtiD grid on
=~&Fq$$ |xTf:@hgHf Rp=abs(rp).^2;
` NC{+A 'CDRb3w}B Rs=abs(rs).^2;
5O9Oi:-!c c0Td a Rn=(Rp+Rs)/2;
XB &-k<C RoXU>a:nS Tp=1-Rp;
xi6Fs, 2S `nO!_3 Ts=1-Rs;
Uv6#d":f; ~%4#R 4& Tn=(Tp+Ts)/2;
89B1\ff &/7 AW(? figure(2)
N~ -N Q -IR9^) subplot(1,2,1);
*R*Tmo" edPnC
{?s plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
3ySP *J5 z,nRw/o legend('R_p','R_s','R_n')
a4
g~'^uC n8zh;vuJ xlabel('\theta_i')
Kfc(GL? BZqb
o `9 ylabel('Amplitude')
lqOv_ q P X](hc= title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
" [=Ee[/ cI3uH1;# axis([0 90 0 1])
yZ aQ{]" GYiUne $ grid on
,\5]n&T;r ,vQkv u z subplot(1,2,2);
J=/| iW (tEW#l'} plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
/^ v4[] &X^~%\F:2 legend('T_p','T_s','T_n')
t=~5 I > ]"^GRFK5 xlabel('\theta_i')
r]'AdJFt J$PE7*NU ylabel('Amplitude')
AKM\1H3U K}O~tff title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
7/(C1II.Q C+}uH:I'L axis([0 90 0 1])
K/A xojo K:P gkc grid on
V Ltb16| Tk/K7h^ 
}6eWdm!B )zt*am; 
A]B D2
