利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
<)H9V-5aZ p8O2Z? \ 1、光疏射向光密
Q(?#'<.# +~$ ]} % clear
7>%8eEc *L^,| close all
Wqw1J=] BX7kO0 j n1=1,n2=1.45;
zwjgE6 E{`fF8]K theta=0:0.1:90;
XNkn|q2 6A- |[(NS a=theta*pi/180;
qR8Lh( "i 2HA:"v8 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
14yv$, \ Gvm9M rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
[RhO$c$[\ g}cq K tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
}&J q}j ~B?y{ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
^hM4j{|&M 7R\<inCQ figure(1)
@qAS*3j O_7| C\] subplot(1,2,1);
]{@-HTt W hy`ziks plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
JK5gQ3C[ %7. 30CA|# legend('r_p','r_s','|r_p|','|r_s|')
hHnYtq h*](a_0 xlabel('\theta_i')
5U$0F$BBp + [mk<pQ ylabel('Amplitude')
Kn1a>fLaJ_ W^l-Y %a/o title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
1oGw4kD^x >|UOz& axis([0 90 -1 1])
fuySN!s }K|oicpUg grid on
LZY"3Jn[nQ /a4{?? #e subplot(1,2,2);
64tvP^kp kt:!
7 plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
F={a;Dvrn f .)O2= legend('t_p','t_s','|t_p|','|t_s|')
$=4QO ^ [@ , xlabel('\theta_i')
cbTm'}R(G a&? :P1$ ylabel('Amplitude')
D*d]aC Q\Vgl(;lX title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
G4" F +%. |yPu!pf l axis([0 90 0 1])
sfl<qD+? xJ .M;SF4 grid on
=t?F6)Q 6Z"X}L,* Rp=abs(rp).^2;
Z,PPu&lmE/ _H@DLhH|= Rs=abs(rs).^2;
6D3B^.r j] j0q &&9/Jj Rn=(Rp+Rs)/2;
H<+T R6k< vnuN6M{ Tp=1-Rp;
U 6)#}
dSHDWu& Ts=1-Rs;
b&U62iq |2A:eI8 ^ Tn=(Tp+Ts)/2;
u=e{]Ax#} )vlhN2 iv figure(2)
0w7DsPdS S,8e lKH4 subplot(1,2,1);
G' 1'/ _lq`a\7e plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
cFX p xskz) kk legend('R_p','R_s','R_n')
MF'JeM;H N?8!3&TiV xlabel('\theta_i')
#GFr`o0$^ <1TA w. ylabel('Amplitude')
#KvlYZ+1 'V>-QD%1 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
uPvEwq*
C CTmT@A{ axis([0 90 0 1])
Dw"\/p:-3 r9XZ(0/p grid on
|DwZ{(R"W +b 6v!7_ subplot(1,2,2);
Q,Eo mt [nh>vqum plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
/x *3}oI E{vbO/|kf legend('T_p','T_s','T_n')
noj0F::m`j l U] nd[x xlabel('\theta_i')
4<v&S2Yq x?<FJ"8"k ylabel('Amplitude')
8 zb /xP> |uJ%5y# title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
*n!J=yS _yT Ed"$
axis([0 90 0 1])
| V(0GB w32y3~ grid on
J[kTlHMD 0*v2y*2V 
J .%IfN q.}CU. dp 
19] E 5'AI
