利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
V_x8
Q+~? 3J3Y t` 1、光疏射向光密
Wx:He8N] H 6E)emFkQ clear
Bn4wr ?@>PKUv{ close all
j;7:aM"BQW +u[^@>_I0 n1=1,n2=1.45;
]jB`"to*} ]B2%\}c theta=0:0.1:90;
vW s#4JoG |7$Q'3V a=theta*pi/180;
qexnsL : Yb_ rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
w$evAPuz^ b_&KL_vo{| rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Fd1 t/B, {0o ,2]o!: tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
o..iT:f;n L=_ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
/S #Z.T~~ w;wgh`ur figure(1)
` 'vNHY h N U.y subplot(1,2,1);
.gJv})Vi 4 }r.g0L plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
/* G-\| o'KBe%@/ legend('r_p','r_s','|r_p|','|r_s|')
J?$`Tnx^ ul& }'jBr xlabel('\theta_i')
)F\^-laMuK {5NE jUu{j ylabel('Amplitude')
Q>yO,H| .5E6 MF title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
*km!<L7Y KX^! t3l6 axis([0 90 -1 1])
cl{kCSZo.z < ht^Ck grid on
5{|tE! E5X#9;U8E" subplot(1,2,2);
#PGpB5vnaA ?~9o2[ plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
AT - Fp(-&,L0fc legend('t_p','t_s','|t_p|','|t_s|')
l|S_10x5 ,ieew` xlabel('\theta_i')
TNF+yj-|X: (v0i]1ly[ ylabel('Amplitude')
y;VmA#k` w?JM;'<AYQ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
y6P-:f/& * (J5M+K\H axis([0 90 0 1])
U}<' [o
V KP xf grid on
QD; :!$Du Reb o.6rG Rp=abs(rp).^2;
v m.%)F#@ Z5 Tu*u= Rs=abs(rs).^2;
bcUa'ZfN< fD#!0^ Rn=(Rp+Rs)/2;
\y0]BH o9+fA H`D Tp=1-Rp;
'p&q}IO `OHdo$Y9 Ts=1-Rs;
:+v4,=fHy :OD-L)Or Tn=(Tp+Ts)/2;
=kBWY9 :$, b" kL)DL1L figure(2)
2!nz>K =GL^tAUJ subplot(1,2,1);
n*m "yp ^TqR0a-* plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
0O|l7mCr%I 4<S=KFT_ legend('R_p','R_s','R_n')
uC5W1LyI $p\ 0/ xlabel('\theta_i')
N 8pzs" X8 x:/]/0 ylabel('Amplitude')
**~1`_7~* \^Z DH title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
8.tp#x,A *|S{%z9> axis([0 90 0 1])
`e }6/~R` Wo,fH Y grid on
<]u]rZc$ B18?)LA subplot(1,2,2);
Df}3^J ~JX S<Uv/pn plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
tREC)+*\ r~; TId} # legend('T_p','T_s','T_n')
9@ 8)ZHf ?dQ#%06mn xlabel('\theta_i')
|K'7BK_^J s%@HchZ 1 ylabel('Amplitude')
7DB!s@"
X~rHNRIU title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
PaBqv] F=V_ ACU axis([0 90 0 1])
m8z414o [Ow rIL grid on
T#=&oy7 `YK% I8 
%s#`Z [8, 3NlG,e'T2 
-3Auo 0
