利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
( tOhuSW _^NaP 1、光疏射向光密
U1 X"UN) )\ceanS clear
l
EzN IQ~qiFCf close all
Ul9^"o .=TXi<8Brw n1=1,n2=1.45;
BZHoRd{E H \U]K!K= theta=0:0.1:90;
@$n
$f kx?Yin8K a=theta*pi/180;
kj[box N 0bM_EC rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
b<~-s sL7a @R50M (@W rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
D|*w6p("z g:a[N%[C tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
'JJKnE zQ S Fd _k9 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
f"P866@oWn aG^E^^Y figure(1)
B K+P cI @qt>& subplot(1,2,1);
T+&fUhSy m2jts(stp plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
Z>bNU A1-,b.Ni legend('r_p','r_s','|r_p|','|r_s|')
*c [w9(fU <tF q^qB xlabel('\theta_i')
?x&}ammid x@Z{5w_a ylabel('Amplitude')
)Bn
}|6` c{I]!y^! title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
#eOHe4Vt 4Ii 5V
c axis([0 90 -1 1])
P>iZ gv hE 5? G; grid on
]zaTX?F: BLt58LYGX subplot(1,2,2);
UtTlJb{-j 1L4-;HYJm plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
b)Da6fp #<b\B qYG legend('t_p','t_s','|t_p|','|t_s|')
g:)iEw> a */aQ+%>jf xlabel('\theta_i')
G&^8)S@1 (9I(e^@] ylabel('Amplitude')
u1 M8nb fEXFnQ# title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
j Db\4QyC #J&3Zds axis([0 90 0 1])
@ewaj! NL$z 4m0 grid on
EA>.SSs! :G@z?ZJ[ Rp=abs(rp).^2;
b.RU%Y#>\ rP7
QW)NF Rs=abs(rs).^2;
ydYsmTr ?#
FYF\P Rn=(Rp+Rs)/2;
dv\bkDF4A ial{A6X Tp=1-Rp;
=bC' >qw} 9gw;MFP)D Ts=1-Rs;
~9FL]qo :Q=y'< Tn=(Tp+Ts)/2;
ZN $%\,< L}rZ1wV 6 figure(2)
HP]5"ziA CYy=f - subplot(1,2,1);
=YgH- { N s0,Z#Z+ plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
: Q("
mP
}<{oh`x legend('R_p','R_s','R_n')
mN]WjfII #W)m({ } xlabel('\theta_i')
]- X6Cl NJsa TBT ylabel('Amplitude')
gjG SI'M0B GxD`M2 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
(@XQ]S}L WyatHC axis([0 90 0 1])
%50)?J=zB AE~a=e\x grid on
)m;*d7l~p W..*!UGl subplot(1,2,2);
LEPLoF3, 2 kVZlt'y plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
"JSIn"/ v[ML=pL legend('T_p','T_s','T_n')
P*iC#w]m xA7~"q&u xlabel('\theta_i')
rIFW1`N}i lH= |Qu ylabel('Amplitude')
o FP8s[B K9Mz4K_ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
C\5G43` 6 ScB:8M axis([0 90 0 1])
} uS0N$4 ]m1p<*0I$ grid on
kR9 7 )}Y Pp!4Ak4TT9 
aBReIK o
"LB
MYZ 
:5& D 6
