利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
L% ?3VW _=;lt O 1、光疏射向光密
m kHcGB!~ o!ycVY$yW clear
Z@b GLS N"rZK/@} close all
{Qr0pjE7R `<^1Ik[g n1=1,n2=1.45;
{kp^@ 3&z.m / theta=0:0.1:90;
K5SP8<. rYw UD7ip a=theta*pi/180;
5D%gDw+" 6 %RN- rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
Y;6%pm $ #Bgq]6G2 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
hpO Uz% hCKx%&[^7 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
hXx. ' 5%`[& ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
W:z!fh- cPm- )/E)i figure(1)
njN]0l{p # cD20t subplot(1,2,1);
fK{[=xMr@ Y94/tjt plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
!@vM@Z" nfbq J legend('r_p','r_s','|r_p|','|r_s|')
@"E{gM@B xMAb=87_
xlabel('\theta_i')
e=%6\&q Lb}
c jI: ylabel('Amplitude')
K0Zq )< 5zBayJh# title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
:u$+lq Lc:DJA axis([0 90 -1 1])
.A2u7*h& J \ V.J/ grid on
9M$N>[og t[- 0/-4 subplot(1,2,2);
,@'M'S 1 KtPq, plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
p~xrl jP$ !(wH}ti legend('t_p','t_s','|t_p|','|t_s|')
]
336Fg T P.;S6i
n xlabel('\theta_i')
@\UoZv( \1p5$0z ylabel('Amplitude')
;y{ VdT VF<C#I title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
\O7Vo<B&D r\-25F<e5 axis([0 90 0 1])
* (4TasQu k-LT'>CWl grid on
+.b@rU6H %V CfcM}5I Rp=abs(rp).^2;
&MlBp I < ?I s ~[2 Rs=abs(rs).^2;
n*oa J<o% %$l^C!qcY Rn=(Rp+Rs)/2;
E(&GZ QE oe5.tkc Tp=1-Rp;
f @c s<x C;j V{sb9c Ts=1-Rs;
dTV:/QM 8zRb )B+ Tn=(Tp+Ts)/2;
Q:x:k+O- h!*++Y?&0 figure(2)
zTa>MzH1-; ?"aj&,q+ subplot(1,2,1);
r12{XW?~ |;-,(509 plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
7Pc0|Z/ Cl H aR legend('R_p','R_s','R_n')
g1uqsqYt ] _/d xlabel('\theta_i')
B9 Dh^9?L 7 \
<4LX ylabel('Amplitude')
1x0 7ua@(v x,Cc$C~YP title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
UEq;}4Bo PSdH9ea axis([0 90 0 1])
4nhe *ip !k&~|_$0@ grid on
8dw]i1t< FNDLqf!j subplot(1,2,2);
MGO.dRy_ _e.b #{=9 plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
~EU[? tH:K6^oR legend('T_p','T_s','T_n')
8i;N|:WdH ?;,Al`/^ xlabel('\theta_i')
6Upg\( { /BEO=8q2 ylabel('Amplitude')
bH 41#B Gc9^Z= title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
7[-jr;v wS2iyrIB axis([0 90 0 1])
Ss u{Lj g[;iVX^1& grid on
ru`;cXa, 34C
^vBp 
F.=u Jdl.! Q ;P ~' 
D ^PsV
