利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
.|-l+ r: c@17 1、光疏射向光密
3`fJzS% O ]> )u+| clear
f2O*8^^Y{Q Y^f 94s:2S close all
eP q13!FC/ -t@y\vZF, n1=1,n2=1.45;
c Pq Dsl3 G`]v_`> theta=0:0.1:90;
Sa;<B:| 11)/] ?/j a=theta*pi/180;
$hjP}- oUX r zM Fof rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Pt6h GSo. aNE9LAms rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
3XeXzPj 4<G? tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
t." g\; m]Fa EQVoE ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
N5 SLF4R1 e 2AN[Ar figure(1)
j+DE|Q&]I cOSxg=~>u subplot(1,2,1);
Iv$:`7|crX ~L %Pz0Gg plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
&W=V%t>Z ;wij}y-6 legend('r_p','r_s','|r_p|','|r_s|')
E?3 0J3 S r0\ cc6 xlabel('\theta_i')
_0'm4?" }>MP{67Dm ylabel('Amplitude')
hLb;5u&!kW B{7Kzwh; title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
]y3pE}R kOs(?= axis([0 90 -1 1])
yicO!:bM )W&o?VRfO grid on
^FP}
qW~;9 J DLT OLG subplot(1,2,2);
$_Y/'IN`k 9[cp7 Rcb plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
{S[I_\3 3B@y &a#& legend('t_p','t_s','|t_p|','|t_s|')
wYZT D*A2h $.H:8^W xlabel('\theta_i')
06 QU "pt+Fe|@c; ylabel('Amplitude')
P!kw;x F-_%>KJS title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
kf<c, 3A S8
:"<B) axis([0 90 0 1])
i 3m3zXt _i5mC,OffN grid on
a% Uw;6|{ ]JOephX2R Rp=abs(rp).^2;
kmryu= HbcOTd)=5 Rs=abs(rs).^2;
!7}IqS s o4$Ott%Wm Rn=(Rp+Rs)/2;
\[:PykS 6SE6AL<b Tp=1-Rp;
<3d;1o vt//)*(.$ Ts=1-Rs;
=WC-Sj{I /QQ8.8=5 Tn=(Tp+Ts)/2;
[+;qWfs B _w 5RK( figure(2)
_@Y"$V]=Vt
CJJD@= subplot(1,2,1);
m9 Ax\lf CdY8 #+"
plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
Mtm/}I }.p<wCPy6 legend('R_p','R_s','R_n')
( ab{F5 Bs[nV}c>> xlabel('\theta_i')
9Gca6e3 /RGNAHtI i ylabel('Amplitude')
F!J J6d53y jk$86ma! title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
zrs<#8!Y_! $(ewk): axis([0 90 0 1])
bp#:UUO%S `-_N@E1'> grid on
p~K9
B-D "
Ya9 ~6 subplot(1,2,2);
h{k_6y m %JtbRs(~q plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
VU|;: v !H:^!z legend('T_p','T_s','T_n')
bLq y!QE (z G.aaz*C xlabel('\theta_i')
$J.T$0pFa IS BV%^la| ylabel('Amplitude')
~]BMrg n Y>B P?l title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
JWROYED {^5?)/< axis([0 90 0 1])
#]9hTa IR - hzj V| grid on
&-%X:~|:X 4 ,G w#@ 
=oX>Ph+ P GMW,*if8p 
MAqLIf<G
