利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
|fhYft .&2Nm&y$ K 1、光疏射向光密
W&`{3 L 1Q>D^yPI[ clear
|';oIYs| $ )E>nr
Z close all
MR`lF-|a| |p\vH#6y+ n1=1,n2=1.45;
{$fd?| 9h 9XEP:}5, theta=0:0.1:90;
u-%|ZSg PRQEk.C a=theta*pi/180;
ZW)_dg 9 #[$zbZ(I>: rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
NL} Q3Vv1. =s5g9n+7 rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
HBp $
|Ta-D++]' tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
,!7\?=G6}v QuWW a|g^. ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
|rr<4>)X (YC{BM} figure(1)
Y~"5HP | PPMAj@B}V subplot(1,2,1);
On*pI37(\ 5R}K8"d plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
Tky P_* v-ZTl4j$ legend('r_p','r_s','|r_p|','|r_s|')
u|{(m_ "H b<E+5;u xlabel('\theta_i')
,6 AnuA @/f'i9?oM` ylabel('Amplitude')
> Sc/E}3 _-c1" Kl title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
Ygs:Ox"[-G EV .F/W h axis([0 90 -1 1])
-Wm'@4bH d$qi. %<kh grid on
$ /#F9>e Z ~`c (7 subplot(1,2,2);
#ni:Bwtl{ Bj *\)lG<
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
WNo7`)Kx 6QM$aLLP? legend('t_p','t_s','|t_p|','|t_s|')
P1qQ )-J f(UB$^4 xlabel('\theta_i')
j {&$_ L"Dos + ylabel('Amplitude')
,Z$!:U s8V:;$ ! title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
T-a [ &qyXi [vw axis([0 90 0 1])
vTsMq>%,< V:<Z grid on
;6}> Shs k3CHv =U{ Rp=abs(rp).^2;
iBh.&K{j S;iJQS Rs=abs(rs).^2;
\&\U&^? @fh:lsw Rn=(Rp+Rs)/2;
rd1EA|T Hj97&C{Q ^ Tp=1-Rp;
{M
^5w '_B;e=v` Ts=1-Rs;
>qS2ha >UnLq:G Tn=(Tp+Ts)/2;
:j&- Lc SQT]' figure(2)
YkF52_^_ 3g87 i r subplot(1,2,1);
~B\O {5W $bFH%EA. plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
hV}C.- 6h lS{ ^*(a legend('R_p','R_s','R_n')
p735i`8 'OI(MuSn xlabel('\theta_i')
V='A;gs GJIZu&C ylabel('Amplitude')
3R<VpN){ FBeo@ title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
6%Pvh- ~_ !CUM* <iV axis([0 90 0 1])
sL], @z8<k nhy:5eSK grid on
-,q
qQf VQ;'SY:` subplot(1,2,2);
WI1DL&*B@<
[L=M=;{4 plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
nQ@<[KNd q{l %k legend('T_p','T_s','T_n')
^G14Z5 . "2 ZuI; w xlabel('\theta_i')
R>hL.+l. yG2rAG_ G& ylabel('Amplitude')
-_BX\iP{ VE )) `? title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
4 9=L9: rN'8,CV axis([0 90 0 1])
C9 j{:& QYps5zcn grid on
Ym?VF{e, K0|:+s@ u 
vG:S(/\> "a-;?S& 
gI00@p:m
