利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
6' *Uo: ] 68GGS`& 1、光疏射向光密
~S_IU">E `rdfROKv clear
898wZ{ 9
= ObtD" close all
JTB~nd> 4Hpu EV8Q n1=1,n2=1.45;
>h( rd1
:E&T}RN theta=0:0.1:90;
yz$1qEII`q U9[A( a=theta*pi/180;
yGG\[I;7 iE gM ~ rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
R}6la.mQ xPJJ
!mY rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
$H@)hY8wA H_=[~mJ tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
<07W&`Dw =yhfL2`aw ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
V >uW|6 NE%yv,B figure(1)
&trh\\I" WCl;#= subplot(1,2,1);
t s\>_/ 5{/Pn%5 plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
PZg]zz=V4 PyxN _agf legend('r_p','r_s','|r_p|','|r_s|')
d#:J\2V"R p}|wO&4h xlabel('\theta_i')
dB/I2uGl> @%G"i:HZ& ylabel('Amplitude')
sH,)e'0 lbU+a$ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
)bU")
m9w
; a axis([0 90 -1 1])
SA n=9MG |A/_Qe|s2 grid on
ZjW| qb
!,!tNs1 K subplot(1,2,2);
WM
)g(i~( ;U3Vows plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
\0T*msYQ 'p+QFT>Ca legend('t_p','t_s','|t_p|','|t_s|')
Pc4c Sw#5 )mVYqlU" xlabel('\theta_i')
-?` l<y( U9sub6w 6 ylabel('Amplitude')
s:F+bG }| Q09~vFBg title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
[ \n.[4gq" |F {E4mg(o axis([0 90 0 1])
cD@lor j g}r5ohqC# grid on
.V:< w~=b `2,_"9Z( Rp=abs(rp).^2;
Ki;5 =) GJfNO- Rs=abs(rs).^2;
Riuv@i^6K ,k' 6<Hw Rn=(Rp+Rs)/2;
-"9)c^KVx M|=$~@9#X Tp=1-Rp;
^F9zS `Yz2 #f;1f8yrN Ts=1-Rs;
F
) ~pw <% uZwk># Tn=(Tp+Ts)/2;
r9*H-V $ T sW6 w figure(2)
.h^Ld,Chj NVZNQ { subplot(1,2,1);
P,=J"%a - =C1Qo#QQ% plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
> .~k?_Of xi= uX xl legend('R_p','R_s','R_n')
D|3QLG y+R$pz X xlabel('\theta_i')
#| E. y^IC \ jdO,-( ylabel('Amplitude')
2dW-WHaM m|FONQ,@D title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
{\Y,UANZ
=H?5fT^
axis([0 90 0 1])
)3 ">%1R bGO_y]Pc grid on
|q)Q <%VS' |]b/5s;> subplot(1,2,2);
io_64K+K < _uv!N plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
30`H
X v@ baG I( Dk legend('T_p','T_s','T_n')
R$`T"C" Y %8QFM xlabel('\theta_i')
K x!|4ya , ~h|L;E" ylabel('Amplitude')
g&5VorGx IWnW(>V title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
:Xr3 3 H} rP{`m axis([0 90 0 1])
y%Ah"UY :PkSX*E[q grid on
nwH|Hs riU 5|z[%x~f 
hV>Ey^Ty 13
p0w 
C,. {y`s'
