利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
_ 9yW; i- #f@}$@ 1、光疏射向光密
m<sCRWa- } _=h]|6t clear
G5hf m- ZZ>F ^t close all
$MqEM~^= os n ,kD* n1=1,n2=1.45;
wEZi eHw vM?,#:5 theta=0:0.1:90;
mWF\h>]|. O{x -9p a=theta*pi/180;
CC)Mws+2 W[f%m0 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
L 8J] X7 ; GEr8_7 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
4k!>JQor s.^+y7$ tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
%$_Y"82 !5
?<QKOe ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
T_[5 ZYy 4C[kj figure(1)
@6 M>x=n5 ]?T,J+S subplot(1,2,1);
{A~3/M%74; %GP`H/H( plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
1HLU
& d ATAH}r& legend('r_p','r_s','|r_p|','|r_s|')
L |8&9F\ 1 F&}e&}c xlabel('\theta_i')
dVMLn4[,MA $9h^tP'CV ylabel('Amplitude')
=zR9^k @~&|BvK% \ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
<qVOd.9c HQSFl=Q axis([0 90 -1 1])
sf
O{.#5< p+228K ;H grid on
`_sKR,LhB F-XMy>9 subplot(1,2,2);
?69E _E n\< uT1n plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
cT@H49#uB )Y`ybAD d3 legend('t_p','t_s','|t_p|','|t_s|')
eM]>" y Ni3@f xlabel('\theta_i')
9VMk ? 9rIv-&7'm ylabel('Amplitude')
J\FLIw4 dGUiMix{N title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
nchpD@'t 0EasPbp axis([0 90 0 1])
[ 'Qh#^p +Ht(_+To1 grid on
']d(m? b8 1cq, Rp=abs(rp).^2;
ix 38|G9U h7 mk< Rs=abs(rs).^2;
zCv)%y KpIY>k Rn=(Rp+Rs)/2;
|"[;0)dw^ (w`_{%T Tp=1-Rp;
nkz<t 1d$wP$ Ts=1-Rs;
?Bb EQr t~$8s G\ Tn=(Tp+Ts)/2;
3BA Q2S} '$VP\Gj. figure(2)
2k m0 ) NTpb subplot(1,2,1);
D&=+PAX 2Ima15^+ F plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
E&Zt<pRf;2 v?}rA %so legend('R_p','R_s','R_n')
i ~)V>x <tm= xlabel('\theta_i')
b ! `6s k8i0`VY5Y ylabel('Amplitude')
[;l;kom rCb$^(w{7 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Jx Kd ~fs}
J axis([0 90 0 1])
PP/#Z ~.M qxcTY|& grid on
BqpJvRJd +U>Y.Y P subplot(1,2,2);
i>C%[dk9 W e*uZ? + plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
lv~ga2 >z =$ T[ legend('T_p','T_s','T_n')
@:@5BCs< )TBm?VMe xlabel('\theta_i')
xZ* B}O{{H > "z`))9 ylabel('Amplitude')
"`mG_qHI[ yTN HM_P title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
wod(P73? LMN`<R(q] axis([0 90 0 1])
- FrK'!\ of > grid on
o)Px d rl2(DA{ 
C0khG9,BL ScQ9p379 
|oQhtk8.
