利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
A{[joo -|ho
8alF 1、光疏射向光密
/k$h2,O"* 3D"2yTM( clear
#y<KO`Es U7)#9qS4 close all
5r *5Co+ $>PXX32 n1=1,n2=1.45;
S'Hb5C2u ne]P -5 0 theta=0:0.1:90;
NUlp4i~Q LW={| 3} a=theta*pi/180;
xJphG )w?DB@Tx rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
laR cEXj
PTU _ <\ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
qXOWCYqs "!ZQ`yl tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
9g7d:zG -/x= `S * ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
|K|[>[?Z/ (=2-*((&(A figure(1)
WpPm|h iKJ-$ x_5 subplot(1,2,1);
6x'F0{U
n[vwwY plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
>@b]t,rrK |fL|tkGEa legend('r_p','r_s','|r_p|','|r_s|')
GI+ x,p aVg~/ xlabel('\theta_i')
: 3J0Q *oby(D"p ylabel('Amplitude')
!"v[\||1 rr@h9bak;g title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
,Wv@D"4? [@VM'@e7 axis([0 90 -1 1])
j W,b"[ /C[Q? grid on
K nl`[Nl VjI=5)+~ subplot(1,2,2);
>>cb0fH5 J?wCqA plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
GI se|[p Q9yIQ{>H[ legend('t_p','t_s','|t_p|','|t_s|')
9QQiIi$74U ~ u)} / xlabel('\theta_i')
!^ y H]v R ]h3a :ic ylabel('Amplitude')
kB)u@`</mV diz= |g=w title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
a2!U9->! GM~Ek] 9C% axis([0 90 0 1])
` !udU,|N Y>/T+ub grid on
=bBV
A0y vGCv J*4! Rp=abs(rp).^2;
afw`Heaa2( |erG cKk Rs=abs(rs).^2;
u ?-|sv* o7r7HmA@ Rn=(Rp+Rs)/2;
-;:.+1 ]\C wa 9 Tp=1-Rp;
>\ 7M f@c 22T \ -g{ Ts=1-Rs;
0h=}BCb+i r4isn^g Tn=(Tp+Ts)/2;
}@y(-7t `SH14A* figure(2)
O"GuVC}B ^Q\Hy\ subplot(1,2,1);
` pY yr/ } Q?a6(4 plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
\{a!Z&df /szwVA legend('R_p','R_s','R_n')
9Q.@RO$%C 45,): U5 xlabel('\theta_i')
Op'&c0l 5|jsv) M+ ylabel('Amplitude')
9+iz+ Y#5v5
title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
`53S[8 Ei 7Oi!1 axis([0 90 0 1])
q'Nafa&a) kz*6%Cg*~ grid on
5SMV3~*P 2<T/N subplot(1,2,2);
rJT YCe1* Eq?U$eE plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
bZ>dr{%%e O'Q,;s`uC legend('T_p','T_s','T_n')
{pc (b WfXwI 'y xlabel('\theta_i')
cy(4g-b]@e r;9 V7C ylabel('Amplitude')
x*sDp3f[* Y?qUO2 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
6oI/*`> ICEyz|
C axis([0 90 0 1])
mKqXB\< Zq~Rkx grid on
%g~&$oZmq Ne)3@? 
~:
fSD0 y`i?Qo3 
AtUt E#K
