利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
\,YF['Qq X.u&4 SH 1、光疏射向光密
o!q3+Pp;} W\JbX<mQ clear
|@{4zoP_N ^_6.*Mvx close all
#S?xR qkc jV 83%%e n1=1,n2=1.45;
E$ B7E@(U "E4;M / theta=0:0.1:90;
aS
]bTYJ' zGtWyXP a=theta*pi/180;
nVkPYeeT q J)[2:.G rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
~s#vP<QHa 6Fp} U rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Q8] lz} >LvQ&fAo tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
74Fv9 uB;\nj5'D ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
<1r#hFU UL {bQi
z figure(1)
${ (c `X xvx5@lx subplot(1,2,1);
d D%Sbb "pQFIV, plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
qa>Z?/w 6N7^`ghTf legend('r_p','r_s','|r_p|','|r_s|')
An cka ii< /!B( xlabel('\theta_i')
-&L(0?*qo 0l.\KF ylabel('Amplitude')
kU*Fif dRl*rP/ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
|wef [|@% wrORyj axis([0 90 -1 1])
~353x%e' 33O O%rWi grid on
ssf.ef$ !CJh6X ! subplot(1,2,2);
S6Er# )k @m#1[n; plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
FLW Q Y, U|[+ M@F_L legend('t_p','t_s','|t_p|','|t_s|')
mpsi{%gA
?^y! }( xlabel('\theta_i')
V:<NQd } "QV{W ylabel('Amplitude')
3Ryae/Nk AeQIsrAHE title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
vdhwFp~Y 8`I/\8;H'p axis([0 90 0 1])
E~[v.3` 0uw3[,I
grid on
"Uk " erhxZ|."P Rp=abs(rp).^2;
-#TF&- Cob<N'. Rs=abs(rs).^2;
g8+Ke'=_ s":\ > Rn=(Rp+Rs)/2;
N:4oVi@Je EZ N38T Tp=1-Rp;
c8R#=^ DD 4hymQ3
g Ts=1-Rs;
oU\Q|mN( >v r! 3 Tn=(Tp+Ts)/2;
|._9;T-Yde 0=WZ 8|R figure(2)
%.NOQ<@W ;usv/8 subplot(1,2,1);
5.]eF$x2 Kuj*U'ed7t plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
-GqMis}c Q&JnF`* legend('R_p','R_s','R_n')
nQ\)~MKd NWN Pq" xlabel('\theta_i')
6Cdc?#& xJ0Q8A ylabel('Amplitude')
1@i 8ASL &efwfnG< title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
{}Q A#:V 7H,p/G?]k axis([0 90 0 1])
.q!U@}k. o3s ME2 grid on
}@ +{;" JQ[~N- subplot(1,2,2);
CO)BF%?B lZ\8$,B) plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
ox}LC, ! 4m qA*c%6S legend('T_p','T_s','T_n')
p!XB\%sv'" /q3 ]AVV xlabel('\theta_i')
qi ;X_\v 96 ozt UK ylabel('Amplitude')
*irYSTA$ N~Kl{" >` title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
t9 Sog~:' z }t{bm axis([0 90 0 1])
L_k'r\L <<ze84 E grid on
GadZ!_.f -0tHc=\u( 
Zy <gA > ne~=^IRB 
BB>R=kt
