利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
B !)9
> @Bkg< 1、光疏射向光密
2'W# x L)7{_s clear
oIj/ V|ByK /u]#dX5 close all
Kt]vTn7!9 P $ h) Y n1=1,n2=1.45;
438r]f?0|{ 9 '|k@i: theta=0:0.1:90;
I}y6ke! xo
^|d3 a=theta*pi/180;
dW5r]D[Cx 7}c[GC )F rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
86qQ "=v "[z/\l8 O rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
^ -~=U^2tC Ha ZV 7 tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
WyQ8}]1 b "=40% j0 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
Ei<+{P(t0 dR,a0+! figure(1)
qOyS8tA.H eo!+UFZbY subplot(1,2,1);
"J}B
lB rfgsas{F plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
TOq xl ,@N.v?p> legend('r_p','r_s','|r_p|','|r_s|')
~Lu,jLKL=[ $'#} f? xlabel('\theta_i')
hD,^mru Rdd9JJsVd ylabel('Amplitude')
/@&uaw ,gkWks l9 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
">fRM= fl &
}"I! axis([0 90 -1 1])
\
sz ](X I;$tBgOWq grid on
!HXsxNe ! ([ v=O# subplot(1,2,2);
QqeF )J[Ady^5 plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
K_N` My 4l/~::y legend('t_p','t_s','|t_p|','|t_s|')
&XhxkN$8 VWCC(YRU|$ xlabel('\theta_i')
>u
.u#d e $ 'a]lR ylabel('Amplitude')
^il gd A>:31C title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
~ :B/`1[m IeIv k 55 axis([0 90 0 1])
d;r,?/C 7:.!R^5H grid on
oJ=u
pnBn- U!nNT== Rp=abs(rp).^2;
|ew:}e: k< t$Ua&w Rs=abs(rs).^2;
:3}K$ =D 5!Xq'| Rn=(Rp+Rs)/2;
<2+FE/3L iW? NxP Tp=1-Rp;
3ZYrNul" ljh,%#95= Ts=1-Rs;
-]Oi/i, { W5 RZsS] Tn=(Tp+Ts)/2;
q1%xk =8 OadGwa\:s figure(2)
-]MZP:s +$G P(Uu, subplot(1,2,1);
j0e1CSE xSjs+Y;Mu plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
j 2J ew o6:p2W legend('R_p','R_s','R_n')
~eGtoEY VDv>I 2% xlabel('\theta_i')
!=+ hU/e b=Q%Jxz? ylabel('Amplitude')
` XbV*{7 EPeV1$ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
ns8I_H F BI^}^#_ axis([0 90 0 1])
^eqq|(<K wQ9fPOm grid on
]<E\J+5K t*!Q9GC_ subplot(1,2,2);
Gh{vExH@5( %X Wb|-= plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
Y<` uq 'V ]3yaIlpD1 legend('T_p','T_s','T_n')
[ Q20c <, c< g{ &YJ xlabel('\theta_i')
0o&7l %Y/ ?|we.{ ylabel('Amplitude')
A j2yAg lV<j?I~?Q title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
,O"zz7 ;jpsH?3g axis([0 90 0 1])
jQ?6I1o nSV
OS6 grid on
[,p[%Dza Q W}N,j$ 
>I.X]<jI '${xZrzmt 
l8Zz Kb-
