利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
l[ P VWM !gKz=-C 1、光疏射向光密
c2,;t)%@ E K*]^0 clear
\H -,^[G3 l )m]<E X close all
RaBq@r*( MB:VACCr n1=1,n2=1.45;
VOY#Y*)g t'.oty= theta=0:0.1:90;
[JzOsi~R 7F;dLd' a=theta*pi/180;
c'XvZNf .C M<
1rQW' rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
M;$LB@h rZt7C(FM$7 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
y. A]un1 V|njgcn d tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
="<S1}. s& gzv =v ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
1vG]-T3VC hHoc>S6^M figure(1)
]%4rL
S
%dN', subplot(1,2,1);
<8SRt-Cr DJu&l plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
&AS<2h B .Y^U Pxf@ legend('r_p','r_s','|r_p|','|r_s|')
U&\2\z3{ VZRM=;V xlabel('\theta_i')
f5droys9 4NN81~v 4 ylabel('Amplitude')
}z/Y
Hv% `"#0 \Wh title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
\B$Q%\- PX ? # cX_ axis([0 90 -1 1])
j. m(Z} \!_:<"nX. grid on
j+E[ [
OG C|elSM subplot(1,2,2);
_s&sA2r<
"?yu^ plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
c/g"/ICs ^&|$&7
legend('t_p','t_s','|t_p|','|t_s|')
v.aSf`K o (OC3 xlabel('\theta_i')
fV` R7m. 0aF&5Lk`y ylabel('Amplitude')
L,QAE)S'a gkN|3^ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
BU|bo") NVF gRJ& axis([0 90 0 1])
OSvv\3= #$qhxYyd grid on
/wj L< <^Hh5kfS' Rp=abs(rp).^2;
Q>FuNdUk ]('isq,P Rs=abs(rs).^2;
?2g\y@ n-cz xq%n Rn=(Rp+Rs)/2;
kdWk{ZT^ <.d0GD`^ Tp=1-Rp;
T]5JsrT $X %w9l e Ts=1-Rs;
NINaOs }
r#by%P Tn=(Tp+Ts)/2;
0 w"&9+kV XG5mfKMt+ figure(2)
C?FUc cI 8bB'[gJ]{ subplot(1,2,1);
~D4%7U"dv T/V 5pYl plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
_Z:WgO]. r9/PmZo4x legend('R_p','R_s','R_n')
Xm7Nr# K4vl#*qn xlabel('\theta_i')
X 9f!F2x N'fE^jqU ylabel('Amplitude')
r lalr+Rf ,qgph^C title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Ld[zOx ;[-dth axis([0 90 0 1])
7W[}7Y qbjLTE= grid on
p ; dH[NW BzqM$F(
L, subplot(1,2,2);
( Ck |RojC /-wAy-W plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
JkazB1h j*m7&wOE legend('T_p','T_s','T_n')
aClXg- &46h!g W xlabel('\theta_i')
-{i;!XE$SR x4cP% {n ylabel('Amplitude')
KZ/U2.{O< 5!r?U title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
cn$E?&- m}] bP axis([0 90 0 1])
+B$ o 8V 9 ve q grid on
0^4uZeW? < @9p|[! 
LlP_`fA TY[1jW~{r 
XSkN9LqZ
