利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
BEWro|]cM .]P2}w)x? 1、光疏射向光密
qI h #~ ;UrK {>B clear
s"i~6})K<$ B oqJ
close all
n$4|P O$X M'q'$)e n1=1,n2=1.45;
qK?$= h. If;R?j0;Q theta=0:0.1:90;
:6Pnie kh3<V'k] a=theta*pi/180;
Ra[ >P _ t,h{+lYU rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
?g3 ]~;# ]9*;;4M g rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
'a#mViPTQ) `4V"s-T' tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
zmiZ] uq $w)~xE5; ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
.%'Z~|K 4 {oUAP1V ^ figure(1)
;wR ' z$8 Z19m@vMsIP subplot(1,2,1);
e3 v5,. F\:(*1C plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
Hm
fXe j)by }} legend('r_p','r_s','|r_p|','|r_s|')
? gSSli[ -W c~ B3E| xlabel('\theta_i')
7J|&U2}c rC BfD ylabel('Amplitude')
6;+jIkkD) m !O;>D title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
26L~X[F ;Lo&}U3F,! axis([0 90 -1 1])
BAV> o |-K U>P|X=) grid on
!^y y0`k6 KKEN'-3 subplot(1,2,2);
I% "'*7 U # ? ?%B plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
f:wd&V r &%.z*q legend('t_p','t_s','|t_p|','|t_s|')
he$XLTmr: (UYF%MA}" xlabel('\theta_i')
G[h(xp?,l +}@ 8p[`) ylabel('Amplitude')
!%PWig- ,*Z. title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
g%a|q~ ) PB 53myDQ axis([0 90 0 1])
"e6|"w@8 lA^+Flh grid on
1J}8sG2` `f9gC3Hk Rp=abs(rp).^2;
2p!"p`b~ ~AZWds(,N Rs=abs(rs).^2;
*5V Xyt2 / nZ;v4 Rn=(Rp+Rs)/2;
G7@ O`N8' 0F:1\9f5 Tp=1-Rp;
xW_yLbE 7N=-Y >$X Ts=1-Rs;
j1-,Sqi ZA(T
Tn=(Tp+Ts)/2;
%o w^dzW "T S figure(2)
{9LWUCpsf jS<_ ) subplot(1,2,1);
P"_$uO( 5x ',I0ih#Ls plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
v [njdP ][;G=oCT legend('R_p','R_s','R_n')
RA O`i>@ 9z>z3,ftN xlabel('\theta_i')
]N{0:Va@D > \Sr{p5KR ylabel('Amplitude')
`ln1$ ArK]0$T title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
)|bC^{kH!l z=7|{ G axis([0 90 0 1])
'gs o'&Uaj [KJ
q grid on
P\nC?!Q%c 58tVx' 1y subplot(1,2,2);
H%F>@(U EZD y+6b plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
od' /% sTRJ:fR legend('T_p','T_s','T_n')
{aY Y85j *QGm/ /b xlabel('\theta_i')
zj{r^D$ 3< Od0J ylabel('Amplitude')
`D$Jv N #/2$+x title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
&opd2 >slm$~rv axis([0 90 0 1])
hr05L<?H kzn[
=P grid on
Z;l`YK^- *hLQ 
NB_ )ZEmF 2d>hi32I 
_PFnh) o
