利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
}ozlED` E l[lUmE 1、光疏射向光密
bg;N B oZd 'b#RfF,7H} clear
s'' ?:
+ //cj$}Rn! close all
.r[b!o^VR e \x=4i n1=1,n2=1.45;
w6DK&@w`'/ fmZ5rmw! theta=0:0.1:90;
wr{03mQHxp d!kiWmw, a=theta*pi/180;
&}wr N(?w hV|pH)Nu{ rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
# TZf\0\! nD6mLNi%a rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
0Ci \( g.OBh_j-v tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
#R$d6N[H Qb}7lm{r ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
le*mr0a W$LaXytmak figure(1)
C c*"cQ e CqX2R:# subplot(1,2,1);
-BUxQ8/, *n mr4Q'v{ plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
$G/ h-6+8 N1$lG?
)+ legend('r_p','r_s','|r_p|','|r_s|')
jqxeO N WmU4~. xlabel('\theta_i')
}&#R-eQT V\^rs41$; ylabel('Amplitude')
w"d~R XH4!|wz title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
>B (%$jG Z , f9V`Pz) axis([0 90 -1 1])
#g[jwl' &No6k~T0:b grid on
Poy ]5 :. >m='#x0>Y subplot(1,2,2);
Sx)b~ * =H6"\`W plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
jqq96hP, tWR>I$O8F legend('t_p','t_s','|t_p|','|t_s|')
)\!_ `ob 'Lu7cb^ xlabel('\theta_i')
R M\it"g 0,+RF "R ylabel('Amplitude')
V5sH:A7GJ h|OqM:J; title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
P#-9{T pQ%~u3 axis([0 90 0 1])
q[p +OpA ; okFm grid on
*sK")Q4N 8 tMfh Rp=abs(rp).^2;
am.}2 QZU x 95[*[ Rs=abs(rs).^2;
{~NiGH Y = Q|_v} Rn=(Rp+Rs)/2;
o C0K!{R*
L,#ij!txS Tp=1-Rp;
0'y9HE'e 2w)[1s[ Ts=1-Rs;
`^HAWo;J ,] HH%/h
Tn=(Tp+Ts)/2;
U-/-aNJ]U l>K z5re^ figure(2)
1{@f:~ v? z5G<h subplot(1,2,1);
w$I<WS{J:Z q\wT[W31@ plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
EIZSV> q#9JJWSs legend('R_p','R_s','R_n')
"@: b' m ,{wA%Oy, xlabel('\theta_i')
&?L
K>QV d]Y-^&]{] ylabel('Amplitude')
oc .H}Eb%Z mlC_E)Ed5 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
fNb2>1 P.^ %8L axis([0 90 0 1])
<Stfqa6FJ ,i KEIxA! grid on
p)l >bC?3 1 .\|,$ subplot(1,2,2);
t\LAotTF/ _ SU%ul plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
^8 dd JU^ {!u legend('T_p','T_s','T_n')
>C,=el M qK}4r5U xlabel('\theta_i')
b0W~*s [4 +$ Q.N{LV ylabel('Amplitude')
Y rq-( & h)G>Sqc title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
')C %CAYW Y r3h=XY axis([0 90 0 1])
W
vh3Y,|3 Gvg)@VNr grid on
EB8=* B8 3I0=^ >A 
paFiu Q D)C^'/8q 
5 9$B
z'LY
