利用菲涅尔公式计算光波在两种介质表面折反射率及折反射能流密度 ME4 Ir
DG1C_hu
i
1、光疏射向光密 ;[9 WB<t
a M I\gCB/
clear ?ra6Lo
u g 47JW
close all S,A\%:Va
Zt G5vdf
n1=1,n2=1.45; WPLM*]6
2@ 4^ 81
theta=0:0.1:90; yv2BbrYyy
N4y$$.uv2
a=theta*pi/180; qo{2 CYG\+
?L.c~w; l
rp=(n2*cos(a)-n1*sqrt(1-(n1/n2*sin(a)).^2))./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); ~{MmUp rS
2' fg
rs=(n1*cos(a)-n2*sqrt(1-(n1/n2*sin(a)).^2))./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2)); 3N%%69JN)
9;,_Q q
tp=2*n1*cos(a)./(n2*cos(a)+n1*sqrt(1-(n1/n2*sin(a)).^2)); Gc6`]7 s
bj)dYj f
ts=2*n1*cos(a)./(n1*cos(a)+n2*sqrt(1-(n1/n2*sin(a)).^2));
wbg_%h:
^@V$'Bk
figure(1) ^}nz^+R
eTi r-7
subplot(1,2,1); ^ ]Mlkd:
7I.7%m,g
plot(theta,rp,'-',theta,rs,'--',theta,abs(rp),':',theta,abs(rs),'-.','LineWidth',2) pi`sx[T@{Z
~]t/|xep
legend('r_p','r_s','|r_p|','|r_s|') > 9KQWeD
@#sBom+K`
xlabel('\theta_i') LZC) vF5
OFS` ?>
ylabel('Amplitude') Mx&
P^#B3
QvB]?D#h
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) )./pS ~
Sw!/ I PO
axis([0 90 -1 1]) _ElA\L4g%
Ya$JX(aUe
grid on 9D
2B 8t"a
8GC (?#Kb
subplot(1,2,2); 9n ][#I)a3
M+R xt.~6
plot(theta,tp,'-',theta,ts,'--',theta,abs(tp),':',theta,abs(ts),'-.','LineWidth',2) 5 $SO
QC+oSb!!?
legend('t_p','t_s','|t_p|','|t_s|') |UbwPL_L
3)SO-Bz\
xlabel('\theta_i') Y>eypfK"
F,B, D^WD
ylabel('Amplitude') b Y6y)l
b.jxkx\nt
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) Mk- C&#'
# f~,8<K
axis([0 90 0 1]) y7 tK> aD}
h Gc')
grid on Gx75EQ2
NG!~<Kx
Rp=abs(rp).^2; &^-quzlZ
RU`m|<
Rs=abs(rs).^2; "D V.%7*^
G~Oj}rn
Rn=(Rp+Rs)/2; imE5 $;
T""y)%
Tp=1-Rp; |
?Js)i
^s'ozCk 0
Ts=1-Rs; nLZT3`@~,
J('p'SlI
Tn=(Tp+Ts)/2; UH8)r
k]*DuVC OX
figure(2) Dim,HPx]d
eZ~^Z8F[6
subplot(1,2,1); >j]*=&,7
,"/<N*vh
plot(theta,Rp,'-',theta,Rs,'--',theta,Rn,':','LineWidth',2) 9gMNS6D'b
l\l\T<wa,
legend('R_p','R_s','R_n') kTQ`$V(>&
,nO:Pxn|
xlabel('\theta_i') wzAp`Zs2Dm
]1]
ylabel('Amplitude') E.Hw|y0_(|
HaeF`gI^Ee
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) 3mWd?!+m=
Q"]C" ?
axis([0 90 0 1]) @:DS/#!
u(Ke S`
grid on htg+V-,
rnxO2
subplot(1,2,2); -0| '{
A0*u(15%
plot(theta,Tp,'-',theta,Ts,'--',theta,Tn,':','LineWidth',2) 8a{FxCBw
Vxif0Bx&/d
legend('T_p','T_s','T_n') Y W"?F y
* wN+Ak q
xlabel('\theta_i') ?Mp)F2'
O 'JH=
'
ylabel('Amplitude') zqAK|jbL
~..h=
title(['n_1=',num2str(n1),',n_2=',num2str(n2)]) BzH7E[R49
]4@z.1Mr
axis([0 90 0 1]) -W>zON|l
FE1dr_i
grid on N i ^pP@('
*>.~ f<V
n15c1=gs
( CY VSO
NWB/N *
