利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
VAkZ@
u3'~ qk}(E#.>F\ 1、光疏射向光密
^X2U
A{ 3{e7j6u\ clear
]RYk Y7>` 5 #jna9Xc close all
om 3$= % : ?_N n1=1,n2=1.45;
Z4S0{:XY `x;8,7W;B theta=0:0.1:90;
.54E*V1 M6n9>aW4 a=theta*pi/180;
Vp3
9`m-W f"XFf@! rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
p #vZYwe=L ^B9rt\,q rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
99Xbp P5 5 i!zh9,i>M tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
\EH: FM}l, d~L`*"/ )[ ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
c*nH = ?h `Ned0P figure(1)
tT7< V{i4 }5n((7@X subplot(1,2,1);
F|._'i+B! !\{&^,y plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
bDdJh}Vz e'<pw^I\ legend('r_p','r_s','|r_p|','|r_s|')
S@suPkQ<> s>sIji xlabel('\theta_i')
`?{Hs+4P5 Sz0M8fYT] ylabel('Amplitude')
mT
N6-V c"aiZ(aP title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
|xI\)V E^ />fP )56* axis([0 90 -1 1])
UA4Q9<>~ ]a%
*$TF grid on
U_a)g
X HgG-r&r!2 subplot(1,2,2);
C]aa^_Ldd- ' 8~cf plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
G~ZDXQ>5CP ]2n&DJu legend('t_p','t_s','|t_p|','|t_s|')
W(*:8}m,p Vv(!Ki} xlabel('\theta_i')
o/I <)sa b6D}GuW ylabel('Amplitude')
=J.)xDx* OwIW;8Z title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
%UJ! (_ G'XlsyaWrb axis([0 90 0 1])
{]ZZ] (_ov _3 grid on
bw M>#@H b5YjhRimS Rp=abs(rp).^2;
SsjO1F ,hK0F3?H> Rs=abs(rs).^2;
D={|&:`L e kr 7f<;rmJ Rn=(Rp+Rs)/2;
P,RCbPC4 8k'UEf`'( Tp=1-Rp;
K5(?6hr; *IzcW6 [9 Ts=1-Rs;
G8`q-B}q =Mq=\T Tn=(Tp+Ts)/2;
p#.B Fy >HnD'y* figure(2)
2F-!SI IS7g{:}=p subplot(1,2,1);
c1wM " " }D u As plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
!B?/6XRUx j' -akXo< legend('R_p','R_s','R_n')
@Z!leya m E66 e4?" xlabel('\theta_i')
Y',s|M1})\ IoNZ' g?d ylabel('Amplitude')
i o
cr .hgc 1 title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
1rC8] M.N q
/|<>s axis([0 90 0 1])
n6WSTh }j TE gog grid on
YP~d1BWvf ~! ]FF}6 subplot(1,2,2);
))!Z2PfD 1o o'\ plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
qB@]$ g=pz&cz;>\ legend('T_p','T_s','T_n')
G 5zZf ~r OM]p"Jd xlabel('\theta_i')
=(*Eh=Pw IGql^,b ylabel('Amplitude')
MLmc]nL= .D^k0V title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
F=B>0Q5 ? $p GG axis([0 90 0 1])
/'E+(Y&:J e }/c`7M grid on
u4$d# 0sA O<f_-n@G| 
PK3)M'[ n_v c}ame 
e @p` -;<
