利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度
/eBcPu"[Vb )+"(7U< 1、光疏射向光密
gk}.L E
]D^zTl3=q clear
=I9 hGj6 a
*bc#!e close all
/GO((v+J -^*8D(j* n1=1,n2=1.45;
p`S~UBcL. Gx |/
Jq theta=0:0.1:90;
P
V9q= KkJE-k*D+w a=theta*pi/180;
; m:I AHT(Z~ C rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
;[7# h8 +3pfBE| rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
w[uw hd )"Wy/P tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2));
m p0 s>R QbYNL9% ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
FC4hvO(/m dCWq~[[ figure(1)
&!*p>Ns)e ;
X/' ujg subplot(1,2,1);
q#0y u"< { !t6&
A plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2)
YCBUc<) B~3qEdoK5` legend('r_p','r_s','|r_p|','|r_s|')
R?] S<Z zB "y^g xlabel('\theta_i')
m.U&O=] 5 =]L#v2@ ylabel('Amplitude')
]Gl5Qf:+z [5]*
Be title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
L @8[. .Pa6HA ! axis([0 90 -1 1])
K14 {c1 %"3tGi : / grid on
i;#AW($+a VKr
oikz@] subplot(1,2,2);
2!a~YT tY?evsVgz plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2)
O.?q8T)n82 =V^8 RlBi legend('t_p','t_s','|t_p|','|t_s|')
?nozB|*>ut A:?w1"7gT xlabel('\theta_i')
"'c
=(P K
st2.Yy ylabel('Amplitude')
c+u) C%g Eqh&<]q title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
oX-h7;SD &wB ?ks axis([0 90 0 1])
4hV~
ir WoWBZ;+ U grid on
iu'r c/=V }}v28"\TA Rp=abs(rp).^2;
cs) z! R(A"6a8* Rs=abs(rs).^2;
v?4MndR y/ah<Y0( Rn=(Rp+Rs)/2;
QsPL^ Ny SG3qNM: g Tp=1-Rp;
M+\LH o(5
( ]bJ Ts=1-Rs;
S]D YEL$ UyOoyyd. Tn=(Tp+Ts)/2;
]m""ga q*{Dy1Tj figure(2)
b %I2ig J_ `\}55n subplot(1,2,1);
O0> ^?dsL [3hOc/]s plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2)
}MV=t7x 9+ :WdiH)Zv legend('R_p','R_s','R_n')
J+8T Ie *m Xs(u xlabel('\theta_i')
\YsYOFc| (yi{<$ U* ylabel('Amplitude')
}Gm/9@oKc 2a O. t title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
J9/}ZD^ (: T\< axis([0 90 0 1])
?!kPW^gD xn 4-^2 grid on
T=->~@5 S%m$LM]NCg subplot(1,2,2);
a(`@u&]WZ qQ UCK plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2)
s:#\U!>0` XqTDLM& legend('T_p','T_s','T_n')
<lwkjt=RV G2}e@L0 xlabel('\theta_i')
/
L/hR4 RgT|^|ZA ylabel('Amplitude')
uvf} 7 7q[a8rUdh title(['n_1=',num2str(n1),',n_2=',num2str(n2)])
V3$Yr"rZ; Q{+&3KXH axis([0 90 0 1])
-:o4|&g<* k\a&4v grid on
)R'~{;z } B @8
]! 
cfd7)(6 u DpCW} 
0 ?
QTi(
