利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
TQ ID-I ~|l>bf 1、光疏射向光密
s )~Q@ze2 ]#3=GFs/ clear
#
?}WQP! l~ F,i n. close all
%S$P+B? sI9~TZ : n1=1,n2=1.45;
a l#yc = 1` theta=0:0.1:90;
t)$>++i SJy:5e?zk a=theta*pi/180;
;M@ /AAZ L.+5`& rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
K< Iz5+oD 2 )4{ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
?=Ma7 y ) ae/+Q8 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
HbZ3QW P } 1X,~y] ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
&*/X*!_HK 6 `X#<#_ & figure(1)
|_!xA/_U'T
<&$!;d8 subplot(1,2,1);
BR"*-$u0; ~3/>;[! plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
x<9 |t( =_.Zv legend('r_p','r_s','|r_p|','|r_s|')
i~ zL,/O8 V ,
)kw{]( xlabel('\theta_i')
jr~76 pX6OhwkTK ylabel('Amplitude')
0r: 8ni%cL <
8WS YZ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
)Bd+jli|s qC`"<R=GX axis([0 90 -1 1])
mP3:Fc _G MEq
()}7P grid on
^t9"!K HYW+,ts' subplot(1,2,2);
%<I0-o fv<($[ 0 plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
< <4G GO o?/N4$&5l legend('t_p','t_s','|t_p|','|t_s|')
N \A)P b>I -4 xlabel('\theta_i')
i"sVk8+o! n# Z6 d` ylabel('Amplitude')
Gh42qar` d3^LalAp title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
8l;0)`PU 0+cRUH9Ew axis([0 90 0 1])
]!h%Jlu _3hCu/BV grid on
{*>$LlL kpK: @ Rp=abs(rp).^2;
qXP)R/~OZ B>AIec\jG Rs=abs(rs).^2;
stoBjDS %Ljc#AVg Rn=(Rp+Rs)/2;
SQa.xLU .^P^lQT]> Tp=1-Rp;
fs~n{z,ja% OuS{ve Ts=1-Rs;
6mM J$FY+ Dzc 4J66 Tn=(Tp+Ts)/2;
%o+bO}/9 X3X~`~bAD figure(2)
9r\8 !R $0iz;!w subplot(1,2,1);
<~X=6 =NyzX&H6 plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
N-K .#5 $T ]1<3\G legend('R_p','R_s','R_n')
<fs2; e`#Gq0}8 xlabel('\theta_i')
Tff7 SEP l<z[)fE{uS ylabel('Amplitude')
,7/un8:%c r/3 !~??x title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
x1mxM#ql +zz9u?2C` axis([0 90 0 1])
dpE\eXoa, k];
<PF grid on
MH{vFA4:, D95 $ subplot(1,2,2);
R!}B^DVt l_^>spF plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
nb0<.ICF%R kDO6:sjR7 legend('T_p','T_s','T_n')
8q_3*++D .F)--% xlabel('\theta_i')
"eh"' Z pPG!{:YT ylabel('Amplitude')
;$[o7Qm5r a, &Kvh title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
69y TGUG3 K< Ct axis([0 90 0 1])
a; /4 ht bp$8hUNYz- grid on
J8? 6yd-7 glk
I9~ 
% ,~; w0 e/@udau 
t#Yyo$9
