利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
5h6c W M8HHyV[AmC 1、光疏射向光密
O${B)C, +#a_Y clear
o+j~~P TCI%Ox| a close all
RC>79e/u< # 2qDn^s n1=1,n2=1.45;
Kk?C x'IYWo
] theta=0:0.1:90;
c&AJFED]< 6mAaFDI,R a=theta*pi/180;
g4-UBDtYt [x\?._> rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
y (w&6: #'&&&_Hu3 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
rE[:j2HF VQ,;~^Td tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
U.,_zEbx, ljw(cUM ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
(o5^@aDr
^D.u figure(1)
`: R7j f 8'}D/4MUr subplot(1,2,1);
aIXN wnq MJDW-KL- plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
.Fe_Z)i> h f0 d*% legend('r_p','r_s','|r_p|','|r_s|')
EO:
VH 2:4:Q[{A xlabel('\theta_i')
UcI;(Va (0W)Jd[ ylabel('Amplitude')
7)Y0D@wg ~As/cd>9 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
W#@Mx I "5VkeIx axis([0 90 -1 1])
cFF'ygJ/ jMV9r-{*+ grid on
lC AD $Ia~ lK 0pr subplot(1,2,2);
lI*uF~ 'D c ?(X (FQ plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
bnYd1 9> |h6 u%t2AY legend('t_p','t_s','|t_p|','|t_s|')
BdQ/k XZu+ % r>v^1Vo xlabel('\theta_i')
alV{| Vf[6 Fk*C 8 ylabel('Amplitude')
&]_2tN=S$ _Q=h3(ZI title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
n=8DC& px>g axis([0 90 0 1])
&o]ic(74c? qQ
T ^d grid on
Fd(o8z8Q 8`GN8 F Rp=abs(rp).^2;
"t!_b ma IL&;2% Rs=abs(rs).^2;
o_@6R"| =_:Mx'7 Rn=(Rp+Rs)/2;
> %B7/l$ Y4j%K~ls Y Tp=1-Rp;
aP}30E*Y ,GTIpPj Ts=1-Rs;
L2}p<?f dZIruZ)x Tn=(Tp+Ts)/2;
l/6$BP U` dc=~EG-_rM figure(2)
3U!#rz" 1 h,m subplot(1,2,1);
iQ#dWxw4 E~K5n2CI plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
z 5q( <C{uodFll legend('R_p','R_s','R_n')
mF`%Z~}b %Xjg/5G - xlabel('\theta_i')
W %_Cda5, Q35jJQ$<` ylabel('Amplitude')
Wa !}$q+ t1rAS.z& title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
^nGKuW7\ s.j6"
Q[W axis([0 90 0 1])
~;D5j ) 9I ,%4~ulKMn grid on
:v o#( hreG5g9{ subplot(1,2,2);
Ds@ nuQ M'>8P6O plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
_<m yM2z U {?#W legend('T_p','T_s','T_n')
RM6*c
. /&!4oBna xlabel('\theta_i')
K1_#Jhz %z "${ zw ylabel('Amplitude')
K!jMW lSK<LytB title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
(>M?
iB | |^+( axis([0 90 0 1])
>(BAIjF
E\ TJ+,G4z grid on
FQqk+P! `Ti?hQm/ 
.X;D I<K /?_5!3K J 
. g8db d
