利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
CdT yUl ^kS44pr\Q 1、光疏射向光密
C*S% aR Ws+Zmpk% clear
K*ZH<@o4 nA!Xb'y& close all
c|kQ 3( 'G.^g}N1 n1=1,n2=1.45;
;+) M~2 = ,P?R
3 theta=0:0.1:90;
-:5]*zVp+- p6VHa$[ a=theta*pi/180;
Oy @vh>RY &H Pzm6.3 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
m4U7{sE CoZXbTq rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
OMr &f8 2MN AY%iT tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
' O\me G(4*e! aZ0 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
E !a5-SrR GM5 s~, figure(1)
L @ b8, . ] S{T subplot(1,2,1);
bt$+l [U^J STC'j1U plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
_` %z vFKt=o$ g legend('r_p','r_s','|r_p|','|r_s|')
?c vXuxCm .ZK|%VGW xlabel('\theta_i')
~}SOd<n)| ; 'o:1{Y ylabel('Amplitude')
C,+ Sv- aZe[Nos title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
kUl:Yj=& - PHVM=: axis([0 90 -1 1])
#>SvYP o'W[v0>
L- grid on
Q7ez?]j 6 'xO^2m+N; subplot(1,2,2);
c8o $WyO SI:+I4i plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
=Vgj=19X( 0FDfB; legend('t_p','t_s','|t_p|','|t_s|')
</K"\EU $/sQatic xlabel('\theta_i')
w3^>{2iqq 7p}.r
J54 ylabel('Amplitude')
6 u3$ .Q pg}9baW? title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
:v_H;UU 6J|E e 1Ez axis([0 90 0 1])
t"0Z=`Wi *):x K;o grid on
{9 >jWNx 5 WR(jl+M Rp=abs(rp).^2;
Pkr0| bs*
RVxlN* Rs=abs(rs).^2;
zSYh\g" I_Q*uH.Y 5 Rn=(Rp+Rs)/2;
M P Ma bK)gB! Tp=1-Rp;
oGzZ.K3 A DU6AlNx Ts=1-Rs;
C+ B `A9 w[>/(R7im Tn=(Tp+Ts)/2;
#+2|ZfCn% $-dz1} figure(2)
q4XS
E, _ ,h h O subplot(1,2,1);
Z4\$h1tl c}nXMA^^ plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
d7x d" 83*"58 legend('R_p','R_s','R_n')
M (,npW jk0Ja@8PK xlabel('\theta_i')
e]\{ Ia + L4_] ylabel('Amplitude')
DrD68$,QN m6P!#=a:l< title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
:/H fMJ \R >!HY axis([0 90 0 1])
|w(@a:2 kw :Mss"L820 grid on
D`=hP( y^ Ybr&z7# 2 subplot(1,2,2);
s=lk K / [ d #@N2 plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
,B>Rc# +tz^ &( legend('T_p','T_s','T_n')
d P(*IOO. h9)QQP P xlabel('\theta_i')
gi;#?gps `[JX}<~i ylabel('Amplitude')
$7{ V+> VWoxi$3v title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Ju$vuEO g ;U f? axis([0 90 0 1])
{`0GAW)q p q`Bg`c grid on
f`X#1w9 Ak!l}d 
J[:3H6%` A*I mruV 
@P:
