利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
\b?O+;5Cj 4|NcWpaV7 1、光疏射向光密
qG0gc\C} 5K00z?kD2V clear
i~tps `3.bux~ close all
=<U'Jtu6' 8},fu3Z n1=1,n2=1.45;
c}w[ T B |S X?X theta=0:0.1:90;
s "l ^v5 Ps~)l#gue a=theta*pi/180;
^@]yiED{g
sVP2$? rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
#J%Fi).^) (ewcj\l4* rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
,62BZyT,T, ?{>5IjL)en tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
y|Y hDO rm,h\ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
= %wBC; 6H:EBj54? figure(1)
/!-ypIY
o
Hrx$>W] subplot(1,2,1);
osTin*T. E9j(%kQ2 plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
X a[lX8$zL q,(hs]\@ legend('r_p','r_s','|r_p|','|r_s|')
,@ t#)HV y |Tv;v1L xlabel('\theta_i')
x 2+ M0 }g G[!<mh4h| ylabel('Amplitude')
AI2 >{V u!cA_, title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
IO?6F@( 'deqF|Iox axis([0 90 -1 1])
:#UN^ "(m} @m" P_1`* grid on
V,:~FufM^ V _pKe~ subplot(1,2,2);
VB{G% !} 5v#_2 Ih plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
KT8]/T`U |C;*GeyS;J legend('t_p','t_s','|t_p|','|t_s|')
ZAMS;e+e i- r y5x xlabel('\theta_i')
G K:pt8= kam \dn04 ylabel('Amplitude')
oOK&+r7 _1P8rc"Dx title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
1 *$-. 0G/_"} @ axis([0 90 0 1])
S31+ j:" ?:sQ]S/Er grid on
yy|F6Pq3` PiwI.c Rp=abs(rp).^2;
@0q* 50 +jX.::UPm Rs=abs(rs).^2;
\+sP<'~M B!z5P" C(~ Rn=(Rp+Rs)/2;
{siIRl2& MM4Eq>F/ Tp=1-Rp;
9~WjCa*,& d!]_n|B@9 Ts=1-Rs;
< 5#}EiT5 &(HIBF'O Tn=(Tp+Ts)/2;
Fs) m;C <3c|S_|L*m figure(2)
9`sIE _%+ `1"Xj ^
YM subplot(1,2,1);
4';(\42 C8.MoFfhe plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
j#E&u*IR rui]_Fn]I legend('R_p','R_s','R_n')
Kzgn h gc <',bqsg[ xlabel('\theta_i')
"+:IA|1wD au 5qbP ylabel('Amplitude')
w"K; e (S H:M;H =0 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
l@ W?qw +cnBE v~y axis([0 90 0 1])
tB7g.)yZb ,BG
L|5?3z grid on
V tr5<:eEx }H<87zH subplot(1,2,2);
|VD}: " tUF,G (< plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
IQ"9# {o ON<X1e U legend('T_p','T_s','T_n')
w*"h#^1z JgY#W1> xlabel('\theta_i')
!*l5%H CERT`W%o ylabel('Amplitude')
:=K <2 <i!7f26r title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
{y
k0Zef_ S
2SJFp axis([0 90 0 1])
@d86l.= Fz4g:8qdA grid on
R
s)Nz< d #$[}JiuL/ 
'u{DFMB-A " ~&d= f0m 
I,xV&j+<
