利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
juI)Do 2_ %, S{9q 1、光疏射向光密
PQUJUs 9&{z?* clear
'e_e*.z3 J%lrXm(l{ close all
RL[F 9g vdigw.=z n1=1,n2=1.45;
ZYt1V "2VJ $)$ _}^.k theta=0:0.1:90;
WcY $=\7 _
^'QHWP a=theta*pi/180;
NLu[<u U* +F.@n_}p-I rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
/|3~LvIt= l[u17,]S rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
tN.$4+ VK!HuO9l tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
0+&K; >f4[OBc ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
5gk Q6& m x3sX=jIW_ figure(1)
Cm]\5}Py `q`ah_ subplot(1,2,1);
>cp v4Pgm RI+Y+z plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
8 llXpe ~7FS'!W,F legend('r_p','r_s','|r_p|','|r_s|')
]~~G<Yh:= M!i*DU+SE xlabel('\theta_i')
m WHyk "l = l:k($%% ylabel('Amplitude')
y} is=h3 B_D0yhh title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
UFLx'VX d jT/P+2hMW axis([0 90 -1 1])
l`uMtv/Wp ~2 nt33" grid on
j5og}P q: /r2S1"(q subplot(1,2,2);
irvd>^&jDC =>kg] plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
%YbL%i|U fwmXIpteK legend('t_p','t_s','|t_p|','|t_s|')
?!<Q8= |ICn/r~ xlabel('\theta_i')
Ya> AI.!K J{\(Y#|rHs ylabel('Amplitude')
vexF|'!}0# ~]].i~EV( title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
s: iBl/N} u+
hRaI;v axis([0 90 0 1])
cNN0-<#c d`C$v j grid on
h?$ J;xn J"@ X>n Rp=abs(rp).^2;
@2mWNYHR*> c##tP*( Rs=abs(rs).^2;
t)uxW
7 iBtG@M Rn=(Rp+Rs)/2;
?3qp?ea 8)1=5 n Tp=1-Rp;
Bo14t*( hW9! Ts=1-Rs;
]N ;n q 4)>UTMF Tn=(Tp+Ts)/2;
Nr+~3:3 p_%,JD figure(2)
ylb)SXBf WPrBK{B`o subplot(1,2,1);
|z"$^|@d? 9+$IulOvk plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
/R?[/`)f& ;MK|l,aIQ legend('R_p','R_s','R_n')
^TjC }Y|M+0 xlabel('\theta_i')
Slj
U=, $} 0!dR2 ylabel('Amplitude')
lrwQ
>N Nd_A8H,&B title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
h oM%|,0 /OB) \{- axis([0 90 0 1])
Iz83T9I& 8DMqjt3B grid on
),y`Iw UjLq[,_! subplot(1,2,2);
mPqK k UZmUYSu; plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
!syyOfu`} 0%%y9;o legend('T_p','T_s','T_n')
7=yjd)Iy9m )4-!]NsV xlabel('\theta_i')
rBU)@I pDG [i.c;'Wy/ ylabel('Amplitude')
K'e,9P{ :SW
vH- ] title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
)V+/@ 4 #Ef! X axis([0 90 0 1])
LR!% iP &/2+'wCp5 grid on
Y~Vc|zM^( 2xiE#l-V2 
0VNpd~G$ !4gHv4v ; 
#8
^b]
