利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
nO;ox*Bk+8 ,t|_Nc
1、光疏射向光密
24f N3 8jiBLZkRf clear
xscR Bx "V,dH%&j close all
h}}7_I9 9QkssI n1=1,n2=1.45;
c~6ywuq+M` Y"dTm;& theta=0:0.1:90;
m-q O yt >f;oY9 {m a=theta*pi/180;
9gZS )MZ
@;[. #hK rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
[Uj,, y.wB V3axwg_ rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
!H @nAz e)i-$0L" tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
]agd Vr^ 5 y p ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
M[e{(iQ: $<f+CtD4 figure(1)
{s?hXB vmW`}FKW subplot(1,2,1);
o/!a7>xO4 V_9\Ax'X plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
f&X M|Bg w!Z,3Yc) legend('r_p','r_s','|r_p|','|r_s|')
?MiMwVR |~=4Z rcCP xlabel('\theta_i')
=J )(=, Kn\$\?u ylabel('Amplitude')
WhV>]B2+" lPz5.(5' title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
[VOw:|Tt "azrc C axis([0 90 -1 1])
eT8h:+k |mz0
] grid on
X<H+Z2d ZaFqGcS~ subplot(1,2,2);
WW~QK2o-@ dw7h@9\ y plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
?\a'; @h `y.i(~^1 legend('t_p','t_s','|t_p|','|t_s|')
QSOJHRl=C @2 SL$0!QA xlabel('\theta_i')
~ o5h}OU" b*< *,Ds/G ylabel('Amplitude')
dr4 m}v. KYFkO~N title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
}iR!uhi# yMf["A vG axis([0 90 0 1])
uTP4r Js !Zk\O grid on
Zg!E}B:z om,=.,|Ld Rp=abs(rp).^2;
bJ6v5YA% ;[-y>qU0 Rs=abs(rs).^2;
jyjK~ !0 ! q*]_1 Rn=(Rp+Rs)/2;
PPNZ(j /?:] f Tp=1-Rp;
1
BVpv7@ lb #`f,r> Ts=1-Rs;
5LxzET"P :*lB86Ly Tn=(Tp+Ts)/2;
bC]GL$ph9* qd?k#Gw& figure(2)
m g4: N a~[]Ye@H subplot(1,2,1);
_l.kb fp@ 'I)E.D oF plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
e@{8 G^o>D 5nG$6H w legend('R_p','R_s','R_n')
C)Hb= -o0~xspF xlabel('\theta_i')
{2%@I~US LgS.%Mn ylabel('Amplitude')
YQ<O .E ?gOZY\[ma title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
V9kL\Ys <l VW; l7 axis([0 90 0 1])
2qUC@d<K K) t+ lJ grid on
B (dq$+4 HzF]hm, subplot(1,2,2);
}c*6|B@f 0sKY;( plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
KcX] g*wy N{6Lvq [8 legend('T_p','T_s','T_n')
zWI C4: *a4nd_! xlabel('\theta_i')
9:l>FoXS h+h` 0(z ylabel('Amplitude')
O
hcPlr ^+ +ec> title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
c o!#. ro4 XA1 axis([0 90 0 1])
9rsty{J8 g&"__~dS-F grid on
NI 136P 3YF*TxKx 
\[w 82%U CCZ]`*wJ 
J m8#M z
