利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
7n,=`0{r foUB/&Ee 1、光疏射向光密
28qlp>U ;Z); k`j clear
#> 6Jsnv1 0Wj,=9q close all
2Z>8ROv^X _L+j6N.h1 n1=1,n2=1.45;
zx5#eMD (67byO{ theta=0:0.1:90;
X;n09 L`CB ZUkM8M$c a=theta*pi/180;
.N7<bt@~) BA[ uO3\4 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
&\%\ "Zh q@g#DP+C rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
/DjsnU~3 H: S<O%f tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
j
WSgO(y w' E(9 gV ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
>#)%/Ti}DU % o 5'M^U figure(1)
6_:I~TTX 5'( T*" subplot(1,2,1);
`~z[Hj=2 f `D( V-4 plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
k* v ${1& bB>.dC legend('r_p','r_s','|r_p|','|r_s|')
aIDv~#l mfG m>U xlabel('\theta_i')
H32o7]lT {Kf5a
m ylabel('Amplitude')
T B- dV'w pKno~jja title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
ltlo$`PR `5e{ec
c7 axis([0 90 -1 1])
|]9L# uq ;yR[w" grid on
U>OAtiq JX cg o subplot(1,2,2);
8+J>jZ @ meT8S9t plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
:t < S q#@r*hl legend('t_p','t_s','|t_p|','|t_s|')
3H0B+F2XQ no W]E}nN xlabel('\theta_i')
z 1]nC]2 {o 2 qY|S ylabel('Amplitude')
k'k}/Hxub 8*x=Fm,Ok title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
jFJ}sX9] R}cNhZC axis([0 90 0 1])
}Z{FPW.QK ! \$V?*p7 grid on
O@4 J=P=w gO)":!_n W Rp=abs(rp).^2;
e#, (a DIw_"$'At Rs=abs(rs).^2;
lx=tOfj8 s2<[@@@q Rn=(Rp+Rs)/2;
6M. |W; Dk>6PBl Tp=1-Rp;
"l9aBBiu +wJ!zab` Ts=1-Rs;
JSi0-S[Y{ N'WC!K.e Tn=(Tp+Ts)/2;
[ `oVMR <e?Eva%t` figure(2)
8#V D u( S1I.l">P subplot(1,2,1);
hxK;f fBctG~CJH plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
n=bdV(?4 r uGeN legend('R_p','R_s','R_n')
* m^\& D4,>g )B xlabel('\theta_i')
yk)]aqic `:;q4zi j; ylabel('Amplitude')
o.Rv<a5.L ]/H6%"CTa title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
gK9d `5 5K 2K'ZkI axis([0 90 0 1])
tV`=o$` ^a_a%ws grid on
*; ]}`r g+4y^x(X@1 subplot(1,2,2);
9i}$245lB Pv/ v=s>X plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
giX[2`^NG | Ia9bg'1U legend('T_p','T_s','T_n')
|Rz.Pt6 vM6W64S xlabel('\theta_i')
cl2_"O M@{#y EP ylabel('Amplitude')
N
UX | U]~@_j title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
]5c(:T F >#x[qX axis([0 90 0 1])
#aU!f"SS T1 $E][@Iv grid on
+q '1P}e (pd $?vRy 
LCH\;07V# l^9gFp~I 
~}}<+ JEEO
