利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
fj&i63 ?e VJS8)oI~ 1、光疏射向光密
2@
Z(P.Gh 'PP #^aI, clear
e#k<d-sf6 x|vqNZ\F close all
O]m+u Ox@P6|m n1=1,n2=1.45;
e zOj+vz ) l#E}Uz theta=0:0.1:90;
e$pMsw 'MJ <wAFy>7 a=theta*pi/180;
Hik3wPnp _l1 NKk rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
:GJ &_YHf &@3H%DP}Ql rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
x)_0OR2lkp $r3i2N-I tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
1Pa tH[T[ y
TDNNK ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
ERN>don2 +k>.Q0n%m figure(1)
;{ [.Zu &%/7E_j7 subplot(1,2,1);
b?'yA Xk +U3m#Y )k plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
NG6& :4! Q6r7.pk"SU legend('r_p','r_s','|r_p|','|r_s|')
RG4 sQ0 ]H1I,`=@ xlabel('\theta_i')
;"m ,:5% >_LDMs[-p ylabel('Amplitude')
qQH]`#P 5;{H &O9Q title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
$ O5UyKI wLH] <k axis([0 90 -1 1])
D~\$~&_]= y? co| grid on
YGsWu7dG L{`J Ru subplot(1,2,2);
\(vY%DL1: Dmu/RD5X: plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
AoI /n4T^ mOlI#5H legend('t_p','t_s','|t_p|','|t_s|')
}_68j8` LoGVwRmoC xlabel('\theta_i')
r#B+(X7LM 9\DQ>V TQ ylabel('Amplitude')
TU
1I} , 'uxX5k/D@t title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
W!&vul5 O7 $hYk axis([0 90 0 1])
'u1?tQ=gmk 8M `#pN^ grid on
Z LD}a:s v7`HQvQEz= Rp=abs(rp).^2;
RSfzRnhmr
7?%k7f Rs=abs(rs).^2;
|1>*;\o- GHeVp/u Rn=(Rp+Rs)/2;
G-.^O,% ]V#M%0:Q82 Tp=1-Rp;
[ n$BRk| heK7pH7;d Ts=1-Rs;
)6J9J+%bi ?28N ^ Tn=(Tp+Ts)/2;
2j[&=R/. UTH_^HAN#G figure(2)
k4 [|'Dk? *NCl fk Z subplot(1,2,1);
'4#}e[e hBX!iukT|{ plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
lh;:M -b9 < "r#:Wr legend('R_p','R_s','R_n')
"PpN0Rr B,=H@[Fj xlabel('\theta_i')
Ch3jxgQY k|H: ylabel('Amplitude')
Yn G_m] :d#VE-e title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
la+[bm< v ]3 l 9:| axis([0 90 0 1])
sU{NHC)5 wLOB}ZMT grid on
sH[ROm e F3,2DD C subplot(1,2,2);
-u8NF_{c ssN6M./6 plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
@0u~ ?!g@ x-?Sn' m legend('T_p','T_s','T_n')
[6XF=L,! f V'ZsJ N xlabel('\theta_i')
F ;2w1S^ I\$X/t +dH ylabel('Amplitude')
#odI EC/ A* /Hj TX title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
j+,d^! "*srx] axis([0 90 0 1])
aD$v2)RR j4R(B grid on
4QIX19{" T"DlT/\ 
<,e+
kL{ Qh'AT o 
`6]%P(#a
