利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
{H s" "/sb l!e8=Ql J 1、光疏射向光密
`84pql, v<bq1QG clear
*0 M#{HQ fVv#| close all
G3&ES3L 9PhdoRE b n1=1,n2=1.45;
tVQ fR*= T.="a2iS2 theta=0:0.1:90;
c/
%5IhX? ElAJR4'{*i a=theta*pi/180;
6' ye-}vD- ^zkTV_,cRp rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
w~ M5)b ep<A d rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
{LTb-CB >|SB]'C | tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
HRi ~ TZ?\ xG0IA 7 ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
{n%-^9b1{& FW& P`Iu figure(1)
9x.vz Pa 2 HFy2 subplot(1,2,1);
t+Mr1e 3P Twpq1 plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
@8C^[fDL ,S;?3? a legend('r_p','r_s','|r_p|','|r_s|')
[ YP8 z~ ,R0@`t1 p xlabel('\theta_i')
* JGm b_ Sh#d& ylabel('Amplitude')
>JS\H6 n"Ec %n title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
ba|x?kz K,tmh1 axis([0 90 -1 1])
%*OKhrM =O1CxsKt6 grid on
&5 /`6-K DU$ ] e1 subplot(1,2,2);
&J^@TgqL^ '[JrP<~^o plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
,0[ 8/)$M cRz7.9-< legend('t_p','t_s','|t_p|','|t_s|')
{@g3AG% Tz58@VY V xlabel('\theta_i')
#0`2wuo
{ CU6rw+Vax ylabel('Amplitude')
t) g1ICt z _! P0` title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
(Z.K3 "m})~va axis([0 90 0 1])
TJ7on.; Svn|vH grid on
?yddr`?W 2j-|.l c Rp=abs(rp).^2;
aGNt?)8WPZ h+zJ"\ Rs=abs(rs).^2;
R|{AIa{} `y0ZFh1>X Rn=(Rp+Rs)/2;
/7|u2!#Ui 8gJ"7,}-' Tp=1-Rp;
JO5~Vj_" +Es3iE @
Ts=1-Rs;
NWwfNb> MR%M[SK1 Tn=(Tp+Ts)/2;
`kyr\+hp N4!YaQQ;} figure(2)
C2AP e{O m W subplot(1,2,1);
cg7NtY W5$jIQ}Bw plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
\% &QIe;:k $ePAsJ legend('R_p','R_s','R_n')
Mp?Ev. /-E>5 w U xlabel('\theta_i')
9+@_ZI- {i~qm4+o ylabel('Amplitude')
^w^cYM, k<"ZNQm$. title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
x*i5g` jx =Z2U axis([0 90 0 1])
+ %eMm.( Cv{rd##Y8 grid on
IyOujdKa 6;dB subplot(1,2,2);
y2\, L " FcA:7 + plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
"pdG%$ S#!PDg legend('T_p','T_s','T_n')
{\`#,[ 9g"a`a?c xlabel('\theta_i')
PQ@(p% PLg`\| ylabel('Amplitude')
H h$D:ZO $&n!j'C: title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
`iv,aQ ' T $GhE axis([0 90 0 1])
Da _g3z 7 -gt V # grid on
n0@e%=H )I nE Qw6q~je 
L2P~moVIi i4'?/UPc 
s=~r. x
