利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
`p|vutk)U Uk6 !Sb 1、光疏射向光密
$SR]7GZ dR :iUw:V clear
%/zZ ~WIf 3wl>a#f close all
v{ Ve sf qT @IY)e n1=1,n2=1.45;
E
ET 2|*} Ghc0{M< theta=0:0.1:90;
T_T{c+,Zd$ p> S/6 [X a=theta*pi/180;
,eCXT=6 t7FQ.E,T rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
x~eEaD5m%J S I5QdX rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
>,Z{wxz J aM!# tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
I1 O?)x~ yw!`1#3. ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
LP /4e` &
jvG]>CS' figure(1)
\
6 a P.DWC'IBN subplot(1,2,1);
v19`7qgR( F9w&!yW: plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
t@M] ec uD5yw #` legend('r_p','r_s','|r_p|','|r_s|')
926oM77 Hc|U @G xlabel('\theta_i')
[rU8
#4.
:`1g{8.+ ylabel('Amplitude')
ypo=y/! +to9].O7y title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
!3# }ZC2 ]M;! ])b$ axis([0 90 -1 1])
Xm'K6JH' Y '1V(5/& grid on
^#se4qQ ,$$$_+m\ subplot(1,2,2);
U:hC! t: %QLYNuG plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
[zEP| 8-YrmP2k legend('t_p','t_s','|t_p|','|t_s|')
v"~I( kf$ W=]" ,< xlabel('\theta_i')
7g+ ] Ct+% ylabel('Amplitude')
Qe.kN dT+_ IQ ~7vk() title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
=A 6O}0z % 6.jh#C axis([0 90 0 1])
rF3]AW( 1Z8oN3 grid on
S'p`ECfVMA -$ z " 74 Rp=abs(rp).^2;
LfXr(2u T?{9Z Rs=abs(rs).^2;
o{W]mr3D ABmDSV5i Rn=(Rp+Rs)/2;
\RyA}P5 S wJ* -K- Tp=1-Rp;
UyK G$6F?3 /,$ \H Ts=1-Rs;
w QB{K 3 ?u!AHSr( Tn=(Tp+Ts)/2;
X>8?p'* G>>u#>0 figure(2)
3eP0v Kg-X]yu*0 subplot(1,2,1);
L b;vrh;A E9 q;>)} plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
8lSn*;S, aZGDtzNG5h legend('R_p','R_s','R_n')
q%Jy>IXt 4,ynt& xlabel('\theta_i')
Al=? j#J6p |ZlT>u ylabel('Amplitude')
u;qMo `- #ic 2ofI title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
y]f^`2L!8> A=]F_ axis([0 90 0 1])
H<7D cwXv s
Xk?.A_D grid on
cGzYW~K @Qjl`SL%O^ subplot(1,2,2);
)\\V
s>9 ,T*_mDVY plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
TM}'XZ& gLMea: legend('T_p','T_s','T_n')
fB,1s}3H n ]O=S2Q xlabel('\theta_i')
=C>`}%XT} EZumJ." ylabel('Amplitude')
1( 7.V-(G aevG<|qP title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
R:B BNzY}f 3H}~eEg, axis([0 90 0 1])
V0
+k3H ^@*zH ?Rx{ grid on
3kqV_Pjg Bx : So6: 
\M"^Oe{Dy? +[8Kl=]L 
K M \+
