利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
xtKU;+# a'. 7)f[g} 1、光疏射向光密
RD jw| V Z:es7<#y clear
] @#w R ^ meU& close all
]qu6/Z )=MK&7 2r n1=1,n2=1.45;
d 2^/ \7 pEn theta=0:0.1:90;
`H$=hr z%iPk'^ a=theta*pi/180;
rm$d v%q lNt xM"G& rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
5h0Hk<N /e*fsQ>M: rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
kqxq'Aq)d iA[o;D# tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
67Qu<9} <- U^% )B I ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
m g@Ol"2 "* N#-=MJF figure(1)
%lGOExV% !$>b}w ' subplot(1,2,1);
:]CL}n$* svb7-.! plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
;Q4,I[?% * fj`+J legend('r_p','r_s','|r_p|','|r_s|')
CQ^(/B^c +"bi]^\z xlabel('\theta_i')
\ym^~ Q| n;$u%2 t2 ylabel('Amplitude')
(
^@i(XQ WVK AA. title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
?%Tx%
dB J YA>Q& axis([0 90 -1 1])
4
2DMmwB $Bl51Vj N grid on
\=mL L|a yjH'< subplot(1,2,2);
r]D U |yzv o"3 plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
#s15AyKz5 Xw< ;)m legend('t_p','t_s','|t_p|','|t_s|')
Dj'?12Onu= &}7R\co3 xlabel('\theta_i')
0GeL">v,:= VBF:MAA ylabel('Amplitude')
fjl 9* ->.9[|lIg title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
#N >66!/V <eS/-W %n6 axis([0 90 0 1])
(_>Su QK J:>o\%sF grid on
d|7LCW+HW Q^nf D
Rp=abs(rp).^2;
i8-Y,&>V v1X[/ \;U Rs=abs(rs).^2;
6
R})KIG CI-za !T Rn=(Rp+Rs)/2;
*AGf'+j*z $@71 w~ y Tp=1-Rp;
9&.md,U ' bHE.E BZ Ts=1-Rs;
rxyeix JJSE@$",\ Tn=(Tp+Ts)/2;
q/ 54=8*h0 -WF((s;<# figure(2)
]4 c+{ r<!nU&FPD: subplot(1,2,1);
' VCuMCV Q1x&Zm1v plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
9X;*GC;d N#"l82^H* legend('R_p','R_s','R_n')
C ~@m6K ,*d8T7T xlabel('\theta_i')
L3xN#W;m7 YW/V}C'> ylabel('Amplitude')
r~ f;g9I PQSmBTs. title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
~M} K]Li UdM2!f axis([0 90 0 1])
at@tS>Dv nQ+5jGP1 grid on
{j]cL !Od JW^ ${4 subplot(1,2,2);
JJ_ Z{ w?|qKO plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
6Z J-oT!. M."/"hV`- legend('T_p','T_s','T_n')
d4\JM 65 )?(Ux1:w) xlabel('\theta_i')
. <z7$lz\ 1v`|mU}i, ylabel('Amplitude')
l7^^Mnk C 5 >0\= title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
z+6PVQ .nrbd#i- axis([0 90 0 1])
0NZ'(qf~9 uK"FopUJ4i grid on
wLqj<ot f)z(9JJL 
\X %FM"r :dULsl$Nz 
Z Se30Rl\
