利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
tpA-IL?KQw n>eDN\5 1、光疏射向光密
a)r["*bTx xWY\,'+Q clear
pO )EYla9 X0j> g^b8 close all
zq$L[ X PPG+~.7 n1=1,n2=1.45;
@ls/3`E/5E G+2fmVB*X theta=0:0.1:90;
.6rbn8h *mj=kJ7(
a=theta*pi/180;
rt*>)GI]b 5K?/-0yG rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
<uF j5. _0^>^he rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
->;2CcpHB K0v S tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
t%^&b'/Z a"O;DYh ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
;q%z\g A havmhS )O figure(1)
B<.\^f uS { )b subplot(1,2,1);
mc2uI-W E+< GsN] plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
St!0MdCH c}S<< LR legend('r_p','r_s','|r_p|','|r_s|')
MK
Sw
A+H8\ew2, xlabel('\theta_i')
)
5I j rZB='(? ylabel('Amplitude')
r~Q E}00@^ 1D [>oK\ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
6/g
82kqpk `w4'DB-R) axis([0 90 -1 1])
,S(Z\[x0 =Sr<d|\O grid on
UB>BVBCt cotySio$ subplot(1,2,2);
->IZZ 5G< /4 K d plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
*a8 <cf mb3aUFxA; legend('t_p','t_s','|t_p|','|t_s|')
L|(U %$ SQ+r'g xlabel('\theta_i')
BL>~~ UB8n,+R ylabel('Amplitude')
qG~6YCqii s%vy^x29 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
jd5kkX8= Q q j9o2 axis([0 90 0 1])
x7gd6"10^ :nl,A c grid on
y eIS} O ?A .ah Rp=abs(rp).^2;
z"Q tP[_m sL\ {.ad5 Rs=abs(rs).^2;
S{MB$JA hiV!/}'7 Rn=(Rp+Rs)/2;
GC r]x ' pH3\X
cn Tp=1-Rp;
tV pXA'"!x U6H3T0# Ts=1-Rs;
M5:*aCN6P e~'z;% O~ Tn=(Tp+Ts)/2;
Tz9 (</y -nUK%a"(D figure(2)
SE i\H$ ! )ryP K"V subplot(1,2,1);
D ZZRu8~ SS _6VE*sI plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
~g#/q~UE DYIp 2-K legend('R_p','R_s','R_n')
{w"Cr0F, ld({1jpX, xlabel('\theta_i')
7042 ?\\= {(Mmv[y ylabel('Amplitude')
#
&5. q;sZwp< title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
\4<|QE Ets6tM` axis([0 90 0 1])
EX , {1^h &IRM<A!8 grid on
ku}`PS0UGd 7\nXJ381 subplot(1,2,2);
;+o6"ky5 D Vg$rm` plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
4 eSFpy 1 $II ~tO legend('T_p','T_s','T_n')
nsV= eFA,xzp xlabel('\theta_i')
Fe:
~M?] 'q}f3u > ylabel('Amplitude')
"j8=%J{ (w(k*b/ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
cmCD}Skk Y8lZ]IB
axis([0 90 0 1])
9Nv?j=*$ =h
~n5wQG grid on
&?xmu204 FQ47j)p; 
?0JNaf fWri7|"0h 
0u( 0*Xl
