利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
VQ^}f/A @3Gr2/a 1、光疏射向光密
Zce /& u_jhmKr~ clear
~ ^*;#[< +{au$v} close all
# b94S?dq J4#rOS n1=1,n2=1.45;
fzjA P7 y B3'-: theta=0:0.1:90;
Cnv?0to2l "g:1br?X,9 a=theta*pi/180;
?<STl-]& qaiR329fx rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
R?i C"s! jW;g{5X rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
)<[)7` A8T8 +M: tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
4KB>O)YNg' raJv$P ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
'J_`CS bP VQ- figure(1)
5F$~ZDu >!W H%J subplot(1,2,1);
OQiyAyX ): 7mK03J plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
x &*2R#Ai x}; sti R legend('r_p','r_s','|r_p|','|r_s|')
h?P-
:E _*fOn@Vwo xlabel('\theta_i')
JVR,Py:%G V,&A ?
Y ylabel('Amplitude')
y& 6 pc D\^\_r): title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
s w+vyBV)r #AF.1;(k axis([0 90 -1 1])
)8%m|v#W ~PyZh5 x grid on
]h%~'8g, 0SV4p . subplot(1,2,2);
7x#."6>Dy ).IK[5Q ` plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
zI" 1.^Trn T\w{&3ONm legend('t_p','t_s','|t_p|','|t_s|')
QC6 :ZxP K)W:@,* xlabel('\theta_i')
d~#:t~
$, vP'#x ylabel('Amplitude')
R\-]t{t` `J %35 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
zKfY0A R LC})aV| axis([0 90 0 1])
3jjV
bm }U@(S>,% grid on
t5 k=ngA [r^WS;9n Rp=abs(rp).^2;
={\9-JJhE }#2I/dn Rs=abs(rs).^2;
rA /T>ZM $S/EIN c Rn=(Rp+Rs)/2;
RMlx[n sq .* &F Tp=1-Rp;
|O{kv}Y Z 3|BB#; Ts=1-Rs;
(BGflb *g "X hk Tn=(Tp+Ts)/2;
soh9Oedml- cUr5x8<W). figure(2)
Lum5Va%0 #6@4c5{2=4 subplot(1,2,1);
4o<'
fY W1ql[ DqE{ plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
[D2< )
xM$AhH legend('R_p','R_s','R_n')
w%$J<Z^-? BBa!l e9P xlabel('\theta_i')
d "25e"(~F Bo0f`EC I ylabel('Amplitude')
6*:U1{Gl) vF9 fXY= title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Qn}M 814cCrr,o axis([0 90 0 1])
#9t3 <H[ GYtp%<<9; grid on
EzU3'x C);I[H4Yfw subplot(1,2,2);
fvRqt)K s |xrnLdng0R plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
iN 1_ T ''}2JJU{ legend('T_p','T_s','T_n')
, 8o
Y(h "] ]aF1 xlabel('\theta_i')
o\X|\nUk d
9]zB-A ylabel('Amplitude')
>aT~ G!y hbm #H7Y title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
nnCz!:9p ?~ <NyJHN% axis([0 90 0 1])
u<EPK*O* 6S^JmYq grid on
=m 6<H Zou;o9Ww 
Z<#h$XUA D&1(qi=x& 
=lL)g"x X
