利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
u>TZt]h8 ))&;}2{ 1、光疏射向光密
A j2OkD d/ bE t& clear
N 5i+3& $OP7l>KZY close all
kAu+zX>S+ L
]w/P| n1=1,n2=1.45;
sFQ^2PwbS _cc#Qlw 7 theta=0:0.1:90;
B<~ NS)w L+C y Qq a=theta*pi/180;
-u? S=h} <
V\I ~; rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
,:??P1 xA9:*>+ > rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
N_pJk2E Fd#m<" tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
]|sAK%/ T' LIrf ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
D5:{fWVsV/
v76D3'8 figure(1)
CE"JS-S? [-)r5Dsdq subplot(1,2,1);
e'=MQ,EWd .N*Pl(<[ plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
2 ho>eRX (sQr X{~ legend('r_p','r_s','|r_p|','|r_s|')
.VkbYK fm^@i;D
xlabel('\theta_i')
#vzEu
)Ul /__PSK ylabel('Amplitude')
l983vKr HU9Sl*/ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
$DZH QH HmZ{L +" axis([0 90 -1 1])
P) uDLFp] r%e KFS grid on
.S%0 YFTjPBV subplot(1,2,2);
rvyr xw%[ 0v]?6wX plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
Gdg )9 t#S<iBAZ legend('t_p','t_s','|t_p|','|t_s|')
8kbY+W%n [y0O{,lI xlabel('\theta_i')
IJJ%$%F/ q| %(3,)ig ylabel('Amplitude')
(4@lKKiU%H gh>>Ibf title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
[lk'x zE }'OHE(s axis([0 90 0 1])
Nk[2nyeO> a"av#Y grid on
=E&OuX-R 0A:n0[V:] Rp=abs(rp).^2;
kqb0>rYa bT6sb#"W Rs=abs(rs).^2;
n9;;x%6 .I g}f9dB,F Rn=(Rp+Rs)/2;
}P3tn ~O 3D[PNW~ Tp=1-Rp;
X1tXqHJF} ^VT1vu
%03 Ts=1-Rs;
v_S4hz6w\ hCYQGx0 Tn=(Tp+Ts)/2;
jEz+1Nl) I5AO?BzJ figure(2)
|BZDhd9<{ 7B3w\ subplot(1,2,1);
B$?qQ|0:= ];+#i"l plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
+fKtG]$ nH !3(X* legend('R_p','R_s','R_n')
GO! uwo: k#}g,0@ xlabel('\theta_i')
d {2 qP*}.Sqk7 ylabel('Amplitude')
kHc<* L_ V ~n!!jM:N title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
BsZ{|,oQnZ iO1nwl !# axis([0 90 0 1])
I /x iT ,}D}oo* grid on
!.5 ),2 \
u+xa{b| subplot(1,2,2);
mV9A {h *B)yy[8j+ plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
e yw'7 &ogt2<1W legend('T_p','T_s','T_n')
2 4\g bv< GgFi9Ffj xlabel('\theta_i')
|I[7,`C ~ b 1^n KB ylabel('Amplitude')
I12KT~z<r ~bA,GfSn0 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
@fn 6<3 &;q<M_< axis([0 90 0 1])
nX%AeDB AT J>S3sP grid on
73NZ:h%= M7y| EB)) 
L0h
G sroGER . 
@6xGJ,s
