利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
7-oH >OF^ 0.}Um 1、光疏射向光密
_l&`*
2d LvJ')HG clear
:%dIX}F ~h~K"GbC? close all
R;P>_ei(LK +$#YW5wy n1=1,n2=1.45;
S^@I 4Z |k wkikGQS theta=0:0.1:90;
i>9/vwe y@;4F n/ a=theta*pi/180;
8 oHyN o }LH>0v_<Y rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
JD^&d~n_ G\\z k rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
BX|+"AeF aW8Bx\q tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
J~5VL |c a g8yWFqE!T ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
W^2Q"c#7F YC:> ) figure(1)
,`/J1(\ nd G4Zs(:a subplot(1,2,1);
T/
CI?sn r sP3?.E plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
/ ;M0tP *'+OA6 legend('r_p','r_s','|r_p|','|r_s|')
?Uzs^rsb XelY?Ph,, xlabel('\theta_i')
V8>%$O
sw ,.` ";='o ylabel('Amplitude')
(Al.hEs' <^&NA<2 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
DT # 1*&- }Po&6^ axis([0 90 -1 1])
q P@4KH} e NIDK:q dR grid on
J!,5HJh1 W'\{8&: ! subplot(1,2,2);
pR8]HNY0 /xjHzva^ w plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
[R j=k)aBm cO2
.gQo' legend('t_p','t_s','|t_p|','|t_s|')
_9E7;ew XljiK8q; % xlabel('\theta_i')
N}wi<P:*) n
5NkjhP~Z ylabel('Amplitude')
dL<okw 06I'#: ] title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
JBdZ] x%ZjGDF m axis([0 90 0 1])
#.0^;M5Nh 8=Di+r grid on
}io9Hk>| w
.l 2 Rp=abs(rp).^2;
|b'fp1< / XpWqL9 s_E Rs=abs(rs).^2;
-IJt( X| F]3iL^v Rn=(Rp+Rs)/2;
|jW82L+!N% pB{QO4q n Tp=1-Rp;
y ";{k+ F#@Mf?#2
Ts=1-Rs;
3|%Q {U E9 {Gaa/{ Tn=(Tp+Ts)/2;
MP?9k )f \;p5Pagx0- figure(2)
z6 bTcs"7h Oh<[8S 7]C subplot(1,2,1);
H E+y1f] ghJ81 plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
IuTZ 2~ pL>Q'{7s3 legend('R_p','R_s','R_n')
xfqgK D> r 4jW =?| xlabel('\theta_i')
l%lkDh!$"
yowvq4e ylabel('Amplitude')
M cE$=Vv UNq !| title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
`!5 ZF@Q>e L #p-AK axis([0 90 0 1])
nCEt*~t9VE :{%6< j grid on
ofl3G
{u - O3^q. subplot(1,2,2);
|)i- c`x T30!'F(*, plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
WA~|:S+ qX`?4"4 legend('T_p','T_s','T_n')
!Pt| Hk dr ;JA2n\iP, xlabel('\theta_i')
?Jm/v%0O p+#J;. ylabel('Amplitude')
'0U+M{ &fq-U5zH title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
{
dw m>a " jmi
"O* axis([0 90 0 1])
:x] gTZ? 9*RfOdnNe grid on
eI8 ^T? }8GCO Y 
P~iu|j E}mnGe 
ItHKpTe r
