利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
5y~[2j B: JL$RBr 1、光疏射向光密
C eg6 o &^ \4&g5vE clear
9(I4x]` RcE%?2l D close all
& 3a+6!L[ %$}iM< n1=1,n2=1.45;
C^~iz
in BdYh: theta=0:0.1:90;
O|/tRkDMP{ &vn9l#\( a=theta*pi/180;
= G_6D AuCVpDH rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
ls5S9R 5 %SE g(< rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
LSs={RD2+p BZBsE
:(F tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
$ S49v ^m7PXY ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
)Qc$UI8L o] j* figure(1)
xWMMHIu g=%& p?1@E subplot(1,2,1);
83 n: h08 E5 dXu5+ye plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
6 q q7: ibq@0CR legend('r_p','r_s','|r_p|','|r_s|')
4>W`XH oMVwId f xlabel('\theta_i')
l1??b
F,
"x~C ylabel('Amplitude')
: ?V; 6IX!9I\sT title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
We ->d |= Dn[1BWM/7 axis([0 90 -1 1])
Dz{e@+>M anvj{ 1 grid on
YJy*OS_& u%pief subplot(1,2,2);
UdG oPzN A>?f bY2n plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
s*IfXv D)_Ei'+*l legend('t_p','t_s','|t_p|','|t_s|')
U,=K_oBAq 2
zy^(%a xlabel('\theta_i')
|}`5< a!6U >4\xcL ylabel('Amplitude')
H@z k8]_P qEAF!iB]L title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
#^ 9;<@M "jHN#} axis([0 90 0 1])
@ K9T )p] 7'|PHQ? S grid on
#RF=a7&F =(ZGaZ} Rp=abs(rp).^2;
3\{Sf /# &Yg/ 08* Rs=abs(rs).^2;
.H"hRYPC? -)oBh Rn=(Rp+Rs)/2;
tX251S L9T u >4 Tp=1-Rp;
Cu:Zn% =CoT{LRQ_ Ts=1-Rs;
K288&D|1WU {#,FlR2 Tn=(Tp+Ts)/2;
z' lNO| nU >-P0wowL figure(2)
q}5A^QX F(<8:`N;G subplot(1,2,1);
i.sq^] j 9QMn%8=j plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
:j ~5(K" akQH+j legend('R_p','R_s','R_n')
#E9['Jn Z 9HX+ sB
M xlabel('\theta_i')
rHlF& ET "|%9xGX|D ylabel('Amplitude')
$1+K}tP LzR iiP^q title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
D.(G 9H ~8*oGG~s axis([0 90 0 1])
7g)3 \C L2'd sO n grid on
1c#'5~ nB $V Wzv4^: subplot(1,2,2);
ImY.HB^& ^d80\PXz plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
M N#C2 qz (x
qA.(F legend('T_p','T_s','T_n')
SbZt\a 8 7X(]r1-+ \ xlabel('\theta_i')
{ yR)} r d|]O<]CG_ ylabel('Amplitude')
aLi_Hrb9 #[y<h3f] title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
<(4#4=ivP |A0$XU{ axis([0 90 0 1])
xXZKj |QLX.. grid on
V=& ,^qZ ,g_onfY 
`mjx4Lb m > (h_j 
<z uE=0P~%
