利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度 ;Kf|a}m -
%5G BMMn
1、光疏射向光密 3CE8+PnT
9NP l]iA)
clear $<cZ<g5)
z6lz*%Yi
close all IR
LPUP
L0kNt
&di
n1=1,n2=1.45; tD]&et
@__;RVQ
theta=0:0.1:90; AD\<}/3U
swlWe}1
a=theta*pi/180; &-fx=gq=
@?m8/t9 .
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); N%f!B"NQ
sAoM=n}!
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); V]]qu:Mh8
`z_7[$\~
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); |y% ].y)
#mhD; .Wg
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); Qu,k
pV6HQ:y1
figure(1) dz|*n 'd
$ rYS
subplot(1,2,1); xLI{= sL
=
Y-Ne6a
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2) Pp_3 n yQ
1fFb 7n~3
legend('r_p','r_s','|r_p|','|r_s|') 1N!g`=}
o0G`Xn
xlabel('\theta_i') <j\;>3Q
uX +<` 3O
ylabel('Amplitude') n[Iu!v\/*
<aaDW
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) ?=!XhU
.
r. b6E% D
axis([0 90 -1 1]) "O[76}I+.q
<cv1$
x ~P
grid on '$XHRS/q]
.crM!{<Y
subplot(1,2,2); (?BgT i\
/i y*3P,`
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2) 5^K#Tj ;2
~H|LWCU)K8
legend('t_p','t_s','|t_p|','|t_s|') lo UwR z
SP*JleQN
xlabel('\theta_i') h
^h-p d
+;* (a3 Gp
ylabel('Amplitude') 0BB @ E(*
8I\er omG
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) Hu9R.[u
; eF4J
axis([0 90 0 1]) P%`R7yk
oPu|Q^I=
grid on V$<G)dwUG5
S"^KJUUc
Rp=abs(rp).^2; ]KK ZbEO
: a q>
Rs=abs(rs).^2; L%K_.!d^
.Spi$>v
Rn=(Rp+Rs)/2; Flujwh@rg
[du>ff
Tp=1-Rp; >3`ctbe
te*Y]-&I|/
Ts=1-Rs; [Ht."VxR
sI RrEea
Tn=(Tp+Ts)/2; :.S41S
H'0*CiHes
figure(2) ]X:
rby$
k<'vP{
subplot(1,2,1); 4 ?@uF[
S`c ]F c
plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2) ?gR\A8:8
22/?JWL>
legend('R_p','R_s','R_n') }1]!#yMfq
C *s0r;
xlabel('\theta_i') UiK+c30FU
-hVv
ylabel('Amplitude') c,+(FQ9
c_z/At;4
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) KBr5bcm4u
kh?#={]Z
axis([0 90 0 1]) e.}3OK
R )d99j^"
grid on K_&c5(-(_
^?6
W<
subplot(1,2,2); g;-+7ViIr
h#@4@x{
plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2) XJJ[F|k~
l<aqiZSY
legend('T_p','T_s','T_n') HhWwc#B
Vgqvv q<S
xlabel('\theta_i') u_4:#~b
#U$ YZ#B
ylabel('Amplitude') /+4^.Q*
%%as>}.
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) 2%5^F i
KO-Zz&2f
axis([0 90 0 1]) ^/%o%J& Hz
7TV>6i+7
grid on tIxhSI^
npO@Haw
d(X\B{
# Z|%0 r_~
yp[<9%Fi
