利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
l5l:'EY> kQRNVdiz 1、光疏射向光密
/<\>j+SC Xv|~1v%s7 clear
%1 9TJn%J$ #(?EL@5 close all
j$4Tot OAEJ?ik n1=1,n2=1.45;
H
:}|UW #O7|&DqF{ theta=0:0.1:90;
A?06fo, {()8 Wr a=theta*pi/180;
aF;QSI o;7!$v>uK rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
RM|<(kq XwOj`N{!H rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
N0,.cd]y` Mmq{]q~At tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Q4JwX=ZVj P#Z$+&)b)s ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
r77?s? t
\kI( G figure(1)
aF1pq O~.A} subplot(1,2,1);
EX7gTf# D"oyl`q plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
8sOM%y9M ]d&6 ?7 !> legend('r_p','r_s','|r_p|','|r_s|')
cxFfAk\,en />S=Y"a/7 xlabel('\theta_i')
~Y<x-)R U<|kA(5 ylabel('Amplitude')
m2AA:u_*j ~9#\+[ d_ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
iqig~fjK~ sa36=:5x- axis([0 90 -1 1])
0%7c?3# 0.4c|-n grid on
RcitW;{|Kg lwIU|T<4 subplot(1,2,2);
!n~p?joJ* kep.+t[ plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
/[`bPKr / Li?;H legend('t_p','t_s','|t_p|','|t_s|')
^|5bK_Z& y+\nj3v6 xlabel('\theta_i')
'oBT*aL m;)[gF ylabel('Amplitude')
3de<H=H' `%EMhk title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
K]yCt~A$ V)V\M6 axis([0 90 0 1])
0&E{[~Pv W]{mEB grid on
Gt{~u^< @jAuSBy Rp=abs(rp).^2;
sqFMO+ g|tnYN Rs=abs(rs).^2;
WBLfxr Ho9 a#9 Rn=(Rp+Rs)/2;
Rb\M63q k)oD Tp=1-Rp;
JL45!+ Q9=X| Ts=1-Rs;
fFudoIC [vV]lWOp' Tn=(Tp+Ts)/2;
DfgqB3U[ N|Mzj|i. figure(2)
=t\HtAXn[ v>CAA"LH subplot(1,2,1);
/DA'p [, usnbGkq plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
`O ?61YUQH [GJ_]w^}j legend('R_p','R_s','R_n')
EH<rUv63 /co^swz xlabel('\theta_i')
_PZGns,u yPL1(i; ylabel('Amplitude')
ZG+8kt!w {>'GE16x title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
eD5.*O me"}1REa axis([0 90 0 1])
Z_Ffiw(p Rqu_[M grid on
ya8MjGo T!u&r subplot(1,2,2);
" )T;3/c %H Pwu & plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
o<Hk/e~ RGvfy/T legend('T_p','T_s','T_n')
MgkeD ?_*X\En*3 xlabel('\theta_i')
w'E&w)Z] <XG&f ylabel('Amplitude')
ZT;$aNy BU],,t\ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
HE#IJB6BS? *\><MXx axis([0 90 0 1])
q9 !)YP+w KPc`5X grid on
CR3<9=Lv> t[/WGF&(R ?T]3I.3
2^ %X)w$}WH "@uKe8r|y