利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
r$={_M$ H ** tMq 1、光疏射向光密
*XhlIQ FFcCoPX_ clear
' 2 r 3E|||3rf close all
H<_Tn$<zH. c]|vg=W n1=1,n2=1.45;
|})s 0TU Hloe7+5UD theta=0:0.1:90;
]H
n:c'aT kzRvLs4xM a=theta*pi/180;
Dz~0 ( NlEyT9 rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
OOSf<I*>
pRobx rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
.A< HM } 3k)xzv%r` tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
$?PI>9g! gt}Atr6>_ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
SF:98#pg Ww-x +U\l figure(1)
g9pKoi|\E F6q= W#~ subplot(1,2,1);
I3nE]OcW@ w"9h_;'C_ plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
KY JI! 1
.]& legend('r_p','r_s','|r_p|','|r_s|')
b \}a
,H7_eVLWR xlabel('\theta_i')
89&9VX^A l ubsL I ylabel('Amplitude')
lm4A%4- db 8?nn4]P title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
-t4:%-wv cn} C I axis([0 90 -1 1])
. :j{d} p} XS &P c grid on
8UjIC4' w PR Ns9^ subplot(1,2,2);
,'(|, f42 R@3HlGuRKw plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
v)du] u<:uL legend('t_p','t_s','|t_p|','|t_s|')
YdCl h2zSOY{su xlabel('\theta_i')
Eq)b=5qrG? U%@PY9# ylabel('Amplitude')
F7=9> , `C>h]H( title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
\ZcI{t'a LOkgeJuWv axis([0 90 0 1])
E^rBs2;9 _E?tVx.6 grid on
R9A:"sJ 2`]c&k;] Rp=abs(rp).^2;
ELwXp|L \2/X$x<?X Rs=abs(rs).^2;
Vl1.]'p_ ] & ]G Rn=(Rp+Rs)/2;
i`@cVYsL Ye On Tp=1-Rp;
V Dn rm* 2*D 2j w Ts=1-Rs;
\5 }PF+)| [ 6VM4l" Tn=(Tp+Ts)/2;
c_qox %@:>hQ2 ; figure(2)
#6@hVR. H,KH}25 subplot(1,2,1);
[Z/P[370 ]CP5s5 plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
|iJz[% RgoF4 g+@ legend('R_p','R_s','R_n')
w-WAgAch vltE2mb xlabel('\theta_i')
auN 8M. -aJ(-Np$f ylabel('Amplitude')
C3 "EZe[R aN"YEL>w title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
x%X3FbF] LF.i0^#J axis([0 90 0 1])
\_.'/<aQ yzfiH4 grid on
7x`$ A ?:vg`m!* subplot(1,2,2);
9Y2u/|!.3 _~Od G plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
z]g#2 xD2 .n IGs'P legend('T_p','T_s','T_n')
FsUH/Y
y '(r?($s xlabel('\theta_i')
# d"M(nt rMG[,:V ylabel('Amplitude')
WuVsW3@ hlyh8=Z6o title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
GX19GI@k \t{4pobo axis([0 90 0 1])
B_Wig2xH0 !pe[H*Cy grid on
FBP # _"z P ` <TO 
2FL_!;p;2E :6C R~p 
vyc<RjS_x
