利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
ub^h&= \S K1mPr^3rC 1、光疏射向光密
}zkHJxZgE Tl(^ clear
}\tdcTMgS QdT}wkX close all
=mS\i6 63 $?Yry. 2 n1=1,n2=1.45;
l~c@^! {@6= Q 6L theta=0:0.1:90;
:o0JY= 5 =XA;[PVx:# a=theta*pi/180;
_ "?.! HW3 }uP\c rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
3h;{!|-3 EYtL_hNp}I rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
7C,&*Ax,9 7mI:| G tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
sD{Wxv B:5Rr}eY+ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
K&vF0*gN3 ah+~y,Gl figure(1)
>zw.GwN | U{7w#>V
. subplot(1,2,1);
]$ L| f!\l g plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
}YB*]<] {@eJtF+2 legend('r_p','r_s','|r_p|','|r_s|')
{IxA)v-` Z, sv9{4r xlabel('\theta_i')
7E!IF>` k.5 u ylabel('Amplitude')
OV m \ v25R_""~ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
iP
=V8g?L ock Te 5U axis([0 90 -1 1])
N >FKy'.gk ]JCvyz
H
grid on
1pt% Kw*@j kDEPs$^ subplot(1,2,2);
I;e=0!9U PH1p2Je plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
fKeT,U`W 0 t Fkd legend('t_p','t_s','|t_p|','|t_s|')
}p}[j t aoT M xlabel('\theta_i')
-&NN51-d\j u?4:H=;> ylabel('Amplitude')
TT2d81I3m "U e. @> title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
1Yq?X: _y),J'W^3u axis([0 90 0 1])
uD}2<$PP cv?06x{ grid on
/M|2 62% 0+?7EL~ Rp=abs(rp).^2;
k[bD\' 4i6q{BeHn Rs=abs(rs).^2;
)k\H@Dy%$ kw#;w=\>R{ Rn=(Rp+Rs)/2;
Wl B 69 5V3R 7 Tp=1-Rp;
G'oG< /A ~ DBcIy? Ts=1-Rs;
ZjgsR|i lt
^GvWg Tn=(Tp+Ts)/2;
$CY~5A `l9 mR&H9 NG figure(2)
v>$'iT~ l j"}*T subplot(1,2,1);
[Yo,*,y31 9X j7~, plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
RZHd9v$ N9jH\0nG legend('R_p','R_s','R_n')
T;L>;E>B x,rlrxI xlabel('\theta_i')
'_GrD>P)- wj,:"ESb4 ylabel('Amplitude')
>d,jKlh^.% Fgp]l2* title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
"/3 db[ U* i{5/$ axis([0 90 0 1])
R$ q;
! C"!gZ8*\!9 grid on
]!I7Y.w6 C1 jHz subplot(1,2,2);
=osv3>&q -n=^U plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
z`!Xh U nSW=LjrO~< legend('T_p','T_s','T_n')
.$qnZWcgG b\3Oyp> xlabel('\theta_i')
$"FQj4%d I(V!Mv8j ylabel('Amplitude')
]~aF2LJ_q L;*ljZ^c title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
P0W*C6&71| ^su<uG<R axis([0 90 0 1])
JBXrFC; l5zS grid on
v*r7Zz6l =T(6#" 
M-1ngI0H; 7s/u(~d) 
XZ sz/#
