利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
L =8rH5 qkC /\![@ 1、光疏射向光密
KjQR$- A]DTUdL clear
0fYj4`4=n
V<j.xd7 close all
d20gf:@BM ?`4+c x}n n1=1,n2=1.45;
T8HF|%I t1%_DPD%W theta=0:0.1:90;
A7n\h-b |M+< m">E a=theta*pi/180;
&cu lbcz o";Z$tAJkC rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
rSJ9 v : WH= EPOR, rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
%wSj%>&-R 4!LCR}K tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
y>aZXa WoBo9aR ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Mz L1 Bh!M FD8 N"p figure(1)
-k"^o!p I hA* " subplot(1,2,1);
;]pJj6J&v >2Kh0r IH plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
PoT`}-9 Q V&D l_ legend('r_p','r_s','|r_p|','|r_s|')
9J?wO9rI X3V'Cy/sy xlabel('\theta_i')
6 C+"`(u%V 8f3vjK' ylabel('Amplitude')
J52
o
g4l :at$ HCaK title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
oDU ;E 48*Do}l] axis([0 90 -1 1])
k0Uyf~p~ )kkh JI*v grid on
afb+GA! Qu]z)";7 subplot(1,2,2);
,Bal F
3'9u# plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
%,k] [V XGkkB legend('t_p','t_s','|t_p|','|t_s|')
T"0,r $3: Xt'sQ} xlabel('\theta_i')
?IGVErnJJC a@0BBihz ylabel('Amplitude')
aM$W*- Y I;1lX
L title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
14"J d\M8 ?|ZTaX6A axis([0 90 0 1])
as>L[jyG/ :G<E^<M\)^ grid on
`8Lo { P ]Ty isaT Rp=abs(rp).^2;
.({smN,B Ey4z.s'-l Rs=abs(rs).^2;
P'O#I}Dmw< 8{Fsm;UsY Rn=(Rp+Rs)/2;
HO' '&hz / 0eYMG+K= Tp=1-Rp;
J:kmqk! @, W vvh Ts=1-Rs;
T0]*{k(FR w&x!,yd; Tn=(Tp+Ts)/2;
l}l Ii8 <bD>m[8, figure(2)
&|` C)6[C E{n:J3_X^d subplot(1,2,1);
+a*^{l}AST jr3ti>,xV plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
&c*^V L\ V,[d66H=N legend('R_p','R_s','R_n')
P(K>=O e~"fn*" xlabel('\theta_i')
d`(@_czdF ?Oc{bF7 ylabel('Amplitude')
3dDX8M? ? mhs$g> title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
tu5*Qp\ m{7(PHpw axis([0 90 0 1])
bl'z<S,
' z1Ov|Q` grid on
51QRM32Y
"k A*Vc# subplot(1,2,2);
UDL
RCS8i A.5i"Ci[ie plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
3ux0 Jr2yT \{EpduwZ legend('T_p','T_s','T_n')
" XT"|KF|D R+7oRXsu xlabel('\theta_i')
5j-]EJb > B>CB3 U ylabel('Amplitude')
CQpCS_M -<_Ww\%8M title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
k`x=D5s\ z-c}NdW axis([0 90 0 1])
|qoKO:B4-[ "hQ_sgz[Z grid on
;q1A*f\:# ":nQgV\ 9 
DU=dLE6-P; _Um d 
{$1J=JbE
