利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度 2@S}x@^
GDLw_usV
1、光疏射向光密 8lQ}-8
%`]+sg[i
clear x/,;:S
Yjoe|
close all oc1BOW z
dN2JOyS
n1=1,n2=1.45; :^7w
sVyV|!K
theta=0:0.1:90; 0F[f%2j
Z],j|rWy6
a=theta*pi/180; LS~at.3zX
]\*^G@HA2
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); .n<vhLDQn
[c{\el9H
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); H07\z1?.K
bq/Aopfr
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); >oW]3)$4S
3aW<FSgP
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); YTTyMn
G9:XEEN
figure(1) .`Rt
@&\Y:aRO%i
subplot(1,2,1); oC5h-4~
HzTmNm)
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2) 7
$y;-[E[
l<6GZ
legend('r_p','r_s','|r_p|','|r_s|') ceUe*}\cr
XPd mz !,b
xlabel('\theta_i') 01r%K@ xX\
x9YQd69
ylabel('Amplitude') 5%}e j)@
$d*9]M4
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) 8w5}9}xF
]oP1c-GEk
axis([0 90 -1 1]) ?i _ACKpw
</ [.1&S+\
grid on C hF~
k%6CkCw
subplot(1,2,2); 2+9VDf2
zsp%Cz7T
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2) l*l*5hA
E@GYl85fI
legend('t_p','t_s','|t_p|','|t_s|') >pF* unC;
UF+Qx/4h0
xlabel('\theta_i') g~DuK|+
3^/w`(-{@
ylabel('Amplitude') M@ U>@x;
?c#s}IH
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) vhBW1/w&F
yCxYFi
axis([0 90 0 1]) E0ED[d,
gGrVpOzBj
grid on 0he3[m}Nr
X.b8qbnq[
Rp=abs(rp).^2; gH//@`6
iVFOOsJ@
Rs=abs(rs).^2; >ai,6!
{;{U@Z
Rn=(Rp+Rs)/2; VM$n|[C~
t'U=K>7
Tp=1-Rp; DP{nvsF
9L>?N:%5
Ts=1-Rs; O=jLZ2os
#
55>?
Tn=(Tp+Ts)/2; BA h'H&;V
YYQvt
figure(2) +(*HDa|
=+iY<~8
subplot(1,2,1); t'eaR-
cQEUHhRg!
plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2) !+SL=xy!{
RRja{*R
legend('R_p','R_s','R_n') L"b&O<No
0x<ASfka
xlabel('\theta_i') {T8;-H0H
lsax.uG5x
ylabel('Amplitude') $9G&
wH>{
KZ 5%q.
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) 'C5id7O&
':n`0+Eh
axis([0 90 0 1]) |S!RQ-CF
o898pg
grid on j:%,lcF
&GLDoLk6[
subplot(1,2,2); ,[;O'g?,g
w-Ph-L/
plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2) %r|sb=(yT
Q+Bl1xl
legend('T_p','T_s','T_n') $9YQ aN%
9Jwd *gevV
xlabel('\theta_i') 3H,x4L5j
wa[L[mw
ylabel('Amplitude') TdPd8ig8{
wa@Rlzij>
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) {8.Zb NEJ
EGl<oxL*R2
axis([0 90 0 1]) "lt <$.
{dF@Vg_n
grid on qxI$F
"w?0f["
%V=%ARP|
:of([e|u6
<{W{
Y\_A>