利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
P_y8[ Y]? s.{nxk. 1、光疏射向光密
<q<kqy5s-R MmL)CT clear
x[m&ILr }z|@X KA# close all
S +mM S #CcC& I
:c n1=1,n2=1.45;
i^I
U)\ 84|oqwZO theta=0:0.1:90;
#y2IHO - W6 y-~ a=theta*pi/180;
Kc,=J?Ob g q`S` rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
mu/GOEZ5 dPx{9Y<FzU rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
+T,Yf/^Fn Q"V S;uh.v tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
G Ch]5\ J =j6rD ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Oh]RIWL m R|;}u;d figure(1)
-w3KBlo ZaKT~f%%z subplot(1,2,1);
UFOUkS
F @ 5tW*:s plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
B*c@w~E w1 5Qqh lK legend('r_p','r_s','|r_p|','|r_s|')
-,Y[`(q k&dLg5O xlabel('\theta_i')
q wd7vYBc, Kb icP< ylabel('Amplitude')
?mME^?x
Mu {%! >0@7 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
UU;U,q OMM5p=2Q axis([0 90 -1 1])
eCfy'US;@3 a"Q> K7K grid on
`rQDX<? !8ch&cr)o+ subplot(1,2,2);
]?"1FSu-8r 1v2pPUH\ plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
S9@2-Oc : l&g5 legend('t_p','t_s','|t_p|','|t_s|')
9s9_a4t5 |OarE 2 xlabel('\theta_i')
K H&o`U(} Q}(D^rGP3 ylabel('Amplitude')
C#3K.0a 1:Dm, d; title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
}u9wD08x G1z0q3< B axis([0 90 0 1])
u
^ }R]:n rfwX:R6,g grid on
L4 ) 2z+-vT% Rp=abs(rp).^2;
AhauNS^"{R x+;"(]# Rs=abs(rs).^2;
>/eV4ma" p*l]I *x'< Rn=(Rp+Rs)/2;
?x-:JME0 # qP W J Tp=1-Rp;
iLR^ V! /GUbc Ts=1-Rs;
ckCb) r_ DwBKqhu Tn=(Tp+Ts)/2;
]Ac&h
aAP yD\[`!sWk figure(2)
9 U!-Zn! c*:H6(u subplot(1,2,1);
q(iM=IeiN \
W.u V[\ plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
( Ymj
W^sH|2g legend('R_p','R_s','R_n')
KH_~DZU*5 ^+b ??K xlabel('\theta_i')
jJU9~5i? !y 7SCz
g ylabel('Amplitude')
)cUFb:D*" ^Ox|q_E
w} title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
%"{j NC? 5L% \rH&N axis([0 90 0 1])
a-(O AzQ_ IN#Z(FMVC grid on
.*acw /ltGSl subplot(1,2,2);
J3P )oM[ gI5" \"T{ plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
:^H#i:4 "T0s7LWp legend('T_p','T_s','T_n')
t.YY?5 l !GL
kAV xlabel('\theta_i')
6'YsSde". 1w|C+m/( ylabel('Amplitude')
mO|YX/> hantGw | title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
0'Y'K6hG` 1GA$nFBVC axis([0 90 0 1])
.*_uXQ
<f+ 9wuZ grid on
PW)Gd +y d> OLnG>
F 
<1xs
ya[e L%sskV( 
^ _C]?D?
