利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
&}*[-z 7@06x+! 1、光疏射向光密
hL&z"_` 7MBz&wE^f clear
1tfm\/V}ho i5:fn@& close all
=G*<WcR WJ/&Ag 1 n1=1,n2=1.45;
ZfIQ Fh> X4 xnr^ theta=0:0.1:90;
E?%rmdyhL! ?HU(0Vgn' a=theta*pi/180;
M`S >Q2{ V>z8 *28S. rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
q?JP\_o: *n}{ )Ef rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
0~"{z >s ' 7eZ,;
x tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
2}hEBw68 o2&mh T ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
9' T
nR[> BK6oW3wD/ figure(1)
ql^n=+U 7)<&,BWc subplot(1,2,1);
'BtvT[KM _18) XR plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
[N.4 i"
Cd rr9N(AoxW legend('r_p','r_s','|r_p|','|r_s|')
k${25*M!3 $xNZ.|al xlabel('\theta_i')
?3 "lI,!0 A"d=,?yE ylabel('Amplitude')
51 l : {0QNqjue title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
}-! 0d*I r m\] axis([0 90 -1 1])
8)/ d8@ f6u<.b grid on
K#t T \ 0.=dOz r subplot(1,2,2);
RMDzPda. ={B%qq plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
d3<7 t 5{L~e>oS9 legend('t_p','t_s','|t_p|','|t_s|')
KZ>cfv- &a >- 0Rq[) xlabel('\theta_i')
4*P#3 B'@V J>TNyVaoQ ylabel('Amplitude')
+9<"Y6 +d>?aqI\A title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
e?,n> T1_O~ < axis([0 90 0 1])
8,7^@[bzXx X @RS
/ grid on
whxTCI V 3 f@@|vZF Rp=abs(rp).^2;
9sJ=Nldq 41Hv)}Yd Rs=abs(rs).^2;
ose(#n4 0 qILb># Rn=(Rp+Rs)/2;
T\?$7$/V z{`K_s %5 Tp=1-Rp;
_wH> h$ E G5 vp(%j
Ts=1-Rs;
dhv?36uE R-LM V Tn=(Tp+Ts)/2;
}IEwGoDwNs <-;/,uu figure(2)
2*Z k^h= g&z8t;@ subplot(1,2,1);
2,aH1Xbex o=J-Ju plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
~I6 N6T Z lg"aB legend('R_p','R_s','R_n')
DK)T2{: 17$'r^t,S xlabel('\theta_i')
,2YZ B*6h{ 7 6H>ST@G| ylabel('Amplitude')
(qglD ' _d4[Olu title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Yw] 7@ v%:VV*MxF axis([0 90 0 1])
A]Bf&+V 0civXZgj grid on
\?Sv O Lq|>n [KY subplot(1,2,2);
!iO2 yp .7]P-]uOZ plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
jr,N+K(@T rk6K0TQ8 legend('T_p','T_s','T_n')
<d @9[]
/~M H]Gh xlabel('\theta_i')
N=A H S jp_|pC' ylabel('Amplitude')
fIl;qGz85 GLgf%A`5/_ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
aaP_^m O {`QA.he. axis([0 90 0 1])
)/?H]o$NU c/Xg ARCO grid on
;$W HTO( D/?Ec\ t 
lD[ 37 U! @/F61Ut 
m>%b4M
