利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
Z4bN|\I (/ Z~0hA[Q 1、光疏射向光密
ObMsnc n )d0&i E`@ clear
)^g}'V=vIr BP*gnXj close all
e +4p__TmZ duV|'ntr n1=1,n2=1.45;
qs 52)$ BalOph4M[ theta=0:0.1:90;
Rm)vY}v hG&RGN_<6+ a=theta*pi/180;
m->
chOu~| F)5QpDmqb
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
8+v6%,K2 8p> %}LX/ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
mkrvWZ jZX fC]+C(*d tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
)Qb1$%r. Fov/?:f$ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
j$& k;S 4
lJ@qhV figure(1)
iN+p>3w^l =EsKFt" subplot(1,2,1);
p5c'gziR X&
O
o1y plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
A{mbL2AxwC = J).(E89 legend('r_p','r_s','|r_p|','|r_s|')
J0mY=vX 2,q^O3F xlabel('\theta_i')
1Qhx$If~ -md2Z0^ Kc ylabel('Amplitude')
;&;coH8` yTkYPx title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
} 9<aX
Y, u)R>ozER axis([0 90 -1 1])
NVeb,Pf I)_072^O grid on
vlp]!7 v .eK1xwhJ subplot(1,2,2);
#x)G2T'? `Ft`8=( plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
) 8st
Rj+}L ~" legend('t_p','t_s','|t_p|','|t_s|')
"g>uNtt~ MvL%*("4b xlabel('\theta_i')
ElqHZ$a? W4|1wd}.t ylabel('Amplitude')
8Tt2T }
Y OA4NXl' title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
{BY`Wu:w H6b omp" axis([0 90 0 1])
<u u1e@P x # BUIi grid on
'[`.&-; g0cC w2S Rp=abs(rp).^2;
NAx( Qi3
2Z7smDJ Rs=abs(rs).^2;
u?Iop/b 08yTTt76t Rn=(Rp+Rs)/2;
C.%iQ x`
zRsT6u Tp=1-Rp;
)amdRc *~prI1e ( Ts=1-Rs;
oel3H 5Nz B/i,QBPF] Tn=(Tp+Ts)/2;
JEU?@J71O e>uV8!u figure(2)
[^1;8Tbk cV&(L]k >` subplot(1,2,1);
7bQ#M )} \P&'4y~PL plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
#++MoW}'g @D[;$YEk legend('R_p','R_s','R_n')
q<,? :g$k nU[RO y5 xlabel('\theta_i')
s-*8= -
Kj$A@~x ylabel('Amplitude')
(ai E!c )
;-AT^ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Vnv<]D
zC xg. d)n axis([0 90 0 1])
F 3,hx 2X]2;W)S; grid on
|+x;18 ]*sXISg1 subplot(1,2,2);
Ij@YOt |+cyb<(V J plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
-t :y y:4 TcO@q ]+S legend('T_p','T_s','T_n')
"MZVwl "E# D
GOc! xlabel('\theta_i')
Zhz.8W [AE]0cO@ ylabel('Amplitude')
\8{C$"F x8b w# title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
q,0o:nI Fg5>CppH axis([0 90 0 1])
-Ww'wH'2 Gob1V grid on
DE."XSni S7E:&E& 
{hZ_f3o D-E30b]e 
s-o0N{b?#'
